title: Evidence for dominant phonon-electron scattering in Weyl semimetal WP$_{2}$

abstract: Topological semimetals have revealed a wide array of novel transport phenomena, including electron hydrodynamics, quantum field theoretic anomalies, and extreme magnetoresistances and mobilities. However, the scattering mechanisms central to these behaviors remain largely unexplored. Here we reveal signatures of significant phonon-electron scattering in the type-II Weyl semimetal WP$_{2}$ via temperature-dependent Raman spectroscopy. Over a large temperature range, we find that the decay rates of the lowest energy $A_{1}$ modes are dominated by phonon-electron rather than phonon-phonon scattering. In conjunction with first-principles calculations, a combined analysis of the momentum, energy, and symmetry-allowed decay paths indicates this results from intraband scattering of the electrons. The excellent agreement with theory further suggests that such results could be true for the acoustic modes. We thus provide evidence for the importance of phonons in the transport properties of topological semimetals and identify specific properties that may contribute to such behavior in other materials.

categories: cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.other cond-mat.str-el


title: Cosmic Ballet III: halo spin evolution in the cosmic web

abstract: We explore the evolution of halo spins in the cosmic web using a very large sample of dark matter haloes in the $\Lambda$CDM Planck-Millennium N-body simulation. We use the NEXUS+ multiscale formalism to identify the hierarchy of filaments and sheets of the cosmic web at several redshifts. We find that at all times the magnitude of halo spins correlates with the web environment, being largest in filaments, and, for the first time, we show that it also correlates with filament thickness as well as the angle between spin-orientation and the spine of the host filament. For example, massive haloes in thick filaments spin faster than their counterparts in thin filaments, while for low-mass haloes the reverse is true. We also have studied the evolution of alignment between halo spin orientations and the preferential axes of filaments and sheets. The alignment varies with halo mass, with the spins of low-mass haloes being predominantly along the filament spine, while those of high-mass haloes being predominantly perpendicular to the filament spine. On average, for all halo masses, halo spins become more perpendicular to the filament spine at later times. At all redshifts, the spin alignment shows a considerable variation with filament thickness, with the halo mass corresponding to the transition from parallel to perpendicular alignment varying by more than one order of magnitude. The environmental dependence of halo spin magnitude shows little evolution for $z\leq2$ and is likely a consequence of the correlations in the initial conditions or high redshift effects


title: The physical origins of low-mass spin bias

abstract: At $z=0$, higher-spin haloes with masses above $\log(\text{M}_{\text{c}}/h^{-1}\text{M}_\odot)\simeq 11.5$ have a higher bias than lower-spin haloes of the same mass. However, this trend is known to invert below this characteristic crossover mass, $\text{M}_{\text{c}}$. In this paper, we measure the redshift evolution and scale dependence of halo spin bias at the low-mass end and demonstrate that the inversion of the signal is entirely produced by the effect of splashback haloes. These low-mass haloes tend to live in the vicinity of significantly more massive haloes, thus sharing their large-scale bias properties. We further show that the location of the redshift-dependent crossover mass scale $\text{M}_{\text{c}}(z)$ is completely determined by the relative abundance of splashbacks in the low- and high-spin subpopulations. Once splashback haloes are removed from the sample, the intrinsic mass dependence of spin bias is recovered. Since splashbacks have been shown to account for some of the assembly bias signal at the low-mass end, our results unveil a specific link between two different secondary bias trends: spin bias and assembly bias.


title: QNEC2 in deformed holographic CFTs

abstract: We use the quantum null energy condition in strongly coupled two-dimensional field theories (QNEC2) as diagnostic tool to study a variety of phase structures, including crossover, second and first order phase transitions. We find a universal QNEC2 constraint for first order phase transitions with kinked entanglement entropy and discuss in general the relation between the QNEC2-inequality and monotonicity of the Casini-Huerta c-function. We then focus on a specific example, the holographic dual of which is modelled by three-dimensional Einstein gravity plus a massive scalar field with one free parameter in the self-interaction potential. We study translation invariant stationary states dual to domain walls and black branes. Depending on the value of the free parameter we find crossover, second and first order phase transitions between such states, and the c-function either flows to zero or to a finite value in the infrared. Strikingly, evaluating QNEC2 for ground state solutions allows to predict the existence of phase transitions at finite temperature.


title: Extended X-ray emission from the z=4.26 radio galaxy 4C 63.20

abstract: We report on deep Chandra X-ray Telescope imaging observations of 4C 63.20, one of the few known radio galaxies at z>3.5. The X-ray counterpart is resolved into a core plus two off-nuclear sources that (combined) account for close to 30% of the total X-ray flux. Their morphology and orientation are consistent with a diffuse, lobe-like nature, albeit compact hotspots cannot be ruled out. The broadband spectral energy distribution of 4C 63.20 can be reproduced with a jet model where the majority of the radio flux can be ascribed to synchrotron emission from the hotspots, whereas the (non-nuclear) X-ray emission is produced via Inverse Compton (IC) off of Cosmic Microwave Background (CMB) photons within the extended lobes. This scenario is broadly consistent with the expectation from highly magnetized lobes in a hotter CMB, and supports the view that IC/CMB may quench less extreme radio lobes at high redshifts.


title: Occurrence of Tachyonic Preheating in the Mixed Higgs-$R^2$ Model

abstract: It has recently been suggested that at the post-inflationary stage of the mixed Higgs-$R^2$ model of inflation efficient particle production can arise from the tachyonic instability of the Higgs field. It might complete the preheating of the Universe if appropriate conditions are satisfied, especially in the Higgs-like regime. In this paper, we study this behavior in more depth, including the conditions for occurrence, analytical estimates for the maximal efficiency, and the necessary degree of fine-tuning among the model parameters to complete preheating by this effect. We find that the parameter sets that cause the most efficient tachyonic instabilities obey simple laws in both the Higgs-like regime and the $R^2$-like regime, respectively. We then estimate the efficiency of this instability. In particular, even in the deep $R^2$-like regime with a small non-minimal coupling, this effect is strong enough to complete preheating although a severe fine-tuning is required among the model parameters. We also estimate how much fine-tuning is needed to complete preheating by this effect. It is shown that the fine-tuning of parameters for the sufficient particle production is at least $ < \mathcal{O}(0.1) $ in the deep Higgs-like regime with a large scalaron mass, while it is more severe $\sim {\cal O}(10^{-4})-{\cal O}(10^{-5})$ in the $R^2$-like regime with a small non-minimal coupling.


title: Equivalence between vortices, twists and chiral gauge fields in Kitaev's honeycomb lattice model

abstract: We demonstrate that $\mathbb{Z}_2$ gauge transformations and lattice deformations in Kitaev's honeycomb lattice model can have the same description in the continuum limit of the model in terms of chiral gauge fields. The chiral gauge fields are coupled to the Majorana fermions that satisfy the Dirac dispersion relation in the non-Abelian sector of the model. For particular values, the effective chiral gauge field becomes equivalent to the $\mathbb{Z}_2$ gauge field, enabling us to associate effective fluxes to lattice deformations. Motivated by this equivalence, we consider Majorana-bounding $\pi$ vortices and Majorana-bounding lattice twists and demonstrate that they are adiabatically connected to each other. This equivalence opens the possibility for novel encoding of Majorana-bounding defects that might be easier to realise in experiments.


title: Modeling the diffusive dynamics of critical fluctuations near the QCD critical point

abstract: The experimental search for the QCD critical point by means of relativistic heavy-ion collisions necessitates the development of dynamical models of fluctuations. In this work we study the fluctuations of the net-baryon density near the critical point. Due to net-baryon number conservation the correct dynamics is given by the fluid dynamical diffusion equation, which we extend by a white noise stochastic term to include intrinsic fluctuations. We quantify finite resolution and finite size effects by comparing our numerical results to analytic expectations for the structure factor and the equal-time correlation function. In small systems the net-baryon number conservation turns out to be quantitatively and qualitatively important, as it introduces anticorrelations at larger distances. Including nonlinear coupling terms in the form of a Ginzburg-Landau free energy functional we observe non-Gaussian fluctuations quantified by the excess kurtosis. We study the dynamical properties of the system close to equilibrium, for a sudden quench in temperature and a Hubble-like temperature evolution. In the real-time dynamical systems we find the important dynamical effects of critical slowing down, weakening of the extremal value and retardation of the fluctuation signal. In this work we establish a set of general tests, which should be met by any model propagating fluctuations, including upcoming $3+1$ dimensional fluctuating fluid dynamics.


title: The distant, galaxy cluster environment of the short GRB 161104A at $z\sim 0.8$ and a comparison to the short GRB host population

abstract: We present optical observations of the Swift short-duration gamma-ray burst (GRB) GRB 161104A and its host galaxy at $z=0.793 \pm 0.003$. We model the multiband photometry and spectroscopy with the stellar population inference code Prospector, and explore the posterior using nested sampling. We find that the mass-weighted age $t_m = 2.12^{+0.23}_{-0.21}$~Gyr, stellar mass $\log{(M/M_\odot)} = 10.21 \pm 0.04$, metallicity $\log{(Z/Z_\odot)} = 0.08^{+0.05}_{-0.06}$, dust extinction $A_V = 0.08^{+0.08}_{-0.05}$ mag, and the star formation rate $\text{SFR} = 9.9 \times 10^{-2} M_\odot$~yr$^{-1}$. These properties, along with a prominent 4000 Angstrom break and optical absorption lines classify this host as an early-type, quiescent galaxy. Using Dark Energy Survey galaxy catalogues, we demonstrate that the host of GRB 161104A resides on the outskirts of a galaxy cluster at $z\approx 0.8$, situated $\approx 1$ Mpc from the likely brightest cluster galaxy. We also present new modeling for 20 additional short GRB hosts ($\approx33\%$ of which are early-type galaxies), finding population medians of $\log(M/M_\odot) = 9.94^{+0.88}_{-0.98}$ and $t_m = 1.07^{+1.98}_{-0.67}$~Gyr ($68\%$ confidence). We further find that the host of GRB 161104A is more distant, less massive, and younger than the four other short GRB hosts known to be associated with galaxy clusters. Cluster short GRBs have faint afterglows, in the lower $\approx 11\%$ ($\approx 30\%$) of observed X-ray (optical) luminosities. We place a lower limit on the fraction of short GRBs in galaxy clusters versus those in the field of $\approx 5-13\%$, consistent with the fraction of stellar mass $\approx 10-20\%$ in galaxy clusters at redshifts $0.1 \leq z \leq 0.8$.


title: Hall diffusion anomaly and transverse Einstein relation

abstract: It is commonly believed that the current response of an electron fluid to a mechanical force (such as an electric field) or to a ``statistical force" (e.g., a gradient of chemical potential) are governed by a single linear transport coefficient - the electric conductivity. We argue that this is not the case in anomalous Hall materials. In particular, we find that transverse (Hall) currents manifest two distinct Hall responses governed by an unconventional {\it transverse} Einstein relation that captures an anomalous relation between the Hall conductivity and the Hall diffusion constant. We give examples of when the Hall diffusion anomaly is prominent, resulting in situations where the transverse diffusion process overwhelms the Hall conductivity and vice versa.


title: Evidence from APOGEE for the presence of a major building block of the halo buried in the inner Galaxy

abstract: We report evidence from APOGEE for the presence of a new metal-poor stellar structure located within $\sim$4~kpc of the Galactic centre. Characterised by a chemical composition resembling those of low mass satellites of the Milky Way, this new inner Galaxy structure (IGS) seems to be chemically and dynamically detached from more metal-rich populations in the inner Galaxy. We conjecture that this structure is associated with an accretion event that likely occurred in the early life of the Milky Way. Comparing the mean elemental abundances of this structure with predictions from cosmological numerical simulations, we estimate that the progenitor system had a stellar mass of $\sim5\times10^8M_\odot$, or approximately twice the mass of the recently discovered Gaia-Enceladus/Sausage system. We find that the accreted:{\it in situ} ratio within our metal-poor ([Fe/H]$<$--0.8) bulge sample is somewhere between 1:3 and 1:2, confirming predictions of cosmological numerical simulations by various groups.


title: A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation

abstract: A precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schr\"odinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from $5 \times 10^{13}\:\rm W/cm^2$ to $5 \times 10^{14}\:\rm W/cm^2$. We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.


title: The Entropy of Ricci Flows with Type-I Scalar Curvature Bounds

abstract: In this paper, we extend the theory of Ricci flows satisfying a Type-I scalar curvature condition at a finite-time singularity. In [Bam16], Bamler showed that a Type-I rescaling procedure will produce a singular shrinking gradient Ricci soliton with singularities of codimension 4. We prove that the entropy of a conjugate heat kernel based at the singular time converges to the soliton entropy of the singular soliton, and use this to characterize the singular set of the Ricci flow solution in terms of a heat kernel density function. This generalizes results previously only known with the stronger assumption of a Type-I curvature bound. We also show that in dimension 4, the singular Ricci soliton is smooth away from finitely many points, which are conical smooth orbifold singularities


title: The Higgs field and early universe cosmology: a (brief) review

abstract: We review and discuss recent work exploring the implications of the Higgs field for early universe cosmology, and vice versa. Depending on the model under consideration the Higgs may be one of a few scalar fields determining the evolution and fate of the Universe, or the Higgs field may be connected to a rich sector of scalar moduli with complicated dynamics. In particular we look at the potential consequences of the Higgs field for inflation and its predictions, for the (meta)stability of the Standard Model vacuum, and for the existence of dynamical selection mechanisms in the landscape.


title: Dark Matter Particle in QCD

abstract: We report on the possibility that the Dark Matter particle is a stable, neutral, as-yet-undiscovered hadron in the standard model. The existence of a compact color-flavor-spin singlet sexaquark (S, uuddss) with mass ~2m_p, is compatible with current knowledge. The S interacts with baryons primarily via a Yukawa interaction of coupling strength alpha_SN, mediated by omega and phi vector mesons having mass ~1 GeV. If it exists, the S is a very attractive DM candidate. The relic abundance of S Dark Matter (SDM) is established when the Universe transitions from the quark-gluon plasma to the hadronic phase at ~150 MeV and is in remarkable agreement with the observed Omega_DM/Omega_b = 5.3+-0.1; this is a no-free-parameters result because the relevant parameters are known from QCD. Survival of this relic abundance to low temperature requires the breakup amplitude gtilde <~ 2 10^-6, comfortably compatible with theory expectations and observational bounds because the breakup amplitude is dynamically suppressed and many orders of magnitude smaller, as we show. The scattering cross section can differ by orders of magnitude from Born approximation, depending on alpha_SN, requiring reanalysis of observational limits. We use direct detection experiments and cosmological constraints to determine the allowed region of alpha_SN. For a range of allowed values, we predict exotic nuclear isotopes at a detectable level with mass offset ~2 amu. The most promising approaches for detecting the sexaquark in accelerator experiments are to search for a long-interaction-length neutral particle component in the central region of relativistic heavy ion collisions or using a beam-dump setup, and to search for evidence of missing particle production characterized by unbalanced baryon number and strangeness using Belle-II or possibly GLUEX at J-Lab.


title: Generative Hierarchical Features from Synthesizing Images

abstract: Generative Adversarial Networks (GANs) have recently advanced image synthesis by learning the underlying distribution of the observed data. However, how the features learned from solving the task of image generation are applicable to other vision tasks remains seldom explored. In this work, we show that learning to synthesize images can bring remarkable hierarchical visual features that are generalizable across a wide range of applications. Specifically, we consider the pre-trained StyleGAN generator as a learned loss function and utilize its layer-wise representation to train a novel hierarchical encoder. The visual feature produced by our encoder, termed as Generative Hierarchical Feature (GH-Feat), has strong transferability to both generative and discriminative tasks, including image editing, image harmonization, image classification, face verification, landmark detection, and layout prediction. Extensive qualitative and quantitative experimental results demonstrate the appealing performance of GH-Feat.


title: Quasi-symmetry groups and many-body scar dynamics

abstract: In quantum systems, a subspace spanned by degenerate eigenvectors of the Hamiltonian may have higher symmetries than those of the Hamiltonian itself. When this enhanced-symmetry group can be generated from local operators, we call it a quasi-symmetry group. When the group is a Lie group, an external field coupled to certain generators of the quasi-symmetry group lifts the degeneracy, and results in exactly periodic dynamics within the degenerate subspace, namely the many-body-scar dynamics (given that Hamiltonian is non-integrable). We provide two related schemes for constructing one-dimensional spin models having on-demand quasi-symmetry groups, with exact periodic evolution of a pre-chosen product or matrix-product state under certain external fields.


title: Towards a complete mass spectrum of type-IIB flux vacua at large complex structure

abstract: The large number of moduli fields arising in a generic string theory compactification makes a complete computation of the low energy effective theory infeasible. A common strategy to solve this problem is to consider Calabi-Yau manifolds with discrete symmetries, which effectively reduce the number of moduli and make the computation of the truncated Effective Field Theory possible. In this approach, however, the couplings (e.g., the masses) of the truncated fields are left undetermined. In the present paper we discuss the tree-level mass spectrum of type-IIB flux compactifications at Large Complex Structure, focusing on models with a reduced one-dimensional complex structure sector. We compute the tree-level spectrum for the dilaton and complex structure moduli, \emph{including the truncated fields}, which can be expressed entirely in terms of the known couplings of the reduced theory. We show that the masses of this set of fields are naturally heavy at vacua consistent with the KKLT construction, and we discuss other phenomenologically interesting scenarios where the spectrum involves fields much lighter than the gravitino. We also derive the probability distribution for the masses on the ensemble of flux vacua, and show that it exhibits universal features independent of the details of the compactification. We check our results on a large sample of flux vacua constructed in an orientifold of the Calabi-Yau $\mathbb{WP}^4_{[1,1,1,1,4]}$. Finally, we also discuss the conditions under which the spectrum derived here could arise in more general compactifications.


title: Integrable semi-discretizations of the Davey-Stewartson system and a $(2+1)$-dimensional Yajima-Oikawa system. II

abstract: This is a continuation of our previous paper arXiv:1904.07924, which is devoted to the construction of integrable semi-discretizations of the Davey-Stewartson system and a $(2+1)$-dimensional Yajima-Oikawa system; in this series of papers, we refer to a discretization of one of the two spatial variables as a semi-discretization. In this paper, we construct an integrable semi-discrete Davey-Stewartson system, which is essentially different from the semi-discrete Davey-Stewartson system proposed in the previous paper arXiv:1904.07924. We first obtain integrable semi-discretizations of the two elementary flows that compose the Davey-Stewartson system by constructing their Lax-pair representations and show that these two elementary flows commute as in the continuous case. Then, we consider a linear combination of the two elementary flows to obtain a new integrable semi-discretization of the Davey-Stewartson system. Using a linear transformation of the continuous independent variables, one of the two elementary Davey-Stewartson flows can be identified with an integrable semi-discretization of the $(2+1)$-dimensional Yajima-Oikawa system proposed in https://link.aps.org/doi/10.1103/PhysRevE.91.062902 .


title: Excellence, F-singularities, and solidity

abstract: An $R$-algebra $S$ is $R$-solid if there exists a nonzero $R$-linear map $S \rightarrow R$. In characteristic $p$, the study of $F$-singularities such as Frobenius splittings implicitly rely on the $R$-solidity of $R^{1/p}$. Following recent results of the first two authors on the Frobenius non-splitting of certain excellent $F$-pure rings, in this paper we use the notion of solidity to systematically study the notion of excellence, with an emphasis on $F$-singularities. We show that for rings $R$ essentially of finite type over complete local rings of characteristic $p$, reducedness implies the $R$-solidity of $R^{1/p}$, $F$-purity implies Frobenius splitting, and $F$-pure regularity implies split $F$-regularity. We demonstrate that Henselizations and completions are not solid, providing obstructions for the $R$-solidity of $R^{1/p}$ for arbitrary excellent rings. This also has negative consequences for the solidity of big Cohen-Macaulay algebras, an important example of which are absolute integral closures of excellent local rings in prime characteristic. We establish a close relationship between the solidity of absolute integral closures and the notion of Japanese rings. Analyzing the Japanese property reveals that Dedekind domains $R$ for which $R^{1/p}$ is $R$-solid are excellent, despite our recent examples of excellent Euclidean domains with no nonzero $p^{-1}$-linear maps. Additionally, we show that while perfect closures are often solid in algebro-geometric situations, there exist locally excellent domains with solid perfect closures whose absolute integral closures are not solid. In an appendix, Karen E. Smith uses the solidity of absolute integral closures to characterize the test ideal for a large class of Gorenstein domains of prime characteristic.


title: Off-resonant coherent electron transport over three nanometers in multi-heme protein bioelectronic junctions

abstract: Multi-heme cytochromes (MHC) are fascinating proteins used by bacterial organisms to shuttle electrons within and between their cells. When placed in a solid state electronic junction, they support temperature-independent currents over several nanometers that are three orders of magnitude higher compared to other redox proteins of comparable size. To gain microscopic insight into their astonishingly high conductivities, we present herein the first current-voltage calculations of its kind, for a MHC sandwiched between two Au(111) electrodes, complemented by photo-emission spectroscopy experiments. We find that conduction proceeds via off-resonant coherent tunneling mediated by a large number of protein valence-band orbitals that are strongly delocalized over heme and protein residues, effectively "gating" the current between the two electrodes. This picture is profoundly different from the dominant electron hopping mechanism supported by the same protein in aqueous solution. Our results imply that current output in MHC junctions could be even further increased in the resonant regime, e.g. by application of a gate voltage, making these proteins extremely interesting for next-generation bionanoelectronic devices.


title: Support Aggregate Analytic Window Function over Large Data by Spilling

abstract: Analytic function, also called window function, is to query the aggregation of data over a sliding window. For example, a simple query over the online stock platform is to return the average price of a stock of the last three days. These functions are commonly used features in SQL databases. They are supported in most of the commercial databases. With the increasing usage of cloud data infra and machine learning technology, the frequency of queries with analytic window functions rises. Some analytic functions only require const space in memory to store the state, such as SUM, AVG, while others require linear space, such as MIN, MAX. When the window is extremely large, the memory space to store the state may be too large. In this case, we need to spill the state to disk, which is a heavy operation. In this paper, we proposed an algorithm to manipulate the state data in the disk to reduce the disk I/O to make spill available and efficienct. We analyze the complexity of the algorithm with different data distribution.


title: Certain subclasses of Spirallike univalent functions related with Pascal distribution series

abstract: The purpose of the present paper is to find the necessary and sufficient conditions for the subclasses of analytic functions associated with Pascal distribution to be in subclasses of spiral-like univalent functions and inclusion relations for such subclasses in the open unit disk $\mathbb{D}.$ Further, we consider the properties of integral operator related to Pascal distribution series. Several corollaries and consequences of the main results are also considered.


title: Restricting Supercuspidal Representations via a Restriction of Data

abstract: Let $F$ be a non-archimedean local field of residual characteristic $p$. Let $\mathbb{G}$ be a reductive group defined over $F$ which splits over a tamely ramified extension and set $G=\mathbb{G}(F)$. We assume that $p$ does not divide the order of the Weyl group of $\mathbb{G}$. Given a closed connected $F$-subgroup $\mathbb{H}$ that contains the derived subgroup of $\mathbb{G}$, we study the restriction to $H$ of an irreducible supercuspidal representation $\pi=\pi_G(\Psi)$ of $G$, where $\Psi$ is a $G$-datum as per the J.K. Yu Construction. We provide a full description of $\pi|_H$ into irreducible components, with multiplicity, via a restriction of data which constructs $H$-data from $\Psi$. Analogously, we define a restriction of Kim-Yu types to study the restriction of irreducible representations of $G$ which are not supercuspidal.


title: Pulsar Timing Signatures of Circumbinary Asteroid Belts

abstract: The gravitational pull of a large number of asteroids perturbs a pulsar's motion to a degree that is detectable through precision timing of millisecond pulsars. The result is a low-frequency, correlated noise process, similar in form to the red timing noise known to affect canonical pulsars, or to the signal expected from a stochastic gravitational-wave background. Motivated by the observed fact that many millisecond pulsars are in binary systems, we describe the ways in which the presence of a binary companion to the pulsar would affect the signal produced by an asteroid belt. The primary effect of the companion is to destabilize the shortest-period orbits, cutting off the high-frequency component of the signal from the asteroid belt. We also discuss the implications of asteroid belts for gravitational-wave search efforts. Compared to the signal from a stochastic gravitational-wave background, asteroid belt noise has a similar frequency and amplitude, and is similarly independent of radio frequency, but is not correlated between different pulsars, which should allow the two kinds of signal to be distinguished.


title: Generalizing Variational Autoencoders with Hierarchical Empirical Bayes

abstract: Variational Autoencoders (VAEs) have experienced recent success as data-generating models by using simple architectures that do not require significant fine-tuning of hyperparameters. However, VAEs are known to suffer from over-regularization which can lead to failure to escape local maxima. This phenomenon, known as posterior collapse, prevents learning a meaningful latent encoding of the data. Recent methods have mitigated this issue by deterministically moment-matching an aggregated posterior distribution to an aggregate prior. However, abandoning a probabilistic framework (and thus relying on point estimates) can both lead to a discontinuous latent space and generate unrealistic samples. Here we present Hierarchical Empirical Bayes Autoencoder (HEBAE), a computationally stable framework for probabilistic generative models. Our key contributions are two-fold. First, we make gains by placing a hierarchical prior over the encoding distribution, enabling us to adaptively balance the trade-off between minimizing the reconstruction loss function and avoiding over-regularization. Second, we show that assuming a general dependency structure between variables in the latent space produces better convergence onto the mean-field assumption for improved posterior inference. Overall, HEBAE is more robust to a wide-range of hyperparameter initializations than an analogous VAE. Using data from MNIST and CelebA, we illustrate the ability of HEBAE to generate higher quality samples based on FID score than existing autoencoder-based approaches.


title: Testing linear inequalities of subgraph statistics

abstract: Property testers are fast randomized algorithms whose task is to distinguish between inputs satisfying some predetermined property ${\cal P}$ and those that are far from satisfying it. Since these algorithms operate by inspecting a small randomly selected portion of the input, the most natural property one would like to be able to test is whether the input does not contain certain forbidden small substructures. In the setting of graphs, such a result was obtained by Alon et al., who proved that for any finite family of graphs ${\cal F}$, the property of being induced ${\cal F}$-free (i.e. not containing an induced copy of any $F \in {\cal F}$) is testable. It is natural to ask if one can go one step further and prove that more elaborate properties involving induced subgraphs are also testable. One such generalization of the result of Alon et al. was formulated by Goldreich and Shinkar who conjectured that for any finite family of graphs ${\cal F}$, and any linear inequality involving the densities of the graphs $F \in {\cal F}$ in the input graph, the property of satisfying this inequality can be tested in a certain restricted model of graph property testing. Our main result in this paper disproves this conjecture in the following strong form: some properties of this type are not testable even in the classical (i.e. unrestricted) model of graph property testing. The proof deviates significantly from prior non-testability results in this area. The main idea is to use a linear inequality relating induced subgraph densities in order to encode the property of being a quasirandom graph.


title: Technicolor coupled models

abstract: When technicolor (TC), QCD, extended technicolor (ETC) and other interactions become coupled through their different Schwinger-Dyson equations, the solution of these equations are modified compared to those of the isolated equations. The change in the self-energies is similar to that obtained in the presence of four-fermion interactions, but without their ad hoc inclusion in the theory. In this case the TC and QCD self-energies decrease logarithmically with the momenta, which allows us to build models where ETC boson masses can be pushed to very high energies, and do not lead to undesirable flavor changing interactions. Viable TC models may be built along this line including a necessary horizontal symmetry. The different fermionic mass scales are dictated by the different strong interactions. Pseudo-Goldstone bosons acquire large masses in this class of models.


title: Statistics for dark matter subhalo searches in gamma rays from a kinematically constrained population model: Fermi-LAT-like telescopes

abstract: Cold dark matter subhalos are expected to populate galaxies in numbers. If dark matter self-annihilates, these objects turn into prime targets for indirect searches, in particular with gamma-ray telescopes. Incidentally, the Fermi-LAT catalog already contains many unidentified sources that might be associated with subhalos. In this paper, we determine the probability for subhalos to be identified as gamma-ray pointlike sources from their predicted distribution properties. We use a semi-analytical model for the Galactic subhalo population, which, in contrast to cosmological simulations, can be made fully consistent with current kinematic constraints in the Milky Way and has no resolution limit. The model incorporates tidal stripping effects from a realistic distribution of baryons in the Milky Way. The same baryonic distribution contributes a diffuse gamma-ray foreground which adds up to that, often neglected in subhalo searches, generated by the smooth dark matter and the unresolved subhalos. This configuration implies a correlation between pointlike subhalo signals and diffuse background. Based on this semi-analytical modeling, we generate mock gamma-ray data assuming an idealized telescope resembling Fermi-LAT and perform a likelihood analysis to estimate the current and future sensitivity to subhalos in the relevant parameter space. We find a number of detectable subhalos of order ${\cal O}(<1)$ for optimistic model parameters and a WIMP mass of 100~GeV, maximized for a cored host halo. This barely provides support to the current interpretation of several Fermi unidentified sources as subhalos. We also find it more likely to detect the smooth Galactic halo itself before subhalos, should dark matter in the GeV-TeV mass range self-annihilate through $s$-wave processes.

categories: astro-ph.HE astro-ph.CO astro-ph.GA hep-ph


title: A coherent likelihood parametrization for doubly robust estimation of a causal effect with missing confounders

abstract: Missing data and confounding are two problems researchers face in observational studies for comparative effectiveness. Williamson et al. (2012) recently proposed a unified approach to handle both issues concurrently using a multiply-robust (MR) methodology under the assumption that confounders are missing at random. Their approach considers a union of models in which any submodel has a parametric component while the remaining models are unrestricted. We show that while their estimating function is MR in theory, the possibility for multiply robust inference is complicated by the fact that parametric models for different components of the union model are not variation independent and therefore the MR property is unlikely to hold in practice. To address this, we propose an alternative transparent parametrization of the likelihood function, which makes explicit the model dependencies between various nuisance functions needed to evaluate the MR efficient score. The proposed method is genuinely doubly-robust (DR) in that it is consistent and asymptotic normal if one of two sets of modeling assumptions holds. We evaluate the performance and doubly robust property of the DR method via a simulation study.


title: Translation Between Waves, wave2wave

abstract: The understanding of sensor data has been greatly improved by advanced deep learning methods with big data. However, available sensor data in the real world are still limited, which is called the opportunistic sensor problem. This paper proposes a new variant of neural machine translation seq2seq to deal with continuous signal waves by introducing the window-based (inverse-) representation to adaptively represent partial shapes of waves and the iterative back-translation model for high-dimensional data. Experimental results are shown for two real-life data: earthquake and activity translation. The performance improvements of one-dimensional data was about 46% in test loss and that of high-dimensional data was about 1625% in perplexity with regard to the original seq2seq.


title: One-loop Graviton Corrections to Conformal Scalars on a de Sitter Background

abstract: We exploit a recent computation of one graviton loop corrections to the self-mass [1] to quantum-correct the field equation for a massless, conformally coupled scalar on a de Sitter background. With the obvious choice for the finite part of the $R^2 \phi^2$ counterterm, we find that neither plane wave mode functions nor the response to a point source acquires large infrared logarithms. However, we do find a decaying logarithmic correction to the mode function and a short distance logarithmic running of the potential in addition to the power-law effect inherited from flat space.


title: NSGANetV2: Evolutionary Multi-Objective Surrogate-Assisted Neural Architecture Search

abstract: In this paper, we propose an efficient NAS algorithm for generating task-specific models that are competitive under multiple competing objectives. It comprises of two surrogates, one at the architecture level to improve sample efficiency and one at the weights level, through a supernet, to improve gradient descent training efficiency. On standard benchmark datasets (C10, C100, ImageNet), the resulting models, dubbed NSGANetV2, either match or outperform models from existing approaches with the search being orders of magnitude more sample efficient. Furthermore, we demonstrate the effectiveness and versatility of the proposed method on six diverse non-standard datasets, e.g. STL-10, Flowers102, Oxford Pets, FGVC Aircrafts etc. In all cases, NSGANetV2s improve the state-of-the-art (under mobile setting), suggesting that NAS can be a viable alternative to conventional transfer learning approaches in handling diverse scenarios such as small-scale or fine-grained datasets. Code is available at https://github.com/mikelzc1990/nsganetv2


title: CACTI: Captcha Avoidance via Client-side TEE Integration

abstract: Preventing abuse of web services by bots is an increasingly important problem, as abusive activities grow in both volume and variety. CAPTCHAs are the most common way for thwarting bot activities. However, they are often ineffective against bots and frustrating for humans. In addition, some recent CAPTCHA techniques diminish user privacy. Meanwhile, client-side Trusted Execution Environments (TEEs) are becoming increasingly widespread (notably, ARM TrustZone and Intel SGX), allowing establishment of trust in a small part (trust anchor or TCB) of client-side hardware. This prompts the question: can a TEE help reduce (or remove entirely) user burden of solving CAPTCHAs? In this paper, we design CACTI: CAPTCHA Avoidance via Client-side TEE Integration. Using client-side TEEs, CACTI allows legitimate clients to generate unforgeable rate-proofs demonstrating how frequently they have performed specific actions. These rate-proofs can be sent to web servers in lieu of solving CAPTCHAs. CACTI provides strong client privacy guarantees, since the information is only sent to the visited website and authenticated using a group signature scheme. Our evaluations show that overall latency of generating and verifying a CACTI rate-proof is less than 0.25 sec, while CACTI's bandwidth overhead is over 98% lower than that of current CAPTCHA systems.


title: Stacking Domain Wall Magnons in Twisted van der Waals Magnets

abstract: Using bilayer CrI$_3$ as an example, we demonstrate that stacking domain walls in van der Waals magnets can host one dimensional (1D) magnon channels, which have lower energies than bulk magnons. Interestingly, some magnon channels are hidden in magnetically homogeneous background and can only be inferred with the knowledge of stacking domain walls. Compared to 1D magnons confined in magnetic domain walls, 1D magnons in stacking domain walls are more stable against external perturbations. We show that the relaxed moir\'e superlattices of small-angle twisted bilayer CrI$_3$ is a natural realization of stacking domain walls and host interconnected moir\'e magnon network. Our work reveals the importance of stacking domain walls in understanding magnetic properties of van der Waals magnets, and extends the scope of stacking engineering to magnetic dynamics.


title: Automatic Story Construction from News Articles in an Online Fashion

abstract: This paper presents a novel story construction system to track the evolution of stories in an online fashion. The proposed system uses a novel sliding window solution, named Inching Window, allowing the processing of each new data element on-the-fly. To assign a new data element into a community in a fast and memory-efficient manner, we apply the modularity maximization idea of Louvain method on-the-fly. As part of the experimental validation, we provide step by step construction of a meaningful news story and support the case with a set of visualizations.


title: Phase Transition Behavior in Knowledge Compilation

abstract: The study of phase transition behaviour in SAT has led to deeper understanding and algorithmic improvements of modern SAT solvers. Motivated by these prior studies of phase transitions in SAT, we seek to study the behaviour of size and compile-time behaviour for random k-CNF formulas in the context of knowledge compilation. We perform a rigorous empirical study and analysis of the size and runtime behavior for different knowledge compilation forms (and their corresponding compilation algorithms): d-DNNFs, SDDs and OBDDs across multiple tools and compilation algorithms. We employ instances generated from the random k-CNF model with varying generation parameters to empirically reason about the expected and median behavior of size and compilation-time for these languages. Our work is similar in spirit to the early work in CSP community on phase transition behavior in SAT/CSP. In a similar spirit, we identify the interesting behavior with respect to different parameters: clause density and solution density, a novel control parameter that we identify for the study of phase transition behavior in the context of knowledge compilation. Furthermore, we summarize our empirical study in terms of two concrete conjectures; a rigorous study of these conjectures will possibly require new theoretical tools.


title: Robust-Adaptive Interval Predictive Control for Linear Uncertain Systems

abstract: We consider the problem of stabilization of a linear system, under state and control constraints, and subject to bounded disturbances and unknown parameters in the state matrix. First, using a simple least square solution and available noisy measurements, the set of admissible values for parameters is evaluated. Second, for the estimated set of parameter values and the corresponding linear interval model of the system, two interval predictors are recalled and an unconstrained stabilizing control is designed that uses the predicted intervals. Third, to guarantee the robust constraint satisfaction, a model predictive control algorithm is developed, which is based on solution of an optimization problem posed for the interval predictor. The conditions for recursive feasibility and asymptotic performance are established. Efficiency of the proposed control framework is illustrated by numeric simulations.


title: Scotogenic dark matter in an orbifold theory of flavor

abstract: We propose a flavour theory in which the family symmetry results naturally from a six-dimensional orbifold compactification. "Diracness" of neutrinos is a consequence of the spacetime dimensionality, and the fact that right-handed neutrinos live in the bulk. Dark matter is incorporated in a scotogenic way, as a result of an auxiliary $\mathbb{Z}_3$ symmetry, and its stability is associated to the conservation of a "dark parity" symmetry. The model leads naturally to a "golden" quark-lepton mass relation.


title: Global gender differences in Wikipedia readership

abstract: Wikipedia represents the largest and most popular source of encyclopedic knowledge in the world today, aiming to provide equal access to information worldwide. From a global online survey of 65,031 readers of Wikipedia and their corresponding reading logs, we present novel evidence of gender differences in Wikipedia readership and how they manifest in records of user behavior. More specifically we report that (1) women are underrepresented among readers of Wikipedia, (2) women view fewer pages per reading session than men do, (3) men and women visit Wikipedia for similar reasons, and (4) men and women exhibit specific topical preferences. Our findings lay the foundation for identifying pathways toward knowledge equity in the usage of online encyclopedic knowledge.


title: Glassy dynamics and memory effects in an intrinsically disordered protein construct

abstract: Glassy, nonexponential relaxations in globular proteins are typically attributed to conformational behaviors that are missing from intrinsically disordered proteins. Yet, we show that single molecules of a disordered-protein construct display two signatures of glassy dynamics, logarithmic relaxations and a Kovacs memory effect, in response to changes in applied tension. We attribute this to the presence of multiple independent local structures in the chain, which we corroborate with a model that correctly predicts the force dependence of the relaxation. The mechanism established here likely applies to other disordered proteins.


title: Collecting Service-Based Maintainability Metrics from RESTful API Descriptions: Static Analysis and Threshold Derivation

abstract: While many maintainability metrics have been explicitly designed for service-based systems, tool-supported approaches to automatically collect these metrics are lacking. Especially in the context of microservices, decentralization and technological heterogeneity may pose challenges for static analysis. We therefore propose the modular and extensible RAMA approach (RESTful API Metric Analyzer) to calculate such metrics from machine-readable interface descriptions of RESTful services. We also provide prototypical tool support, the RAMA CLI, which currently parses the formats OpenAPI, RAML, and WADL and calculates 10 structural service-based metrics proposed in scientific literature. To make RAMA measurement results more actionable, we additionally designed a repeatable benchmark for quartile-based threshold ranges (green, yellow, orange, red). In an exemplary run, we derived thresholds for all RAMA CLI metrics from the interface descriptions of 1,737 publicly available RESTful APIs. Researchers and practitioners can use RAMA to evaluate the maintainability of RESTful services or to support the empirical evaluation of new service interface metrics.


title: Hermite functions and Fourier series

abstract: Using normalized Hermite functions, we construct bases in the space of square integrable functions on the unit circle ($L^2(\mathcal C)$) and in $l_2(\mathbb Z)$, which are related to each other by means of the Fourier transform and the discrete Fourier transform. These relations are unitary. The construction of orthonormal bases requires the use of the Gramm--Schmidt method. On both spaces, we have provided ladder operators with the same properties as the ladder operators for the one-dimensional quantum oscillator. These operators are linear combinations of some multiplication- and differentiation-like operators that, when applied to periodic functions, preserve periodicity. Finally, we have constructed riggings for both $L^2(\mathcal C)$ and $l_2(\mathbb Z)$, so that all the mentioned operators are continuous.


title: Fusing Concurrent Orthogonal Wide-aperture Sonar Images for Dense Underwater 3D Reconstruction

abstract: We propose a novel approach to handling the ambiguity in elevation angle associated with the observations of a forward looking multi-beam imaging sonar, and the challenges it poses for performing an accurate 3D reconstruction. We utilize a pair of sonars with orthogonal axes of uncertainty to independently observe the same points in the environment from two different perspectives, and associate these observations. Using these concurrent observations, we can create a dense, fully defined point cloud at every time-step to aid in reconstructing the 3D geometry of underwater scenes. We will evaluate our method in the context of the current state of the art, for which strong assumptions on object geometry limit applicability to generalized 3D scenes. We will discuss results from laboratory tests that quantitatively benchmark our algorithm's reconstruction capabilities, and results from a real-world, tidal river basin which qualitatively demonstrate our ability to reconstruct a cluttered field of underwater objects.


title: PDO-eConvs: Partial Differential Operator Based Equivariant Convolutions

abstract: Recent research has shown that incorporating equivariance into neural network architectures is very helpful, and there have been some works investigating the equivariance of networks under group actions. However, as digital images and feature maps are on the discrete meshgrid, corresponding equivariance-preserving transformation groups are very limited. In this work, we deal with this issue from the connection between convolutions and partial differential operators (PDOs). In theory, assuming inputs to be smooth, we transform PDOs and propose a system which is equivariant to a much more general continuous group, the $n$-dimension Euclidean group. In implementation, we discretize the system using the numerical schemes of PDOs, deriving approximately equivariant convolutions (PDO-eConvs). Theoretically, the approximation error of PDO-eConvs is of the quadratic order. It is the first time that the error analysis is provided when the equivariance is approximate. Extensive experiments on rotated MNIST and natural image classification show that PDO-eConvs perform competitively yet use parameters much more efficiently. Particularly, compared with Wide ResNets, our methods result in better results using only 12.6% parameters.


title: Exploring student facility with "goes like'' reasoning in introductory physics

abstract: Covariational reasoning -- reasoning about how changes in one quantity relate to changes in another quantity -- has been examined extensively in mathematics education research. Little research has been done, however, on covariational reasoning in introductory physics contexts. We explore one aspect of covariational reasoning: ``goes like'' reasoning. ``Goes like'' reasoning refers to ways physicists relate two quantities through a simplified function. For example, physicists often say that ``the electric field goes like one over r squared.'' While this reasoning mode is used regularly by physicists and physics instructors, how students make sense of and use it remains unclear. We present evidence from reasoning inventory items which indicate that many students are sense making with tools from prior math instruction, that could be developed into expert ``goes like'' thinking with direct instruction. Recommendations for further work in characterizing student sense making as a foundation for future development of instruction are made.


title: Volumetric heating of nanowire arrays to keV temperatures using kilojoule-scale petawatt laser interactions

abstract: We present picosecond-resolution streaked K-shell spectra from 400 nm-diameter nickel nanowire arrays, demonstrating the ability to generate large volumes of high energy density plasma when combined with the longer pulses typical of the largest short pulse lasers. After irradiating the wire array with 100 J, 600 fs ultra-high-contrast laser pulses focussed to $>10^{20}$ W/cm$^{2}$ at the Orion laser facility, we combine atomic kinetics modeling of the streaked spectra with 2D collisional particle-in-cell simulations to describe the evolution of material conditions within these samples for the first time. We observe a three-fold enhancement of helium-like emission compared to a flat foil in a near-solid-density plasma sustaining keV temperatures for tens of picoseconds, the result of strong electric return currents heating the wires and causing them to explode and collide.


title: Assessment of a symmetry preserving JFNK method for atmospheric convection

abstract: Numerical simulations of nonhydrostatic atmospheric flow, based on linearly decoupled semi-implicit or fully-implicit techniques, usually solve linear systems by a pre-conditioned Krylov method without preserving the skew-symmetry of convective operators. We propose to perform atmospheric simulations in such a fully-implicit manner that the difference operators preserve both the skew-symmetry and the tightly nonlinear coupling of the differential operators. We demonstrate that a symmetry-preserving Jacobian-free Newton-Krylov~(JFNK) method mimics a balance between convective transport and turbulence dissipation. We present a wavelet method as an effective symmetry preserving discretization technique. The symmetry-preserving JFNK method for solving equations of nonhydrostatic atmospheric flows has been examined using two benchmark simulations of penetrative convection -- a) dry thermals rising in a neutrally stratified and stably stratified environment, and b) urban heat island circulations for effects of the surface heat flux $H_0$ varying in the range of $25 \le H_0 \le 930$Wm$^{-2}$.The results show that an eddy viscosity model provides the necessary dissipation of the subgrid-scale modes, while the symmetry-preserving JFNK method provides the conservation of mass and energy at a satisfactory level. Comparisons of the results from a laboratory experiment of heat island circulation and a field measurement of potential temperature also suggest the modelling accuracy of the present symmetry-preserving JFNK framework.


title: Randomized Automatic Differentiation

abstract: The successes of deep learning, variational inference, and many other fields have been aided by specialized implementations of reverse-mode automatic differentiation (AD) to compute gradients of mega-dimensional objectives. The AD techniques underlying these tools were designed to compute exact gradients to numerical precision, but modern machine learning models are almost always trained with stochastic gradient descent. Why spend computation and memory on exact (minibatch) gradients only to use them for stochastic optimization? We develop a general framework and approach for randomized automatic differentiation (RAD), which can allow unbiased gradient estimates to be computed with reduced memory in return for variance. We examine limitations of the general approach, and argue that we must leverage problem specific structure to realize benefits. We develop RAD techniques for a variety of simple neural network architectures, and show that for a fixed memory budget, RAD converges in fewer iterations than using a small batch size for feedforward networks, and in a similar number for recurrent networks. We also show that RAD can be applied to scientific computing, and use it to develop a low-memory stochastic gradient method for optimizing the control parameters of a linear reaction-diffusion PDE representing a fission reactor.


title: Sample Shape and Boundary Dependence of Measured Transverse Thermal Properties

abstract: Despite increased interest in thermal Hall measurements for the analysis of insulating quantum materials, there remains large uncertainty in such measurements due to contact misalignment. In this paper we propose that sample geometry and uncertain boundary conditions may account for uncertainty in the measurement of $D_{xy}$ or $\kappa_{xy}$ as well. By running simple simulations in an open source finite-element solver, we demonstrate that measured $D_{xy}$ in a thermal Hall bar can be changed by a factor of order unity in samples with similar width and length. This geometric corrective factor depends on the distinction between uniform heat flow and constant temperature boundary couplings to a bath. Sample geometry and boundary conditions can be accounted for through simulation or by using rectangular samples to make thermal Hall measurements more reliable and reproducible when the amplitude of $\kappa_{xy}$ is important. Finally, we propose a contactless optical method for measuring $D_{xy}$ which is insensitive to the longitudinal diffusivity pollution caused by contact misalignment.


title: Topology of many-body edge and extended quantum states in an open spin chain: 1/3--plateau, Kosterlitz-Thouless transition, and Luttinger liquid

abstract: Quantum many-body edge and extended magnon excitations from the 1/3 -- plateau of the anisotropic Heisenberg model on an open AB$_2$ chain in a magnetic field $h$ are unveiled using the density matrix renormalization group and exact diagonalization. By tuning both the anisotropy and $h$ in the rich phase diagram, the edge states penetrate in the bulk, whose gap closes in a symmetry-protected topological Kosterlitz-Thouless transition. Also, we witness the squeezed chain effect, the breaking of the edge states degeneracy, and a topological change of the excitations from gapped magnons with quadratic long-wavelength dispersion to a linear spinon dispersion in the Luttinger liquid gapless phase as the anisotropy $\lambda$ approaches the critical point from the $\lambda>0$ side of the phase diagram.


title: The Historical Impact of Anthropogenic Climate Change on Global Agricultural Productivity

abstract: Agricultural research has fostered productivity growth, but the historical influence of anthropogenic climate change on that growth has not been quantified. We develop a robust econometric model of weather effects on global agricultural total factor productivity (TFP) and combine this model with counterfactual climate scenarios to evaluate impacts of past climate trends on TFP. Our baseline model indicates that anthropogenic climate change has reduced global agricultural TFP by about 21% since 1961, a slowdown that is equivalent to losing the last 9 years of productivity growth. The effect is substantially more severe (a reduction of ~30-33%) in warmer regions such as Africa and Latin America and the Caribbean. We also find that global agriculture has grown more vulnerable to ongoing climate change.


title: Integrative Analysis for COVID-19 Patient Outcome Prediction

abstract: While image analysis of chest computed tomography (CT) for COVID-19 diagnosis has been intensively studied, little work has been performed for image-based patient outcome prediction. Management of high-risk patients with early intervention is a key to lower the fatality rate of COVID-19 pneumonia, as a majority of patients recover naturally. Therefore, an accurate prediction of disease progression with baseline imaging at the time of the initial presentation can help in patient management. In lieu of only size and volume information of pulmonary abnormalities and features through deep learning based image segmentation, here we combine radiomics of lung opacities and non-imaging features from demographic data, vital signs, and laboratory findings to predict need for intensive care unit (ICU) admission. To our knowledge, this is the first study that uses holistic information of a patient including both imaging and non-imaging data for outcome prediction. The proposed methods were thoroughly evaluated on datasets separately collected from three hospitals, one in the United States, one in Iran, and another in Italy, with a total 295 patients with reverse transcription polymerase chain reaction (RT-PCR) assay positive COVID-19 pneumonia. Our experimental results demonstrate that adding non-imaging features can significantly improve the performance of prediction to achieve AUC up to 0.884 and sensitivity as high as 96.1%, which can be valuable to provide clinical decision support in managing COVID-19 patients. Our methods may also be applied to other lung diseases including but not limited to community acquired pneumonia. The source code of our work is available at https://github.com/DIAL-RPI/COVID19-ICUPrediction.


title: Bayesian Few-Shot Classification with One-vs-Each P\'olya-Gamma Augmented Gaussian Processes

abstract: Few-shot classification (FSC), the task of adapting a classifier to unseen classes given a small labeled dataset, is an important step on the path toward human-like machine learning. Bayesian methods are well-suited to tackling the fundamental issue of overfitting in the few-shot scenario because they allow practitioners to specify prior beliefs and update those beliefs in light of observed data. Contemporary approaches to Bayesian few-shot classification maintain a posterior distribution over model parameters, which is slow and requires storage that scales with model size. Instead, we propose a Gaussian process classifier based on a novel combination of P\'olya-Gamma augmentation and the one-vs-each softmax approximation that allows us to efficiently marginalize over functions rather than model parameters. We demonstrate improved accuracy and uncertainty quantification on both standard few-shot classification benchmarks and few-shot domain transfer tasks.


title: The arrow of time (second law) as a randomness-driven emergent property of large systems

abstract: The arrow of time is an irreversible phenomenon for a system of particles undergoing reversible dynamics. Since the time of Boltzmann to this day, the arrow of time has led to debate and research. However, the enormous growth of nanotechnology and associated experimental techniques has brought the arrow of time at the forefront because of its practical implications. Using simulations of one-dimensional diffusion of a system of particles, we show that the arrow of time is an emergent property of a large system. We show that the recurrence time for a system of particles to return to its original configuration grows rapidly as the number of particles grows. Based on the simulations, we have provided the expressions for recurrence times for classical particles, Fermions, and Bosons. A system of Bosons has the shortest recurrence time, whereas a system of classical particles has the longest recurrence time. The underlying distribution around the mean recurrence time is Poisson-distributed for Bosons and Gaussian-distributed for Fermions and classical particles. The probabilistic approach to encode dynamics enables testing processes other than diffusion and quantify their effects on the recurrence time.


title: Visual Testing of GUIs by Abstraction

abstract: Ensuring the correct visual appearance of graphical user interfaces (GUIs) is important because visual bugs can cause substantial losses for businesses. An application might behave functionally correct in an automated test, but visual bugs can make the GUI effectively unusable for the user. Most of today's approaches for visual testing are pixel-based and tend to have flaws that are characteristic for image differencing. For instance, minor and unimportant visual changes often cause false positives, which confuse the user with unnecessary error reports. Our idea is to introduce an abstract GUI state (AGS), where we define structural relations to identify relevant GUI changes and ignore those that are unimportant from the user's point of view. In addition, we explore several strategies to address the GUI element identification problem in terms of AGS. This allows us to provide rich diagnostic information that help the user to better interpret changes. Based on the principles of golden master testing, we can support a fully-automated approach to visual testing by using the AGS. We have implemented our approach to visually test web pages and our experiments show that we are able to reliably detect GUI changes.


title: Non-Markov Policies to Reduce Sequential Failures in Robot Bin Picking

abstract: A new generation of automated bin picking systems using deep learning is evolving to support increasing demand for e-commerce. To accommodate a wide variety of products, many automated systems include multiple gripper types and/or tool changers. However, for some objects, sequential grasp failures are common: when a computed grasp fails to lift and remove the object, the bin is often left unchanged; as the sensor input is consistent, the system retries the same grasp over and over, resulting in a significant reduction in mean successful picks per hour (MPPH). Based on an empirical study of sequential failures, we characterize a class of "sequential failure objects" (SFOs) -- objects prone to sequential failures based on a novel taxonomy. We then propose three non-Markov picking policies that incorporate memory of past failures to modify subsequent actions. Simulation experiments on SFO models and the EGAD dataset suggest that the non-Markov policies significantly outperform the Markov policy in terms of the sequential failure rate and MPPH. In physical experiments on 50 heaps of 12 SFOs the most effective Non-Markov policy increased MPPH over the Dex-Net Markov policy by 107%.


title: Curl-curl conforming elements on tetrahedra

abstract: In [24], we proposed H(curl^2)-conforming elements on both a triangle and a rectangle. This family of elements provides a brand new method to solve the quad-curl problem in 2 dimensions. In this paper, we turn our focus to 3 dimensions and construct an H(curl^2)-conforming tetrahedral finite element. The newly proposed element has been proved to have the optimal interpolation error estimate. Having tetrahedral elements, we can solve the quad-curl problem in any Lipschitz domain by conforming finite element method. We also provide several numerical examples of using our element to solve the quad-curl problem. The results of the numerical experiments show the effectiveness and correctness of our element.


title: The quantum non-linear $\sigma$-model RG flow and integrability in wormhole geometries

abstract: The target space of the non-linear $\sigma$-model is a Riemannian manifold. Although it can be any Riemannian metric, there are certain physically interesting geometries which are worth to study. Here, we numerically evolve the time-symmetric foliations of a family of spherically symmetric asymptotically flat wormholes under the $1$-loop renormalization group flow of the non-linear $\sigma$-model, the Ricci flow, and under the $2$-loop approximation, RG-2 flow. We rely over some theorems adapted from the compact case for studying the evolution of different wormhole types, specially those with high curvature zones. Some metrics expand and others contract at the beginning of the flow, however, all metrics pinch-off at a certain time. This is related with the fact that the flow does not converge to a fixed point when its starting geometry is the spatial sections of a Morris-Thorne wormhole, and therefore the corresponding non-linear $\sigma$-model is non-integrable/renormalizable. We present a numerical study of the evolution of wormhole singularities in three dimensions extending the theoretical estimations. Finally, we compute the evolution of the Hamilton's entropy and the Brown-York energy.


title: A new model for the effective thermal conductivity of polycrystalline solids

abstract: A new model for predicting the effective thermal conductivity of polycrystalline materials is presented. In contrast to existing models, our new model is based on the thin-interface description of grain boundaries (GBs) and treats GBs as an autonomous phase with its own thermal conductivity. This treatment allows the model to account for the effects of segregation/doping, interface structural/phase transition, and interface decohesion on the effective thermal conductivity. finite-element simulations were performed to validate the new model. The simulations proved its higher accuracy compared to existing models.


title: Thermodynamical Description of Hot, Rapidly Rotating Neutron Stars, Protoneutron Stars, and Neutron Star Merger Remnants

abstract: The prediction of the equation of state of hot, dense nuclear matter is one of the most complicated and interesting problems in nuclear astrophysics. At the same time, knowledge of it is the basic ingredient for some of the most interesting studies. In the present work, we concentrate our study on the construction of the equation of state of hot, dense nuclear matter, related mainly to the interior of the neutron star. We employ a theoretical nuclear model, which includes momentum-dependent interaction among the nucleons, along with the \textit{state-of-the-art} microscopic calculations. Thermal effects are introduced in a self-consistent way, and a set of isothermal and isentropic equations of state are predicted. The predicted equations of state are used in order to acquire and to extend the knowledge of the thermal effect on both nonrotating and rapidly rotating with the Kepler frequency neutron stars. The simultaneous study of thermal and rotation effects provides useful information on some of the most important quantities, including the mass (gravitational and baryon) and radius, the Kepler frequency and Kerr parameter, the moment of inertia, etc. These quantities are directly related to studies of protoneutron stars and mainly the hot and rapidly rotating remnant of a binary neutron star merger. Data from the late observations of binary neutron star mergers and the present study may offer useful tools for their investigation and help in providing possible constraints on the equation of state of nuclear matter.


title: Electrical control of spin relaxation anisotropy during drift transport in a two-dimensional electron gas

abstract: Spin relaxation was studied in a two-dimensional electron gas confined in a wide GaAs quantum well. Recently, the control of the spin relaxation anisotropy by diffusive motion was first shown in D. Iizasa et al., arXiv:2006.08253 (2020). Here, we demonstrate electrical control by drift transport in a system with two-subbands occupied. The combined effect of in-plane and gate voltages was investigated using time-resolved Kerr rotation. The measured relaxation time present strong anisotropy with respect to the transport direction. For an in-plane accelerating electric field along $\left[110\right]$, the lifetime was strongly suppressed irrespective of the applied gate voltage. Remarkably, for transport along $\left[1\bar{1}0\right]$, the data shows spin lifetime that was gate-dependent and longer than in the $\left[110\right]$ direction regardless of the in-plane voltage. In agreement, independent results of anisotropic spin precession frequencies are also presented. Nevertheless, the long spin lifetime, strong anisotropy and drift response seen in the data are beyond the existing models for spin drift and diffusion.


title: Ultracompact HII regions with extended emission: The case of G43.89-0.78 and its molecular environment

abstract: The Karl Jansky Very Large Array (VLA), Owens Valley Radio Observatory (OVRO), Atacama Large Millimetric Array (ALMA), and the infrared \textit{Spitzer} observatories, are powerful facilities to study massive star formation regions and related objects such as ultra--compact (UC) \hii regions, molecular clumps, and cores. We used these telescopes to study the \uchiir G43.89--0.78. The morphological study at arcminute scales using NVSS and \textit{Spitzer} data shows that this region is similar to those observed in the \textit{ bubble--like} structures revealed by \textit{Spitzer} observations. With this result, and including a physical characterization based on 3.6 cm data, we suggest G43.89--0.78 be classified as an \uchiir with Extended Emission because it meets the operational definition given in this paper comparing radio continuum data at 3.6 and 20~cm. For the ultra-compact component, we use VLA data to obtain physical parameters at 3.6~cm confirming this region as an \uchii region. Using ALMA observations, we detect the presence of a dense ($2.6\times10^7$ cm$^{-3}$) and small ($\sim$ 2.0\arcsec; 0.08 pc) molecular clump with a mass of 220 M$_{\odot}$ and average kinetic temperature of 21~K, located near to the \uchii region. In this clump, catalogued as G43.890--0.784, water masers also exist, possibly tracing a bipolar outflow. We discover in this vicinity two additional clumps which we label as G43.899--0.786 (T$_d$ = 50 K; M = 11 M$_{\odot}$), and G43.888--0.787 (T$_d$ = 50 K; M = 15 M$_{\odot}$).


title: Faraday-waves contact-line shear gradient induces streaming flow and tracers' self-organization: from rotating rings to spiral galaxy-like patterns

abstract: In this letter we experimentally demonstrate self-organization of small tracers under the action of longitudinal Faraday waves in a narrow container. We observe a steady current formation dividing the interface in small cells given by the symmetries of the Faraday wave. These streaming currents are rotating in each cell and their circulation increases with wave amplitude. This streaming flow drives the tracers to form patterns, whose shapes depend on the Faraday wave's amplitude: from low to high amplitudes we find dispersed tracers, a narrow rotating ring and a spiral galaxy-like pattern. We first describe the main pattern features, and characterize the wave and tracers' motion. We then show experimentally that the main source of the streaming flow comes from the time and spatial dependent shear at the wall contact line, created by the Faraday wave itself. We end by presenting a 2D model that considers the minimal ingredients present in the Faraday experiment, namely the stationary circulation, the stretching component due to the oscillatory wave and a steady converging field, which combined produce the observed self-organized patterns.


title: The large magnetocaloric effect and refrigerant capacity in nanocrystalline/ amorphous Gd$_3$Ni/Gd$_{65}$Ni$_{35}$ composite microwires

abstract: We report on a novel class of nanocrystalline/amorphous Gd$_3$Ni/Gd$_{65}$Ni$_{35}$ composite microwires, which was created directly by melt-extraction through controlled solidification. X-ray diffraction (XRD) and transmission electron microscopy (TEM) confirmed the formation of a biphase nanocrystalline/amorphous structure in these wires. Magnetic and magnetocaloric experiments indicate the large magnetic entropy change (-$\Delta$SM ~9.64 J/kg K) and the large refrigerant capacity (RC ~742.1 J/kg) around the Curie temperature of ~120 K for a field change of 5 T. These values are ~1.5 times larger relative to its bulk counterpart, and are superior to other candidate materials being considered for active magnetic refrigeration in the liquid nitrogen temperature range.


title: Bouquets of curves in surfaces

abstract: We characterise when a set of simple closed curves in an orientable surface forms a bouquet, in terms of relations between the corresponding Dehn twists.


title: A Survey of Algorithms for Geodesic Paths and Distances

abstract: Numerical computation of shortest paths or geodesics on curved domains, as well as the associated geodesic distance, arises in a broad range of applications across digital geometry processing, scientific computing, computer graphics, and computer vision. Relative to Euclidean distance computation, these tasks are complicated by the influence of curvature on the behavior of shortest paths, as well as the fact that the representation of the domain may itself be approximate. In spite of the difficulty of this problem, recent literature has developed a wide variety of sophisticated methods that enable rapid queries of geodesic information, even on relatively large models. This survey reviews the major categories of approaches to the computation of geodesic paths and distances, highlighting common themes and opportunities for future improvement.


title: The Timing Behavior of the Central Compact Object Pulsar 1E 1207.4-5209

abstract: We present 20 years of timing observations for 1E 1207.4-5209, the central compact object in supernova remnant PKS 1209-51/52, to follow up on our detection of an unexpected timing glitch in its spin-down. Using new XMM-Newton and NICER observations of 1E 1207.4-5209, we now find that the phase ephemeris can be well modelled by either two small glitches, or extreme timing noise. The implied magnitudes of the frequency glitches are Delta f/f = (9+\-2)E-10 and Delta f/f = (3.7+/-0.7)E-10, at epochs 2010.9 and 2014.4, respectively. The updated timing solutions also rule out our previous suggestion of a large glitch in the frequency derivative fdot. No other canonical pulsar with such a small spin-down rate (fdot = -1.2E-16 Hz/s) or surface dipole magnetic field strength (B_s = 9.8E10 G) has been observed to glitch; the glitch activity parameter of 1E 1207.4-5209 is larger than that of more energetic pulsars. Alternative parameterizations that do not involve glitches can fit the data, but they have timing residuals or a second frequency derivative fddot that are orders of magnitude larger than in pulsars with similar spin-down parameters. These timing properties of 1E 1207.4-5209 further motivate the leading theory of central compact objects, that an initial B-field of normal strength was buried in the neutron star crust by fallback of supernova ejecta, suppressing the surface dipole field. The slow reemergence of the buried field may be involved in triggering glitches or excess timing noise.


title: Treatment Effects with Targeting Instruments

abstract: Multivalued treatments are commonplace in applications. We explore the use of discrete-valued instruments to control for selection bias in this setting. Our discussion revolves around the concept of targeting instruments: which instruments target which treatments. It allows us to establish conditions under which counterfactual averages and treatment effects are point- or partially-identified for composite complier groups. We illustrate the usefulness of our framework by applying it to data from the Head Start Impact Study. Under a plausible positive selection assumption, we derive informative bounds that suggest less beneficial effects of Head Start expansions than the parametric estimates of Kline and Walters (2016).


title: Finite Cell Method for functionally graded materials based on V-models and homogenized microstructures

abstract: This paper proposes an extension of the finite cell method (FCM) to V-rep models, a novel geometric framework for volumetric representations. This combination of an embedded domain approach (FCM) and a new modeling framework (V-rep) forms the basis for an efficient and accurate simulation of mechanical artifacts, which are not only characterized by complex shapes but also by their non-standard interior structure. These objects gain more and more interest in the context of the new design opportunities opened by additive manufacturing, mainly when graded or micro-structured material is applied. Two different types of functionally graded materials (FGM) are considered: The first one, multi-material FGM, is described using the V-rep models' inherent property to assign different properties throughout the interior of a domain. The second, single-material FGM -- which is heterogeneously micro-structured -- characterizes the effective material behavior of representative volume elements by homogenization and performs large-scale simulations using the embedded domain approach.


title: Conformer-Kernel with Query Term Independence for Document Retrieval

abstract: The Transformer-Kernel (TK) model has demonstrated strong reranking performance on the TREC Deep Learning benchmark---and can be considered to be an efficient (but slightly less effective) alternative to BERT-based ranking models. In this work, we extend the TK architecture to the full retrieval setting by incorporating the query term independence assumption. Furthermore, to reduce the memory complexity of the Transformer layers with respect to the input sequence length, we propose a new Conformer layer. We show that the Conformer's GPU memory requirement scales linearly with input sequence length, making it a more viable option when ranking long documents. Finally, we demonstrate that incorporating explicit term matching signal into the model can be particularly useful in the full retrieval setting. We present preliminary results from our work in this paper.


title: Simulation of XANES spectroscopy and the calculation of total energies for N-heterocyclic carbenes on Au(111)

abstract: It has recently been demonstrated that N-heterocyclic carbenes (NHCs) form self-assembled monolayers (SAMs) on metal surfaces. Consequently, it is important to both characterize and understand their binding modes to fully exploit NHCs in functional surface systems. To assist with this effort, we have performed {\it first-principles} total energy calculations for NHCs on Au(111) and simulations of X-ray absorption near edge structure (XANES). The NHCs we have considered are N,N-dimethyl-, N,N-diethyl-, N,N-diisopropylbenzimidazolylidene ($^B$NHC$^X$, with X=Me, Et, and iPr, respectively) and the bis-$^B$NHC$^X$ complexes with Au derived from these molecules. We present a comprehensive analysis of the energetic stability of both the $^B$NHC$^X$ and the complexes on Au(111) and, for the former, examine the role of the wing group in determining the attachment geometry. Further structural characterization is performed by calculating the nitrogen K-edge X-ray absorption spectra. Our simulated XANES results give insight into (i) the relationship between the $^B$NHC$^X$/Au geometry and the N($1s$) $\rightarrow$ $\pi^\ast/\sigma^\ast$, pre-edge/near-edge, absorption intensities, and (ii) the contributions of the molecular deformation and molecule-surface electronic interaction to the XANES spectrum. Our simulations are compared with recent experimental results.


title: How different age groups responded to the COVID-19 pandemic in terms of mobility behaviors: a case study of the United States

abstract: The rapid spread of COVID-19 has affected thousands of people from different socio-demographic groups all over the country. A decisive step in preventing or slowing the outbreak is the use of mobility interventions, such as government stay-at-home orders. However, different socio-demographic groups might have different responses to these orders and regulations. In this paper, we attempt to fill the current gap in the literature by examining how different communities with different age groups performed social distancing by following orders such as the national emergency declaration on March 13, as well as how fast they started changing their behavior after the regulations were imposed. For this purpose, we calculated the behavior changes of people in different mobility metrics, such as percentage of people staying home during the study period (March, April, and May 2020), in different age groups in comparison to the days before the pandemic (January and February 2020), by utilizing anonymized and privacy-protected mobile device data. Our study indicates that senior communities outperformed younger communities in terms of their behavior change. Senior communities not only had a faster response to the outbreak in comparison to young communities, they also had better performance consistency during the pandemic.


title: Tunable transverse spin-motion coupling for quantum information processing

abstract: Laser-controlled entanglement between atomic qubits (`spins') and collective motion in trapped ion Coulomb crystals requires conditional momentum transfer from the laser. Since the spin-dependent force is derived from a spatial gradient in the spin-light interaction, this force is typically longitudinal -- parallel and proportional to the average laser $k$-vector (or two beams' $k$-vector difference), which constrains both the direction and relative magnitude of the accessible spin-motion coupling. Here, we show how momentum can also be transferred perpendicular to a single laser beam due to the gradient in its transverse profile. By controlling the transverse gradient at the position of the ion through beam shaping, the relative strength of the sidebands and carrier can be tuned to optimize the desired interaction and suppress undesired, off-resonant effects that can degrade gate fidelity. We also discuss how this effect may already be playing an unappreciated role in recent experiments.


title: Conditional calibration for false discovery rate control under dependence

abstract: We introduce a new class of methods for finite-sample false discovery rate (FDR) control in multiple testing problems with dependent test statistics where the dependence is fully or partially known. Our approach separately calibrates a data-dependent p-value rejection threshold for each hypothesis, relaxing or tightening the threshold as appropriate to target exact FDR control. In addition to our general framework we propose a concrete algorithm, the dependence-adjusted Benjamini-Hochberg (dBH) procedure, which adaptively thresholds the q-value for each hypothesis. Under positive regression dependence the dBH procedure uniformly dominates the standard BH procedure, and in general it uniformly dominates the Benjamini-Yekutieli (BY) procedure (also known as BH with log correction). Simulations and real data examples illustrate power gains over competing approaches to FDR control under dependence.


title: Subgroups of simple groups are as diverse as possible

abstract: For a finite group $G$, let $\sigma(G)$ be the number of subgroups of $G$ and $\sigma_\iota(G)$ the number of isomorphism types of subgroups of $G$. Let $L=L_r(p^e)$ denote a simple group of Lie type, rank $r$, over a field of order $p^e$ and characteristic $p$. If $r\neq 1$, $L\not\cong {^2 B_2}(2^{1+2m})$, then there are constants $c,d$, dependent on the Lie type, such that as $re$ grows $$p^{(c-o(1))r^4e^2}\leq\sigma_{\iota}(L_r(p^e))\leq\sigma(L_r(p^e)) \leq p^{(d+o(1))r^4e^2}.$$ For type $A$, $c=d=1/64$. For other classical groups $1/64\leq c\leq d\leq 1/4$. For exceptional and twisted groups $1/2^{100}\leq c\leq d\leq 1/4$. Furthermore, $$2^{(1/36-o(1))k^2)}\leq\sigma_{\iota}(\mathrm{Alt}_k)\leq \sigma(\mathrm{Alt}_k)\leq 24^{(1/6+o(1))k^2}.$$ For abelian and sporadic simple groups $G$, $\sigma_{\iota}(G),\sigma(G)\in O(1)$. In general these bounds are best possible amongst groups of the same orders. Thus with the exception of finite simple groups with bounded ranks and field degrees, the subgroups of finite simple groups are as diverse as possible.


title: Fundamental quantum limits in ellipsometry

abstract: We establish the ultimate limits that quantum theory imposes on the accuracy attainable in optical ellipsometry. We show that the standard quantum limit, as usual reached when the incident light is in a coherent state, can be surpassed with the use of appropriate squeezed states and, for tailored beams, even pushed to the ultimate Heisenberg limit.


title: Zero-Error Tracking for Autonomous Vehicles through Epsilon-Trajectory Generation

abstract: This paper presents a control method and trajectory planner for vehicles with first-order nonholonomic constraints that guarantee asymptotic convergence to a time-indexed trajectory. To overcome the nonholonomic constraint, a fixed point in front of the vehicle can be controlled to track a desired trajectory, albeit with a steady-state error. To eliminate steady state error, a sufficiently smooth trajectory is reformulated for the new reference point such that, when tracking the new trajectory, the vehicle asymptotically converges to the original trajectory. The resulting zero-error tracking law is demonstrated through a novel framework for creating time-indexed Clothoids. The Clothoids can be planned to pass through arbitrary waypoints using traditional methods yet result in trajectories that can be followed with zero steady-state error. The results of the control method and planner are illustrated in simulation wherein zero-error tracking is demonstrated.


title: Unlocking the Potential of Deep Counterfactual Value Networks

abstract: Deep counterfactual value networks combined with continual resolving provide a way to conduct depth-limited search in imperfect-information games. However, since their introduction in the DeepStack poker AI, deep counterfactual value networks have not seen widespread adoption. In this paper we introduce several improvements to deep counterfactual value networks, as well as counterfactual regret minimization, and analyze the effects of each change. We combined these improvements to create the poker AI Supremus. We show that while a reimplementation of DeepStack loses head-to-head against the strong benchmark agent Slumbot, Supremus successfully beats Slumbot by an extremely large margin and also achieves a lower exploitability than DeepStack against a local best response. Together, these results show that with our key improvements, deep counterfactual value networks can achieve state-of-the-art performance.


title: Magnetic Reconnection During the Post-Impulsive Phase of a Long-Duration Solar Flare: Bi-Directional Outflows as a Cause of Microwave and X-ray Bursts

abstract: Magnetic reconnection plays a crucial role in powering solar flares, production of energetic particles, and plasma heating. However, where the magnetic reconnections occur, how and where the released magnetic energy is transported, and how it is converted to other forms remain unclear. Here we report recurring bi-directional plasma outflows located within a large-scale plasma sheet observed in extreme ultraviolet emission and scattered white light during the post-impulsive gradual phase of the X8.2 solar flare on 2017 September 10. Each of the bi-directional outflows originates in the plasma sheet from a discrete site, identified as a magnetic reconnection site. These reconnection sites reside at very low altitudes ($< 180$ Mm, or 0.26 $R_{\odot}$) above the top of the flare arcade, a distance only $<3\%$ of the total length of a plasma sheet that extends to at least 10 $R_{\odot}$. Each arrival of sunward outflows at the looptop region appears to coincide with an impulsive microwave and X-ray burst dominated by a hot source (10-20 MK) at the looptop, which is immediately followed by a nonthermal microwave burst located in the loopleg region. We propose that the reconnection outflows transport the magnetic energy released at localized magnetic reconnection sites outward in the form of kinetic energy flux and/or electromagnetic Poynting flux. The sunward-directed energy flux induces particle acceleration and plasma heating in the post-flare arcades, observed as the hot and nonthermal flare emissions.


title: There are no $\sigma$-finite absolutely continuous invariant measures for multicritical circle maps

abstract: It is well-known that every multicritical circle map without periodic orbits admits a unique invariant Borel probability measure which is purely singular with respect to Lebesgue measure. Can such a map leave invariant an infinite, $\sigma$-finite invariant measure which is absolutely continuous with respect to Lebesgue measure? In this paper, using an old criterion due to Katznelson, we show that the answer to this question is no.


title: PanRep: Graph neural networks for extracting universal node embeddings in heterogeneous graphs

abstract: Learning unsupervised node embeddings facilitates several downstream tasks such as node classification and link prediction. A node embedding is universal if it is designed to be used by and benefit various downstream tasks. This work introduces PanRep, a graph neural network (GNN) model, for unsupervised learning of universal node representations for heterogenous graphs. PanRep consists of a GNN encoder that obtains node embeddings and four decoders, each capturing different topological and node feature properties. Abiding to these properties the novel unsupervised framework learns universal embeddings applicable to different downstream tasks. PanRep can be furthered fine-tuned to account for possible limited labels. In this operational setting PanRep is considered as a pretrained model for extracting node embeddings of heterogenous graph data. PanRep outperforms all unsupervised and certain supervised methods in node classification and link prediction, especially when the labeled data for the supervised methods is small. PanRep-FT (with fine-tuning) outperforms all other supervised approaches, which corroborates the merits of pretraining models. Finally, we apply PanRep-FT for discovering novel drugs for Covid-19. We showcase the advantage of universal embeddings in drug repurposing and identify several drugs used in clinical trials as possible drug candidates.


title: Density Functional Theory-based Quantum Mechanics/Coarse-grained Molecular Mechanics: Theory and Implementation

abstract: Quantum mechanics/molecular mechanics (QM/MM) is a standard computational tool for describing chemical reactivity in systems with many degrees of freedom, including polymers, enzymes, and reacting molecules in complex solvents. However, QM/MM is less suitable for systems with complex MM dynamics due to associated long relaxation times, the high computational cost of QM energy evaluations, and expensive long-range electrostatics. Recently, a systematic coarse-graining of the MM part was proposed to overcome these QM/MM limitations in the form of the quantum mechanics/coarse-grained molecular mechanics (QM/CG-MM) approach. Herein, we recast QM/CG-MM in the density functional theory formalism and, by employing the force-matching variational principle, access the method performance for two model systems: QM CCl4 in the MM CCl4 liquid and the reaction of tert-butyl hypochlorite with the benzyl radical in the MM CCl4 solvent. We find that DFT-QM/CG-MM accurately reproduces DFT-QM/MM radial distribution functions and 3-body correlations between QM and CG-MM subsystems. The free energy profile of the reaction is also described well, with an error < 1-2 kcal/mol. DFT-QM/CG-MM is a general, systematic, and computationally efficient approach to include chemical reactivity in coarse-grained molecular models.


title: On $q$-analogs of zeta functions associated with a pair of $q$-analogs of Bernoulli numbers and polynomials

abstract: In this paper, we use two different approaches to introduce $q$-analogs of Riemann's zeta function and prove that their values at even integers are related to the $q$-Bernoulli and $q$ Euler's numbers introduced by Ismail and Mansour [Analysis and Applications, {\bf{17}}, 6, 2019, 853--895].