title: Disorder-Free Localization and Many-Body Quantum Scars from Magnetic Frustration

abstract: The concept of geometrical frustration has led to rich insights into condensed matter physics, especially as a mechansim to produce exotic low energy states of matter. Here we show that frustration provides a natural vehicle to generate models exhibiting anomalous thermalization of various types within high energy states. We consider three classes of non-integrable frustrated spin models: (I) systems with local conserved quantities where the number of symmetry sectors grows exponentially with the system size but more slowly than the Hilbert space dimension, (II) systems with exact eigenstates that are singlet coverings, and (III) flat band systems hosting magnon crystals. We argue that several 1D and 2D models from class (I) exhibit disorder-free localization in high energy states so that information propagation is dynamically inhibited on length scales greater than a few lattice spacings. We further show that models of class (II) and (III) exhibit quantum many-body scars -- eigenstates of non-integrable Hamiltonians with finite energy density and anomalously low entanglement entropy. Our results demonstrate that magnetic frustration supplies a means to systematically construct classes of non-integrable models exhibiting anomalous thermalization in mid-spectrum states.


title: White Dwarfs in the Era of the LSST and its Synergies with Space-Based Missions

abstract: With the imminent start of the Legacy Survey for Space and Time (LSST) on the Vera C. Rubin Observatory, and several new space telescopes expected to begin operations later in this decade, both time domain and wide-field astronomy are on the threshold of a new era. In this paper, we use a new, multi-component model for the distribution of white dwarfs (WDs) in our Galaxy to simulate the WD populations in four upcoming wide-field surveys (i.e., LSST, Euclid, the Roman Space Telescope and CASTOR) and use the resulting samples to explore some representative WD science cases. Our results confirm that LSST will provide a wealth of information for Galactic WDs, detecting more than 150 million WDs at the final depth of its stacked, 10-year survey. Within this sample, nearly 300,000 objects will have 5$\sigma$ parallax measurements and nearly 7 million will have 5$\sigma$ proper motion measurements, allowing the detection of the turn-off in the halo WD luminosity function and the discovery of more than 200,000 ZZ Ceti stars. The wide wavelength coverage that will be possible by combining LSST data with observations from Euclid, and/or the Roman Space Telescope, will also discover more than 3,500 WDs with debris disks, highlighting the advantages of combining data between the ground- and space-based missions.


title: Dust size and spatial distributions in debris discs: predictions for exozodiacal dust dragged in from an exo-Kuiper belt

abstract: The SEDs of some nearby stars show mid-infrared excesses from warm habitable zone dust, known as exozodiacal dust. This dust may originate in collisions in a planetesimal belt before being dragged inwards. This paper presents an analytical model for the size distribution of particles at different radial locations in such a scenario, considering evolution due to destructive collisions and Poynting-Robertson (P-R) drag. Results from more accurate but computationally expensive numerical simulations of this process are used to validate the model and fit its free parameters. The model predicts 11 $\mu$m excesses ($R_{11}$) for discs with a range of dust masses and planetesimal belt radii using realistic grain properties. We show that P-R drag should produce exozodiacal dust levels detectable with the Large Binocular Telescope Interferometer (LBTI) ($R_{11} > 0.1\%$) in systems with known outer belts; non-detection may indicate dust depletion, e.g. by an intervening planet. We also find that LBTI could detect exozodiacal dust dragged in from a belt too faint to detect at far-infrared wavelengths, with fractional luminosity $f\sim 10^{-7}$ and radius $\sim 10-80$ au. Application to systems observed with LBTI shows that P-R drag can likely explain most (5/9) of the exozodiacal dust detections in systems with known outer belts; two systems ($\beta$ Uma and $\eta$ Corvi) with bright exozodi may be due to exocomets. We suggest that the three systems with exozodiacal dust detections but no known belt may have cold planetesimal belts too faint to be detectable in the far-infrared. Even systems without outer belt detections could have exozodiacal dust levels $R_{11} > 0.04\%$ which are problematic for exo-Earth imaging.


title: Clues on the history of early-type galaxies from SDSS spectra and GALEX photometry

abstract: Stellar population studies of early-type galaxies (ETGs) based on their optical stellar continuum suggest that these are quiescent systems. However, emission lines and ultraviolet photometry reveal a diverse population. We use a new version of the STARLIGHT spectral synthesis code and state-of-the-art stellar population models to simultaneously fit SDSS spectra and GALEX photometry for a sample of 3453 galaxies at $z < 0.1$ with $NUV-r > 5$ that are classified as elliptical by Galaxy Zoo. We reproduce $FUV$ magnitudes for 80 per cent of UV upturn galaxies selected using criteria from the literature, suggesting that additional stellar population ingredients such as binaries and extreme horizontal branch stars may have a limited contribution to the UV upturn. The addition of ultraviolet data leads to a broadening of the distributions of mean stellar ages, metallicities and attenuation. Stellar populations younger than $1\,$Gyr are required to reproduce the ultraviolet emission in 17 per cent of our sample. These systems represent 43 per cent of the sample at $5<NUV-r<5.5$ and span the same stellar mass range as other ETGs in our sample. ETGs with young stellar components have larger $H\alpha$ equivalent widths ($W_{H\alpha}$) and larger dust attenuation. Emission line ratios and $W_{H\alpha}$ indicate that the ionising source in these systems is a mixture of young and old stellar populations. Their young stellar populations are metal-poor, especially for high-mass galaxies, indicating recent star formation associated with rejuvenation events triggered by external processes, such as minor mergers.


title: Coarse graining spin foam quantum gravity -- a review

abstract: In quantum gravity, we envision renormalization as the key tool for bridging the gap between microscopic models and observable scales. For spin foam quantum gravity, which is defined on a discretisation akin to lattice gauge theories, the goal is to derive an effective theory on a coarser discretisation from the dynamics on the finer one, coarse graining the system in the process and thus relating physics at different scales. In this review I will discuss the motivation for studying renormalization in spin foam quantum gravity, e.g. to restore diffeomorphism symmetry, and explain how to define renormalization in a background independent setting by formulating it in terms of boundary data. I will motivate the importance of the boundary data by studying coarse graining of a concrete example and extending this to the spin foam setting. This will naturally lead me to the methods currently used for renormalizing spin foam quantum gravity, such as tensor network renormalization, and a discussion of recent results. I will conclude with an overview of future prospects and research directions.


title: Scaling relations of Fuzzy Dark Matter haloes I: individual systems in their cosmological environment

abstract: Dark matter models involving a very light bosonic particle, generally known as Fuzzy Dark Matter (FDM), have been recently attracting great interest in the cosmology community, as their wave-like phenomenology would simultaneously explain the longstanding mis-detection of a dark matter particle and help easing the small-scale issues related to the standard Cold Dark Matter (CDM) scenario. With the present work, we initiate a series of papers aiming at investigating the evolution of FDM structures in a cosmological framework performed with our N-body code AX-GADGET, detailing for the first time in the literature how the actual scaling relations between solitonic cores and host haloes properties are significantly affected by the dynamical state, morphology and merger history of the individual systems. In particular, in this first paper we confirm the ability of AX-GADGET to correctly reproduce the typical FDM solitonic core and we employ it to study the non-linear evolution of eight FDM haloes in their cosmological context through the zoom-in simulation approach. We find that the scaling relations identified in previous works for isolated systems are generally modified for haloes evolving in a realistic cosmological environment, and appear to be valid only as a limit for the most relaxed and spherically symmetric systems.


title: Local starburst conditions and formation of GRB 980425 / SN 1998bw within a collisional ring

abstract: We present the first spatially resolved study of molecular gas in the vicinity of a Gamma Ray Burst, using CO(2-1) emission line observations with the Atacama Large Millimetre Array (ALMA) at ~50 pc scales. The host galaxy of GRB 980425 contains a ring of high column density HI gas which is likely to have formed due to a collision between the GRB host and its companion galaxy, within which the GRB is located. We detect eleven molecular gas clumps in the galaxy, seven of which are within the gas ring. The clump closest to the GRB position is at a projected separation of ~280 pc. Although it is plausible that the GRB progenitor was ejected from clusters formed in this clump, we argue that the in situ formation of the GRB progenitor is the most likely scenario. We measure the molecular gas masses of the clumps and find them to be sufficient for forming massive star clusters. The molecular gas depletion times of the clumps show a variation of ~2 dex, comparable with the large variation in depletion times found in starburst galaxies in the nearby Universe. This demonstrates the presence of starburst modes of star formation on local scales in the galaxy, even while the galaxy as a whole cannot be categorised as a starburst based on its global properties. Our findings suggest that the progenitor of GRB 9802425 was originated in a young massive star cluster formed in the starburst mode of star formation.


title: Experimental realization of high-fidelity teleportation via non-Markovian open quantum system

abstract: Open quantum systems and study of decoherence are important for our fundamental understanding of quantum physical phenomena. For practical purposes, there exists a large number of quantum protocols exploiting quantum resources, e.g. entanglement, which allows to go beyond what is possible to achieve by classical means. We combine concepts from open quantum systems and quantum information science, and give a proof-of-principle experimental demonstration -- with teleportation -- that it is possible to implement efficiently a quantum protocol via non-Markovian open system. The results show that, at the time of implementation of the protocol, it is not necessary to have the quantum resource in the degree of freedom used for the basic protocol -- as long as there exists some other degree of freedom, or environment of an open system, which contains useful resources. The experiment is based on a pair of photons, where their polarizations act as open system qubits and frequencies as their environments -- while the path degree of freedom of one of the photons represents the state of Alice's qubit to be teleported to Bob's polarization qubit.


title: ALMA resolves molecular clouds in the metal poor Magellanic Bridge A

abstract: (Abridged)We characterize gas and dust emission in Magellanic Bridge A, which has the highest 870$\mu$m excess of emission found in single dish surveys. Using the ALMA telescope, we mapped the Magellanic Bridge A molecular cloud with sub-parsec resolution, in 1.3 mm continuum and CO(2-1) line emission. We also map the cloud in 870$\mu$m continuum and CO(2-1) line emission at ~6 pc resolution with APEX. We combine the ALMA and APEX CO(2-1) line cubes to study the molecular gas emission. Magellanic Bridge A breaks up into two distinct molecular clouds in dust and CO(2-1) emission, which we call North and South. Dust emission in the North source, according to our best parameters from fitting the far-infrarred fluxes, is ~3 K colder than in the South source in correspondence to its less developed star formation. Both dust sources present large submillimeter excesses in LABOCA data: according to our best fits the excess over the modified blackbody (MBB) fit to the Spitzer/Herschel continuum are ~7 and ~3 for the North and South sources respectively. Nonetheless, we do not detect the corresponding 1.3 mm continuum with ALMA. Our limits are compatible with the extrapolation of the MBB fits and therefore we cannot independently confirm the excess at this longer wavelength. The CO(2-1) emission is in two parsec-sized clouds with virial masses around 400 and 700 Mo each. Their volume densities are ~700-2600 cm$^{-3}$, larger than typical bulk densities of Galactic molecular clouds. The CO-to-H2 conversion factor is 6.5 and 15 M$_{\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ for the North and South clouds, calculated using their respective virial masses and CO(2-1) luminosities. Gas mass estimates from our MBB fits to dust emission yields masses $M\sim1.3\times10^3$ M$_{\odot}$ and $2.9\times10^3$ M$_{\odot}$ for North and South respectively, a factor ~4 larger than the virial masses we infer from CO.


title: Existence and Instability of Novel Hairy Black Holes in Shift-symmetric Horndeski Theories

abstract: Shift-symmetric Horndeski theories admit an interesting class of Schwarzschild black hole solutions exhibiting time-dependent scalar hair. By making use of Lema\^{i}tre coordinates, we analyze perturbations around these types of black holes, and demonstrate that scalar perturbations around black hole backgrounds inevitably have gradient instabilities. Taken together with previously established results, this newly-discovered instability rules out black holes with time-dependent scalar hair in Horndeski theories.


title: Optimal control of mean field equations with monotone coefficients and applications in neuroscience

abstract: We are interested in the optimal control problem associated with certain quadratic cost functionals depending on the solution $X=X^\alpha$ of the stochastic mean-field type evolution equation in $\mathbb R^d$ $dX_t=b(t,X_t,\mathcal L(X_t),\alpha_t)dt+\sigma(t,X_t,\mathcal L(X_t),\alpha_t)dW_t,$ $X_0\sim \mu$ given, under assumptions that enclose a sytem of FitzHugh-Nagumo neuron networks, and where for practical purposes the control $\alpha_t$ is deterministic. To do so, we assume that we are given a drift coefficient that satisfies a one-sided Lipshitz condition, and that the dynamics is subject to a (convex) level set constraint of the form $\pi(X_t)\leq0$. The mathematical treatment we propose follows the lines of the recent monograph of Carmona and Delarue for similar control problems with Lipshitz coefficients. After addressing the existence of minimizers via a martingale approach, we show a maximum principle and then numerically investigate a gradient algorithm for the approximation of the optimal control.


title: The VIMOS Ultra-Deep Survey: the Ly$\alpha$ emission line morphology at $2 < z < 6$

abstract: The Lyman-$\alpha$ (Ly$\alpha$) emission line has been ubiquitously used to confirm and study high redshift galaxies. We report on the line morphology as seen in the 2D spectra from the VIMOS Ultra Deep Survey in a sample of 914 Ly$\alpha$ emitters from a parent sample of 4192 star-forming galaxies at $2<z_\mathrm{spec}\lesssim6$. The study of the spatial extent of Ly$\alpha$ emission provides insight into the escape of Ly$\alpha$ photons from galaxies. We classify the line emission as either non-existent, coincident, projected spatial offset, or extended with respect to the observed 2D UV continuum emission. The line emitters in our sample are classified as ~45% coincident, ~24% extended and ~11% offset emitters. For galaxies with detected UV continuum, we show that extended Ly$\alpha$ emitters (LAEs) correspond to the highest equivalent width galaxies (with an average $W_\mathrm{Ly\alpha}\sim-22${\AA}). This means that this class of objects is the most common in narrow-band selected samples, which usually select high equivalent width LAEs, $<-20${\AA}. Extended Ly$\alpha$ emitters are found to be less massive, less star-forming, with lower dust content, and smaller UV continuum sizes ($r_{50}\sim0.9$kpc) of all the classes considered here. We also find that galaxies with larger UV-sizes have lower fractions of Ly$\alpha$ emitters. By stacking the spectra per emitter class we find that the weaker Ly$\alpha$ emitters have stronger low ionization inter-stellar medium (ISM) absorption lines. Interestingly, we find that galaxies with Ly$\alpha$ offset emission (median separation of $1.1_{-0.8}^{+1.3}$kpc from UV continuum) show similar velocity offsets in the ISM as those with no visible emission (and different from other Ly$\alpha$ emitting classes). This class of objects may hint at episodes of gas accretion, bright offset clumps, or on-going merging activity into the larger galaxies.


title: Dualities for three-dimensional $\mathcal{N} = 2$ $SU(N_c)$ chiral adjoint SQCD

abstract: We study dualities for 3d $\mathcal{N} = 2$ $SU(N_c)$ SQCD at Chern-Simons level $k$ in presence of an adjoint with polynomial superpotential. The dualities are dubbed chiral because there is a different amount of fundamentals $N_f$ and antifundamentals $N_a$. We build a complete classification of such dualities in terms of $ |N_f - N_a| $ and $k$. The classification is obtained by studying the flow from the non-chiral case, and we corroborate our proposals by matching the three-sphere partition functions. Finally, we revisit the case of $SU(N_c)$ SQCD without the adjoint, comparing our results with previous literature.


title: Quasi-probability distributions in Loop Quantum Cosmology

abstract: In this paper, we introduce a complete family of parametrized quasi-probability distributions in phase space and their corresponding Weyl quantization maps with the aim to generalize the recently developed Wigner-Weyl formalism within the Loop Quantum Cosmology program (LQC). In particular, we intend to define those quasi-distributions for states valued on the Bohr compactification of the real line in such a way that they are labeled by a parameter that accounts for the ordering ambiguity corresponding to non-commutative quantum operators. Hence, we notice that the projections of the parametrized quasi-probability distributions result in marginal probability densities which are invariant under any ordering prescription. We also note that, in opposition to the standard Schr\"odinger representation, for an arbitrary character the quasi-distributions determine a positive function independently of the ordering. Further, by judiciously implementing a parametric-ordered Weyl quantization map for LQG, we are able to recover in a simple manner the relevant cases of the standard, anti-standard, and Weyl symmetric orderings, respectively. We expect that our results may serve to analyze several fundamental aspects within the LQC program, in special those related to coherence, squeezed states, and the convergence of operators, as extensively analyzed in the quantum optics and in the quantum information frameworks.


title: Structure of Submetries

abstract: We investigate the geometric and topological structure of equidistant decompositions of Riemannian manifolds.


title: Minimally coupled scalar fields as imperfect fluids

abstract: We revisit the issue of the fluid description of minimally coupled scalar fields. While in a cosmological setup the interpretation of a time-evolving scalar field as a perfect fluid is well-understood, the situation is more intricate when the scalar field is static, but has a spatial gradient, a situation motivated by black hole perturbations in scalar-tensor theories. Then the scalar field is interpreted as either a particular imperfect fluid of type I or a superposition of a pair of leftgoing (incoming) and rightgoing (outgoing) null dusts with a perfect fluid. Finally, when the scalar gradient is null, it is equivalent to an imperfect fluid of type II, degenerating into null dust when the energy conditions are imposed. We also propose the suitable action in terms of the fluid pressure components for each case and discuss the variational principle for a generic class of minimally coupled scalar fields.


title: Debiasing Distributed Second Order Optimization with Surrogate Sketching and Scaled Regularization

abstract: In distributed second order optimization, a standard strategy is to average many local estimates, each of which is based on a small sketch or batch of the data. However, the local estimates on each machine are typically biased, relative to the full solution on all of the data, and this can limit the effectiveness of averaging. Here, we introduce a new technique for debiasing the local estimates, which leads to both theoretical and empirical improvements in the convergence rate of distributed second order methods. Our technique has two novel components: (1) modifying standard sketching techniques to obtain what we call a surrogate sketch; and (2) carefully scaling the global regularization parameter for local computations. Our surrogate sketches are based on determinantal point processes, a family of distributions for which the bias of an estimate of the inverse Hessian can be computed exactly. Based on this computation, we show that when the objective being minimized is $l_2$-regularized with parameter $\lambda$ and individual machines are each given a sketch of size $m$, then to eliminate the bias, local estimates should be computed using a shrunk regularization parameter given by $\lambda^{\prime}=\lambda\cdot(1-\frac{d_{\lambda}}{m})$, where $d_{\lambda}$ is the $\lambda$-effective dimension of the Hessian (or, for quadratic problems, the data matrix).


title: An efficient Kinetic Monte Carlo to study analyte capture by a nanopore: Transients, boundary conditions and time-dependent fields

abstract: To better understand the capture process by a nanopore, we introduce an efficient Kinetic Monte Carlo (KMC) algorithm that can simulate long times and large system sizes by mapping the dynamic of a point-like particle in a 3D spherically symmetric system onto the 1D biased random walk. Our algorithm recovers the steady-state analytical solution and allows us to study time-dependent processes such as transients. Simulation results show that the steady-state depletion zone near pore is barely larger than the pore radius and narrows at higher field intensities; as a result, the time to reach steady-state is much smaller than the time required to empty a zone of the size of the capture radius $\lambda_e$. When the sample reservoir has a finite size, a second depletion region propagates inward from the outer wall, and the capture rate starts decreasing when it reaches the capture radius $\lambda_e$. We also note that the flatness of the electric field near the pore, which is often neglected, induces a traffic jam that can increase the transient time by several orders of magnitude. Finally, we propose a new proof-of-concept scheme to separate two analytes of the same mobility but different diffusion coefficients using time-varying fields.


title: On the arithmetic of Pad\'e approximants to the exponential function

abstract: The $(u,v)$-Pad\'e approximation to a function $f$ is the (unique, up to scaling) rational approximation $f(x) = P(x)/Q(x) + O(x^{u+v+1})$, where $P$ has degree $u$ and $Q$ has degree $v$. Motivated by recent work of Molin, Pazuki, and Rabarison, we study the arithmetic of the Pad\'e approximants of the exponential polynomials. By viewing the approximants as certain Generalized Laguerre Polynomials, we determine the Galois groups of the diagonal approximants and prove some special cases of irreducibility.


title: A priori and a posteriori error estimates for the quad-curl eigenvalue problem

abstract: In this paper, we propose a new family of H(curl^2)-conforming elements for the quad-curl eigenvalue problem in 2D. The accuracy of this family is one order higher than that in [32]. We prove a priori and a posteriori error estimates. The a priori estimate of the eigenvalue with a convergence order 2(s-1) is obtained if the eigenvector u\in H^{s+1}(\Omega). For the a posteriori estimate, by analyzing the associated source problem, we obtain lower and upper bounds for the eigenvector in an energy norm and an upper bound for the eigenvalues. Numerical examples are presented for validation.


title: Detailed 3D structure of OrionA in dust with Gaia DR2

abstract: The unprecedented astrometry from Gaia DR2 provides us with an opportunity to study in detail molecular clouds in the solar neighbourhood. Extracting the wealth of information in these data remains a challenge, however. We have further improved our Gaussian Processes-based, three-dimensional dust mapping technique to allow us to study molecular clouds in more detail. These improvements include a significantly better scaling of the computational cost with the number of stars, and taking into account distance uncertainties to individual stars. Using Gaia DR2 astrometry together with 2MASS and WISE photometry for 30 000 stars, we infer the distribution of dust out to 600 pc in the direction of the Orion A molecular cloud. We identify a bubble-like structure in front of Orion A, centred at a distance of about 350 pc from the Sun. The main Orion A structure is visible at slightly larger distances, and we clearly see a tail extending over 100 pc that is curved and slightly inclined to the line-of-sight. The location of our foreground structure coincides with 5-10 Myr old stellar populations, suggesting a star formation episode that predates that of the Orion Nebula Cluster itself. We identify also the main structure of the Orion B molecular cloud, and in addition discover a background component to this at a distance of about 460 pc from the Sun. Finally, we associate our dust components at different distances with the plane-of-the-sky magnetic field orientation as mapped by Planck. This provides valuable information for modelling the magnetic field in 3D around star forming regions.


title: Meta-Learning Stationary Stochastic Process Prediction with Convolutional Neural Processes

abstract: Stationary stochastic processes (SPs) are a key component of many probabilistic models, such as those for off-the-grid spatio-temporal data. They enable the statistical symmetry of underlying physical phenomena to be leveraged, thereby aiding generalization. Prediction in such models can be viewed as a translation equivariant map from observed data sets to predictive SPs, emphasizing the intimate relationship between stationarity and equivariance. Building on this, we propose the Convolutional Neural Process (ConvNP), which endows Neural Processes (NPs) with translation equivariance and extends convolutional conditional NPs to allow for dependencies in the predictive distribution. The latter enables ConvNPs to be deployed in settings which require coherent samples, such as Thompson sampling or conditional image completion. Moreover, we propose a new maximum-likelihood objective to replace the standard ELBO objective in NPs, which conceptually simplifies the framework and empirically improves performance. We demonstrate the strong performance and generalization capabilities of ConvNPs on 1D regression, image completion, and various tasks with real-world spatio-temporal data.


title: A Renormalization Group Approach to Spontaneous Stochasticity

abstract: We develop a theoretical approach to ``spontaneous stochasticity'' in classical dynamical systems that are nearly singular and weakly perturbed by noise. This phenomenon is associated to a breakdown in uniqueness of solutions for fixed initial data and underlies many fundamental effects of turbulence (unpredictability, anomalous dissipation, enhanced mixing). Based upon analogy with statistical-mechanical critical points at zero temperature, we elaborate a renormalization group (RG) theory that determines the universal statistics obtained for sufficiently long times after the precise initial data are ``forgotten''. We apply our RG method to solve exactly the ``minimal model'' of spontaneous stochasticity given by a 1D singular ODE. Generalizing prior results for the infinite-Reynolds limit of our model, we obtain the RG fixed points that characterize the spontaneous statistics in the near-singular, weak-noise limit, determine the exact domain of attraction of each fixed point, and derive the universal approach to the fixed points as a singular large-deviations scaling, distinct from that obtained by the standard saddle-point approximation to stochastic path-integrals in the zero-noise limit. We present also numerical simulation results that verify our analytical predictions, propose possible experimental realizations of the ``minimal model'', and discuss more generally current empirical evidence for ubiquitous spontaneous stochasticity in Nature. Our RG method can be applied to more complex, realistic systems and some future applications are briefly outlined.

categories: cond-mat.stat-mech math-ph math.MP nlin.CD physics.ao-ph physics.flu-dyn


title: Multi-agent Planning for thermalling gliders using multi level graph-search

abstract: This paper solves a path planning problem for a group of gliders. The gliders are tasked with visiting a set of interest points. The gliders have limited range but are able to increase their range by visiting special points called thermals. The problem addressed in this paper is of path planning for the gliders such that, the total number of interest points visited by the gliders is maximized. This is referred to as the multi-agent problem. The problem is solved by first decomposing it into several single-agent problems. In a single-agent problem a set of interest points are allocated to a single glider. This problem is solved by planning a path which maximizes the number of visited interest points from the allocated set. This is achieved through a uniform cost graph search, as shown in our earlier work. The multi-agent problem now consists of determining the best allocation (of interest points) for each glider. Two ways are presented of solving this problem, a brute force search approach as shown in earlier work and a Branch\&Bound type graph search. The Branch&Bound approach is the main contribution of the paper. This approach is proven to be optimal and shown to be faster than the brute force search using simulations.


title: Clustering of Electromagnetic Showers and Particle Interactions with Graph Neural Networks in Liquid Argon Time Projection Chambers Data

abstract: Liquid Argon Time Projection Chambers (LArTPCs) are a class of detectors that produce high resolution images of charged particles within their sensitive volume. In these images, the clustering of distinct particles into superstructures is of central importance to the current and future neutrino physics program. Electromagnetic (EM) activity typically exhibits spatially detached fragments of varying morphology and orientation that are challenging to efficiently assemble using traditional algorithms. Similarly, particles that are spatially removed from each other in the detector may originate from a common interaction. Graph Neural Networks (GNNs) were developed in recent years to find correlations between objects embedded in an arbitrary space. The Graph Particle Aggregator (GrapPA) first leverages GNNs to predict the adjacency matrix of EM shower fragments and to identify the origin of showers, i.e. primary fragments. On the PILArNet public LArTPC simulation dataset, the algorithm achieves achieves a shower clustering accuracy characterized by a mean adjusted Rand index (ARI) of 97.8 % and a primary identification accuracy of 99.8 %. It yields a relative shower energy resolution of $(4.1+1.4/\sqrt{E (\text{GeV})})\,\%$ and a shower direction resolution of $(2.1/\sqrt{E(\text{GeV})})^{\circ}$. The optimized algorithm is then applied to the related task of clustering particle instances into interactions and yields a mean ARI of 99.2 % for an interaction density of $\sim\mathcal{O}(1)\,m^{-3}$.


title: The unbounded denominator conjecture for the noncongruence subgroups of index $7$

abstract: We study modular forms for the minimal index noncongruence subgroups of the modular group. Our main theorem is a proof of the unbounded denominator conjecture for these groups, and we also provide a study of the Fourier coefficients of Eisenstein series for one of these minimal groups.


title: Metric Dimension of Hamming Graphs and Applications to Computational Biology

abstract: Genetic sequencing has become an increasingly affordable and accessible source of genomic data in computational biology. This data is often represented as $k$-mers, i.e., strings of some fixed length $k$ with symbols chosen from a reference alphabet. In contrast, some of the most effective and well-studied machine learning algorithms require numerical representations of the data. The concept of metric dimension of the so-called Hamming graphs presents a promising way to address this issue. A subset of vertices in a graph is said to be resolving when the distances to those vertices uniquely characterize every vertex in the graph. The metric dimension of a graph is the size of a smallest resolving subset of vertices. Finding the metric dimension of a general graph is a challenging problem, NP-complete in fact. Recently, an efficient algorithm for finding resolving sets in Hamming graphs has been proposed, which suffices to uniquely embed $k$-mers into a real vector space. Since the dimension of the embedding is the cardinality of the associated resolving set, determining whether or not a node can be removed from a resolving set while keeping it resolving is of great interest. This can be quite challenging for large graphs since only a brute-force approach is known for checking whether a set is a resolving set or not. In this thesis, we characterize resolvability of Hamming graphs in terms of a linear system over a finite domain: a set of nodes is resolving if and only if the linear system has only a trivial solution over said domain. We can represent the domain as the roots of a polynomial system so the apparatus of Gr\"obner bases comes in handy to determine, whether or not a set of nodes is resolving. As proof of concept, we study the resolvability of Hamming graphs associated with octapeptides i.e. proteins sequences of length eight.


title: Algebraic curves with automorphism groups of large prime order

abstract: Let $\mathcal{X}$ be an algebraic curve of genus $g$ defined over an algebraically closed field $K$ of characteristic $p \geq 0$, and $q$ a prime dividing $|\mbox{Aut}(\mathcal{X})|$. We say that $\mathcal{X}$ is a $q$-curve. Homma proved that either $q \leq g+1$ or $q = 2g+1$, and classified $(2g+1)$-curves. In this note, we classify $(g+1)$-curves, and fully characterize the automorphism groups of $q$-curves for $q= 2g+1, g+1$. We also give some partial results on $q$-curves for $q = g, g-1$.


title: Improved Cauchy-characteristic evolution system for high-precision numerical relativity waveforms

abstract: We present several improvements to the Cauchy-characteristic evolution procedure that generates high-fidelity gravitational waveforms at $\mathcal{I}^+$ from numerical relativity simulations. Cauchy-characteristic evolution combines an interior solution of the Einstein field equations based on Cauchy slices with an exterior solution based on null slices that extend to $\mathcal{I}^+$. The foundation of our improved algorithm is a comprehensive method of handling the gauge transformations between the arbitrarily specified coordinates of the interior Cauchy evolution and the unique (up to BMS transformations) Bondi-Sachs coordinate system of the exterior characteristic evolution. We present a reformulated set of characteristic evolution equations better adapted to numerical implementation. In addition, we develop a method to ensure that the angular coordinates used in the volume during the characteristic evolution are asymptotically inertial. This provides a direct route to an expanded set of waveform outputs and is guaranteed to avoid pure-gauge logarithmic dependence that has caused trouble for previous spectral implementations of the characteristic evolution equations. We construct a set of Weyl scalars compatible with the Bondi-like coordinate systems used in characteristic evolution, and determine simple, easily implemented forms for the asymptotic Weyl scalars in our suggested set of coordinates.


title: Lack of evidence for a substantial rate of templated mutagenesis in B cell diversification

abstract: B cell receptor sequences diversify through mutations introduced by purpose-built cellular machinery. A recent paper has concluded that a "templated mutagenesis" process is a major contributor to somatic hypermutation, and therefore immunoglobulin diversification, in mice and humans. In this proposed process, mutations in the immunoglobulin locus are introduced by copying short segments from other immunoglobulin genes. If true, this would overturn decades of research on B cell diversification, and would require a complete re-write of computational methods to analyze B cell data for these species. In this paper, we re-evaluate the templated mutagenesis hypothesis. By applying the original inferential method using potential donor templates absent from B cell genomes, we obtain estimates of the methods's false positive rates. We find false positive rates of templated mutagenesis in murine and human immunoglobulin loci that are similar to or even higher than the original rate inferences, and by considering the bases used in substitution we find evidence that if templated mutagenesis occurs, it is at a low rate. We also show that the statistically significant results in the original paper can easily result from a slight misspecification of the null model.


title: Ideal Free Dispersal under General Spatial heterogeneity and Time Periodicity

abstract: A population is said to have an ideal free distribution in a spatially heterogeneous but temporally constant environment if each of its members have chosen a fixed spatial location in a way that optimizes its individual fitness, allowing for the effects of crowding. In this paper, we extend the idea of individual fitness associated with a specific location in space to account for the full path that an individual organism takes in space and time over a periodic cycle, and extend the mathematical formulation of an ideal free distribution to general time periodic environments. We find that, as in many other cases, populations using dispersal strategies that can produce a generalized ideal free distribution have a competitive advantage relative to populations using strategies that do not produce an ideal free distribution. A sharp criterion on the environmental functions is found to be necessary and sufficient for such ideal free distribution to be feasible. In the case the criterion is met, we showed that there exist dispersal strategies that can be identified as producing a time-periodic version of an ideal free distribution, and such strategies are evolutionarily steady and are neighborhood invaders from the viewpoint of adaptive dynamics. Our results extend previous works in which the environments are either temporally constant, or temporally periodic but the total carrying capacity is temporally constant.


title: X-ray Properties of Young Radio Quasars at z > 4.5

abstract: We present a comprehensive analysis of Chandra X-ray observations of 15 young radio quasars at redshifts $4.5 < z < 5.0$. All sources are detected in the $0.5-7.0$ keV energy band. Emission spectra are extracted, and the average photon index for the sample is measured to be $1.5\pm0.1$. Unabsorbed rest-frame $2-10$ keV luminosities are found to range between $(0.5-23.2) \times 10^{45}$ erg s$^{-1}$. The optical-X-ray power-law spectral index $\alpha_{ox}$ is calculated for each source using optical/UV data available in the literature. The $\alpha_{ox}$-UV relationship is compared with other quasar surveys, and an anticorrelation is observed that agrees with independent estimates. Rest-frame radio and X-ray luminosities are established for the sample, and a correlation between the luminosities is detected. These multiwavelength results reinforce a lack of spectral evolution for quasars over a broad redshift range. We additionally identify three quasars from our multiwavelength analysis that are statistically significant outliers, with one source being a Compton-thick candidate in the early universe, and discuss each in detail.


title: Molecular theory of Langevin dynamics for active self-diffusiophoretic colloids

abstract: Active colloidal particles that are propelled by a self-diffusiophoretic mechanism are often described by Langevin equations that are either postulated on physical grounds or derived using the methods of fluctuating hydrodynamics. While these descriptions are appropriate for colloids of micrometric and larger size, they will break down for very small active particles. A fully microscopic derivation of Langevin equations for self-diffusiophoretic particles powered by chemical reactions catalyzed asymmetrically by the colloid is given in this paper. The derivation provides microscopic expressions for the translational and rotational friction tensors, as well as reaction rate coefficients appearing in the Langevin equations. The diffusiophoretic force and torque are expressed in terms of nonequilibrium averages of fluid fields that satisfy generalized transport equations. The results provide a description of active motion on small scales where descriptions in terms of coarse grained continuum fluid equations combined with boundary conditions that account for the presence of the colloid may not be appropriate.


title: Decoding asymptomatic COVID-19 infection and transmission

abstract: Coronavirus disease 2019 (COVID-19) is a continuously devastating public health and the world economy. One of the major challenges in controlling the COVID-19 outbreak is its asymptomatic infection and transmission, which are elusive and defenseless in most situations. The pathogenicity and virulence of asymptomatic COVID-19 remain mysterious. Based on the genotyping of 20656 Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genome isolates, we reveal that asymptomatic infection is linked to SARS-CoV-2 11083G>T mutation, i.e., leucine (L) to phenylalanine (F) substitution at the residue 37 (L37F) of nonstructure protein 6 (NSP6). By analyzing the distribution of 11083G>T in various countries, we unveil that 11083G>T may correlate with the hypotoxicity of SARS-CoV-2. Moreover, we show a global decaying tendency of the 11083G>T mutation ratio indicating that 11083G>T hinders SARS-CoV-2 transmission capacity. Sequence alignment found both NSP6 and residue 37 neighborhoods are relatively conservative over a few coronaviral species, indicating their importance in regulating host cell autophagy to undermine innate cellular defense against viral infection. Using machine learning and topological data analysis, we demonstrate that mutation L37F has made NSP6 energetically less stable. The rigidity and flexibility index and several network models suggest that mutation L37F may have compromised the NSP6 function, leading to a relatively weak SARS-CoV subtype. This assessment is a good agreement with our genotyping of SARS-CoV-2 evolution and transmission across various countries and regions over the past few months.


title: Convexity of the weighted Mabuchi functional and the uniqueness of weighted extremal metrics

abstract: We prove the uniqueness, up to a pull-back by an element of a suitable subgroup of complex automorphisms, of the weighted extremal K\"ahler metrics on a compact K\"ahler manifold introduced in our previous work. This extends a result by Berman--Berndtsson and Chen--Paun--Zeng in the extremal K\"ahler case. Furthermore, we show that a weighted extremal K\"ahler metric is a global minimum of a suitable weighted version of the modified Mabuchi energy. This implies a suitable notion of weighted K-semistability of a K\"ahler manifold admitting a weighted extremal K\"ahler metric.


title: Spectral Methods for Ranking with Scarce Data

abstract: Given a number of pairwise preferences of items, a common task is to rank all the items. Examples include pairwise movie ratings, New Yorker cartoon caption contests, and many other consumer preferences tasks. What these settings have in common is two-fold: a scarcity of data (it may be costly to get comparisons for all the pairs of items) and additional feature information about the items (e.g., movie genre, director, and cast). In this paper we modify a popular and well studied method, RankCentrality for rank aggregation to account for few comparisons and that incorporates additional feature information. This method returns meaningful rankings even under scarce comparisons. Using diffusion based methods, we incorporate feature information that outperforms state-of-the-art methods in practice. We also provide improved sample complexity for RankCentrality in a variety of sampling schemes.


title: Distinguishing Black Hole and Naked Singularity in MOG via Inertial Frame Dragging Effect

abstract: We analyze the generalized spin precession of a test gyroscope around a stationary spacetime i.e. for Kerr-MOG black hole~(BH) in scalar-tensor-vector gravity or modified gravity~(MOG). A detailed study of generalized spin frequency has been done for \emph{non} extremal Kerr-MOG BH, \emph{extremal} Kerr-MOG BH and \emph{naked singularity~(NS)} in comparison to non-extremal BH, extremal BH, and NS of Kerr spacetime. The generalized spin frequency that {we have} computed could be expressed in terms of {the} BH mass parameter, the angular momentum parameter, and the MOG parameter. Moreover, we differentiate the non-extremal BH, extremal BH, and NS via computation of the said precession frequency. The Lense-Thirring~(LT) frequency {can} obtain from generalized spin frequency by taking the limit as $\Omega=0$ i. e. {when the} angular frequency is set to zero limit. Furthermore, we compute the LT frequency for various {values of} angular coordinates i.e. starting from polar to {the} equatorial plane. We show that the LT frequency diverges at the horizon for extremal BH. Finally, we study the accretion disk physics by computing three epicyclic frequencies namely the Keplerian frequency, {the} radial epicyclic frequency and {the} vertical epicyclic frequency. We also compute the periastron frequency and nodal frequency. With the aid of these frequency profiles, {one} can distinguish three compact objects i. e. \emph{non-extremal BH, extremal BH} {and} \emph{NS}.


title: Efficient Neural Network Deployment for Microcontroller

abstract: Edge computing for neural networks is getting important especially for low power applications and offline devices. TensorFlow Lite and PyTorch Mobile were released for this purpose. But they mainly support mobile devices instead of microcontroller level yet. Microcontroller support is an emerging area now. There are many approaches to reduce network size and compute load like pruning, binarization and layer manipulation i.e. operator reordering. This paper is going to explore and generalize convolution neural network deployment for microcontrollers with two novel optimization proposals offering memory saving and compute efficiency in 2D convolutions as well as fully connected layers. The first one is in-place max-pooling, if the stride is greater than or equal to pooling kernel size. The second optimization is to use ping-pong buffers between layers to reduce memory consumption significantly. The memory savings and performance will be compared with CMSIS-NN framework developed for ARM Cortex-M CPUs. The final purpose is to develop a tool consuming PyTorch model with trained network weights, and it turns into an optimized inference engine(forward pass) in C/C++ for low memory(kilobyte level) and limited computing capable microcontrollers.


title: Networks with Growth and Preferential Attachment: Modeling and Applications

abstract: In this article we presented a brief study of the main network models with growth and preferential attachment. Such models are interesting because they present several characteristics of real systems. We started with the classical model proposed by Barabasi and Albert: nodes are added to the network connecting preferably to other nodes that are more connected. We also presented models that consider more representative elements from social perspectives, such as the homophily between the vertices or the fitness that each node has to build connections. Furthermore, we showed a version of these models including the Euclidean distance between the nodes as a preferential attachment rule. Our objective is to investigate the basic properties of these networks as distribution of connectivity, degree correlation, shortest path, cluster coefficient and how these characteristics are affected by the preferential attachment rules. Finally, we also provided a comparison of these synthetic networks with real ones. We found that characteristics as homophily, fitness and geographic distance are significant preferential attachment rules to modeling real networks. These rules can change the degree distribution form of these synthetic network models and make them more suitable to model real networks.


title: Uncertainty Prediction for Deep Sequential Regression Using Meta Models

abstract: Generating high quality uncertainty estimates for sequential regression, particularly deep recurrent networks, remains a challenging and open problem. Existing approaches often make restrictive assumptions (such as stationarity) yet still perform poorly in practice, particularly in presence of real world non-stationary signals and drift. This paper describes a flexible method that can generate symmetric and asymmetric uncertainty estimates, makes no assumptions about stationarity, and outperforms competitive baselines on both drift and non drift scenarios. This work helps make sequential regression more effective and practical for use in real-world applications, and is a powerful new addition to the modeling toolbox for sequential uncertainty quantification in general.


title: Polychromatic soliton molecules

abstract: Soliton molecules, bound states composed of interacting fundamental solitons, exhibit remarkable resemblance with chemical compounds and phenomena in quantum mechanics. Whereas optical molecules composed of two or more temporally locked solitons have been observed in a variety of platforms, soliton molecules formed by bound solitons at different frequencies have only recently been theoretically proposed. Here, we report the observation of polychromatic soliton molecules within a mode-locked laser cavity, achieving the desired dispersion by implementing a spectral pulse-shaper. This system supports two or more coincident solitons with different frequencies, but common group velocities. Our results open new directions of exploration in nonlinear dynamics within systems with complex dispersion and offer a convenient platform for optical analogies to mutual trapping and spectral tunneling in quantum mechanics.


title: Algebraic invariants of orbit configuration spaces in genus zero associated to finite groups

abstract: We consider orbit configuration spaces associated to finite groups acting freely by orientation preserving homeomorphisms on the $2$-sphere minus a finite number of points. Such action is equivalent to a homography action of a finite subgroup $G\subset \mathrm{PGL}(\mathbb{C}^2)$ on the complex projective line $\mathbb{P}^1$ minus a finite set $Z$ stable under $G$. We compute the cohomology ring and the Poincar\'e series of the orbit configuration space $C_n^G(\mathbb{P}^1 \setminus Z)$. This can be seen as a generalization of the work of Arnold for the classical configuration space $C_n(\mathbb{C})$ ($(G,Z)=(\{1\},\infty$)). It follows from the work that $C_n^G(\mathbb{P}^1\setminus Z)$ is formal in the sense of rational homotopy theory. We also prove the existence of an LCS formula relating the Poincar\'e series of $C_n^G(\mathbb{P}^1\setminus Z)$ to the ranks of quotients of successive terms of the lower central series of the fundamental group of $C_n^G(\mathbb{P}^1 \setminus Z)$. The successive quotients correspond to homogenous elements of graded Lie algebras introduced by the author in an earlier work. Such formula is also known for classical configuration spaces of $\mathbb{C}$, where fundamental groups are Artin braid groups and the ranks correspond to dimensions of homogenous elements of the Kohno-Drinfeld Lie algebras.


title: Spatially offset black holes in the Horizon-AGN simulation and comparison to observations

abstract: We study the displacements between the centres of galaxies and their supermassive black holes (BHs) in the cosmological hydrodynamical simulation Horizon-AGN, and in a variety of observations from the literature. The BHs in Horizon-AGN feel a sub-grid dynamical friction force, sourced by the surrounding gas, which prevents recoiling BHs being ejected from the galaxy. We find that i) the fraction of spatially offset BHs increases with cosmic time, ii) BHs live on prograde orbits in the plane of the galaxy with an orbital radius that decays with time but stalls near $z=0$, and iii) the magnitudes of offsets from the galaxy centres are substantially larger in the simulation than in observations. We attribute the stalling of the infall and excessive offset magnitudes to the fact that dynamical friction from stars and dark matter is not modelled in the simulation, and hence provide a way to improve the black hole dynamics of future simulations.


title: On the involution fixity of simple groups

abstract: Let $G$ be a finite permutation group of degree $n$ and let ${\rm ifix}(G)$ be the involution fixity of $G$, which is the maximum number of fixed points of an involution. In this paper we study the involution fixity of almost simple primitive groups whose socle $T$ is an alternating or sporadic group; our main result classifies the groups of this form with ${\rm ifix}(T) \leqslant n^{4/9}$. This builds on earlier work of Burness and Thomas, who studied the case where $T$ is an exceptional group of Lie type, and it strengthens the bound ${\rm ifix}(T) > n^{1/6}$ (with prescribed exceptions), which was proved by Liebeck and Shalev in 2015. A similar result for classical groups will be established in a sequel.


title: Impact of evidence-based flipped or active-engagement non-flipped courses on student performance in introductory physics

abstract: We describe the impact of physics education research-based pedagogical techniques in flipped and active-engagement non-flipped courses on student performance on validated conceptual surveys. We compare student performance in courses that make significant use of evidence-based active engagement (EBAE) strategies with courses that primarily use lecture-based (LB) instruction. All courses had large enrollment and often had 100-200 students. The analysis of data for validated conceptual surveys presented here includes data from large numbers of students from two-semester sequences of introductory algebra-based and calculus-based introductory physics courses. The conceptul surveys used to assess student learning in the first and second semester courses were the Force Concept Inventory and the Conceptual Survey of Electricity and Magnetism, respectively. In the research discussed here, the performance of students in EBAE courses at a particular level is compared with LB courses in two situations: (i) the same instructor taught two courses, one of which was a flipped course involving EBAE methods and the other an LB course, while the homework, recitations, and final exams were kept the same; (ii) student perforamnce in all of the EBAE courses taught by different instructors was averaged and compared with LB courses of the same type also averaged over different instructors. In all classes, we find that students in courses that make significant use of active-engagement strategies, on average, outperformed students in courses using primarily LB instruction of the same type on conceptual surveys even though there was no statistically significant difference on the pretest before instruction. We also discuss correlation between the performance on the validated conceptual surveys and the final exam, which typically placed a heavy weight on quantitative problem solving.


title: Decoder-free Robustness Disentanglement without (Additional) Supervision

abstract: Adversarial Training (AT) is proposed to alleviate the adversarial vulnerability of machine learning models by extracting only robust features from the input, which, however, inevitably leads to severe accuracy reduction as it discards the non-robust yet useful features. This motivates us to preserve both robust and non-robust features and separate them with disentangled representation learning. Our proposed Adversarial Asymmetric Training (AAT) algorithm can reliably disentangle robust and non-robust representations without additional supervision on robustness. Empirical results show our method does not only successfully preserve accuracy by combining two representations, but also achieve much better disentanglement than previous work.


title: Efficient computation and analysis of distributional Shapley values

abstract: Distributional data Shapley value (DShapley) has recently been proposed as a principled framework to quantify the contribution of individual datum in machine learning. DShapley develops the foundational game theory concept of Shapley values into a statistical framework and can be applied to identify data points that are useful (or harmful) to a learning algorithm. Estimating DShapley is computationally expensive, however, and this can be a major challenge to using it in practice. Moreover, there has been little mathematical analyses of how this value depends on data characteristics. In this paper, we derive the first analytic expressions for DShapley for the canonical problems of linear regression, binary classification, and non-parametric density estimation. These analytic forms provide new algorithms to estimate DShapley that are several orders of magnitude faster than previous state-of-the-art methods. Furthermore, our formulas are directly interpretable and provide quantitative insights into how the value varies for different types of data. We demonstrate the practical efficacy of our approach on multiple real and synthetic datasets.


title: Amplitudes of minima in dynamic conductance spectra of the SNS Andreev contact

abstract: Despite several theoretical approaches describing multiple Andreev reflections (MAR) effect in superconductor-normal metal-superconductor (SNS) junction are elaborated, the problem of comprehensive and adequate description of MAR is highly actual. In particular, a broadening parameter $\Gamma$ is still unaccounted at all, whereas a ballistic condition (the mean free path for inelastic scattering $l$ to the barrier width $d$ ratio) is considered only in the framework of K\"{u}mmel, Gunsenheimer, and Nikolsky (KGN), as well as Gunsenheimer-Zaikin approaches, for an isotropic case and fully-transparent constriction. Nonetheless, an influence of $l/d$ ratio to the dynamic conductance spectrum ($dI/dV$) features remains disregarded, thus being one of the aims of the current work. Our numerical calculations in the framework of an extended KGN approach develop the $l/d$ variation to determine both the number of the Andreev features and their amplitudes in the $dI/dV$ spectrum. We show, in the spectrum of a diffusive SNS junction ($l/d \rightarrow 1$) a suppression of the Andreev excess current, dramatic change in the current voltage $I(V)$-curve slope at low bias, with only the main harmonic at $eV=2\Delta$ bias voltage remains well-distinguished in the $dI/dV$-spectrum. Additionally, we attempt to make a first-ever comparison between experimental data for the high-transparency SNS junctions (more than $ 85~\%$) and theoretical predictions. As a result, we calculate the temperature dependences of amplitudes and areas of Andreev features within the extended KGN approach, which qualitatively agrees with our experimental data obtained using a "break-junction" technique.


title: A Bayesian Multilingual Document Model for Zero-shot Topic Identification and Discovery

abstract: In this paper, we present a Bayesian multilingual document model for learning language-independent document embeddings. The model is an extension of BaySMM [Kesiraju et al 2020] to the multilingual scenario. It learns to represent the document embeddings in the form of Gaussian distributions, thereby encoding the uncertainty in its covariance. We propagate the learned uncertainties through linear classifiers that benefit zero-shot cross-lingual topic identification. Our experiments on 17 languages show that the proposed multilingual Bayesian document model performs competitively, when compared to other systems based on large-scale neural networks (LASER, XLM-R, mUSE) on 8 high-resource languages, and outperforms these systems on 9 mid-resource languages. We revisit cross-lingual topic identification in zero-shot settings by taking a deeper dive into current datasets, baseline systems and the languages covered. We identify shortcomings in the existing evaluation protocol (MLDoc dataset), and propose a robust alternative scheme, while also extending the cross-lingual experimental setup to 17 languages. Finally, we consolidate the observations from all our experiments, and discuss points that can potentially benefit the future research works in applications relying on cross-lingual transfers.


title: A New ECDF Two-Sample Test Statistic

abstract: Empirical cumulative distribution functions (ECDFs) have been used to test the hypothesis that two samples come from the same distribution since the seminal contribution by Kolmogorov and Smirnov. This paper describes a statistic which is usable under the same conditions as Kolmogorov-Smirnov, but provides more power than other extant tests in that vein. I demonstrate a valid (conservative) procedure for producing finite-sample p-values. I outline the close relationship between this statistic and its two main predecessors. I also provide a public R package (CRAN: twosamples [2018]) implementing the testing procedure in $O(N\log(N))$ time with $O(N)$ memory. Using the package's functions, I perform several simulation studies showing the power improvements.


title: Bound States in Continuum and Zero-Index Metamaterials: A Review

abstract: Bound states in the continuum (BICs) are waves that remain localized even though they coexist with a continuous spectrum of radiating waves that can carry energy away. Their very existence defies conventional wisdom. Although BICs were first proposed in quantum mechanics, they are a general wave phenomenon and have since been identified in electromagnetic waves, acoustic waves in air, water waves and elastic waves in solids. These states have been studied in a wide range of material systems, such as piezoelectric materials, dielectric photonic crystals, optical waveguides and fibers, quantum dots, graphene and topological insulators. In this Review, we describe recent developments in this field with an emphasis on the physical mechanisms that lead to BICs across seemingly very different materials and types of waves. We also discuss experimental realizations, existing applications and directions for future work. At last, we present our recent effort to design a novel type of silicon (Si) based mu-Near-Zero (uNZ) metamaterials that reduces radiative loss through destructive interference of multiple loss channels, resulting in a bound state in the continuum. This design consists of Si pillars array. By adjusting the unit-cell dimensions including radius, pitch and height, we can eliminate the out-of-plane radiation.


title: Thermodynamics and the Evolution of Stochastic Populations

abstract: The appeal of thermodynamics to problems outside physics is undeniable, as is the growing recognition of its apparent universality, yet in the absence of a rigorous formalism divorced from the peculiarities of molecular systems all attempts to generalize thermodynamics remain qualitative and heuristic at best. In this paper we formulate a probabilistic theory of thermodynamics and and set the basis for its application to generic stochastic populations.


title: Optimal evolutionary decision-making to store immune memory

abstract: The adaptive immune system provides a diverse set of molecules that can mount specific responses against a multitude of pathogens. Memory is a key feature of adaptive immunity, which allows organisms to respond more readily upon re-infections. However, differentiation of memory cells is still one of the least understood cell fate decisions. Here, we introduce a mathematical framework to characterize optimal strategies to store memory to maximize the utility of immune response over an organism's lifetime. We show that memory production should be actively regulated to balance between affinity and cross-reactivity of immune receptors for an effective protection against evolving pathogens. Moreover, we predict that specificity of memory should depend on the organism's lifespan, and shorter-lived organisms with fewer pathogenic encounters should store more cross-reactive memory. Our framework provides a baseline to gauge the efficacy of immune memory in light of an organism's coevolutionary history with pathogens.


title: Generic and Mod p Kazhdan-Lusztig Theory for GL_2

abstract: Let $F$ be a non-archimedean local field with residue field $\mathbb{F}_q$ and let $G = GL_{2/F}$. Let $\mathbf{q}$ be an indeterminate and let $H^{(1)}(\mathbf{q})$ be the generic pro-p Iwahori-Hecke algebra of the group $G(F)$. Let $V_{\widehat{G}}$ be the Vinberg monoid of the dual group $\widehat{G}$. We establish a generic version for $H^{(1)}(\mathbf{q})$ of the Kazhdan-Lusztig-Ginzburg spherical representation, the Bernstein map and the Satake isomorphism. We define the flag variety for the monoid $V_{\widehat{G}}$ and establish the characteristic map in its equivariant K-theory. These generic constructions recover the classical ones after the specialization $\mathbf{q} = q \in \mathbb{C}$. At $\mathbf{q} = q = 0 \in\overline{\mathbb{F}}_q$, the spherical map provides a dual parametrization of all the irreducible $H^{(1)}_{\overline{\mathbb{F}}_q}(0)$-modules.


title: Pulsating chromosphere of classical Cepheids. Calcium infrared triplet and H$\alpha$ profile variations

abstract: It has been shown recently that the infrared emission of Cepheids, which is constant over the pulsation cycle, might be due to a pulsating shell of ionized gas of about 15\% of the stellar radius, which could be attributed to the chromospheric activity of Cepheids. The aim of this paper is to investigate the dynamical structure of the chromosphere of Cepheids along the pulsation cycle and quantify its size. We present H$\alpha$ and Calcium Near InfraRed triplet (Ca IR) profile variations using high-resolution spectroscopy with the UVES spectrograph of a sample of 24 Cepheids with a good period coverage from $\approx$ 3 to 60 days. After a qualitative analysis of the spectral lines profiles, we quantify the Van Hoof effect (velocity gradient between the H$\alpha$ and Ca IR) as a function of the period of the Cepheids. Then, we use the Schwarzschild mechanism (a line doubling due to a shock wave) to quantify the size of the chromosphere. We find a significant Van Hoof effect for Cepheids with period larger than $P=10$ days, in particular H$\alpha$ lines are delayed with a velocity gradient up to $\Delta v \approx$30 km/s compared to Ca IR. We find that the size of the chromosphere of long-period Cepheids is of at least $\approx$ 50\% of the stellar radius, which is consistent at first order with the size of the shell made of ionized gas previously found from the analysis of infrared excess. Last, for most of the long-period Cepheids in the sample, we report a motionless absorption feature in the H$\alpha$ line that we attribute to a circumstellar envelope that surrounds the chromosphere. Analyzing the Ca~IR lines of Cepheids is of importance to potentially unbias the period-luminosity relation from their infrared excess, particularly in the context of forthcoming observations from the Radial Velocity Spectrometer (RVS) on board \textit{Gaia}, that could be sensitive to their chromosphere.


title: Modular categories with transitive Galois actions

abstract: In this paper, we study modular categories whose Galois group actions on their simple objects are transitive. We show that such modular categories admit unique factorization into prime transitive factors. The representations of $SL_2(\mathbb{Z})$ associated with transitive modular categories are proven to be minimal and irreducible. Together with the Verlinde formula, we characterize prime transitive modular categories as the Galois conjugates of the adjoint subcategory of the quantum group modular category $\mathcal{C}(\mathfrak{sl}_2,p-2)$ for some prime $p > 3$. As a consequence, we completely classify transitive modular categories. Transitivity of super-modular categories can be similarly defined. A unique factorization of any transitive super-modular category into s-simple transitive factors is obtained, and the split transitive super-modular categories are completely classified.


title: Lecture Notes on Control System Theory and Design

abstract: This is a collection of the lecture notes of the three authors for a first-year graduate course on control system theory and design (ECE 515 , formerly ECE 415) at the ECE Department of the University of Illinois at Urbana-Champaign. This is a fundamental course on the modern theory of dynamical systems and their control, and builds on a first-level course in control that emphasizes frequency-domain methods (such as the course ECE 486 , formerly ECE 386, at UIUC ). The emphasis in this graduate course is on state space techniques, and it encompasses modeling , analysis (of structural properties of systems, such as stability, controllability, and observability), synthesis (of observers/compensators and controllers) subject to design specifications, and optimization . Accordingly, this set of lecture notes is organized in four parts, with each part dealing with one of the issues identified above. Concentration is on linear systems , with nonlinear systems covered only in some specific contexts, such as stability and dynamic optimization. Both continuous-time and discrete-time systems are covered, with the former, however, in much greater depth than the latter. The main objective of this course is to teach the student some fundamental principles within a solid conceptual framework, that will enable her/him to design feedback loops compatible with the information available on the "states" of the system to be controlled, and by taking into account considerations such as stability, performance, energy conservation, and even robustness. A second objective is to familiarize her/him with the available modern computational, simulation, and general software tools that facilitate the design of effective feedback loops


title: Recovery of Returning Halley-Type Comet 12P/Pons-Brooks With the Lowell Discovery Telescope

abstract: We report the recovery of returning Halley-type comet 12P/Pons-Brooks using the 4.3 m Lowell Discovery Telescope, at a heliocentric distance of 11.89 au. Comparative analysis with a dust model suggests that the comet may have been active since $\sim30$ au from the Sun. We derive a nucleus radius of $17\pm6$ km from the nucleus photometry, though this number is likely an overestimation due to the contamination from dust and gas. Continuing monitoring is encouraged in anticipation of the comet's forthcoming perihelion in 2024 April.


title: Low-Power Object Counting with Hierarchical Neural Networks

abstract: Deep Neural Networks (DNNs) can achieve state-of-the-art accuracy in many computer vision tasks, such as object counting. Object counting takes two inputs: an image and an object query and reports the number of occurrences of the queried object. To achieve high accuracy on such tasks, DNNs require billions of operations, making them difficult to deploy on resource-constrained, low-power devices. Prior work shows that a significant number of DNN operations are redundant and can be eliminated without affecting the accuracy. To reduce these redundancies, we propose a hierarchical DNN architecture for object counting. This architecture uses a Region Proposal Network (RPN) to propose regions-of-interest (RoIs) that may contain the queried objects. A hierarchical classifier then efficiently finds the RoIs that actually contain the queried objects. The hierarchy contains groups of visually similar object categories. Small DNNs are used at each node of the hierarchy to classify between these groups. The RoIs are incrementally processed by the hierarchical classifier. If the object in an RoI is in the same group as the queried object, then the next DNN in the hierarchy processes the RoI further; otherwise, the RoI is discarded. By using a few small DNNs to process each image, this method reduces the memory requirement, inference time, energy consumption, and number of operations with negligible accuracy loss when compared with the existing object counters.


title: Epidemiology of exposure to mixtures: we cant be casual about causality when using or testing methods

abstract: Background: There is increasing interest in approaches for analyzing the effect of exposure mixtures on health. A key issue is how to simultaneously analyze often highly collinear components of the mixture, which can create problems such as confounding by co-exposure and co-exposure amplification bias (CAB). Evaluation of novel mixtures methods, typically using synthetic data, is critical to their ultimate utility. Objectives: This paper aims to answer two questions. How do causal models inform the interpretation of statistical models and the creation of synthetic data used to test them? Are novel mixtures methods susceptible to CAB? Methods: We use directed acyclic graphs (DAGs) and linear models to derive closed form solutions for model parameters to examine how underlying causal assumptions affect the interpretation of model results. Results: The same beta coefficients estimated by a statistical model can have different interpretations depending on the assumed causal structure. Similarly, the method used to simulate data can have implications for the underlying DAG (and vice versa), and therefore the identification of the parameter being estimated with an analytic approach. We demonstrate that methods that can reproduce results of linear regression, such as Bayesian kernel machine regression and the new quantile g-computation approach, will be subject to CAB. However, under some conditions, estimates of an overall effect of the mixture is not subject to CAB and even has reduced uncontrolled bias. Discussion: Just as DAGs encode a priori subject matter knowledge allowing identification of variable control needed to block analytic bias, we recommend explicitly identifying DAGs underlying synthetic data created to test statistical mixtures approaches. Estimates of the total effect of a mixture is an important but relatively underexplored topic that warrants further investigation.


title: Abundance Analysis of New $r$-Process-Enhanced Stars from the HESP-GOMPA Survey

abstract: We present a study on the detailed chemical abundances of five new relatively bright $r$-process-enhanced stars that were initially observed as part of the SDSS/MARVELS pre-survey. These stars were selected, on the basis of their metallicities and carbon abundances, among a total of 60 stars, for high-resolution spectroscopic follow-up as part of the HESP-GOMPA survey (Hanle Echelle SPectrograph -- Galactic survey Of Metal Poor stArs). Here we discuss the three new $r$-I and two new $r$-II stars found in this survey. We have carried out a detailed abundance analysis for each of these stars, at a resolving power of $R \sim 30,000$, and compare our results to the existing literature. We could measure three of the first $r$-process-peak elements (Sr, Y and Zr) in all five stars, while Ba, Ce, Nd, Sm, Eu, and Dy could be detected among the second $r$-process-peak elements. Thorium could also be detected in one of the targets, which is found to be an actinide-boost star. We have carried out a comparative study among the sub-populations of the $r$-process-enhanced stars and other stars of the Milky Way halo population to constrain the origin of this class of objects. These bright $r$-process-enhanced stars provide an excellent opportunity to study the nucleosynthesis history of this population in great detail, and shed light on their chemical-enrichment histories.


title: Discriminating between Neutron Stars and Black Holes with Imperfect Knowledge of the Maximum Neutron Star Mass

abstract: Although gravitational-wave signals from exceptional low-mass compact binary coalescences, like GW170817, may carry matter signatures that differentiate the source from a binary black hole system, only one out of every eight events detected by the current Advanced LIGO and Virgo observatories are likely to have signal-to-noise ratios large enough to measure matter effects, even if they are present. Nonetheless, the systems' component masses will generally be constrained precisely. Constructing an explicit mixture model for the total rate density of merging compact objects, we develop a hierarchical Bayesian analysis to classify gravitational-wave sources according to the posterior odds that their component masses are drawn from different subpopulations. Accounting for current uncertainty in the maximum neutron star mass, and adopting different reasonable models for the total rate density, we examine two recent events from the LIGO-Virgo Collaboration's third observing run, GW190425 and GW190814. For population models with no overlap between the neutron star and black hole mass distributions, we typically find that there is a $\gtrsim 70\%$ chance that GW190425 was a binary neutron star merger rather than a neutron-star--black-hole merger. On the other hand, we find that there is a $\lesssim 6\%$ chance that GW190814 involved a slowly spinning neutron star, regardless of our assumed population model.


title: Phonons, Q-dependent Kondo spin fluctuations, and 4$\textit{f}$/phonon resonance in YbAl$_3$

abstract: The intermediate valence (IV) compound YbAl$_3$ exhibits nonintegral valence (Yb 4$f^{14-n_f}$ (5d6s)$^z$ where z = 2+n$_f$ = 2.75) in a moderately heavy (m* = 20-30me) ground state with a large Kondo temperature (T$_K$ ~ 500-600K). We have measured the magnetic fluctuations and the phonon spectra on single crystals of this material by time-of-flight inelastic neutron scattering (INS) and inelastic x-ray scattering (IXS). We find that at low temperature, the Kondo-scale spin fluctuations have a momentum (Q) dependence similar to that seen recently in the IV compound CePd$_3$ and which can be attributed to particle-hole excitations in a coherent itinerant 4$f$ correlated ground state. The Q-dependence disappears as the temperature is raised towards room temperature and the 4$f$ electron band states become increasingly incoherent. The measured phonons can be described adequately by a calculation based on standard DFT+$U$ density functional theory, without recourse to considering 4$f$ correlations dynamically. A low temperature magnetic peak observed in the neutron scattering at ~ 30meV shows dispersion identical to an optic phonon branch. This 4$f$/phonon resonance disappears for T > 150K. The phonons appear to remain unaffected by the resonance. We discuss several possibilities for the origin of this unusual excitation, including the idea that it arises from the large amplitude beating of the light Al atoms against the heavy Yb atoms, resulting in a dynamic 4$f$/3$p$ hybridization.


title: Smt-Switch: a solver-agnostic C++ API for SMT Solving

abstract: This extended abstract describes work in progress on Smt-Switch, an open-source, solver-agnostic API for SMT solving. Smt-Switch provides an abstract interface, which can be implemented by different SMT solvers. Smt-Switch provides simple, uniform, and high-performance access to SMT solving for applications in areas such as automated reasoning, planning, and formal verification. The interface allows the user to create, traverse, and manipulate terms, as well as to dynamically dispatch queries to different underlying SMT solvers.


title: LSTFCoDel: CoDel with LSTF-Style Priority Queuing

abstract: Congestion control is vastly important in computer networks. Arising naturally from the bursty nature of Internet traffic, congestion plagues not only the network edge, but also the network core. Many remedies have been proposed to fight congestion; active queue management (AQM) is one such proposal. AQM seeks to prevent congestion by actively avoiding it. Some queuing disciplines such as Random Early Detection (RED) will prematurely drop a random packet (with some probability) when the queue nears capacity to signal the sender to back off. However, RED utilizes queue length as a mechanism to indicate congestion. On the other hand, the Controlled Delay (CoDel) queuing discipline uses queuing delay as an indication of congestion. The problem with both RED and CoDel are that they indiscriminately treat all packets the same. Normally implemented using a FIFO queue, CoDel simply enqueues and dequeues packets in a first-come, first-served manner. Priority queuing can be carefully utilized to selectively service packets utilizing the very same metric CoDel uses for AQM, queuing delay. That said, Least Slack Time First (LSTF), a multi-processor scheduling algorithm employs priority scheduling, which coincidentally, is also based on delay. In the context of computer networks LSTF can be applied in the control plane or in the data plane. At the control plane, LSTF functions across the entire network, but in doing so requires all intermediary routers to implement it; LSTF also requires support at the packet level in terms of a slack entry. Within the data plane, LSTF can be implemented as a queuing mechanism based on delay spent in the router (just like CoDel AQM). This paper applies data plane level LSTF to CoDel AQM to enable delay-based packet classification within the confines of the CoDel AQM algorithm.


title: Improved bounds for noisy group testing with constant tests per item

abstract: The group testing problem is concerned with identifying a small set of infected individuals in a large population. At our disposal is a testing procedure that allows us to test several individuals together. In an idealized setting, a test is positive if and only if at least one infected individual is included and negative otherwise. Significant progress was made in recent years towards understanding the information-theoretic and algorithmic properties in this noiseless setting. In this paper, we consider a noisy variant of group testing where test results are flipped with certain probability, including the realistic scenario where sensitivity and specificity can take arbitrary values. Using a test design where each individual is assigned to a fixed number of tests, we derive explicit algorithmic bounds for two commonly considered inference algorithms and thereby naturally extend the results of Scarlett \& Cevher (2016) and Scarlett \& Johnson (2020). We provide improved performance guarantees for the efficient algorithms in these noisy group testing models -- indeed, for a large set of parameter choices the bounds provided in the paper are the strongest currently proved.


title: DATE: Defense Against TEmperature Side-Channel Attacks in DVFS Enabled MPSoCs

abstract: Given the constant rise in utilizing embedded devices in daily life, side channels remain a challenge to information flow control and security in such systems. One such important security flaw could be exploited through temperature side-channel attacks, where heat dissipation and propagation from the processing elements are observed over time in order to deduce security flaws. In our proposed methodology, DATE: Defense Against TEmperature side-channel attacks, we propose a novel approach of reducing spatial and temporal thermal gradient, which makes the system more secure against temperature side-channel attacks, and at the same time increases the reliability of the device in terms of lifespan. In this paper, we have also introduced a new metric, Thermal-Security-in-Multi-Processors (TSMP), which is capable of quantifying the security against temperature side-channel attacks on computing systems, and DATE is evaluated to be 139.24% more secure at the most for certain applications than the state-of-the-art, while reducing thermal cycle by 67.42% at the most.


title: Generative embeddings of brain collective dynamics using variational autoencoders

abstract: We consider the problem of encoding pairwise correlations between coupled dynamical systems in a low-dimensional latent space based on few distinct observations. We used variational autoencoders (VAE) to embed temporal correlations between coupled nonlinear oscillators that model brain states in the wake-sleep cycle into a two-dimensional manifold. Training a VAE with samples generated using two different parameter combinations resulted in an embedding that represented the whole repertoire of collective dynamics, as well as the topology of the underlying connectivity network. We first followed this approach to infer the trajectory of brain states measured from wakefulness to deep sleep from the two endpoints of this trajectory; next, we showed that the same architecture was capable of representing the pairwise correlations of generic Landau-Stuart oscillators coupled by complex network topology


title: Detecting Ongoing Events Using Contextual Word and Sentence Embeddings

abstract: This paper introduces the Ongoing Event Detection (OED) task, which is a specific Event Detection task where the goal is to detect ongoing event mentions only, as opposed to historical, future, hypothetical, or other forms or events that are neither fresh nor current. Any application that needs to extract structured information about ongoing events from unstructured texts can take advantage of an OED system. The main contribution of this paper are the following: (1) it introduces the OED task along with a dataset manually labeled for the task; (2) it presents the design and implementation of an RNN model for the task that uses BERT embeddings to define contextual word and contextual sentence embeddings as attributes, which to the best of our knowledge were never used before for detecting ongoing events in news; (3) it presents an extensive empirical evaluation that includes (i) the exploration of different architectures and hyperparameters, (ii) an ablation test to study the impact of each attribute, and (iii) a comparison with a replication of a state-of-the-art model. The results offer several insights into the importance of contextual embeddings and indicate that the proposed approach is effective in the OED task, outperforming the baseline models.


title: Deep reinforcement learning driven inspection and maintenance planning under incomplete information and constraints

abstract: Determination of inspection and maintenance policies for minimizing long-term risks and costs in deteriorating engineering environments constitutes a complex optimization problem. Major computational challenges include the (i) curse of dimensionality, due to exponential scaling of state/action set cardinalities with the number of components; (ii) curse of history, related to exponentially growing decision-trees with the number of decision-steps; (iii) presence of state uncertainties, induced by inherent environment stochasticity and variability of inspection/monitoring measurements; (iv) presence of constraints, pertaining to stochastic long-term limitations, due to resource scarcity and other infeasible/undesirable system responses. In this work, these challenges are addressed within a joint framework of constrained Partially Observable Markov Decision Processes (POMDP) and multi-agent Deep Reinforcement Learning (DRL). POMDPs optimally tackle (ii)-(iii), combining stochastic dynamic programming with Bayesian inference principles. Multi-agent DRL addresses (i), through deep function parametrizations and decentralized control assumptions. Challenge (iv) is herein handled through proper state augmentation and Lagrangian relaxation, with emphasis on life-cycle risk-based constraints and budget limitations. The underlying algorithmic steps are provided, and the proposed framework is found to outperform well-established policy baselines and facilitate adept prescription of inspection and intervention actions, in cases where decisions must be made in the most resource- and risk-aware manner.


title: D-NetPAD: An Explainable and Interpretable Iris Presentation Attack Detector

abstract: An iris recognition system is vulnerable to presentation attacks, or PAs, where an adversary presents artifacts such as printed eyes, plastic eyes, or cosmetic contact lenses to circumvent the system. In this work, we propose an effective and robust iris PA detector called D-NetPAD based on the DenseNet convolutional neural network architecture. It demonstrates generalizability across PA artifacts, sensors and datasets. Experiments conducted on a proprietary dataset and a publicly available dataset (LivDet-2017) substantiate the effectiveness of the proposed method for iris PA detection. The proposed method results in a true detection rate of 98.58\% at a false detection rate of 0.2\% on the proprietary dataset and outperfoms state-of-the-art methods on the LivDet-2017 dataset. We visualize intermediate feature distributions and fixation heatmaps using t-SNE plots and Grad-CAM, respectively, in order to explain the performance of D-NetPAD. Further, we conduct a frequency analysis to explain the nature of features being extracted by the network. The source code and trained model are available at https://github.com/iPRoBe-lab/D-NetPAD.


title: WattScale: A Data-driven Approach for Energy Efficiency Analytics of Buildings at Scale

abstract: Buildings consume over 40% of the total energy in modern societies, and improving their energy efficiency can significantly reduce our energy footprint. In this paper, we present \texttt{WattScale}, a data-driven approach to identify the least energy-efficient buildings from a large population of buildings in a city or a region. Unlike previous methods such as least-squares that use point estimates, \texttt{WattScale} uses Bayesian inference to capture the stochasticity in the daily energy usage by estimating the distribution of parameters that affect a building. Further, it compares them with similar homes in a given population. \texttt{WattScale} also incorporates a fault detection algorithm to identify the underlying causes of energy inefficiency. We validate our approach using ground truth data from different geographical locations, which showcases its applicability in various settings. \texttt{WattScale} has two execution modes -- (i) individual, and (ii) region-based, which we highlight using two case studies. For the individual execution mode, we present results from a city containing >10,000 buildings and show that more than half of the buildings are inefficient in one way or another indicating a significant potential from energy improvement measures. Additionally, we provide probable cause of inefficiency and find that 41\%, 23.73\%, and 0.51\% homes have poor building envelope, heating, and cooling system faults, respectively. For the region-based execution mode, we show that \texttt{WattScale} can be extended to millions of homes in the US due to the recent availability of representative energy datasets.


title: Deep Interactive Learning: An Efficient Labeling Approach for Deep Learning-Based Osteosarcoma Treatment Response Assessment

abstract: Osteosarcoma is the most common malignant primary bone tumor. Standard treatment includes pre-operative chemotherapy followed by surgical resection. The response to treatment as measured by ratio of necrotic tumor area to overall tumor area is a known prognostic factor for overall survival. This assessment is currently done manually by pathologists by looking at glass slides under the microscope which may not be reproducible due to its subjective nature. Convolutional neural networks (CNNs) can be used for automated segmentation of viable and necrotic tumor on osteosarcoma whole slide images. One bottleneck for supervised learning is that large amounts of accurate annotations are required for training which is a time-consuming and expensive process. In this paper, we describe Deep Interactive Learning (DIaL) as an efficient labeling approach for training CNNs. After an initial labeling step is done, annotators only need to correct mislabeled regions from previous segmentation predictions to improve the CNN model until the satisfactory predictions are achieved. Our experiments show that our CNN model trained by only 7 hours of annotation using DIaL can successfully estimate ratios of necrosis within expected inter-observer variation rate for non-standardized manual surgical pathology task.


title: Metric SYZ conjecture and non-archimedean geometry

abstract: We show that assuming a conjecture in non-archimedean geometry, then a metric formulation of the SYZ conjecture can be proved in large generality.


title: Trace Densities and Algebraic Index Theorems for Sheaves of Formal Cherednik Algebras

abstract: We show how a novel construction of the sheaf of Cherednik algebras on a quotient orbifold Y=X/G by virtue of formal geometry in author's prior work leads to results for the sheaf of Cherednik algebra which until recently were viewed as intractable. First, for every orbit type stratum in $X$, we define a trace density map for the Hochschild chain complex of the sheaf of Cherednik algebras, which generalises the standard Engeli-Felder's trace density construction for the sheaf of differential operators. Second, by means of the newly obtained trace density maps, we prove an isomorphism in the derived category of complexes of $\mathbb C_Y[[\hbar]]$-modules which computes the hypercohomology of the Hochschild chain complex $\mathcal{C}_{\bullet}$ of the sheaf of formal Cherednik algebras. We show that this hypercohomology is isomorphic to the Chen-Ruan cohomology of the orbifold $X/G$ with values in the ring of formal power series $\mathbb C[[\hbar]]$. We infer that the Hochschild chain complex of the sheaf of skew group algebras $\mathcal{D}_X\rtimes G$ has a well-defined Euler characteristic which is proportional to the topological Euler characteristic of $X/G$. Finally, we prove an algebraic index theorem.


title: Posterior Adaptation With New Priors

abstract: Classification approaches based on the direct estimation and analysis of posterior probabilities will degrade if the original class priors begin to change. We prove that a unique (up to scale) solution is possible to recover the data likelihoods for a test example from its original class posteriors and dataset priors. Given the recovered likelihoods and a set of new priors, the posteriors can be re-computed using Bayes' Rule to reflect the influence of the new priors. The method is simple to compute and allows a dynamic update of the original posteriors.


title: An Algebraic Approach for the Stability Analysis of BLDC Motor Controllers

abstract: This paper presents an algebraic technique to compute the maximum time-delay that can be accepted in the control loop of a Brushless DC Motor (BLDCM) speed controller before the closed loop response becomes unstable. Using a recently proposed time-delay stability analysis methodology, we derive accurate stability conditions for the BLDCM speed controller. The results of applying the new method show that tuning the PI controller for very fast response in the order of magnitude of the BLDCM mechanical time constant cause the time-delay to significantly affect the system stability.


title: Learn Faster and Forget Slower via Fast and Stable Task Adaptation

abstract: Training Deep Neural Networks (DNNs) is still highly time-consuming and compute-intensive. It has been shown that adapting a pretrained model may significantly accelerate this process. With a focus on classification, we show that current fine-tuning techniques make the pretrained models catastrophically forget the transferred knowledge even before anything about the new task is learned. Such rapid knowledge loss undermines the merits of transfer learning and may result in a much slower convergence rate compared to when the maximum amount of knowledge is exploited. We investigate the source of this problem from different perspectives and to alleviate it, introduce Fast And Stable Task-adaptation (FAST), an easy to apply fine-tuning algorithm. The paper provides a novel geometric perspective on how the loss landscape of source and target tasks are linked in different transfer learning strategies. We empirically show that compared to prevailing fine-tuning practices, FAST learns the target task faster and forgets the source task slower.


title: Bertrand's postulate and the existence of finite fields

abstract: This is an expository note discussing how the Erdos--Ramanujan proof of Bertrand's postulate may be adapted to show the existence of finite fields.