title: SNIa-Cosmology Analysis Results from Simulated LSST Images: from Difference Imaging to Constraints on Dark Energy

abstract: The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ${\sim}10^6$ transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration (DESC) Data Challenge (DC2) simulation for the Wide-Fast-Deep (WFD) survey area. DC2 is the first large-scale (300 deg$^2$) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ${\sim}$15 deg$^2$ of DC2 over a 5-year time-span in which artificial point-sources from Type Ia Supernovae (SNIa) light curves have been overlaid onto the images. We measure the detection efficiency as a function of Signal-to-Noise Ratio (SNR) and find a $50\%$ efficiency at $\rm{SNR}=5.8$. The magnitude limits for each filter are: $u=23.66$, $g=24.69$, $r=24.06$, $i=23.45$, $z=22.54$, $y=21.62$ $\rm{mag}$. The artifact contamination is $\sim90\%$ of detections, corresponding to $\sim1000$ artifacts/deg$^2$ in $g$ band, and falling to 300 per deg$^2$ in $y$ band. The photometry has biases $<1\%$ for magnitudes $19.5 < m <23$. Our DIA performance on simulated images is similar to that of the Dark Energy Survey pipeline applied to real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of SNR, redshift, and fitted light-curve properties. Applying a realistic SNIa-cosmology analysis for redshifts $z<1$, we recover the input cosmology parameters to within statistical uncertainties.


title: Higher-Order Coverage Errors of Batching Methods via Edgeworth Expansions on $t$-Statistics

abstract: While batching methods have been widely used in simulation and statistics, it is open regarding their higher-order coverage behaviors and whether one variant is better than the others in this regard. We develop techniques to obtain higher-order coverage errors for batching methods by building Edgeworth-type expansions on $t$-statistics. The coefficients in these expansions are intricate analytically, but we provide algorithms to estimate the coefficients of the $n^{-1}$ error term via Monte Carlo simulation. We provide insights on the effect of the number of batches on the coverage error where we demonstrate generally non-monotonic relations. We also compare different batching methods both theoretically and numerically, and argue that none of the methods is uniformly better than the others in terms of coverage. However, when the number of batches is large, sectioned jackknife has the best coverage among all.


title: Hypersurfaces, Geodesics and Isoperimetric Inequalities in Cartan-Hadamard Manifolds

abstract: We prove an inequality for submanifolds of Cartan-Hadamard manifolds, which relates the geometry of a submanifold to the measure of the geodesics in the ambient space which it intersects. For hypersurfaces, this gives an extension of Banchoff and Pohl's isoperimetric inequality to spaces of non-positive curvature. We also prove a modified version of Croke's isoperimetric inequality for hypersurfaces immersed in Cartan-Hadamard manifolds and a sharp, quantitative version of an isoperimetric inequality of Yau in spaces of negative curvature. We discuss the relationship between these results, and we develop several facts about the spaces of geodesics in Cartan-Hadamards manifold that may be of independent interest.


title: Celestial IR divergences in general most-subleading-color gluon and gravity amplitudes

abstract: Gluon amplitudes at most-subleading order in the $1/N$ expansion share a remarkable simplicity with graviton amplitudes: collinear divergences are completely absent in both and, as a consequence, their full IR behavior arises from soft gluon/graviton exchange among the external states. In this paper we study the effect of all-loop IR divergences of celestial most-subleading color gluon amplitudes and their similarities with the celestial gravity case. In particular, a simple celestial exponentiation formula for the dipole part can be written. We also analize how this exponentiation is modified by non-dipole contributions. Finally we also show that, in the Regge limit, the soft factor satisfies the Knizhnik-Zamolodchikov equation hinting at the possibility that, in this limit, an effective Wess-Zumino-Witten model would describe the dynamics of the infrared sector.


title: Silicon Photonic Architecture for Training Deep Neural Networks with Direct Feedback Alignment

abstract: There has been growing interest in using photonic processors for performing neural network inference operations; however, these networks are currently trained using standard digital electronics. Here, we propose on-chip training of neural networks enabled by a CMOS-compatible silicon photonic architecture to harness the potential for massively parallel, efficient, and fast data operations. Our scheme employs the direct feedback alignment training algorithm, which trains neural networks using error feedback rather than error backpropagation, and can operate at speeds of trillions of multiply-accumulate (MAC) operations per second while consuming less than one picojoule per MAC operation. The photonic architecture exploits parallelized matrix-vector multiplications using arrays of microring resonators for processing multi-channel analog signals along single waveguide buses to calculate the gradient vector for each neural network layer in situ. We also experimentally demonstrate training deep neural networks with the MNIST dataset using on-chip MAC operation results. Our novel approach for efficient, ultra-fast neural network training showcases photonics as a promising platform for executing AI applications.


title: Hierarchical Clustering: New Bounds and Objective

abstract: Hierarchical Clustering has been studied and used extensively as a method for analysis of data. More recently, Dasgupta [2016] defined a precise objective function. Given a set of $n$ data points with a weight function $w_{i,j}$ for each two items $i$ and $j$ denoting their similarity/dis-similarity, the goal is to build a recursive (tree like) partitioning of the data points (items) into successively smaller clusters. He defined a cost function for a tree $T$ to be $Cost(T) = \sum_{i,j \in [n]} \big(w_{i,j} \times |T_{i,j}| \big)$ where $T_{i,j}$ is the subtree rooted at the least common ancestor of $i$ and $j$ and presented the first approximation algorithm for such clustering. Then Moseley and Wang [2017] considered the dual of Dasgupta's objective function for similarity-based weights and showed that both random partitioning and average linkage have approximation ratio $1/3$ which has been improved in a series of works to $0.585$ [Alon et al. 2020]. Later Cohen-Addad et al. [2019] considered the same objective function as Dasgupta's but for dissimilarity-based metrics, called $Rev(T)$. It is shown that both random partitioning and average linkage have ratio $2/3$ which has been only slightly improved to $0.667078$ [Charikar et al. SODA2020]. Our first main result is to consider $Rev(T)$ and present a more delicate algorithm and careful analysis that achieves approximation $0.71604$. We also introduce a new objective function for dissimilarity-based clustering. For any tree $T$, let $H_{i,j}$ be the number of $i$ and $j$'s common ancestors. Intuitively, items that are similar are expected to remain within the same cluster as deep as possible. So, for dissimilarity-based metrics, we suggest the cost of each tree $T$, which we want to minimize, to be $Cost_H(T) = \sum_{i,j \in [n]} \big(w_{i,j} \times H_{i,j} \big)$. We present a $1.3977$-approximation for this objective.


title: Perfectly absorbed and emitted currents by complex potentials in nonlinear media

abstract: Recently it was demonstrated that the concept of a spectral singularity (SS) can be generalized to waves propagating in nonlinear media, like matter waves or electromagnetic waves in Kerr media. The corresponding solutions represent nonlinear currents sustained by a localized linear complex potential in a nonlinear Schr\"odinger equation. A key feature allowing the nonlinear generalization of a SS is a possibility to reduce a nonlinear current to the linear limit, where a SS has the unambiguous definition. In the meantime, known examples of nonlinear modes bifurcating from linear spectral singularities are few and belong to the specific class of constant-amplitude waves. Here we propose to extend the class of nonlinear SSs by incorporating solutions whose amplitudes are inhomogeneous. We show that the continuation from the linear limit requires a deformation of the complex potential, and this deformation is not unique. Examples include the deformation preserving the gain-and-loss distribution and the deformation preserving geometry of the potential. For the case example of a rectangular potential, we demonstrate that the nonlinear currents can be divided into two types: solutions of the first type bifurcate from the linear spectral singularities, and solutions of the second type cannot be reduced to the linear limit.


title: Generalized active information: extensions to unbounded domains

abstract: In the last three decades, several measures of complexity have been proposed. Up to this point, most of such measures have only been developed for finite spaces. In these scenarios the baseline distribution is uniform. This makes sense because, among other things, the uniform distribution is the measure of maximum entropy over the relevant space. Active information traditionally assumes a finite interval universe of discourse but can be extended to other cases where maximum entropy is defined. Illustrating this is the purpose of this paper. Disequilibrium from maximum entropy, measured as active information, can be evaluated from baselines with unbounded support.

categories: cs.IT math.IT math.ST physics.data-an stat.TH


title: Elastic three-dimensional metaframe for selective wave filtering and polarization control

abstract: We experimentally achieve selective wave filtering and polarization control in a three-dimensional elastic frame embedding local resonators. By connecting multi-resonating elements to a frame structure, a complete low-frequency, subwavelength bandgap with strong selective filtering properties is obtained. Theory and experiments demonstrate the metaframe capability to selectively stop transverse waves while allowing longitudinal wave propagation, as in 'fluid-like' elasticity. This peculiar behaviour, together with the complete bandgap structure, may open opportunities in the context of wave control, envisaging concurrent applications for three-dimensional filters and elastic wave polarizers.


title: Flatee: Federated Learning Across Trusted Execution Environments

abstract: Federated learning allows us to distributively train a machine learning model where multiple parties share local model parameters without sharing private data. However, parameter exchange may still leak information. Several approaches have been proposed to overcome this, based on multi-party computation, fully homomorphic encryption, etc.; many of these protocols are slow and impractical for real-world use as they involve a large number of cryptographic operations. In this paper, we propose the use of Trusted Execution Environments (TEE), which provide a platform for isolated execution of code and handling of data, for this purpose. We describe Flatee, an efficient privacy-preserving federated learning framework across TEEs, which considerably reduces training and communication time. Our framework can handle malicious parties (we do not natively solve adversarial data poisoning, though we describe a preliminary approach to handle this).


title: HybridQ: A Hybrid Simulator for Quantum Circuits

abstract: Developing state-of-the-art classical simulators of quantum circuits is of utmost importance to test and evaluate early quantum technology and understand the true potential of full-blown error-corrected quantum computers. In the past few years, multiple theoretical and numerical advances have continuously pushed the boundary of what is classically simulable, hence the development of a plethora of tools which are often limited to a specific purpose or designed for a particular hardware (e.g. CPUs vs. GPUs). Moreover, such tools are typically developed using tailored languages and syntax, which makes it hard to compare results from, and create hybrid approaches using, different simulation techniques. To support unified and optimized use of these techniques across platforms, we developed HybridQ, a highly extensible platform designed to provide a common framework to integrate multiple state-of-the-art techniques to run on a variety of hardware. The philosophy behind its development has been driven by three main pillars: "Easy to Use", "Easy to Extend", and "Use the Best Available Technology". The powerful tools of HybridQ allow users to manipulate, develop, and extend noiseless and noisy circuits for different hardware architectures. HybridQ supports large-scale high-performance computing (HPC) simulations, automatically balancing workload among different processor nodes and enabling the use of multiple backends to maximize parallel efficiency. Everything is then glued together by a simple and expressive language that allows seamless switching from one technique to another as well as from one hardware to the next, without the need to write lengthy translations, thus greatly simplifying the development of new hybrid algorithms and techniques.


title: A null test of the Cosmological Principle with BAO measurements

abstract: The assumption of a Friedmann-Lema\^itre-Roberton-Walker (FLRW) Universe is one of the fundamental pillars of modern Cosmology. It is crucial to confront it against cosmological observations in order to confirm (or rule out) the validity of the $\Lambda$CDM paradigm. We perform a consistency test of the FLRW metric using measurements of both radial and transversal modes of Baryonic Acoustic Oscillations in a model-independent fashion using Gaussian Processes. The impact of different prior assumptions for the sound horizon scale and the Hubble Constant is also explored. We find a mild deviation of the FLRW hypothesis ($\gtrsim 3\sigma$ confidence level) at lower redshift ranges ($0.1<z<0.3$) under the assumption the SH0ES $H_0$ prior, which does not occur for the $H_0$ prior by Planck. We also find that different reconstruction kernels reduce the statistical significance of such a deviation. At higher reshift ranges, the null test are able to confirm the FLRW assumption within $3\sigma$ confidence level regardless of the priors.


title: Radiation hydrodynamics modeling of kilonovae with SNEC

abstract: We develop a method to compute synthetic kilonova light curves that combines numerical relativity simulations of neutron star mergers and the $\texttt{SNEC}$ radiation-hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the $\texttt{SNEC}$ results with those of two semi-analytic light curve models. We apply our code to the calculation of color light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models which neglect non-LTE effects and possible additional energy injection from the merger remnant, and to the need to go beyond the assumption of spherical symmetry adopted in this work.


title: Sampling from multimodal distributions using tempered Hamiltonian transitions

abstract: Hamiltonian Monte Carlo (HMC) methods are widely used to draw samples from unnormalized target densities due to high efficiency and favorable scalability with respect to increasing space dimensions. However, HMC struggles when the target distribution is multimodal, because the maximum increase in the potential energy function (i.e., the negative log density function) along the simulated path is bounded by the initial kinetic energy, which follows a half of the $\chi_d^2$ distribution, where d is the space dimension. In this paper, we develop a Hamiltonian Monte Carlo method where the constructed paths can travel across high potential energy barriers. This method does not require the modes of the target distribution to be known in advance. Our approach enables frequent jumps between the isolated modes of the target density by continuously varying the mass of the simulated particle while the Hamiltonian path is constructed. Thus, this method can be considered as a combination of HMC and the tempered transitions method. Compared to other tempering methods, our method has a distinctive advantage in the Gibbs sampler settings, where the target distribution changes at each step. We develop a practical tuning strategy for our method and demonstrate that it can construct globally mixing Markov chains targeting high-dimensional, multimodal distributions, using mixtures of normals and a sensor network localization problem.


title: Constraints on leptoquarks from lepton-flavour-violating tau-lepton processes

abstract: Leptoquarks are ubiquitous in several extensions of the Standard Model and seem to be able to accommodate the universality-violation-driven $B$-meson-decay anomalies and the $(g-2)_\mu$ discrepancy interpreted as deviations from the Standard Model predictions. In addition, the search for lepton-flavour violation in the charged sector is, at present, a major research program that could also be facilitated by the dynamics generated by leptoquarks. In this article, we consider a rather wide framework of both scalar and vector leptoquarks as the generators of lepton-flavour violation in processes involving the tau lepton. We single out its couplings to leptoquarks, thus breaking universality in the lepton sector, and we integrate out leptoquarks at tree level, generating the corresponding dimension-6 operators of the Standard Model Effective Field Theory. In Ref. [1] we obtained model-independent bounds on the Wilson coefficients of those operators contributing to lepton-flavour-violating hadron tau decays and $\ell$-$\tau$ conversion in nuclei, with $\ell=e,\mu$. Hence, we use those results to translate the bounds into the couplings of leptoquarks to the Standard Model fermions.


title: Elliptic hypergeometric function and $6j$-symbols for the SL(2,$\mathbb{C}$) group

abstract: We show that the complex hypergeometric function describing $6j$-symbols for $SL(2,\mathbb{C})$ group is a special degeneration of the $V$-function -- an elliptic analogue of the Euler-Gauss $_2F_1$ hypergeometric function. For this function, we derive mixed difference-recurrence relations as limiting forms of the elliptic hypergeometric equation and some symmetry transformations. At the intermediate steps of computations, there emerge a function describing the $6j$-symbols for the Faddeev modular double and the corresponding difference equations and symmetry transformations.


title: Uniform Fourier restriction for convex curves

abstract: We extend the estimates for maximal Fourier restriction operators proved by M\"{u}ller, Ricci, and Wright in \cite{MR3960255} and Ramos in \cite{MR4055940} to the case of arbitrary convex curves in the plane, with constants uniform in the curve. The improvement over M\"{u}ller, Ricci, and Wright and Ramos is given by the removal of the $\mathcal{C}^2$ regularity condition on the curve. This requires the choice of an appropriate measure for each curve, that is suggested by an affine invariant construction of Oberlin in \cite{MR1960918}. As corollaries, we obtain a uniform Fourier restriction theorem for arbitrary convex curves, and a result on the Lebesgue points of the Fourier transform on the curve.


title: Cooperative multi-agent reinforcement learning for high-dimensional nonequilibrium control

abstract: Experimental advances enabling high-resolution external control create new opportunities to produce materials with exotic properties. In this work, we investigate how a multi-agent reinforcement learning approach can be used to design external control protocols for self-assembly. We find that a fully decentralized approach performs remarkably well even with a "coarse" level of external control. More importantly, we see that a partially decentralized approach, where we include information about the local environment allows us to better control our system towards some target distribution. We explain this by analyzing our approach as a partially-observed Markov decision process. With a partially decentralized approach, the agent is able to act more presciently, both by preventing the formation of undesirable structures and by better stabilizing target structures as compared to a fully decentralized approach.


title: Towards the Finite-$N$ Twisted Holography from the Geometry of Phase Space

abstract: We study ring of functions on the (classical and quantized) phase space of $2d$ BF theory with the gauge group $\mathrm{GL}_N$ coupled to a $1d$ quantum mechanics with global symmetry $\mathrm{GL}_K$. These functions are gauge-invariant local observables of the coupled system, which appears as boundary side of an early example of twisted holography involving stacks of $N$ D2-branes and $K$ D4-branes intersecting along a real line. We first construct the classical phase space of this system and describe its ring of functions and their large-$N$ limit. We next compute the Hilbert series of these algebras for finite-$N$ and also in the large-$N$ limit. We then study the quantization of this phase space and the deformation quantization of its ring of functions, elaborate its relation to the Yangian, and construct its co-product. Finally, we identify these quantized algebras with the quantized Coulomb-branch algebras of certain $3d$ $\mathcal{N}=4$ quiver gauge theories.


title: Continuity, Uniqueness and Long-Term Behavior of Nash Flows Over Time

abstract: We consider a dynamic model of traffic that has received a lot of attention in the past few years. Users control infinitesimal flow particles aiming to travel from an origin to a destination as quickly as possible. Flow patterns vary over time, and congestion effects are modeled via queues, which form whenever the inflow into a link exceeds its capacity. Despite lots of interest, some very basic questions remain open in this model. We resolve a number of them: - We show uniqueness of journey times in equilibria. - We show continuity of equilibria: small perturbations to the instance or to the traffic situation at some moment cannot lead to wildly different equilibrium evolutions. - We demonstrate that, assuming constant inflow into the network at the source, equilibria always settle down into a "steady state" in which the behavior extends forever in a linear fashion. One of our main conceptual contributions is to show that the answer to the first two questions, on uniqueness and continuity, are intimately connected to the third. To resolve the third question, we substantially extend the approach of Cominetti et al., who show a steady-state result in the regime where the input flow rate is smaller than the network capacity.


title: FIXP-membership via Convex Optimization: Games, Cakes, and Markets

abstract: We introduce a new technique for proving membership of problems in FIXP - the class capturing the complexity of computing a fixed-point of an algebraic circuit. Our technique constructs a "pseudogate" which can be used as a black box when building FIXP circuits. This pseudogate, which we term the "OPT-gate", can solve most convex optimization problems. Using the OPT-gate, we prove new FIXP-membership results, and we generalize and simplify several known results from the literature on fair division, game theory and competitive markets. In particular, we prove complexity results for two classic problems: computing a market equilibrium in the Arrow-Debreu model with general concave utilities is in FIXP, and computing an envy-free division of a cake with very general valuations is FIXP-complete. We further showcase the wide applicability of our technique, by using it to obtain simplified proofs and extensions of known FIXP-membership results for equilibrium computation for various types of strategic games, as well as the pseudomarket mechanism of Hylland and Zeckhauser.


title: On the Search for Multicenter AdS Black Holes from M-theory

abstract: We study the effective potentials for various probe branes surrounding AdS$_4$ black holes with massive halos in consistent truncations of M-theory on the Sasaki-Einstein$_7$ manifolds $Q^{111}$ and $M^{111}$. These probes are either M2 branes extended in spacetime or "particle-like" probes such as internally wrapped M2 branes and, upon reduction to type IIA string theory, D6 branes corresponding to baryon operators in the dual Chern-Simons theory. We find both global and local minima of the potential outside the horizon, indicating the existence of stable and metastable multicenter AdS black holes in the extreme mass ratio regime, at fixed temperature and charges. For the planar case, we also find an instability towards nucleation of spacetime-filling M2 branes. With this analysis, we address some open questions on the holographic description of glassy phases of matter.


title: Robust Eigenvectors of Symmetric Tensors

abstract: The tensor power method generalizes the matrix power method to higher order arrays, or tensors. Like in the matrix case, the fixed points of the tensor power method are the eigenvectors of the tensor. While every real symmetric matrix has an eigendecomposition, the vectors generating a symmetric decomposition of a real symmetric tensor are not always eigenvectors of the tensor. In this paper we show that whenever an eigenvector is a generator of the symmetric decomposition of a symmetric tensor, then (if the order of the tensor is sufficiently high) this eigenvector is robust, i.e., it is an attracting fixed point of the tensor power method. We exhibit new classes of symmetric tensors whose symmetric decomposition consists of eigenvectors. Generalizing orthogonally decomposable tensors, we consider equiangular tight frame decomposable and equiangular set decomposable tensors. Our main result implies that such tensors can be decomposed using the tensor power method.


title: Multimodal Virtual Point 3D Detection

abstract: Lidar-based sensing drives current autonomous vehicles. Despite rapid progress, current Lidar sensors still lag two decades behind traditional color cameras in terms of resolution and cost. For autonomous driving, this means that large objects close to the sensors are easily visible, but far-away or small objects comprise only one measurement or two. This is an issue, especially when these objects turn out to be driving hazards. On the other hand, these same objects are clearly visible in onboard RGB sensors. In this work, we present an approach to seamlessly fuse RGB sensors into Lidar-based 3D recognition. Our approach takes a set of 2D detections to generate dense 3D virtual points to augment an otherwise sparse 3D point cloud. These virtual points naturally integrate into any standard Lidar-based 3D detectors along with regular Lidar measurements. The resulting multi-modal detector is simple and effective. Experimental results on the large-scale nuScenes dataset show that our framework improves a strong CenterPoint baseline by a significant 6.6 mAP, and outperforms competing fusion approaches. Code and more visualizations are available at https://tianweiy.github.io/mvp/


title: Multi-band light curves from eccentric accreting supermassive black hole binaries

abstract: We use long-run, high-resolution hydrodynamics simulations to compute the multi-wavelength light curves (LCs) from thermal disk emission around accreting equal-mass supermassive black hole (BH) binaries, with a focus on revealing binary eccentricity. LCs are obtained by modeling the disk thermodynamics with an adiabatic equation of state, a local blackbody cooling prescription, and corrections to approximate the effects of radiation pressure. We find that modulation of multi-band LCs on the orbital time scale are generally in-phase (to within $\sim\,$2% of a binary orbital period), but they contain pulse substructure in the time domain that is not necessarily reflected in BH accretion rates $\dot M$. We thus predict that binary-hosting AGN will exhibit highly correlated, in-phase, periodic brightness modulations in their low-energy disk emission. However, detectability of these modulations in multi-wavelength observing campaigns could be seriously compromised because observed stochastic variability in AGNs typically has a higher amplitude than our proposed signal. It is possible that observations over temporal baselines of many binary periods may make the signal more prominent, but this would need to be analyzed carefully. If jet emission is predicted by $\dot{M}$, then we predict a weaker correlation with low-energy disk emission due to the differing sub-peak structure. For the binary parameters we explore, we show that LC variability due to hydrodynamics likely dominates Doppler brightening for all equal-mass binaries with disk Mach numbers $\lesssim 20$. A promising signature of eccentricity is weak or absent "lump" periodicity. We find hints that a significant lag exists between $\dot{M}$ and low-energy disk emission for circular binaries, but they are in-phase for eccentric binaries, which might explain some "orphan" blazar flares with no $\gamma$-ray counterpart.


title: Search for oscillations of fundamental constants using molecular spectroscopy

abstract: A possible implication of an ultralight dark matter (UDM) field interacting wibeginth the Standard Model (SM) degrees of freedom is oscillations of fundamental constants. Here, we establish direct experimental bounds on the coupling of an oscillating UDM field to the up, down, and strange quarks and to the gluons, for oscillation frequencies between 10 Hz and 10^8 Hz. We employ spectroscopic experiments that take advantage of the dependence of molecular transition frequencies on the nuclear masses. Our results apply to previously unexplored frequency bands, and improve on existing bounds at frequencies > 5 MHz. We identify a sector of UDM - SM coupling space where the bounds from Equivalence Principle tests may be challenged by next-generation experiments of the present kind.


title: Signatures of Primordial Gravitational Waves on the Large-Scale Structure of the Universe

abstract: We study the generation and evolution of second-order energy-density perturbations arising from primordial gravitational waves. Such "tensor-induced scalar modes" approximately evolve as standard linear matter perturbations and may leave observable signatures in the Large-Scale Structure of the Universe. We study the imprint on the matter power-spectrum of some primordial models which predict a large gravitational-wave signal at high frequencies. This novel mechanism in principle allows us to constrain or detect primordial gravitational waves by looking at specific features in the matter or galaxy power-spectrum, thereby allowing to probe them on a range of scales unexplored so far.


title: Guided Sampling-based Evolutionary Deep Neural Network for Intelligent Fault Diagnosis

abstract: The diagnostic performance of most of the deep learning models is greatly affected by the selection of model architecture and hyperparameters. Manual selection of model architecture is not feasible as training and evaluating the different architectures of deep learning models is a time-consuming process. Therefore, we have proposed a novel framework of evolutionary deep neural network which uses policy gradient to guide the evolution of DNN architecture towards maximum diagnostic accuracy. We have formulated a policy gradient-based controller which generates an action to sample the new model architecture at every generation such that the optimality is obtained quickly. The fitness of the best model obtained is used as a reward to update the policy parameters. Also, the best model obtained is transferred to the next generation for quick model evaluation in the NSGA-II evolutionary framework. Thus, the algorithm gets the benefits of fast non-dominated sorting as well as quick model evaluation. The effectiveness of the proposed framework has been validated on three datasets: the Air Compressor dataset, Case Western Reserve University dataset, and Paderborn university dataset.


title: Performance vs Persistence : Assess the alpha to identify outperformers

abstract: The number of pension funds has multiplied exponentially over the last decade. Active portfolio management requires a precise analysis of the performance drivers. Several risk and performance attribution metrics have been developed since the 70s to guide investors in their investment choices. Based on the study made by Fama and French (2010) we reproduce the experiment they had carried out in order to complete their work using additionnal features. Throughout this study we draw a parallel between the results obtained by Fama and French (2010) with the 3-factor model. The aim of this paper is to assess the usefulness of two additional factors in the analysis of the persistence of alphas. We also look at the quality of the manager through his investment choices in order to generate alpha considering the environment in which he operates.


title: How to create an artificial magnetosphere for Mars

abstract: If humanity is ever to consider substantial, long-term colonization of Mars, the resources needed are going to be extensive. For a long-term human presence on Mars to be established, serious thought would need to be given to terraforming the planet. One major requirement for such terraforming is having the protection of a planetary magnetic field which Mars currently does not have. In this article we explore comprehensively for the first time, the practical and engineering challenges that affect the feasibility of creating an artificial magnetic field capable of encompassing Mars. This includes the concerns that define the design, where to locate the magnetic field generator and possible construction strategies. The rationale here is not to justify the need for a planetary magnetosphere but to put figures on the practicalities so as to be able to weigh the pros and cons of the different engineering approaches. The optimum solution proposed is completely novel, although inspired by natural situations and fusion plasma techniques. The solution with the lowest power, assembly and mass is to create an artificial charged particle ring (similar in form to a "radiation belt"), around the planet possibly formed by ejecting matter from one of the moons of Mars (in fashion similar to that that forms the Io-Jupiter plasma torus), but using electromagnetic and plasma waves to drive a net current in the ring(s) that results in an overall magnetic field. With a new era of space exploration underway, this is the time to start thinking about these new and bold future concepts and to begin filling strategic knowledge gaps. Furthermore, the principles explored here are also applicable to smaller scale objects like manned spacecraft, space stations or moon bases, which would benefit from the creation of protective mini-magnetospheres.

categories: physics.space-ph astro-ph.EP astro-ph.IM physics.pop-ph


title: Learning Generalized Gumbel-max Causal Mechanisms

abstract: To perform counterfactual reasoning in Structural Causal Models (SCMs), one needs to know the causal mechanisms, which provide factorizations of conditional distributions into noise sources and deterministic functions mapping realizations of noise to samples. Unfortunately, the causal mechanism is not uniquely identified by data that can be gathered by observing and interacting with the world, so there remains the question of how to choose causal mechanisms. In recent work, Oberst & Sontag (2019) propose Gumbel-max SCMs, which use Gumbel-max reparameterizations as the causal mechanism due to an intuitively appealing counterfactual stability property. In this work, we instead argue for choosing a causal mechanism that is best under a quantitative criteria such as minimizing variance when estimating counterfactual treatment effects. We propose a parameterized family of causal mechanisms that generalize Gumbel-max. We show that they can be trained to minimize counterfactual effect variance and other losses on a distribution of queries of interest, yielding lower variance estimates of counterfactual treatment effect than fixed alternatives, also generalizing to queries not seen at training time.


title: DriverGym: Democratising Reinforcement Learning for Autonomous Driving

abstract: Despite promising progress in reinforcement learning (RL), developing algorithms for autonomous driving (AD) remains challenging: one of the critical issues being the absence of an open-source platform capable of training and effectively validating the RL policies on real-world data. We propose DriverGym, an open-source OpenAI Gym-compatible environment specifically tailored for developing RL algorithms for autonomous driving. DriverGym provides access to more than 1000 hours of expert logged data and also supports reactive and data-driven agent behavior. The performance of an RL policy can be easily validated on real-world data using our extensive and flexible closed-loop evaluation protocol. In this work, we also provide behavior cloning baselines using supervised learning and RL, trained in DriverGym. We make DriverGym code, as well as all the baselines publicly available to further stimulate development from the community.


title: Impact of loss functions on the performance of a deep neural network designed to restore low-dose digital mammography

abstract: Digital mammography is still the most common imaging tool for breast cancer screening. Although the benefits of using digital mammography for cancer screening outweigh the risks associated with the x-ray exposure, the radiation dose must be kept as low as possible while maintaining the diagnostic utility of the generated images, thus minimizing patient risks. Many studies investigated the feasibility of dose reduction by restoring low-dose images using deep neural networks. In these cases, choosing the appropriate training database and loss function is crucial and impacts the quality of the results. In this work, a modification of the ResNet architecture, with hierarchical skip connections, is proposed to restore low-dose digital mammography. We compared the restored images to the standard full-dose images. Moreover, we evaluated the performance of several loss functions for this task. For training purposes, we extracted 256,000 image patches from a dataset of 400 images of retrospective clinical mammography exams, where different dose levels were simulated to generate low and standard-dose pairs. To validate the network in a real scenario, a physical anthropomorphic breast phantom was used to acquire real low-dose and standard full-dose images in a commercially avaliable mammography system, which were then processed through our trained model. An analytical restoration model for low-dose digital mammography, previously presented, was used as a benchmark in this work. Objective assessment was performed through the signal-to-noise ratio (SNR) and mean normalized squared error (MNSE), decomposed into residual noise and bias. Results showed that the perceptual loss function (PL4) is able to achieve virtually the same noise levels of a full-dose acquisition, while resulting in smaller signal bias compared to other loss functions.


title: UV Spectropolarimetry with Polstar: Protoplanetary Disks

abstract: Polstar is a proposed NASA MIDEX mission that would feature a high resolution UV spectropolarimeter capable of measure all four Stokes parameters onboard a 60cm telescope. The mission would pioneer the field of time-domain UV spectropolarimetry. Time domain UV spectropolarimetry offers the best resource to determine the geometry and physical conditions of protoplanetary disks from the stellar surface to <5 AU. We detail two key objectives that a dedicated time domain UV spectropolarimetry survey, such as that enabled by Polstar, could achieve: 1) Test the hypothesis that magneto-accretion operating in young planet-forming disks around lower-mass stars transitions to boundary layer accretion in planet-forming disks around higher mass stars; and 2) Discriminate whether transient events in the innermost regions of planet-forming disks of intermediate mass stars are caused by inner disk mis-alignments or from stellar or disk emissions.


title: Analysis of chronic diseases progression using stochastic modeling

abstract: This book handles the fatty liver disease from the bio-statistical point of view . It discusses the disease process in the simple general form of health-disease-death multi-states model . Continuous Time Markov Chains are used to estimate the rate transition matrix utilizing the MLE and Quasi-Newton formula , once obtained , the probability transition matrix can be estimated by exponentiation of the rate matrix . The probability transition matrix can also be obtained by solving the forward Kolmogorov differential equations , which yields more stable solution than exponentiation of rate matrix. The disease process is expanded in 9 states model to explain the transition among the detailed stages of the disease process , in more elaborate form. The probability transition matrix is used to estimate the number of patients in each stage , this matrix along with the rate transition matrix , both are used to estimate life expectancy of patients is each stage. These statistical indices are of great value as they can help the health policy makers and medical insurance managers to allocate the resources for investigating and treating patients in different stages of the disease . This method is of a high potential value to be used in longitudinal studies conducted by the pharmaceutical companies to evaluate the effect of anti-fibrotic drugs used to treat patients within the initial stages of fibrosis . Poisson regression model is also used to relate the high risk covariates such as type 2 diabetes, hypercholesterolemia , obesity and hypertension with the rate of progression and evolution of the stages of the disease over time. The general model , the expanded model and the model with covariates are illustrated by artificial hypothetical examples to demonstrate the mathematical statistical indices .


title: High-accuracy numerical models of Brownian thermal noise in thin mirror coatings

abstract: Brownian coating thermal noise in detector test masses is limiting the sensitivity of current gravitational-wave detectors on Earth. Therefore, accurate numerical models can inform the ongoing effort to minimize Brownian coating thermal noise in current and future gravitational-wave detectors. Such numerical models typically require significant computational resources and time, and often involve closed-source commercial codes. In contrast, open-source codes give complete visibility and control of the simulated physics, enable direct assessment of the numerical accuracy, and support the reproducibility of results. In this article, we use the open-source SpECTRE numerical relativity code and adopt a novel discontinuous Galerkin numerical method to model Brownian coating thermal noise. We demonstrate that SpECTRE achieves significantly higher accuracy than a previous approach at a fraction of the computational cost. Furthermore, we numerically model Brownian coating thermal noise in multiple sub-wavelength crystalline coating layers for the first time. Our new numerical method has the potential to enable fast exploration of realistic mirror configurations, and hence to guide the search for optimal mirror geometries, beam shapes and coating materials for gravitational-wave detectors.

categories: astro-ph.IM gr-qc physics.comp-ph physics.ins-det


title: Convolutional Nets Versus Vision Transformers for Diabetic Foot Ulcer Classification

abstract: This paper compares well-established Convolutional Neural Networks (CNNs) to recently introduced Vision Transformers for the task of Diabetic Foot Ulcer Classification, in the context of the DFUC 2021 Grand-Challenge, in which this work attained the first position. Comprehensive experiments demonstrate that modern CNNs are still capable of outperforming Transformers in a low-data regime, likely owing to their ability for better exploiting spatial correlations. In addition, we empirically demonstrate that the recent Sharpness-Aware Minimization (SAM) optimization algorithm considerably improves the generalization capability of both kinds of models. Our results demonstrate that for this task, the combination of CNNs and the SAM optimization process results in superior performance than any other of the considered approaches.


title: Simulation of energy transport by dark matter scattering in stars

abstract: Asymmetric dark matter (ADM) that is captured in stars can act as an efficient conductor of heat. Small ADM-induced changes in a star's temperature gradient are known to alter neutrino fluxes and asteroseismological signatures, erase convective cores and modify a star's main sequence lifetime. The Sun's proximity to us makes it an ideal laboratory for studying these effects. However, the two formalisms commonly used to parametrize such heat transport were developed over 30 years ago, and calibrated with a single set of simulations. What's more, both are based on assumptions that break down at the Knudsen transition, where heat transport is maximized. We construct a Monte Carlo simulation to exactly solve the Boltzmann collision equation, determining the steady-state distribution and luminosity carried in stars by ADM with cross sections that depend on velocity and momentum. We find that, although the established (Gould and Raffelt) formalism based on local thermal equilibrium does well for constant cross sections, the isothermal (Spergel and Press) method actually performs better across all models with a simple, universal rescaling function. Based on simulation results, we provide recommendations on the parametrization of DM heat transport in stellar evolution models.


title: Observational limits on the early-time dust mass in SN1987A

abstract: In recent years, dust masses of a few tenths of a solar mass have been found in the expanding ejecta of a number of core-collapse supernovae. How dust forms in such quantities remains poorly understood; theories of dust formation predict lower total masses and much faster formation rates than observations imply. One suggestion to reconcile observations and theory was made by Dwek et al. (2019), who proposed that the dust forms very rapidly, and because of its optical depth, is not initially observationally detectable, only being gradually revealed as the ejecta expand. Observational dust masses at early times would then only be lower limits. Using a large grid of radiative transfer models covering dust masses from 10$^{-4}$ to 1M$_\odot$ to calculate both the spectral energy distribution and the emission line profiles from clumpy dust shells, we show that this cannot be the case. Some clump distributions allow dust masses of $\sim$0.01 M$_\odot$ to be concealed in clumps and still predict an SED consistent with the observations. However, these geometries predict emission line profiles which are inconsistent with the observations. Similarly, clump geometries which reproduce the observed emission line profiles with dust masses >0.01 M$_\odot$ do not reproduce the SED. However, models with $\sim$10$^{-3}$ M$_\odot$ of amorphous carbon can reproduce both the SED and the emission line profiles. We conclude that no large masses of dust can be hidden from view in the ejecta of SN 1987A at early epochs, and that the majority of dust must thus have formed at epochs >1000 days.


title: Weiss variation for general boundaries

abstract: The Weiss variation of the Einstein-Hilbert action with an appropriate boundary term has been studied for general boundary surfaces; the boundary surfaces can be spacelike, timelike, or null. To achieve this we introduce an auxiliary reference connection and find that the resulting Weiss variation yields the Einstein equations as expected, with additional boundary contributions. Among these boundary contributions, we obtain the dynamical variable and the associated conjugate momentum, irrespective of the spacelike, timelike or, null nature of the boundary surface. We also arrive at the generally non-vanishing covariant generalization of the Einstein energy-momentum pseudotensor. We study this tensor in the Schwarzschild geometry and find that the pseudotensorial ambiguities translate into ambiguities in the choice of coordinates on the reference geometry. Moreover, we show that from the Weiss variation, one can formally derive a gravitational Schr{\"o}dinger equation, which may, despite ambiguities in the definition of the Hamiltonian, be useful as a tool for studying the problem of time in quantum general relativity. Implications have been discussed.


title: Dirac Radiative Neutrino Mass with Modular Symmetry and Leptogenesis

abstract: Minimalistic Dirac radiative neutrino mass model based on modular symmetry is proposed. We predict maximum number of observables possible including neutrino mass splittings, neutrino mass scale, lepton mixing angles, and Dirac phases in the leptonic sector with minimum number of input parameters possible. Model is capable of accommodating multicomponent dark matter, thanks to the $\emph{R}-$parity and accidental scotogenic $\mathbb{Z}_2$ discrete symmetry. Furthermore, even-though neutrinos are Dirac in our model, matter-antimatter asymmetry of the Universe is achieved via neutrinogenesis mechanism. Phenomenology of the dark sector including various dark matter candidates is briefly discussed.


title: RR Lyrae-based Distances for 39 Nearby Dwarf Galaxies Calibrated to Gaia eDR3

abstract: We provide uniform RR Lyrae-based distances to 39 dwarf galaxies in and around the Local Group. We determine distances based on a Bayesian hierarchical model that uses periods and magnitudes of published RR Lyrae in dwarf galaxies and is anchored to well-studied Milky Way (MW) RR Lyrae with spectroscopic metallicities and Gaia eDR3 parallaxes. Gaia eDR3 parallaxes for the anchor sample are a factor of 2, on average, more precise than DR2 parallaxes, and allow for a much better constrained period-luminosity-metallicity relation. While $\sim75$% of our distances are within 1-$\sigma$ of recent literature RR Lyrae distances, our distances are also $\sim2$-$3$ times more precise than literature distances, on average. On average, our distances are $\sim0.05$ mag closer than literature distances, as well as $\sim0.06$ mag closer than distances derived using a theoretical period-luminosity-metallicity relation. These discrepancies are largely due to our eDR3 parallax anchor. We show that the Hipparcos-anchored RR Lyrae distance scale of Carretta et al. (2000) over-predicts distances to MW RR Lyrae by $\sim0.05$ mag. The largest uncertainties in our distances are (i) the lack of direct metallicity measurements for RR Lyrae and (ii) the heterogeneity of published RR Lyrae photometry. We provide simple formulae to place new dwarf galaxies with RR Lyrae on a common distance scale with this work. We provide public code that can easily incorporate additional galaxies and data from future surveys, providing a versatile framework for cartography of the local Universe with RR Lyrae.


title: Constraining the cosmic-ray ionization rate and their spectrum with NIR spectroscopy of dense clouds -- A test-bed for JWST

abstract: Low-energy cosmic-rays (CRs) control the thermo-chemical state and the coupling between gas and magnetic fields in dense molecular clouds, the sites of star-formation. However, current estimates of the low-energy CR spectrum ($E \lesssim 1$ GeV) and the associated CR ionization rate are highly uncertain. We apply, for the first time, a new method for constraining the CR ionization rate and the CR spectral shape using H$_2$ rovibrational lines from cold molecular clouds. Using the MMIRS instrument on the MMT, we obtained deep near-infrared (NIR) spectra in six positions within four dense cores, G150, G157, G163, G198, with column densities $N_{\rm H_2} \approx 10^{22}$ cm$^{-2}$. We derive 3$\sigma$ upper limits on the H$_2$ $(1-0)$S(0) line (2.22 $\mu$m) brightness in the range $I = 5.9 \times 10^{-8}$ to $1.2 \times 10^{-7}$ erg cm$^{-2}$ s$^{-1}$ sr$^{-1}$ for the different targets. Using both an analytic model and a numerical model of CR propagation, we convert these into upper limits on the CR ionization rate in the clouds' interior, $\zeta = 1.5$ to $3.6 \times 10^{-16}$ s$^{-1}$, and lower limits on the low-energy spectral slope of interstellar CR protons, $\alpha = -0.97$ to $-0.79$. We show that while MMT was unable to detect the H$_2$ lines due to high atmospheric noise, JWST/NIRSpec will be able to efficiently detect the CR-excited H$_2$ lines, making it the ideal method for constraining the otherwise elusive low-energy CRs, shedding light on the sources and propagation modes of CRs.


title: Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

abstract: Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}\mu$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretisation of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following $\propto \Dres^{-1/2}$, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy.


title: MS-LaTTE: A Dataset of Where and When To-do Tasks are Completed

abstract: Tasks are a fundamental unit of work in the daily lives of people, who are increasingly using digital means to keep track of, organize, triage and act on them. These digital tools -- such as task management applications -- provide a unique opportunity to study and understand tasks and their connection to the real world, and through intelligent assistance, help people be more productive. By logging signals such as text, timestamp information, and social connectivity graphs, an increasingly rich and detailed picture of how tasks are created and organized, what makes them important, and who acts on them, can be progressively developed. Yet the context around actual task completion remains fuzzy, due to the basic disconnect between actions taken in the real world and telemetry recorded in the digital world. Thus, in this paper we compile and release a novel, real-life, large-scale dataset called MS-LaTTE that captures two core aspects of the context surrounding task completion: location and time. We describe our annotation framework and conduct a number of analyses on the data that were collected, demonstrating that it captures intuitive contextual properties for common tasks. Finally, we test the dataset on the two problems of predicting spatial and temporal task co-occurrence, concluding that predictors for co-location and co-time are both learnable, with a BERT fine-tuned model outperforming several other baselines. The MS-LaTTE dataset provides an opportunity to tackle many new modeling challenges in contextual task understanding and we hope that its release will spur future research in task intelligence more broadly.


title: 4d N=2 supergravity observables from Nekrasov-like partition functions

abstract: We reinterpret the OSV formula for the on-shell action/entropy function of asymptotically flat BPS black holes as a fixed point formula that is formally equivalent to a recent gluing proposal for asymptotically AdS$_4$ black holes. This prompts a conjecture that the complete perturbative answer for the most general gravitational building block of 4d $\mathcal{N}=2$ supergravity at a single fixed point takes the form of a Nekrasov-like partition function with equivariant parameters related to the higher-derivative expansion of the prepotential. In turn this leads to a simple localization-like proposal for a set of supersymmetric partition functions in (UV completed) 4d $\mathcal{N}=2$ supergravity theories. The conjecture is shown to be in agreement with a number of available results for different BPS backgrounds with both Minkowski and AdS asymptotics. In particular, it follows that the OSV formula comes from the unrefined limit of the general expression including only the so-called $\mathbb{W}$ tower of higher derivatives, while the on-shell action of pure (Euclidean) AdS$_4$ with round S$^3$ boundary comes from the NS limit that includes only the $\mathbb{T}$ tower. Backgrounds preserving less supersymmetry, such as the under-rotating black holes in flat space, the holographic squashed S$^3$, and the static/rotating twisted and non-twisted Kerr-Newman-like black holes in AdS$_4$ lead to a more general refined version of the corresponding gravitational blocks as dictated by the supersymmetric gluing rules.


title: Triaxial Orbit-based Dynamical Modeling of Galaxies with Supermassive Black Holes and an Application to Massive Elliptical Galaxy NGC 1453

abstract: Most stellar-dynamical determinations of the masses of nearby supermassive black holes (SMBHs) have been obtained with the orbit superposition technique under the assumption of axisymmetry. However, few galaxies -- in particular massive early-type galaxies -- obey exact axisymmetry. Here we present a revised orbit superposition code and a new approach for dynamically determining the intrinsic shapes and mass parameters of triaxial galaxies based on spatially-resolved stellar kinematic data. The triaxial TriOS code described here corrects an error in the original van den Bosch et al. (2008) code that gives rise to incorrect projections for most orbits in triaxial models and can significantly impact parameter search results. The revised code also contains significant improvements in orbit sampling, mass constraints, and run time. Furthermore, we introduce two new parameter searching strategies -- a new set of triaxial shape parameters and a novel grid-free sampling technique -- that together lead to a remarkable gain in efficiency in locating the best-fit model. We apply the updated code and search method to NGC 1453, a fast-rotating massive elliptical galaxy. A full 6D parameter search finds $p=b/a=0.933^{+0.014}_{-0.015}$ and $q=c/a=0.779\pm0.012$ for the intrinsic axis ratios and $T=0.33\pm0.06$ for the triaxiality parameter. Despite the deviations from axisymmetry, the best-fit SMBH mass, stellar mass-to-light ratio, and dark matter enclosed mass for NGC 1453 are consistent with the axisymmetric results. More comparisons between axisymmetric and triaxial modeling are needed before drawing general conclusions.


title: Revised Stellar Parameters for V471 Tau, A Post-common Envelope Binary in the Hyades

abstract: V471 Tau is a post-common-envelope binary consisting of an eclipsing DA white dwarf and a K-type main-sequence star in the Hyades star cluster. We analyzed publicly available photometry and spectroscopy of V471 Tau to revise the stellar and orbital parameters of the system. We used archival K2 photometry, archival Hubble Space Telescope spectroscopy, and published radial-velocity measurements of the K-type star. Employing Gaussian processes to fit for rotational modulation of the system flux by the main-sequence star, we recovered the transits of the white dwarf in front of the main-sequence star for the first time. The transits are shallower than would be expected from purely geometric occultations owing to gravitational microlensing during transit, which places an additional constraint on the white-dwarf mass. Our revised mass and radius for the main-sequence star is consistent with single-star evolutionary models given the age and metallicity of the Hyades. However, as noted previously in the literature, the white dwarf is too massive and too hot to be the result of single-star evolution given the age of the Hyades, and may be the product of a merger scenario. We independently estimate the conditions of the system at the time of common envelope that would result in the measured orbital parameters today.


title: Path Verification for Dynamic Indirect Illumination

abstract: In this paper we present a technique that improves rendering performance for real-time scenes with ray traced lighting in the presence of dynamic lights and objects. In particular we verify photon paths from the previous frame against dynamic objects in the current frame, and show how most photon paths are still valid. When using area lights, we use a data structure to store light distribution that tracks light paths allowing photons to be reused when the light source is moving in the scene. We also show that by reusing paths when the error in the reflected energy is below a threshold value, even more paths can be reused. We apply this technique to Indirect Illumination using a screen space photon splatting rendering engine. By reusing photon paths and applying our error threshold, our method can reduce the number of rays traced by up to 5x, and improve performance by up to 2x.


title: Improving Experience Replay through Modeling of Similar Transitions' Sets

abstract: In this work, we propose and evaluate a new reinforcement learning method, COMPact Experience Replay (COMPER), which uses temporal difference learning with predicted target values based on recurrence over sets of similar transitions, and a new approach for experience replay based on two transitions memories. Our objective is to reduce the required number of experiences to agent training regarding the total accumulated rewarding in the long run. Its relevance to reinforcement learning is related to the small number of observations that it needs to achieve results similar to that obtained by relevant methods in the literature, that generally demand millions of video frames to train an agent on the Atari 2600 games. We report detailed results from five training trials of COMPER for just 100,000 frames and about 25,000 iterations with a small experiences memory on eight challenging games of Arcade Learning Environment (ALE). We also present results for a DQN agent with the same experimental protocol on the same games set as the baseline. To verify the performance of COMPER on approximating a good policy from a smaller number of observations, we also compare its results with that obtained from millions of frames presented on the benchmark of ALE.


title: Explainable AI for Psychological Profiling from Digital Footprints: A Case Study of Big Five Personality Predictions from Spending Data

abstract: Every step we take in the digital world leaves behind a record of our behavior; a digital footprint. Research has suggested that algorithms can translate these digital footprints into accurate estimates of psychological characteristics, including personality traits, mental health or intelligence. The mechanisms by which AI generates these insights, however, often remain opaque. In this paper, we show how Explainable AI (XAI) can help domain experts and data subjects validate, question, and improve models that classify psychological traits from digital footprints. We elaborate on two popular XAI methods (rule extraction and counterfactual explanations) in the context of Big Five personality predictions (traits and facets) from financial transactions data (N = 6,408). First, we demonstrate how global rule extraction sheds light on the spending patterns identified by the model as most predictive for personality, and discuss how these rules can be used to explain, validate, and improve the model. Second, we implement local rule extraction to show that individuals are assigned to personality classes because of their unique financial behavior, and that there exists a positive link between the model's prediction confidence and the number of features that contributed to the prediction. Our experiments highlight the importance of both global and local XAI methods. By better understanding how predictive models work in general as well as how they derive an outcome for a particular person, XAI promotes accountability in a world in which AI impacts the lives of billions of people around the world.


title: Active information requirements for fixation on the Wright-Fisher model of population genetics

abstract: In the context of population genetics, active information can be extended to measure the change of information of a given event (e.g., fixation of an allele) from a neutral model in which only genetic drift is taken into account to a non-neutral model that includes other sources of frequency variation (e.g., selection and mutation). In this paper we illustrate active information in population genetics through the Wright-Fisher model.

categories: cs.IT math.IT math.PR physics.bio-ph q-bio.PE


title: Optimal fine-scale structures in compliance minimization for a uniaxial load in three space dimensions

abstract: We consider the shape and topology optimization problem to design a structure that minimizes a weighted sum of material consumption and (linearly) elastic compliance under a fixed given boundary load. As is well-known, this problem is in general not well-posed since its solution typically requires the use of infinitesimally fine microstructure. Therefore we examine the effect of singularly perturbing the problem by adding the structure perimeter to the cost. For a uniaxial and a shear load in two space dimensions, corresponding energy scaling laws were already derived in the literature. This work now derives the scaling law for the case of a uniaxial load in three space dimensions, which can be considered the simplest three-dimensional setting. In essence, it is expected (and confirmed in this article) that for a uniaxial load the compliance behaves almost like the dissipation in a scalar flux problem so that lower bounds from pattern analysis in superconductors can directly be applied. The upper bounds though require nontrivial modifications of the constructions known from superconductors. Those become necessary since in elasticity one has the additional constraint of torque balance.


title: On fiber bundles and quaternionic slice regular functions

abstract: The papers \cite{O1,O2} are the first works to apply the theory of fiber bundles in the study of the quaternionic slice regular functions. The main goal of the present work is to extend the results given in \cite{O1}, where the quaternionic right linear space of quaternionic slice regular functions was presented as the base space of a fiber bundle. When the quaternionic right linear space of quaternionic slice regular functions is associated to certain domains then this paper shows that the elements of total space, given in \cite{O1}, are defined from a pair of harmonic functions and a pair of orthogonal vectors. Simplifying the computations presented in \cite{O1}, where each element of the total space is formed by two pair of conjugate harmonic functions and a pair of orthogonal unit vectors. This work also gives some interpretations of the behavior of the zero sets of some quaternionic slice regular polynomials in terms of the theory of fiber bundles.


title: Microcanonical Action and the Entropy of Hawking Radiation

abstract: The island formula -- an extremization prescription for generalized entropy -- is known to result in a unitary Page curve for the entropy of Hawking radiation. This semi-classical entropy formula has been derived for Jackiw-Teitelboim (JT) gravity coupled to conformal matter using the "replica trick" to evaluate the Euclidean path integral. Alternatively, for eternal Anti-de Sitter black holes, we derive the extremization of generalized entropy from minimizing the microcanonical action of an entanglement wedge. The on-shell action is minus the entropy and arises in the saddle-point approximation of the (nonreplicated) microcanonical path integral. When the black hole is coupled to a bath, islands emerge from maximizing the entropy at fixed energy, consistent with the island formula. Our method applies to JT gravity as well as other two-dimensional dilaton gravity theories.


title: Visual Intelligence through Human Interaction

abstract: Over the last decade, Computer Vision, the branch of Artificial Intelligence aimed at understanding the visual world, has evolved from simply recognizing objects in images to describing pictures, answering questions about images, aiding robots maneuver around physical spaces and even generating novel visual content. As these tasks and applications have modernized, so too has the reliance on more data, either for model training or for evaluation. In this chapter, we demonstrate that novel interaction strategies can enable new forms of data collection and evaluation for Computer Vision. First, we present a crowdsourcing interface for speeding up paid data collection by an order of magnitude, feeding the data-hungry nature of modern vision models. Second, we explore a method to increase volunteer contributions using automated social interventions. Third, we develop a system to ensure human evaluation of generative vision models are reliable, affordable and grounded in psychophysics theory. We conclude with future opportunities for Human-Computer Interaction to aid Computer Vision.


title: Elastic Properties of Confined Fluids from Molecular Modeling to Ultrasonic Experiments on Porous Solids

abstract: Fluids confined in nanopores are ubiquitous in nature and technology. In recent years, the interest in confined fluids has grown, driven by research on unconventional hydrocarbon resources -- shale gas and shale oil, much of which are confined in nanopores. When fluids are confined in nanopores, many of their properties differ from those of the same fluid in the bulk. These properties include density, freezing point, transport coefficients, thermal expansion coefficient, and elastic properties. The elastic moduli of a fluid confined in the pores contribute to the overall elasticity of the fluid-saturated porous medium and determine the speed at which elastic waves traverse through the medium. Wave propagation in fluid-saturated porous media is pivotal for geophysics, as elastic waves are used for characterization of formations and rock samples. In this paper, we present a comprehensive review of experimental works on wave propagation in fluid-saturated nanoporous media, as well as theoretical works focused on calculation of compressibility of fluids in confinement. We discuss models that bridge the gap between experiments and theory, revealing a number of open questions that are both fundamental and applied in nature. While some results were demonstrated both experimentally and theoretically (e.g. the pressure dependence of compressibility of fluids), others were theoretically predicted, but not verified in experiments (e.g. linear scaling of modulus with the pore size). Therefore, there is a demand for the combined experimental-modeling studies on porous samples with various characteristic pore sizes. The extension of molecular simulation studies from simple model fluids to the more complex molecular fluids is another open area of practical interest.

categories: cond-mat.soft physics.app-ph physics.chem-ph physics.geo-ph


title: First-principles investigations of the electronic, magnetic and thermoelectric properties of VTiRhZ (Z= Al, Ga, In) Quaternary Heusler Alloys

abstract: Calculations using density functional theory (DFT) were performed to investigate the structural, dynamical, mechanical, electronic, magnetic, and thermoelectric properties of VTiRhZ (Z = Al, Ga, In) alloys. The most stable structure of these alloys was found to be the type-I configuration. Using GGA-PBE functional, VTiRhGa, and VTiRhIn alloys are predicted as half-metallic ferromagnets with a 100% spin-polarization and a total magnetic moment of 3{\mu}B, which is promising for spintronic applications. The thermoelectric properties and lattice thermal conductivity of VTiRhZ alloys were obtained using the Boltzmann transport theory within the constant relaxation time and Slack equation, respectively. The figure-of-merit (ZT) values of VTiRhAl, VTiRhGa, and VTiRhIn alloys were found to be 0.96, 0.88 and 0.64, respectively, which are promising for future thermoelectric applications.


title: Offense Detection in Dravidian Languages using Code-Mixing Index based Focal Loss

abstract: Over the past decade, we have seen exponential growth in online content fueled by social media platforms. Data generation of this scale comes with the caveat of insurmountable offensive content in it. The complexity of identifying offensive content is exacerbated by the usage of multiple modalities (image, language, etc.), code-mixed language and more. Moreover, even after careful sampling and annotation of offensive content, there will always exist a significant class imbalance between offensive and non-offensive content. In this paper, we introduce a novel Code-Mixing Index (CMI) based focal loss which circumvents two challenges (1) code-mixing in languages (2) class imbalance problem for Dravidian language offense detection. We also replace the conventional dot product-based classifier with the cosine-based classifier which results in a boost in performance. Further, we use multilingual models that help transfer characteristics learnt across languages to work effectively with low resourced languages. It is also important to note that our model handles instances of mixed script (say usage of Latin and Dravidian-Tamil script) as well. To summarize, our model can handle offensive language detection in a low-resource, class imbalanced, multilingual and code-mixed setting.


title: Positive periodic solutions for systems of impulsive delay differential equations

abstract: A class of periodic differential $n$-dimensional systems with patch structure with (possibly infinite) delay and nonlinear impulses is considered. These systems incorporate very general nonlinearities and impulses whose signs may vary. Criteria for the existence of at least one positive periodic solution are presented, extending and improving previous ones established for the scalar case. Applications to systems inspired in mathematical biology models, such as impulsive hematopoiesis and Nicholson-type systems, are also included.


title: Regular Path Query Evaluation Sharing a Reduced Transitive Closure Based on Graph Reduction

abstract: Regular path queries (RPQs) find pairs of vertices of paths satisfying given regular expressions on an edge-labeled, directed multigraph. When evaluating an RPQ, the evaluation of a Kleene closure (i.e., Kleene plus or Kleene star) is very expensive. Furthermore, when multiple RPQs include a Kleene closure as a common sub-query, repeated evaluations of the common sub-query cause serious performance degradation. In this paper, we present a novel concept of RPQ-based graph reduction, which significantly simplifies the original graph through edge-level and vertex-level reductions. Interestingly, RPQ-based graph reduction can replace the evaluation of the Kleene closure on the large original graph to that of the transitive closure to the small reduced graph. We then propose a reduced transitive closure (RTC) as a lightweight structure for efficiently sharing the result of a Kleene closure. We also present an RPQ evaluation algorithm, RTCSharing, which treats each clause in the disjunctive normal form of the given RPQ as a batch unit. If the batch units include a Kleene closure as a common sub-query, we share the lightweight RTC instead of the heavyweight result of the Kleene closure. RPQ-based graph reduction further enables us to formally represent the result of an RPQ including a Kleene closure as a relational algebra expression including the RTC. Through the formal expression, we optimize the evaluation of the batch unit by eliminating useless and redundant operations of the previous method. Experiments show that RTCSharing improves the performance significantly by up to 73.86 times compared with existing methods in terms of query response time.


title: Explicit Factorization of a Categorical Center

abstract: Given a braided fusion category $C$, it is well known that the natural map $C \boxtimes C^{bop} \to Z(C)$ from the square of $C$ to the (Drinfeld) categorical center $Z(C)$ is an equivalence if and only if $C$ is modular. However, it is not clear how to construct the inverse and the natural isomorphisms. In this work, we provide an explicit construction using insights from a specific quantum field theory, and explore how the equivalence fails for the degenerate cases.


title: Neural optimal feedback control with local learning rules

abstract: A major problem in motor control is understanding how the brain plans and executes proper movements in the face of delayed and noisy stimuli. A prominent framework for addressing such control problems is Optimal Feedback Control (OFC). OFC generates control actions that optimize behaviorally relevant criteria by integrating noisy sensory stimuli and the predictions of an internal model using the Kalman filter or its extensions. However, a satisfactory neural model of Kalman filtering and control is lacking because existing proposals have the following limitations: not considering the delay of sensory feedback, training in alternating phases, and requiring knowledge of the noise covariance matrices, as well as that of systems dynamics. Moreover, the majority of these studies considered Kalman filtering in isolation, and not jointly with control. To address these shortcomings, we introduce a novel online algorithm which combines adaptive Kalman filtering with a model free control approach (i.e., policy gradient algorithm). We implement this algorithm in a biologically plausible neural network with local synaptic plasticity rules. This network performs system identification and Kalman filtering, without the need for multiple phases with distinct update rules or the knowledge of the noise covariances. It can perform state estimation with delayed sensory feedback, with the help of an internal model. It learns the control policy without requiring any knowledge of the dynamics, thus avoiding the need for weight transport. In this way, our implementation of OFC solves the credit assignment problem needed to produce the appropriate sensory-motor control in the presence of stimulus delay.


title: Performance Analysis over Correlated/Independent Fisher-Snedecor F Fading Multi-User Channels

abstract: In this paper, we investigate the impact of correlated fading on the performance of wireless multiple access channels (MAC) in the presence and absence of side information (SI) at transmitters, where the fading coefficients are modeled according to the Fisher-Snedecor F distribution. Specifically, we represent two scenarios: (i) clean MAC (i.e, without SI at transmitters), (ii) doubly dirty MAC (i.e., with the non-causally known SI at transmitters). For both system models, we derive the closedform expressions of the outage probability (OP) as well as the average capacity (AC) under independent fading conditions. Besides, exploiting copula theory, we obtain the exact analytical expressions for the OP and the AC under positive dependence fading conditions in both considered models. Finally, the validity of the analytical results is illustrated numerically.


title: Morse Index Bound for Minimal Torus

abstract: The min-max construction of minimal spheres using harmonic replacement is introduced by Colding and Minicozzi and generalized by Zhou to conformal harmonic torus. We prove that the Morse index of the min-max conformal harmonic torus is bounded by one.


title: ARC Nav -- A 3D Navigation Stack for Autonomous Robots

abstract: Popular navigation stacks implemented on top of open-source frameworks such as ROS(Robot Operating System) and ROS2 represent the robot workspace using a discretized 2D occupancy grid. This method, while requiring less computation, restricts the use of such navigation stacks to wheeled robots navigating on flat surfaces. In this paper, we present a navigation stack that uses a volumetric representation of the robot workspace, and hence can be extended to aerial and legged robots navigating through uneven terrain. Additionally, we present a new sampling-based motion planning algorithm which introduces a bi-directional approach to the Batch Informed Trees (BIT*) motion planning algorithm, whilst wrapping it with a strategy switching approach in order to reduce the initial time taken to find a path, in addition to the time taken to find the shortest path.


title: A Simple and Fast Baseline for Tuning Large XGBoost Models

abstract: XGBoost, a scalable tree boosting algorithm, has proven effective for many prediction tasks of practical interest, especially using tabular datasets. Hyperparameter tuning can further improve the predictive performance, but unlike neural networks, full-batch training of many models on large datasets can be time consuming. Owing to the discovery that (i) there is a strong linear relation between dataset size & training time, (ii) XGBoost models satisfy the ranking hypothesis, and (iii) lower-fidelity models can discover promising hyperparameter configurations, we show that uniform subsampling makes for a simple yet fast baseline to speed up the tuning of large XGBoost models using multi-fidelity hyperparameter optimization with data subsets as the fidelity dimension. We demonstrate the effectiveness of this baseline on large-scale tabular datasets ranging from $15-70\mathrm{GB}$ in size.


title: Action2video: Generating Videos of Human 3D Actions

abstract: We aim to tackle the interesting yet challenging problem of generating videos of diverse and natural human motions from prescribed action categories. The key issue lies in the ability to synthesize multiple distinct motion sequences that are realistic in their visual appearances. It is achieved in this paper by a two-step process that maintains internal 3D pose and shape representations, action2motion and motion2video. Action2motion stochastically generates plausible 3D pose sequences of a prescribed action category, which are processed and rendered by motion2video to form 2D videos. Specifically, the Lie algebraic theory is engaged in representing natural human motions following the physical law of human kinematics; a temporal variational auto-encoder (VAE) is developed that encourages diversity of output motions. Moreover, given an additional input image of a clothed human character, an entire pipeline is proposed to extract his/her 3D detailed shape, and to render in videos the plausible motions from different views. This is realized by improving existing methods to extract 3D human shapes and textures from single 2D images, rigging, animating, and rendering to form 2D videos of human motions. It also necessitates the curation and reannotation of 3D human motion datasets for training purpose. Thorough empirical experiments including ablation study, qualitative and quantitative evaluations manifest the applicability of our approach, and demonstrate its competitiveness in addressing related tasks, where components of our approach are compared favorably to the state-of-the-arts.


title: The unitary Cuntz semigroup on the classification of non-simple C*-algebras

abstract: This paper argues that the unitary Cuntz semigroup, introduced in [10] and termed Cu$_1$, contains crucial information regarding the classification of non-simple C$^*$-algebras. We exhibit two (non-simple) C$^*$-algebras that agree on their Cuntz semigroups, termed Cu, and their K$_1$-groups and yet disagree at level of their unitary Cuntz semigroups. In the process, we establish that the unitary Cuntz semigroup contains rigorously more information about non-simple C$^*$-algebras than Cu and K$_1$ alone.


title: Improving blazar redshift constraints with the edge of the Ly$\alpha$ forest: 1ES 1553+113 and implications for observations of the WHIM

abstract: Blazars are some of the brightest UV and X-ray sources in the sky and are valuable probes of the elusive warm-hot intergalactic medium (WHIM; $T{\simeq} 10^5-10^7$ K). However, many of the brightest blazars$-$called BL Lac objects such as 1ES1553+113$-$have quasi-featureless spectra and poorly constrained redshifts. Here, we significantly improve the precision of indirect redshift constraints for blazars based on the edge of the $\rm{H\,I}$ Ly$\alpha$ forest observed in their UV spectra. We develop a robust technique to constrain the redshift of a $z<0.5$ AGN or blazar with a $1\sigma$ uncertainty of ${\approx}0.01$ using only the position of its highest-redshift Ly$\alpha$ absorber with $\log N_{\rm{H\,I}}/{\rm cm^{-2}} > 12.6$. We use a large sample of 192 AGN/QSOs at $0.01\lesssim z\lesssim0.45$ that have high-quality COS FUV spectra to characterize the intrinsic scatter in the gap between the AGN redshift and the edge of their Ly$\alpha$ forest. We present new COS NUV data for 1ES1553+113 and confirm its redshift of $z=0.433$ using our technique. We apply our Ly$\alpha$-forest-based redshift estimation technique to nine additional blazars with archival ${\it HST}$ UV spectra, most of which are key targets for future X-ray missions. Our inferred redshift constraints improve estimates for two BL Lacs (1ES1118+424 and S50716+714) and are consistent with previous estimates for the rest. Our results emphasize the need to obtain further UV spectra of bright blazars, of which many have uncertain redshifts, in order to maximize the scientific value of future X-ray WHIM observations that will improve our understanding of galaxy evolution.


title: Generalized Nested Rollout Policy Adaptation with Dynamic Bias for Vehicle Routing

abstract: In this paper we present an extension of the Nested Rollout Policy Adaptation algorithm (NRPA), namely the Generalized Nested Rollout Policy Adaptation (GNRPA), as well as its use for solving some instances of the Vehicle Routing Problem. We detail some results obtained on the Solomon instances set which is a conventional benchmark for the Vehicle Routing Problem (VRP). We show that on all instances, GNRPA performs better than NRPA. On some instances, it performs better than the Google OR Tool module dedicated to VRP.


title: Hierarchical Bayesian Bandits

abstract: Meta-, multi-task, and federated learning can be all viewed as solving similar tasks, drawn from a distribution that reflects task similarities. We provide a unified view of all these problems, as learning to act in a hierarchical Bayesian bandit. We propose and analyze a natural hierarchical Thompson sampling algorithm (HierTS) for this class of problems. Our regret bounds hold for many variants of the problems, including when the tasks are solved sequentially or in parallel; and show that the regret decreases with a more informative prior. Our proofs rely on a novel total variance decomposition that can be applied beyond our models. Our theory is complemented by experiments, which show that the hierarchy helps with knowledge sharing among the tasks. This confirms that hierarchical Bayesian bandits are a universal and statistically-efficient tool for learning to act with similar bandit tasks.


title: Quantum teleportation via a two-qubit Heisenberg XXX chain with x-component of Dzyaloshinskii-Moriya interaction

abstract: In this paper, we investigate the thermal entanglement and teleportation of a thermally mixed entangled via a two-qubit Heisenberg XXX chain with the x-component $D_x$ of the Dzyaloshinskii-Moriya interaction. Our findings suggest that temperature $T$, spin coupling constant J, and the x-components $D_x$ may influence entanglement of the output states and, consequently, the possibilities of teleportation protocols. Furthermore, these results indicate that the entanglement of the output states requires a low-temperature regime, weak Dzyaloshinskii-Moriya interaction, or an antiferromagnetic chain. Finally, the channel becomes entangled, making the teleportation protocol conceivable and feasible.


title: Solving A System Of Linear Equations By Randomized Orthogonal Projections

abstract: Randomization has shown catalyzing effects in linear algebra with promising perspectives for tackling computational challenges in large-scale problems. For solving a system of linear equations, we demonstrate the convergence of a broad class of algorithms that at each step pick a subset of $n$ equations at random and update the iterate with the orthogonal projection to the subspace those equations represent. We identify, in this context, a specific degree-$n$ polynomial that non-linearly transforms the singular values of the system towards equalization. This transformation to singular values and the corresponding condition number then characterizes the expected convergence rate of iterations. As a consequence, our results specify the convergence rate of the stochastic gradient descent algorithm, in terms of the mini-batch size $n$, when used for solving systems of linear equations.


title: A Cost-Effective Upgrade Path for the Fermilab Accelerator Complex

abstract: The Fermilab Proton Improvement Plan II, or PIP-II, would enable the world's most intense high-energy neutrino beam and would help scientists search for rare particle physics processes. The PIP-II goal is to deliver 1.2 MW of proton beam power from the Fermilab Main Injector, over the energy range 60 - 120 GeV, at the start of operation of the LBNF/DUNE program. PIP-II provides a variety of upgrade paths to higher beam power from the Main Injector, as demanded by the neutrino science program and as recommended by the 2014 P5 report. Delivering more than 2 MW to the LBNF target in the future will require a replacement of the existing Booster. This report outlines a cost-effective Booster replacement option and an upgrade path for the Fermilab Accelerator Complex to attain 2.4-MW beam power on the LBNF target, as well as to retain the capability to provide 8-GeV proton beams to the existing Fermilab Muon Campus via the existing Recycler ring. Its cost-effectiveness is achieved by: (1) using a small-diameter metallic vacuum chamber in the Booster replacement and (2) reusing the existing Recycler ring. Reusing the Recycler ring may be of particular advantage since it is presently employed to deliver 8-GeV beams to the Muon Campus experiments. The present concept also retains such a capability.


title: Nanometer-scale lateral p-n junctions in graphene/$\alpha$-RuCl$_3$ heterostructures

abstract: The ability to create high-quality lateral p-n junctions at nanometer length scales is essential for the next generation of two-dimensional (2D) electronic and plasmonic devices. Using a charge-transfer heterostructure consisting of graphene on $\alpha$-RuCl$_3$, we conduct a proof-of-concept study demonstrating the existence of intrinsic nanoscale lateral p-n junctions in the vicinity of graphene nanobubbles. Our multi-pronged experimental approach incorporates scanning tunneling microscopy (STM) and spectroscopy (STS) and scattering-type scanning near-field optical microscopy ($\textit{s}$-SNOM) in order to simultaneously probe both the electronic and optical responses of nanobubble p-n junctions. Our STM and STS results reveal that p-n junctions with a band offset of more than 0.6 eV can be achieved over lateral length scale of less than 3 nm, giving rise to a staggering effective in-plane field in excess of 10$^8$ V/m. Concurrent $\textit{s}$-SNOM measurements confirm the utility of these nano-junctions in plasmonically-active media, and validate the use of a point-scatterer formalism for modeling surface plasmon polaritons (SPPs). Model $\textit{ab initio}$ density functional theory (DFT) calculations corroborate our experimental data and reveal a combination of sub-angstrom and few-angstrom decay processes dictating the dependence of charge transfer on layer separation. Our study provides experimental and conceptual foundations for the use of charge-transfer interfaces such as graphene/$\alpha$-RuCl$_3$ to generate p-n nano-junctions.


title: Contrastive Feature Loss for Image Prediction

abstract: Training supervised image synthesis models requires a critic to compare two images: the ground truth to the result. Yet, this basic functionality remains an open problem. A popular line of approaches uses the L1 (mean absolute error) loss, either in the pixel or the feature space of pretrained deep networks. However, we observe that these losses tend to produce overly blurry and grey images, and other techniques such as GANs need to be employed to fight these artifacts. In this work, we introduce an information theory based approach to measuring similarity between two images. We argue that a good reconstruction should have high mutual information with the ground truth. This view enables learning a lightweight critic to "calibrate" a feature space in a contrastive manner, such that reconstructions of corresponding spatial patches are brought together, while other patches are repulsed. We show that our formulation immediately boosts the perceptual realism of output images when used as a drop-in replacement for the L1 loss, with or without an additional GAN loss.


title: Quantum corrected thermodynamics of nonlinearly charged BTZ black holes in massive gravity's rainbow

abstract: In this paper, we consider three-dimensional massive gravity's rainbow and obtain black hole solutions in three different cases of Born-Infeld, logarithmic, and exponential theories of nonlinear electrodynamics. We discuss the horizon structure and geometrical properties. Then, we study thermodynamics of these models by considering the first-order quantum correction effects, which appears as a logarithmic term in the black hole entropy. We discuss such effects on the black hole stability and phase transitions. We find that, due to the quantum corrections, the second-order phase transition happens in Born-Infeld and logarithmic models. We obtain the modified first law of black hole thermodynamics in the presence of logarithmic corrections.


title: Whirling interlayer fields as a source of stable topological order in moir\'e CrI3

abstract: The moir\'e engineering of two-dimensional magnets opens unprecedented opportunities to design novel magnetic states with promises for spintronic device applications. The possibility of stabilizing skyrmions in these materials without chiral spin-orbit couplings or dipolar interactions is yet to be explored. Here, we investigate the formation and control of ground state topological spin textures (TSTs) in moir\'e CrI3 using stochastic Landau-Lifshitz-Gilbert simulations. We unveil the emergence of interlayer vortex and antivortex Heisenberg exchange fields, stabilizing spontaneous and field-assisted ground state TSTs with various topologies. The developed study accounts for the full bilayer spin dynamics, thermal fluctuations, and intrinsic spin-orbit couplings. By examining the effect of the Kitaev interaction and the next nearest-neighbor Dzyaloshinskii-Moriya interaction, we propose the latter as the unique spin-orbit coupling mechanism compatible with experiments on monolayer and twisted CrI3. Our findings contribute to the current knowledge about moir\'e skyrmionics and uncover the nature of spin-orbit coupling in CrI3.