title: Bose-Einstein condensation for two dimensional bosons in the Gross-Pitaevskii regime

abstract: We consider systems of N bosons trapped on the two-dimensional unit torus, in the Gross-Pitaevskii regime, where the scattering length of the repulsive interaction is exponentially small in the number of particles. We show that low-energy states exhibit complete Bose-Einstein condensation, with almost optimal bounds on the number of orthogonal excitations.

categories: math-ph cond-mat.quant-gas math.AP math.MP


title: Wetting of a Solid Surface by Active Matter

abstract: A lattice model is used to study repulsive active particles at a planar surface. A rejection-free Kinetic Monte Carlo method is employed to characterize the wetting behaviour. The model predicts a mobility induced phase separation of active particles, and the bulk coexistence of dense liquid-like and dilute vapour-like steady states is determined. An "ensemble", with a varying number of particles, analogous to a grand canonical ensemble in equilibrium, is introduced. The formation and growth of the liquid film between the solid surface and the vapour phase is investigated. For all of the activities considered, the thickness of the adsorbed film exhibits a diverging behaviour as the system is brought towards coexistence from the vapour side, suggesting a complete wetting scenario along the full coexistence curve.


title: A 2-component Camassa-Holm equation, Euler-Bernoulli Beam Problem and Non-Commutative Continued Fractions

abstract: A new approach to the Euler-Bernoulli beam based on an inhomogeneous matrix string problem is presented. Three ramifications of the approach are developed: (1) motivated by an analogy with the Camassa-Holm equation a class of isospectral deformations of the beam problem is formulated; (2) a reformulation of the matrix string problem in terms of a certain compact operator is used to obtain basic spectral properties of the inhomogeneous matrix string problem with Dirichlet boundary conditions; (3) the inverse problem is solved for the special case of a discrete Euler-Bernoulli beam. The solution involves a non-commutative generalization of Stieltjes' continued fractions, leading to the inverse formulas expressed in terms of ratios of Hankel-like determinants.


title: Predictive risk estimation for the Expectation Maximization algorithm with Poisson data

abstract: In this work, we introduce a novel estimator of the predictive risk with Poisson data, when the loss function is the Kullback-Leibler divergence, in order to define a regularization parameter's choice rule for the Expectation Maximization (EM) algorithm. To this aim, we prove a Poisson counterpart of the Stein's Lemma for Gaussian variables, and from this result we derive the proposed estimator showing its analogies with the well-known Stein's Unbiased Risk Estimator valid for a quadratic loss. We prove that the proposed estimator is asymptotically unbiased with increasing number of measured counts, under certain mild conditions on the regularization method. We show that these conditions are satisfied by the EM algorithm and then we apply this estimator to select its optimal reconstruction. We present some numerical tests in the case of image deconvolution, comparing the performances of the proposed estimator with other methods available in the literature, both in the inverse crime and non-inverse crime setting.


title: Solsticial Hadley Cell ascending edge theory from supercriticality

abstract: How far the Hadley circulation's ascending branch extends into the summer hemisphere is a fundamental but incompletely understood characteristic of Earth's climate. Here, we present a predictive, analytical theory for this ascending edge latitude based on the extent of supercritical forcing. Supercriticality sets the minimum extent of a large-scale circulation based on the angular momentum and absolute vorticity distributions of the hypothetical state were the circulation absent. We explicitly simulate this latitude-by-latitude radiative-convective equilibrium (RCE) state. Its depth-averaged temperature profile is suitably captured by a simple analytical approximation that increases linearly with $\sin\varphi$, where $\varphi$ is latitude, from the winter to the summer pole. This, in turn, yields a one-third power-law scaling of the supercritical forcing extent with the thermal Rossby number. In moist and dry idealized GCM simulations under solsticial forcing performed with a wide range of planetary rotation rates, the ascending edge latitudes largely behave according to this scaling.


title: Vital role of anisotropy in cubic chiral skyrmion hosts

abstract: The impact of magnetic anisotropy on the skyrmion lattice (SkL) state in cubic chiral magnets has been overlooked for long, partly because a semi-quantitative description of the thermodynamically stable SkL phase pocket forming near the Curie temperature could be achieved without invoking anisotropy effects. However, there has been a range of phenomena reported recently in these materials, such as the formation of low-temperature tilted conical and SkL states as well as temperature-induced transformations of lattice geometry in metastable SkL states, where anisotropy was suspected to play a key role. To settle this issue on experimental basis, we quantified the cubic anisotropy in a series of CoZnMn-type cubic chiral magnets. We found that the strength of anisotropy is highly enhanced towards low temperatures in all the compounds, moreover, not only the magnitude but also the character of cubic anisotropy drastically varies upon changing the Co/Mn ratio. We correlate these changes with temperature- and composition-induced variations of the helical modulation vectors, the anharmonicity and structural rearrangements of the metastable SkLs and the spin relaxation rates. Similar systematic studies on magnetic anisotropy may not only pave the way for a quantitative and unified description of the stable and metastable modulated spin textures in cubic chiral magnets but would also help exploring further topological spin textures in this large class of skyrmion hosts.


title: Isotope Shifts in the Metastable a$^5$F and Excited y$^5$G$^\circ$ Terms of Atomic Titanium

abstract: We measure and analyze the isotope shifts the multiplet of transitions between the metastable a$^5$F and excited y$^5$G$^\circ$ terms of neutral titanium by probing a titanium vapor in a hollow cathode lamp using saturated absorption spectroscopy. We resolve the five $J\to J+1$ and the four $J\to J$ transitions within the multiplet for each of the the three $I=0$ stable isotopes ($^{46}$Ti, $^{48}$Ti, and $^{50}$Ti). The isotope shifts on these transitions allow us to determine the isotope-dependent variation in the fine-structure splitting of the a$^5$F and y$^5$G$^\circ$ levels themselves. Combined with existing knowledge of the nuclear charge radii of titanium nuclei, we derive the specific mass and field shifts, which arise from correlated electronic motion and electronic density at the nucleus respectively, and further observe a strong $J$-dependent variation in each. Our results yield insight into the electronic and nuclear structure of transition metal atoms like titanium, and also characterize optical transitions that may allow for optical manipulation of ultracold gases of transition metal species.


title: Unified mapmaking for an anisotropic stochastic gravitational wave background

abstract: A stochastic gravitational wave background (SGWB), created by the superposition of signals from unresolved astrophysical sources, may be detected in the next few years. Estimating the variation of intensity across the sky can, therefore, play a key role in improving our understanding of astrophysical models. Skymaps have been produced for all the data-taking runs of the advanced ground-based interferometric detectors. While these maps are being produced in pixel and SpH bases, to probe, respectively, localized and diffuse astrophysical and cosmological sources, with algorithms that employ cross-correlation as the common strategy, the underlying algebra and numerical implementation remain different. As a consequence, there was a need for producing skymaps in both bases in those analyses. We show that these manifestly redundant methods could indeed be unified to a single analysis that can probe very different scales and demonstrate it by applying them on real data. We first develop the algebra to show that the results in two different bases are easily transformable. We then incorporate both the schemes in the now-standard analysis pipeline for anisotropic SGWB, PyStoch. This will enable SGWB anisotropy searches in SpH basis also to take full advantage of integrated \hpx tools and makes it computationally feasible to perform the search in every frequency bin. We, however, follow a different approach for direct estimation of the SpH moments. We show that the results obtained from these different methods match very well; the differences are less than $0.1$\% for the SpH moments and less than $0.01$\% for the Fisher information matrices. Thus we conclude that a single skymap will be sufficient to describe the anisotropies in a stochastic background. The multiple capabilities of PyStoch will be useful for estimating and constraining various measures that characterize an anisotropic background.


title: Transformers for One-Shot Visual Imitation

abstract: Humans are able to seamlessly visually imitate others, by inferring their intentions and using past experience to achieve the same end goal. In other words, we can parse complex semantic knowledge from raw video and efficiently translate that into concrete motor control. Is it possible to give a robot this same capability? Prior research in robot imitation learning has created agents which can acquire diverse skills from expert human operators. However, expanding these techniques to work with a single positive example during test time is still an open challenge. Apart from control, the difficulty stems from mismatches between the demonstrator and robot domains. For example, objects may be placed in different locations (e.g. kitchen layouts are different in every house). Additionally, the demonstration may come from an agent with different morphology and physical appearance (e.g. human), so one-to-one action correspondences are not available. This paper investigates techniques which allow robots to partially bridge these domain gaps, using their past experience. A neural network is trained to mimic ground truth robot actions given context video from another agent, and must generalize to unseen task instances when prompted with new videos during test time. We hypothesize that our policy representations must be both context driven and dynamics aware in order to perform these tasks. These assumptions are baked into the neural network using the Transformers attention mechanism and a self-supervised inverse dynamics loss. Finally, we experimentally determine that our method accomplishes a $\sim 2$x improvement in terms of task success rate over prior baselines in a suite of one-shot manipulation tasks.


title: Higher-Derivative Supergravity, Wrapped M5-branes, and Theories of Class $\mathcal{R}$

abstract: We study the interplay between four-derivative 4d gauged supergravity, holography, wrapped M5-branes, and theories of class $\mathcal{R}$. Using results from Chern-Simons theory on hyperbolic three-manifolds and the 3d-3d correspondence we are able to constrain the two independent coefficients in the four-derivative supergravity Lagrangian. This in turn allows us to calculate the subleading terms in the large-$N$ expansion of supersymmetric partition functions for an infinite class of three-dimensional $\mathcal{N}=2$ SCFTs of class $\mathcal{R}$. We also determine the leading correction to the Bekenstein-Hawking entropy of asymptotically AdS$_4$ black holes arising from wrapped M5-branes. In addition, we propose and test some conjectures about the perturbative partition function of Chern-Simons theory with complexified ADE gauge groups on closed hyperbolic three-manifolds.


title: Bianchi II and VII$_{h=0}$ Models Revisited Via The Euclidean-Signature Semi Classical Method

abstract: We apply in a novel fashion a modified semi-classical method to the Bianchi II and VII$_{h=0}$ models when a cosmological constant, aligned electromagnetic field and stiff matter are present. Additionally we study the non-commutative quantum Bianchi II models when an aligned electromagnetic field is included. Through the use of the Euclidean-signature semi classical method we find a plethora of new solutions to these model's corresponding Lorentzian signature Wheeler DeWitt equations which we can interpret qualitatively. These new solutions for the aforementioned models involving matter sources reveal some potentially interesting effects that should be chronicled as possible phenomena that a toy model of quantum gravity can induce on the evolution of a quantum universe. Furthermore we find 'excited' states which behave differently from the 'excited' states of the Bianchi IX and Taub models that were previously uncovered using this method. By comparing and contrasting the 'excited' states given by these models we help facilitate a better understating of what constitutes an 'excited' state solution of the Wheeler Dewitt equation. Our results further show the utility of the Euclidean-signature semi classical method for tackling Lorentzian signature problems without having to invoke a Wick rotation. This feature of not needing to apply a Wick rotation makes this method potentially very useful for tackling a variety of problems in bosonic relativistic field theory and quantum gravity.


title: A survey on practical adversarial examples for malware classifiers

abstract: Machine learning based solutions have been very helpful in solving problems that deal with immense amounts of data, such as malware detection and classification. However, deep neural networks have been found to be vulnerable to adversarial examples, or inputs that have been purposefully perturbed to result in an incorrect label. Researchers have shown that this vulnerability can be exploited to create evasive malware samples. However, many proposed attacks do not generate an executable and instead generate a feature vector. To fully understand the impact of adversarial examples on malware detection, we review practical attacks against malware classifiers that generate executable adversarial malware examples. We also discuss current challenges in this area of research, as well as suggestions for improvement and future research directions.


title: Beyond the Fr\"ohlich Hamiltonian: Path integral treatment of large polarons in anharmonic solids

abstract: The properties of an electron in a typical solid are modified by the interaction with the crystal ions, leading to the formation of a quasiparticle: the polaron. Such polarons are often described using the Fr\"ohlich Hamiltonian, which assumes the underlying lattice phonons to be harmonic. However, this approximation is invalid in several interesting materials, including the recently discovered high-pressure hydrides which superconduct at temperatures above $200$K. In this paper, we show that Fr\"ohlich theory can be extended to eliminate this problem. We derive four additional terms in the Fr\"ohlich Hamiltonian to account for anharmonicity up to third order. We calculate the energy and effective mass of the new polaron, using both perturbation theory and Feynman's path integral formalism. It is shown that the anharmonic terms lead to significant additional trapping of the electron. The derived Hamiltonian is well-suited for analytical calculations, due to its simplicity and since the number of model parameters is low. Since it is a direct extension of the Fr\"ohlich Hamiltonian, it can readily be used to investigate the effect of anharmonicity on other polaron properties, such as the optical conductivity and the formation of bipolarons.


title: Small extensions of abelian ordered groups

abstract: Let $\Gamma$ be a totally ordered group. We use Hahn's embedding theorem to construct a totally ordered set $\Gamma\subset \Gamma_{\operatorname{sme}}$ which classifies small extensions of $\Gamma$. This small-extensions closure $\Gamma_{\operatorname{sme}}$ is complete and plays a crucial role in the description of equivalence classes of valuations on the polynomial ring $K[x]$ over a field $K$.


title: The Vulnerability of the Neural Networks Against Adversarial Examples in Deep Learning Algorithms

abstract: With further development in the fields of computer vision, network security, natural language processing and so on so forth, deep learning technology gradually exposed certain security risks. The existing deep learning algorithms cannot effectively describe the essential characteristics of data, making the algorithm unable to give the correct result in the face of malicious input. Based on current security threats faced by deep learning, this paper introduces the problem of adversarial examples in deep learning, sorts out the existing attack and defense methods of the black box and white box, and classifies them. It briefly describes the application of some adversarial examples in different scenarios in recent years, compares several defense technologies of adversarial examples, and finally summarizes the problems in this research field and prospects for its future development. This paper introduces the common white box attack methods in detail, and further compares the similarities and differences between the attack of the black and white box. Correspondingly, the author also introduces the defense methods, and analyzes the performance of these methods against the black and white box attack.


title: Dyonic Objects and Tensor Network Representation

abstract: Motivated by particle phyiscs results, we investigate certain dyonic solutions in arbitrary dimensions. Concretely, we study the stringy constructions of such objects from concrete compactifications. Then we elaborate their tensor network realizations using multistate particle formalism.


title: The Impact of Text Presentation on Translator Performance

abstract: Widely used computer-aided translation (CAT) tools divide documents into segments such as sentences and arrange them in a side-by-side, spreadsheet-like view. We present the first controlled evaluation of these design choices on translator performance, measuring speed and accuracy in three experimental text processing tasks. We find significant evidence that sentence-by-sentence presentation enables faster text reproduction and within-sentence error identification compared to unsegmented text, and that a top-and-bottom arrangement of source and target sentences enables faster text reproduction compared to a side-by-side arrangement. For revision, on the other hand, our results suggest that presenting unsegmented text results in the highest accuracy and time efficiency. Our findings have direct implications for best practices in designing CAT tools.


title: Probing episodic accretion with chemistry: CALYPSO observations of IRAM 04191+1522. Results from the CALYPSO IRAM-PdBI survey

abstract: Context. The process of mass accretion in the earliest phases of star formation is still not fully understood: Does the accretion rate smoothly decline with the age of the protostar or are there short, intermittent accretion bursts? Aims. We aim to explore whether or not the observed C$^{18}$O and N$_2$H$^+$ emission pattern towards the VeLLO IRAM 04191+1522 can be understood in the framework of a scenario where the emission is chemically tracing a past accretion burst. Methods.We used high-angular-resolution Plateau de Bure Interferometer (PdBI) observations of C$^{18}$O and N$_2$H$^+$ towards IRAM 04191+1522 that were obtained as part of the CALYPSO IRAM Large Program. We model these observations using a chemical code with a time-dependent physical structure coupled with a radiative transfer module, where we allow for variations in the source luminosity. Results. We find that the N$_2$H$^+$ line emission shows a central hole, while the C$^{18}$O emission is compact. The morphology of these two lines cannot be reproduced with a constant luminosity model based on the present-day internal luminosity (0.08 L$_{sun}$). However, the N$_2$H$^+$ peaks are consistent with a constant-luminosity model of 12 L$_{sun}$. Using a model with time-dependent temperature and density profiles, we show that the observed N$_2$H$^+$ peak emission could indeed be caused by a past accretion burst. Such a burst should have occurred a couple of hundred years ago. Conclusions. We suggest that an accretion burst occurred in IRAM 04191+1522 in the recent past. If such bursts are common and sufficiently long in VeLLOs, they could lead to higher accretion onto the central object than their luminosity suggests. For IRAM 04191 in particular, our results yield an estimated final mass of 0.2 - 0.25 M$_{sun}$ by the end of the Class 0 phase, which would make this object a low-mass star rather than a brown dwarf.


title: GRCNN: Graph Recognition Convolutional Neural Network for Synthesizing Programs from Flow Charts

abstract: Program synthesis is the task to automatically generate programs based on user specification. In this paper, we present a framework that synthesizes programs from flow charts that serve as accurate and intuitive specifications. In order doing so, we propose a deep neural network called GRCNN that recognizes graph structure from its image. GRCNN is trained end-to-end, which can predict edge and node information of the flow chart simultaneously. Experiments show that the accuracy rate to synthesize a program is 66.4%, and the accuracy rates to recognize edge and nodes are 94.1% and 67.9%, respectively. On average, it takes about 60 milliseconds to synthesize a program.


title: T-Weyl calculus

abstract: Let $\left( W,\sigma \right) $ be a symplectic vector space and let $% T:W\rightarrow W$ be a linear map that satisfies a certain condition of non-degeneracy. We define the Schur multiplier $\omega _{\sigma ,T}$ on $W$. To this multiplier we associate a $\omega _{\sigma ,T}$-representation and and we build the $T$-Weyl calculus, $\mathrm{Op}_{\sigma ,T}$, whose properties are are systematically studied further.


title: The Photospheric Temperatures of Betelgeuse during the Great Dimming of 2019/2020: No New Dust Required

abstract: The processes that shape the extended atmospheres of red supergiants (RSGs), heat their chromospheres, create molecular reservoirs, drive mass loss, and create dust remain poorly understood. Betelgeuse's V-band "Great Dimming" event of 2019 September /2020 February and its subsequent rapid brightening provides a rare opportunity to study these phenomena. Two different explanations have emerged to explain the dimming; new dust appeared in our line of sight attenuating the photospheric light, or a large portion of the photosphere had cooled. Here we present five years of Wing three-filter (A, B, and C band) TiO and near-IR photometry obtained at the Wasatonic Observatory. These reveal that parts of the photosphere had a mean effective temperature $(T_{\rm eff}$) significantly lower than that found by (Levesque & Massey 2020). Synthetic photometry from MARCS -model photospheres and spectra reveal that the V band, TiO index, and C-band photometry, and previously reported 4000-6800 Angstrom spectra can be quantitatively reproduced if there are multiple photospheric components, as hinted at by VLT-SPHERE images (Montarges et al. 2020). If the cooler component has $\Delta T_{\rm eff} \ge 250$ K cooler than 3650 K, then no new dust is required to explain the available empirical constraints. A coincidence of the dominant short- ($\sim 430$ day) and long-period ($\sim 5.8$ yr) V-band variations occurred near the time of deep minimum (Guinan et al. 2019). This is in tandem with the strong correlation of V mag and photospheric radial velocities, recently reported by Dupree et al. (2020b). These suggest that the cooling of a large fraction of the visible star has a dynamic origin related to the photospheric motions, perhaps arising from pulsation or large-scale convective motions.


title: X-ray redshifts for obscured AGN: a case study in the J1030 deep field

abstract: We present a procedure to constrain the redshifts of obscured ($N_H > 10^{22}$ cm$^{-2}$) Active Galactic Nuclei (AGN) based on low-count statistics X-ray spectra, which can be adopted when photometric and/or spectroscopic redshifts are unavailable or difficult to obtain. We selected a sample of 54 obscured AGN candidates on the basis of their X-ray hardness ratio, $HR>-0.1$, in the Chandra deep field ($\sim$479 ks, 335 arcmin$^2$) around the $z=6.3$ QSO SDSS J1030+0524. The sample has a median value of $\approx80$ net counts in the 0.5-7 keV energy band. We estimate reliable X-ray redshift solutions taking advantage of the main features in obscured AGN spectra, like the Fe 6.4 keV K$\mathrm{\alpha}$ emission line, the 7.1 keV Fe absorption edge and the photoelectric absorption cut-off. The significance of such features is investigated through spectral simulations, and the derived X-ray redshift solutions are then compared with photometric redshifts. Both photometric and X-ray redshifts are derived for 33 sources. When multiple solutions are derived by any method, we find that combining the redshift solutions of the two techniques improves the rms by a factor of two. Using our redshift estimates ($0.1\lesssim z \lesssim 4$), we derived absorbing column densities in the range $\sim 10^{22}-10^{24}$ cm$^{-2}$ and absorption-corrected, 2-10 keV rest-frame luminosities between $\sim 10^{42}$ and $10^{45}$ erg s$^{-1}$, with median values of $N_H = 1.7 \times 10^{23}$ cm$^{-2}$ and $L_{\mathrm{2-10\, keV}} = 8.3\times10^{43}$ erg s$^{-1}$, respectively. Our results suggest that the adopted procedure can be applied to current and future X-ray surveys, for sources detected only in the X-rays or that have uncertain photometric or single-line spectroscopic redshifts.


title: Dynamics of market states and risk assessment

abstract: Previous research explored various conditions of financial markets based on the similarity of correlation structures and classified as market states. We introduce modifications to previous selection criteria for these market states, mainly due to increased attention to the transition matrix between the states. Clustering and thus market states are fixed by the optimization of two parameters -- number of clusters and noise suppression, but in similar conditions, we give preference to the clustering which avoids large jumps in the transition matrix. We found statistically significant results applying this model to the SP 500 and Nikkei 225 markets for the pre-COVID-19 pandemic era (2006-2019). Retaining the epoch length of 20 trading days but reducing the shift of the epoch to a single trading day we are led to the concept of a trajectory of the market in the space of correlation matrices. We may visualize these states after dimensional scaling to two or three dimensions. This approach, using dynamics, improves the options of risk assessment, opens the door to dynamical treatments of markets (e.g. hedging), and shows noise suppression in a new light.


title: Dirichlet Pruning for Neural Network Compression

abstract: We introduce Dirichlet pruning, a novel post-processing technique to transform a large neural network model into a compressed one. Dirichlet pruning is a form of structured pruning that assigns the Dirichlet distribution over each layer's channels in convolutional layers (or neurons in fully-connected layers) and estimates the parameters of the distribution over these units using variational inference. The learned distribution allows us to remove unimportant units, resulting in a compact architecture containing only crucial features for a task at hand. The number of newly introduced Dirichlet parameters is only linear in the number of channels, which allows for rapid training, requiring as little as one epoch to converge. We perform extensive experiments, in particular on larger architectures such as VGG and ResNet (45% and 58% compression rate, respectively) where our method achieves the state-of-the-art compression performance and provides interpretable features as a by-product.


title: Thermal noise and mechanical loss of SiO$_2$/Ta$_2$O$_5$ optical coatings at cryogenic temperatures

abstract: Mechanical loss of dielectric mirror coatings sets fundamental limits for both gravitational wave detectors and cavity-stabilized optical local oscillators for atomic clocks. Two approaches are used to determine the mechanical loss: ringdown measurements of the coating quality factor and direct measurement of the coating thermal noise. Here we report a systematic study of the mirror thermal noise from room temperature to 4 K by operating reference cavities at these temperatures. The directly measured thermal noise is used to extract the corresponding mechanical loss for SiO$_2$/Ta$_2$O$_5$ coatings, which are compared with previously reported values.


title: Physics-constrained Deep Learning of Multi-zone Building Thermal Dynamics

abstract: We present a physics-constrained control-oriented deep learning method for modeling building thermal dynamics. The proposed method is based on the systematic encoding of physics-based prior knowledge into a structured recurrent neural architecture. Specifically, our method incorporates structural priors from traditional physics-based building modeling into the neural network thermal dynamics model structure. Further, we leverage penalty methods to provide inequality constraints, thereby bounding predictions within physically realistic and safe operating ranges. Observing that stable eigenvalues accurately characterize the dissipativeness of the system, we additionally use a constrained matrix parameterization based on the Perron-Frobenius theorem to bound the dominant eigenvalues of the building thermal model parameter matrices. We demonstrate the proposed data-driven modeling approach's effectiveness and physical interpretability on a dataset obtained from a real-world office building with 20 thermal zones. Using only 10 days' measurements for training, we demonstrate generalization over 20 consecutive days, significantly improving the accuracy compared to prior state-of-the-art results reported in the literature.


title: Maximum sampled conditional likelihood for informative subsampling

abstract: Subsampling is a computationally effective approach to extract information from massive data sets when computing resources are limited. After a subsample is taken from the full data, most available methods use an inverse probability weighted (IPW) objective function to estimate the model parameters. The IPW estimator does not fully utilize the information in the selected subsample. In this paper, we propose to use the maximum sampled conditional likelihood estimator (MSCLE) based on the sampled data. We established the asymptotic normality of the MSCLE and prove that its asymptotic variance covariance matrix is the smallest among a class of asymptotically unbiased estimators, including the IPW estimator. We further discuss the asymptotic results with the L-optimal subsampling probabilities and illustrate the estimation procedure with generalized linear models. Numerical experiments are provided to evaluate the practical performance of the proposed method.

categories: math.ST cs.IT cs.LG math.IT stat.CO stat.ME stat.TH


title: Linear Dilation-Erosion Perceptron for Binary Classification

abstract: In this work, we briefly revise the reduced dilation-erosion perceptron (r-DEP) models for binary classification tasks. Then, we present the so-called linear dilation-erosion perceptron (l-DEP), in which a linear transformation is applied before the application of the morphological operators. Furthermore, we propose to train the l-DEP classifier by minimizing a regularized hinge-loss function subject to concave-convex restrictions. A simple example is given for illustrative purposes.


title: The Non Relativistic Interiors of Ultra-Relativistic Explosions: Extension to the Blandford McKee Solutions

abstract: The hydrodynamics of an ultrarelativistic flow, enclosed by a strong shock wave, are described by the well known Blandford-McKee solutions in spherical geometry. These solutions, however, become inaccurate at a distance $\sim R/2$ behind the shock wave, where $R$ is the shock radius, as the flow approaches Newtonian velocities. In this work we find a new self-similar solution which is an extension to the Blandford-McKee solutions, and which describes the interior part of the blast wave, where the flow reaches mildly relativistic to Newtonian velocities. We find that the velocity profile of the internal part of the flow does not depend on the value of the shock Lorentz factor, $\Gamma$, and is accurate from $r=0$ down to a distance of $R/\Gamma^2$ behind the shock. Despite the fact that the shock wave is in causal contact with the entire flow behind it, a singular point appears in the equations. Nevertheless, the solution is not required to pass through the singular point: for ambient density that decreases slowly enough, $\rho \propto r^{-k}$ with $k<\frac{1}{2}(5-\sqrt{10})\cong0.92$, a secondary shock wave forms with an inflow towards the origin.


title: Solving high-dimensional parameter inference: marginal posterior densities & Moment Networks

abstract: High-dimensional probability density estimation for inference suffers from the "curse of dimensionality". For many physical inference problems, the full posterior distribution is unwieldy and seldom used in practice. Instead, we propose direct estimation of lower-dimensional marginal distributions, bypassing high-dimensional density estimation or high-dimensional Markov chain Monte Carlo (MCMC) sampling. By evaluating the two-dimensional marginal posteriors we can unveil the full-dimensional parameter covariance structure. We additionally propose constructing a simple hierarchy of fast neural regression models, called Moment Networks, that compute increasing moments of any desired lower-dimensional marginal posterior density; these reproduce exact results from analytic posteriors and those obtained from Masked Autoregressive Flows. We demonstrate marginal posterior density estimation using high-dimensional LIGO-like gravitational wave time series and describe applications for problems of fundamental cosmology.


title: Neural networks as optimal estimators to marginalize over baryonic effects

abstract: Many different studies have shown that a wealth of cosmological information resides on small, non-linear scales. Unfortunately, there are two challenges to overcome to utilize that information. First, we do not know the optimal estimator that will allow us to retrieve the maximum information. Second, baryonic effects impact that regime significantly and in a poorly understood manner. Ideally, we would like to use an estimator that extracts the maximum cosmological information while marginalizing over baryonic effects. In this work we show that neural networks can achieve that. We made use of data where the maximum amount of cosmological information is known: power spectra and 2D Gaussian density fields. We also contaminate the data with simplified baryonic effects and train neural networks to predict the value of the cosmological parameters. For this data, we show that neural networks can 1) extract the maximum available cosmological information, 2) marginalize over baryonic effects, and 3) extract cosmological information that is buried in the regime dominated by baryonic physics. We also show that neural networks learn the priors of the data they are trained on. We conclude that a promising strategy to maximize the scientific return of cosmological experiments is to train neural networks on state-of-the-art numerical simulations with different strengths and implementations of baryonic effects.


title: Constraints on Axions from Cosmic Distance Measurements

abstract: Axion couplings to photons could induce photon-axion conversion in the presence of magnetic fields in the Universe. This conversion could impact various cosmic distance measurements, such as luminosity distances to type Ia supernovae and angular distances to galaxy clusters, in different ways. In this paper we consider different combinations of the most up-to-date distance measurements to constrain the axion-photon coupling. Employing the conservative cell magnetic field model for the magnetic fields in the intergalactic medium (IGM) and ignoring the conversion in the intracluster medium (ICM), we find the upper bounds on axion-photon couplings to be around $5 \times 10^{-12}$ (nG/$B$) $\sqrt{\mathrm{Mpc/s}}$ GeV$^{-1}$ for axion masses $m_a$ below $10^{-13}$ eV, where $B$ is the strength of the IGM magnetic field, and $s$ is the comoving size of the magnetic domains. When including the conversion in the ICM, the upper bound is lowered and could reach $5 \times 10^{-13}\, $GeV$^{-1}$ for $m_a < 5 \times 10^{-12}$ eV. While this stronger bound depends on the ICM modeling, it is independent of the strength of the IGM magnetic field, for which there is no direct evidence yet. These constraints could be placed on firmer footing with an enhanced understanding and control of the astrophysical uncertainties associated with the IGM and ICM. All the bounds are determined by the shape of the Hubble rate as a function of redshift reconstructable from various distance measurements, and insensitive to today's Hubble rate, of which there is a tension between early and late cosmological measurements. As an appendix, we discuss the model building challenges of the use of photon-axion conversion to make type Ia supernovae brighter to alleviate the Hubble problem/crisis.


title: ${\cal N}=2$ Supersymmetric Partially Massless Fields and Non-Unitary Superconformal Representations

abstract: We find and classify the simplest ${\cal N}=2$ SUSY multiplets on AdS$_4$ which contain partially massless fields. We do this by studying representations of the ${\cal N}=2$, $d=3$ superconformal algebra of the boundary, including new shortening conditions that arise in the non-unitary regime. Unlike the ${\cal N}=1$ case, the simplest ${\cal N}=2$ multiplet containing a partially massless spin-2 is short, containing several exotic fields. More generally, we argue that ${\cal N}=2$ supersymmetry allows for short multiplets that contain partially massless spin-$s$ particles of depth $t=s-2$.


title: Heavy Axion Opportunities at the DUNE Near Detector

abstract: While the QCD axion is often considered to be necessarily light ($\lesssim$ eV), recent work has opened a viable and interesting parameter space for heavy axions, which solve both the Strong CP and the axion Quality Problems. These well-motivated heavy axions, as well as the generic axion-like-particles, call for explorations in the GeV mass realm at collider and beam dump environments. The primary upcoming neutrino experiment, Deep Underground Neutrino Experiment (DUNE), is simultaneously also a powerful beam dump experiment, enabled by its multipurpose Near Detector (ND) complex. In this study, we show with detailed analyses that the DUNE ND has a unique sensitivity to heavy axions for masses between $20$ MeV and $2$ GeV, complementary to other future experiments.


title: The influence of infall on the properties of protoplanetary discs

abstract: We perform a population synthesis of protoplanetary discs including infall with a total of $50\,000$ simulations using a 1D vertically integrated viscous evolution code, studying a large parameter space in final stellar mass. Initial conditions and infall locations are chosen based on the results from a radiation-hydrodynamic population synthesis of circumstellar discs. We also consider a different infall prescription based on a magnetohydrodynamic (MHD) collapse simulation in order to assess the influence of magnetic fields on disc formation. The duration of the infall phase is chosen to produce a stellar mass distribution in agreement with the observationally determined stellar initial mass function. We find that protoplanetary discs are very massive early in their lives. When averaged over the entire stellar population, the discs have masses of $\sim 0.3$ and $0.1\,\mathrm{M_\odot}$ for systems based on hydrodynamic or MHD initial conditions, respectively. In systems with final stellar mass $\sim 1\,\mathrm{M_\odot}$, we find disc masses of $\sim 0.7\,\mathrm{M_\odot}$ for the `hydro' case and $\sim 0.2\,\mathrm{M_\odot}$ for the `MHD' case at the end of the infall phase. Furthermore, the inferred total disc lifetimes are long, $\approx 5-7\,\mathrm{Myr}$ on average, despite our choice of a high value of $10^{-2}$ for the background viscosity $\alpha$-parameter. In addition, fragmentation is common in systems that are simulated using hydrodynamic cloud collapse, with more fragments of larger mass formed in more massive systems. In contrast, if disc formation is limited by magnetic fields, fragmentation is suppressed entirely.


title: Modeling in-ice radio propagation with parabolic equation methods

abstract: We investigate the use of parabolic equation (PE) methods for solving radio-wave propagation in polar ice. PE methods provide an approximate solution to Maxwell's equations, in contrast to full-field solutions such as finite-difference-time-domain (FDTD) methods, yet provide a more complete model of propagation than simple geometric ray-tracing (RT) methods that are the current state of the art for simulating in-ice radio detection of neutrino-induced cascades. PE are more computationally efficient than FDTD methods, and more flexible than RT methods, allowing for the inclusion of diffractive effects, and modeling of propagation in regions that cannot be modeled with geometric methods. We present a new PE approximation suited to the in-ice case. We conclude that current ray-tracing methods may be too simplistic in their treatment of ice properties, and their continued use could overestimate experimental sensitivity for in-ice neutrino detection experiments. We discuss the implications for current in-ice Askaryan-type detectors and for the upcoming Radar Echo Telescope; two families of experiments for which these results are most relevant. We suggest that PE methods be investigated further for in-ice radio applications.


title: An Amusing Look at the Host of the Periodic Nuclear Transient ASASSN-14ko Reveals a Second AGN

abstract: We present Multi-Unit Spectroscopic Explorer (MUSE) integral-field spectroscopy of ESO 253$-$G003, which hosts a known Active Galactic Nucleus (AGN) and the periodic nuclear transient ASASSN-14ko, observed as part of the All-weather MUse Supernova Integral-field of Nearby Galaxies (AMUSING) survey. The MUSE observations reveal that the inner region hosts two AGN separated by $1.4\pm0.1~\rm{kpc}$ ($\approx 1.\!\!^{\prime\prime}7$). The brighter nucleus has asymmetric broad, permitted emission-line profiles and is associated with the archival AGN designation. The fainter nucleus does not have a broad emission-line component but exhibits other AGN characteristics, including $v_{\rm{FWHM}}\approx 700~\rm{km}~\rm{s}^{-1}$ forbidden line emission, $\log_{10}(\rm{[OIII]}/\rm{H}\beta) \approx 1.1$, and high excitation potential emission lines such as [Fe$~$VII]$~\lambda6086$ and He$~$II$~\lambda4686$. The host galaxy exhibits a disturbed morphology with large kpc-scale tidal features, potential outflows from both nuclei, and a likely superbubble. A circular relativistic disk model cannot reproduce the asymmetric broad emission-line profiles in the brighter nucleus, but two non-axisymmetric disk models provide good fits to the broad emission-line profiles: an elliptical disk model and a circular disk + spiral arm model. Implications for the periodic nuclear transient ASASSN-14ko are discussed.


title: Primordial black hole dark matter in dilaton-extended two-field Starobinsky inflation

abstract: We investigate the production of primordial black holes and their contribution to the presently observed dark matter in a dilaton two-field extension of Starobinsky's quadratic $f(R)$ model of inflation. The model features a multi-field amplification mechanism which leads to the generation of a sharp peak in the inflationary power spectrum at small wavelengths responsible for the production of primordial black holes. This mechanism is significantly different from single-field models and requires a stochastic treatment during an intermediate phase of the inflationary dynamics. We find that the model leads to a successful phase of effective single-field Starobinsky inflation for wavelengths probed by the cosmic microwave background radiation and explains the observed cold dark matter content in the Universe by the formation of primordial black holes.


title: Floquet Majorana bound states in voltage-biased Planar Josephson Junctions

abstract: We study a planar Josephson junction under an applied DC voltage bias in the presence of an in-plane magnetic field. Upon tuning the bias voltage across the junction, $V_{{\rm J}}$, the two ends of the junction are shown to simultaneously host both $0$- and $\pi$-Majorana modes. These modes can be probed using either a Scanning-Tunneling-Microscopy measurement or through resonant Andreev tunneling from a lead coupled to the junction. While these modes are mostly bound to the junction's ends, they can hybridize with the bulk by absorbing or emitting photons. We analyze this process both numerically and analytically, demonstrating that it can become negligible under typical experimental conditions. Transport signatures of the $0$- and $\pi$-Majorana states are shown to be robust to moderate disorder.


title: Unveiling the rarest morphologies of the LOFAR Two-metre Sky Survey radio source population with self-organised maps

abstract: The Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) is a low-frequency radio continuum survey of the Northern sky at an unparalleled resolution and sensitivity. In order to fully exploit this huge dataset and those produced by the Square Kilometre Array in the next decade, automated methods in machine learning and data-mining will be increasingly essential both for morphological classifications and for identifying optical counterparts to the radio sources. Using self-organising maps (SOMs), a form of unsupervised machine learning, we created a dimensionality reduction of the radio morphologies for the $\sim$25k extended radio continuum sources in the LoTSS first data release, which is only $\sim$2 percent of the final LoTSS survey. We made use of \textsc{PINK}, a code which extends the SOM algorithm with rotation and flipping invariance, increasing its suitability and effectiveness for training on astronomical sources. After training, the SOMs can be used for a wide range of science exploitation and we present an illustration of their potential by finding an arbitrary number of morphologically rare sources in our training data (424 square degrees) and subsequently in an area of the sky ($\sim$5300 square degrees) outside the training data. Objects found in this way span a wide range of morphological and physical categories: extended jets of radio active galactic nuclei, diffuse cluster haloes and relics, and nearby spiral galaxies. Finally, to enable accessible, interactive, and intuitive data exploration, we showcase the LOFAR-PyBDSF Visualisation Tool, which allows users to explore the LoTSS dataset through the trained SOMs.


title: TDCOSMO VI: Distance Measurements in Time-delay Cosmography under the Mass-sheet transformation

abstract: Time-delay cosmography with gravitationally lensed quasars plays an important role in anchoring the absolute distance scale and hence measuring the Hubble constant, $H_{0}$, independent of traditional distance ladder methodology. A current potential limitation of time delay distance measurements is the "mass-sheet transformation" (MST) which leaves the lensed imaging unchanged but changes the distance measurements and the derived value of $H_0$. In this work we show that the standard method of addressing the MST in time delay cosmography, through a combination of high-resolution imaging and the measurement of the stellar velocity dispersion of the lensing galaxy, depends on the assumption that the ratio, $D_{\rm s}/D_{\rm ds}$, of angular diameter distances to the background quasar and between the lensing galaxy and the quasar can be constrained. This is typically achieved through the assumption of a particular cosmological model. Previous work (TDCOSMO IV) addressed the mass-sheet degeneracy and derived $H_{0}$ under the assumption of $\Lambda$CDM model. In this paper we show that the mass sheet degeneracy can be broken without relying on a specific cosmological model by combining lensing with relative distance indicators such as supernovae type Ia and baryon acoustic oscillations, which constrain the shape of the expansion history and hence $D_{\rm s}/D_{\rm ds}$. With this approach, we demonstrate that the mass-sheet degeneracy can be constrained in a cosmological-model-independent way, and hence model-independent distance measurements in time-delay cosmography under mass-sheet transformations can be obtained.


title: The bi-conical vector model at $1/N$

abstract: We study finite $N$ aspects of the $O(m)\times O(N-m)$ vector model with quartic interactions in general $2\leq d \leq 6$ spacetime dimensions. This model has recently been shown to display the phenomenon of persistent symmetry breaking at a perturbative Wilson-Fisher-like fixed point in $d=4-\epsilon$ dimensions. The large rank limit of the bi-conical model displays a conformal manifold and a moduli space of vacua. We find a set of three double trace scalar operators that are respectively irrelevant, relevant and marginal deformations of the conformal manifold in general $d$. We calculate the anomalous dimensions of the single and multi-trace scalar operators to the first sub-leading order in the large rank expansion. The anomalous dimension of the marginal operator does not vanish in general, indicating that the conformal manifold is lifted at finite $N$. In the case of equal ranks we are able to derive explicitly the scaling dimensions of various operators as functions of only $d$.


title: The Characteristic Momentum of Radiatively Cooling Energy-Driven Galactic Winds

abstract: Energy injection by supernovae may drive hot supersonic galactic winds in rapidly star-forming galaxies, driving metal-enriched gas into the circumgalactic medium and potentially accelerating cool gas. If sufficiently mass-loaded, such flows become radiative within the wind-driving region, reducing the overall mass outflow rate from the host galaxy. We show that this sets a maximum on the total outflow momentum for hot energy-driven winds. For a spherical wind of Solar metallicity driven by continuous star formation, $\dot p_\rm{max} \simeq 1.9\times10^4\ M_\odot/\rm{yr\ km/s}\ (\alpha/0.9)^{0.86}(R_\star/300\ \rm{pc})^{0.14}(\dot M_\star/20\ M_\odot/\rm{yr})^{0.86}$, where $\alpha$ is the fraction of supernova energy that thermalizes the wind, and $\dot M_\star$ and $R_\star$ are the star formation rate and radius of the wind-driving region. This maximum momentum for hot winds can also apply to cool, ionized outflows that are typically observed in starburst galaxies, if the hot wind undergoes bulk radiative cooling or if the hot wind transfers mass and momentum to cool clouds within the flow. We show that requiring the hot wind to undergo single-phase cooling on large scales sets a minimum on the total outflow momentum rate. These maximum and minimum outflow momenta have similar values, setting a characteristic momentum rate of hot galactic winds that can become radiative on large scales. We find that most observations of photoionized outflow wind momentum fall below the theoretical maximum and thus may be signatures of cooling hot flows. On the other hand, many systems fall below the minimum momentum required for bulk cooling, indicating that perhaps the cool material observed has instead been entrained in or mixed with the hot flow.


title: Signatures of global symmetry violation in relative entropies and replica wormholes

abstract: It is widely believed that exact global symmetries do not exist in theories that admit quantum black holes. Here we propose a way to quantify the degree of global symmetry violation in the Hawking radiation of a black hole by using certain relative entropies. While the violations of global symmetry that we consider are non-perturbative effects, they nevertheless give $\mathcal{O}(1)$ contributions to the relative entropy after the Page time. Furthermore, using "island" formulas, these relative entropies can be computed within semi-classical gravity, which we demonstrate with explicit examples. These formulas give a rather precise operational sense to the statement that a global charge thrown into an old black hole will be lost after a scrambling time. The relative entropies considered here may also be computed using a replica trick. At integer replica index, the global symmetry violating effects manifest themselves as charge flowing through the replica wormhole.


title: Towards NNGP-guided Neural Architecture Search

abstract: The predictions of wide Bayesian neural networks are described by a Gaussian process, known as the Neural Network Gaussian Process (NNGP). Analytic forms for NNGP kernels are known for many models, but computing the exact kernel for convolutional architectures is prohibitively expensive. One can obtain effective approximations of these kernels through Monte-Carlo estimation using finite networks at initialization. Monte-Carlo NNGP inference is orders-of-magnitude cheaper in FLOPs compared to gradient descent training when the dataset size is small. Since NNGP inference provides a cheap measure of performance of a network architecture, we investigate its potential as a signal for neural architecture search (NAS). We compute the NNGP performance of approximately 423k networks in the NAS-bench 101 dataset on CIFAR-10 and compare its utility against conventional performance measures obtained by shortened gradient-based training. We carry out a similar analysis on 10k randomly sampled networks in the mobile neural architecture search (MNAS) space for ImageNet. We discover comparative advantages of NNGP-based metrics, and discuss potential applications. In particular, we propose that NNGP performance is an inexpensive signal independent of metrics obtained from training that can either be used for reducing big search spaces, or improving training-based performance measures.


title: Understanding College Students' Phone Call Behaviors Towards a Sustainable Mobile Health and Wellbeing Solution

abstract: During the transition from high school to on-campus college life, a student leaves home and starts facing enormous life changes, including meeting new people, more responsibilities, being away from family, and academic challenges. These recent changes lead to an elevation of stress and anxiety, affecting a student's health and wellbeing. With the help of smartphones and their rich collection of sensors, we can continuously monitor various factors that affect students' behavioral patterns, such as communication behaviors associated with their health, wellbeing, and academic success. In this work, we try to assess college students' communication patterns (in terms of phone call duration and frequency) that vary across various geographical contexts (e.g., dormitories, classes, dining) during different times (e.g., epochs of a day, days of a week) using visualization techniques. Findings from this work will help foster the design and delivery of smartphone-based health interventions; thereby, help the students adapt to the changes in life.


title: Subleading BMS Charges and The Lorentz Group

abstract: The extended BMS group includes supertranslation, dual supertranslation and Lorentz transformations. The generators of these symmetries can be classified according to their parity into "electric" and "magnetic" types. Using a multipole expansion of gravitational sources in the Newman-Penrose formalism, we associate each one of these charges with a particular moment. At leading order, the "electric" and "magnetic" monopole moments correspond respectively to supertranslations and dual supertranslations. At the first subleading order, the "electric" and "magnetic" dipole moments correspond respectively to boosts and rotations, therefore comprising the entire Lorentz group. Electric-magnetic type of duality then rotates the boost and rotation generators into each other.


title: Objective trajectories in hybrid classical-quantum dynamics

abstract: Consistent dynamics which couples classical and quantum degrees of freedom exists, provided it is stochastic. This dynamics is linear in the hybrid state, completely positive and trace preserving. One application of this is to study the back-reaction of quantum fields on space-time which does not suffer from the pathologies of the semi-classical equations. Here we introduce several toy models in which to study hybrid classical-quantum evolution, including a qubit coupled to a particle in a potential, and a quantum harmonic oscillator coupled to a classical one. We present an unravelling approach to calculate the dynamics, and provide code to numerically simulate it. Unlike the purely quantum case, the trajectories (or histories) of this unravelling can be unique, conditioned on the classical degrees of freedom for discrete realisations of the dynamics, when different jumps in the classical degrees of freedom are accompanied by the action of unique operators on the quantum system. As a result, the ``measurement postulate'' of quantum theory is not needed; quantum systems become classical because they interact with a fundamentally classical field.


title: Constraints on ultralight scalar bosons within black hole spin measurements from LIGO-Virgo's GWTC-2

abstract: Clouds of ultralight bosons - such as axions - can form around a rapidly spinning black hole, if the black hole radius is comparable to the bosons' wavelength. The cloud rapidly extracts angular momentum from the black hole, and reduces it to a characteristic value that depends on the boson's mass as well as on the black hole mass and spin. Therefore, a measurement of a black hole mass and spin can be used to reveal or exclude the existence of such bosons. Using the black holes released by LIGO and Virgo in their GWTC-2, we perform a simultaneous measurement of the black hole spin distribution at formation and the mass of the scalar boson. We find that the data strongly disfavors the existence of scalar bosons in the mass range between $1.3\times10^{-13}~\mathrm{eV}$ and $2.7\times10^{-13}~\mathrm{eV}$ for a decay constant $f_a\gtrsim 10^{14}~\mathrm{GeV}$. The statistical evidence is mostly driven by the two {binary black holes} systems GW190412 and GW190517, which host rapidly spinning black holes. The region where bosons are excluded narrows down if these two systems merged shortly ($\sim 10^5$ years) after the black holes formed.


title: Isospectral twirling and quantum chaos

abstract: We show that the most important measures of quantum chaos like frame potentials, scrambling, Loschmidt echo, and out-of-time-order correlators (OTOCs) can be described by the unified framework of the isospectral twirling, namely the Haar average of a $k$-fold unitary channel. We show that such measures can then be always cast in the form of an expectation value of the isospectral twirling. In literature, quantum chaos is investigated sometimes through the spectrum and some other times through the eigenvectors of the Hamiltonian generating the dynamics. We show that, by exploiting random matrix theory, these measures of quantum chaos clearly distinguish the finite time profiles of probes to quantum chaos corresponding to chaotic spectra given by the Gaussian Unitary Ensemble (GUE) from the integrable spectra given by Poisson distribution and the Gaussian Diagonal Ensemble (GDE). On the other hand, we show that the asymptotic values do depend on the eigenvectors of the Hamiltonian. We see that the isospectral twirling of Hamiltonians with eigenvectors stabilizer states does not possess chaotic features, unlike those Hamiltonians whose eigenvectors are taken from the Haar measure. As an example, OTOCs obtained with Clifford resources decay to higher values compared with universal resources. Finally, we show a crossover in the OTOC behavior between a class of integrable models and quantum chaos.


title: SDSS-IV MaNGA: Refining Strong Line Diagnostic Classifications Using Spatially Resolved Gas Dynamics

abstract: We use the statistical power of the MaNGA integral-field spectroscopic galaxy survey to improve the definition of strong line diagnostic boundaries used to classify gas ionization properties in galaxies. We detect line emission from 3.6 million spaxels distributed across 7400 individual galaxies spanning a wide range of stellar masses, star formation rates, and morphological types, and find that the gas-phase velocity dispersion sigma_HAlpha correlates strongly with traditional optical emission line ratios such as [S II]/HAlpha, [N II]/HAlpha, [O I]/HAlpha, and [O III]/HBeta. Spaxels whose line ratios are most consistent with ionization by galactic HII regions exhibit a narrow range of dynamically cold line of sight velocity distributions (LOSVDs) peaked around 25 km/s corresponding to a galactic thin disk, while those consistent with ionization by active galactic nuclei (AGN) and low-ionization emission-line regions (LI(N)ERs) have significantly broader LOSVDs extending to 200 km/s. Star-forming, AGN, and LI(N)ER regions are additionally well separated from each other in terms of their stellar velocity dispersion, stellar population age, HAlpha equivalent width, and typical radius within a given galaxy. We use our observations to revise the traditional emission line diagnostic classifications so that they reliably identify distinct dynamical samples both in two-dimensional representations of the diagnostic line ratio space and in a multi-dimensional space that accounts for the complex folding of the star forming model surface. By comparing the MaNGA observations to the SDSS single-fiber galaxy sample we note that the latter is systematically biased against young, low metallicity star-forming regions that lie outside of the 3 arcsec fiber footprint.


title: OGLE-2007-BLG-224L: Confirmation of Terrestrial Parallax

abstract: We present limits on the lens flux of OGLE-2007-BLG-224 based on MagAO imaging taken seven years after the microlensing event. At the time of the observations, the lens should have been separated from the microlensing source by 292 mas. However, no new sources are detected with MagAO. We place an upper limit on the lens flux of $H>20.57$. This measurement supports the conclusion of Gould et al. (2009) that the lens in this event should be a brown dwarf. This is the first test of a prediction based on the terrestrial microlens parallax effect and the first AO confirmation of a sub-stellar/dark microlens.


title: Football tracking networks: Beyond event-based connectivity

abstract: We propose using Network Science as a complementary tool to analyze player and team behavior during a football match. Specifically, we introduce four kinds of networks based on different ways of interaction between players. Our approach's main novelty is to use tracking datasets to create football tracking networks, instead of constructing and analyzing the traditional networks based on events. In this way, we are able to capture player interactions that go beyond passes and introduce the concepts of (a) Ball Flow Networks, (b) Marking Networks, (c) Signed Proximity Networks and (d) Functional Coordination Networks. After defining the methodology for creating each kind of network, we show some examples using tracking datasets from four different matches of LaLiga Santander. Finally, we discuss some of the applications, limitations, and further improvements of football tracking networks.


title: GANMEX: One-vs-One Attributions Guided by GAN-based Counterfactual Explanation Baselines

abstract: Attribution methods have been shown as promising approaches for identifying key features that led to learned model predictions. While most existing attribution methods rely on a baseline input for performing feature perturbations, limited research has been conducted to address the baseline selection issues. Poor choices of baselines limit the ability of one-vs-one (1-vs-1) explanations for multi-class classifiers, which means the attribution methods were not able to explain why an input belongs to its original class but not the other specified target class. 1-vs-1 explanation is crucial when certain classes are more similar than others, e.g. two bird types among multiple animals, by focusing on key differentiating features rather than shared features across classes. In this paper, we present GAN-based Model EXplainability (GANMEX), a novel approach applying Generative Adversarial Networks (GAN) by incorporating the to-be-explained classifier as part of the adversarial networks. Our approach effectively selects the counterfactual baseline as the closest realistic sample belong to the target class, which allows attribution methods to provide true 1-vs-1 explanations. We showed that GANMEX baselines improved the saliency maps and led to stronger performance on perturbation-based evaluation metrics over the existing baselines. Existing attribution results are known for being insensitive to model randomization, and we demonstrated that GANMEX baselines led to better outcome under the cascading randomization of the model.


title: Collision in the interior of wormhole

abstract: The Schwarzschild wormhole has been interpreted as an entangled state. If Alice and Bob fall into each of the black hole, they can meet in the interior. We interpret this meeting in terms of the quantum circuit that prepares the entangled state. Alice and Bob create growing perturbations in the circuit, and we argue that the overlap of these perturbations represents their meeting. We compare the gravity picture with circuit analysis, and identify the post-collision region as the region storing the gates that are not affected by any of the perturbations.


title: Estimating heating times in periodically driven quantum many-body systems via avoided crossing spectroscopy

abstract: Periodic driving of a quantum (or classical) many-body system can alter the systems properties significantly and therefore has emerged as a promising way to engineer exotic quantum phases, such as topological insulators and discrete time crystals. A major limitation in such setups, is that generally interacting, driven systems will heat up over time and lose the desired properties. Understanding the relevant time scales is thus an important topic in the field and so far, there have only been few approaches to determine heating times for a concrete system quantitatively, and in a computationally efficient way. In this article we propose a new approach, based on building the heating rate from microscopic processes, encoded in avoided level crossings of the Floquet propagator. We develop a method able to resolve individual crossings and show how to construct the heating rate based on these. The method is closely related to the Fermi Golden Rule approach for weak drives, but can go beyond it, since it captures non-perturbative effects by construction. This enables our method to be applicable in scenarios such as the heating time of discrete time crystals or frequency dependent couplings, which are very relevant for Floquet engineering, where previously no efficient methods for estimating heating times were available.


title: Extremal Eigenvalues Of The Conformal Laplacian Under Sire-Xu Normalization

abstract: Let $(M^n,g)$ be a closed Riemannian manifold of dimension $n\ge 3$. We study the variational properties of the $k$-th eigenvalue functional $\tilde g\in[g] \mapsto \lambda_k(L_{\tilde g})$ under a non-volume normalization proposed by Sire-Xu. We discuss necessary conditions for the existence of extremal eigenvalues under such normalization. Also, we discuss the general existence problem when $k=1$.


title: Policing Chronic and Temporary Hot Spots of Violent Crime: A Controlled Field Experiment

abstract: Hot-spot-based policing programs aim to deter crime through increased proactive patrols at high-crime locations. While most hot spot programs target easily identified chronic hot spots, we introduce models for predicting temporary hot spots to address effectiveness and equity objectives for crime prevention, and present findings from a crossover experiment evaluating application of hot spot predictions to prevent serious violent crime in Pittsburgh, PA. Over a 12-month experimental period, the Pittsburgh Bureau of Police assigned uniformed patrol officers to weekly predicted chronic and temporary hot spots of serious violent crimes comprising 0.5 percent of the city's area. We find statistically and practically significant reductions in serious violent crime counts within treatment hot spots as compared to control hot spots, with an overall reduction of 25.3 percent in the FBI-classified Part 1 Violent (P1V) crimes of homicide, rape, robbery, and aggravated assault, and a 39.7 percent reduction of African-American and other non-white victims of P1V crimes. We find that temporary hot spots increase spatial dispersion of patrols and have a greater percentage reduction in P1V crimes than chronic hot spots but fewer total number of crimes prevented. Only foot patrols, not car patrols, had statistically significant crime reductions in hot spots. We find no evidence of crime displacement; instead, we find weakly statistically significant spillover of crime prevention benefits to adjacent areas. In addition, we find no evidence that the community-oriented hot spot patrols produced over-policing arrests of minority or other populations.


title: "Ears" formation in Supernova Remnants: Overhearing an interaction history with bipolar circumstellar structures

abstract: A characteristic feature that is frequently met in nearby supernova remnants (SNRs) is the existence of two antisymmetric, local protrusions that are projected as two "ears" in the morphology of the nebula. In this work, we present a novel scenario for the "ear" formation process according to which the two lobes are formed through the interaction of the SNR with a bipolar circumstellar medium (CSM) that was surrounding the explosion center. We conduct two dimensional hydrodynamic simulations and we show that the SNR shock breakout from the bipolar CSM triggers the inflation of two opposite protrusions at the equator of the remnant that retain their size and shape for several hundreds up to a few thousand years of the SNR evolution. We run a set of models by varying the supernova (SN) and CSM properties and we demonstrate that the extracted results reveal a good agreement with the observables, regarding the "ears" sizes, lifespan, morphology and kinematics. We discuss the plausibility of our model in nature and we suggest that the most likely progenitors of the "ear-carrying" SNRs are the Luminous Blue Variables or the Red/Yellow Supergiants for the SNRs resulted by core collapse SN events and the symbiotic binaries or the planetary nebulae for the remnants formed by Type Ia SNe. Finally, we compare our model with other "ear" formation models of the literature and we show that there are distinctive differences among them, concerning the "ears" orientation and the phase in which the "ear" formation process occurs.


title: A WLAV-based Robust Hybrid State Estimation using Circuit-theoretic Approach

abstract: For reliable and secure power grid operation, AC state-estimation (ACSE) must provide certain guarantees of convergence while being resilient against bad-data. This paper develops a circuit-theoretic weighted least absolute value (WLAV) based hybrid ACSE that satisfies these needs to overcome some of the limitations of existing ACSE methods. Hybrid refers to the inclusion of RTU and PMU measurement data, and the use of the LAV objective function enables automatic rejection of bad data while providing clear identification of suspicious measurements from the sparse residual vector. Taking advantage of linear construction of the measurement models in circuit-theoretic approach, the proposed hybrid SE is formulated as a LP problem with guaranteed convergence. To address efficiency, we further develop problem-specific heuristics for fast convergence. To validate the efficacy of the proposed approach, we run ACSE on large cases and compare the results against WLS-based algorithms. We further demonstrate the advantages of our solution methodology over standard commercial LP solvers through comparison of runtime and convergence performance.


title: Machine Learning Based Path Planning for Improved Rover Navigation (Pre-Print Version)

abstract: Enhanced AutoNav (ENav), the baseline surface navigation software for NASA's Perseverance rover, sorts a list of candidate paths for the rover to traverse, then uses the Approximate Clearance Evaluation (ACE) algorithm to evaluate whether the most highly ranked paths are safe. ACE is crucial for maintaining the safety of the rover, but is computationally expensive. If the most promising candidates in the list of paths are all found to be infeasible, ENav must continue to search the list and run time-consuming ACE evaluations until a feasible path is found. In this paper, we present two heuristics that, given a terrain heightmap around the rover, produce cost estimates that more effectively rank the candidate paths before ACE evaluation. The first heuristic uses Sobel operators and convolution to incorporate the cost of traversing high-gradient terrain. The second heuristic uses a machine learning (ML) model to predict areas that will be deemed untraversable by ACE. We used physics simulations to collect training data for the ML model and to run Monte Carlo trials to quantify navigation performance across a variety of terrains with various slopes and rock distributions. Compared to ENav's baseline performance, integrating the heuristics can lead to a significant reduction in ACE evaluations and average computation time per planning cycle, increase path efficiency, and maintain or improve the rate of successful traverses. This strategy of targeting specific bottlenecks with ML while maintaining the original ACE safety checks provides an example of how ML can be infused into planetary science missions and other safety-critical software.


title: Discovering alignment relations with Graph Convolutional Networks: a biomedical case study

abstract: Knowledge graphs are freely aggregated, published, and edited in the Web of data, and thus may overlap. Hence, a key task resides in aligning (or matching) their content. This task encompasses the identification, within an aggregated knowledge graph, of nodes that are equivalent, more specific, or weakly related. In this article, we propose to match nodes within a knowledge graph by (i) learning node embeddings with Graph Convolutional Networks such that similar nodes have low distances in the embedding space, and (ii) clustering nodes based on their embeddings, in order to suggest alignment relations between nodes of a same cluster. We conducted experiments with this approach on the real world application of aligning knowledge in the field of pharmacogenomics, which motivated our study. We particularly investigated the interplay between domain knowledge and GCN models with the two following focuses. First, we applied inference rules associated with domain knowledge, independently or combined, before learning node embeddings, and we measured the improvements in matching results. Second, while our GCN model is agnostic to the exact alignment relations (e.g., equivalence, weak similarity), we observed that distances in the embedding space are coherent with the ``strength'' of these different relations (e.g., smaller distances for equivalences), letting us considering clustering and distances in the embedding space as a means to suggest alignment relations in our case study.


title: BeyondPlanck II. CMB map-making through Gibbs sampling

abstract: We present a Gibbs sampling solution to the map-making problem for CMB measurements, building on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps; noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods, and allows us for the first time to bring the destriping baseline length to a single sample. We apply the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we may compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction, and this methodology forms the conceptual basis for the map-making algorithm employed in the BeyondPlanck framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.


title: Effective field theory analysis of composite higgsino-like and wino-like thermal relic dark matter

abstract: We study the effective field theory (including operators up to dimension five) of models in which dark matter is composite, consisting of either an electroweak doublet Dirac fermion (`higgsino-like dark matter') or an electroweak triplet Majorana fermion (`wino-like dark matter'). Some of the dimension-five operators in the former case cause mass splittings between the neutralino and chargino states, leading to a depleted rate of coannihilations and viable thermal relic dark matter with masses of the order of tens to hundreds of GeV rather than the usual pure higgsino thermal relic mass of 1 TeV. No such effects are found in the latter case (where the usual thermal relic mass is 3 TeV). Other operators, present for both wino- and higgsino-like dark matter, correspond to inelastic electromagnetic dipole moment interactions and annihilation through these can lead to viable models with dark matter masses up by an order of magnitude compared to the usual values.


title: Solving Puzzles in Deformed JT Gravity: Phase Transitions and Non-Perturbative Effects

abstract: Recent work has shown that certain deformations of the scalar potential in Jackiw-Teitelboim gravity can be written as double-scaled matrix models. However, some of the deformations exhibit an apparent breakdown of unitarity in the form of a negative spectral density at disc order. We show here that the source of the problem is the presence of a multi-valued solution of the leading order matrix model string equation. While for a class of deformations we fix the problem by identifying a first order phase transition, for others we show that the theory is both perturbatively and non-perturbatively inconsistent. Aspects of the phase structure of the deformations are mapped out, using methods known to supply a non-perturbative definition of undeformed JT gravity. Some features are in qualitative agreement with a semi-classical analysis of the phase structure of two-dimensional black holes in these deformed theories.


title: Role of impurity clusters for the current-driven motion of magnetic Skyrmions

abstract: We study how impurities influence the current-induced dynamics of magnetic Skyrmions moving in a racetrack geometry. For this, we solve numerically the generalized Landau-Lifshitz-Gilbert equation extended by the current-induced spin transfer torque. In particular, we investigate two classes of impurities, non-conducting and magnetic impurities. The former are magnetically rigid objects and yield to an inhomogeneous current density over the racetrack which we determine separately by solving the fundamental electrostatic equations. In contrast, magnetic impurities leave the applied current density homogeneous throughout the stripe. Depending on parameters, we observe four different scenarios of Skyrmion motions in the presence of disorder, the Skyrmion decay, the pinning, the creation of additional Skyrmions, and ordinary Skyrmion passage. We calculate and discuss phase diagrams in dependence of the impurity concentration and radii of the impurities.


title: Search for the rare decay of the W boson into a pion and a photon in proton-proton collisions at $\sqrt{s} = $ 13 TeV

abstract: A search is performed for the rare decay W$^\pm\to\pi^\pm\gamma$ in proton-proton collisions at $\sqrt{s} =$ 13 TeV. Data corresponding to an on W integrated luminosity of 137 fb$^{-1}$ were collected during 2016 to 2018 with the CMS detector. This analysis exploits a novel search strategy based on W boson production in top quark pair events. An inclusive search for the W$^\pm\to\pi^\pm\gamma$ decay is not optimal at the LHC because of the high trigger thresholds. Instead, a trigger selection is exploited in which the W boson originating from one of the top quarks is used to tag the event in a leptonic decay. The W boson emerging from the other top quark is used to search for the W$^\pm\to\pi^\pm\gamma$ signature. Such decays are characterized by an isolated track pointing to a large energy deposit, and by an isolated photon of large transverse momentum. The presence of b quark jets reduces the background from the hadronization of light-flavor quarks and gluons. The W$^\pm\to\pi^\pm\gamma$ decay is not observed. An upper exclusion limit is set to this branching fraction, corresponding to 1.50 $\times$ 10$^{-5}$ at 95% confidence level, whereas the expected upper limit exclusion limit is 0.85 $^{+0.52}_{-0.29}$ $\times$ 10$^{-5}$.


title: Gelfand-Tsetlin modules: canonicity and calculations

abstract: In this paper, we give a more down-to-earth introduction to the connection between Gelfand-Tsetlin modules over $\mathfrak{gl}_n$ and diagrammatic KLRW algebras, and develop some of its consequences. In addition to a new proof of this description of the category Gelfand-Tsetlin modules appearing in earlier work, we show three new results of independent interest: (1) we show that every simple Gelfand-Tsetlin module is a canonical module in the sense of Early, Mazorchuk and Vishnyakova, and characterize when two maximal ideals have isomorphic canonical modules, (2) we show that the dimensions of Gelfand-Tsetlin weight spaces in simple modules can be computed using an appropriate modification of Leclerc's algorithm for computing dual canonical bases, and (3) we construct a basis of the Verma modules of $\mathfrak{sl}_n$ which consists of generalized eigenvectors for the Gelfand-Tsetlin subalgebra. Furthermore, we present computations of multiplicities and Gelfand-Kirillov dimensions for all integral Gelfand-Tsetlin modules in ranks 3 and 4; unfortunately, for ranks $>4$, our computers are not adequate to perform these computations.


title: Modeling Interplanetary Expansion and Deformation of CMEs with ANTEATR-PARADE I: Relative Contribution of Different Forces

abstract: Coronal Mass Ejections (CMEs) are key drivers of space weather activity but most predictions have been limited to the expected arrival time of a CME, rather than the internal properties that affect the severity of an impact. Many properties, such as the magnetic field density and mass density, follow conservation laws and vary systematically with changes in the size of a CME. We present ANTEATR-PARADE, the newest version of the ANTEATR arrival time model, which now includes physics-driven changes in the size and shape of both the CME's central axis and its cross section. Internal magnetic and thermal and external drag forces affect the acceleration of the CME in different directions, inducing asymmetries between the radial and perpendicular directions. These improvements should lead to more realistic CME velocities, both bulk and expansion, sizes and shapes, and internal properties. We present the model details, an initial illustration of the general behavior, and a study of the relative importance of the different forces. The model shows a pancaking of both the cross section and central axis of the CME so that their radial extent becomes smaller than their extent in the perpendicular direction. We find that the initial velocities, drag, any form of cross section expansion, and the precise form of thermal expansion have strong effects. The results are less sensitive to axial forces and the specific form of the cross section expansion.


title: swdpwr: A SAS Macro and An R Package for Power Calculation in Stepped Wedge Cluster Randomized Trials

abstract: Background and objective: The stepped wedge cluster randomized trial is a study design increasingly used for public health intervention evaluations. Most previous literature focuses on power calculations for this particular type of cluster randomized trials for continuous outcomes, along with an approximation to this approach for binary outcomes. Although not accurate for binary outcomes, it has been widely used. To improve the approximation for binary outcomes, two new methods for stepped wedge designs (SWDs) of binary outcomes have recently been published. However, these new methods have not been implemented in publicly available software. The objective of this paper is to present power calculation software for SWDs in various settings for both continuous and binary outcomes. Methods: We have developed a SAS macro %swdpwr and an R package swdpwr for power calculation in SWDs. Different scenarios including cross-sectional and cohort designs, binary and continuous outcomes, marginal and conditional models, three link functions, with and without time effects are accommodated in this software. Results: swdpwr provides an efficient tool to support investigators in the design and analysis of stepped wedge cluster randomized trails. swdpwr addresses the implementation gap between newly proposed methodology and their application to obtain more accurate power calculations in SWDs. Conclusions: This user-friendly software makes the new methods more accessible and incorporates as many variations as currently available, which were not supported in other related packages. swdpwr is implemented under two platforms: SAS and R, satisfying the needs of investigators from various backgrounds.


title: On the Formalism of The Screening Paradox

abstract: Bayes' Theorem imposes inevitable limitations on the accuracy of screening tests by tying the test's predictive value to the disease prevalence. The aforementioned limitation is independent of the adequacy and make-up of the test and thus implies inherent Bayesian limitations to the screening process itself. As per the WHO's $Wilson-Jungner$ criteria, one of the prerequisite steps before undertaking screening is to ensure that a treatment for the condition screened exists. However, in so doing, a paradox, henceforth termed the screening paradox, ensues. If a disease process is screened for and subsequently treated, its prevalence would drop in the population, which as per Bayes' theorem, would make the tests' predictive value drop in return. Put another way, a very powerful screening test would, by performing and succeeding at the very task it was developed to do, paradoxically reduce its ability to correctly identify individuals with the disease it screens for in the future. Where $J$ is Youden's statistic (sensitivity [$a$] + specificity [$b$] - 1), and $\phi$ is the prevalence, the ratio of positive predictive values at subsequent time $k$, $\rho(\phi_{k})$, over the original $\rho(\phi_{0})$ at $t_0$ is given by: $\zeta(\phi_{0},k) = \frac{\rho(\phi_{k})}{\rho(\phi_{0})} =\frac{\phi_k(1-b)+J\phi_0\phi_k}{\phi_0(1-b)+J\phi_0\phi_k}$ In this manuscript, we explore the mathematical model which formalizes said screening paradox and explore its implications for population level screening programs. In particular, we define the number of positive test iterations (PTI) needed to reverse the effects of the paradox as follows: $n_{i\phi_e}=\left\lceil\frac{ln\left[\frac{\omega\phi_e\phi_k-\omega\phi_e}{\omega\phi_e\phi_k-\phi_k}\right]}{2ln\omega}\right\rceil$ where $\omega$ is the square root of the positive likelihood ratio (LR+).


title: FastPathology: An open-source platform for deep learning-based research and decision support in digital pathology

abstract: Deep convolutional neural networks (CNNs) are the current state-of-the-art for digital analysis of histopathological images. The large size of whole-slide microscopy images (WSIs) requires advanced memory handling to read, display and process these images. There are several open-source platforms for working with WSIs, but few support deployment of CNN models. These applications use third-party solutions for inference, making them less user-friendly and unsuitable for high-performance image analysis. To make deployment of CNNs user-friendly and feasible on low-end machines, we have developed a new platform, FastPathology, using the FAST framework and C++. It minimizes memory usage for reading and processing WSIs, deployment of CNN models, and real-time interactive visualization of results. Runtime experiments were conducted on four different use cases, using different architectures, inference engines, hardware configurations and operating systems. Memory usage for reading, visualizing, zooming and panning a WSI were measured, using FastPathology and three existing platforms. FastPathology performed similarly in terms of memory to the other C++ based application, while using considerably less than the two Java-based platforms. The choice of neural network model, inference engine, hardware and processors influenced runtime considerably. Thus, FastPathology includes all steps needed for efficient visualization and processing of WSIs in a single application, including inference of CNNs with real-time display of the results. Source code, binary releases and test data can be found online on GitHub at https://github.com/SINTEFMedtek/FAST-Pathology/.


title: Quasiperiodicity, band topology, and moir\'e graphene

abstract: A number of moir\'e graphene systems have nearly flat topological bands where electron motion is strongly correlated. Though microscopically these systems are only quasiperiodic, they can typically be treated as translation invariant to an excellent approximation. Here we reconsider this question for magic angle twisted bilayer graphene that is nearly aligned with a hexagonal boron nitride(h-BN) substrate. We carefully study the effect of the periodic potential induced by h-BN on the low energy physics. The combination of this potential and the moir\'e lattice produced by the twisted graphene generates a quasi-periodic term that depends on the alignment angle between h-BN and the moir\'e graphene. We find that the alignment angle has a significant impact on both the band gap near charge neutrality and the behavior of electrical transport. We also introduce and study toy models to illustrate how a quasi-periodic potential can give rise to localization and change in transport properties of topological bands.


title: Clifford operators in SU(N)1; N not odd prime

abstract: Farinholt gives a characterization of Clifford operators for qudits; d both odd and even. In this comment it is shown that the necessary gates for the construction of Clifford operators; N both odd and even, are obtained directly from operations that appear in SU(N)1. A witness for W3 states in SU(2)1 is discussed. See e.g. [1-4].


title: Elastic turbulence generates anomalous flow resistance in porous media

abstract: Diverse processes rely on the viscous flow of polymer solutions through porous media. In many cases, the macroscopic flow resistance abruptly increases above a threshold flow rate in a porous medium---but not in bulk solution. The reason why has been a puzzle for over half a century. Here, by directly visualizing the flow in a transparent 3D porous medium, we demonstrate that this anomalous increase is due to the onset of an elastic instability. We establish that the energy dissipated by the unstable flow fluctuations, which vary across pores, generates the anomalous increase in flow resistance through the entire medium. Thus, by linking the pore-scale onset of unstable flow to macroscopic transport, our work provides generally-applicable guidelines for predicting and controlling polymer solution flows.

categories: physics.flu-dyn cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.soft nlin.CD


title: Unsupervised Video Representation Learning by Bidirectional Feature Prediction

abstract: This paper introduces a novel method for self-supervised video representation learning via feature prediction. In contrast to the previous methods that focus on future feature prediction, we argue that a supervisory signal arising from unobserved past frames is complementary to one that originates from the future frames. The rationale behind our method is to encourage the network to explore the temporal structure of videos by distinguishing between future and past given present observations. We train our model in a contrastive learning framework, where joint encoding of future and past provides us with a comprehensive set of temporal hard negatives via swapping. We empirically show that utilizing both signals enriches the learned representations for the downstream task of action recognition. It outperforms independent prediction of future and past.


title: Dilaton-gravity, deformations of the minimal string, and matrix models

abstract: A large class of two-dimensional dilaton-gravity theories in asymptotically AdS$_2$ spacetimes are holographically dual to a matrix integral, interpreted as an ensemble average over Hamiltonians. Viewing these theories as Jackiw-Teitelboim gravity with a gas of defects, we extend this duality to a broader class of dilaton potentials compared to previous work by including conical defects with small deficit angles. In order to do this we show that these theories are equal to the large $p$ limit of a natural deformation of the $(2,p)$ minimal string theory.


title: Recovery of nonlinear terms for reaction diffusion equations from boundary measurements

abstract: We consider the inverse problem of determining a general semilinear term appearing in nonlinear parabolic equations. For this purpose, we derive a new criterion that allows to prove global recovery of some general class of semilinear terms from lateral boundary measurements of solutions of the equation with initial condition fixed at zero. More precisely, we prove, for what seems to be the first time, the unique and stable recovery of general semilinear terms depending on time and space variables independently of the solution of the nonlinear equation from the knowledge of the parabolic Dirichlet-to-Neumann map associated with the solution of the equation with initial condition fixed at zero. Our approach is based on the second linearization of the inverse problem under consideration.


title: Multicomponent CASSCF Revisited: Large Active Spaces are Needed for Qualitatively Accurate Protonic Densities

abstract: Multicomponent methods seek to treat select nuclei, typically protons, fully quantum mechanically and equivalent to the electrons of a chemical system. In such methods, it is well known that due to the neglect of electron-proton correlation, a Hartree-Fock (HF) description of the electron-proton interaction catastrophically fails leading to qualitatively incorrect protonic properties. In single-component quantum chemistry, the qualitative failure of HF is normally indicative of the need for multireference methods such as complete active space self-consistent field (CASSCF). While a multicomponent CASSCF method was implemented nearly twenty years ago, it is only able to perform calculations with very small active spaces (~105 multicomponent configurations). Therefore, in order to extend the realm of applicability of the multicomponent CASSCF method, this study derives and implements a new two-step multicomponent CASSCF method that uses multicomponent heat-bath configuration interaction for the configuration interaction step, enabling calculations with very large active spaces (up to 16 electrons in 48 orbitals). We find that large electronic active spaces are needed to obtain qualitatively accurate protonic densities for the HCN and FHF- molecules. Additionally, the multicomponent CASSCF method implemented here should have further applications for double-well protonic potentials and systems that are inherently electronically multireference.