title: Plan Arithmetic: Compositional Plan Vectors for Multi-Task Control

abstract: Autonomous agents situated in real-world environments must be able to master large repertoires of skills. While a single short skill can be learned quickly, it would be impractical to learn every task independently. Instead, the agent should share knowledge across behaviors such that each task can be learned efficiently, and such that the resulting model can generalize to new tasks, especially ones that are compositions or subsets of tasks seen previously. A policy conditioned on a goal or demonstration has the potential to share knowledge between tasks if it sees enough diversity of inputs. However, these methods may not generalize to a more complex task at test time. We introduce compositional plan vectors (CPVs) to enable a policy to perform compositions of tasks without additional supervision. CPVs represent trajectories as the sum of the subtasks within them. We show that CPVs can be learned within a one-shot imitation learning framework without any additional supervision or information about task hierarchy, and enable a demonstration-conditioned policy to generalize to tasks that sequence twice as many skills as the tasks seen during training. Analogously to embeddings such as word2vec in NLP, CPVs can also support simple arithmetic operations -- for example, we can add the CPVs for two different tasks to command an agent to compose both tasks, without any additional training.


title: Are Out-of-Distribution Detection Methods Effective on Large-Scale Datasets?

abstract: Supervised classification methods often assume the train and test data distributions are the same and that all classes in the test set are present in the training set. However, deployed classifiers often require the ability to recognize inputs from outside the training set as unknowns. This problem has been studied under multiple paradigms including out-of-distribution detection and open set recognition. For convolutional neural networks, there have been two major approaches: 1) inference methods to separate knowns from unknowns and 2) feature space regularization strategies to improve model robustness to outlier inputs. There has been little effort to explore the relationship between the two approaches and directly compare performance on anything other than small-scale datasets that have at most 100 categories. Using ImageNet-1K and Places-434, we identify novel combinations of regularization and specialized inference methods that perform best across multiple outlier detection problems of increasing difficulty level. We found that input perturbation and temperature scaling yield the best performance on large scale datasets regardless of the feature space regularization strategy. Improving the feature space by regularizing against a background class can be helpful if an appropriate background class can be found, but this is impractical for large scale image classification datasets.


title: Dark-Energy Instabilities induced by Gravitational Waves

abstract: We point out that dark-energy perturbations may become unstable in the presence of a gravitational wave of sufficiently large amplitude. We study this effect for the cubic Horndeski operator (braiding), proportional to $\alpha_{\rm B}$. The scalar that describes dark-energy fluctuations features ghost and/or gradient instabilities for gravitational-wave amplitudes that are produced by typical binary systems. Taking into account the populations of binary systems, we conclude that the instability is triggered in the whole Universe for $|\alpha_{\rm B} |\gtrsim 10^{-2}$, i.e. when the modification of gravity is sizeable. The instability is triggered by massive black-hole binaries down to frequencies corresponding to $10^{10}$ km: the instability is thus robust, unless new physics enters on even longer wavelengths. The fate of the instability and the subsequent time-evolution of the system depend on the UV completion, so that the theory may end up in a state very different from the original one. The same kind of instability is present in beyond-Horndeski theories for $|\alpha_{\rm H}| \gtrsim 10^{-20}$. In conclusion, the only dark-energy theories with sizeable cosmological effects that avoid these problems are $k$-essence models, with a possible conformal coupling with matter.


title: A systematic study of stacked square nets: from Dirac fermions to material realizations

abstract: Non-symmorphic symmetries protect Dirac line nodes in square net materials. This phenomenon has been most prominently observed in ZrSiS. Here, we systematically study the symmetry-protected nodal fermions that result from different ways of embedding the square net into a larger unit cell. Surprisingly, we find that a nonsymmorphic space group is not a necessary condition for a filling enforced semimetal: symmorphic space groups can also host nodal fermions that are enforced by band folding and electron count, that is, a combination of a particular structural motif combined with electron filling. We apply the results of this symmetry analysis to define an algorithm, which we utilize to find square net materials with nodal fermions in specific symmorphic space groups. We highlight one result of this search, the compound ThGeSe, which has not been discussed before in the context of nodal fermions. Finally, we discuss how band folding can impose constraints on band connectivity beyond the connectivity of single elementary band representations.


title: Effects of the Great American Solar Eclipse on the lower ionosphere observed with VLF waves

abstract: The altitude of the ionospheric lower layer (D-region) is highly influenced by the solar UV flux affecting in turn, the propagation of Very Low Frequency (VLF) signals inside the waveguide formed between this layer and the Earth surface. A rapid change of the solar irradiance, as during a solar eclipse, can help to understand the details of the energy transfer of the solar radiation onto the ionospheric D-layer. Using the "Latin American VLF Network" (LAVNet-Mex) receiver station in Mexico City, Mexico, we detected the phase and amplitude changes of the VLF signals transmitted by the NDK station at 25.2 kHz in North Dakota, USA during the August 21, 2017, solar eclipse. As the Sunlight was eclipsed, the rate of ionization in the ionosphere (D-region) was reduced and the effective reflection height increased, causing a considerable drop of the phase and amplitude of the observed VLF waves. The corresponding waveguide path is 3007.15 km long and crossed almost perpendicularly the total eclipse path. Circumstantially, at the time of the total eclipse, a C3 flare took place allowing us to isolate the flare flux from the background flux of a large portion of the disk. In this work, we report the observations and present a new model of the ionospheric effects of the eclipse and flare. The model is based on a detailed setup of the degree of Moon shadow that affects the entire Great Circle Path (GCP). During the eclipse, the maximum phase variation was -63.36$^{\circ}$ at 18:05 UT which, according to our model, accounts for a maximum increase of the reflection height of 9.3 km.

categories: physics.space-ph astro-ph.EP astro-ph.SR physics.ao-ph physics.plasm-ph


title: The Soft Quark Sudakov

abstract: There has been recent interest in understanding the all loop structure of the subleading power soft and collinear limits, with the goal of achieving a systematic resummation of subleading power infrared logarithms. Most of this work has focused on subleading power corrections to soft gluon emission, whose form is strongly constrained by symmetries. In this paper we initiate a study of the all loop structure of soft fermion emission. In $\mathcal{N}=1$ QCD we perform an operator based factorization and resummation of the associated infrared logarithms, and prove that they exponentiate into a Sudakov due to their relation to soft gluon emission. We verify this result through explicit calculation to $\mathcal{O}(\alpha_s^3)$. We show that in QCD, this simple Sudakov exponentiation is violated by endpoint contributions proportional to $(C_A-C_F)^n$ which contribute at leading logarithmic order. Combining our $\mathcal{N}=1$ result and our calculation of the endpoint contributions to $\mathcal{O}(\alpha_s^3)$, we conjecture a result for the soft quark Sudakov in QCD, a new all orders function first appearing at subleading power, and give evidence for its universality. Our result, which is expressed in terms of combinations of cusp anomalous dimensions in different color representations, takes an intriguingly simple form and also exhibits interesting similarities to results for large-x logarithms in the off diagonal splitting functions.


title: Anomaly Matching in the Symmetry Broken Phase: Domain Walls, CPT, and the Smith Isomorphism

abstract: Symmetries in Quantum Field Theory may have 't Hooft anomalies. If the symmetry is unbroken in the vacuum, the anomaly implies a nontrivial low-energy limit, such as gapless modes or a topological field theory. If the symmetry is spontaneously broken, for the continuous case, the anomaly implies low-energy theorems about certain couplings of the Goldstone modes. Here we study the case of spontaneously broken discrete symmetries, such as Z/2 and T. Symmetry breaking leads to domain walls, and the physics of the domain walls is constrained by the anomaly. We investigate how the physics of the domain walls leads to a matching of the original discrete anomaly. We analyze the symmetry structure on the domain wall, which requires a careful analysis of some properties of the unbreakable CPT symmetry. We demonstrate the general results on some examples and we explain in detail the mod 4 periodic structure that arises in the Z/2 and T case. This gives a physical interpretation for the Smith isomorphism, which we also extend to more general abelian groups. We show that via symmetry breaking and the analysis of the physics on the wall, the computations of certain discrete anomalies are greatly simplified. Using these results we perform new consistency checks on the infrared phases of 2+1 dimensional QCD.


title: The possible role of stellar mergers for the formation of multiple stellar populations in globular clusters

abstract: Many possible scenarios for the formation of multiple stellar populations (MSP) in globular clusters (GCs) have been discussed so far, including the involvement of asymptotic giant branch stars, fast rotating main sequence stars, very massive main sequence stars and mass-transferring massive binaries based on stellar evolution modelling. But self-consistent, dynamical simulations of very young GCs are usually not considered. In this work, we perform direct $N$-body modelling such systems with total masses up to $3.2\times10^5$ M$_\odot$, taking into account the observationally constrained primordial binary properties, and discuss the stellar-mergers driven both by binary stellar evolution and dynamical evolution of GCs. The occurrence of stellar mergers is enhanced significantly in binary-rich clusters such that stars forming from the gas polluted by mergers-driven ejection/winds would appears as MSPs. We thus emphasize that stellar mergers can be an important process that connects MSP formation with star cluster dynamics, and that multiple MSP formation channels can naturally work together. The scenario studied here, also in view of a possible top-heavy IMF, may be particularly relevant for explaining the high mass fraction of MSPs (the mass budget problem) and the absence of MSPs in young and low-mass star clusters.


title: Simulating the interstellar medium of galaxies with radiative transfer, non-equilibrium thermochemistry, and dust

abstract: We present a novel framework to self-consistently model the effects of radiation fields, dust physics and molecular chemistry (H$_2$) in the interstellar medium (ISM) of galaxies. The model combines a state-of-the-art radiation hydrodynamics module with a non-equilibrium thermochemistry module that accounts for H$_2$ coupled to a realistic dust formation and destruction model, all integrated into the new stellar feedback framework SMUGGLE. We test this model on high-resolution isolated Milky-Way (MW) simulations. We show that photoheating from young stars makes stellar feedback more efficient, but this effect is quite modest in low gas surface density galaxies like the MW. The multi-phase structure of the ISM, however, is highly dependent on the strength of the interstellar radiation field. We are also able to predict the distribution of H$_2$, that allow us to match the molecular Kennicutt-Schmidt (KS) relation, without calibrating for it. We show that the dust distribution is a complex function of density, temperature and ionization state of the gas which cannot be reproduced by simple scaling relations often used in the literature. Our model is only able to match the observed dust temperature distribution if radiation from the old stellar population is considered, implying that these stars have a non-negligible contribution to dust heating in the ISM. Our state-of-the-art model is well-suited for performing next generation cosmological galaxy formation simulations, which will be able to predict a wide range of resolved ($\sim 10$ pc) properties of galaxies.


title: HyDRA-H: Simultaneous Hybrid Retrieval of Exoplanetary Emission Spectra

abstract: High-resolution Doppler spectroscopy has been used to detect several chemical species in exoplanetary atmospheres. Such detections have traditionally relied on cross correlation of observed spectra against spectral model templates, an approach that is successful for detecting chemical species but not optimised for constraining abundances. Recent work has explored ways to perform atmospheric retrievals on high-resolution spectra (HRS) and combine them with retrievals routinely performed for low-resolution spectra (LRS) by developing a mapping from the cross correlation function to a likelihood metric. We build upon previous studies and report HyDRA-H, a hybrid retrieval code for simultaneous analysis of low- and high- resolution thermal emission spectra of exoplanets in a fully Bayesian approach. We demonstrate HyDRA-H on the hot Jupiter HD 209458b as a case study. We validate our HRS retrieval capability by confirming previous results and report a simultaneous hybrid retrieval using both HRS and LRS data. The LRS data span the HST WFC3 (1.1-1.7 $\mu$m) and Spitzer photometry (IRAC 3.6-8 $\mu$m) bands, while the HRS data were obtained with CRIRES on VLT at 2.3 $\mu$m. The constraints on the composition and temperature profiles for the hybrid retrieval are more stringent than retrievals with either LRS or HRS datasets individually. We retrieve abundances of $\log(\mathrm{H_2O)} = -4.11^{+0.91}_{-0.30}$ and $\log(\mathrm{CO}) = {-2.16}^{+0.99}_{-0.47}$, and $\mathrm{C/O} = 0.99^{+0.01}_{-0.02}$, consistent with previous works. We constrain the photospheric temperature to be $1498^{+216}_{-57}$ K, consistent with the equilibrium temperature. Our results demonstrate the significant advantages of hybrid retrievals by combining strengths of both HRS and LRS observations which probe complementary aspects of exoplanetary atmospheres.


title: Quantum aspects of "hydrodynamic" transport from weak electron-impurity scattering

abstract: Recent experimental observations of apparently hydrodynamic electronic transport have generated much excitement. However, the understanding of the observed non-local transport (whirlpool) effects and parabolic (Poiseuille-like) current profiles has largely been motivated by a phenomenological analogy to classical fluids. This is due to difficulty in incorporating strong correlations in quantum mechanical calculation of transport, which has been the primary angle for interpreting the apparently hydrodynamic transport. Here we demonstrate that even free fermion systems, in the presence of (inevitable) disorder, exhibit non-local conductivity effects such as those observed in experiment because of the fermionic system's long-range entangled nature. On the basis of explicit calculations of the conductivity at finite wavevector, $\sigma({\bf q})$, for selected weakly disordered free fermion systems, we propose experimental strategies for demonstrating distinctive quantum effects in non-local transport at odds with the expectations of classical kinetic theory. Our results imply that the observation of whirlpools or other "hydrodynamic" effects does not guarantee the dominance of electron-electron scattering over electron-impurity scattering.


title: Relativistic kinetic gases as direct sources of gravity

abstract: We propose a new model for the description of a gravitating multi particle system, viewed as a kinetic gas. The properties of the, colliding or non-colliding, particles are encoded into a so called one-particle distribution function, which is a density on the space of allowed particle positions and velocities, i.e. on the tangent bundle of the spacetime manifold. We argue that an appropriate theory of gravity, describing the gravitational field generated by a kinetic gas, must also be modeled on the tangent bundle. The most natural mathematical framework for this task is Finsler spacetime geometry. Following this line of argumentation, we construct a coupling between the kinetic gas and a recently proposed Finsler geometric extension of general relativity. Additionally, we explicitly show how the coordinate invariance of the action of the kinetic gas leads to a novel formulation of conservation of the energy-momentum distribution of the gas on the tangent bundle.


title: The Dragonfly Wide Field Survey. I. Telescope, Survey Design and Data Characterization

abstract: We present a description of the Dragonfly Wide Field Survey (DWFS), a deep photometric survey of a wide area of sky. The DWFS covers 330 $\mathrm{deg}^2$ in the equatorial GAMA fields and the Stripe 82 fields in the SDSS $g$ and $r$ bands. It is carried out with the 48-lens Dragonfly Telephoto Array, a telescope that is optimized for the detection of low surface brightness emission. The main goal of the survey is to study the dwarf galaxy population beyond the Local Group. In this paper, we describe the survey design and show early results. We reach $1\sigma$ depths of $\mu_g\approx 31$ mag arcsec$^{-2}$ on arcminute scales and show that Milky Way satellites such as Sextans, Bootes, and Ursa Major should be detectable out to $D\gtrsim 10$ Mpc. We also provide an overview of the elements and operation of the 48-lens Dragonfly telescope and a detailed description of its data reduction pipeline. The pipeline is fully automated, with individual frames subjected to a rigorous series of quality tests. The sky subtraction is performed in two stages, ensuring that emission features with spatial scales up to $\sim 0.^{\circ}9 \times 0.^{\circ}6$ are preserved. The DWFS provides unparalleled sensitivity to low surface brightness features on arcminute scales.


title: Decorated $\mathbb{Z}_{2}$ Symmetry Defects and Their Time-Reversal Anomalies

abstract: We discuss an isomorphism between the possible anomalies of $(d+1)$-dimensional quantum field theories with $\mathbb{Z}_{2}$ unitary global symmetry, and those of $d$-dimensional quantum field theories with time-reversal symmetry $\mathsf{T}$. This correspondence is an instance of symmetry defect decoration. The worldvolume of a $\mathbb{Z}_{2}$ symmetry defect is naturally invariant under $\mathsf{T},$ and bulk $\mathbb{Z}_{2}$ anomalies descend to $\mathsf{T}$ anomalies on these defects. We illustrate this correspondence in detail for $(1+1)d$ bosonic systems where the bulk $\mathbb{Z}_{2}$ anomaly leads to a Kramers degeneracy in the symmetry defect Hilbert space, and exhibit examples. We also discuss $(1+1)d$ fermion systems protected by $\mathbb{Z}_{2}$ global symmetry where interactions lead to a $\mathbb{Z}_{8}$ classification of anomalies. Under the correspondence, this is directly related to the $\mathbb{Z}_{8}$ classification of $(0+1)d$ fermions protected by $\mathsf{T}$. Finally, we consider $(3+1)d$ bosonic systems with $\mathbb{Z}_{2}$ symmetry where the possible anomalies are classified by $\mathbb{Z}_{2}\times \mathbb{Z}_{2}$. We construct topological field theories realizing these anomalies and show that their associated symmetry defects support anyons that can be either fermions or Kramers doublets.


title: A Potential Mechanism for Inflation from Swampland Conjectures

abstract: Inflation is the currently accepted paradigm for the beginnings of the Universe. To explain the observed almost scale invariant spectrum of density perturbations with only a slight spectral tilt, inflation must have been "slow roll", that is with a potential with sufficiently small slope. While the origin of inflationary structure is intrinsically quantum mechanical, gravity gets treated semiclassically within inflationary models. Recent work, in terms of the so called de-Sitter swampland conjecture, has called into question whether slow roll inflation is consistent with a complete theory of quantum gravity in the presence of a positive vacuum energy density, which is a key ingredient in the inflationary paradigm. In this work, we show that, in fact, if we understand this conjecture correctly and with another swampland conjecture, the so-called distance conjecture, involved we get a potential mechanism for slow roll inflation and we argue that here fine-tuning is not a technical problem.


title: Thermalization and its absence within Krylov subspaces of a constrained Hamiltonian

abstract: We study the quantum dynamics of a simple translation invariant, center-of-mass (CoM) preserving model of interacting fermions in one dimension (1D), which arises in multiple experimentally realizable contexts. We show that this model naturally displays the phenomenology associated with fractonic systems, wherein single charges can only move by emitting dipoles. This allows us to demonstrate the rich Krylov fractured structure of this model, whose Hilbert space shatters into exponentially many dynamically disconnected subspaces. Focusing on exponentially large Krylov subspaces, we show that these can be either be integrable or non-integrable, thereby establishing the notion of Krylov-restricted thermalization. We analytically find a tower of integrable Krylov subspaces of this Hamiltonian, all of which map onto spin-1/2 XX models of various system sizes. We also discuss the physics of the non-integrable subspaces, where we show evidence for weak Eigenstate Thermalization Hypothesis (ETH) restricted to each non-integrable Krylov subspace. Further, we show that constraints in some of the thermal Krylov subspaces cause the long-time expectation values of local operators to deviate from behavior typically expected from translation-invariant systems. Finally, we show using a Schrieffer-Wolff transformation that such models naturally appear as effective Hamiltonians in the large electric field limit of the interacting Wannier-Stark problem, and comment on connections of our work with the phenomenon of Bloch many-body localization.


title: Inverse Compton Signatures of Gamma-Ray Burst Afterglows

abstract: The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~GeV) remains uncertain. The recent detection of sub-TeV emission from GRB~190114C by MAGIC raises further debate on what powers the very high-energy (VHE; $\gtrsim 300$GeV) emission. Here, we explore the inverse Compton scenario as a candidate for the HE and VHE emissions, considering two sources of seed photons for scattering: synchrotron photons from the blast wave (synchrotron self-Compton or SSC) and isotropic photon fields external to the blast wave (external Compton). For each case, we compute the multi-wavelength afterglow spectra and light curves. We find that SSC will dominate particle cooling and the GeV emission, unless a dense ambient infrared photon field, typical of star-forming regions, is present. Additionally, considering the extragalactic background light attenuation, we discuss the detectability of VHE afterglows by existing and future gamma-ray instruments for a wide range of model parameters. Studying GRB~190114C, we find that its afterglow emission in the \fermi-LAT band is synchrotron-dominated.The late-time \fermi-LAT measurement (i.e., $t\sim 10^4$~s), and the MAGIC observation also set an upper limit on the energy density of a putative external infrared photon field (i.e. $\lesssim 3\times 10^{-9}\,{\rm erg\,cm^{-3}}$), making the inverse Compton dominant in the sub-TeV energies.


title: Electron-Positron Annihilation Freeze-Out in the Early Universe

abstract: Electron-positron annihilation largely occurs in local thermal and chemical equilibrium after the neutrinos fall out of thermal equilibrium and during the Big Bang Nucleosynthesis (BBN) epoch. The effects of this process are evident in BBN yields as well as the relativistic degrees of freedom. We self-consistently calculate the collision integral for electron-positron creation and annihilation using the Klein-Nishina amplitude and appropriate statistical factors for Fermi-blocking and Bose-enhancement. Our calculations suggest that this annihilation freezes out when the photon-electron-positron-baryon plasma temperature is approximately 16 keV, after which its rate drops below the Hubble rate. In the temperature regime near 16 keV, we break the assumption of chemical equilibrium between the electrons, positrons, and photons to independently calculate the evolution of the chemical potentials of the electrons and positrons while computing the associated collision integrals at every time step. We find that the electron and positron chemical potentials deviate from the case with chemical equilibrium. While our results do not affect the interpretation of precision cosmological measurements in elucidating the standard cosmological model, these out of equilibrium effects may be important for testing physics beyond the standard model.


title: The Cosmological Bootstrap: Weight-Shifting Operators and Scalar Seeds

abstract: A key insight of the bootstrap approach to cosmological correlations is the fact that all correlators of slow-roll inflation can be reduced to a unique building block---the four-point function of conformally coupled scalars, arising from the exchange of a massive scalar. Correlators corresponding to the exchange of particles with spin are then obtained by applying a spin-raising operator to the scalar-exchange solution. Similarly, the correlators of massless external fields can be derived by acting with a suitable weight-raising operator. In this paper, we present a systematic and highly streamlined derivation of these operators (and their generalizations) using tools of conformal field theory. Our results greatly simplify the theoretical foundations of the cosmological bootstrap program.


title: The Myriad Uses of Instantons

abstract: In quantum chromodynamics (QCD), the role which topologically non-trivial configurations play in splitting the singlet pseudo-Goldstone meson, the $\eta^\prime$, from the octet is familiar. In addition, such configurations contribute to other processes which violate the axial $U(1)_A$ symmetry. While the nature of topological fluctuations in the confined phase is still unsettled, at temperatures well above that for the chiral phase transition, they can be described by a dilute gas of instantons. We show that instantons of arbitrary topological charge $Q$ generate anomalous interactions between $2 N_f |Q|$ quarks, which for $Q = 1$ make the $\eta^\prime$ heavy. For two flavors we compute an anomalous quartic meson coupling and discuss its implications for the phenomenology of the chiral phase transition. A dilute instanton gas suggests that for cold, dense quarks, instantons do not evaporate until very high densities, when the baryon chemical potential is $\gtrsim 2$ GeV.


title: Instanton Resummation and the Weak Gravity Conjecture

abstract: We develop methods for resummation of instanton lattice series. Using these tools, we investigate the consequences of the Weak Gravity Conjecture for large-field axion inflation. We find that the Sublattice Weak Gravity Conjecture implies a constraint on the volume of the axion fundamental domain. However, we also identify conditions under which alignment and clockwork constructions, and a new variant of N-flation that we devise, can evade this constraint. We conclude that some classes of low-energy effective theories of large-field axion inflation are consistent with the strongest proposed form of the Weak Gravity Conjecture, while others are not.


title: uGMRT search for cold gas at z~1-1.4 towards red quasars

abstract: We present results from our search for HI 21-cm and OH 18-cm absorption at z~1-1.4 towards red quasars showing strong MgII absorption using uGMRT. The quasars J1501+1822 and J1521+5508 show multiple strong associated MgII absorption at z~1.1 and signature of reddening in their optical spectra. We report the detection of HI 21-cm absorption towards J1521+5508 at the systemic redshift of the quasar, with N(HI) =(1.2 +/- 0.2) x 10^{20} cm^{-2} for spin temperature of 100 K and unit covering factor. The HI 21-cm absorption is offset from the blueshifted strong MgII absorbers by >~1500 km/s. We do not detect HI 21-cm absorption at the redshift of the associated MgII absorption and at the systemic redshift towards J1501+1822. We argue that lack of one-to-one correspondence between MgII and HI 21-cm absorption could be related with clumpiness of the neutral gas and the radio and optical sightlines probing different volume of the gas. We find that the presence of strong associated MgII absorption and reddening in the optical spectra of the quasars lead to an increased detection rate of associated HI 21-cm absorption at z>=1. We also report non-detection of intervening OH absorption ([OH]/[HI] <= (1-4) x 10^{-8}) at z = 1.3 towards two red quasars, J0850+5159 and J0852+3435, which show strong metal and HI 21-cm absorption and the 2175 A dust extinction bump.


abstract: Solving goal-oriented tasks is an important but challenging problem in reinforcement learning (RL). For such tasks, the rewards are often sparse, making it difficult to learn a policy effectively. To tackle this difficulty, we propose a new approach called Policy Continuation with Hindsight Inverse Dynamics (PCHID). This approach learns from Hindsight Inverse Dynamics based on Hindsight Experience Replay, enabling the learning process in a self-imitated manner and thus can be trained with supervised learning. This work also extends it to multi-step settings with Policy Continuation. The proposed method is general, which can work in isolation or be combined with other on-policy and off-policy algorithms. On two multi-goal tasks GridWorld and FetchReach, PCHID significantly improves the sample efficiency as well as the final performance.


title: Unsupervised Star Galaxy Classification with Cascade Variational Auto-Encoder

abstract: The increasing amount of data in astronomy provides great challenges for machine learning research. Previously, supervised learning methods achieved satisfactory recognition accuracy for the star-galaxy classification task, based on manually labeled data set. In this work, we propose a novel unsupervised approach for the star-galaxy recognition task, namely Cascade Variational Auto-Encoder (CasVAE). Our empirical results show our method outperforms the baseline model in both accuracy and stability.


title: Observational Properties of Field UDGs: Colours and Number Densities

abstract: While much of the focus around Ultra-Diffuse Galaxies (UDGs) has been given to those in galaxy groups and clusters, relatively little is known about them in less-dense environments. These isolated UDGs provide fundamental insights into UDG formation because environmentally driven evolution and survivability play less of a role in determining their physical and observable properties. We have recently conducted a statistical analysis of UDGs in the field using a new catalogue of sources detected in the deep Kilo-Degree Survey (KiDS) and Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) optical imaging surveys. Using an empirical model to assess our contamination from interloping sources, we show that a scenario in which cluster-like quiescent UDGs occupy a large fraction of the field UDG population is unlikely, with most being significantly bluer and some showing signs of localised star formation. We estimate an upper-limit on the total field abundance of UDGs of 8$\pm$3$\times10^{-3}$cMpc$^{-3}$ within our selection range. The mass formation efficiency of UDGs implied by this upper-limit is similar to what is measured in groups and clusters, meaning that secular formation channels may significantly contribute to the overall UDG population.


title: How mergers magnetise massive stars

abstract: Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars (<1.5 solar masses) are generated in-situ by a dynamo process in their turbulent, convective envelopes. Unlike such stars, intermediate-mass and high-mass stars (>1.5 solar masses; referred to as "massive" stars here) have relatively quiet, radiative envelopes where a solar-like dynamo cannot operate. However, about 10% of them, including 78 Vir, have strong, large-scale surface magnetic fields whose origin has remained a major mystery. The massive star $\tau$ Sco is a prominent member of this group and appears to be surprisingly young compared to other presumably coeval members of the Upper Scorpius association. Here, we present the first 3D magneto-hydrodynamical simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain $\tau$ Sco's magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. These stars are promising progenitors of those neutron stars that host the strongest magnetic fields in the Universe, so-called magnetars, and that may give rise to some of the enigmatic fast radio bursts. Strong magnetic fields affect the explosions of core-collapse supernovae and, moreover, those magnetic stars that have rapidly-rotating cores at the end of their lives might provide the right conditions to power long-duration gamma-ray bursts and super-luminous supernovae.


title: Intensity mapping: a new window into the cosmos

abstract: The technique of intensity mapping (IM) has emerged as a powerful tool to explore the universe at $z < 6$. IM measures the integrated emission from sources over a broad range of frequencies, unlocking significantly more information than traditional galaxy surveys. Astrophysical uncertainties, however, constitute an important systematic in our attempts to constrain cosmology with IM. I describe an innovative approach which allows us to fully utilize our current knowledge of astrophysics in order to develop cosmological forecasts from IM. This framework can be used to exploit synergies with other complementary surveys, thereby opening up the fascinating possibility of constraining physics beyond $\Lambda$CDM from future IM observations.


title: Brill-Noether and existence of semistable sheaves on del Pezzo surfaces

abstract: Let $X$ be a del Pezzo surface. When the degree of $X$ is at least 4, we compute the cohomology of a general sheaf in the moduli space of Gieseker semistable sheaves. We also classify the Chern characters for which the general sheaf in the moduli space is non-special, i.e. has at most one nonzero cohomology group. Our results hold for arbitrary polarizations, slope semistability, and semi-exceptional moduli spaces. When the degree of $X$ is at least 3, we further show our construction of certain vector bundles implies the existence of stable and semistable sheaves with respect to the anti-canonical polarization.


title: Moir\'e quantum chemistry: charge transfer in transition metal dichalcogenide superlattices

abstract: Transition metal dichalcogenide (TMD) bilayers have recently emerged as a robust and tunable moir\'e system for studying and designing correlated electron physics. In this work, by combining large-scale first principle calculation and continuum model approach, we provide an electronic structure theory that maps long-period heterobilayer TMD superlattices onto diatomic crystals with cations and anions. We find that the interplay between moir\'e potential and Coulomb interaction leads to filling-dependent charge transfer between MM and MX regions several nanometers apart. We show that the insulating state at half-filling found in recent experiments on WSe$_2$/WS$_2$ is a charge-transfer insulator rather than a Mott-Hubbard insulator. Our work reveals the richness of simplicity in moir\'e quantum chemistry.


title: Anomalous Magnetic Moments from Asymptotic Safety

abstract: The measurements of the muon and electron anomalous magnetic moments hint at physics beyond the standard model. We show why and how models inspired by asymptotic safety can explain deviations from standard model predictions naturally. Our setup features an enlarged scalar sector and Yukawa couplings between leptons and new vector-like fermions. Using the complete two-loop running of couplings, we observe a well-behaved high energy limit of models including a stabilization of the Higgs. We find that a manifest breaking of lepton universality beyond standard model Yukawas is not necessary to explain the muon and electron anomalies. We further predict the tau anomalous magnetic moment, and new particles in the TeV energy range whose signatures at colliders are indicated. With small CP phases, the electron EDM can be as large as the present bound.


title: LaplacianNet: Learning on 3D Meshes with Laplacian Encoding and Pooling

abstract: 3D models are commonly used in computer vision and graphics. With the wider availability of mesh data, an efficient and intrinsic deep learning approach to processing 3D meshes is in great need. Unlike images, 3D meshes have irregular connectivity, requiring careful design to capture relations in the data. To utilize the topology information while staying robust under different triangulation, we propose to encode mesh connectivity using Laplacian spectral analysis, along with Mesh Pooling Blocks (MPBs) that can split the surface domain into local pooling patches and aggregate global information among them. We build a mesh hierarchy from fine to coarse using Laplacian spectral clustering, which is flexible under isometric transformation. Inside the MPBs there are pooling layers to collect local information and multi-layer perceptrons to compute vertex features with increasing complexity. To obtain the relationships among different clusters, we introduce a Correlation Net to compute a correlation matrix, which can aggregate the features globally by matrix multiplication with cluster features. Our network architecture is flexible enough to be used on meshes with different numbers of vertices. We conduct several experiments including shape segmentation and classification, and our LaplacianNet outperforms state-of-the-art algorithms for these tasks on ShapeNet and COSEG datasets.


title: Infinite Switch Simulated Tempering in Force (FISST)

abstract: Many proteins in cells are capable of sensing and responding to piconewton scale forces, a regime in which conformational changes are small but significant for biological processes. In order to efficiently and effectively sample the response of these proteins to small forces, enhanced sampling techniques will be required. In this work, we derive, implement, and evaluate an efficient method to simultaneously sample the result of applying any constant pulling force within a specified range to a molecular system of interest. We start from Simulated Tempering in Force, whereby force is applied as a linear bias on a collective variable to the system's Hamiltonian, and the coefficient is taken as a continuous auxiliary degree of freedom. We derive a formula for an average collective-variable-dependent force, which depends on a set of weights, learned on-the-fly throughout a simulation, that reflect the limit where force varies infinitely quickly. These weights can then be used to retroactively compute averages of any observable at any force within the specified range. This technique is based on recent work deriving similar equations for Infinite Switch Simulated Tempering in Temperature, that showed the infinite switch limit is the most efficient for sampling. Here, we demonstrate that our method accurately and simultaneously samples molecular systems at all forces within a user defined force range, and show how it can serve as an enhanced sampling tool for cases where the pulling direction destabilizes states of low free-energy at zero-force. This method is implemented in, and will be freely-distributed with, the PLUMED open-source sampling library, and hence can be readily applied to problems using a wide range of molecular dynamics software packages.


title: Whittaker functions from motivic Chern classes

abstract: We prove a `motivic' analogue of the Weyl character formula, computing the Euler characteristic of a line bundle on a generalized flag manifold $G/B$ multiplied either by a motivic Chern class of a Schubert cell, or a Segre analogue of it. The result, given in terms of Demazure-Lusztig (D-L) operators, recovers formulas found by Brubaker, Bump and Licata for the Iwahori-Whittaker functions of the principal series representation of a $p$-adic group. In particular, we obtain a new proof of the classical Casselman-Shalika formula for the spherical Whittaker function. The proofs are based on localization in equivariant K theory, and require a geometric interpretation of how the Hecke dual (or inverse) of a D-L operator acts on the class of a point. We prove that the Hecke dual operators give Grothendieck-Serre dual classes of the motivic classes, a result which might be of independent interest. In an Appendix joint with Dave Anderson we show that if the line bundle is trivial, we recover a generalization of a classical formula by Kostant, Macdonald, Shapiro and Steinberg for the Poincar{\'e} polynomial of $G/B$; the generalization we consider is due to Aky{\i}ld{\i}z and Carrell and replaces $G/B$ by any smooth Schubert variety.


title: The factorization bias in the Van der Meer method: Run 2 experiences at the CMS experiment

abstract: The luminosity measurement of the CMS experiment at the CERN LHC is calibrated with Van der Meer (VdM) scans. A bias occurs in the VdM method due to the assumption of transversely factorizable proton densities of the LHC beams. Here, the different methods applied in Run 2 to estimate the size of the factorization bias are reported. The beam-imaging method reconstructs the transverse proton densities from beam-imaging scans. Additional methods exploit offset scans, analyze the evolution of the measured luminous region, and evaluate diagonal scans.


title: Spectral properties of kernel matrices in the flat limit

abstract: Kernel matrices are of central importance to many applied fields. In this manuscript, we focus on spectral properties of kernel matrices in the so-called "flat limit", which occurs when points are close together relative to the scale of the kernel. We establish asymptotic expressions for the determinants of the kernel matrices, which we then leverage to obtain asymptotic expressions for the main terms of the eigenvalues. Analyticity of the eigenprojectors yields expressions for limiting eigenvectors, which are strongly tied to discrete orthogonal polynomials. Both smooth and finitely smooth kernels are covered, with stronger results available in the finite smoothness case.


title: High resolution time- and angle-resolved photoemission spectroscopy with 11 eV laser pulses

abstract: Performing time and angle resolved photoemission spectroscopy (tr-ARPES) at high momenta necessitates extreme ultraviolet laser pulses, which are typically produced via high harmonic generation (HHG). Despite recent advances, HHG-based setups still require large pulse energies (hundreds of $\mu$J to mJ) and their energy resolution is limited to tens of meV. Here, we present a novel 11 eV tr-ARPES setup that generates a flux of $5\times10^{10}$ photons/s and achieves an unprecedented energy resolution of 16 meV. It can be operated at high repetition rates (up to 250 kHz) while using input pulse energies down to 3 $\mu$J. We demonstrate these unique capabilities by simultaneously capturing the energy and momentum resolved dynamics in two well-separated momentum space regions of a charge density wave material ErTe$_3$. This novel setup offers opportunity to study the non-equilibrium band structure of solids with exceptional energy and time resolutions at high repetition rates.


title: The magnon-mediated attraction between two holes doped in a CuO$_2$ layer

abstract: Using a realistic multi-band model for two holes doped into a CuO$_2$ layer, we devise a method to turn off the magnon-mediated interaction between the holes. This allows us to verify that this interaction is attractive, and therefore could indeed be (part of) the superconducting glue. We derive its analytical expression and show that it consists of a novel kind of pair-hopping+spin-exchange terms. Its coupling constant is fitted from the ground state energy obtained with variational exact diagonalization, and it faithfully reproduces the effect of the magnon-mediated attraction in the entire Brillouin zone. For realistic parameter values, this effective interaction is borderline strong enough to bind the holes into preformed pairs.


title: Comment on "Convergence towards asymptotic state in 1-D mappings: A scaling investigation"

abstract: Nonequilibrium phase transitions are characterized by the so-called critical exponents, each of which is related to a different observable. Systems that share the same set of values for these exponents also share the same universality class. Thus, it is important that the exponents are named and treated in a standardized framework. In this comment, we reinterpret the exponents obtained in [Phys Lett A 379:1246-12 (2015)] for the logistic and cubic maps in order to correctly state the universality class of their bifurcations, since these maps may describe the mean-field solution of stochastic spreading processes.

categories: nlin.AO cond-mat.dis-nn cond-mat.stat-mech physics.bio-ph


title: On the Stability of Self-similar Blow-up for $C^{1,\alpha}$ Solutions to the Incompressible Euler Equations on $\mathbb{R}^3$

abstract: We study the stability of recently constructed self-similar blow-up solutions to the incompressible Euler equation. A consequence of our work is the existence of finite-energy $C^{1,\alpha}$ solutions that become singular in finite time in a locally self-similar manner. As a corollary, we also observe that the Beale-Kato-Majda criterion cannot be improved in the class of $C^{1,\alpha}$ solutions.


title: Unsupervised inference approach to facial attractiveness

abstract: The perception of facial beauty is a complex phenomenon depending on many, detailed and global facial features influencing each other. In the machine learning community this problem is typically tackled as a problem of supervised inference. However, it has been conjectured that this approach does not capture the complexity of the phenomenon. A recent original experiment (Ib\'a\~nez-Berganza et al., Scientific Reports 9, 8364, 2019) allowed different human subjects to navigate the face-space and ``sculpt'' their preferred modification of a reference facial portrait. Here we present an unsupervised inference study of the set of sculpted facial vectors in that experiment. We first infer minimal, interpretable, and faithful probabilistic models (through Maximum Entropy and artificial neural networks) of the preferred facial variations, that capture the origin of the observed inter-subject diversity in the sculpted faces. The application of such generative models to the supervised classification of the gender of the sculpting subjects, reveals an astonishingly high prediction accuracy. This result suggests that much relevant information regarding the subjects may influence (and be elicited from) her/his facial preference criteria, in agreement with the multiple motive theory of attractiveness proposed in previous works.


title: A Primal-dual weak Galerkin finite element method for linear convection equations in non-divergence form

abstract: A new primal-dual weak Galerkin (PD-WG) finite element method was developed and analyzed in this article for first-order linear convection equations in non-divergence form. The PD-WG method results in a symmetric discrete system involving not only the original equation for the primal variable, but also the dual/adjoint equation for the dual variable (also known as Lagrangian multiplier). Optimal order of error estimates in various discrete Sobolev norms are derived for the numerical solutions arising from the PD-WG scheme. Numerical results are produced and reported to illustrate the accuracy and effectiveness of the new PD-WG method.


title: Regularity dimensions: quantifying doubling and uniform perfectness

abstract: We study the upper and lower regularity dimensions in relation to the notions of doubling and uniformly perfect. These two regularity properties are closely related which is quantified thanks to the regularity dimensions. The regularity dimensions of pushforward measures onto graphs of Brownian motion are calculated, similarly for pushforwards with respect to quasisymmetric homeomorphisms. We finish by introducing an application to Diophantine approximation in the setting of Kleinian groups.


title: Fill in the Blanks: Imputing Missing Sentences for Larger-Context Neural Machine Translation

abstract: Most neural machine translation systems still translate sentences in isolation. To make further progress, a promising line of research additionally considers the surrounding context in order to provide the model potentially missing source-side information, as well as to maintain a coherent output. One difficulty in training such larger-context (i.e. document-level) machine translation systems is that context may be missing from many parallel examples. To circumvent this issue, two-stage approaches, in which sentence-level translations are post-edited in context, have recently been proposed. In this paper, we instead consider the viability of filling in the missing context. In particular, we consider three distinct approaches to generate the missing context: using random contexts, applying a copy heuristic or generating it with a language model. In particular, the copy heuristic significantly helps with lexical coherence, while using completely random contexts hurts performance on many long-distance linguistic phenomena. We also validate the usefulness of tagged back-translation. In addition to improving BLEU scores as expected, using back-translated data helps larger-context machine translation systems to better capture long-range phenomena.


title: A Neural Topic-Attention Model for Medical Term Abbreviation Disambiguation

abstract: Automated analysis of clinical notes is attracting increasing attention. However, there has not been much work on medical term abbreviation disambiguation. Such abbreviations are abundant, and highly ambiguous, in clinical documents. One of the main obstacles is the lack of large scale, balance labeled data sets. To address the issue, we propose a few-shot learning approach to take advantage of limited labeled data. Specifically, a neural topic-attention model is applied to learn improved contextualized sentence representations for medical term abbreviation disambiguation. Another vital issue is that the existing scarce annotations are noisy and missing. We re-examine and correct an existing dataset for training and collect a test set to evaluate the models fairly especially for rare senses. We train our model on the training set which contains 30 abbreviation terms as categories (on average, 479 samples and 3.24 classes in each term) selected from a public abbreviation disambiguation dataset, and then test on a manually-created balanced dataset (each class in each term has 15 samples). We show that enhancing the sentence representation with topic information improves the performance on small-scale unbalanced training datasets by a large margin, compared to a number of baseline models.


title: Quantum dephasing induced by non-Markovian random telegraph noise

abstract: We theoretically study the dynamical dephasing of a quantum two level system interacting with an environment exhibiting non-Markovian random telegraph fluctuations. The time evolution of the conditional probability of the environmental noise is governed by a generalized master equation depending on the environmental memory effect. The expression of the dephasing factor is derived exactly which is closely associated with the memory kernel in the generalized master equation for the conditional probability of the environmental noise. In terms of three important types memory kernels, we discuss the quantum dephasing dynamics of the system and the non-Markovian character exhibiting in the dynamical dephasing induced by non-Markovian random telegraph noise. We show that the dynamical dephasing of the quantum system does not always exhibit non-Markovian character which results from that the non-Markovian character in the dephasing dynamics depends both on the environmental non-Markovian character and the interaction between the system and environment. In addition, the dynamical dephasing of the quantum system can be modulated by the external modulation frequency of the environment. This result is significant to quantum information processing and helpful for further understanding non-Markovian dynamics of open quantum systems.


title: Probabilistic Detection and Estimation of Conic Sections from Noisy Data

abstract: Inferring unknown conic sections on the basis of noisy data is a challenging problem with applications in computer vision. A major limitation of the currently available methods for conic sections is that estimation methods rely on the underlying shape of the conics (being known to be ellipse, parabola or hyperbola). A general purpose Bayesian hierarchical model is proposed for conic sections and corresponding estimation method based on noisy data is shown to work even when the specific nature of the conic section is unknown. The model, thus, provides probabilistic detection of the underlying conic section and inference about the associated parameters of the conic section. Through extensive simulation studies where the true conics may not be known, the methodology is demonstrated to have practical and methodological advantages relative to many existing techniques. In addition, the proposed method provides probabilistic measures of uncertainty of the estimated parameters. Furthermore, we observe high fidelity to the true conics even in challenging situations, such as data arising from partial conics in arbitrarily rotated and non-standard form, and where a visual inspection is unable to correctly identify the type of conic section underlying the data.


title: Implicit multirate GARK methods

abstract: This work considers multirate generalized-structure additively partitioned Runge-Kutta (MrGARK) methods for solving stiff systems of ordinary differential equations (ODEs) with multiple time scales. These methods treat different partitions of the system with different timesteps for a more targeted and efficient solution compared to monolithic single rate approaches. With implicit methods used across all partitions, methods must find a balance between stability and the cost of solving nonlinear equations for the stages. In order to characterize this important trade-off, we explore multirate coupling strategies, problems for assessing linear stability, and techniques to efficiently implement Newton iterations for stage equations. Unlike much of the existing multirate stability analysis which is limited in scope to particular methods, we present general statements on stability and describe fundamental limitations for certain types of multirate schemes. New implicit multirate methods up to fourth order are derived, and their accuracy and efficiency properties are verified with numerical tests.


title: Contextual Text Denoising with Masked Language Models

abstract: Recently, with the help of deep learning models, significant advances have been made in different Natural Language Processing (NLP) tasks. Unfortunately, state-of-the-art models are vulnerable to noisy texts. We propose a new contextual text denoising algorithm based on the ready-to-use masked language model. The proposed algorithm does not require retraining of the model and can be integrated into any NLP system without additional training on paired cleaning training data. We evaluate our method under synthetic noise and natural noise and show that the proposed algorithm can use context information to correct noise text and improve the performance of noisy inputs in several downstream tasks.


title: Duality and Stability in Complex Multiagent State-Dependent Network Dynamics

abstract: Despite significant progress on stability analysis of conventional multiagent networked systems with weakly coupled state-network dynamics, most of the existing results have shortcomings in addressing multiagent systems with highly coupled state-network dynamics. Motivated by numerous applications of such dynamics, in our previous work [1], we initiated a new direction for stability analysis of such systems that uses a sequential optimization framework. Building upon that, in this paper, we extend our results by providing another angle on multiagent network dynamics from a duality perspective, which allows us to view the network structure as dual variables of a constrained nonlinear program. Leveraging that idea, we show that the evolution of the coupled state-network multiagent dynamics can be viewed as iterates of a primal-dual algorithm for a static constrained optimization/saddle-point problem. This view bridges the Lyapunov stability of state-dependent network dynamics and frequently used optimization techniques such as block coordinated descent, mirror descent, the Newton method, and the subgradient method. As a result, we develop a systematic framework for analyzing the Lyapunov stability of state-dependent network dynamics using techniques from nonlinear optimization. Finally, we support our theoretical results through numerical simulations from social science.


title: ALMA evidence for ram pressure compression and stripping of molecular gas in the Virgo cluster galaxy NGC 4402

abstract: High resolution (1'' $\times$ 2'') ALMA CO(2-1) observations of the ram pressure stripped galaxy NGC 4402 in the Virgo cluster show some of the clearest evidence yet for the impacts of ram pressure on the molecular ISM of a galaxy. The eastern side of the galaxy at $r \sim 4.5$ kpc, upon which ram pressure is incident, has a large (width $\sim$1 kpc, height $\sim$1 kpc above the disk midplane) extraplanar plume of molecular gas and dust. Molecular gas in the plume region shows distinct non-circular motions in the direction of the ram pressure; the kinematic offset of up to 60 km s$^{-1}$ is consistent with acceleration by ram pressure. We also detect a small amount of gas in clouds below the plume that are spatially and kinematically distinct from the surrounding medium, and appear to be decoupled from the stripped ISM. We propose that diffuse molecular gas is directly stripped but GMC density gas is not directly stripped, and so decouples from lower density stripped gas. However, GMCs become effectively stripped on short timescales. We also find morphological and kinematic signatures of ram pressure compression of molecular gas in a region of intense star formation on the leading side at $r \sim 3.5$ kpc. We propose that the compressed and stripped zones represent different evolutionary stages of the ram pressure interaction, and that feedback from star formation in the compressed zone facilitates the effective stripping of GMCs by making the gas cycle rapidly to a lower density diffuse state.


title: Statistics of TASEP with three merging characteristics

abstract: In this paper we consider the totally asymmetric simple exclusion process, with non-random initial condition having three regions of constant densities of particles. From left to right, the densities of the three regions are increasing. Consequently, there are three characteristics which meet, i.e. two shocks merge. We study the particle fluctuations at this merging point and show that they are given by a product of three (properly scaled) GOE Tracy-Widom distribution functions. We work directly in TASEP without relying on the connection to last passage percolation.


title: Building an Application Independent Natural Language Interface

abstract: Traditional approaches to building natural language (NL) interfaces typically use a semantic parser to parse the user command and convert it to a logical form, which is then translated to an executable action in an application. However, it is still challenging for a semantic parser to correctly parse natural language. For a different domain, the parser may need to be retrained or tuned, and a new translator also needs to be written to convert the logical forms to executable actions. In this work, we propose a novel and application independent approach to building NL interfaces that does not need a semantic parser or a translator. It is based on natural language to natural language matching and learning, where the representation of each action and each user command are both in natural language. To perform a user intended action, the system only needs to match the user command with the correct action representation, and then execute the corresponding action. The system also interactively learns new (paraphrased) commands for actions to expand the action representations over time. Our experimental results show the effectiveness of the proposed approach.


title: Twin Higgs with Exact $Z_2$

abstract: We present a novel mechanism for realistic electroweak symmetry breaking in Twin Higgs/neutral naturalness models where the $Z_2$ exchange symmetry can remain exactly unbroken. The exchange symmetry in the Yukawa sector will be implemented as an "N-trigonometric parity" $\sin N \frac{h}{f} \leftrightarrow \cos N \frac{h}{f}$. The Yukawa couplings will be suppressed leading to an N-suppressed Higgs quadratic term, without significantly affecting the quartic. We present a concrete implementation of this idea for general (odd) values of N using maximal symmetry, and a realistic benchmark model for $N=3$. We find that the tuning in the resulting Higgs potential is negligible, and also show that two-loop N-suppression violating gauge contributions can be sufficiently small. The Higgs potential and its couplings in top sector are different from other neutral naturalness models, which are the main predictions of our model and can be tested in colliders.


title: Higher-Form Symmetries, Bethe Vacua, and the 3d-3d Correspondence

abstract: By incorporating higher-form symmetries, we propose a refined definition of the theories obtained by compactification of the 6d $(2,0)$ theory on a three-manifold $M_3$. This generalization is applicable to both the 3d $\mathcal{N}=2$ and $\mathcal{N}=1$ supersymmetric reductions. An observable that is sensitive to the higher-form symmetries is the Witten index, which can be computed by counting solutions to a set of Bethe equations that are determined by $M_3$. This is carried out in detail for $M_3$ a Seifert manifold, where we compute a refined version of the Witten index. In the context of the 3d-3d correspondence, we complement this analysis in the dual topological theory, and determine the refined counting of flat connections on $M_3$, which matches the Witten index computation that takes the higher-form symmetries into account.


title: Towards Generalizable Neuro-Symbolic Systems for Commonsense Question Answering

abstract: Non-extractive commonsense QA remains a challenging AI task, as it requires systems to reason about, synthesize, and gather disparate pieces of information, in order to generate responses to queries. Recent approaches on such tasks show increased performance, only when models are either pre-trained with additional information or when domain-specific heuristics are used, without any special consideration regarding the knowledge resource type. In this paper, we perform a survey of recent commonsense QA methods and we provide a systematic analysis of popular knowledge resources and knowledge-integration methods, across benchmarks from multiple commonsense datasets. Our results and analysis show that attention-based injection seems to be a preferable choice for knowledge integration and that the degree of domain overlap, between knowledge bases and datasets, plays a crucial role in determining model success.


title: A Need for Dedicated Outreach Expertise and Online Programming: Astro2020 Science White Paper

abstract: Maximizing the public impact of astronomy projects in the next decade requires NSF-funded centers to support the development of online, mobile-friendly outreach and education activities. EPO teams with astronomy, education, and web development expertise should be in place to build accessible programs at scale and support astronomers doing outreach.


title: Weak inverse problem of calculus of variations for geodesic mappings and relation to harmonic maps

abstract: In this paper, we study the relation between geodesic and harmonic mappings. Harmonic mappings are defined between Riemannian manifolds as critical points of the energy functional, on the other hand, geodesic mappings are defined in a more general setting (manifolds with affine connections). Using the well-established formalism of calculus of variations on fibred manifolds we solve the weak inverse problem for the equation of geodesic mappings and get a variational equation which is a consequence of the geodesic mappings equation. For the connection on the target manifold we get the expected result, that it is a metric connection. However, we find that the connection on the source manifold need not be metric. The interesting result is that the metric which induces the connection on the target manifold can change between fibres and these changes are related to the connection on the source manifold. These results hint onto a possibility for a more general structure on the fibred manifold, than usually assumed.


title: Conditional Density Estimation, Latent Variable Discovery and Optimal Transport

abstract: A framework is proposed that addresses both conditional density estimation and latent variable discovery. The objective function maximizes explanation of variability in the data, achieved through the optimal transport barycenter generalized to a collection of conditional distributions indexed by a covariate --either given or latent-- in any suitable space. Theoretical results establish the existence of barycenters, a minimax formulation of optimal transport maps, and a general characterization of variability via the optimal transport cost. This framework leads to a family of non-parametric neural network-based algorithms, the BaryNet, with a supervised version that estimates conditional distributions and an unsupervised version that assigns latent variables. The efficacy of BaryNets is demonstrated by tests on both artificial and real-world data sets. A parallel drawn between autoencoders and the barycenter framework leads to the Barycentric autoencoder algorithm (BAE).


title: Measure rigidity of Anosov flows via the factorization method

abstract: Using the factorization method, a method pioneered by Eskin and Mirzakhani in their groundbreaking work about measure classification over the moduli space of translation surfaces, we show that generalized $u$-Gibbs states over quantitatively non-integrable Anosov systems are absolutely continuous with respect to the whole unstable manifold.


title: A long-lived stop with freeze-in and freeze-out dark matter in the hidden sector

abstract: In extended supersymmetric models with a hidden sector the lightest $R$-parity odd particle can reside in the hidden sector and act as dark matter. We consider the case when the hidden sector has ultraweak interactions with the visible sector. An interesting phenomenon arises if the LSP of the visible sector is charged in which case it will decay to the hidden sector dark matter. Due to the ultraweak interactions, the LSP of the visible sector will be long-lived decaying outside the detector after leaving a track inside. We investigate this possibility in the framework of a $U(1)_X$-extended MSSM/SUGRA model with a small gauge kinetic mixing and mass mixing between the $U(1)_X$ and $U(1)_Y$ where $U(1)_Y$ is the gauge group of the hypercharge. Specifically we investigate the case when the LSP of MSSM is a stop which decays into the hidden sector dark matter and has a lifetime long enough to traverse the LHC detector without decay. It is shown that such a particle can be detected at the HL-LHC and HE-LHC as an $R$-hadron which will look like a slow moving muon with a large transverse momentum $p_T$ and so can be detected by the track it leaves in the inner tracker and in the muon spectrometer. Further, due to the ultraweak couplings between the hidden sector and the MSSM fields, the dark matter particle has a relic density arising from a combination of the freeze-out and freeze-in mechanisms. It is found that even for the ultraweak or feeble interactions the freeze-out contribution relative to freeze-in contribution to the relic density is substantial to dominant, varying between 30\% to 74\% for the model points considered. It is subdominant to freeze-in for relatively small stop masses with relatively larger stop annihilation cross-sections and the dominant contribution to the relic density for relatively large stop masses and relatively smaller stop annihilation cross-sections.


title: Superconvergence of differential structure for finite element methods on perturbed surface meshes

abstract: Superconvergence of differential structure on discretized surfaces is studied in this paper. The newly introduced geometric supercloseness provides us with a fundamental tool to prove the superconvergence of gradient recovery on deviated surfaces. An algorithmic framework for gradient recovery without exact geometric information is introduced. Several numerical examples are documented to validate the theoretical results.


title: Indirect methods in nuclear astrophysics with relativistic radioactive beams

abstract: Reactions with radioactive nuclear beams at relativistic energies have opened new doors to clarify the mechanisms of stellar evolution and cataclysmic events involving stars and during the big bang epoch. Numerous nuclear reactions of astrophysical interest cannot be assessed directly in laboratory experiments. Ironically, some of the information needed to describe such reactions, at extremely low energies (e.g., keVs), can only be studied on Earth by using relativistic collisions between heavy ions at GeV energies. In this contribution, we make a short review of experiments with relativistic radioactive beams and of the theoretical methods needed to understand the physics of stars, adding to the knowledge inferred from astronomical observations. We continue by introducing a more detailed description of how the use of relativistic radioactive beams can help to solve astrophysical puzzles and several successful experimental methods. State-of-the-art theories are discussed at some length with the purpose of helping us understand the experimental results reported. The review is not complete and we have focused most of it to traditional methods aiming at the determination of the equation of state of symmetric and asymmetric nuclear matter and the role of the symmetry energy. Whenever possible, under the limitations of our present understanding of experimental data and theory, we try to pinpoint the information still missing to further understand how stars evolve, explode, and how their internal structure might be.


title: Explainable Prediction of Adverse Outcomes Using Clinical Notes

abstract: Clinical notes contain a large amount of clinically valuable information that is ignored in many clinical decision support systems due to the difficulty that comes with mining that information. Recent work has found success leveraging deep learning models for the prediction of clinical outcomes using clinical notes. However, these models fail to provide clinically relevant and interpretable information that clinicians can utilize for informed clinical care. In this work, we augment a popular convolutional model with an attention mechanism and apply it to unstructured clinical notes for the prediction of ICU readmission and mortality. We find that the addition of the attention mechanism leads to competitive performance while allowing for the straightforward interpretation of predictions. We develop clear visualizations to present important spans of text for both individual predictions and high-risk cohorts. We then conduct a qualitative analysis and demonstrate that our model is consistently attending to clinically meaningful portions of the narrative for all of the outcomes that we explore.


title: Robust and Computationally-Efficient Anomaly Detection using Powers-of-Two Networks

abstract: Robust and computationally efficient anomaly detection in videos is a problem in video surveillance systems. We propose a technique to increase robustness and reduce computational complexity in a Convolutional Neural Network (CNN) based anomaly detector that utilizes the optical flow information of video data. We reduce the complexity of the network by denoising the intermediate layer outputs of the CNN and by using powers-of-two weights, which replaces the computationally expensive multiplication operations with bit-shift operations. Denoising operation during inference forces small valued intermediate layer outputs to zero. The number of zeros in the network significantly increases as a result of denoising, we can implement the CNN about 10% faster than a comparable network while detecting all the anomalies in the testing set. It turns out that denoising operation also provides robustness because the contribution of small intermediate values to the final result is negligible. During training we also generate motion vector images by a Generative Adversarial Network (GAN) to improve the robustness of the overall system. We experimentally observe that the resulting system is robust to background motion.


title: On the $\zeta_3$-Pell equation

abstract: Let $K = \mathbb{Q}(\zeta_3)$, where $\zeta_3$ is a primitive root of unity. In this paper we study the distribution of integers $\alpha \in \mathcal{O}_K$ for which the norm equation $N_{K(\sqrt[3]{\alpha})/K}(\mathbf{x}) = \zeta_3$ is solvable for integers $\mathbf{x} \in \mathcal{O}_{K(\sqrt[3]{\alpha})}$. The analogous question for $\zeta_2 = -1$ is the well-known negative Pell equation. We also address the natural generalization of Stevenhagen's conjecture on the negative Pell equation in this setting.


title: Comparison of Different Spike Sorting Subtechniques Based on Rat Brain Basolateral Amygdala Neuronal Activity

abstract: Developing electrophysiological recordings of brain neuronal activity and their analysis provide a basis for exploring the structure of brain function and nervous system investigation. The recorded signals are typically a combination of spikes and noise. High amounts of background noise and possibility of electric signaling recording from several neurons adjacent to the recording site have led scientists to develop neuronal signal processing tools such as spike sorting to facilitate brain data analysis. Spike sorting plays a pivotal role in understanding the electrophysiological activity of neuronal networks. This process prepares recorded data for interpretations of neurons interactions and understanding the overall structure of brain functions. Spike sorting consists of three steps: spike detection, feature extraction, and spike clustering. There are several methods to implement each of spike sorting steps. This paper provides a systematic comparison of various spike sorting sub-techniques applied to real extracellularly recorded data from a rat brain basolateral amygdala. An efficient sorted data resulted from careful choice of spike sorting sub-methods leads to better interpretation of the brain structures connectivity under different conditions, which is a very sensitive concept in diagnosis and treatment of neurological disorders. Here, spike detection is performed by appropriate choice of threshold level via three different approaches. Feature extraction is done through PCA and Kernel PCA methods, which Kernel PCA outperforms. We have applied four different algorithms for spike clustering including K-means, Fuzzy C-means, Bayesian and Fuzzy maximum likelihood estimation. As one requirement of most clustering algorithms, optimal number of clusters is achieved through validity indices for each method. Finally, the sorting results are evaluated using inter-spike interval histograms.


title: Notes on the complex Sachdev-Ye-Kitaev model

abstract: We describe numerous properties of the Sachdev-Ye-Kitaev model for complex fermions with $N\gg 1$ flavors and a global U(1) charge. We provide a general definition of the charge in the $(G,\Sigma)$ formalism, and compute its universal relation to the infrared asymmetry of the Green function. The same relation is obtained by a renormalization theory. The conserved charge contributes a compact scalar field to the effective action, from which we derive the many-body density of states and extract the charge compressibility. We compute the latter via three distinct numerical methods and obtain consistent results. Finally, we present a two dimensional bulk picture with free Dirac fermions for the zero temperature entropy.


title: Invariant-mass spectroscopy of $^{14}$O excited states

abstract: Excited states in $^{14}$O have been investigated both experimentally and theoretically. Experimentally, these states were produced via neutron-knockout reactions with a fast $^{15}$O beam and the invariant-mass technique was employed to isolate the 1$p$ and 2$p$ decay channels and determine their branching ratios. The spectrum of excited states was also calculated with the Shell Model Embedded in the Continuum that treats bound and scattering states in a unified model. By comparing energies, widths and decay branching patterns, spin and parity assignments for all experimentally observed levels below 8 MeV are made. This includes the location of the second 2$^{+}$ state that we find is in near degeneracy with the third 0$^{+}$ state. An interesting case of sequential 2$p$ decay through a pair of degenerate $^{13}$N excited states with opposite parities was found where the interference between the two sequential decay pathways produces an unusual relative-angle distribution between the emitted protons.


title: Bayesian inference for high-dimensional nonstationary Gaussian processes

abstract: In spite of the diverse literature on nonstationary spatial modeling and approximate Gaussian process (GP) methods, there are no general approaches for conducting fully Bayesian inference for moderately sized nonstationary spatial data sets on a personal laptop. For statisticians and data scientists who wish to learn about spatially-referenced data and conduct posterior inference and prediction with appropriate uncertainty quantification, the lack of such approaches and corresponding software is a significant limitation. In this paper, we develop methodology for implementing formal Bayesian inference for a general class of nonstationary GPs. Our novel approach uses pre-existing frameworks for characterizing nonstationarity in a new way that is applicable for small to moderately sized data sets via modern GP likelihood approximations. Posterior sampling is implemented using flexible MCMC methods, with nonstationary posterior prediction conducted as a post-processing step. We demonstrate our novel methods on two data sets, ranging from several hundred to several thousand locations, and compare our methodology with related statistical methods that provide off-the-shelf software. All of our methods are implemented in the freely available BayesNSGP software package for R.


title: Continuity of Lyapunov exponents for non-uniformly fiber-bunched cocycles

abstract: We provide conditions which imply the continuity of the Lyapunov exponents for non-uniformly fiber-bunched cocycles in $SL(2,\mathbb{R})$. The main theorem is an extension of the result of Backes, Brown and Butler and gives a partial answer to a conjecture of Marcelo Viana.


title: CALC2.0: Combining Appearance, Semantic and Geometric Information for Robust and Efficient Visual Loop Closure

abstract: Traditional attempts for loop closure detection typically use hand-crafted features, relying on geometric and visual information only, whereas more modern approaches tend to use semantic, appearance or geometric features extracted from deep convolutional neural networks (CNNs). While these approaches are successful in many applications, they do not utilize all of the information that a monocular image provides, and many of them, particularly the deep-learning based methods, require user-chosen thresholding to actually close loops -- which may impact generality in practical applications. In this work, we address these issues by extracting all three modes of information from a custom deep CNN trained specifically for the task of place recognition. Our network is built upon a combination of a semantic segmentator, Variational Autoencoder (VAE) and triplet embedding network. The network is trained to construct a global feature space to describe both the visual appearance and semantic layout of an image. Then local keypoints are extracted from maximally-activated regions of low-level convolutional feature maps, and keypoint descriptors are extracted from these feature maps in a novel way that incorporates ideas from successful hand-crafted features. These keypoints are matched globally for loop closure candidates, and then used as a final geometric check to refute false positives. As a result, the proposed loop closure detection system requires no touchy thresholding, and is highly robust to false positives -- achieving better precision-recall curves than the state-of-the-art NetVLAD, and with real-time speeds.


title: End-to-end Microphone Permutation and Number Invariant Multi-channel Speech Separation

abstract: An important problem in ad-hoc microphone speech separation is how to guarantee the robustness of a system with respect to the locations and numbers of microphones. The former requires the system to be invariant to different indexing of the microphones with the same locations, while the latter requires the system to be able to process inputs with varying dimensions. Conventional optimization-based beamforming techniques satisfy these requirements by definition, while for deep learning-based end-to-end systems those constraints are not fully addressed. In this paper, we propose transform-average-concatenate (TAC), a simple design paradigm for channel permutation and number invariant multi-channel speech separation. Based on the filter-and-sum network (FaSNet), a recently proposed end-to-end time-domain beamforming system, we show how TAC significantly improves the separation performance across various numbers of microphones in noisy reverberant separation tasks with ad-hoc arrays. Moreover, we show that TAC also significantly improves the separation performance with fixed geometry array configuration, further proving the effectiveness of the proposed paradigm in the general problem of multi-microphone speech separation.


title: Use of R-trees to improve reconstruction time in pixel trackers

abstract: Computing time is becoming a key issue for tracking algorithms both online and off-line. Programming using adequate data structures can largely improve the efficiency of the reconstruction in terms of time response. We propose using one such data structure, called R-tree, that performs a fast, flexible and custom spatial indexing of the hits based on a neighbourhood organisation. The overhead required to prepare the data structure shows to be largely compensated by the efficiency in the search of hits that are candidate to belong to the same track when events present a large number of hits. The study, including different indexing approaches, is performed for a generic pixel tracker largely inspired in the upgrade of the LHCb vertex locator with a backwards reconstruction algorithm of the cellular automaton type.


title: Sample Complexity of Learning Mixtures of Sparse Linear Regressions

abstract: In the problem of learning mixtures of linear regressions, the goal is to learn a collection of signal vectors from a sequence of (possibly noisy) linear measurements, where each measurement is evaluated on an unknown signal drawn uniformly from this collection. This setting is quite expressive and has been studied both in terms of practical applications and for the sake of establishing theoretical guarantees. In this paper, we consider the case where the signal vectors are sparse; this generalizes the popular compressed sensing paradigm. We improve upon the state-of-the-art results as follows: In the noisy case, we resolve an open question of Yin et al. (IEEE Transactions on Information Theory, 2019) by showing how to handle collections of more than two vectors and present the first robust reconstruction algorithm, i.e., if the signals are not perfectly sparse, we still learn a good sparse approximation of the signals. In the noiseless case, as well as in the noisy case, we show how to circumvent the need for a restrictive assumption required in the previous work. Our techniques are quite different from those in the previous work: for the noiseless case, we rely on a property of sparse polynomials and for the noisy case, we provide new connections to learning Gaussian mixtures and use ideas from the theory of error-correcting codes.


title: Investigating Resistance of Deep Learning-based IDS against Adversaries using min-max Optimization

abstract: With the growth of adversarial attacks against machine learning models, several concerns have emerged about potential vulnerabilities in designing deep neural network-based intrusion detection systems (IDS). In this paper, we study the resilience of deep learning-based intrusion detection systems against adversarial attacks. We apply the min-max (or saddle-point) approach to train intrusion detection systems against adversarial attack samples in NSW-NB 15 dataset. We have the max approach for generating adversarial samples that achieves maximum loss and attack deep neural networks. On the other side, we utilize the existing min approach [2] [9] as a defense strategy to optimize intrusion detection systems that minimize the loss of the incorporated adversarial samples during the adversarial training. We study and measure the effectiveness of the adversarial attack methods as well as the resistance of the adversarially trained models against such attacks. We find that the adversarial attack methods that were designed in binary domains can be used in continuous domains and exhibit different misclassification levels. We finally show that principal component analysis (PCA) based feature reduction can boost the robustness in intrusion detection system (IDS) using a deep neural network (DNN).


title: Asteroseismic ages of red-giant stars from grid-based modelling: the impact of systematics in effective temperature or metallicity

abstract: An increasingly popular method to determine stellar ages of red-giant stars for the purpose of Galactic archaeology is asteroseismic grid-based modelling (GBM). In asteroseismic GBM of red-giant stars with solar-like oscillations the large frequency separation ($\Delta\nu$) and the frequency of maximum oscillation power ($\nu_{\rm max}$) are commonly used asteroseismic observables, in addition to the usual spectroscopic parameters effective temperature ($T_{\rm eff}$) and metallicity ([Fe/H]). The precision with which $\Delta\nu$ and $\nu_{\rm max}$ can be determined largely depends on the length of the timeseries data (assuming the stars are bright enough that the oscillations can be detected). The question that we aim to answer here is: with what precision should [Fe/H] and $T_{\rm eff}$ be obtained to derive stellar ages of red-giant stars through asteroseismic GBM given the precision of $\Delta\nu$ and $\nu_{\rm max}$ that we can expect from the TESS data? In most cases presented here the value of $\sigma_{\rm age}$ shows only weak, or no, dependence on the uncertainties and biases in [Fe/H] and $T_{\rm eff}$ ingested in this study. Only in case a large ($>100$~K) $T_{\rm eff}$ biases is accompanied by a small (50 K) uncertainties in $T_{\rm eff}$ on the red-giant branch, $\Delta_{\rm age}$ is of the order of or larger than $\sigma_{\rm age}$ and in these cases we expect an impact of this bias on the resulting age estimate.


title: Evaluating a Semi-Autonomous Brain-Computer Interface Based on Conformal Geometric Algebra and Artificial Vision

abstract: In this paper, we evaluate a semi-autonomous brain-computer interface (BCI) for manipulation tasks. In such system, the user controls a robotic arm through motor imagery commands. In traditional process-control BCI systems, the user has to provide those commands continuously in order manipulate the effector of the robot step-by-step, which results in a tiresome process for simple tasks such as pick and replace an item from a surface. Here, we take a semi-autonomous approach based on a conformal geometric algebra model that solves the inverse kinematics of the robot on the fly, then the user only has to decide on the start of the movement and the final position of the effector (goal-selection approach). Under these conditions, we implemented pick-and-place tasks with a disk as an item and two target areas placed on the table at arbitrary positions. An artificial vision (AV) algorithm was used to obtain the positions of the items expressed in the robot frame through images captured with a webcam. Then, the AV algorithm is integrated to the inverse kinematics model to perform the manipulation tasks. As proof-of-concept, different users were trained to control the pick-and-place tasks through the process-control and semi-autonomous goal-selection approaches, so that the performance of both schemes could be compared. Our results show the superiority in performance of the semi-autonomous approach, as well as evidence of less mental fatigue with it.


title: Spatially Coupled Generalized LDPC Codes: Asymptotic Analysis and Finite Length Scaling

abstract: Generalized low-density parity-check (GLDPC) codes are a class of LDPC codes in which the standard single parity check (SPC) constraints are replaced by constraints defined by a linear block code. These stronger constraints typically result in improved error floor performance, due to better minimum distance and trapping set properties, at a cost of some increased decoding complexity. In this paper, we study spatially coupled generalized low-density parity-check (SC-GLDPC) codes and present a comprehensive analysis of these codes, including: (1) an iterative decoding threshold analysis of SC-GLDPC code ensembles demonstrating capacity approaching thresholds via the threshold saturation effect; (2) an asymptotic analysis of the minimum distance and free distance properties of SC-GLDPC code ensembles, demonstrating that the ensembles are asymptotically good; and (3) an analysis of the finite-length scaling behavior of both GLDPC block codes and SC-GLDPC codes based on a peeling decoder (PD) operating on a binary erasure channel (BEC). Results are compared to GLDPC block codes, and the advantages and disadvantages of SC-GLDPC codes are discussed.


title: Phase Field Benchmark Problems for Dendritic Growth and Linear Elasticity

abstract: We present the second set of benchmark problems for phase field models that are being jointly developed by the Center for Hierarchical Materials Design (CHiMaD) and the National Institute of Standards and Technology (NIST) along with input from other members in the phase field community. As the integrated computational materials engineering (ICME) approach to materials design has gained traction, there is an increasing need for quantitative phase field results. New algorithms and numerical implementations increase computational capabilities, necessitating standard problems to evaluate their impact on simulated microstructure evolution as well as their computational performance. We propose one benchmark problem for solidification and dendritic growth in a single-component system, and one problem for linear elasticity via the shape evolution of an elastically constrained precipitate. We demonstrate the utility and sensitivity of the benchmark problems by comparing the results of 1) dendritic growth simulations performed with different time integrators and 2) elastically constrained precipitate simulations with different precipitate sizes, initial conditions, and elastic moduli. These numerical benchmark problems will provide a consistent basis for evaluating different algorithms, both existing and those to be developed in the future, for accuracy and computational efficiency when applied to simulate physics often incorporated in phase field models.


title: Alexa, Who Am I Speaking To? Understanding Users' Ability to Identify Third-Party Apps on Amazon Alexa

abstract: Many Internet of Things (IoT) devices have voice user interfaces (VUIs). One of the most popular VUIs is Amazon's Alexa, which supports more than 47,000 third-party applications ("skills"). We study how Alexa's integration of these skills may confuse users. Our survey of 237 participants found that users do not understand that skills are often operated by third parties, that they often confuse third-party skills with native Alexa functions, and that they are unaware of the functions that the native Alexa system supports. Surprisingly, users who interact with Alexa more frequently are more likely to conclude that a third-party skill is native Alexa functionality. The potential for misunderstanding creates new security and privacy risks: attackers can develop third-party skills that operate without users' knowledge or masquerade as native Alexa functions. To mitigate this threat, we make design recommendations to help users distinguish native and third-party skills.


title: Linear Response Based Parameter Estimation in the Presence of Model Error

abstract: Recently, we proposed a method to estimate parameters of stochastic dynamics based on the linear response statistics. The method rests upon a nonlinear least-squares problem that takes into account the response properties that stem from the Fluctuation-Dissipation Theory. In this article, we address an important issue that arises in the presence of model error. In particular, when the equilibrium density function is high dimensional and non-Gaussian, and in some cases, is unknown, the linear response statistics are inaccessible. We show that this issue can be resolved by fitting the imperfect model to appropriate marginal linear response statistics that can be approximated using the available data and parametric or nonparametric models. The effectiveness of the parameter estimation approach is demonstrated in the context of molecular dynamical models (Langevin dynamics) with a non-uniform temperature profile, where the modeling error is due to coarse-graining, and a PDE (non-Langevin dynamics) that exhibits spatiotemporal chaos, where the model error arises from a severe spectral truncation. In these examples, we show how the imperfect models, the Langevin equation with parameters estimated using the proposed scheme, can predict the nonlinear response statistics of the underlying dynamics under admissible external disturbances.


title: Finslerian geometrization of quantum mechanics in the hydrodynamical representation

abstract: We consider a Finslerian type geometrization of the non-relativistic quantum mechanics in its hydrodynamical (Madelung) formulation, by also taking into account the effects of the presence of the electromagnetic fields on the particle motion. In the Madelung representation the Schr\"{o}dinger equation can be reformulated as the classical continuity and Euler equations of classical fluid mechanics in the presence of a quantum potential, representing the quantum hydrodynamical evolution equations. The equation of particle motion can then be obtained from a Lagrangian similar to its classical counterpart. After the reparametrization of the Lagrangian it turns out that the Finsler metric describing the geometric properties of quantum hydrodynamics is a Kropina metric. We present and discuss in detail the metric and the geodesic equations describing the geometric properties of the quantum motion in the presence of electromagnetic fields. As an application of the obtained formalism we consider the Zermelo navigation problem in a quantum hydrodynamical system, whose solution is given by a Kropina metric. The case of the Finsler geometrization of the quantum hydrodynamical motion of spinless particles in the absence of electromagnetic interactions is considered in detail, and the Zermelo navigation problem for this case is also discussed.