title: Many-body quantum teleportation via operator spreading in the traversable wormhole protocol

abstract: By leveraging shared entanglement between a pair of qubits, one can teleport a quantum state from one particle to another. Recent advances have uncovered an intrinsically many-body generalization of quantum teleportation, with an elegant and surprising connection to gravity. In particular, the teleportation of quantum information relies on many-body dynamics, which originate from strongly-interacting systems that are holographically dual to gravity; from the gravitational perspective, such quantum teleportation can be understood as the transmission of information through a traversable wormhole. Here, we propose and analyze a new mechanism for many-body quantum teleportation -- dubbed peaked-size teleportation. Intriguingly, peaked-size teleportation utilizes precisely the same type of quantum circuit as traversable wormhole teleportation, yet has a completely distinct microscopic origin: it relies upon the spreading of local operators under generic thermalizing dynamics and not gravitational physics. We demonstrate the ubiquity of peaked-size teleportation, both analytically and numerically, across a diverse landscape of physical systems, including random unitary circuits, the Sachdev-Ye-Kitaev model (at high temperatures), one-dimensional spin chains and a bulk theory of gravity with stringy corrections. Our results pave the way towards using many-body quantum teleportation as a powerful experimental tool for: (i) characterizing the size distributions of operators in strongly-correlated systems and (ii) distinguishing between generic and intrinsically gravitational scrambling dynamics. To this end, we provide a detailed experimental blueprint for realizing many-body quantum teleportation in both trapped ions and Rydberg atom arrays; effects of decoherence and experimental imperfections are analyzed.

categories: quant-ph cond-mat.quant-gas gr-qc hep-th


title: Zero sum sets in abelian groups

abstract: The distribution of cardinalities of zero-sum sets in abelian groups is completely determined. A complex summation involving the M\"obius function is given for the general abelian group, while in many special cases, including the case of elementary abelian groups, solved earlier by Li and Wan, it has a compact form. The proof involves two different M\"obius transforms, on positive integers and on set partitions.


title: An even lighter QCD axion

abstract: We explore whether the axion which solves the strong CP problem can naturally be much lighter than the canonical QCD axion. The $Z_\mathcal{N}$ symmetry proposed by Hook, with $\mathcal{N}$ mirror and degenerate worlds coexisting in Nature and linked by the axion field, is considered in terms of generic effective axion couplings. We show that the total potential is safely approximated by a single cosine in the large $\mathcal{N}$ limit, and we determine the analytical formula for the exponentially suppressed axion mass. The resulting universal enhancement of all axion interactions relative to those of the canonical QCD axion has a strong impact on the prospects of axion-like particle experiments such as ALPS II, IAXO and many others. The finite density axion potential is also analyzed and we show that the $Z_\mathcal{N}$ asymmetric background of high-density stellar environments sets already significant model-independent constraints: $3\le\mathcal{N}\lesssim47$ for an axion scale $f_a\lesssim 2.4\times10^{15}$ GeV, with tantalizing discovery prospects for any value of $f_a$ and down to $\mathcal{N}\sim9$ with future neutron star and gravitational wave data, down to the ultra-light mass region. In addition, two specific ultraviolet $Z_\mathcal{N}$ completions are developed: a composite axion one and a KSVZ-like model with improved Peccei-Quinn quality.


title: Long-distance entanglement of purification and reflected entropy in conformal field theory

abstract: Quantifying entanglement properties of mixed states in quantum field theory via entanglement of purification and reflected entropy is a new and challenging subject. In this work, we study both quantities for two spherical subregions far away from each other in the vacuum of a conformal field theory in any number of dimensions. Using lattice techniques, we find an elementary proof that the decay of both, the entanglement of purification and reflected entropy, is enhanced with respect to the mutual information behaviour by a logarithm of the distance between the subregions. In the case of the Ising spin chain at criticality and the related free fermion conformal field theory, we compute also the overall coefficients numerically for the both quantities of interest.

categories: hep-th cond-mat.stat-mech cond-mat.str-el quant-ph


title: Topological magnons for thermal Hall transport in frustrated magnets with bond-dependent interactions

abstract: Thermal transport in topologically-ordered phases of matter provides valuable insights as it can detect the charge-neutral quasiparticles that would not directly couple to electromagnetic probes. An important example is edge heat transport of Majorana fermions in a chiral spin liquid, which leads to a half-quantized thermal Hall conductivity. This signature is precisely what has recently been measured in $\alpha$-RuCl$_3$ under external magnetic fields. The plateau-like behavior of the half-quantized thermal Hall conductivity as a function of external magnetic field, and the peculiar sign change depending on the magnetic field orientation, has been proposed as strong evidence for the non-Abelian Kitaev spin liquid. Alternatively, for in-plane magnetic fields, it was theoretically shown that such a sign structure can also arise from topological magnons in the field-polarized state. In this work, we investigate the full implications of topological magnons as heat carriers on thermal transport measurements. We first prove analytically that for any commensurate order with a finite magnetic unit cell, reversing the field direction leads to a sign change in the magnon thermal Hall conductivity in two-dimensional systems. We verify this proof numerically with nontrivial magnetic orders as well as the field-polarized state in Kitaev magnets subjected to an in-plane field. In the case of a tilted magnetic field, whereby there exist both finite in-plane and out-of-plane field components, we find that the plateau-like behavior of the thermal Hall conductivity and the sign change upon reversing the in-plane component of the magnetic field arise in the partially-polarized state, as long as the in-plane field contribution to the Zeeman energy is significant. While these results are consistent with the experimental observations, we comment on other aspects requiring investigation in future studies.


title: Improved gravitational radiation time-scales II: spin-orbit contributions and environmental perturbations

abstract: Peters' formula is an analytical estimate of the time-scale of gravitational wave (GW)-induced coalescence of binary systems. It is used in countless applications, where the convenience of a simple formula outweighs the need for precision. However, many promising sources of the Laser Interferometer Space Antenna (LISA), such as supermassive black hole binaries and extreme mass-ratio inspirals (EMRIs), are expected to enter the LISA band with highly eccentric ($e \gtrsim$ 0.9) and highly relativistic orbits. These are exactly the two limits in which Peters' estimate performs the worst. In this work, we expand upon previous results and give simple analytical fits to quantify how the inspiral time-scale is affected by the relative 1.5 post-Newtonian (PN) hereditary fluxes and spin-orbit couplings. We discuss several cases that demand a more accurate GW time-scale. We show how this can have a major influence on quantities that are relevant for LISA event-rate estimates, such as the EMRI critical semi-major axis. We further discuss two types of environmental perturbations that can play a role in the inspiral phase: the gravitational interaction with a third massive body and the energy loss due to dynamical friction and torques from a surrounding gas medium ubiquitous in galactic nuclei. With the aid of PN corrections to the time-scale in vacuum, we find simple analytical expressions for the regions of phase space in which environmental perturbations are of comparable strength to the effects of any particular PN order, being able to qualitatively reproduce the results of much more sophisticated analyses.


title: Dynamic electron correlations with charge order wavelength along all directions in the copper oxide plane

abstract: In strongly correlated systems the strength of Coulomb interactions between electrons, relative to their kinetic energy, plays a central role in determining their emergent quantum mechanical phases. We perform resonant x-ray scattering on Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$, a prototypical cuprate superconductor, to probe electronic correlations within the CuO$_2$ plane. We discover a dynamic quasi-circular pattern in the $x$-$y$ scattering plane with a radius that matches the wave vector magnitude of the well-known static charge order. Along with doping- and temperature-dependent measurements, our experiments reveal a picture of charge order competing with superconductivity where short-range domains along $x$ and $y$ can dynamically rotate into any other in-plane direction. This quasi-circular spectrum, a hallmark of Brazovskii-type fluctuations, has immediate consequences to our understanding of rotational and translational symmetry breaking in the cuprates. We discuss how the combination of short- and long-range Coulomb interactions results in an effective non-monotonic potential that may determine the quasi-circular pattern.


title: ${\tt bajes}$: Bayesian inference of multimessenger astrophysical data, methods and application to gravitational-waves

abstract: We present ${\tt bajes}$, a parallel and lightweight framework for Bayesian inference of multimessenger transients. ${\tt bajes}$ is a Python modular package with minimal dependencies on external libraries adaptable to the majority of the Bayesian models and to various sampling methods. We describe the general workflow and the parameter estimation pipeline for compact-binary-coalescence gravitational-wave transients. The latter is validated against injections of binary black hole and binary neutron star waveforms, including confidence interval tests that demonstrates the inference is well-calibrated. Binary neutron star postmerger injections are also studied using a network of five detectors made of LIGO, Virgo, KAGRA and Einstein Telescope. Postmerger signals will be detectable for sources at ${\lesssim}80\,$Mpc, with Einstein Telescope contributing over 90\% of the total signal-to-noise ratio. As a full scale application, we re-analyze the GWTC-1 black hole transients using the effective-one-body ${\tt TEOBResumS}$ approximant, and reproduce selected results with other approximants. ${\tt bajes}$ inferences are consistent with previous results; the direct comparison of ${\tt bajes}$ and ${\tt bilby}$ analyses of GW150914 shows a maximum Jensen-Shannon divergence of $5.2{\times}10^{-4}$. GW170817 is re-analyzed using ${\tt TaylorF2}$ with 5.5PN point-mass and 7.5PN tides, ${\tt TEOBResumSPA}$, and ${\tt IMRPhenomPv2\_NRTidal}$ with different cutoff-frequencies of $1024\,$Hz and $2048\,$Hz. We find that the former choice minimizes systematics on the reduced tidal parameter, while a larger amount of tidal information is gained with the latter choice. ${\tt bajes}$ can perform these analyses in about 1~day using 128 CPUs.


title: Self-consistent determination of proton and nuclear PDFs at the Electron Ion Collider

abstract: We quantify the impact of unpolarized lepton-proton and lepton-nucleus inclusive deep-inelastic scattering (DIS) cross section measurements from the future Electron-Ion Collider (EIC) on the proton and nuclear parton distribution functions (PDFs). To this purpose we include neutral- and charged-current DIS pseudodata in a self-consistent set of proton and nuclear global PDF determinations based on the NNPDF methodology. We demonstrate that the EIC measurements will reduce the uncertainty of the light quark PDFs of the proton at large values of the momentum fraction $x$, and, more significantly, of the quark and gluon PDFs of heavy nuclei, especially at small and large $x$. We illustrate the implications of the improved precision of nuclear PDFs for the interaction of ultra-high energy cosmic neutrinos with matter.


title: Hamiltonian reconstruction as metric for variational studies

abstract: Variational approaches are among the most powerful modern techniques to approximately solve quantum many-body problems. These encompass both variational states based on tensor or neural networks, and parameterized quantum circuits in variational quantum eigensolvers. However, self-consistent evaluation of the quality of variational wavefunctions is a notoriously hard task. Using a recently developed Hamiltonian reconstruction method, we propose a multi-faceted approach to evaluating the quality of neural-network based wavefunctions. Specifically, we consider convolutional neural network (CNN) and restricted Boltzmann machine (RBM) states trained on a square lattice spin-1/2 $J_1$-$J_2$ Heisenberg model. We find that the reconstructed Hamiltonians are typically less frustrated, and have easy-axis anisotropy near the high frustration point. Furthermore, the reconstructed Hamiltonians suppress quantum fluctuations in the large $J_2$ limit. Our results highlight the critical importance of the wavefunction's symmetry. Moreover, the multi-faceted insight from the Hamiltonian reconstruction reveals that a variational wave function can fail to capture the true ground state through suppression of quantum fluctuations.


title: Measurements of Partial Branching Fractions of Inclusive $B \to X_u \, \ell^+\, \nu_{\ell}$ Decays with Hadronic Tagging

abstract: We present measurements of partial branching fractions of inclusive semileptonic $B \to X_u \, \ell^+\, \nu_{\ell}$ decays using the full Belle data set of 711 fb$^{-1}$ of integrated luminosity at the $\Upsilon(4S)$ resonance and for $\ell = e, \mu$. Inclusive semileptonic $B \to X_u \, \ell^+\, \nu_{\ell}$ decays are CKM suppressed and measurements are complicated by the large background from CKM-favored $B \to X_c \, \ell^+\, \nu_{\ell}$ transitions, which have a similar signature. Using machine learning techniques, we reduce this and other backgrounds effectively, whilst retaining access to a large fraction of the $B \to X_u \, \ell^+\, \nu_{\ell}$ phase space and high signal efficiency. We measure partial branching fractions in three phase-space regions covering about $31\%$ to $86\%$ of the accessible $B \to X_u \, \ell^+\, \nu_{\ell}$ phase space. The most inclusive measurement corresponds to the phase space with lepton energies of $E_\ell^B > 1 $ GeV, and we obtain $\Delta \mathcal{B}(B \to X_u \ell^+ \, \nu_\ell) = \left( 1.59 \pm 0.07 \pm 0.16 \right) \times 10^{-3}$ from a two-dimensional fit of the hadronic mass spectrum and the four-momentum-transfer squared distribution, with the uncertainties denoting the statistical and systematic error. We find $\left| V_{ub} \right| = \left( 4.10 \pm 0.09 \pm 0.22 \pm 0.15 \right) \times 10^{-3}$ from an average of four calculations for the partial decay rate with the third uncertainty denoting the average theory error. This value is higher but compatible with the determination from exclusive semileptonic decays within 1.3 standard deviations. In addition, we report charmless inclusive partial branching fractions separately for $B^+$ and $B^0$ mesons as well as for electron and muon final states. No isospin breaking or lepton flavor universality violating effects are observed.


title: Security in Quantum Cryptography

abstract: Quantum cryptography exploits principles of quantum physics for the secure processing of information. A prominent example is secure communication, i.e., the task of transmitting confidential messages from one location to another. The cryptographic requirement here is that the transmitted messages remain inaccessible to anyone other than the designated recipients, even if the communication channel is untrusted. In classical cryptography, this can usually only be guaranteed under computational hardness assumptions, e.g., that factoring large integers is infeasible. In contrast, the security of quantum cryptography relies entirely on the laws of quantum mechanics. Here we review this physical notion of security, focusing on quantum key distribution and secure communication.


title: MICROSCOPE's constraint on a short-range fifth force

abstract: The MICROSCOPE experiment was designed to test the weak equivalence principle in space, by comparing the low-frequency dynamics of cylindrical "free-falling" test masses controlled by electrostatic forces. We use data taken during technical sessions aimed at estimating the electrostatic stiffness of MICROSCOPE's sensors to constrain a short-range Yukawa deviation from Newtonian gravity. We take advantage of the fact that in the limit of small displacements, the gravitational interaction (both Newtonian and Yukawa-like) between nested cylinders is linear, and thus simply characterised by a stiffness. By measuring the total stiffness of the forces acting on a test mass as it moves, and comparing it with the theoretical electrostatic stiffness (expected to dominate), it is a priori possible to infer constraints on the Yukawa potential parameters. However, we find that measurement uncertainties are dominated by the gold wires used to control the electric charge of the test masses, though their related stiffness is indeed smaller than the expected electrostatic stiffness. Moreover, we find a non-zero unaccounted for stiffness that depends on the instrument's electric configuration, hinting at the presence of patch-field effects. Added to significant uncertainties on the electrostatic model, they only allow for poor constraints on the Yukawa potential. This is not surprising, as MICROSCOPE was not designed for this measurement, but this analysis is the first step to new experimental searches for non-Newtonian gravity in space.


title: Constraints on chameleon gravity from the measurement of the electrostatic stiffness of the MICROSCOPE mission accelerometers

abstract: This article is dedicated to the use the MICROSCOPE mission's data to test chameleon theory of gravity. We take advantage of the technical sessions aimed to characterize the electrostatic stiffness of MICROSCOPE's instrument intrinsic to its capacitive measurement system. Any discrepancy between the expected and measured stiffness may result from unaccounted-for contributors, i.e. extra-forces. This work considers the case of chameleon gravity as a possible contributor. It was previously shown that in situations similar to these measurement sessions, a chameleon fifth force appears and acts as a stiffness for small displacements. The magnitude of this new component of the stiffness is computed over the chameleon's parameter space. It allows us to derive constraints by excluding any force inconsistent with the MICROSCOPE data. As expected --since MICROSCOPE was not designed for the purpose of such an analysis--, these new bounds are not competitive with state-of-the-art constraints, but they could be improved by a better estimation of all effects at play in these sessions. Hence our work illustrates this novel technique as a new way of constraining fifth forces.


title: Beams with propagation-invariant transverse polarization pattern

abstract: A wide class of nonuniformly totally polarized beams is introduced that preserve their transverse polarization pattern during paraxial propagation. They are obtained as suitable combinations of Gaussian modes and find applications in polarimetric techniques that use a single input beam for the determination of the Mueller matrix of a homogeneous sample. The class also includes beams that present all possible polarization states across their transverse section (Full-Poincar\`e beams). An example of such beams and its use in polarimetry is discussed in detail. The requirement of invariance of the polarization pattern can be limited to the propagation in the far field, in which case less restrictive conditions are found and a wider class of beams is obtained.


title: Multiband charge density wave exposed in a transition metal dichalcogenide

abstract: In the presence of multiple bands, well-known electronic instabilities may acquire new complexity. While multiband superconductivity is the subject of extensive studies, the possibility of multiband charge density waves (CDWs) has been largely ignored so far. Here, combining energy dependent scanning tunnelling microscopy (STM) topography with a simple model of the charge modulations and a self-consistent calculation of the CDW gap, we find evidence for a multiband CDW in 2H-NbSe$_2$. This CDW not only involves the opening of a gap on the inner band around the K-point, but also on the outer band. This leads to spatially out-of-phase charge modulations from electrons on these two bands, which we detect through a characteristic energy dependence of the CDW contrast in STM images.


title: Decay in norm of transfer operators for semiflows

abstract: We establish exponential decay in H\"older norm of transfer operators applied to smooth observables of uniformly and nonuniformly expanding semiflows with exponential decay of correlations.


title: Gesture Recognition in Robotic Surgery: a Review

abstract: Objective: Surgical activity recognition is a fundamental step in computer-assisted interventions. This paper reviews the state-of-the-art in methods for automatic recognition of fine-grained gestures in robotic surgery focusing on recent data-driven approaches and outlines the open questions and future research directions. Methods: An article search was performed on 5 bibliographic databases with the following search terms: robotic, robot-assisted, JIGSAWS, surgery, surgical, gesture, fine-grained, surgeme, action, trajectory, segmentation, recognition, parsing. Selected articles were classified based on the level of supervision required for training and divided into different groups representing major frameworks for time series analysis and data modelling. Results: A total of 52 articles were reviewed. The research field is showing rapid expansion, with the majority of articles published in the last 4 years. Deep-learning-based temporal models with discriminative feature extraction and multi-modal data integration have demonstrated promising results on small surgical datasets. Currently, unsupervised methods perform significantly less well than the supervised approaches. Conclusion: The development of large and diverse open-source datasets of annotated demonstrations is essential for development and validation of robust solutions for surgical gesture recognition. While new strategies for discriminative feature extraction and knowledge transfer, or unsupervised and semi-supervised approaches, can mitigate the need for data and labels, they have not yet been demonstrated to achieve comparable performance. Important future research directions include detection and forecast of gesture-specific errors and anomalies. Significance: This paper is a comprehensive and structured analysis of surgical gesture recognition methods aiming to summarize the status of this rapidly evolving field.


title: High $H_0$ Values from CMB E-mode Data: A Clue for Resolving the Hubble Tension?

abstract: The E-mode (EE) CMB power spectra measured by Planck, ACTPol, and SPTpol constrain the Hubble constant to be $70.0\pm2.7$, $72.4^{+3.9}_{-4.8}$, and $73.1^{+3.3}_{-3.9}$ km s$^{-1}$ Mpc$^{-1}$ within the standard $\Lambda$CDM model (posterior mean and central 68% interval bounds). These values are higher than the constraints from the Planck temperature (TT) power spectrum, and consistent with the Cepheid-supernova distance ladder measurement $H_0=73.2\pm1.3$ km s$^{-1}$ Mpc$^{-1}$. If this preference for a higher value was strengthened in a joint analysis it could provide an intriguing hint at the resolution of the Hubble disagreement. We show, however, that combining the Planck, ACTPol, and SPTpol EE likelihoods yields $H_0=68.7\pm1.3$ km s$^{-1}$ Mpc$^{-1}$, $2.4\sigma$ lower than the distance ladder measurement. This is due to different degeneracy directions across the full parameter space, particularly involving the baryon density, $\Omega_bh^2$, and scalar tilt, $n_s$, arising from sensitivity to different multipole ranges. We show that the E-mode $\Lambda$CDM constraints are consistent across the different experiments within $1.4\sigma$, and with the Planck TT results at $0.8\sigma$. Combining the Planck, ACTPol, and SPTpol EE data constrains the phenomenological lensing amplitude, $A_L=0.89\pm0.10$, consistent with the expected value of unity.


title: You Only Query Once: Effective Black Box Adversarial Attacks with Minimal Repeated Queries

abstract: Researchers have repeatedly shown that it is possible to craft adversarial attacks on deep classifiers (small perturbations that significantly change the class label), even in the "black-box" setting where one only has query access to the classifier. However, all prior work in the black-box setting attacks the classifier by repeatedly querying the same image with minor modifications, usually thousands of times or more, making it easy for defenders to detect an ensuing attack. In this work, we instead show that it is possible to craft (universal) adversarial perturbations in the black-box setting by querying a sequence of different images only once. This attack prevents detection from high number of similar queries and produces a perturbation that causes misclassification when applied to any input to the classifier. In experiments, we show that attacks that adhere to this restriction can produce untargeted adversarial perturbations that fool the vast majority of MNIST and CIFAR-10 classifier inputs, as well as in excess of $60-70\%$ of inputs on ImageNet classifiers. In the targeted setting, we exhibit targeted black-box universal attacks on ImageNet classifiers with success rates above $20\%$ when only allowed one query per image, and $66\%$ when allowed two queries per image.


title: Optimistic Policy Iteration for MDPs with Acyclic Transient State Structure

abstract: We consider Markov Decision Processes (MDPs) in which every stationary policy induces the same graph structure for the underlying Markov chain and further, the graph has the following property: if we replace each recurrent class by a node, then the resulting graph is acyclic. For such MDPs, we prove the convergence of the stochastic dynamics associated with a version of optimistic policy iteration (OPI), suggested in Tsitsiklis (2002), in which the values associated with all the nodes visited during each iteration of the OPI are updated.


title: FA2PbBr4: synthesis, structure and unusual optical properties of two polymorphs of formamidinium-based layered (110) hybrid perovskite

abstract: Small cations such as guanidinium and cesium can act as templating cations to form low dimensional phases (2D, 1D, 0D) in the case of excess of organic halides. However, such phases with the widely used formamidinium (FA+) cation have not been reported so far. In this study, we discovered two novel low dimensional phases with a composition of FA2PbBr4 and investigated the prerequisites of their formation upon crystallization of FABr-excessive solutions of FAPbBr3. We found that both phases have the structure of (110) layered perovskite but are represented by two different polymorphs with eclipsed and staggered arrangement of adjacent layers. It was shown that FA2PbBr4 phases usually exist in a labile equilibrium with FAPbBr3 3D perovskite and can form composites with it. The optical properties of both polymorphs were comprehensively studied by means of absorption spectroscopy, diffuse reflection spectroscopy and photoluminescence spectroscopy. DFT calculations were applied to investigate the band structure of the FA2PbBr4 and to corroborate the conclusions on their optoelectronic properties. As a result, we found that FA2PbBr4 phases irradiated by UV can exhibit effective green photoluminescence due to a transfer of excitation energy to defective states or 3D perovskite inclusions.


title: Electroluminescence yield in pure krypton

abstract: The krypton electroluminescence yield was studied, at room temperature, as a function of electric field in the gas scintillation gap. A large area avalanche photodiode has been used to allow the simultaneous detection of the electroluminescence pulses as well as the direct interaction of x-rays, the latter being used as a reference for the calculation of the number of charge carriers produced by the electroluminescence pulses and, thus, the determination of the number of photons impinging the photodiode. An amplification parameter of 113 photons per kV per drifting electron and a scintillation threshold of 2.7 Td ( 0.7 kV/cm/bar at 293 K ) was obtained, in good agreement with the simulation data reported in the literature. On the other hand, the ionisation threshold in krypton was found to be around 13.5 Td (3.4 kV/cm/bar), less than what had been obtained by the most recent simulation work-package. The krypton amplification parameter is about 80% and 140% of those measured for xenon and argon, respectively. The electroluminescence yield in krypton is of great importance for modeling krypton-based double-phase or high-pressure gas detectors, which may be used in future rare event detection experiments.


title: On the functional equation for classical orthogonal polynomials on lattices

abstract: Necessary and sufficient conditions for the regularity of solutions of the functional equation appearing in the theory of classical orthogonal polynomials on lattices are stated. Moreover, the functional Rodrigues formula and a closed formula for the recurrence coefficients are presented.


title: Dynamic imaging using a deep generative SToRM (Gen-SToRM) model

abstract: We introduce a generative smoothness regularization on manifolds (SToRM) model for the recovery of dynamic image data from highly undersampled measurements. The model assumes that the images in the dataset are non-linear mappings of low-dimensional latent vectors. We use the deep convolutional neural network (CNN) to represent the non-linear transformation. The parameters of the generator as well as the low-dimensional latent vectors are jointly estimated only from the undersampled measurements. This approach is different from traditional CNN approaches that require extensive fully sampled training data. We penalize the norm of the gradients of the non-linear mapping to constrain the manifold to be smooth, while temporal gradients of the latent vectors are penalized to obtain a smoothly varying time-series. The proposed scheme brings in the spatial regularization provided by the convolutional network. The main benefit of the proposed scheme is the improvement in image quality and the orders-of-magnitude reduction in memory demand compared to traditional manifold models. To minimize the computational complexity of the algorithm, we introduce an efficient progressive training-in-time approach and an approximate cost function. These approaches speed up the image reconstructions and offers better reconstruction performance.


title: MF-Net: Compute-In-Memory SRAM for Multibit Precision Inference using Memory-immersed Data Conversion and Multiplication-free Operators

abstract: We propose a co-design approach for compute-in-memory inference for deep neural networks (DNN). We use multiplication-free function approximators based on ell_1 norm along with a co-adapted processing array and compute flow. Using the approach, we overcame many deficiencies in the current art of in-SRAM DNN processing such as the need for digital-to-analog converters (DACs) at each operating SRAM row/column, the need for high precision analog-to-digital converters (ADCs), limited support for multi-bit precision weights, and limited vector-scale parallelism. Our co-adapted implementation seamlessly extends to multi-bit precision weights, it doesn't require DACs, and it easily extends to higher vector-scale parallelism. We also propose an SRAM-immersed successive approximation ADC (SA-ADC), where we exploit the parasitic capacitance of bit lines of SRAM array as a capacitive DAC. Since the dominant area overhead in SA-ADC comes due to its capacitive DAC, by exploiting the intrinsic parasitic of SRAM array, our approach allows low area implementation of within-SRAM SA-ADC. Our 8$\times$62 SRAM macro, which requires a 5-bit ADC, achieves $\sim$105 tera operations per second per Watt (TOPS/W) with 8-bit input/weight processing at 45 nm CMOS.


title: Facilitating Knowledge Sharing from Domain Experts to Data Scientists for Building NLP Models

abstract: Data scientists face a steep learning curve in understanding a new domain for which they want to build machine learning (ML) models. While input from domain experts could offer valuable help, such input is often limited, expensive, and generally not in a form readily consumable by a model development pipeline. In this paper, we propose Ziva, a framework to guide domain experts in sharing essential domain knowledge to data scientists for building NLP models. With Ziva, experts are able to distill and share their domain knowledge using domain concept extractors and five types of label justification over a representative data sample. The design of Ziva is informed by preliminary interviews with data scientists, in order to understand current practices of domain knowledge acquisition process for ML development projects. To assess our design, we run a mix-method case-study to evaluate how Ziva can facilitate interaction of domain experts and data scientists. Our results highlight that (1) domain experts are able to use Ziva to provide rich domain knowledge, while maintaining low mental load and stress levels; and (2) data scientists find Ziva's output helpful for learning essential information about the domain, offering scalability of information, and lowering the burden on domain experts to share knowledge. We conclude this work by experimenting with building NLP models using the Ziva output by our case study.


title: Monte-Carlo Applications for Partially Polarised Inverse External-Compton Scattering (MAPPIES) -- I. Description of the code and First Results

abstract: The radiation mechanisms responsible for the multiwavelength emission from relativistic jet sources are poorly understood. The modelling of the spectral energy distributions (SEDs) and light curves alone is not adequate to distinguish between existing models. Polarisation in the $X$-ray and $\gamma$-ray regime of these sources may provide new and unique information about the jet physics and radiation mechanisms. Several upcoming projects will be able to deliver polarimetric measurements of the brightest $X$-ray sources, including active galactic nuclei (AGN) jets and $\gamma$-ray bursts (GRBs). This article describes the development of a new Monte-Carlo code -- MAPPIES (Monte-Carlo Applications for Partially Polarised Inverse External-Compton Scattering) -- for polarisation-dependent Compton scattering in relativistic jet sources. Generic results for Compton polarisation in the Thomson and Klein-Nishina regimes are presented.


title: Path-dependent Hamilton-Jacobi equations with super-quadratic growth in the gradient and the vanishing viscosity method

abstract: The non-exponential Schilder-type theorem in Backhoff-Veraguas, Lacker and Tangpi [Ann. Appl. Probab., 30 (2020), pp. 1321-1367] is expressed as a convergence result for path-dependent partial differential equations with appropriate notions of generalized solutions. This entails a non-Markovian counterpart to the vanishing viscosity method. We show uniqueness of maximal subsolutions for path-dependent viscous Hamilton-Jacobi equations related to convex super-quadratic backward stochastic differential equations. We establish well-posedness for the Hamilton-Jacobi-Bellman equation associated to a Bolza problem of the calculus of variations with path-dependent terminal cost. In particular, uniqueness among lower semi-continuous solutions holds and state constraints are admitted.


title: Monte-Carlo Applications for Partially Polarized Inverse External-Compton Scattering (MAPPIES) II -- Application to the UV/Soft X-ray Excess in Blazar Spectra

abstract: The spectral energy distributions (SEDs) of some blazars exhibit an ultraviolet (UV) and/or soft X-ray excess, which can be modelled with different radiation mechanisms. Polarization measurements of the UV/X-ray emission from blazars may provide new and unique information about the astrophysical environment of blazar jets, and could thus help to distinguish between different emission scenarios. In this paper, a new Monte-Carlo code -- MAPPIES (Monte-Carlo Applications for Partially Polarized Inverse External-Compton Scattering) -- for polarization-dependent Compton scattering is used to simulate the polarization signatures in a model where the UV/soft X-ray excess arises from the bulk Compton process. Predictions of the expected polarization signatures of Compton emission from the soft X-ray excess in the SED of AO 0235+164, and the UV excess in the SED of 3C 279 are made for upcoming and proposed polarimetry missions.


title: A Search for correlations between turbulence and star formation in LITTLE THINGS dwarf irregular galaxies

abstract: Turbulence has the potential for creating gas density enhancements that initiate cloud and star formation (SF), and it can be generated locally by SF. To study the connection between turbulence and SF, we looked for relationships between SF traced by FUV images, and gas turbulence traced by kinetic energy density (KED) and velocity dispersion ($v_{disp}$) in the LITTLE THINGS sample of nearby dIrr galaxies. We performed 2D cross-correlations between FUV and KED images, measured cross-correlations in annuli to produce correlation coefficients as a function of radius, and determined the cumulative distribution function of the cross correlation value. We also plotted on a pixel-by-pixel basis the locally excess KED, $v_{disp}$, and HI mass surface density, $\Sigma_{\rm HI}$, as determined from the respective values with the radial profiles subtracted, versus the excess SF rate density $\Sigma_{\rm SFR}$, for all regions with positive excess $\Sigma_{\rm SFR}$. We found that $\Sigma_{\rm SFR}$ and KED are poorly correlated. The excess KED associated with SF implies a $\sim0.5$% efficiency for supernova energy to pump local HI turbulence on the scale of resolution here, which is a factor of $\sim2$ too small for all of the turbulence on a galactic scale. The excess $v_{disp}$ in SF regions is also small, only $\sim0.37$ km s$^{-1}$. The local excess in $\Sigma_{\rm HI}$ corresponding to an excess in $\Sigma_{\rm SFR}$ is consistent with an HI consumption time of $\sim1.6$ Gyr in the inner parts of the galaxies. The similarity between this timescale and the consumption time for CO implies that CO-dark molecular gas has comparable mass to HI in the inner disks.


title: MeV astronomy with Herwig?

abstract: GeV-scale dark matter is an increasingly attractive target for direct detection, indirect detection, and collider searches. Especially for masses in the MeV to GeV range, indirect detection is expected to give a leading constraint. In that range, dark matter annihilations into hadronic final states produce a challenging zoo of light hadronic resonances. With an update of Herwig7, we provide a modeling of these processes and study energy spectra from annihilation through a vector mediator. We cover dark matter masses between 250 MeV and 5 GeV and include an error estimate. This opens up new opportunities for fully studying indirect detection signals for hidden dark matter sectors with vector mediators in the sub-GeV range.


title: Example of an Highly Branching CD Space

abstract: Ketterer and Rajala showed an example of metric measure space, satisfying the measure contraction property $MCP(0,3)$, that has different topological dimensions at different regions of the space. In this article I propose a refinement of that example, which satisfies the $CD(0,\infty)$ condition, proving the non-constancy of topological dimension for CD spaces. This example also shows that the weak curvature dimension bound, in the sense of Lott-Sturm-Villani, is not sufficient to deduce any reasonable non-branching condition. Moreover, it allows to answer to some open question proposed by Schultz, about strict curvature dimension bounds and their stability with respect to the measured Gromov Hausdorff convergence.


title: Error estimates for the Smagorinsky turbulence model: enhanced stability through scale separation and numerical stabilization

abstract: In the present work we show some results on the effect of the Smagorinsky model on the stability of the associated perturbation equation. We show that in the presence of a spectral gap, such that the flow can be decomposed in a large scale with moderate gradient and a small amplitude fine scale with arbitratry gradient, the Smagorinsky model admits stability estimates for perturbations, with exponential growth depending only on the large scale gradient. We then show in the context of stabilized finite element methods that the same result carries over to the approximation and that in this context, for suitably chosen finite element spaces the Smagorinsky model acts as a stabilizer yielding close to optimal error estimates in the $L^2$-norm for smooth flows in the pre-asymptotic high Reynolds number regime.


title: Multi-Wavelength Photometry Derived from Monochromatic Kepler Data

abstract: The Kepler mission has provided a wealth of data, revealing new insights in time-domain astronomy. However, Kepler's single band-pass has limited studies to a single wavelength. In this work we build a data-driven, pixel-level model for the Pixel Response Function (PRF) of Kepler targets, modeling the image data from the spacecraft. Our model is sufficiently flexible to capture known detector effects, such as non-linearity, intra-pixel sensitivity variations, and focus change. In theory, the shape of the Kepler PRF should also be weakly wavelength dependent, due to optical chromatic aberration and wavelength dependent detector response functions. We are able to identify these predicted shape changes to the PRF using the residuals between Kepler data and our model. In this work, we show that these PRF changes correspond to wavelength variability in Kepler targets using a small sample of eclipsing binaries. Using our model, we demonstrate that pixel-level light curves of eclipsing binaries show variable eclipse depths, ellipsoidal modulation and limb darkening. These changes at the pixel level are consistent with multi-wavelength photometry. Our work suggests each pixel in the Kepler data of a single target has a different effective wavelength, ranging from $\approx$ 550-750 $nm$. In this proof of concept, we demonstrate our model, and discuss possible use cases for the wavelength dependent Pixel Response Function of Kepler. These use cases include characterizing variable systems, and vetting exoplanet discoveries at the pixel level. The chromatic PRF of Kepler is due to weak wavelength dependence in the optical systems and detector of the telescope, and similar chromatic PRFs are expected in other similar telescopes, notably the NASA TESS telescope.


title: A Statistical Analysis of the Nulling Pulsar Population

abstract: Approximately 8% of the $\sim$2800 known pulsars exhibit "nulling," a temporary broadband cessation of normal pulsar emission. Nulling behaviour can be coarsely quantified by the nulling fraction, which describes the percentage of time a given pulsar will be found in a null state. In this paper, we perform the most thorough statistical analysis thus far of the properties of 141 known nulling pulsars. We find weak, non-linear correlations between nulling fraction and pulse width, as well as nulling fraction and spin period which could be attributed to selection effects. We also further investigate a recently-hypothesized gap at 40% nulling fraction. While a local minimum does exist in the distribution, we cannot confirm a consistent and unique break in the distribution when we investigate with univariate and multivariate clustering methods, nor can we prove the existence of two statistically distinct populations about this minimum. Using the same methods, we find that nulling pulsars are a statistically different population from normal, radio, non-nulling pulsars, which has never been quantitatively verified. In addition, we summarize the findings of the prior nulling pulsar statistics literature, which are notoriously contradictory. This study, in context, furthers the idea that nulling fraction alone does not contain enough information to describe the behaviour of a nulling pulsar and that other parameters such as null lengths and null randomness, in addition to a better understanding of selection effects, are required to fully understand this phenomenon.


title: Recovery of Power Flow to Critical Infrastructures using Mode-dependent Droop-based Inverters

abstract: Recovery of power flow to critical infrastructures, after grid failure, is a crucial need arising in scenarios that are increasingly becoming more frequent. This article proposes a power transition and recovery strategy by proposing a mode-dependent droop control-based inverters. The control strategy of inverters achieves the following objectives 1) regulate the output active and reactive power by the droop-based inverters to a desired value while operating in on-grid mode 2) seamless transition and recovery of power flow injections into the critical loads in the network by inverters operating in off-grid mode after the main grid fails; 3) require minimal information of grid/network status and conditions for the mode transition of droop control. A framework for assessing the stability of the system and to guide the choice of parameters for controllers is developed using control-oriented modeling. A comprehensive controller hardware-in-the-loop-based real-time simulation study on a test-system based on the realistic electrical network of M-Health Fairview, University of Minnesota Medical Center, corroborates the efficacy of the proposed controller strategy.


title: Model Adaptation for Image Reconstruction using Generalized Stein's Unbiased Risk Estimator

abstract: Deep learning image reconstruction algorithms often suffer from model mismatches when the acquisition scheme differs significantly from the forward model used during training. We introduce a Generalized Stein's Unbiased Risk Estimate (GSURE) loss metric to adapt the network to the measured k-space data and minimize model misfit impact. Unlike current methods that rely on the mean square error in kspace, the proposed metric accounts for noise in the measurements. This makes the approach less vulnerable to overfitting, thus offering improved reconstruction quality compared to schemes that rely on mean-square error. This approach may be useful to rapidly adapt pre-trained models to new acquisition settings (e.g., multi-site) and different contrasts than training data


title: A Metric Stability Result for the Very Strict CD Condition

abstract: In (Calc.Var.PDE 2018) Schultz generalized the work of Rajala and Sturm (Calc.Var.PDE 2014), proving that a weak non-branching condition holds in the more general setting of very strict CD spaces. Anyway, similar to what happens for the strong CD condition, the very strict CD condition seems not to be stable with respect to the measured Gromov Hausdorff convergence. In this article I prove a stability result for the very strict CD condition, assuming some metric requirements on the converging sequence and on the limit space. The proof relies on the notions of \textit{consistent geodesic flow} and \textit{consistent plan selection}, which allow to treat separately the static and the dynamic part of a Wasserstein geodesic. As an application, I prove that the metric measure space $\mathbb R^N$ equipped with a crystalline norm and with the Lebesgue measure satisfies the very strict $\ CD(0,\infty)$ condition.


title: Extension of boundary homeomorphisms to mappings of finite distortion

abstract: We provide sufficient conditions so that a homeomorphism of the real line or of the circle admits an extension to a mapping of finite distortion in the upper half-plane or the disk, respectively. Moreover, we can ensure that the quasiconformal dilatation of the extension satisfies certain integrability conditions, such as $p$-integrability or exponential integrability. Mappings satisfying the latter integrability condition are also known as David homeomorphisms. Our extension operator is the same as the one used by Beurling and Ahlfors in their celebrated work. We prove an optimal bound for the quasiconformal dilatation of the Beurling--Ahlfors extension of a homeomorphism of the real line, in terms of its symmetric distortion function. More specifically, the quasiconformal dilatation is bounded above by an average of the symmetric distortion function and below by the symmetric distortion function itself. As a consequence, the quasiconformal dilatation of the Beurling--Ahlfors extension of a homeomorphism of the real line is (sub)exponentially integrable, is $p$-integrable, or has a $BMO$ majorant if and only if the symmetric distortion is (sub)exponentially integrable, is $p$-integrable, or has a $BMO$ majorant, respectively. These theorems are all new and reconcile several sufficient extension conditions that have been established in the past.


title: Sequential prediction under log-loss and misspecification

abstract: We consider the question of sequential prediction under the log-loss in terms of cumulative regret. Namely, given a hypothesis class of distributions, learner sequentially predicts the (distribution of the) next letter in sequence and its performance is compared to the baseline of the best constant predictor from the hypothesis class. The well-specified case corresponds to an additional assumption that the data-generating distribution belongs to the hypothesis class as well. Here we present results in the more general misspecified case. Due to special properties of the log-loss, the same problem arises in the context of competitive-optimality in density estimation, and model selection. For the $d$-dimensional Gaussian location hypothesis class, we show that cumulative regrets in the well-specified and misspecified cases asymptotically coincide. In other words, we provide an $o(1)$ characterization of the distribution-free (or PAC) regret in this case -- the first such result as far as we know. We recall that the worst-case (or individual-sequence) regret in this case is larger by an additive constant ${d\over 2} + o(1)$. Surprisingly, neither the traditional Bayesian estimators, nor the Shtarkov's normalized maximum likelihood achieve the PAC regret and our estimator requires special "robustification" against heavy-tailed data. In addition, we show two general results for misspecified regret: the existence and uniqueness of the optimal estimator, and the bound sandwiching the misspecified regret between well-specified regrets with (asymptotically) close hypotheses classes.


title: Dynamics of Laterally Propagating Flames in X-ray Bursts. II. Realistic Burning & Rotation

abstract: We continue to investigate two-dimensional laterally propagating flames in type I X-ray bursts using fully compressible hydrodynamics simulations. In the current study we relax previous approximations where we artificially boosted the flames. We now use more physically realistic reaction rates, thermal conductivities, and rotation rates, exploring the effects of neutron star rotation rate and thermal structure on the flame. We find that at lower rotation rates the flame becomes harder to ignite, whereas at higher rotation rates the nuclear burning is enhanced by increased confinement from the Coriolis force and the flame propagates steadily. At higher crustal temperatures, the flame moves more quickly and accelerates as it propagates through the atmosphere. If the temperature is too high, instead of a flame propagating across the surface the entire atmosphere burns steadily. All of the software used for these simulations is freely available.


title: Phase curve and variability analysis of WASP-12b using TESS photometry

abstract: We analyse Sector 20 TESS photometry of the ultra-hot Jupiter WASP-12b, and extract its phase curve to study the planet's atmospheric properties. We successfully recover the phase curve with an amplitude of 549 $\pm$ 62 ppm, and a secondary eclipse depth of 609$^{+74}_{-73}$ ppm. The peak of the phase curve is shifted by 0.049 $\pm$ 0.015 in phase, implying that the brightest spot in the atmosphere is shifted from the substellar point towards the planet's evening terminator. Assuming zero albedo, the eclipse depth infers a day-side brightness temperature of 3128$^{+64}_{-68}$ K. No significant detection of flux from the night-side is found at 60 $\pm$ 97 ppm, implying a night-side brightness temperature of $<$2529 K (1-$\sigma$). We do not detect any significant variability in the light from the planet over the $\sim$27 days of the TESS observations. Finally, we note that an ephemeris model taking orbital decay into account provides a significantly better fit than a constant-period model.


title: Poincar\'{e}-Bendixson Limit Sets in Multi-Agent Learning

abstract: A key challenge of evolutionary game theory and multi-agent learning is to characterize the limit behavior of game dynamics. Whereas convergence is often a property of learning algorithms in games satisfying a particular reward structure (e.g., zero-sum games), even basic learning models, such as the replicator dynamics, are not guaranteed to converge for general payoffs. Worse yet, chaotic behavior is possible even in rather simple games, such as variants of the Rock-Paper-Scissors game. Although chaotic behavior in learning dynamics can be precluded by the celebrated Poincar\'e-Bendixson theorem, it is only applicable to low-dimensional settings. Are there other characteristics of a game that can force regularity in the limit sets of learning? We show that behavior consistent with the Poincar\'e-Bendixson theorem (limit cycles, but no chaotic attractor) can follow purely from the topological structure of the interaction graph, even for high-dimensional settings with an arbitrary number of players and arbitrary payoff matrices. We prove our result for a wide class of follow-the-regularized leader (FoReL) dynamics, which generalize replicator dynamics, for binary games characterized interaction graphs where the payoffs of each player are only affected by one other player (i.e., interaction graphs of indegree one). Since chaos occurs already in games with only two players and three strategies, this class of non-chaotic games may be considered maximal. Moreover, we provide simple conditions under which such behavior translates into efficiency guarantees, implying that FoReL learning achieves time-averaged sum of payoffs at least as good as that of a Nash equilibrium, thereby connecting the topology of the dynamics to social-welfare analysis.


title: Fast spectroscopic mapping of two-dimensional quantum materials

abstract: Spectroscopic mapping refers to the massive recording of spectra whilst varying an additional degree of freedom, such as: magnetic field, location, temperature, or charge carrier concentration. As this involves two serial tasks, spectroscopic mapping can become excruciatingly slow. We demonstrate exponentially faster mapping through our combination of sparse sampling and parallel spectroscopy. We exemplify our concept using quasiparticle interference imaging of Au(111) and Bi2Sr2CaCu2O8 (Bi2212), as two well-known model systems. Our method is accessible, straightforward to implement with existing scanning tunneling microscopes, and can be easily extended to enhance gate or field-mapping spectroscopy. In view of a possible four orders of magnitude speed advantage, it is setting the stage to fundamentally promote the discovery of novel quantum materials.


title: Lower Bounds on Information Requirements for Causal Network Inference

abstract: Recovery of the causal structure of dynamic networks from noisy measurements has long been a problem of intense interest across many areas of science and engineering. Many algorithms have been proposed, but there is no work that compares the performance of the algorithms to converse bounds in a non-asymptotic setting. As a step to address this problem, this paper gives lower bounds on the error probability for causal network support recovery in a linear Gaussian setting. The bounds are based on the use of the Bhattacharyya coefficient for binary hypothesis testing problems with mixture probability distributions. Comparison of the bounds and the performance achieved by two representative recovery algorithms are given for sparse random networks based on the Erd\H{o}s-R\'enyi model.


title: A Graphical Interpretation of Circumgalactic Precipitation

abstract: Both observations and recent numerical simulations of the circumgalactic medium (CGM) support the hypothesis that a self-regulating feedback loop suspends the gas density of the ambient CGM close to the galaxy in a state with a ratio of cooling time to freefall time >10. This limiting ratio is thought to arise because circumgalactic gas becomes increasingly susceptible to multiphase condensation as the ratio declines. If the timescale ratio gets too small, then cold clouds precipitate out of the CGM, rain into the galaxy, and fuel energetic feedback that raises the ambient cooling time. The astrophysical origin of this so-called precipitation limit is not simple but is critical to understanding the CGM and its role in galaxy evolution. This paper therefore attempts to interpret its origin as simply as possible, relying mainly on conceptual reasoning and schematic diagrams. It illustrates how the precipitation limit can depend on both the global configuration of a galactic atmosphere and the degree to which dynamical disturbances drive CGM perturbations. It also frames some tests of the precipitation hypothesis that can be applied to both CGM observations and numerical simulations of galaxy evolution.


title: An Information-Theoretic Law Governing Human Multi-Task Navigation Decisions

abstract: To better understand the process by which humans make navigation decisions when tasked with multiple stopovers, we analyze motion data captured from shoppers in a grocery store. We discover several trends in the data that are consistent with a noisy decision making process for the order of item retrieval, and decompose a shopping trip into a sequence of discrete choices about the next item to retrieve. Our analysis reveals that the likelihood of inverting any two items in the order is monotonically bound to the entropy of the pair-wise ordering task. Based on this analysis, we propose a noisy distance estimation model for predicting the order of item retrieval given a shopping list. We show that our model theoretically reproduces the entropy law seen in the data with high accuracy, and in practice matches the trends in the data when used to simulate the same shopping lists. Our approach has direct applications to improving simulations of human navigation in retail and other settings.


title: Statistical Inference after Kernel Ridge Regression Imputation under item nonresponse

abstract: Imputation is a popular technique for handling missing data. We consider a nonparametric approach to imputation using the kernel ridge regression technique and propose consistent variance estimation. The proposed variance estimator is based on a linearization approach which employs the entropy method to estimate the density ratio. The root-n consistency of the imputation estimator is established when a Sobolev space is utilized in the kernel ridge regression imputation, which enables us to develop the proposed variance estimator. Synthetic data experiments are presented to confirm our theory.


title: Debt Representation in UTXO Blockchains

abstract: We provide a UTXO model of blockchain transactions that is able to represent both credit and debt on the same blockchain. Ordinarily, the UTXO model is solely used to represent credit and the representation of credit and debit together is achieved using the account model because of its support for balances. However, the UTXO model provides superior privacy, safety, and scalability when compared to the account model. In this work, we introduce a UTXO model that has the flexibility of balances with the usual benefits of the UTXO model. This model extends the conventional UTXO model, which represents credits as unmatched outputs, by representing debts as unmatched inputs. We apply our model to solving the problem of transparency in reverse mortgage markets, in which some transparency is necessary for a healthy market but complete transparency leads to adverse outcomes. Here the pseudonymous properties of the UTXO model protect the privacy of loan recipients while still allowing an aggregate view of the loan market. We present a prototype of our implementation in Tendermint and discuss the design and its benefits.


title: Modifications to the signal from a gravitational wave event due to a surrounding shell of matter

abstract: In previous work, we established theoretical results concerning the effect of matter shells surrounding a gravitational wave (GW) source, and we now apply these results to astrophysical scenarios. Firstly, it is shown that GW echoes that are claimed to be present in LIGO data of certain events, could not have been caused by a matter shell. However, it is also shown that there are scenarios in which matter shells could make modifications of order a few percent to a GW signal; these scenarios include binary black hole mergers, binary neutron star mergers, and core collapse supernovae.


title: Evidence for a high-energy tail in the gamma-ray spectra of globular clusters

abstract: Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes--curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars--have long been unclear. To address this, we search for evidence of inverse Compton emission in 8-year $\textit{Fermi}$-LAT data from the directions of 157 Milky Way globular clusters. We find a mildly statistically significant (3.8$\sigma$) correlation between the measured globular cluster gamma-ray luminosities and their photon field energy densities. However, this may also be explained by a hidden correlation between the photon field densities and the stellar encounter rates of globular clusters. Analysed $\textit{in toto}$, we demonstrate that the gamma-ray emission of globular clusters can be resolved spectrally into two components: i) an exponentially cut-off power law and ii) a pure power law. The latter component--which we uncover at a significance of 8.2$\sigma$--has a power index of 2.79 $\pm$ 0.25. It is most naturally interpreted as inverse Compton emission by cosmic-ray electrons and positrons injected by millisecond pulsars. We find the luminosity of this power-law component is comparable to, or slightly smaller than, the luminosity of the curved component, suggesting the fraction of millisecond pulsar spin-down luminosity into relativistic leptons is similar to the fraction of the spin-down luminosity into prompt magnetospheric radiation.


title: Neural 3D Clothes Retargeting from a Single Image

abstract: In this paper, we present a method of clothes retargeting; generating the potential poses and deformations of a given 3D clothing template model to fit onto a person in a single RGB image. The problem is fundamentally ill-posed as attaining the ground truth data is impossible, i.e., images of people wearing the different 3D clothing template model at exact same pose. We address this challenge by utilizing large-scale synthetic data generated from physical simulation, allowing us to map 2D dense body pose to 3D clothing deformation. With the simulated data, we propose a semi-supervised learning framework that validates the physical plausibility of the 3D deformation by matching with the prescribed body-to-cloth contact points and clothing silhouette to fit onto the unlabeled real images. A new neural clothes retargeting network (CRNet) is designed to integrate the semi-supervised retargeting task in an end-to-end fashion. In our evaluation, we show that our method can predict the realistic 3D pose and deformation field needed for retargeting clothes models in real-world examples.


title: Recurrent Localization Networks applied to the Lippmann-Schwinger Equation

abstract: The bulk of computational approaches for modeling physical systems in materials science derive from either analytical (i.e. physics based) or data-driven (i.e. machine-learning based) origins. In order to combine the strengths of these two approaches, we advance a novel machine learning approach for solving equations of the generalized Lippmann-Schwinger (L-S) type. In this paradigm, a given problem is converted into an equivalent L-S equation and solved as an optimization problem, where the optimization procedure is calibrated to the problem at hand. As part of a learning-based loop unrolling, we use a recurrent convolutional neural network to iteratively solve the governing equations for a field of interest. This architecture leverages the generalizability and computational efficiency of machine learning approaches, but also permits a physics-based interpretation. We demonstrate our learning approach on the two-phase elastic localization problem, where it achieves excellent accuracy on the predictions of the local (i.e., voxel-level) elastic strains. Since numerous governing equations can be converted into an equivalent L-S form, the proposed architecture has potential applications across a range of multiscale materials phenomena.

categories: physics.comp-ph cond-mat.dis-nn cond-mat.mtrl-sci cs.LG


title: Choquet-Sugeno-like operator based on relation and conditional aggregation operators

abstract: We introduce a~\textit{Choquet-Sugeno-like operator} generalizing many operators for bounded functions and monotone measures from the literature, e.g., Sugeno-like operator, Lov\'{a}sz and Owen measure extensions, $\rF$-decomposition integral with respect to a~partition decomposition system, and others. The new operator is based on the concepts of dependence relation and conditional aggregation operators, but it does not depend on $t$-level sets. We also provide conditions for which the Choquet-Sugeno-like operator coincides with some Choquet-like integrals defined on finite spaces and appeared recently in the literature, e.g. reverse Choquet integral, $d$-Choquet integral, $\rF$-based discrete Choquet-like integral, some version of $C_{\rF_1\rF_2}$-integral, $\mathrm{C}\mathrm{C}$-integrals (or Choquet-like Copula-based integral) and discrete inclusion-exclusion integral. Some basic properties of the Choquet-Sugeno-like operator are studied.


title: Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data

abstract: Radio interferometric (RI) data are noisy under-sampled spatial Fourier components of the unknown radio sky affected by direction-dependent antenna gains. Failure to model these antenna gains accurately results in a radio sky estimate with limited fidelity and resolution. The RI inverse problem has been recently addressed via a joint calibration and imaging approach which consists in solving a non-convex minimisation task, involving suitable priors for the DDEs, namely temporal and spatial smoothness, and sparsity for the unknown radio map via an $\ell_1$-norm prior, in the context of realistic RI simulations. Building on these developments, we propose to promote sparsity of the radio map via a log-sum prior, enforcing sparsity more strongly than the $\ell_1$-norm. The resulting minimisation task is addressed via a sequence of non-convex minimisation tasks composed of re-weighted $\ell_1$ image priors, which are solved approximately. We demonstrate the efficiency of the approach on RI observations of the celebrated radio galaxy Cygnus~A obtained with the Karl G. Jansky Very Large Array at X, C, and S bands. More precisely, we showcase that the approach enhances data fidelity significantly while achieving high resolution high dynamic range radio maps, confirming the suitability of the priors considered for the unknown DDEs and radio image. As a clear qualitative indication of the high fidelity achieved by the data and the proposed approach, we report the detection of three background sources in the vicinity of Cyg~A, at S band.


title: R-process-rich stellar streams in the Milky Way

abstract: We present high-resolution Magellan/MIKE spectra of 22 bright ($9<V<13.5$) metal-poor stars ($-3.18<\mbox{[Fe/H]}<-1.37$) in three different stellar streams, the Helmi debris stream, the Helmi trail stream, and the $\omega$ Centauri progenitor stream. We augment our Helmi debris sample with results for ten stars by Roederer et al. 2010 (arXiv:1001.1745), for a total of 32 stars. Detailed chemical abundances of light elements as well as heavy neutron-capture elements have been determined for our 22 stars. All three streams contain carbon-enhanced stars. For 13 stars, neutron-capture element lines were detectable and they all show signatures in agreement with the scaled solar $r$-process pattern, albeit with a large spread of $-0.5<\mbox{[Eu/Fe]}<+1.3$. Eight of these stars show an additional small $s$-process contribution superposed onto their $r$-process pattern. This could be discerned because of the relatively high $S/N$ of the spectra given that the stars are close by in the halo. Our results suggest that the progenitors of these streams experienced one or more $r$-process events, such as a neutron star merger or another prolific $r$-process source, early on that widely enriched these host systems before their accretion by the Milky Way. The small $s$-process contribution suggests the presence of AGB stars and associated local (inhomogeneous) enrichment as part of the ongoing chemical evolution by low mass stars. Stars in stellar streams may thus be a promising avenue for studying the detailed history of large dwarf galaxies and their role in halo assembly with easily accessible targets for high-quality spectra of many stars.


title: Multi-Block Sparse Functional Principal Components Analysis for Longitudinal Microbiome Multi-Omics Data

abstract: Microbiome researchers often need to model the temporal dynamics of multiple complex, nonlinear outcome trajectories simultaneously. This motivates our development of multivariate Sparse Functional Principal Components Analysis (mSFPCA), extending existing SFPCA methods to simultaneously characterize multiple temporal trajectories and their inter-relationships. As with existing SFPCA methods, the mSFPCA algorithm characterizes each trajectory as a smooth mean plus a weighted combination of the smooth major modes of variation about the mean, where the weights are given by the component scores for each subject. Unlike existing SFPCA methods, the mSFPCA algorithm allows estimation of multiple trajectories simultaneously, such that the component scores, which are constrained to be independent within a particular outcome for identifiability, may be arbitrarily correlated with component scores for other outcomes. A Cholesky decomposition is used to estimate the component score covariance matrix efficiently and guarantee positive semi-definiteness given these constraints. Mutual information is used to assess the strength of marginal and conditional temporal associations across outcome trajectories. Importantly, we implement mSFPCA as a Bayesian algorithm using R and stan, enabling easy use of packages such as PSIS-LOO for model selection and graphical posterior predictive checks to assess the validity of mSFPCA models. Although we focus on application of mSFPCA to microbiome data in this paper, the mSFPCA model is of general utility and can be used in a wide range of real-world applications.


title: The central parsec of NGC 3783: a rotating broad emission line region, asymmetric hot dust structure, and compact coronal line region

abstract: Using VLTI/GRAVITY and SINFONI data, we investigate the sub-pc gas and dust structure around the nearby type 1 AGN hosted by NGC 3783. The K-band coverage of GRAVITY uniquely allows a simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-IR continuum emitting hot dust, and the size of the coronal line region (CLR). We find the BLR probed through broad Br$\gamma$ emission is well described by a rotating, thick disk with a radial distribution of clouds peaking in the inner region. In our BLR model the physical mean radius of 16 light days is nearly twice the 10 day time lag that would be measured, which matches very well the 10 day time lag that has been measured by reverberation mapping. We measure a hot dust FWHM size of 0.74 mas (0.14 pc) and further reconstruct an image of the hot dust which reveals a faint (5% of the total flux) offset cloud which we interpret as an accreting cloud heated by the central AGN. Finally, we directly measure the FWHM size of the nuclear CLR as traced by the [CaVIII] and narrow Br$\gamma$ line. We find a FWHM size of 2.2 mas (0.4 pc), fully in line with the expectation of the CLR located between the BLR and narrow line region. Combining all of these measurements together with larger scale near-IR integral field unit and mid-IR interferometry data, we are able to comprehensively map the structure and dynamics of gas and dust from 0.01--100 pc.


title: Bochner-Riesz Means Convergence of Prolate Spheroidal Series and Their Extensions

abstract: In this paper, we study the $L^p$-Bochner-Riesz mean summability problem related to the spectrum of some particular Sturm-Liouville operators in the weighted $L^p([a,b],\omega).$ Our purpose is to establish suitable conditions under which the Bochner-Riesz expansion of a function $f \in L^p([a,b],\omega)$,$1<p<\infty$, in two generalisations of Slepian's basis, converges to $f$ in $L^p([a,b],\omega)$.


title: Network-centric indicators for fragility in global financial indices

abstract: Over the last two decades, financial systems have been studied and analysed from the perspective of complex networks, where the nodes and edges in the network represent the various financial components and the strengths of correlations between them. Here, we adopt a similar network-based approach to analyse the daily closing prices of 69 global financial market indices across 65 countries over a period of 2000-2014. We study the correlations among the indices by constructing threshold networks superimposed over minimum spanning trees at different time frames. We investigate the effect of critical events in financial markets (crashes and bubbles) on the interactions among the indices by performing both static and dynamic analyses of the correlations. We compare and contrast the structures of these networks during periods of crashes and bubbles, with respect to the normal periods in the market. In addition, we study the temporal evolution of traditional market indicators, various global network measures and the recently developed edge-based curvature measures. We show that network-centric measures can be extremely useful in monitoring the fragility in the global financial market indices.


title: Ampere-class Bright Field Emission Cathode Operated at 100 MV/m

abstract: High current bright sources are needed to power the next generation of compact rf and microwave systems. A major requirement is that such a source could be sustainably operated at high frequencies, well above 1 GHz, and high gradients, well above 100 MV/m. Field emission sources offer simplicity and scalability in a high frequency era of the injector design, but the output rf cycle charge and high gradient operation remain a great and largely unaddressed challenge. Here, a field emission cathode based on ultra-nano-crystalline diamond or UNCD, an efficient planar field emission material, was tested at 100 MV/m in an L-band injector. A very high charge of 38 pC per rf cycle (300 nC per rf pulse corresponding to rf pulse current of 120 mA) was demonstrated. This operating condition revealed a space charge dominated emission from the cathode and revealed a condition under which the 1D Child Langmuir limit was surpassed. Finally, a beam brightness of $\sim 10^{14}-10^{15}$ A/m$^2\times$rad$^2$ was estimated.


title: Strong Anharmonic and Quantum Effects in Pm-3n AlH3 Under High Pressure: A First-Principles Study

abstract: Motivated by the absence of experimental superconductivity in the metallic Pm-3n phase of AlH3 despite the predictions, we reanalyze its vibrational and superconducting properties at pressures above 99 GPa making use of first-principles techniques. In our calculations based on the self-consistent harmonic approximation method that treats anharmonicity beyond perturbation theory, we predict a strong anharmonic correction to the phonon spectra and demonstrate that the superconducting critical temperatures predicted in previous calculations based on the harmonic approximation are strongly suppressed by anharmonicity. The electron-phonon coupling concentrates on the lowest-energy hydrogen-character optical modes at the X point of the Brillouin zone. As a consequence of the strong anharmonic enhancement of their frequency, the electron-phonon coupling is suppressed by at least a 30%. The suppression in {\lambda} makes Tc smaller than 4.2 K above 120 GPa, which is well consistent with the experimental evidence. Our results underline that metal hydrides with hydrogen atoms in interstitial sites are subject to huge anharmonic effects.


title: Moving Pearl vortices in thin-film superconductors

abstract: The magnetic field $h_z$ of a moving Pearl vortex in a superconducting thin-film in $(x,y)$ plane is studied with the help of time-dependent London equation. It is found that for a vortex at the origin moving in $+x$ direction, $h_z(x,y)$ is suppressed in front of the vortex, $x>0$, and enhanced behind ($x<0$). The distribution asymmetry is proportional to the velocity and to the conductivity of normal quasiparticles. The vortex self-energy and the interaction of two moving vortices are evaluated.


title: Water Production Rate of C/2020 F3 (NEOWISE) from SOHO/SWAN over Its Active Apparition

abstract: C/2020 F3 (NEOWISE) was discovered in images from the Near Earth Object program of the Wide-Field Infrared Survey Explorer (NEOWISE) taken on 27 March 2020 and has become the Great Comet of 2020. The Solar Wind ANisotropies (SWAN) camera on the Solar and Heliospheric Observatory (SOHO) spacecraft, located in a halo orbit around the Earth-Sun L1 Lagrange point, makes daily full-sky images of hydrogen Lyman-alpha. Water production rates were determined from the SWAN hydrogen Lyman-alpha brightness and spatial distribution of the comet measured over a 4-month period of time on either side of the comet's perihelion on 3 July 2020. The water production rate in s^-1 was moderately asymmetric around perihelion and varied with the heliocentric distance, r, in au as (6.9+/-0.5) x 10^28 r^-2.5+/-0.2 and (10.1+/-0.5) x 10^28 r^-3.5+/-0.1 before and after perihelion, respectively. This is consistent with the comet having been through the planetary region of the solar system on one or more previous apparitions. A water production rates as large as 5.27 x 10^30 s^-1 were determined shortly after perihelion, once the comet was outside the solar avoidance area of SWAN, when the comet was 0.324 au from the Sun.


title: RecSSD: Near Data Processing for Solid State Drive Based Recommendation Inference

abstract: Neural personalized recommendation models are used across a wide variety of datacenter applications including search, social media, and entertainment. State-of-the-art models comprise large embedding tables that have billions of parameters requiring large memory capacities. Unfortunately, large and fast DRAM-based memories levy high infrastructure costs. Conventional SSD-based storage solutions offer an order of magnitude larger capacity, but have worse read latency and bandwidth, degrading inference performance. RecSSD is a near data processing based SSD memory system customized for neural recommendation inference that reduces end-to-end model inference latency by 2X compared to using COTS SSDs across eight industry-representative models.


title: Scalable creation of silicon-vacancy color centers in diamond by ion implantation through a 1-$\mu$m pinhole

abstract: The controlled creation of quantum emitters in diamond represents a major research effort in the fabrication of single-photon devices. Here, we present the scalable production of silicon-vacancy (SiV) color centers in single-crystal diamond by ion implantation. The lateral position of the SiV is spatially controlled by a 1-$\mu$m pinhole placed in front of the sample, which can be moved nanometer precise using a piezo stage. The initial implantation position is controlled by monitoring the ion beam position with a camera. Hereby, silicon ions are implanted at the desired spots in an area comparable to the diffraction limit. We discuss the role of ions scattered by the pinhole and the activation yield of the SiV color centers for the creation of single quantum emitters.


title: Scalable Voltage Control using Structure-Driven Hierarchical Deep Reinforcement Learning

abstract: This paper presents a novel hierarchical deep reinforcement learning (DRL) based design for the voltage control of power grids. DRL agents are trained for fast, and adaptive selection of control actions such that the voltage recovery criterion can be met following disturbances. Existing voltage control techniques suffer from the issues of speed of operation, optimal coordination between different locations, and scalability. We exploit the area-wise division structure of the power system to propose a hierarchical DRL design that can be scaled to the larger grid models. We employ an enhanced augmented random search algorithm that is tailored for the voltage control problem in a two-level architecture. We train area-wise decentralized RL agents to compute lower-level policies for the individual areas, and concurrently train a higher-level DRL agent that uses the updates of the lower-level policies to efficiently coordinate the control actions taken by the lower-level agents. Numerical experiments on the IEEE benchmark 39-bus model with 3 areas demonstrate the advantages and various intricacies of the proposed hierarchical approach.


title: Prospects of gravitational-waves detections from common-envelope evolution with LISA

abstract: Understanding common envelope (CE) evolution is an outstanding problem in binary evolution. Although the CE phase is not driven by gravitational-wave (GW) emission, the in-spiraling binary emits GWs that passively trace the CE dynamics. Detecting this GW signal would provide direct insight into the gas-driven physics. Even a non-detection might offer invaluable constraints. We investigate the prospects of detection of a Galactic CE by LISA. While the dynamical phase of the CE is likely sufficiently loud for detection, it is short and thus rare. We focus instead on the self-regulated phase that proceeds on a thermal timescale. Based on population synthesis calculations and the (unknown) signal duration in the LISA band, we expect $\sim 0.1-100$ sources in the Galaxy during the mission duration. We map the GW observable parameter space of frequency $f_\mathrm{GW}$ and its derivative $\dot f_\mathrm{GW}$ remaining agnostic on the specifics of the inspiral, and find that signals with $\mathrm{SNR}>10$ are possible if the CE stalls at separations such that $f_\mathrm{GW}\gtrsim2\times10^{-3}\,\mathrm{Hz}$. We investigate the possibility of misidentifying the signal with other known sources. If the second derivative $\ddot f_\mathrm{GW}$ can also be measured, the signal can be distinguished from other sources using a GW braking-index. Alternatively, coupling LISA with electromagnetic observations of peculiar red giant stars and/or infrared and optical transients might allow for the disentangling of a Galactic CE from other Galactic and extra-galactic GW sources.


title: CAMBI: Contrast-aware Multiscale Banding Index

abstract: Banding artifacts are artificially-introduced contours arising from the quantization of a smooth region in a video. Despite the advent of recent higher quality video systems with more efficient codecs, these artifacts remain conspicuous, especially on larger displays. In this work, a comprehensive subjective study is performed to understand the dependence of the banding visibility on encoding parameters and dithering. We subsequently develop a simple and intuitive no-reference banding index called CAMBI (Contrast-aware Multiscale Banding Index) which uses insights from Contrast Sensitivity Function in the Human Visual System to predict banding visibility. CAMBI correlates well with subjective perception of banding while using only a few visually-motivated hyperparameters.


title: The sign problem in density matrix quantum Monte Carlo

abstract: Density matrix quantum Monte Carlo (DMQMC) is a recently-developed method for stochastically sampling the $N$-particle thermal density matrix to obtain exact-on-average energies for model and \emph{ab initio} systems. We report a systematic numerical study of the sign problem in DMQMC based on simulations of atomic and molecular systems. In DMQMC, the density matrix is written in an outer product basis of Slater determinants and has a size of space which is the square of the number of Slater determinants. In principle this means DMQMC needs to sample a space which scales in the system size, $N$, as $\mathcal{O}[(\exp(N))^2]$. In practice, there is a system-dependent critical walker population ($N_c$) which must be exceeded in order to remove the sign problem, and this imposes limitations by way of storage and computer time. We establish that $N_c$ for DMQMC is the square of $N_c$ for FCIQMC. By contrast, the minimum $N_c$ in the interaction picture modification of DMQMC (IP-DMQMC) only is directly proportionate to the $N_c$ for FCIQMC. We find that this comes from the asymmetric propagation of IP-DMQMC compared to the symmetric propagation of canonical DMQMC. An asymmetric mode of propagation is prohibitively expensive for DMQMC because it has a much greater stochastic error. Finally, we find that the equivalence between IP-DMQMC and FCIQMC seems to extend to the initiator approximation, which is often required to study larger basis sets and other systems. This suggests IP-DMQMC offers a way to ameliorate the cost of moving between a Slater determinant space and an outer product basis.


title: Phase shifting, dispersion variation and defocusing suppression in wave breaking

abstract: We present an investigation of the fundamental physical processes involved in deep water wave breaking. Our motivation is to identify the underlying reason causing the deficiency of the eddy viscosity breaking model (EVBM) in predicting surface elevation for strongly nonlinear waves. Owing to the limitation of experimental methods in the provision of high-resolution flow information, we propose a numerical methodology by developing an EVBM enclosed standalone fully-nonlinear quasi-potential (FNP) flow model and a coupled FNP plus Navier-Stokes flow model. The numerical models were firstly verified with a wave train subject to modulational instability, then used to simulate a series of broad-banded focusing wave trains under non-, moderate- and strong-breaking conditions. A systematic analysis was carried out to investigate the discrepancies of numerical solutions produced by the two models in surface elevation and other important physical properties. It is found that EVBM predicts accurately the energy dissipated by breaking and the amplitude spectrum of free waves in terms of magnitude, but fails to capture accurately breaking induced phase shifting. The shift of phase grows with breaking intensity and is especially strong for high wavenumber components. This is identified as a cause of the upshift of wave dispersion relation, which increases the frequencies of large wavenumber components. Such a variation drives large-wavenumber components to propagate at nearly the same speed, which is significantly higher than the linear dispersion levels. This suppresses the instant dispersive spreading of harmonics after the focal point, prolonging the lifespan of focused waves and expanding their propagation space.


title: Settling the Sharp Reconstruction Thresholds of Random Graph Matching

abstract: This paper studies the problem of recovering the hidden vertex correspondence between two edge-correlated random graphs. We focus on the Gaussian model where the two graphs are complete graphs with correlated Gaussian weights and the Erd\H{o}s-R\'enyi model where the two graphs are subsampled from a common parent Erd\H{o}s-R\'enyi graph $\mathcal{G}(n,p)$. For dense graphs with $p=n^{-o(1)}$, we prove that there exists a sharp threshold, above which one can correctly match all but a vanishing fraction of vertices and below which correctly matching any positive fraction is impossible, a phenomenon known as the "all-or-nothing" phase transition. Even more strikingly, in the Gaussian setting, above the threshold all vertices can be exactly matched with high probability. In contrast, for sparse Erd\H{o}s-R\'enyi graphs with $p=n^{-\Theta(1)}$, we show that the all-or-nothing phenomenon no longer holds and we determine the thresholds up to a constant factor. Along the way, we also derive the sharp threshold for exact recovery, sharpening the existing results in Erd\H{o}s-R\'enyi graphs. The proof of the negative results builds upon a tight characterization of the mutual information based on the truncated second-moment computation and an "area theorem" that relates the mutual information to the integral of the reconstruction error. The positive results follows from a tight analysis of the maximum likelihood estimator that takes into account the cycle structure of the induced permutation on the edges.


title: Reduced operator inference for nonlinear partial differential equations

abstract: We present a new scientific machine learning method that learns from data a computationally inexpensive surrogate model for predicting the evolution of a system governed by a time-dependent nonlinear partial differential equation (PDE), an enabling technology for many computational algorithms used in engineering settings. Our formulation generalizes to the function space PDE setting the Operator Inference method previously developed in [B. Peherstorfer and K. Willcox, Data-driven operator inference for non-intrusive projection-based model reduction, Computer Methods in Applied Mechanics and Engineering, 306 (2016)] for systems governed by ordinary differential equations. The method brings together two main elements. First, ideas from projection-based model reduction are used to explicitly parametrize the learned model by low-dimensional polynomial operators which reflect the known form of the governing PDE. Second, supervised machine learning tools are used to infer from data the reduced operators of this physics-informed parametrization. For systems whose governing PDEs contain more general (non-polynomial) nonlinearities, the learned model performance can be improved through the use of lifting variable transformations, which expose polynomial structure in the PDE. The proposed method is demonstrated on two examples: a heat equation model problem that demonstrates the benefits of the function space formulation in terms of consistency with the underlying continuous truth, and a three-dimensional combustion simulation with over 18 million degrees of freedom, for which the learned reduced models achieve accurate predictions with a dimension reduction of five orders of magnitude and model runtime reduction of up to nine orders of magnitude.


title: A linearized framework and a new benchmark for model selection for fine-tuning

abstract: Fine-tuning from a collection of models pre-trained on different domains (a "model zoo") is emerging as a technique to improve test accuracy in the low-data regime. However, model selection, i.e. how to pre-select the right model to fine-tune from a model zoo without performing any training, remains an open topic. We use a linearized framework to approximate fine-tuning, and introduce two new baselines for model selection -- Label-Gradient and Label-Feature Correlation. Since all model selection algorithms in the literature have been tested on different use-cases and never compared directly, we introduce a new comprehensive benchmark for model selection comprising of: i) A model zoo of single and multi-domain models, and ii) Many target tasks. Our benchmark highlights accuracy gain with model zoo compared to fine-tuning Imagenet models. We show our model selection baseline can select optimal models to fine-tune in few selections and has the highest ranking correlation to fine-tuning accuracy compared to existing algorithms.


title: Long-Range Seasonal Forecasting of 2m-Temperature with Machine Learning

abstract: A significant challenge in seasonal climate prediction is whether a prediction can beat climatology. We hereby present results from two data-driven models - a convolutional (CNN) and a recurrent (RNN) neural network - that predict 2 m temperature out to 52 weeks for six geographically-diverse locations. The motivation for testing the two classes of ML models is to allow the CNN to leverage information related to teleconnections and the RNN to leverage long-term historical temporal signals. The ML models boast improved accuracy of long-range temperature forecasts up to a lead time of 30 weeks for PCC and up 52 weeks for RMSESS, however only for select locations. Further iteration is required to ensure the ML models have value beyond regions where the climatology has a noticeably reduced correlation skill, namely the tropics.


title: Challenges in Automated Debiasing for Toxic Language Detection

abstract: Biased associations have been a challenge in the development of classifiers for detecting toxic language, hindering both fairness and accuracy. As potential solutions, we investigate recently introduced debiasing methods for text classification datasets and models, as applied to toxic language detection. Our focus is on lexical (e.g., swear words, slurs, identity mentions) and dialectal markers (specifically African American English). Our comprehensive experiments establish that existing methods are limited in their ability to prevent biased behavior in current toxicity detectors. We then propose an automatic, dialect-aware data correction method, as a proof-of-concept. Despite the use of synthetic labels, this method reduces dialectal associations with toxicity. Overall, our findings show that debiasing a model trained on biased toxic language data is not as effective as simply relabeling the data to remove existing biases.


title: Are Work Zones and Connected Automated Vehicles Ready for a Harmonious Coexistence? A Scoping Review and Research Agenda

abstract: The recent advent of connected and automated vehicles (CAVs) is expected to transform the transportation system. CAV technologies are being developed rapidly and they are foreseen to penetrate the market at a rapid pace. On the other hand, work zones (WZs) have become common areas on highway systems as a result of the increasing construction and maintenance activities. The near future will therefore bring the coexistence of CAVs and WZs which makes their interaction inevitable. WZs expose all vehicles to a sudden and complex geometric change in the roadway environment, something that may challenge many of CAV navigation capabilities. WZs however also impose a space contraction resulting in adverse traffic impacts, something that legitimately calls for benefiting from the highly efficient CAV functions. CAVs should be able to reliably traverse WZ geometry and WZs should benefit from CAV intelligent functions. This paper reviews the state-of-the-art and the key concepts, opportunities, and challenges of deploying CAV systems at WZs. The reviewed subjects include traffic performance and behaviour, technologies and infrastructure, and regulatory considerations. Eighteen CAV mobility, safety, and environmental concepts and functions were distributed over the WZ area which was subdivided into five segments: further upstream, approach area, queuing area, WZ activity, and termination area. In addition, among other topics reviewed and discussed are detection of WZ features, smart traffic control devices, various technologies at connected WZs, cross-border harmonization, liability, insurance, and privacy. The paper also provides a research agenda with a list of research needs supported by experts rating and inputs. The paper aims to provide a bird eye view, but with necessary details that can benefit researchers, practitioners, and transportation agencies.