title: Warping of Radar Data into Camera Image for Cross-Modal Supervision in Automotive Applications

abstract: We present an approach to automatically generate semantic labels for real recordings of automotive range-Doppler (RD) radar spectra. Such labels are required when training a neural network for object recognition from radar data. The automatic labeling approach rests on the simultaneous recording of camera and lidar data in addition to the radar spectrum. By warping radar spectra into the camera image, state-of-the-art object recognition algorithms can be applied to label relevant objects, such as cars, in the camera image. The warping operation is designed to be fully differentiable, which allows backpropagating the gradient computed on the camera image through the warping operation to the neural network operating on the radar data. As the warping operation relies on accurate scene flow estimation, we further propose a novel scene flow estimation algorithm which exploits information from camera, lidar and radar sensors. The proposed scene flow estimation approach is compared against a state-of-the-art scene flow algorithm, and it outperforms it by approximately 30% w.r.t. mean average error. The feasibility of the overall framework for automatic label generation for RD spectra is verified by evaluating the performance of neural networks trained with the proposed framework for Direction-of-Arrival estimation.


title: Optimal dimension dependence of the Metropolis-Adjusted Langevin Algorithm

abstract: Conventional wisdom in the sampling literature, backed by a popular diffusion scaling limit, suggests that the mixing time of the Metropolis-Adjusted Langevin Algorithm (MALA) scales as $O(d^{1/3})$, where $d$ is the dimension. However, the diffusion scaling limit requires stringent assumptions on the target distribution and is asymptotic in nature. In contrast, the best known non-asymptotic mixing time bound for MALA on the class of log-smooth and strongly log-concave distributions is $O(d)$. In this work, we establish that the mixing time of MALA on this class of target distributions is $\widetilde\Theta(d^{1/2})$ under a warm start. Our upper bound proof introduces a new technique based on a projection characterization of the Metropolis adjustment which reduces the study of MALA to the well-studied discretization analysis of the Langevin SDE and bypasses direct computation of the acceptance probability.


title: Oriented expressions of graph properties

abstract: Several graph properties are characterized as the class of graphs that admit an orientation avoiding finitely many oriented structures. For instance, if $F_k$ is the set of homomorphic images of the directed path on $k+1$ vertices, then a graph is $k$-colourable if and only if it admits an orientation with no induced oriented graph in $F_k$. There is a fundamental question underlying this kind of characterizations: given a graph property, $\mathcal{P}$, is there a finite set of oriented graphs, $F$, such that a graph belongs to $\mathcal{P}$ if and only if it admits an orientation with no induced oriented graph in $F$? We address this question by exhibiting necessary conditions upon certain graph classes to admit such a characterization. Consequently, we exhibit an uncountable family of hereditary classes, for which no such finite set exists. In particular, the class of graphs with no holes of prime length belongs to this family.


title: Selecting velocity models using Bayesian Information Criterion

abstract: We present a strategy for selecting the values of elasticity parameters by comparing walk-away vertical seismic profiling data with a multilayered model in the context of Bayesian Information Criterion. We consider $P$-wave traveltimes and assume elliptical velocity dependence. The Bayesian Information Criterion approach requires two steps of optimization. In the first step, we find the signal trajectory and, in the second step, we find media parameters by minimizing the misfit between the model and data.


title: Design of a non-linear human breast diffusion coil

abstract: Diffusion-weighted imaging (DWI) in the female breast is a magnetic resonance imaging (MRI) technique, which may potentially replace or effectively complement both x-ray mammography and contrast-enhanced MR mammography. To further improve specificity of DWI in the breast, stronger and faster diffusion weighting is advantageous. With the target application in mind, the gradient coil is designed on an irregularly shaped semi-open current-carrying surface. Due to the coil former closely fitting the non-spherical target region, non-linear encoding fields become particularly advantageous for achieving locally exceptionally high gradient strengths. As the breast tissue has a predominantly isotropic cellular microstructure, the direction of the diffusion-weighting gradient may be allowed to vary within the target volume. However, due to the quadratic dependency of the b-factor on the gradient strength, variation of the gradient magnitude needs to be carefully controlled. To achieve the above design goals the corresponding multi-objective optimization problem is reformulated as a constrained optimization, allowing for flexible and precise control of the coil properties. A novel constraint is proposed limiting the gradient magnitude variation within every slice while allowing for variations both the direction of the gradient within the slice and the magnitude across the slices. The above innovations enable the design of a unilateral coil for diffusion weighting in the female breast with the local gradient strengths exceeding 1~T/m with highly homogeneous diffusion weighting for imaging in the coronal slice orientation.


title: Data segmentation algorithms: Univariate mean change and beyond

abstract: Data segmentation a.k.a. multiple change point analysis has received considerable attention due to its importance in time series analysis and signal processing, with applications in a variety of fields including natural and social sciences, medicine, engineering and finance. In the first part of this survey, we review the existing literature on the canonical data segmentation problem which aims at detecting and localising multiple change points in the mean of univariate time series. We provide an overview of popular methodologies on their computational complexity and theoretical properties. In particular, our theoretical discussion focuses on the separation rate relating to which change points are detectable by a given procedure, and the localisation rate quantifying the precision of corresponding change point estimators, and we distinguish between whether a homogeneous or multiscale viewpoint has been adopted in their derivation. We further highlight that the latter viewpoint provides the most general setting for investigating the optimality of data segmentation algorithms. Arguably, the canonical segmentation problem has been the most popular framework to propose new data segmentation algorithms and study their efficiency in the last decades. In the second part of this survey, we motivate the importance of attaining an in-depth understanding of strengths and weaknesses of methodologies for the change point problem in a simpler, univariate setting, as a stepping stone for the development of methodologies for more complex problems. We illustrate this with a range of examples showcasing the connections between complex distributional changes and those in the mean. We also discuss extensions towards high-dimensional change point problems where we demonstrate that the challenges arising from high dimensionality are orthogonal to those in dealing with multiple change points.


title: A note on Griffiths' conjecture about the positivity of Chern-Weil forms

abstract: Let $ (E,h) $ be a Griffiths semipositive Hermitian holomorphic vector bundle of rank $ 3 $ over a complex manifold. In this paper, we prove the positivity of the characteristic differential form $ c_1(E,h) \wedge c_2(E,h) - c_3(E,h) $, thus providing a new evidence towards a conjecture by Griffiths about the positivity of the Schur polynomials in the Chern forms of Griffiths semipositive vector bundles. As a consequence, we establish a new chain of inequalities between Chern forms. Moreover, we point out how to obtain the positivity of the second Chern form $ c_2(E,h) $ in any rank, starting from the well-known positivity of such form if $ (E,h) $ is just Griffiths positive of rank $ 2 $. The final part of the paper gives an overview on the state of the art of Griffiths' conjecture, collecting several remarks and open questions.


title: Operator learning for predicting multiscale bubble growth dynamics

abstract: Simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long term objective to simplify multiscale modeling by avoiding the fragile and time-consuming "hand-shaking" interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh-Plesset (R-P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.


title: Propagation of Strongly Hybridized Edge Modes in Plasmonic Waveguides

abstract: In our work we investigate the propagation of optical modes in nanoscale hybrid plasmonic waveguides. Frequency domain Maxwell equations based simulations are implemented to study properties of mixed modes in 3D. The results of our analysis are applied to design an advanced hybrid waveguide with useful characteristics. It is demonstrated here how to couple two edge plasmons with a waveguide mode to form a hybrid waveguide well tuned for applications in the field of electro-optical transistors (modulators). A mixed polarization state of edge plasmons allows them to couple with a horizontally polarized (TE) optical mode (which is not possible for usual surface plasmons) and open the way to polarization independent plasmonic modulators.


title: Permutation groups with restricted stabilizers

abstract: Fix a positive integer $d$ and let $\Gamma_d$ be the class of finite groups without sections isomorphic to the alternating group $A_d$. The groups in $\Gamma_d$ were studied by Babai, Cameron and P\'{a}lfy in the 1980s and they determined bounds on the order of a primitive permutation group with this property, which have found a wide range of applications. Subsequently, results on the base sizes of such groups were also obtained. In this paper we replace the structural conditions on the group by restrictions on its point stabilizers, and we obtain similar, and sometimes stronger conclusions. For example, we prove that there is a linear function $f$ such that the base size of any finite primitive group with point stabilizers in $\Gamma_d$ is at most $f(d)$. This generalizes a recent result of the first author on primitive groups with solvable point stabilizers. For non-affine primitive groups we obtain stronger results, assuming only that stabilizers of $c$ points lie in $\Gamma_d$. We also show that if $G$ is any permutation group of degree $n$ whose $c$-point stabilizers lie in $\Gamma_d$, then $|G| \leqslant ((1+o_c(1))d/e)^{n-1}$. This asymptotically extends and improves a $d^{n-1}$ upper bound on $|G|$ obtained by Babai, Cameron and P\'{a}lfy assuming $G \in \Gamma_d$.


title: Filamentary collapse flow in molecular clouds

abstract: We present idealized numerical simulations of prestellar gravitational collapse of a moderate initial filamentary perturbation with an additional central ellipsoidal enhancement (a core) considering a uniform, and a stratified background, the latter representing flattened clouds. Both simulations maintain the filamentary structure during the collapse, developing a hierarchical accretion flow from the cloud to the plane; from there to the filament, and from the filament to the core. The flow changes direction smoothly at every step of the hierarchy, with no density divergence nor a shock developing at the filament's axis during the studied prestellar evolution. The flow drives accretion onto the central core and drains material from the filament, slowing down the growth of the latter. As a consequence, the ratio of the central density of the core to the filament density increases in time, diverging at the time of singularity formation in the core. The stratified simulation produces the best match for observed Plummer-like radial column density profiles of filaments, while the uniform simulation does not produce a flat central density profile. This result supports recent suggestions that MCs may be preferentially flattened structures. We examine the possibility that the filamentary flow might approach a quasi-stationary regime in which the radial accretion onto the filament is balanced by the longitudinal accretion onto the core. A simple argument suggests that such a stationary state may be an attractor for the system. Our simulations, do not attain this stationary stage, but appear to be approaching it during the prestellar stage.


title: Multiclass Spinal Cord Tumor Segmentation on MRI with Deep Learning

abstract: Spinal cord tumors lead to neurological morbidity and mortality. Being able to obtain morphometric quantification (size, location, growth rate) of the tumor, edema, and cavity can result in improved monitoring and treatment planning. Such quantification requires the segmentation of these structures into three separate classes. However, manual segmentation of 3-dimensional structures is time-consuming and tedious, motivating the development of automated methods. Here, we tailor a model adapted to the spinal cord tumor segmentation task. Data were obtained from 343 patients using gadolinium-enhanced T1-weighted and T2-weighted MRI scans with cervical, thoracic, and/or lumbar coverage. The dataset includes the three most common intramedullary spinal cord tumor types: astrocytomas, ependymomas, and hemangioblastomas. The proposed approach is a cascaded architecture with U-Net-based models that segments tumors in a two-stage process: locate and label. The model first finds the spinal cord and generates bounding box coordinates. The images are cropped according to this output, leading to a reduced field of view, which mitigates class imbalance. The tumor is then segmented. The segmentation of the tumor, cavity, and edema (as a single class) reached 76.7 $\pm$ 1.5% of Dice score and the segmentation of tumors alone reached 61.8 $\pm$ 4.0% Dice score. The true positive detection rate was above 87% for tumor, edema, and cavity. To the best of our knowledge, this is the first fully automatic deep learning model for spinal cord tumor segmentation. The multiclass segmentation pipeline is available in the Spinal Cord Toolbox (https://spinalcordtoolbox.com/). It can be run with custom data on a regular computer within seconds.


title: Focal Frequency Loss for Image Reconstruction and Synthesis

abstract: Image reconstruction and synthesis have witnessed remarkable progress thanks to the development of generative models. Nonetheless, gaps could still exist between the real and generated images, especially in the frequency domain. In this study, we show that narrowing gaps in the frequency domain can ameliorate image reconstruction and synthesis quality further. We propose a novel focal frequency loss, which allows a model to adaptively focus on frequency components that are hard to synthesize by down-weighting the easy ones. This objective function is complementary to existing spatial losses, offering great impedance against the loss of important frequency information due to the inherent bias of neural networks. We demonstrate the versatility and effectiveness of focal frequency loss to improve popular models, such as VAE, pix2pix, and SPADE, in both perceptual quality and quantitative performance. We further show its potential on StyleGAN2.


title: A Multi-Epoch Study of Radio Continuum Emission from MassiveProtostars

abstract: We report the results of the Jansky Very Large Array (JVLA) observation of five massive protostars at 6 and 22.2 GHz. The aim of the study was to compare their current fluxes and positions with previous observations to search for evidence of variability. Most of the observed sources present the morphologies of a thermal core, hosting the protostar and exhibiting no proper motion, and associated non-thermal radio lobes that are characterised by proper motions and located away from the thermal core. Some of the protostars drive jets whose lobes have dissimilar displacement vectors, implying precession of the jets or the presence of multiple jet drivers. The jets of the protostars were found to have proper motions that lie in the range 170$ \leq v \leq $650 kms$^{-1}$, and precessions of periods 40$\leq p \leq$50 years and angles 2$\leq \alpha \leq$10$^\circ$, assuming that their velocities $v=$500 kms$^{-1}$. The core of one of the sources, S255 NIRS 3, which was in outburst at the time of our observations, showed a significant change in flux compared to the other sources. Its spectral index decreased during the outburst, consistent with the model of an expanding gas bubble. Modelling the emission of the outburst as that of a new non-thermal lobe that is emerging from a thermal core whose emission enshrouds that of the lobe also has the potential to account for the increase in flux and a decrease in the spectral index of the source's outburst.


title: High-speed generation of vector beams through random spatial multiplexing

abstract: Complex vector modes have become topical of late due to their fascinating properties and the many applications they have found across a broad variety of research fields. Even though such modes can be generated in a wide variety of ways, digital holography stands out as one of the most flexible and versatile. Along this line, Digital Micromirror Devices (DMDs) have gain popularity in recent time due to their high refresh rates, which allows the generation of vector modes at kHz rates. Nonetheless, most techniques are limited either by the diversity of vector modes that can be generated or by the speed at which they can be switched. Here we propose a technique based on the concept of random encoding, which allows the generation of arbitrary vector beams at speeds limited only by the refresh rate of the DMD. Our technique will be of great relevance in research fields such as optical communications, laser material processing and optical manipulation, amongst others.


title: Dark Energy Survey Year 3 Results: Deep Field Optical + Near-Infrared Images and Catalogue

abstract: We describe the Dark Energy Survey (DES) Deep Fields, a set of images and associated multi-wavelength catalogue ($ugrizJHKs$) built from Dark Energy Camera (DECam) and Visible and Infrared Survey Telescope for Astronomy (VISTA) data. The DES Deep Fields comprise 11 fields (10 DES supernova fields plus COSMOS), with a total area of $\sim30~$ square degrees in $ugriz$ bands and reaching a maximum $i$-band depth of 26.75 (AB, $10\sigma$, 2 arcsec). We present a catalogue for the DES 3-year cosmology analysis of those four fields with full 8-band coverage, totalling $5.88~$ sq. deg. after masking. Numbering $2.8~$million objects ($1.6~$million post masking), our catalogue is drawn from images coadded to consistent depths of $r=25.7, i=25, z=24.3$ mag. We use a new model-fitting code, built upon established methods, to deblend sources and ensure consistent colours across the $u$-band to $Ks$-band wavelength range. We further detail the tight control we maintain over the point-spread function modelling required for the model fitting, astrometry and consistency of photometry between the four fields. The catalogue allows us to perform a careful star-galaxy separation and produces excellent photometric redshift performance (${\rm NMAD} = 0.023$ at $i<23$). The Deep-Fields catalogue will be made available as part of the cosmology data products release, following the completion of the DES 3-year weak lensing and galaxy clustering cosmology work.


title: Dark Energy Survey Year 3 Results: Measuring the Survey Transfer Function with Balrog

abstract: We describe an updated calibration and diagnostic framework, Balrog, used to directly sample the selection and photometric biases of the Dark Energy Survey's (DES) Year 3 (Y3) dataset. We systematically inject onto the single-epoch images of a random 20% subset of the DES footprint an ensemble of nearly 30 million realistic galaxy models derived from DES Deep Field observations. These augmented images are analyzed in parallel with the original data to automatically inherit measurement systematics that are often too difficult to capture with traditional generative models. The resulting object catalog is a Monte Carlo sampling of the DES transfer function and is used as a powerful diagnostic and calibration tool for a variety of DES Y3 science, particularly for the calibration of the photometric redshifts of distant "source" galaxies and magnification biases of nearer "lens" galaxies. The recovered Balrog injections are shown to closely match the photometric property distributions of the Y3 GOLD catalog, particularly in color, and capture the number density fluctuations from observing conditions of the real data within 1% for a typical galaxy sample. We find that Y3 colors are extremely well calibrated, typically within ~1-8 millimagnitudes, but for a small subset of objects we detect significant magnitude biases correlated with large overestimates of the injected object size due to proximity effects and blending. We discuss approaches to extend the current methodology to capture more aspects of the transfer function and reach full coverage of the survey footprint for future analyses.


title: Dark Energy Survey Year 3 Results: Calibration of Lens Sample Redshift Distributions using Clustering Redshifts with BOSS/eBOSS

abstract: We present clustering redshift measurements for Dark Energy Survey (DES) lens sample galaxies to be used in weak gravitational lensing and galaxy clustering studies. To perform this measurement, we cross-correlate with spectroscopic galaxies from the Baryon Acoustic Oscillation Survey (BOSS) and its extension, eBOSS. We validate our methodology in simulations, including a new technique to calibrate systematic errors due to the galaxy clustering bias, finding our method to be generally unbiased in calibrating the mean redshift. We apply our method to the data, and estimate the redshift distribution for eleven different photometrically-selected bins. We find general agreement between clustering redshift and photometric redshift estimates, with differences on the inferred mean redshift to be below $|\Delta z|=0.01$ in most of the bins. We also test a method to calibrate a width parameter for redshift distributions, which we found necessary to use for some of our samples. Our typical uncertainties on the mean redshift ranged from 0.003 to 0.008, while our uncertainties on the width ranged from 4 to 9\%. We discuss how these results calibrate the photometric redshift distributions used in companion DES Year 3 Results papers.


title: 5d and 4d SCFTs: Canonical Singularities, Trinions and S-Dualities

abstract: Canonical threefold singularities in M-theory and Type IIB string theory give rise to superconformal field theories (SCFTs) in 5d and 4d, respectively. In this paper, we study canonical hypersurface singularities whose resolutions contain residual terminal singularities and/or 3-cycles. We focus on a certain class of `trinion' singularities which exhibit these properties. In Type IIB, they give rise to 4d $\mathcal{N}=2$ SCFTs that we call $D_p^b(G)$-trinions, which are marginal gaugings of three SCFTs with $G$ flavor symmetry. In order to understand the 5d physics of these trinion singularities in M-theory, we reduce these 4d and 5d SCFTs to 3d $\mathcal{N}=4$ theories, thus determining the electric and magnetic quivers (or, more generally, quiverines). In M-theory, residual terminal singularities give rise to free sectors of massless hypermultiplets, which often are discretely gauged. These free sectors appear as `ugly' components of the magnetic quiver of the 5d SCFT. The 3-cycles in the crepant resolution also give rise to free hypermultiplets, but their physics is more subtle, and their presence renders the magnetic quiver `bad'. We propose a way to redeem the badness of these quivers using a class $\mathcal{S}$ realization. We also discover new S-dualities between different $D_p^b(G)$-trinions. For instance, a certain $E_8$ gauging of the $E_8$ Minahan-Nemeschansky theory is S-dual to an $E_8$-shaped Lagrangian quiver SCFT.


title: Constructing Turing complete Euler flows in dimension $3$

abstract: Can every physical system simulate any Turing machine? This is a classical problem which is intimately connected with the undecidability of certain physical phenomena. Concerning fluid flows, Moore asked in [15] if hydrodynamics is capable of performing computations. More recently, Tao launched a programme based on the Turing completeness of the Euler equations to address the blow up problem in the Navier-Stokes equations. In this direction, the undecidability of some physical systems has been studied in recent years, from the quantum gap problem [7] to quantum field theories [11]. To the best of our knowledge, the existence of undecidable particle paths of 3D fluid flows has remained an elusive open problem since Moore's works in the early 1990's. In this article we construct a Turing complete stationary Euler flow on a Riemannian $S^3$ and speculate on its implications concerning Tao's approach to the blow up problem in the Navier-Stokes equations.


title: Effects of latitude-dependent gravity wave source variations on the middle and upper atmosphere

abstract: Atmospheric gravity waves (GWs) are generated in the lower atmosphere by various weather phenomena. They propagate upward, carry energy and momentum to higher altitudes, and appreciably influence the general circulation upon depositing them in the middle and upper atmosphere. We use a three-dimensional first-principle general circulation model (GCM) with an implemented nonlinear whole atmosphere GW parameterization to study the global climatology of wave activity and produced effects at altitudes up to the upper thermosphere. The numerical experiments were guided by the GW momentum fluxes and temperature variances as measured in 2010 by the SABER (Sounding of the Atmosphere using Broadband Emission Radiometry) instrument onboard NASA's TIMED (Thermosphere Ionosphere Mesosphere Energetics Dynamics) satellite. This includes the latitudinal dependence and magnitude of GW activity in the lower stratosphere for the boreal summer season. The modeling results were compared to the SABER temperature and total absolute momentum flux, and Upper Atmosphere Research Satellite (UARS) data in the mesosphere and lower thermosphere. Simulations suggest that, in order to reproduce the observed circulation and wave activity in the middle atmosphere, smaller than the measured GW fluxes have to be used at the source level in the lower atmosphere. This is because observations contain a broader spectrum of GWs, while parameterizations capture only a portion relevant to the middle and upper atmosphere dynamics. Accounting for the latitudinal variations of the source appreciably improves simulations.

categories: physics.space-ph astro-ph.EP physics.ao-ph physics.comp-ph physics.flu-dyn


title: Tractability frontiers in probabilistic team semantics and existential second-order logic over the reals

abstract: Probabilistic team semantics is a framework for logical analysis of probabilistic dependencies. Our focus is on the axiomatizability, complexity, and expressivity of probabilistic inclusion logic and its extensions. We identify a natural fragment of existential second-order logic with additive real arithmetic that captures exactly the expressivity of probabilistic inclusion logic. We furthermore relate these formalisms to linear programming, and doing so obtain PTIME data complexity for the logics. Moreover, on finite structures, we show that the full existential second-order logic with additive real arithmetic can only express NP properties. Lastly, we present a sound and complete axiomatization for probabilistic inclusion logic at the atomic level.


title: Approximation Limitations of Pure Dynamic Programming

abstract: We prove the first, even super-polynomial, lower bounds on the size of tropical (min,+) and (max,+) circuits approximating given optimization problems. Many classical dynamic programming (DP) algorithms for optimization problems are pure in that they only use the basic min, max, + operations in their recursion equations. Tropical circuits constitute a rigorous mathematical model for this class of algorithms. An algorithmic consequence of our lower bounds for tropical circuits is that the approximation powers of pure DP algorithms and greedy algorithms are incomparable. That pure DP algorithms can hardly beat greedy in approximation, is long known. New in this consequence is that also the converse holds.


title: Soft topological lattice wheels

abstract: We investigate the large deformation and extreme load-management capabilities of a soft topologically polarized kagome lattice mapped to a cylindrical domain through the problem of a lattice wheel rolling on an irregular surface. We test the surface-lattice interaction experimentally by subjecting a 3D-printed topological lattice wheel prototype to localized and distributed boundary loads. This investigation reveals a dichotomy in the force transfer between the two loading scenarios, whereby localized loads are absorbed with limited stress penetration into the bulk and small force transfer to the wheel axle, compared to distributed loads. Through numerical simulations, we compare the lattice wheel against a baseline solid wheel to highlight the unique stress management opportunities offered by the lattice configuration. These findings promote the design of rolling objects enabled by topological mechanics, in which a surplus of softness, activated by local asperities, can coexist with a globally stiff response to distributed loads that ensures satisfactory load-bearing capabilities.


title: Smoothness of cohomology sheaves of stacks of shtukas

abstract: We prove, for all reductive groups, that the cohomology sheaves with compact support of stacks of shtukas are ind-smooth over $(X - N)^I$ and that their geometric generic fibers are equiped with an action of $\text{Weil}(X - N, \overline{\eta})^I$. Our method does not use any compactification of stacks of shtukas.


title: General Theory of Large D Membranes Consistent with Second Law of Thermodynamics

abstract: We write down the most general membrane equations dual to black holes for a general class of gravity theories, up to sub-leading order in $1/D$ in large $D$ limit. We derive a "minimal" entropy current which satisfies a local form of second law from these membrane equations. We find that consistency with second law requires the membrane equations to satisfy certain constraints. We find additional constraints on the membrane equations from the existence of membrane solutions dual to stationary black holes. Finally we observe a tension between second law and matching with Wald entropy for dual stationary black hole configurations, for the minimal entropy current. We propose a simple modification of the membrane entropy current so that it satisfies second law and also the stationary membrane entropy matches the Wald entropy.


title: Enabling Secure and Effective Biomedical Data Sharing through Cyberinfrastructure Gateways

abstract: Dynaswap project reports on developing a coherently integrated and trustworthy holistic secure workflow protection architecture for cyberinfrastructures which can be used on virtual machines deployed through cyberinfrastructure (CI) services such as JetStream. This service creates a user-friendly cloud environment designed to give researchers access to interactive computing and data analysis resources on demand. The Dynaswap cybersecurity architecture supports roles, role hierarchies, and data hierarchies, as well as dynamic changes of roles and hierarchical relations within the scientific infrastructure. Dynaswap combines existing cutting-edge security frameworks (including an Authentication Authorization-Accounting framework, Multi-Factor Authentication, Secure Digital Provenance, and Blockchain) with advanced security tools (e.g., Biometric-Capsule, Cryptography-based Hierarchical Access Control, and Dual-level Key Management). The CI is being validated in life-science research environments and in the education settings of Health Informatics.


title: The potential role of binary neutron star merger afterglows in multimessenger cosmology

abstract: Binary neutron star mergers offer a new and independent means of measuring the Hubble constant $H_0$ by combining the gravitational-wave inferred source luminosity distance with its redshift obtained from electromagnetic follow-up. This method is limited by intrinsic degeneracy between the system distance and orbital inclination in the gravitational-wave signal. Observing the afterglow counterpart to a merger can further constrain the inclination angle, allowing this degeneracy to be partially lifted and improving the measurement of $H_0$. In the case of the binary neutron star merger GW170817, afterglow light-curve and imagery modeling thus allowed to improve the $H_0$ measurement by a factor of 3. However, systematic access to afterglow data is far from guaranteed. In fact, though each one allows a leap in $H_0$ precision, these afterglow counterparts should prove rare in forthcoming multimessenger campaigns. We combine models for emission and detection of gravitational-wave and electromagnetic radiation from binary neutron star mergers with realistic population models and estimates for afterglow inclination angle constraints. Using these models, we quantify how fast $H_0$ will be narrowed-down by successive multimessenger events with and without the afterglow. We find that, because of its rareness and though it greatly refines angle estimates, the afterglow counterpart should not significantly contribute to the measurement of $H_0$ in the long run.


title: Two-photon Above Threshold Ionization of Helium

abstract: Multiphoton ionization provides a clear window into the nature of electron correlations in the helium atom. In the present study, the final state energy range extends up to the region near the $N=2$ and $N=3$ ionization thresholds, where two-photon ionization proceeds via continuum intermediate states above the lowest threshold. Our calculations are performed using multichannel quantum defect theory (MQDT) and the streamlined R-matrix method. The sum and integration over all intermediate states in the two-photon ionization amplitude is evaluated using the inhomogeneous R-matrix method developed by Robicheaux and Gao. The seamless connection of that method with MQDT allows us to present high resolution spectra of the final state Rydberg resonances. Our analysis classifies the resonances above the $N=2$ threshold in terms of their group theory quantum numbers. Their dominant decay channels are found to obey the previously conjectured propensity rule far more weakly for these even parity states than was observed for the odd-parity states relevant to single photon ionization.


title: Sorting Can Exponentially Speed Up Pure Dynamic Programming

abstract: Many discrete minimization problems, including various versions of the shortest path problem, can be efficiently solved by dynamic programming (DP) algorithms that are "pure" in that they only perform basic operations, as min, max, +, but no conditional branchings via if-then-else in their recursion equations. It is known that any pure (min,+) DP algorithm solving the minimum weight spanning tree problem on undirected n-vertex graphs must perform at least $2^{\Omega(\sqrt{n})}$ operations. We show that this problem can be solved by a pure (min,max,+) DP algorithm performing only $O(n^3)$ operations. The algorithm is essentially a (min,max) algorithm: addition operations are only used to output the final values. The presence of both min and max operations means that now DP algorithms can sort: this explains the title of the paper.


title: Cohorting to isolate asymptomatic spreaders: An agent-based simulation study on the Mumbai Suburban Railway

abstract: The Mumbai Suburban Railways, \emph{locals}, are a key transit infrastructure of the city and is crucial for resuming normal economic activity. To reduce disease transmission, policymakers can enforce reduced crowding and mandate wearing of masks. \emph{Cohorting} -- forming groups of travelers that always travel together, is an additional policy to reduce disease transmission on \textit{locals} without severe restrictions. Cohorting allows us to: ($i$) form traveler bubbles, thereby decreasing the number of distinct interactions over time; ($ii$) potentially quarantine an entire cohort if a single case is detected, making contact tracing more efficient, and ($iii$) target cohorts for testing and early detection of symptomatic as well as asymptomatic cases. Studying impact of cohorts using compartmental models is challenging because of the ensuing representational complexity. Agent-based models provide a natural way to represent cohorts along with the representation of the cohort members with the larger social network. This paper describes a novel multi-scale agent-based model to study the impact of cohorting strategies on COVID-19 dynamics in Mumbai. We achieve this by modeling the Mumbai urban region using a detailed agent-based model comprising of 12.4 million agents. Individual cohorts and their inter-cohort interactions as they travel on locals are modeled using local mean field approximations. The resulting multi-scale model in conjunction with a detailed disease transmission and intervention simulator is used to assess various cohorting strategies. The results provide a quantitative trade-off between cohort size and its impact on disease dynamics and well being. The results show that cohorts can provide significant benefit in terms of reduced transmission without significantly impacting ridership and or economic \& social activity.


title: Existence of Kazdan-Warner equation with sign-changing prescribed function

abstract: In this paper, we study the following Kazdan-Warner equation with sign-changing prescribed function $h$ \begin{align*} -\Delta u=8\pi\left(\frac{he^{u}}{\int_{\Sigma}he^{u}}-1\right) \end{align*} on a closed Riemann surface whose area is equal to one. The solutions are the critical points of the functional $J_{8\pi}$ which is defined by \begin{align*} J_{8\pi}(u)=\frac{1}{16\pi}\int_{\Sigma}|\nabla u|^2+\int_{\Sigma}u-\ln\left|\int_{\Sigma}he^{u}\right|,\quad u\in H^1\left(\Sigma\right). \end{align*} We prove the existence of minimizer of $J_{8\pi}$ by assuming \begin{equation*} \Delta \ln h^++8\pi-2K>0 \end{equation*}at each maximum point of $2\ln h^++A$, where $K$ is the Gaussian curvature, $h^+$ is the positive part of $h$ and $A$ is the regular part of the Green function. This generalizes the existence result of Ding, Jost, Li and Wang [Asian J. Math. 1(1997), 230-248] to the sign-changing prescribed function case. We are also interested in the blow-up behavior of a sequence $u_{\varepsilon}$ of critical points of $J_{8\pi-\varepsilon}$ with $\int_{\Sigma}he^{u_{\varepsilon}}=1, \lim\limits_{\varepsilon\searrow 0}J_{8\pi-\varepsilon}\left(u_{\varepsilon}\right)<\infty$ and obtain the following identity during the blow-up process \begin{equation*} -\varepsilon=\frac{16\pi}{(8\pi-\varepsilon)h(p_\varepsilon)}\left[\Delta \ln h(p_\varepsilon)+8\pi-2K(p_\varepsilon)\right]\lambda_{\varepsilon}e^{-\lambda_{\varepsilon}}+O\left(e^{-\lambda_{\varepsilon}}\right), \end{equation*}where $p_\varepsilon$ and $\lambda_\varepsilon$ are the maximum point and maximum value of $u_\varepsilon$, respectively. Moreover, $p_{\varepsilon}$ converges to the blow-up point which is a critical point of the function $2\ln h^{+}+A$.


title: The L\'evy-Rosenzweig-Porter random matrix ensemble

abstract: In this paper we consider an extension of the Rosenzweig-Porter (RP) model, the L\'evy-RP (L-RP) model, in which the off-diagonal matrix elements are broadly distributed, providing a more realistic benchmark to develop an effective description of non-ergodic extended (NEE) states in interacting many-body disordered systems. We put forward a simple, general, and intuitive argument that allows one to unveil the multifractal structure of the mini-bands in the local spectrum when hybridization is due to anomalously large transition amplitudes in the tails of the distribution. The idea is that the energy spreading of the mini-bands can be determined self-consistently by requiring that the maximum of the matrix elements between a site $i$ and the other $N^{D_1}$ sites of the support set is of the same order of the Thouless energy itself $N^{D_1 - 1}$. This argument yields the fractal dimensions that characterize the statistics of the multifractal wave-functions in the NEE phase, as well as the whole phase diagram of the L-RP ensemble. Its predictions are confirmed both analytically, by a thorough investigation of the self-consistent equation for the local density of states obtained using the cavity approach, and numerically, via extensive exact diagonalizations.


title: Point-Based Value Iteration and Approximately Optimal Dynamic Sensor Selection for Linear-Gaussian Processes

abstract: The problem of synthesizing an optimal sensor selection policy is pertinent to a variety of engineering applications ranging from event detection to autonomous navigation. We consider such a synthesis problem over an infinite time horizon with a discounted cost criterion. We formulate this problem in terms of a value iteration over the continuous space of covariance matrices. To obtain a computationally tractable solution, we subsequently formulate an approximate sensor selection problem, which is solvable through a point-based value iteration over a finite "mesh" of covariance matrices with a user-defined bounded trace. We provide theoretical guarantees bounding the suboptimality of the sensor selection policies synthesized through this method and provide numerical examples comparing them to known results.


title: EQ-Net: A Unified Deep Learning Framework for Log-Likelihood Ratio Estimation and Quantization

abstract: In this work, we introduce EQ-Net: the first holistic framework that solves both the tasks of log-likelihood ratio (LLR) estimation and quantization using a data-driven method. We motivate our approach with theoretical insights on two practical estimation algorithms at the ends of the complexity spectrum and reveal a connection between the complexity of an algorithm and the information bottleneck method: simpler algorithms admit smaller bottlenecks when representing their solution. This motivates us to propose a two-stage algorithm that uses LLR compression as a pretext task for estimation and is focused on low-latency, high-performance implementations via deep neural networks. We carry out extensive experimental evaluation and demonstrate that our single architecture achieves state-of-the-art results on both tasks when compared to previous methods, with gains in quantization efficiency as high as $20\%$ and reduced estimation latency by up to $60\%$ when measured on general purpose and graphical processing units (GPU). In particular, our approach reduces the GPU inference latency by more than two times in several multiple-input multiple-output (MIMO) configurations. Finally, we demonstrate that our scheme is robust to distributional shifts and retains a significant part of its performance when evaluated on 5G channel models, as well as channel estimation errors.


title: Automating Surgical Peg Transfer: Calibration with Deep Learning Can Exceed Speed, Accuracy, and Consistency of Humans

abstract: Peg transfer is a well-known surgical training task in the Fundamentals of Laparoscopic Surgery (FLS). While human sur-geons teleoperate robots such as the da Vinci to perform this task with high speed and accuracy, it is challenging to automate. This paper presents a novel system and control method using a da Vinci Research Kit (dVRK) surgical robot and a Zivid depth sensor, and a human subjects study comparing performance on three variants of the peg-transfer task: unilateral, bilateral without handovers, and bilateral with handovers. The system combines 3D printing, depth sensing, and deep learning for calibration with a new analytic inverse kinematics model and a time-minimized motion controller. In a controlled study of 3384 peg transfer trials performed by the system, an expert surgical resident, and 9 volunteers, results suggest that the system achieves accuracy on par with the experienced surgical resident and is significantly faster and more consistent than the surgical resident and volunteers. The system also exhibits the highest consistency and lowest collision rate. To our knowledge, this is the first autonomous system to achieve superhuman performance on a standardized surgical task.


title: Searching optimal conditions for quantum gates application with the new 3-body F\"orster resonances in Rb and Cs Rydberg atoms

abstract: Three body resonant interactions between Rydberg atoms are considered in order to perform few-body quantum gates. So far, the resonances found in cesium or rubidium atoms relied on an adjacent two-body resonance which ceases to exist for principal quantum numbers above $n \simeq 40$. We have proposed recently a new class of 3-body interaction resonances in alkali-metal Rydberg atoms [P. Cheinet \textit{et al.}, Quant. Elect. \textbf{50}, 213 (2020)], which circumvienes this limit. We investigate here the relative strength between this new class of 3-body interaction resonance and quasi-forbidden 2-body interaction resonances in rubidium and cesium Rydberg atoms. We then identify the best case scenario for detecting and using this 3-body interaction.


title: Improved Algorithm for Computing the Maximum-volume Bichromatic Separating Box

abstract: We consider the problem of computing the largest-area bichromatic separating box among a set of $n$ red points and a set of $m$ blue points in three dimensions. Currently, the best-known algorithm to solve this problem takes $O(m^2 (m + n))$ time and $O(m + n)$ space. In this paper, we come up with an improved algorithm for the problem, which takes $O(m^2 + n)$ time.


title: Dark Matter and Nature of Electroweak Phase Transition with an Inert Doublet

abstract: We provide a comprehensive and up-to-date analysis of the prospects to realize Dark Matter (DM) in the Inert Doublet Model, while simultaneously enhancing the Electroweak Phase Transition (EWPhT) such as to allow for electroweak baryogenesis. Instead of focusing on certain aspects or mass hierarchies, we perform extensive, yet fine-grained, parameter space scans, where we analyze the nature of the EWPhT in both the light and the heavy DM regions, confronting it with the amount of DM potentially residing in the lightest inert-doublet state. Thereby, we point out a viable region where a non-trivial two-step EWPhT can appear, without being in conflict with direct-detection bounds, which could leave interesting imprints in gravitational wave signatures. We propose new benchmarks with this feature as well as update benchmarks with a strong first-order transition in the light of new XENON1T limits. Moreover, taking into account these latest bounds as well as relevant collider constraints, we envisage a new region for light DM with a small mass splitting, lifting the usual assumption of exact degeneracy of the new non-DM scalars, such as to avoid collider bounds while providing a fair DM abundance over a rather large DM mass range. This follows from a detailed survey of the impact of co-annihilations on the abundance, dissecting the various channels.


title: R\'enyi free energy and variational approximations to thermal states

abstract: We propose the construction of thermodynamic ensembles that minimize the R\'enyi free energy, as an alternative to Gibbs states. For large systems, the local properties of these R\'enyi ensembles coincide with those of thermal equilibrium, and they can be used as approximations to thermal states. We provide algorithms to find tensor network approximations to the 2-R\'enyi ensemble. In particular, a matrix-product-state representation can be found by using gradient-based optimization on Riemannian manifolds, or via a non-linear evolution which yields the desired state as a fixed point. We analyze the performance of the algorithms and the properties of the ensembles on one-dimensional spin chains.


title: Exceptional points of degeneracy in traveling wave tubes

abstract: Traveling wave tube (TWT) is a powerful vacuum electronic device used to amplify radio-frequency (RF) signals with numerous applications, including radar, television and telephone satellite communications. TWT design in a nutshell comprises of a pencil-like electron beam (e-beam) in vacuum interacting with guiding it slow-wave structure (SWS). In our studies here the e-beam is represented by one-dimensional electron flow and SWS is represented by a transmission line (TL). The interaction between the e-beam and the TL is modeled by an analytic theory that generalizes the well-known Pierce model by taking into account the so-called space-charge effects particularly electron-to-electron repulsion (debunching). Many important aspects of the analytic theory of TWTs have been already analyzed in our monograph on the subject. The focus of the studies here is on degeneracies of the TWT dispersion relations particularly on exceptional points of degeneracy and their applications. The term exceptional point of degeneracy (EPD) refers to the property of the relevant matrix to have nontrivial Jordan block structure. Using special parameterization particularly suited to chosen EPD we derive exact formulas for the relevant Jordan basis including the eigenvectors and the so-called root vector associated with the Jordan block. Based on these studies we develop constructive approach to sensing of small signals.


title: Homochirality: a prerequisite or consequence of life?

abstract: Many of the building blocks of life such as amino acids and nucleotides are chiral, i.e., different from their mirror image. Contemporary life selects and synthesizes only one of two possible handednesses. In an abiotic environment, however, there are usually equally many left- and right-handed molecules. If homochirality was a prerequisite of life, there must have been physical or chemical circumstances that led to the selection of a certain preference. Conversely, if it was a consequence of life, we must identify possible pathways for accomplishing a transition from a racemic to a homochiral chemistry. After a discussion of the observational evidence, I will review ideas where homochirality of any handedness could emerge as a consequence of the first polymerization events of nucleotides in an emerging RNA world. These mechanisms are not limited to nucleotides, but can also occur for peptides, as a precursor to the RNA world. The question of homochirality is, in this sense, intimately tied to the origin of life. Future Mars missions may be able to detect biomolecules of extant or extinct life. We will therefore also discuss possible experimental setups for determining the chirality of primitive life forms in situ on Mars.


title: A note on Bridgeland stability conditions and Catalan numbers

abstract: In this short note, we describe a problem in algebraic geometry where the solution involves Catalan numbers. More specifically, we consider the derived category of coherent sheaves on an elliptic surface, and the action of its autoequivalence group on its Bridgeland stability manifold. In solving an equation involving this group action, the generating function of Catalan numbers arises, allowing us to use asymptotic estimates of Catalan numbers to arrive at a bound for the solution set.


title: New aspects of Argyres--Douglas theories and their dimensional reduction

abstract: Argyres-Douglas (AD) theories constitute an infinite class of superconformal field theories in four dimensions with a number of interesting properties. We study several new aspects of AD theories engineered in $A$-type class $\mathcal{S}$ with one irregular puncture of Type I or Type II and also a regular puncture. These include conformal manifolds, structures of the Higgs branch, as well as the three dimensional gauge theories coming from the reduction on a circle. The latter admit a description in terms of a linear quiver with unitary and special unitary gauge groups, along with a number of twisted hypermultiplets. The origin of these twisted hypermultiplets is explained from the four dimensional perspective. We also propose the three dimensional mirror theories for such linear quivers. These provide explicit descriptions of the magnetic quivers of all AD theories in question in terms of quiver diagrams with unitary gauge groups, together with a collection of free hypermultiplets. A number of quiver gauge theories presented in this paper are new and have not been studied elsewhere in the literature.


title: Explaining cosmic ray antimatter with secondaries from old supernova remnants

abstract: Despite significant efforts over the past decade, the origin of the cosmic ray positron excess has still not been unambiguously established. A popular class of candidate sources are pulsars or pulsar wind nebulae but these cannot also account for the observed hard spectrum of cosmic ray antiprotons. We revisit the alternative possibility that the observed high-energy positrons are secondaries created by spallation in supernova remnants during the diffusive shock acceleration of the primary cosmic rays, which are further accelerated by the same shocks. The resulting source spectrum of positrons at high energies is then naturally harder than that of the primaries, as is the spectrum of other secondaries such as antiprotons. We present the first comprehensive investigation of the full parameter space of this model -- both the source parameters as well as those governing galactic transport. Various parameterisations of the cross sections for the production of positrons and antiprotons are considered, and the uncertainty in the model parameters discussed. We obtain an excellent fit to recent precision measurements by AMS-02 of cosmic ray protons, helium, positrons and antiprotons, as well as of various primary and secondary nuclei. This model thus provides an economical explanation of the spectra of all secondary species -- from a single well-motivated population of sources.


title: Deep manifold learning reveals hidden dynamics of proteasome autoregulation

abstract: The 2.5-MDa 26S proteasome maintains proteostasis and regulates myriad cellular processes. How polyubiquitylated substrate interactions regulate proteasome activity is not understood. Here we introduce a deep manifold learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy (cryo-EM) reconstructions of nonequilibrium conformational continuum and reconstitutes hidden dynamics of proteasome autoregulation in the act of substrate degradation. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of free-energy landscapes, which directs 3D clustering via an energy-based particle-voting algorithm. In blind assessments using simulated heterogeneous cryo-EM datasets, AlphaCryo4D achieved 3D classification accuracy three times that of conventional method and reconstructed continuous conformational changes of a 130-kDa protein at sub-3-angstrom resolution. By using AlphaCryo4D to analyze a single experimental cryo-EM dataset, we identified 64 conformers of the substrate-bound human 26S proteasome, revealing conformational entanglement of two regulatory particles in the doubly capped holoenzymes and their energetic differences with singly capped ones. Novel ubiquitin-binding sites are discovered on the RPN2, RPN10 and Alpha5 subunits to remodel polyubiquitin chains for deubiquitylation and recycle. Importantly, AlphaCryo4D choreographs single-nucleotide-exchange dynamics of proteasomal AAA-ATPase motor during translocation initiation, which upregulates proteolytic activity by allosterically promoting nucleophilic attack. Our systemic analysis illuminates a grand hierarchical allostery for proteasome autoregulation.


title: Two-Dimensional Multifunctional Materials from Endohedral Fullerenes

abstract: A new multifunctional 2D material is theoretically predicted based on systematic ab-initio calculations and model simulations for the honeycomb lattice of endohedral fullerene W@C28 molecules. It has structural bistability, ferroelectricity, multiple magnetic phases, and excellent valley characters and can be easily functionalized by the proximity effect with magnetic isolators such as MnTiO3. Furthermore, we may also manipulate the valley Hall and spin transport properties by selectively switch a few W@C28 molecules to the metastable phase. These findings pave a new way in integrate different functions in a single 2D material for technological innovations.


title: Adaptive Mirror Descent Methods for Convex Programming Problems with delta-subgradients

abstract: We propose some adaptive mirror descent dethods for convex programming problems with delta-subgradients and prove some theoretical results.


title: On the extension of Muckenhoupt weights in metric spaces

abstract: A theorem by Wolff states that weights defined on a measurable subset of $\mathbb{R}^n$ and satisfying a Muckenhoupt-type condition can be extended into the whole space as Muckenhoupt weights of the same class. We give a complete and self-contained proof of this theorem generalized into metric measure spaces supporting a doubling measure. Related to the extension problem, we also show estimates for Muckenhoupt weights on Whitney chains in the metric setting.


title: Landau equation for self-gravitating classical and quantum particles: Application to dark matter

abstract: We develop the kinetic theory of classical and quantum particles (fermions and bosons) in gravitational interaction. The kinetic theory of quantum particles may have applications in the context of dark matter. For simplicity, we consider an infinite and spatially homogeneous system (or make a local approximation) and neglect collective effects. This leads to the quantum Landau equation derived heuristically in [Chavanis, Physica A 332, 89 (2004)]. We establish its main properties: conservation laws, $H$-theorem, equilibrium state, relaxation time, quantum diffusion and friction coefficients, quantum Rosenbluth potentials, self-consistent evolution, (thermal) bath approximation, quantum Fokker-Planck equation, quantum King model... For bosonic particles, the Landau equation can describe the process of Bose-Einstein condensation. We discuss the relation of our study with the works of [Levkov et al., Phys. Rev. Lett. 121, 151301 (2018); Bar-Or et al., Astrophys. J. 871, 28 (2019)] on fuzzy dark matter halos and the formation of Bose stars and solitons.


title: Limit Theorems for Fr\'echet Mean Sets

abstract: For $1\le p \le \infty$, the Fr\'echet $p$-mean of a probability measure on a metric space is an important notion of central tendency that generalizes the usual notions in the real line of mean ($p=2$) and median ($p=1$). In this work we prove a collection of limit theorems for Fr\'echet means and related objects, which, in general, constitute a sequence of random closed sets. On the one hand, we show that many limit theorems (a strong law of large numbers, an ergodic theorem, and a large deviations principle) can be simply descended from analogous theorems on the space of probability measures via purely topological considerations. On the other hand, we provide the first sufficient conditions for the strong law of large numbers to hold in a $T_2$ topology (in particular, the Fell topology), and we show that this condition is necessary in some special cases. We also discuss statistical and computational implications of the results herein.


title: On weighted $L^p$-Hardy inequality on domains in $\mathbb{R}^n$

abstract: We consider weighted $L^p$-Hardy inequalities involving the distance to the boundary of a domain in the $n$-dimensional Euclidean space with nonempty boundary. Using criticality theory, we give an alternative proof of the following result of F.~G.~Avkhadiev (2006) Theorem: Let $\Omega \subsetneqq \mathbb{R}^n$, $n\geq 2$, be an arbitrary domain, $1<p<\infty$ and $\alpha + p>n$. Let $\mathrm{d}_\Omega(x) =\mathrm{dist}(x,\partial \Omega )$ denote the distance of a point $x\in \Omega$ to $\partial \Omega$. Then the following Hardy-type inequality holds $$ \int_{\Omega }\frac{|\nabla \varphi |^p}{\mathrm{d}_\Omega^{\alpha}}\,\mathrm{d}x \geq \left( \frac{\alpha +p-n}{p}\right)^p \int_{\Omega }\frac{|\varphi|^p}{\mathrm{d}_\Omega^{p+\alpha}}\,\mathrm{d}x \qquad \forall \varphi\in C^{\infty }_c(\Omega),$$ and the lower bound constant $\left( \frac{\alpha +p-n}{p}\right)^p$ is sharp.


title: Credit Freezes, Equilibrium Multiplicity, and Optimal Bailouts in Financial Networks

abstract: We analyze how interdependencies between organizations in financial networks can lead to multiple possible equilibrium outcomes. A multiplicity arises if and only if there exists a certain type of dependency cycle in the network that allows for self-fulfilling chains of defaults. We provide necessary and sufficient conditions for banks' solvency in any equilibrium. Building on these conditions, we characterize the minimum bailout payments needed to ensure systemic solvency, as well as how solvency can be ensured by guaranteeing a specific set of debt payments. Bailout injections needed to eliminate self-fulfilling cycles of defaults (credit freezes) are fully recoverable, while those needed to prevent cascading defaults outside of cycles are not. We show that the minimum bailout problem is computationally hard, but provide an upper bound on optimal payments and show that the problem has intuitive solutions in specific network structures such as those with disjoint cycles or a core-periphery structure.


title: Towards Fair Personalization by Avoiding Feedback Loops

abstract: Self-reinforcing feedback loops are both cause and effect of over and/or under-presentation of some content in interactive recommender systems. This leads to erroneous user preference estimates, namely, overestimation of over-presented content while violating the right to be presented of each alternative, contrary of which we define as a fair system. We consider two models that explicitly incorporate, or ignore the systematic and limited exposure to alternatives. By simulations, we demonstrate that ignoring the systematic presentations overestimates promoted options and underestimates censored alternatives. Simply conditioning on the limited exposure is a remedy for these biases.


title: A Phase Field Crystal Method for Multilayer Graphene Structure

abstract: Bilayer graphene has been a subject of intense study in recent years. We extend a structural phase field crystal method to include an external potential from adjacent layer(s), which is generated by the corresponding phase field and changes over time. Moreover, multiple layers can be added into the structure. Using the thickness of the boundaries between different stacking variants of the bilayer structure as the key parameter, we quantify the strength of the adjacent layer potential by comparing with atomistic simulation results. We then test the multiple graphene structures, including bilayers, triple layers, up to 6 layers. We find that besides the initial conditions, the way of new layers added into the structure will also affect the layout of the atomic configuration. We believe tour results can help understanding the mechanism of graphene structure consists of more than one layer.


title: High-Entropy Rare Earth Tetraborides

abstract: Six high-entropy rare earth tetraborides of the tetragonal UB4-prototyped structure have been successfully synthesized for the first time. The specimens are prepared from elemental precursors via high-energy ball mill and in-situ reactive spark plasma sintering. The sintered specimens are >98% in relative densities without detectable oxide impurities (albeit the presence of minor hexaborides in some compositions). No detectable secondary phase is observed in the composition (Y$_{0.2}$Nd$_{0.2}$Sm$_{0.2}$Gd$_{0.2}$Tb$_{0.2}$)B$_{4}$, which is proven homogeneous at both microscale and nanoscale. The Vickers microhardness are determined to be ~13-15 GPa at a standard indentation load of 9.8 N. A scientifically interesting observation is represented by the anisotropic lattice distortion from the rule-of-mixture averages. This work expands the family of high-entropy ceramics via fabricating a new class of high-entropy borides with a unique tetragonal quasi-layered crystal structure.


title: Decoding Grain Boundary Thermodynamics in High-Entropy Alloys in a 5D Space: Coupled Segregation and Disordering

abstract: Grain boundaries (GBs) can critically influence the microstructural evolution and various materials properties. However, a fundamental understanding of GBs in high-entropy alloys (HEAs) is lacking because of the complex couplings of the segregations of multiple elements and interfacial disordering, which can generate new phenomena and challenge the classical theories. Here, by combining large-scale atomistic simulations and machine learning, we demonstrate the feasibility of predicting the GB properties as functions of four independent compositional degrees of freedoms and temperature in a 5D space. Subsequently, GB counterparts to bulk phase diagrams are constructed for the first time for quinary HEAs. A data-driven discovery further reveals new coupled segregation and disordering effects in HEAs. Notably, an analysis of a large dataset discovers a critical compensation temperature at which the segregations of all elements virtually vanish simultaneously. While the machine learning model can predict GB properties via a black-box approach, a surrogate data-based analytical model (DBAM) is constructed to provide more physics insights and better transferability, with good accuracies. This study not only provides a new paradigm enabling prediction of GB properties in a 5D space, but also uncovers new GB segregation phenomena in HEAs beyond the classical GB segregation models.


title: Single-Phase Duodenary High-Entropy Fluorite/Pyrochlore Oxides with an Order-Disorder Transition

abstract: Improved thermomechanical properties have been reported for various high-entropy oxides containing typically five metal cations. This study further investigates a series of duodenary (11 metals + oxygen) high-entropy oxides by mixing different fractions of a five-cation fluorite-structured niobate and a seven-cation pyrochlore (both containing Yb) with matching lattice constants. Nine compositions of duodenary high-entropy oxides have been examined. All of them exhibit single high-entropy phases of either disordered fluorite or ordered pyrochlore structure. An order-disorder transition (ODT) is evident with changing composition, accompanied by a reduction in thermal conductivity (k). In comparison with the ODT criteria developed from ternary oxides, these duodenary oxides are more prone to disorder, but the ODT is still controlled by similar factors (but at different thresholds). Interestingly, there are abrupt increases in Young's modulus (E) at low mixing concentrations near both endmembers. The E/k ratios are increased, in comparison with both endmembers. This study suggests a new route to tailor high-entropy ceramics via controlling cation ordering vs. disordering.


title: Integrity bases for cubic nonlinear magnetostriction

abstract: A so-called smart material is a material that is the seat of one or more multiphysical coupling. One of the key points in the development of the constitutive laws of these materials, either at the local or at the global scale, is to formulate a free energy density (or enthalpy) from vectors, tensors, at a given order and for a class of given symmetry, depending on the symmetry classes of the crystal constituting the material or the symmetry of the representative volume element. This article takes as a support of study the stress and magnetization couple ($\sigma$, m) involved in the phenomena of magnetoelastic coupling in a cubic symmetry medium. Several studies indeed show a non-monotonic sensitivity of the magnetic susceptibility and magnetostriction of certain soft magnetic materials under stress. Modeling such a phenomenon requires the introduction of a second order stress term in the Gibbs free energy density. A polynomial formulation in the two variables stress and magnetization is preferred over a tensorial formulation. For a given material symmetry class, this allows to express more easily the free energy density at any bi-degree in $\sigma$ and m (i.e. at any constitutive tensors order for the so-called tensorial formulation). A rigorous and systematic method is essential to obtain the high-degree magneto-mechanical coupling terms and to build a free energy density function at any order which is invariant by the action of the cubic (octahedral) group. For that aim, theoretical and computer tools in Invariant Theory, that allow for the mathematical description of cubic nonlinear magneto-elasticity, are introduced. Minimal integrity bases of the invariant algebra for the pair (m, $\sigma$), under the proper (orientation-preserving) and the full cubic groups, are then proposed. The minimal integrity basis for the proper cubic group is constituted of 60 invariants, while the minimal integrity basis for the full cubic group (the one of interest for magneto-elasticity) is made up of 30 invariants. These invariants are formulated in a (coordinate free) intrinsic manner, using a generalized cross product to write some of them. The counting of independent invariants of a given multi-degree in (m, $\sigma$) is performed. It is shown accordingly that it is possible to list without error all the material parameters useful for the description of the coupled magnetoelastic behavior from the integrity basis. The technique is applied to derive general expressions $\Psi$ $\star$ ($\sigma$, m) of the free energy density at the magnetic domains scale exhibiting cubic symmetry. The classic results for an isotropic medium are recovered.


title: Flat-Combining-Based Persistent Data Structures for Non-Volatile Memory

abstract: Flat combining (FC) is a synchronization paradigm in which a single thread, holding a global lock, collects requests by multiple threads for accessing a concurrent data structure and applies their combined requests to it. Although FC is sequential, it significantly reduces synchronization overheads and cache invalidations and thus often provides better performance than that of lock-free implementations. The recent emergence of non-volatile memory (NVM) technologies increases the interest in the development of persistent concurrent objects. These are objects that are able to recover from system failures and ensure consistency by retaining their state in NVM and fixing it, if required, upon recovery. Of particular interest are detectable objects that, in addition to ensuring consistency, allow recovery code to infer if a failed operation took effect before the crash and, if it did, obtain its response. In this work, we present the first FC-based persistent object implementations. Specifically, we introduce a detectable FC-based implementation of a concurrent LIFO stack, a concurrent FIFO queue, and a double-ended queue. Our empirical evaluation establishes that due to flat combining, the novel implementations require a much smaller number of costly persistence instructions than competing algorithms and are therefore able to significantly outperform them.


title: 3d Convex Contact Forms And The Ruelle Invariant

abstract: Let $X \subset \mathbb{R}^4$ be a convex domain with smooth boundary $Y$. We use a relation between the extrinsic curvature of $Y$ and the Ruelle invariant $\text{Ru}(Y)$ of the natural Reeb flow on $Y$ to prove that there exist constants $C > c > 0$ independent of $Y$ such that \[c < \frac{\text{Ru}(Y)^2}{\text{vol}(X)} \cdot \text{sys}(Y) < C\] Here $\text{sys}(Y)$ is the systolic ratio, i.e. the square of the minimal period of a closed Reeb orbit of $Y$ divided by twice the volume of $X$. We then construct dynamically convex contact forms on $S^3$ that violate this bound using methods of Abbondandolo-Bramham-Hryniewicz-Salom\~{a}o. These are the first examples of dynamically convex contact $3$-spheres that are not strictly contactomorphic to a convex boundary $Y$.


title: Lunar neutrinos

abstract: Cosmic rays bombard the lunar surface producing mesons, which attenuate inside the regolith. They get slower and decay weakly into mostly sub-GeV neutrinos leaving the surface. Thus the Moon shines in neutrinos. Here we calculate spectra of low energy neutrinos, which exhibit bright features potentially recognisable above isotropic neutrino background in the direction towards the Moon. Their observation, though a very challenging task for future neutrino large volume experiments, would make the Moon the nearest astrophysical source for which the concept of multimessenger astronomy works truly. Remarkably, some features of the lunar neutrino flux are sensitive to the surface mass density of the Moon.


title: A Multimodal Framework for the Detection of Hateful Memes

abstract: An increasingly common expression of online hate speech is multimodal in nature and comes in the form of memes. Designing systems to automatically detect hateful content is of paramount importance if we are to mitigate its undesirable effects on the society at large. The detection of multimodal hate speech is an intrinsically difficult and open problem: memes convey a message using both images and text and, hence, require multimodal reasoning and joint visual and language understanding. In this work, we seek to advance this line of research and develop a multimodal framework for the detection of hateful memes. We improve the performance of existing multimodal approaches beyond simple fine-tuning and, among others, show the effectiveness of upsampling of contrastive examples to encourage multimodality and ensemble learning based on cross-validation to improve robustness. We furthermore analyze model misclassifications and discuss a number of hypothesis-driven augmentations and their effects on performance, presenting important implications for future research in the field. Our best approach comprises an ensemble of UNITER-based models and achieves an AUROC score of 80.53, placing us 4th on phase 2 of the 2020 Hateful Memes Challenge organized by Facebook.


title: Energy efficient spark plasma sintering: Breaking the threshold of large dimension tooling energy consumption

abstract: An energy efficient spark plasma sintering method enabling the densification of large size samples assisted by very low electric current levels is developed. In this method, the electric current is concentrated in the graphite foils around the sample. High energy dissipation is then achieved in this area enabling the heating and full densification of large (alumina) parts ({{\O}} 40 mm) at relatively low currents (800 A). The electrothermal mechanical simulation reveals that the electric current needed to heat the large samples is 70 % lower in the energy efficient configuration compared to the traditional configuration. The presence of thermal and densification gradients is also revealed for the larger size samples. Potential solutions for this problem are discussed. The experiments confirm the possibility of full densification (96-99 %) of large alumina samples. This approach allows using small (and low cost) SPS devices (generally limited to 10-15 mm samples) for large size samples (40-50 mm). The developed technique enables also an optimized energy consumption by large scale SPS systems.


title: Actuation of cylindrical nematic elastomer balloons

abstract: Nematic elastomers are programmable soft materials that display large, reversible and predictable deformation under an external stimulus such as a change in temperature or light. While much of the work in the field has focused on actuation from flat sheets, recent advances in 3D printing and other methods of directed synthesis have motivated the study of actuation of curved shells. Snap-through buckling has been a topic of particular interest. In this work, we present theoretical calculations to motivate another mode of actuation that combines programmable soft materials as well as instabilities associated with large deformation. Specifically, we analyze the deformation of a cylindrical shell of a patterned nematic elastomer under pressure, show that it can undergo an enormous change of volume with changing temperature and suggest its application as a pump with extremely high ejection fraction.


title: The free flight step analysis in the Monte Carlo Method. Test case: charge transport in graphene

abstract: Monte Carlo methods are now a definitively well established methods for studying the transport of charges in semiconductors. A lot of drawbacks are however present about the correct reconstruction of the distribution function, which has to be bounded between $0$ and $1$. The standard Ensemble Monte Carlo method is not able to properly include the Pauli exclusion principle and the numerical distribution functions exceed the maximum value. A lot of attempts were made to solve the problem and recently a new DSMC scheme has been developed to overcome the problem. All the improvements are generally based on some approximations or computational techniques involving the distribution function and the scattering terms. Only an attempt was devoted to opportunely change the formulation of the free flight step but there is still a lack in the literature regarding a both theoretical and computational analysis of corrected Monte Carlo methods based on the free flights. In this paper we would like to fulfill such a gap by analyzing deeper the role of the free flight step into a Monte Carlo simulation, its physical properties and the conditions that allow us to preserve the physical consistency of the method.


title: Post-Newtonian limit of generalized symmetric teleparallel gravity

abstract: In this article we analyze the post-Newtonian approximation of a generalization of the symmetric teleparallel gravity with the help of the parameterized post-Newtonian (PPN) formalism. This class of theories is based on a free function of the five independent quadratic contractions of the non-metricity tensor. By calculating the PPN metric of these theories, we can restrict the Taylor coefficients of the free function with the help of the PPN parameters and their observational bounds. We find two families of theories whose PPN parameters are identical to those of general relativity, and thus in full agreement with observations. For three further families, we find that only the PPN parameters $\beta$ and $\gamma$ deviate, but can be brought arbitrarily close to their general relativity values by an appropriate choice of the Lagrangian, so that also these families contain candidate theories which agree with observations. The remaining theories either possess no well-defined solution of the post-Newtonian field equations, or possess a post-Newtonian limit which exceeds the form assumed in the PPN formalism.


title: Arnol'd Tongues in Oscillator Systems with Nonuniform Spatial Driving

abstract: Nonlinear oscillator systems are ubiquitous in biology and physics, and their control is a practical problem in many experimental systems. Here we study this problem in the context of the two models of spatially-coupled oscillators: the complex Ginzburg-Landau equation (CGLE) and a generalization of the CGLE in which oscillators are coupled through an external medium (emCGLE). We focus on external control drives that vary in both space and time. We find that the spatial distribution of the drive signal controls the frequency ranges over which oscillators synchronize to the drive and that boundary conditions strongly influence synchronization to external drives for the CGLE. Our calculations also show that the emCGLE has a low density regime in which a broad range of frequencies can be synchronized for low drive amplitudes. We study the bifurcation structure of these models and find that they are very similar to results for the driven Kuramoto model, a system with no spatial structure. We conclude by discussing the implications of our results for controlling coupled oscillator systems such as the social amoebae \emph{Dictyostelium} and populations of BZ catalytic particles using spatially structured external drives.


title: Training data-efficient image transformers & distillation through attention

abstract: Recently, neural networks purely based on attention were shown to address image understanding tasks such as image classification. However, these visual transformers are pre-trained with hundreds of millions of images using an expensive infrastructure, thereby limiting their adoption. In this work, we produce a competitive convolution-free transformer by training on Imagenet only. We train them on a single computer in less than 3 days. Our reference vision transformer (86M parameters) achieves top-1 accuracy of 83.1% (single-crop evaluation) on ImageNet with no external data. More importantly, we introduce a teacher-student strategy specific to transformers. It relies on a distillation token ensuring that the student learns from the teacher through attention. We show the interest of this token-based distillation, especially when using a convnet as a teacher. This leads us to report results competitive with convnets for both Imagenet (where we obtain up to 85.2% accuracy) and when transferring to other tasks. We share our code and models.


title: Resetting with stochastic return through linear confining potential

abstract: We consider motion of an overdamped Brownian particle subject to stochastic resetting in one dimension. In contrast to the usual setting where the particle is instantaneously reset to a preferred location (say, the origin), here we consider a finite time resetting process facilitated by an external linear potential $V(x)=\lambda|x|~ (\lambda>0)$. When resetting occurs, the trap is switched on and the particle experiences a force $-\partial_x V(x)$ which helps the particle to return to the resetting location. The trap is switched off as soon as the particle makes a first passage to the origin. Subsequently, the particle resumes its free diffusion motion and the process keeps repeating. In this set-up, the system attains a non-equilibrium steady state. We study the relaxation to this steady state by analytically computing the position distribution of the particle at all time and then analysing this distribution using the spectral properties of the corresponding Fokker-Planck operator. As seen for the instantaneous resetting problem, we observe a `cone spreading' relaxation with travelling fronts such that there is an inner core region around the resetting point that reaches the steady state, while the region outside the core still grows ballistically with time. In addition to the unusual relaxation phenomena, we compute the large deviation functions associated to the corresponding probability density and find that the large deviation functions describe a dynamical transition similar to what is seen previously in case of instantaneous resetting. Notably, our method, based on spectral properties, complements the existing renewal formalism and reveals the intricate mathematical structure responsible for such relaxation phenomena. We verify our analytical results against extensive numerical simulations.


title: Parabolicity conjecture of $F$-isocrystals

abstract: In this article we prove Crew's parabolicity conjecture of $F$-isocrystals. For this purpose, we introduce and study the notion of $\dagger$-hull of a sub-$F$-isocrystal. On the way, we prove a new Lefschetz theorem for overconvergent $F$-isocrystals.


title: Exploring Instance-Level Uncertainty for Medical Detection

abstract: The ability of deep learning to predict with uncertainty is recognized as key for its adoption in clinical routines. Moreover, performance gain has been enabled by modelling uncertainty according to empirical evidence. While previous work has widely discussed the uncertainty estimation in segmentation and classification tasks, its application on bounding-box-based detection has been limited, mainly due to the challenge of bounding box aligning. In this work, we explore to augment a 2.5D detection CNN with two different bounding-box-level (or instance-level) uncertainty estimates, i.e., predictive variance and Monte Carlo (MC) sample variance. Experiments are conducted for lung nodule detection on LUNA16 dataset, a task where significant semantic ambiguities can exist between nodules and non-nodules. Results show that our method improves the evaluating score from 84.57% to 88.86% by utilizing a combination of both types of variances. Moreover, we show the generated uncertainty enables superior operating points compared to using the probability threshold only, and can further boost the performance to 89.52%. Example nodule detections are visualized to further illustrate the advantages of our method.


title: Lifespan Functors and Natural Dualities in Persistent Homology

abstract: We introduce lifespan functors, which are endofunctors on the category of persistence modules that filter out intervals from barcodes according to their boundedness properties. They can be used to classify injective and projective objects in the category of barcodes and the category of pointwise finite-dimensional persistence modules. They also naturally appear in duality results for absolute and relative versions of persistent (co)homology, generalizing previous results in terms of barcodes. Due to their functoriality, we can apply these results to morphisms in persistent homology that are induced by morphisms between filtrations. This lays the groundwork for the efficient computation of barcodes for images, kernels, and cokernels of such morphisms.


title: Decoherence-free mechanism to protect long-range entanglement against decoherence

abstract: Quantum spin models with variable-range interactions can exhibit certain quantum characteristics that a short-ranged model cannot possess. By considering the quantum XYZ model whose interaction strength between different sites varies either exponentially or polynomially, we report the creation of long-range entanglement in dynamics both in the absence and presence of system-bath interactions. Specifically, during closed dynamics, we determine a parameter regime from which the system should start its evolution so that the resulting state after quench can produce a high time-averaged entanglement having low fluctuations. Both in the exponential and power-law decays, it occurs when the magnetic field is weak and the interactions in the z-direction are nonvanishing. When part of the system interacts with the bath repeatedly or is attached to a collection of harmonic oscillators along with dephasing noise in the z-direction, we observe that long-range entanglement of the subparts which are not attached with the environment remains constant with time in the beginning of the evolution, known as freezing of entanglement, thereby demonstrating a method to protect long-range entanglement. We find that the frozen entanglement content in any length and the time up to which freezing occurs called the freezing terminal to follow a complementary relation for all ranges of interactions. However, we find that for a fixed range of entanglement, there exists a critical value of interaction length which leads to the maximum freezing terminal.


title: Evaluating structural edge importance in temporal networks

abstract: To monitor risk in temporal financial networks, we need to understand how individual behaviours affect the global evolution of networks. Here we define a structural importance metric - which we denote as $l_e$ - for the edges of a network. The metric is based on perturbing the adjacency matrix and observing the resultant change in its largest eigenvalues. We then propose a model of network evolution where this metric controls the probabilities of subsequent edge changes. We show using synthetic data how the parameters of the model are related to the capability of predicting whether an edge will change from its value of $l_e$. We then estimate the model parameters associated with five real financial and social networks, and we study their predictability. These methods have application in financial regulation whereby it is important to understand how individual changes to financial networks will impact their global behaviour. It also provides fundamental insights into spectral predictability in networks, and it demonstrates how spectral perturbations can be a useful tool in understanding the interplay between micro and macro features of networks.


title: Vid2Actor: Free-viewpoint Animatable Person Synthesis from Video in the Wild

abstract: Given an "in-the-wild" video of a person, we reconstruct an animatable model of the person in the video. The output model can be rendered in any body pose to any camera view, via the learned controls, without explicit 3D mesh reconstruction. At the core of our method is a volumetric 3D human representation reconstructed with a deep network trained on input video, enabling novel pose/view synthesis. Our method is an advance over GAN-based image-to-image translation since it allows image synthesis for any pose and camera via the internal 3D representation, while at the same time it does not require a pre-rigged model or ground truth meshes for training, as in mesh-based learning. Experiments validate the design choices and yield results on synthetic data and on real videos of diverse people performing unconstrained activities (e.g. dancing or playing tennis). Finally, we demonstrate motion re-targeting and bullet-time rendering with the learned models.


title: Light-matter interactions near photonic Weyl points

abstract: Weyl photons appear when two three-dimensional photonic bands with linear dispersion are degenerate at a single momentum point, labeled as Weyl point. These points have remarkable properties such as being robust topological monopoles of Berry curvature as well as an associated vanishing density of states. In this work, we report on a systematic theoretical study of the quantum optical consequences of such Weyl photons. First, we analyze the dynamics of a single quantum emitter coupled to a Weyl photonic bath as a function of its detuning with respect to the Weyl point and study the corrresponding emission patterns, using both perturbative and exact treatments. Our calculations show an asymmetric dynamical behavior when the emitter is detuned away from the Weyl frequency, as well as different regimes of highly collimated emission, which ultimately translate in a variety of directional collective decays. Besides, we find that the incorporation of staggered mass and hopping terms in the bath Hamiltonian both enriches the observed phenomenology and increases the tunability of the interaction. Finally, we analyze the competition between the coherent and dissipative components of the dynamics for the case of two emitters and derive the conditions under which an effective interacting spin model description is valid.


title: NBIHT: An Efficient Algorithm for 1-bit Compressed Sensing with Optimal Error Decay Rate

abstract: The Binary Iterative Hard Thresholding (BIHT) algorithm is a popular reconstruction method for one-bit compressed sensing due to its simplicity and fast empirical convergence. There have been several works about BIHT but a theoretical understanding of the corresponding approximation error and convergence rate still remains open. This paper shows that the normalized version of BIHT (NBHIT) achieves an approximation error rate optimal up to logarithmic factors. More precisely, using $m$ one-bit measurements of an $s$-sparse vector $x$, we prove that the approximation error of NBIHT is of order $O \left(1 \over m \right)$ up to logarithmic factors, which matches the information-theoretic lower bound $\Omega \left(1 \over m \right)$ proved by Jacques, Laska, Boufounos, and Baraniuk in 2013. To our knowledge, this is the first theoretical analysis of a BIHT-type algorithm that explains the optimal rate of error decay empirically observed in the literature. This also makes NBIHT the first provable computationally-efficient one-bit compressed sensing algorithm that breaks the inverse square root error decay rate $O \left(1 \over m^{1/2} \right)$.


title: SDSS-IV/MaNGA: Can impulsive gaseous inflows explain steep oxygen abundance profiles \& anomalously-low-metallicity regions?

abstract: Gaseous inflows are necessary suppliers of galaxies' star-forming fuel, but are difficult to characterize at the survey scale. We use integral-field spectroscopic measurements of gas-phase metallicity and single-dish radio measurements of total atomic gas mass to estimate the magnitude and frequency of gaseous inflows incident on star-forming galaxies. We reveal a mutual correlation between steep oxygen abundance profiles between $0.25-1.5 R_e$, increased variability of metallicity between $1.25-1.75 R_e$, and elevated HI content at fixed total galaxy stellar mass. Employing a simple but intuitive inflow model, we find that galaxies with total stellar mass less than $10^{10.1} {\rm M_{\odot}}$ have local oxygen abundance profiles consistent with reinvigoration by inflows. Approximately 10-25\% of low-mass galaxies possess signatures of recent accretion, with estimated typical enhancements of approximately 10-90\% in local gas mass surface density. Higher-mass galaxies have limited evidence for such inflows. The large diversity of HI mass implies that inflow-associated gas ought to reside far from the star-forming disk. We therefore propose that a combination of high HI mass, steep metallicity profile between $0.25-1.5 R_e$, and wide metallicity distribution function between $1.25 - 1.75 R_e$ be employed to target possible hosts of inflowing gas for high-resolution radio follow-up.


title: Uniform, localized asymptotics for sub-Riemannian heat kernels and diffusions

abstract: We show how Molchanov's method provides a systematic approach to determining the small-time asymptotics of the heat kernel on a sub-Riemannian manifold away from any abnormal minimizers. The expansion is closely connected to the structure of the minimizing geodesics between two points. If the normal form of the exponential map at the minimal geodesics between two points is sufficiently explicit, a complete asymptotic expansion of the heat kernel can, in principle, be given. (But one can also exhibit metrics for which the exponential map is quite degenerate, so that even giving the leading term of the expansion seems doubtful.) In a different direction, we have uniform bounds on the heat kernel and its derivatives over any compact set with no abnormals, in an inherently local way. The method extends naturally to logarithmic derivatives of the heat kernel, which are closely related to the law of large numbers for the corresponding bridge process. This allows for a general treatment of the small-time behavior of the bridge process on the cut locus, as well as the determination of the limiting measure of the bridge process if the normal form of the exponential map at the minimal geodesics between two points is sufficiently explicit. Further, we give an expression for the leading term of the logarithmic derivative of the heat kernel, of any order, as a cumulant of geometrically natural random variables and a characterization of the cut locus in terms of the blow-up of the logarithmic second derivative of the heat kernel. This method also provides uniform, local estimates of the logarithmic derivatives.


title: Stahl-Totik Regularity for Dirac Operators

abstract: We develop a theory of regularity for Dirac operators with uniformly locally square-integrable operator data. This is motivated by Stahl--Totik regularity for orthogonal polynomials and by recent developments for continuum Schr\"odinger operators, but contains significant new phenomena. We prove that the symmetric Martin function at $\infty$ for the complement of the essential spectrum has the two-term asymptotic expansion $\Im \left( z - \frac{b}{2 z}\right) + o(\frac 1z)$ as $z \to i \infty$, which is seen as a thickness statement for the essential spectrum. The constant $b$ plays the role of a renormalized Robin constant and enters a universal inequality involving the lower average $L^2$-norm of the operator data. However, we show that regularity of Dirac operators is not precisely characterized by a single scalar equality involving $b$ and is instead characterized by a family of equalities. This work also contains a sharp Combes--Thomas estimate (root asymptotics of eigensolutions), a study of zero counting measures, and applications to ergodic and decaying operator data.