title: Cosmic ray streaming in the turbulent interstellar medium

abstract: We study the streaming instability of GeV$-100~$GeV cosmic rays (CRs) and its damping in the turbulent interstellar medium (ISM). We find that the damping of streaming instability is dominated by ion-neutral collisional damping in weakly ionized molecular clouds, turbulent damping in the highly ionized warm medium, and nonlinear Landau damping in the Galactic halo. Only in the Galactic halo, is the streaming speed of CRs close to the Alfv\'{e}n speed. Alfv\'{e}nic turbulence plays an important role in both suppressing the streaming instability and regulating the diffusion of streaming CRs via magnetic field line tangling, with the effective mean free path of streaming CRs in the observer frame determined by the Alfv\'{e}nic scale in super-Alfv\'{e}nic turbulence. The resulting diffusion coefficient is sensitive to Alfv\'{e}n Mach number, which has a large range of values in the multi-phase ISM. Super-Alfv\'{e}nic turbulence contributes to additional confinement of streaming CRs, irrespective of the dominant damping mechanism.


title: The fate of Parisi-Sourlas supersymmetry in Random Field models

abstract: By the Parisi-Sourlas conjecture, the critical point of a theory with random field (RF) disorder is described by a supersymmeric (SUSY) conformal field theory (CFT), related to a $d-2$ dimensional CFT without SUSY. Numerical studies indicate that this is true for the RF $\phi^3$ model but not for RF $\phi^4$ model in $d<5$ dimensions. Here we argue that the SUSY fixed point is not reached because of new relevant SUSY-breaking interactions. We use perturbative renormalization group in a judiciously chosen field basis, allowing systematic exploration of the space of interactions. Our computations agree with the numerical results for both cubic and quartic potential.


title: Probing ultra-light axions with the 21-cm Signal during Cosmic Dawn

abstract: Ultra-light axions (ULAs) are a promising and intriguing set of dark-matter candidates. We study the prospects to use forthcoming measurements of 21-cm fluctuations from cosmic dawn to probe ULAs. We focus in particular on the velocity acoustic oscillations (VAOs) in the large-scale 21-cm power spectrum, features imprinted by the long-wavelength ($k\sim0.1\,{\rm Mpc}^{-1}$) modulation, by dark-matter--baryon relative velocities, of the small-scale ($k\sim 10-10^3\, {\rm Mpc}^{-1}$) power required to produce the stars that heat the neutral hydrogen. Damping of small-scale power by ULAs reduces the star-formation rate at cosmic dawn which then leads to a reduced VAO amplitude. Accounting for different assumptions for feedback and foregrounds, experiments like HERA may be sensitive to ULAs with masses up to $m_{\alpha}\approx 10^{-18}\text{eV}$, two decades of mass higher than current constraints.


title: Dynamical Friction, Buoyancy and Core-Stalling -- I. A Non-perturbative Orbit-based Analysis

abstract: We examine the origin of dynamical friction using a non-perturbative, orbit-based approach. Unlike the standard perturbative approach, in which dynamical friction arises from the LBK torque due to pure resonances, this alternative, complementary view nicely illustrates how a massive perturber significantly changes the energies and angular momenta of field particles on near-resonant orbits, with friction arising from an imbalance between particles that gain energy and those that lose energy. We treat dynamical friction in a spherical host system as a restricted three-body problem. This treatment is applicable in the `slow' regime, in which the perturber sinks slowly and the standard perturbative framework fails due to the onset of non-linearities. Hence it is especially suited to investigate the origin of core-stalling: the cessation of dynamical friction in central constant-density cores. We identify three different families of near-co-rotation-resonant orbits that dominate the contribution to dynamical friction. Their relative contribution is governed by the Lagrange points (fixed points in the co-rotating frame). In particular, one of the three families, which we call Pac-Man orbits because of their appearance in the co-rotating frame, is unique to cored density distributions. When the perturber reaches a central core, a bifurcation of the Lagrange points drastically changes the orbital make-up, with Pac-Man orbits becoming dominant. In addition, due to relatively small gradients in the distribution function inside a core, the net torque from these Pac-Man orbits becomes positive (enhancing), thereby effectuating a dynamical buoyancy. We argue that core stalling occurs where this buoyancy is balanced by friction.

categories: astro-ph.GA astro-ph.EP nlin.CD physics.class-ph


title: Analytic Modeling of Synchrotron-Self-Compton Spectra: Application to GRB 190114C

abstract: Observations of TeV emission from early gamma-ray burst (GRB) afterglows revealed the long sought for inverse Compton (IC) upscattering of the lower energy synchrotron. However, it turned out that the long hoped for ability to easily interpret the synchrotron-self-Compton (SSC) spectra didn't materialize. The TeV emission is in the Klein-Nishina (KN) regime and the simple Thomson regime SSC spectrum is modified, complicating the scene. We describe here a methodology to obtain an analytic approximation to an observed spectrum and infer the conditions at the emitting region. The methodology is general and can be used in any such source. As a test case, we apply it to the observations of GRB 190114C. We find that the procedure of fitting the model parameters using the analytic SSC spectrum suffers from some generic problems. However, at the same time, it conveniently gives a useful insight into the conditions that shape the spectrum. Once we introduce a correction to the standard KN approximation, the best fit solution is consistent with the one found in detailed numerical simulations. As in the numerical analysis, we find a family of solutions that provide a good approximation to the data and satisfy roughly $B\propto \Gamma^{-3}$ between the magnetic field and the bulk Lorentz factor, and we provide a tentative explanation why such a family arises.


title: Non-zero momentum requires long-range entanglement

abstract: We show that a quantum state in a lattice spin (boson) system must be long-range entangled if it has non-zero lattice momentum, i.e. if it is an eigenstate of the translation symmetry with eigenvalue $e^{iP}\neq1$. Equivalently, any state that can be connected with a non-zero momentum state through a finite-depth local unitary transformation must also be long-range entangled. The statement can also be generalized to fermion systems. Some non-trivial consequences follow immediately from our theorem: (1) several different types of Lieb-Schultz-Mattis-Oshikawa-Hastings (LSMOH) theorems, including a previously unknown version involving only a discrete $\mathbb{Z}_n$ symmetry, can be derived in a simple manner from our result; (2) a gapped topological order (in space dimension $d>1$) must weakly break translation symmetry if one of its ground states on torus has nontrivial momentum - this generalizes the familiar physics of Tao-Thouless; (3) our result provides further evidence of the "smoothness" assumption widely used in the classification of crystalline symmetry-protected topological (cSPT) phases.

categories: cond-mat.str-el cond-mat.mes-hall hep-th quant-ph


title: HEALPix Alchemy: Fast All-Sky Geometry and Image Arithmetic in a Relational Database for Multimessenger Astronomy Brokers

abstract: Efficient searches for electromagnetic counterparts to gravitational wave, high-energy neutrino, and gamma-ray burst events demand rapid processing of image arithmetic and geometry set operations in a database to cross-match galaxy catalogs, observation footprints, and all-sky images. Here we introduce HEALPix Alchemy, an open-source, pure Python implementation of a set of methods that enables rapid all-sky geometry calculations. HEALPix Alchemy is built upon HEALPix, a spatial indexing strategy that is widely used in astronomical databases as well as the native format of LIGO-Virgo-KAGRA gravitational-wave sky localization maps. Our approach leverages new multirange types built into the PostgreSQL 14 database engine. This enables fast all-sky queries against probabilistic multimessenger event localizations and telescope survey footprints. Questions such as "What are the galaxies contained within the 90% credible region of an event?" and "What is the rank-ordered list of the fields within an observing footprint with the highest probability of containing the event?" can be performed in less than a few seconds on commodity hardware using off-the-shelf cloud-managed database implementations without server-side database extensions. Common queries scale roughly linearly with the number of telescope pointings. As the number of fields grows into the hundreds or thousands, HEALPix Alchemy is orders of magnitude faster than other implementations. HEALPix Alchemy is now used as the spatial geometry engine within SkyPortal, which forms the basis of the Zwicky Transient Facility transient marshal, called Fritz.


title: Stripped-Envelope Stars in Different Metallicity Environments I. Evolutionary Phases, Classification and Populations

abstract: Massive stars that become stripped of their hydrogen envelope through binary interaction or winds can be observed either as Wolf-Rayet stars, if they have optically thick winds, or as transparent-wind stripped-envelope stars. We approximate their evolution through evolutionary models of single helium stars, and compute detailed model grids in the initial mass range 1.5 to 70 M$_{\odot}$ for metallicities between 0.01 and 0.04, from core helium ignition until core collapse. Throughout their lifetime, some stellar models expose the ashes of helium burning. We propose that models that have nitrogen-rich envelopes are candidate WN stars, while models with a carbon-rich surface are candidate WC stars during core helium burning, and WO stars afterwards. We measure metallicity dependance of the total lifetime of our models and the duration of their evolutionary phases. We propose an analytic estimate of the wind optical depth to distinguish models of Wolf-Rayet stars from transparent-wind stripped-envelope stars, and find that the luminosity ranges at which WN, WC and WO type stars can exist is a strong function of metallicity. We find that all carbon-rich models produced in our grids have optically thick winds and match the luminosity distribution of observed populations. We construct population models and predict the numbers of transparent-wind stripped-envelope stars and Wolf-Rayet stars, and derive their number ratios at different metallicities. We find that as metallicity increases, the number of transparent-wind stripped-envelope stars decreases and the number of Wolf-Rayet stars increases. At high metallicities WC and WO type stars become more common. We apply our population models to nearby galaxies, and find that populations are more sensitive to the transition luminosity between Wolf-Rayet stars and transparent-wind helium stars than to the metallicity dependent mass loss rates.


title: Red-Giant Branch Stellar Cores as Macroscopic Dark Matter Detectors

abstract: We show that macroscopic dark matter (DM) impacts on the degenerate helium cores of red-giant branch (RGB) stars can ignite helium fusion via DM-baryon elastic scattering. The onset of helium burning leads to a characteristic drop in luminosity and rise in temperature that marks the transition to a horizontal branch star. We show that such impacts can alter the RGB luminosity function of globular clusters (GCs), focusing in particular on the GC M15. Using models of M15 stars constructed with the stellar simulation code MESA, we compute the expected DM-ignition event rates and the theoretical RGB luminosity functions under the null and signal hypotheses. We constrain DM with masses $10^{17}\ {\rm g} \lesssim m_{\chi} \lesssim 10^{20}\ \rm{g}$ and geometric cross sections $10^2\ {\rm cm}^2 \lesssim \sigma_{\chi n} \lesssim 10^{7}\ \rm{ cm}^2 $ assuming that the DM in M15 is sourced by the background Milky Way halo. We also place more stringent constraints assuming that M15 formed in a DM subhalo that survives today.

categories: hep-ph astro-ph.CO astro-ph.GA astro-ph.HE astro-ph.SR


title: Probe for bound states of SU(3) fermions and colour deconfinement

abstract: Fermionic artificial matter realized with cold atoms grants access to an unprecedented degree of control on sophisticated many-body effects with an enhanced flexibility of the operating conditions. We consider three-component fermions with attractive interactions to study the formation of complex bound states whose nature goes beyond the standard fermion pairing occurring in quantum materials. Such systems display clear analogies with quark matter. Here, we address the nature of the bound states of a three-component fermionic system in a ring-shaped trap through the persistent current. In this way, we demonstrate that we can distinguish between color superfluid and trionic bound states. By analyzing finite temperature effects, we show how finite temperature can lead to the deconfinement of bound states. For weak interactions the deconfinement occurs because of scattering states. In this regime, the deconfinement depends on the trade-off between interactions and thermal fluctuations temperature. For strong interactions the features of the persistent current result from the properties of a suitable gas of bound states.

categories: cond-mat.quant-gas cond-mat.mes-hall hep-th nucl-th quant-ph


title: Impacts of Nuclear Reaction Rate Uncertainties on Sub-Chandrasekhar-Mass White Dwarf Detonations

abstract: The precise nature of Type Ia supernova (SN Ia) progenitors remains a mystery, but the detonation of a sub-Chandrasekhar-mass white dwarf (WD) has become a promising candidate. There is a growing body of work suggesting that the carbon core detonation of a sub-Chandrasekhar-mass WD can be triggered by the detonation of a helium shell accreted from a companion WD, through either inward shock convergence near the center or direct edge-lit detonation. This "double-detonation" SN Ia can be triggered by a small helium shell and is therefore well approximated by the detonation of a bare carbon-oxygen white dwarf (C/O WD). The impacts of uncertainties in experimentally and theoretically determined nuclear reaction rates on nucleosynthesis in the detonations of sub-Chandrasekhar-mass WDs have not yet been fully explored. We investigate the sensitivity of this model to nuclear reaction rate uncertainties to better constrain the nucleosynthetic yields resulting from these phenomena and identify the nuclear reaction rates whose uncertainties have the most significant impacts on nucleosynthesis. We find that the chemical abundances, and in particular those of the iron-group elements, are relatively insensitive to these nuclear reaction rate uncertainties.


title: Effective-one-body multipolar waveforms for eccentric binary black holes with non-precessing spins

abstract: We construct an inspiral-merger-ringdown eccentric gravitational-wave (GW) model for binary black holes with non-precessing spins within the effective-one-body formalism. This waveform model, SEOBNRv4EHM, extends the accurate quasi-circular SEOBNRv4HM model to eccentric binaries by including recently computed eccentric corrections up to 2PN order in the gravitational waveform modes, notably the $(l,|m|)=(2,2),(2,1),(3,3),(4,4),(5,5)$ multipoles. The waveform model reproduces the zero eccentricity limit with an accuracy comparable to the underlying quasi-circular model, with the unfaithfulness of $\lesssim1\%$ against quasi-circular numerical-relativity (NR) simulations. When compared against 28 public eccentric NR simulations from the Simulating eXtreme Spacetimes catalog with initial orbital eccentricities up to $e\simeq0.3$ and dimensionless spin magnitudes up to $+0.7$, the model provides unfaithfulness $<1\%$, showing that both the $(2,|2|)$-modes and the higher-order modes are reliably described without calibration to NR datasets in the eccentric sector. The waveform model SEOBNRv4EHM is able to qualitatively reproduce the phenomenology of dynamical captures, and can be extended to include spin-precession effects. It can be employed for upcoming observing runs with the LIGO-Virgo-KAGRA detectors and used to re-analyze existing GW catalogs to infer the eccentricity parameters for binaries with $e\lesssim0.3$ (at 20 Hz or lower) and spins up to $\lesssim 0.9-0.95$. The latter is a promising region of the parameter space where some astrophysical formation scenarios of binaries predict mild eccentricity in the ground-based detectors' bandwidth. Assessing the accuracy and robustness of the eccentric waveform model SEOBNRv4EHM for larger eccentricities and spins will require comparisons with, and, likely, calibration to eccentric NR waveforms in a larger region of the parameter space.


title: Controlled Cue Generation for Play Scripts

abstract: In this paper, we use a large-scale play scripts dataset to propose the novel task of theatrical cue generation from dialogues. Using over one million lines of dialogue and cues, we approach the problem of cue generation as a controlled text generation task, and show how cues can be used to enhance the impact of dialogue using a language model conditioned on a dialogue/cue discriminator. In addition, we explore the use of topic keywords and emotions for controlled text generation. Extensive quantitative and qualitative experiments show that language models can be successfully used to generate plausible and attribute-controlled texts in highly specialised domains such as play scripts. Supporting materials can be found at: https://catlab-team.github.io/cuegen.


title: Index theorems, generalized Hall currents and topology for gapless defect fermions

abstract: We show how the index of the fermion operator from the Euclidean action can be used to uncover the existence of gapless modes living on defects (such as edges and vortices) in topological insulators and superconductors. The 1-loop Feynman diagram that computes the index reveals an analog of the Quantum Hall current flowing on and off the defect -- even in systems without conserved currents or chiral anomalies -- and makes explicit the interplay between topology in momentum and coordinate space. We provide several explicit examples.

categories: cond-mat.mes-hall cond-mat.str-el hep-lat hep-th


title: On the space of null geodesics of a spacetime: the compact case, Engel geometry and retrievability

abstract: We compute the contact manifold of null geodesics of the family of spacetimes $\{(\mathbb{S}^2\times\mathbb{S}^1, g_\circ-\frac{d^2}{c^2}dt^2)\}_{d,c\in\mathbb{N}^+\text{ coprime}}$, with $g_\circ$ the round metric on $\mathbb{S}^2$ and $t$ the $\mathbb{S}^1$-coordinate. We find that these are the lens spaces $L(2c,1)$ together with the pushforward of the canonical contact structure on $ST\mathbb{S}^2\cong L(2,1)$ under the natural projection $L(2,1)\to L(2c,1)$. We extend this computation to $Z\times \mathbb{S}^1$ for $Z$ a Zoll manifold. On the other hand, motivated by these examples, we show how Engel geometry can be used to describe the manifold of null geodesics of a certain class of three-dimensional spacetimes, by considering the Cartan deprolongation of their Lorentz prolongation. We characterize the three-dimensional contact manifolds that are contactomorphic to the space of null geodesics of a spacetime. The characterization consists in the existence of an overlying Engel manifold with a certain foliation and, in this case, we also retrieve the spacetime.


title: Nonperturbative Negative Geometries: Amplitudes at Strong Coupling and the Amplituhedron

abstract: The amplituhedron determines scattering amplitudes in planar ${\cal N}=4$ super Yang-Mills by a single "positive geometry" in the space of kinematic and loop variables. We study a closely related definition of the amplituhedron for the simplest case of four-particle scattering, given as a sum over complementary "negative geometries", which provides a natural geometric understanding of the exponentiation of infrared (IR) divergences, as well as a new geometric definition of an IR finite observable ${\cal F}(g,z)$ - dually interpreted as the expectation value of the null polygonal Wilson loop with a single Lagrangian insertion - which is directly determined by these negative geometries. This provides a long-sought direct link between canonical forms for positive (negative) geometries, and a completely IR finite post-loop-integration observable depending on a single kinematical variable $z$, from which the cusp anomalous dimension $\Gamma_{\rm cusp}(g)$ can also be straightforwardly obtained. We study an especially simple class of negative geometries at all loop orders, associated with a "tree" structure in the negativity conditions, for which the contributions to ${\cal F}(g,z)$ and $\Gamma_{\rm cusp}$ can easily be determined by an interesting non-linear differential equation immediately following from the combinatorics of negative geometries. This lets us compute these "tree" contributions to ${\cal F}(g,z)$ and $\Gamma_{\rm cusp}$ for all values of the 't Hooft coupling. The result for $\Gamma_{\rm cusp}$ remarkably shares all main qualitative characteristics of the known exact results obtained using integrability.


title: Spinning cylinders in general relativity: a canonical form for the Lewis metrics of the Weyl class

abstract: In the main article [CQG 38 (2021) 055003], a new "canonical" form for the Lewis metrics of the Weyl class has been obtained, depending only on three parameters -- Komar mass and angular momentum per unit length, plus the angle deficit -- corresponding to a coordinate system fixed to the "distant stars" and an everywhere timelike Killing vector field. Such form evinces the local but non-global static character of the spacetime, and striking parallelisms with the electromagnetic analogue. We discuss here its generality, main physical features and important limits (the Levi-Civita static cylinder, and spinning cosmic strings). We contrast it on geometric and physical grounds with the Kerr spacetime -- as an example of a metric which is locally non-static.


title: Designing $\mathbb{Z}_2$ and $\mathbb{Z}_2 \times \mathbb{Z}_2$ topological orders in networks of Majorana bound states

abstract: Topological orders have been intrinsically identified in a class of systems such as fractional quantum Hall states and spin liquids. Accessing such states often requires extreme conditions such as low temperatures, high magnetic fields, pure samples, etc. Another approach would be to engineer the topological orders in systems with more accessible ingredients. In this work, we present networks of Majorana bound states, which are currently accessible in semiconductor nanowires proximitized to conventional superconductors, and show that the effective low-energy theory is topologically ordered. We first demonstrate the main principles in a lattice made of Kitaev superconducting chains comprising both spin species. The lattice is coupled to free magnetic moments through the Kondo interaction. We then show that at the weak coupling limit, effective ring spin interactions are induced between magnetic moments with a topological order enjoying a local $\mathbb{Z}_2 \times \mathbb{Z}_2$ gauge symmetry. We then show that the same topological order and also the $\mathbb{Z}_2$ one can be engineered in architecture patterns of semiconductor nanowires hosting Majorana bound states. The basic blocks of patterns are the time-reversal Majorana Cooper boxes coupled to each other by metallic leads, and the Majorana states are allowed to tunnel to quantum dots sitting on the vertices of the lattices. In the limit of strong onsite Coulomb interactions, where the charge fluctuations are suppressed, the magnetic moments of dots on the square and honeycomb lattices are described by topologically ordered spin models with underlying $\mathbb{Z}_2$ and $\mathbb{Z}_2 \times \mathbb{Z}_2$ gauge symmetries, respectively. Finally, we show that the latter topological order can also be realized in a network of purely Majorana zero modes in the absence of coupling to quantum dots.


title: Volume-law entanglement entropy of typical pure quantum states

abstract: The entanglement entropy of subsystems of typical eigenstates of quantum many-body Hamiltonians has been recently conjectured to be a diagnostic of quantum chaos and integrability. In quantum chaotic systems it has been found to behave as in typical pure states, while in integrable systems it has been found to behave as in typical pure Gaussian states. In this tutorial, we provide a pedagogical introduction to known results about the entanglement entropy of subsystems of typical pure states and of typical pure Gaussian states. They both exhibit a leading term that scales with the volume of the subsystem, when smaller than one half of the volume of the system, but the prefactor of the volume law is fundamentally different. It is constant (and maximal) for typical pure states, and it depends on the ratio between the volume of the subsystem and of the entire system for typical pure Gaussian states. Since particle-number conservation plays an important role in many physical Hamiltonians, we discuss its effect on typical pure states and on typical pure Gaussian states. We prove that while the behavior of the leading volume-law terms does not change qualitatively, the nature of the subleading terms can change. In particular, subleading corrections can appear that depend on the square root of the volume of the subsystem. We unveil the origin of those corrections. Finally, we discuss the connection between the entanglement entropy of typical pure states and analytical results obtained in the context of random matrix theory, as well as numerical results obtained for physical Hamiltonians.

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


title: The Contribution of AGN Accretion Disks to Hydrogen Reionization

abstract: We examine the contribution of high-redshift (z>6) active galactic nuclei (AGNs) to cosmic hydrogen reionization, by tracing the growth and ionizing output of the first generation of supermassive black holes (SMBHs). Our calculations are anchored to the observed population of z~6 quasars, and trace back the evolving spectral energy distributions (SEDs) of the accretion flows that power these early AGNs and consider a variety of growth histories, including super-Eddington accretion. Compared to a fixed-shape SED, the evolving thin disks can produce ionizing radiation that is higher by up to ~80%. Across a variety of SMBH growth scenarios, the contribution of AGNs to reionization is limited to late epochs (z<7), and remains sub-dominant compared to star-forming galaxies. This conclusion holds irrespective of the (still unknown) space density of low-luminosity z=6 AGNs, and for growth scenarios that allow super-Eddington accretion. The contribution of AGNs to reionization can extend to earlier epochs (z>8) in scenarios with relative slow SMBH mass growth, i.e., for low accretion rates and/or high spins. We finally demonstrate that our framework can reproduce the observed quasar proximity zone sizes, and that compact proximity zones around z=6 quasars can be explained by the late onset of super-Eddington accretion.


title: AdS$_5$ Holography and Higher-Derivative Supergravity

abstract: We study four-derivative corrections to 5d $\mathcal{N}=2$ minimal gauged supergravity using tools from conformal supergravity. There are two supersymmetric invariants at the four-derivative order and we show explicitly how to write down the action in the Poincar\'{e} frame in terms of their coefficients. We apply our results in the context of holography where the four-derivative terms in the supergravity action translate into corrections to the conformal anomalies of the dual 4d $\mathcal{N}=1$ SCFT. We use these results to study corrections to the entropy of the supersymmetric black string solutions of minimal supergravity. We also discuss explicit embeddings of the 5d supergravity in string and M-theory and use them to fix the higher-derivative coefficients in terms of the microscopic parameters.


title: One-phase Free Boundary Problems on RCD Metric Measure Spaces

abstract: In this paper, we consider a vector-valued one-phase Bernoulli-type free boundary problem on a metric measure space $(X,d,\mu)$ with Riemannian curvature-dimension condition $RCD(K,N)$. We first prove the existence and the local Lipschitz regularity of the solutions, provided that the space $X$ is non collapsed, i.e. $\mu$ is the $N$-dimensional Hausdorff measure of $X$. And then we show that the free boundary of the solutions is an $(N-1)$-dimensional topological manifold away from a relatively closed subset of Hausdorff dimension $\leqslant N-3$.


title: Meterstick: Benchmarking Performance Variability in Cloud and Self-hosted Minecraft-like Games Extended Technical Report

abstract: Due to increasing popularity and strict performance requirements, online games have become a workload of interest for the performance engineering community. One of the most popular types of online games is the Minecraft-like Game (MLG), in which players can terraform the environment. The most popular MLG, Minecraft, provides not only entertainment, but also educational support and social interaction, to over 130 million people world-wide. MLGs currently support their many players by replicating isolated instances that support each only up to a few hundred players under favorable conditions. In practice, as we show here, the real upper limit of supported players can be much lower. In this work, we posit that performance variability is a key cause for the lack of scalability in MLGs. We propose a novel operational model for MLGs and use it to design the first benchmark that focuses on MLG performance variability, defining specialized workloads, metrics, and processes. We conduct real-world benchmarking of MLGs and find environment-based workloads and cloud deployment to be significant sources of performance variability: peak-latency degrades sharply to 20.7 times the arithmetic mean, and exceeds by a factor of 7.4 the performance requirements. We derive actionable insights for game-developers, game-operators, and other stakeholders to tame performance variability.


title: Expanding the time domain of multiple populations: evidences of nitrogen variations in the ~1.5 Gyr old star cluster NGC 1783

abstract: We present the result of a detailed analysis of HST UV and optical deep images of the massive and young (~1.5 Gyr) stellar cluster NGC 1783 in the Large Magellanic Cloud. This system does not show evidences of multiple populations (MPs) along the red giant branch (RGB) stars. However, we find that the cluster main-sequence (MS) shows evidence of a significant broadening (50% larger than what is expected from photometric errors) along with hints of possible bimodality in the MP sensitive (F343N - F438W, F438W) color-magnitude-diagram (CMD). Such an effect is observed in all color combinations including the F343N filter, while it is not found in the optical CMDs. This observational evidence suggests we might have found light-element chemical abundance variations along the MS of NGC 1783, which represents the first detection of MPs in a system younger than 2 Gyr. A comparison with isochrones including MP-like abundances shows that the observed broadening is compatible with a N abundance enhancement of Delta([N/Fe]) ~0.3. Our analysis also confirms previous results about the lack of MPs along the cluster RGB. However, we find that the apparent disagreement between the results found on the MS and the RGB is compatible with the mixing effects linked to the first dredge-up. This study provides new key information about the MP phenomenon and suggests that star clusters form in a similar way at any cosmic age.


title: Experimental Higher-Order Interference in a Nonlinear Triple Slit

abstract: Interference between two waves is a well-known concept in physics, and its generalization to more than two waves is straight-forward. The order of interference is defined as the number of paths that interfere in a manner that cannot be reduced to patterns of a lower order. In practice, second-order interference means that in, say, a triple-slit experiment, the interference pattern when all three slits are open can be predicted from the interference patterns between all possible pairs of slits. Quantum mechanics is often said to only exhibit second-order interference. However, this is only true under specific assumptions, typically single-particles undergoing linear evolution. Here we experimentally show that nonlinear evolution can in fact lead to higher-order interference. The higher-order interference in our experiment has a simple quantum mechanical description; namely, optical coherent states interacting in a nonlinear medium. Our work shows that nonlinear evolution could open a loophole for experiments attempting to verify Born's rule by ruling out higher-order interference.


title: A stability bound on the $T$-linear resistivity of conventional metals

abstract: Perturbative considerations account for the properties of conventional metals, including the range of temperatures where the transport scattering rate is $1/\tau_\text{tr} = 2\pi \lambda T$, where $\lambda$ is a dimensionless strength of the electron-phonon coupling. The fact that measured values satisfy $\lambda \lesssim 1$ has been noted in the context of a possible "Planckian" bound on transport. However, since the electron-phonon scattering is quasi-elastic in this regime, no such Planckian considerations can be relevant. We present and analyze Monte Carlo results on the Holstein model which show that a different sort of bound is at play: a "stability" bound on $\lambda$ consistent with metallic transport. We conjecture that a qualitatively similar bound on the strength of residual interactions, which is often stronger than Planckian, may apply to metals more generally.


title: Bath deformations, islands and holographic complexity

abstract: Considering a doubly holographic model, we study the evolution of holographic subregion complexity corresponding to deformations of bath state by a relevant scalar operator, which corresponds to a renormalization group flow from the AdS-Schwarzschild to the Kasner universe in the bulk. The subregion complexity shows a discontinuous jump at Page time at a fixed perturbation, where the discontinuity depends solely on the system's parameters. We show that the amount of discontinuity decreases with the perturbation as well as with the scaling dimension of the relevant scalar operator.


title: Astrophysical Distance Scale IV. Preliminary Zero-Point Calibration of the JAGB Method in the HST/WFC3-IR Broad J-Band (F110W) Filter

abstract: We present an absolute calibration of the J-region Asymptotic Giant Branch (JAGB) method using published photometry of resolved stars in 20 nearby galaxies observed with HST using the WFC3-IR camera and the F110W (Broad J-Band) filter. True distance moduli for each of the galaxies are based on the Tip of the Red Giant Branch (TRGB) method as uniformly determined by Dalcanton et al. (2012). From a composite color-magnitude diagram composed of over 6 million stars, leading to a sample of 453 JAGB stars in these galaxies, we find M_{F110W}{JAGB} = -5.77 +/- 0.02 mag(statistical error on the mean). The external scatter seen in a comparison of the individual TRGB and the JAGB moduli is +/-0.081 mag (or 4% in distance). Some of this scatter can be attributed to small-number statistics arising from the sparse JAGB populations found in the generally low-luminosity galaxies that comprise the particular sample studied here. However, if this inter-method scatter is shared equitably between the JAGB and TRGB methods that implies that each are good to +/-0.06 mag, or better than 3% in distance.


title: A gradient sampling method with complexity guarantees for Lipschitz functions in high and low dimensions

abstract: Zhang et al. introduced a novel modification of Goldstein's classical subgradient method, with an efficiency guarantee of $O(\varepsilon^{-4})$ for minimizing Lipschitz functions. Their work, however, makes use of a nonstandard subgradient oracle model and requires the function to be directionally differentiable. In this paper, we show that both of these assumptions can be dropped by simply adding a small random perturbation in each step of their algorithm. The resulting method works on any Lipschitz function whose value and gradient can be evaluated at points of differentiability. We additionally present a new cutting plane algorithm that achieves better efficiency in low dimensions: $O(d\varepsilon^{-3})$ for Lipschitz functions and $O(d\varepsilon^{-2})$ for those that are weakly convex.


title: Structural and thermodynamics properties of pure phase alkanes, monoamides and alkane/monoamide mixtures with an ab initio based force-field model

abstract: A polarizable force-field (FF) model for short- and long-alkane chains and amide derivatives was constructed based solely on accurate quantum chemical (QC) calculations. First, the FF model accuracy was accessed by performing molecular dynamics (MD) simulations to calculate liquid-phase thermodynamic and structural properties for alkanes, for which experimental data are available. Second, The FF was then used to perform molecular dynamics simulations to calculate thermodynamic, structural and excess properties of monoamide/dodecane mixtures, namely DEHiBA/dodecane and DEHBA/dodecane. Aggregation phenomena appear for both types of mixtures and monoamide pure phases. A detailed structural analysis revealed, at small monoamide mole fraction the formation of dimers, while trimerization at larger monoamide concentrations and in their pure phases. Analysis of the relative orientation of the dimers have also been performed and showed a small difference for both phases.


title: Reading M87's DNA: A Double Helix revealing a large scale Helical Magnetic Field

abstract: We present unprecedented high fidelity radio images of the M87 jet. We analyzed Jansky Very Large Array (VLA) broadband, full polarization, radio data from 4 to 18 GHz. The observations were taken with the most extended configuration (A configuration), which allow the study of the emission of the jet up to kpc scales with a linear resolution $\sim$10 pc. The high sensitivity and resolution of our data allow to resolve the jet width. We confirm a double-helix morphology of the jet material between $\sim$300 pc and $\sim$1 kpc. We found a gradient of the polarization degree with a minimum at the projected axis and maxima at the jet edges, and a gradient in the Faraday depth with opposite signs at the jet edges. We also found that the behavior of the polarization properties along the wide range of frequencies is consistent with internal Faraday depolarization. All these characteristics strongly support the presence of a helical magnetic field in the M87 jet up to 1 kpc from the central black hole although the jet is most likely particle dominated at these large scales. Therefore, we propose a plausible scenario in which the helical configuration of the magnetic field has been maintained to large scales thanks to the presence of Kelvin-Helmholtz instabilities.


title: Discovering of the L ligand impact on luminescence enhancement of Eu(Dibenzoylmethane)3.Lx complexes employing transient absorption spectroscopy

abstract: The effect of luminescent enhancement under exchange of the auxiliary ligand in Europium(III) tris(1,3-diphenyl-1,3-propanedionato) monohydrate was investigated by steady-state and time-resolved transient absorption spectroscopy. The excited state relaxation dynamics of this complex was analyzed through a comparison of the experimental data obtained for several model compounds, namely Eu(DBM)3.NH3, Eu(DBM)3.EDA, Eu(DBM)3.Phen, Al(DBM)3 and dibenzoyl methane (DBM) in various solutions and polymer matrices. The results show there is no linear relationship between enhancement of the emission quantum yield and the luminescent lifetime, which suggests that the auxiliary ligand reduces the rate of nonradiative relaxation of the lanthanide ion, but also affects the excited state energy transfer from ligand to metal ion. Transient absorption data shows a clear correlation between the efficiency of the energy transfer and the degree of triplet state population expressed by an amplification of the signal for its excited state absorption band on going from Eu(DBM)3.H2O to the Eu(DBM).L complex. The results show that this auxiliary ligand exchange acts as a switch turning the intersystem crossing on or off as a competitive pathway for excited state relaxation of the europium(III) complexes.


title: A Generalization of Lanchester's Model of Warfare

abstract: The classical Lanchester's model is shortly reviewed and analysed, with particular attention to the critical issues that intrinsically arise from the mathematical formalization of the problem. We then generalize a particular version of such a model describing the dynamics of warfare when three or more armies are involved in the conflict. Several numerical simulations are provided.


title: Control-Oriented Modeling of Pipe Flow through Intersecting Pipe Geometries

abstract: We present control-oriented models for transient dynamics of isothermal one-dimensional gas flow through multiple pipes in series and intersecting pipe geometries. These composite models subsume algebraic constraints that would otherwise appear due to boundary conditions, so that our linear state-space models are well-suited for model-based control design for gas flow in pipe networks with non-trivial geometries.


title: Pursuing the Precision Study for Color Glass Condensate in Forward Hadron Productions

abstract: With the tremendous accomplishments of RHIC and the LHC experiments and the advent of the future Electron-Ion Collider on the horizon, the quest for compelling evidence of the color glass condensate (CGC) has become one of the most aspiring goals in the high energy Quantum Chromodynamics research. Pursuing this question requires developing the precision test of the CGC formalism. By systematically implementing the threshold resummation, we significantly improve the stability of the next-to-leading-order calculation in CGC for forward rapidity hadron productions in $pp$ and $pA$ collisions, especially in the high $p_T$ region, and obtain reliable descriptions of all existing data measured at RHIC and the LHC across all $p_T$ regions. Consequently, this technique can pave the way for the precision studies of the CGC next-to-leading-order predictions by confronting them with a large amount of precise data.


title: General memory kernels and further corrections to the variational path integral approach for the Bogoliubov-Fr\"{o}hlich Hamiltonian

abstract: The celebrated variational path integral approach to the polaron problem shows remarkable discrepancies with diagrammatic Monte Carlo for the Bogoliubov-Fr\"{o}hlich Hamiltonian which describes an impurity weakly coupled to a Bose condensed atomic gas. It has been shown both via a renormalization group approach and by the method of correlated Gaussian wavefunctions that the model has a subtle UV divergence caused by quantum fluctuations, which are not captured within Feynman's approach. In this work we address the issues with Feynman's approach and show that by extending the model action to a more general form, and by considering higher order corrections beyond the Jensen-Feynman inequality, a good agreement with diagrammatic Monte Carlo can be obtained.


title: Medial packing and elastic asymmetry stabilize the double-gyroid in block copolymers

abstract: Triply-periodic networks are among the most complex and functionally valuable self-assembled morphologies, yet they form in nearly every class of biological and synthetic soft matter building blocks. In contrast to simpler assembly motifs -- spheres, cylinders, layers -- TPN assemblies require molecules to occupy variable local domain shapes, confounding attempts to understand their formation. Here, we examine the double-gyroid (DG) network phase of block copolymer (BCP) melts, a prototypical soft self-assembly system, by using a geometric formulation of the strong stretching theory (SST) of BCP melts. The theory establishes the direct link between molecular BCP packing, thermodynamics of melt assembly and the {\it medial map}, a generic geometric measure of the center of complex shapes. We show that "medial packing" is essential for thermodynamic stability of DG in strongly-segregated melts, reconciling a long-standing contradiction between infinite- and finite-segregation theories, corroborating our SST predictions at finite-segregation via self-consistent field calculations. Additionally, we find a previously unrecognized non-monotonic dependence of DG stability on the elastic asymmetry, the comparative entropic stiffness of matrix-forming to tubular-network forming blocks. The composition window of stable DG -- intermediate to competitor lamellar and columnar phases -- widens both for large and small elastic asymmetry, seemingly overturning the heuristic view that packing frustration is localized to the tubular domains. This study demonstrates utility of geometric optimization of {\it medial tesselations} for understanding in soft-molecular assembly. As such, the particular medial-based approach deployed here is readily generalizable to study packing frustration far beyond the case of DG morphologies in neat BCP assemblies.


title: Exploring Latent Dimensions of Crowd-sourced Creativity

abstract: Recently, the discovery of interpretable directions in the latent spaces of pre-trained GANs has become a popular topic. While existing works mostly consider directions for semantic image manipulations, we focus on an abstract property: creativity. Can we manipulate an image to be more or less creative? We build our work on the largest AI-based creativity platform, Artbreeder, where users can generate images using pre-trained GAN models. We explore the latent dimensions of images generated on this platform and present a novel framework for manipulating images to make them more creative. Our code and dataset are available at http://github.com/catlab-team/latentcreative.


title: The Brain Tumor Sequence Registration (BraTS-Reg) Challenge: Establishing Correspondence Between Pre-Operative and Follow-up MRI Scans of Diffuse Glioma Patients

abstract: Registration of longitudinal brain MRI scans containing pathologies is challenging due to dramatic changes in tissue appearance. Although there has been progress in developing general-purpose medical image registration techniques, they have not yet attained the requisite precision and reliability for this task, highlighting its inherent complexity. Here we describe the Brain Tumor Sequence Registration (BraTS-Reg) challenge, as the first public benchmark environment for deformable registration algorithms focusing on estimating correspondences between pre-operative and follow-up scans of the same patient diagnosed with a diffuse brain glioma. The BraTS-Reg data comprise de-identified multi-institutional multi-parametric MRI (mpMRI) scans, curated for size and resolution according to a canonical anatomical template, and divided into training, validation, and testing sets. Clinical experts annotated ground truth (GT) landmark points of anatomical locations distinct across the temporal domain. Quantitative evaluation and ranking were based on the Median Euclidean Error (MEE), Robustness, and the determinant of the Jacobian of the displacement field. The top-ranked methodologies yielded similar performance across all evaluation metrics and shared several methodological commonalities, including pre-alignment, deep neural networks, inverse consistency analysis, and test-time instance optimization per-case basis as a post-processing step. The top-ranked method attained the MEE at or below that of the inter-rater variability for approximately 60% of the evaluated landmarks, underscoring the scope for further accuracy and robustness improvements, especially relative to human experts. The aim of BraTS-Reg is to continue to serve as an active resource for research, with the data and online evaluation tools accessible at https://bratsreg.github.io/.


title: Almost all subgeneric third-order Chow decompositions are identifiable

abstract: For real and complex homogeneous cubic polyomials in $n+1$ variables, we prove that the Chow variety of products of linear forms is generically complex identifiable for all ranks up to the generic rank minus two. By integrating fundamental results of [Oeding, Hyperdeterminants of polynomials, Adv. Math., 2012], [Casarotti and Mella, From non defectivity to identifiability, J. Eur. Math. Soc., 2021], and [Torrance and Vannieuwenhoven, All secant varieties of the Chow variety are nondefective for cubics and quaternary forms, Trans. Amer. Math. Soc., 2021] the proof is reduced to only those cases in up to $103$ variables. These remaining cases are proved using the Hessian criterion for tangential weak defectivity from [Chiantini, Ottaviani, and Vannieuwenhoven, An algorithm for generic and low-rank specific identifiability of complex tensors, SIAM J. Matrix Anal. Appl., 2014]. We also establish that the smooth loci of the real and complex Chow varieties are immersed minimal submanifolds in their usual ambient spaces.


title: Public Release and Validation of SPEC CPU2017 PinPoints

abstract: Phase-based statistical sampling methods such as SimPoints have proven to be effective at dramatically reducing the long time for architectural simulators to run large workloads such as SPEC CPU2017. However, generating and validating them is a long and tenuous process. While checkpoints of program phases, or "pinballs", of SPEC CPU2017 have been collected by other researchers and shared with the research community, they are outdated and produce errors when used with the latest versions of the Sniper architectural simulator. To facilitate our own research as well as contribute to the community, we collect and validate our own pinballs for the SPEC CPU2017 SPECspeed suite and release them to the public domain. In this work we document our methodology, the hardware and software details of the collection process, and our validation results. In terms of CPI, our pinballs have an average error rate of 12% when compared with the native whole-program benchmark execution.


title: Optical setup for a piston-cylinder type pressure cell: a double-volume approach

abstract: Measurement of the absolute value of the applied pressure in high-pressure muon and neutron experiments is a complicated task. It requires both the presence of a calibration material inside the sample volume as well as additional time for refining the response of the calibrant. Here we describe the use of optical calibrants for precise determination of the pressure value inside the piston-cylinder clamp cells. Utilizing the concept of separate volumes for the sample and the optical media, a new setup for conducting in-situ pressure measurements has been successfully tested. Pressures in both the `sample' and the `optical' volumes were proved to be the same within experimental accuracy. The use of SrB$_{4}$O$_7$:(0.01 Sm$^{2+}$, 0.03 Eu$^{2+}$) as a pressure calibrant allows for a high accuracy of pressure determination by considering up to eight fluorescence lines.


title: Sylvester-Cayley vector partitions algorithm and the Gaussian polynomials

abstract: We extend an algorithm suggested in 1858 by Sylvester and implemented in 1860 by Cayley for a problem of double partitions and apply it to derivation of explicit expressions for coefficients of the Gaussian polynomials through convolution of restricted partition functions.


title: Implementing a Category-Theoretic Framework for Typed Abstract Syntax

abstract: In previous work ("From signatures to monads in UniMath"), we described a category-theoretic construction of abstract syntax from a signature, mechanized in the UniMath library based on the Coq proof assistant. In the present work, we describe what was necessary to generalize that work to account for simply-typed languages. First, some definitions had to be generalized to account for the natural appearance of non-endofunctors in the simply-typed case. As it turns out, in many cases our mechanized results carried over to the generalized definitions without any code change. Second, an existing mechanized library on $\omega$-cocontinuous functors had to be extended by constructions and theorems necessary for constructing multi-sorted syntax. Third, the theoretical framework for the semantical signatures had to be generalized from a monoidal to a bicategorical setting, again to account for non-endofunctors arising in the typed case. This uses actions of endofunctors on functors with given source, and the corresponding notion of strong functors between actions, all formalized in UniMath using a recently developed library of bicategory theory. We explain what needed to be done to plug all of these ingredients together, modularly. The main result of our work is a general construction that, when fed with a signature for a simply-typed language, returns an implementation of that language together with suitable boilerplate code, in particular, a certified monadic substitution operation.


title: Long-Term Variation of Helioseismic Far-Side Images and What Causes It

abstract: A new time--distance far-side imaging technique was recently developed by utilizing multiple multi-skip acoustic waves. The measurement procedure is applied to 11 years of Doppler observations from the Solar Dynamics Observatory / Helioseismic and Magnetic Imager, and over 8000 far-side images of the Sun have been obtained with a 12-hour temporal cadence. The mean travel-time shifts in these images unsurprisingly vary with the solar cycle. However, the temporal variation does not show good correlations with the magnetic activity in their respective northern or southern hemisphere, but show very good anti-correlation with the global-scale magnetic activity. We investigate four possible causes of this travel-time variation. Our analysis demonstrates that the acoustic waves that are used for mapping the Sun's far side experience surface reflections around the globe, where they may interact with surface or near-surface magnetic field, and carry travel-time deficits with them. The mean far-side travel-time shifts from these acoustic waves therefore vary in phase with the Sun's magnetic activity.


title: On The Reliability Of Machine Learning Applications In Manufacturing Environments

abstract: The increasing deployment of advanced digital technologies such as Internet of Things (IoT) devices and Cyber-Physical Systems (CPS) in industrial environments is enabling the productive use of machine learning (ML) algorithms in the manufacturing domain. As ML applications transcend from research to productive use in real-world industrial environments, the question of reliability arises. Since the majority of ML models are trained and evaluated on static datasets, continuous online monitoring of their performance is required to build reliable systems. Furthermore, concept and sensor drift can lead to degrading accuracy of the algorithm over time, thus compromising safety, acceptance and economics if undetected and not properly addressed. In this work, we exemplarily highlight the severity of the issue on a publicly available industrial dataset which was recorded over the course of 36 months and explain possible sources of drift. We assess the robustness of ML algorithms commonly used in manufacturing and show, that the accuracy strongly declines with increasing drift for all tested algorithms. We further investigate how uncertainty estimation may be leveraged for online performance estimation as well as drift detection as a first step towards continually learning applications. The results indicate, that ensemble algorithms like random forests show the least decay of confidence calibration under drift.


title: What's Behind the Elephant's Trunk? Identifying Young Stellar Objects on the Outskirts of IC1396

abstract: Empirically, the estimated lifetime of a typical protoplanetary disk is $<5-10$ Myr. However, the disk lifetimes required to produce a variety of observed exoplanetary systems may exceed this timescale. Some hypothesize that this inconsistency is due to estimating disk fractions at the cores of clusters, where radiation fields external to a star-disk system can photoevaporate the disk. To test this, we have observed a field on the western outskirts of the IC 1396 star-forming region with \textit{XMM-Newton} to identify new Class III YSO cluster members. Our X-ray sample is complete for YSOs down to $1.8\,M_{\odot}$. We use a subset of these X-ray sources that have near- and mid-infrared counterparts to determine the disk fraction for this field. We find that the fraction of X-ray-detected cluster members that host disks in the field we observe is $17_{-7}^{+10}\%$ (1$\sigma$), comparable with the $29_{-3}^{+4}\%$ found in an adjacent field centered on the cometary globule IC 1396A. We re-evaluate YSO identifications in the IC 1396A field using \textit{Gaia} parallaxes compared to previous color-cut-only identifications, finding that incorporating independent distance measurements provides key additional constraints. Given the existence of at least one massive star producing an external radiation field in the cluster core, the lack of statistically significant difference in disk fraction in each observed field suggests that disk lifetimes remain consistent as a function of distance from the cluster core.


title: Event-guided Deblurring of Unknown Exposure Time Videos

abstract: Motion deblurring is a highly ill-posed problem due to the loss of motion information in the blur degradation process. Since event cameras can capture apparent motion with a high temporal resolution, several attempts have explored the potential of events for guiding deblurring. These methods generally assume that the exposure time is the same as the reciprocal of the video frame rate. However, this is not true in real situations, and the exposure time might be unknown and dynamically varies depending on the video shooting environment(e.g., illumination condition). In this paper, we address the event-guided motion deblurring assuming dynamically variable unknown exposure time of the frame-based camera. To this end, we first derive a new formulation for event-guided motion deblurring by considering the exposure and readout time in the video frame acquisition process. We then propose a novel end-to-end learning framework for event-guided motion deblurring. In particular, we design a novel Exposure Time-based Event Selection(ETES) module to selectively use event features by estimating the cross-modal correlation between the features from blurred frames and the events. Moreover, we propose a feature fusion module to fuse the selected features from events and blur frames effectively. We conduct extensive experiments on various datasets and demonstrate that our method achieves state-of-the-art performance.


title: Analyzing a Caching Model

abstract: Machine Learning has been successfully applied in systems applications such as memory prefetching and caching, where learned models have been shown to outperform heuristics. However, the lack of understanding the inner workings of these models -- interpretability -- remains a major obstacle for adoption in real-world deployments. Understanding a model's behavior can help system administrators and developers gain confidence in the model, understand risks, and debug unexpected behavior in production. Interpretability for models used in computer systems poses a particular challenge: Unlike ML models trained on images or text, the input domain (e.g., memory access patterns, program counters) is not immediately interpretable. A major challenge is therefore to explain the model in terms of concepts that are approachable to a human practitioner. By analyzing a state-of-the-art caching model, we provide evidence that the model has learned concepts beyond simple statistics that can be leveraged for explanations. Our work provides a first step towards explanability of system ML models and highlights both promises and challenges of this emerging research area.


title: Factorization and pseudofactorization of weighted graphs

abstract: For unweighted graphs, finding isometric embeddings is closely related to decompositions of $G$ into Cartesian products of smaller graphs. When $G$ is isomorphic to a Cartesian graph product, we call the factors of this product a factorization of $G$. When $G$ is isomorphic to an isometric subgraph of a Cartesian graph product, we call those factors a pseudofactorization of $G$. Prior work has shown that an unweighted graph's pseudofactorization can be used to generate a canonical isometric embedding into a product of the smallest possible pseudofactors. However, for arbitrary weighted graphs, which represent a richer variety of metric spaces, methods for finding isometric embeddings or determining their existence remain elusive, and indeed pseudofactorization and factorization have not previously been extended to this context. In this work, we address the problem of finding the factorization and pseudofactorization of a weighted graph $G$, where $G$ satisfies the property that every edge constitutes a shortest path between its endpoints. We term such graphs minimal graphs, noting that every graph can be made minimal by removing edges not affecting its path metric. We generalize pseudofactorization and factorization to minimal graphs and develop new proof techniques that extend the previously proposed algorithms due to Graham and Winkler [Graham and Winkler, '85] and Feder [Feder, '92] for pseudofactorization and factorization of unweighted graphs. We show that any $m$-edge, $n$-vertex graph with positive integer edge weights can be factored in $O(m^2)$ time, plus the time to find all pairs shortest paths (APSP) distances in a weighted graph, resulting in an overall running time of $O(m^2+n^2\log\log n)$ time. We also show that a pseudofactorization for such a graph can be computed in $O(mn)$ time, plus the time to solve APSP, resulting in an $O(mn+n^2\log\log n)$ running time.


title: Biomorphic propulsion system diving thunniform robotic fish

abstract: A biomorphic propulsion systemfor underwater robotic fish is presented. The system is based on a combination of an elastic chord with a tail fin fixed on it. The tail fin is connected to servomotor by two symmetric movable thrusts simulating muscle contraction. The propulsion system provides oscillatory tail movement with controllable amplitude and frequency. Tail oscillations results in translational movement of the robotic fish implementing the thunniform principle of locomotion. The shape of the body of the robotic fish and the tail fin were designed using computational model simulating virtual body in an aquatic medium. A prototype of robotic fish device was constructed and tested in experimental conditions. Dependencies of fish velocity on the amplitude and frequency of tail oscillations were analyzed.


title: Second generation star formation in globular clusters of different masses

abstract: By means of three-dimensional hydrodynamical simulations, we investigate the formation of second generation (SG) stars in young globular clusters of different masses. We consider clusters with a first generation of asymptotic giant branch (AGB) stars with mass 10^5 and 10^6 Msun moving at constant velocity through a uniform gas with density 10^(-24) and 10^(-23) g cm^(-3). Our setup is designed to reproduce the encounter of a young cluster with a reservoir of dense gas, e. g. during its orbital motion in the host galaxy. In the low-density models, as a result of the cooling AGB ejecta which collect in the centre, weakly perturbed by the external ram pressure, a compact central He-rich SG stellar component is formed on a timescale which decreases with increasing initial cluster mass. Our high-density models are subject to stronger ram pressure, which prevents the accumulation of the most He-rich AGB ejecta in the cluster centre. As a result, the SG is more extended and less He-enhanced than in the low-density models. By combining our results with previous simulations, we are able to study relevant, cluster-related scaling relations across a dynamical range of two orders of magnitude in mass (from 10^5 Msun to 10^7 Msun). In agreement with current observationally-based estimates, we find positive correlations between the SG-to-total number ratio and maximum He enhancement in SG stars as a function of the initial cluster mass.


title: On Control Schemes of Voltage Source Converters

abstract: This paper discusses some aspects of control schemes for voltage source converters under abnormal conditions. The control schemes are developed specifically for the situations when one or more system parameters vary significantly to the extent that the system becomes unstable with a conventional controller. The paper will present some recent works on control of grid-interactive converters for parameter variations in weak grid. The paper will also discuss some methods under abnormal dc-bus variation. Finally, the paper focuses on the control schemes suitable for ac-dc converters that addresses input frequency, voltage and load variation. All the discussed control schemes have shown robust performance under abnormal conditions.


title: Isometric Hamming embeddings of weighted graphs

abstract: A mapping $\alpha : V(G) \to V(H)$ from the vertex set of one graph $G$ to another graph $H$ is an isometric embedding if the shortest path distance between any two vertices in $G$ equals the distance between their images in $H$. Here, we consider isometric embeddings of a weighted graph $G$ into unweighted Hamming graphs, called Hamming embeddings, when $G$ satisfies the property that every edge is a shortest path between its endpoints. Using a Cartesian product decomposition of $G$ called its pseudofactorization, we show that every Hamming embedding of $G$ may be partitioned into Hamming embeddings for each irreducible pseudofactor graph of $G$, which we call its canonical partition. This implies that $G$ permits a Hamming embedding if and only if each of its irreducible pseudofactors is Hamming embeddable. This result extends prior work on unweighted graphs that showed that an unweighted graph permits a Hamming embedding if and only if each irreducible pseudofactor is a complete graph. When a graph $G$ has nontrivial pseudofactors, determining whether $G$ has a Hamming embedding can be simplified to checking embeddability of two or more smaller graphs.


title: Finiteness for self-dual classes in integral variations of Hodge structure

abstract: We generalize the finiteness theorem for the locus of Hodge classes with fixed self-intersection number, due to Cattani, Deligne, and Kaplan, from Hodge classes to self-dual classes. The proof uses the definability of period mappings in the o-minimal structure $\mathbb{R}_{\mathrm{an},\exp}$.


title: A Computational Model of Force within the Ligaments and Tendons in Progressive Collapsing Foot Deformity

abstract: Progressive collapsing foot deformity results from degeneration of the ligaments and the posterior tibial tendon (PTT). Our understanding of the relationship between the failures of them remains incomplete. We sought to improve this understanding through computational modeling of force in these soft tissues. The impact of PTT and ligament tears on force changes in the remaining ligaments was investigated using a previously validated finite element model of the foot. The ability of the PTT to restore foot alignment in a collapsed foot was evaluated by increasing the PTT force in a foot with attenuated ligaments and comparing the alignment angles to the intact foot. Rupture of any one of the ligaments led to overloading the remaining ligaments, except for the plantar naviculocuneiform, first plantar tarsometatarsal, and spring ligaments, where removing one led to unloading the other two. The attenuation of the plantar fascia, long plantar, short plantar, and spring ligaments significantly overloaded the deltoid and talocalcaneal ligaments. Isolated PTT rupture had no effect on foot alignment, but did increase the force in the deltoid and spring ligaments. Increasing the force within the PTT to 30% of body weight was effective at restoring foot alignment during quiet stance, primarily through reducing hindfoot valgus and forefoot abduction as opposed to improving arch height. The attenuation of any one ligament often leads to overload of the remaining ligaments which may lead to progressive degeneration. The PTT can maintain alignment in the collapsing foot, but at an increased load which could lead to its injury. Early intervention, in the form of ligament repair or reconstruction, might be used to prevent the progression of deformity. Moreover, strengthening the PTT through therapeutic exercise might improve its ability to restore foot alignment.


title: ELF: Exact-Lipschitz Based Universal Density Approximator Flow

abstract: Normalizing flows have grown more popular over the last few years; however, they continue to be computationally expensive, making them difficult to be accepted into the broader machine learning community. In this paper, we introduce a simple one-dimensional one-layer network that has closed form Lipschitz constants; using this, we introduce a new Exact-Lipschitz Flow (ELF) that combines the ease of sampling from residual flows with the strong performance of autoregressive flows. Further, we show that ELF is provably a universal density approximator, more computationally and parameter efficient compared to a multitude of other flows, and achieves state-of-the-art performance on multiple large-scale datasets.


title: Coupled Hubbard ladders at weak coupling: Pairing and spin excitations

abstract: The Hubbard model provides a simple framework in which one can study how certain aspects of the electronic structure of strongly interacting systems can be tuned to optimize the superconducting pairing correlations and how these changes affect the mechanisms giving rise to them. Here we use a weak-coupling random phase approximation to study a two-dimensional Hubbard model with a unidirectional modulation of the hopping amplitudes as the system evolves from the uniform square lattice to an array of weakly coupled two-leg ladders. We find that the pairing correlations retain their dominant $d_{x^2-y^2}$-wave like structure and that they are significantly enhanced for a slightly modulated lattice. This enhancement is traced backed to an increase in the strength of the spin-fluctuation pairing interacting due to favorable Fermi surface nesting in the modulated system. We then use a random-phase approximation BCS framework to examine the evolution of the neutron resonance in the superconducting state. We find that it changes only weakly for moderate modulations, but breaks up into two distinct resonances at incommensurate wave-vectors in the limit of weakly coupled ladders.


title: Designing weighted and multiplex networks for deep learning user geolocation in Twitter

abstract: Predicting the geographical location of users of social media like Twitter has found several applications in health surveillance, emergency monitoring, content personalization, and social studies in general. In this work we contribute to the research in this area by designing and evaluating new methods based on the literature of weighted multigraphs combined with state-of-the-art deep learning techniques. The explored methods depart from a similar underlying structure (that of an extended mention and/or follower network) but use different information processing strategies, e.g., information diffusion through transductive and inductive algorithms -- RGCNs and GraphSAGE, respectively -- and node embeddings with Node2vec+. These graphs are then combined with attention mechanisms to incorporate the users' text view into the models. We assess the performance of each of these methods and compare them to baseline models in the publicly available Twitter-US dataset; we also make a new dataset available based on a large Twitter capture in Latin America. Finally, our work discusses the limitations and validity of the comparisons among methods in the context of different label definitions and metrics.


title: Small deformation theory for a magnetic droplet in a rotating field

abstract: A three dimensional small deformation theory is developed to examine the motion of a magnetic droplet in a uniform rotating magnetic field. The equations describing the droplet's shape evolution are derived using two different approaches - a phenomenological equation for the tensor describing the anisotropy of the droplet, and the hydrodynamic solution using perturbation theory. We get a system of ordinary differential equations for the parameters describing the droplet's shape, which we further analyze for the particular case when the droplet's elongation is in the plane of the rotating field. The qualitative behavior of this system is governed by a single dimensionless quantity $\tau\omega$ - the product of the characteristic relaxation time of small perturbations and the angular frequency of the rotating magnetic field. Values of $\tau\omega$ determine whether the droplet's equilibrium will be closer to an oblate or a prolate shape, as well as whether it's shape will undergo oscillations as it settles to this equilibrium.We show that for small deformations, the droplet pseudo-rotates in the rotating magnetic field - its long axis follows the field, which is reminiscent of a rotation, nevertheless the torque exerted on the surrounding fluid is zero. We compare the analytic results with a boundary element simulation to determine their accuracy and the limits of the small deformation theory.


title: Rationality of Fano threefolds with terminal Gorenstein singularities, II

abstract: We classify nonrational Fano threefolds $X$ with terminal Gorenstein singularities such that $\mathrm{\rk}\, \mathrm{\Pic}(X)=1$, $(-K_X)^3\ge 8$, and $\mathrm{\rk}\, \mathrm{\Cl}(X)\le 2$.


title: Optimizing the expected maximum of two linear functions defined on a multivariate Gaussian distribution

abstract: We study stochastic optimization problems with objective function given by the expectation of the maximum of two linear functions defined on the component random variables of a multivariate Gaussian distribution. We consider random variables that are arbitrarily correlated, and we show that the problem is NP-hard even if the space of feasible solutions is unconstrained. We exploit a closed-form expression for the objective function from the literature to construct a cutting-plane algorithm that can be seen as an extension of the integer L-shaped method for a highly nonlinear function, which includes the evaluation of the c.d.f and p.d.f of a standard normal random variable with decision variables as part of the arguments. To exhibit the model's applicability, we consider two featured applications. The first is daily fantasy sports, where the algorithm identifies entries with positive returns during the 2018-2019 National Football League season. The second is a special case of makespan minimization for two parallel machines and jobs with uncertain processing times; for the special case where the jobs are uncorrelated, we prove the equivalence between its deterministic and stochastic versions and show that our algorithm can deliver a constant-factor approximation guarantee for the problem. The results of our computational evaluation involving synthetic and real-world data suggest that our discretization and upper bounding techniques lead to significant computational improvements and that the proposed algorithm outperforms sub-optimal solutions approaches.


title: Generalizations of Loday's assembly maps for Lawvere's algebraic theories

abstract: Loday's assembly maps approximate the K-theory of group rings by the K-theory of the coefficient ring and the corresponding homology of the group. We present a generalization that places both ingredients on the same footing. Building on Elmendorf--Mandell's multiplicativity results and our earlier work, we show that the K-theory of Lawvere theories is lax monoidal. This result makes it possible to present our theory in a user-friendly way without using higher categorical language. It also allows us to extend the idea to new contexts and set up a non-abelian interpolation scheme, raising novel questions. Numerous examples illustrate the scope of our extension.


title: Double cohomology of moment-angle complexes

abstract: We put a cochain complex structure ${CH}^*(\mathcal Z_K)$ on the cohomology of a moment-angle complex $\mathcal Z_K$ and call the resulting cohomology the double cohomology, ${HH}^*(\mathcal Z_K)$. We give three equivalent definitions for the differential, and compute ${HH}^*(\mathcal Z_K)$ for a family of simplicial complexes containing clique complexes of chordal graphs.


title: On the gap between deterministic and probabilistic Lyapunov exponents for continuous-time linear systems

abstract: Consider a non-autonomous continuous-time linear system in which the time-dependent matrix determining the dynamics is piecewise constant and takes finitely many values $A_1, \dotsc, A_N$. This paper studies the equality cases between the maximal Lyapunov exponent associated with the set of matrices $\{A_1, \dotsc, A_N\}$, on the one hand, and the corresponding ones for piecewise deterministic Markov processes with modes $A_1, \dotsc, A_N$, on the other hand. A fundamental step in this study consists in establishing a result of independent interest, namely, that any sequence of Markov processes associated with the matrices $A_1,\dotsc, A_N$ converges, up to extracting a subsequence, to a Markov process associated with a suitable convex combination of those matrices.


title: On permutation quadrinomials from Niho exponents in characteristic two

abstract: In a recent paper Zheng et al. characterized the coefficients of $ f(x) = x + a_1x^{s_1(2m-1)+1} + a_2x^{s_2(2m-1)+1} + a_3x^{s_3(2m-1)+1} $ over $\mathbb{ F}_2^{2m}$ that lead $ f(x) $ to be a permutation of $\mathbb{ F}_2^{2m}$ for $ (s_1, s_2, s_3) = (\frac{1}{4},1, \frac{3}{4})$. They left open the question whether those conditions were also necessary. In this paper we give a positive answer to that question.


title: Optimizing over an ensemble of neural networks

abstract: We study optimization problems where the objective function is modeled through feedforward neural networks with rectified linear unit (ReLU) activation. Recent literature has explored the use of a single neural network to model either uncertain or complex elements within an objective function. However, it is well known that ensembles of neural networks produce more stable predictions and have better generalizability than models with single neural networks, which motivates the investigation of ensembles of neural networks rather than single neural networks in decision-making pipelines. We study how to incorporate a neural network ensemble as the objective function of an optimization model and explore computational approaches for the ensuing problem. We present a mixed-integer linear program based on existing popular big-M formulations for optimizing over a single neural network. We develop a two-phase approach for our model that combines preprocessing procedures to tighten bounds for critical neurons in the neural networks with a Lagrangian relaxation-based branch-and-bound approach. Experimental evaluations of our solution methods suggest that using ensembles of neural networks yields more stable and higher quality solutions, compared to single neural networks, and that our optimization algorithm outperforms (the adaption of) a state-of-the-art approach in terms of computational time and optimality gaps.


title: Mechanical annealing and memories in a disordered solid

abstract: Shearing a disordered or amorphous solid for many cycles with a constant strain amplitude can anneal it, relaxing a sample to a steady state that encodes a memory of that amplitude. This steady state also features a remarkable stability to amplitude variations that allows one to read the memory. Here we shed new light on both annealing and memory, by considering how to mechanically anneal a sample to have as little memory content as possible. In experiments, we show that a "ring-down" protocol reaches a comparable steady state, but with no discernible memories and minimal structural anisotropy. We introduce a method to characterize the population of rearrangements within a sample, and show how it connects with the response to amplitude variation and the size of annealing steps. These techniques can be generalized to other forms of glassy matter and a wide array of disordered solids, especially those that yield by flowing homogeneously.


title: A Torelli Theorem for Graph Isomorphisms

abstract: It is known that isomorphisms of graph Jacobians induce cyclic bijections on the associated graphs. We characterize when such cyclic bijections can be strengthened to graph isomorphisms, in terms of an easily computed divisor. The result refines tools used in algebraic geometry to examine the fibers of the compactified Torelli map.


title: Slowing Down for Performance and Energy: An OS-Centric Study in Network Driven Workloads

abstract: This paper studies three fundamental aspects of an OS that impact the performance and energy efficiency of network processing: 1) batching, 2) processor energy settings, and 3) the logic and instructions of the OS networking paths. A network device's interrupt delay feature is used to induce batching and processor frequency is manipulated to control the speed of instruction execution. A baremetal library OS is used to explore OS path specialization. This study shows how careful use of batching and interrupt delay results in 2X energy and performance improvements across different workloads. Surprisingly, we find polling can be made energy efficient and can result in gains up to 11X over baseline Linux. We developed a methodology and a set of tools to collect system data in order to understand how energy is impacted at a fine-grained granularity. This paper identifies a number of other novel findings that have implications in OS design for networked applications and suggests a path forward to consider energy as a focal point of systems research.


title: Event Based Time-Vectors for auditory features extraction: a neuromorphic approach for low power audio recognition

abstract: In recent years tremendous efforts have been done to advance the state of the art for Natural Language Processing (NLP) and audio recognition. However, these efforts often translated in increased power consumption and memory requirements for bigger and more complex models. These solutions falls short of the constraints of IoT devices which need low power, low memory efficient computation, and therefore they fail to meet the growing demand of efficient edge computing. Neuromorphic systems have proved to be excellent candidates for low-power low-latency computation in a multitude of applications. For this reason we present a neuromorphic architecture, capable of unsupervised auditory feature recognition. We then validate the network on a subset of Google's Speech Commands dataset.


title: Gravelamps: Gravitational Wave Lensing Mass Profile Model Selection

abstract: We present the package Gravelamps which is designed to analyse lensed gravitational wave signals in order to constrain the mass density profile of the lensing object. Gravelamps does this via parameter estimation using the framework of bilby, which enables estimation of both the lens and the source parameters. The package can be used to study both microlensing and macrolensing cases -- where the lensing mass distribution is described by a point mass and extended mass density profile respectively -- and allows the user to easily and freely switch between the full wave optics and approximate geometric optics description. The performance of Gravelamps is demonstrated via simulated analysis of both mircolensing and macrolensing events, illustrating its capability for both parameter estimation and model selection in the wave optics and hybrid environments. To further demonstrate the utility of the package, the real gravitational-wave event GW170809 was analysed using Gravelamps; this event was found to yield no strong evidence supporting the lensing hypothesis, consistent with previously published results.


title: PantheonRL: A MARL Library for Dynamic Training Interactions

abstract: We present PantheonRL, a multiagent reinforcement learning software package for dynamic training interactions such as round-robin, adaptive, and ad-hoc training. Our package is designed around flexible agent objects that can be easily configured to support different training interactions, and handles fully general multiagent environments with mixed rewards and n agents. Built on top of StableBaselines3, our package works directly with existing powerful deep RL algorithms. Finally, PantheonRL comes with an intuitive yet functional web user interface for configuring experiments and launching multiple asynchronous jobs. Our package can be found at https://github.com/Stanford-ILIAD/PantheonRL.


title: Identifying Marginal Treatment Effects in the Presence of Sample Selection

abstract: This article presents identification results for the marginal treatment effect (MTE) when there is sample selection. We show that the MTE is partially identified for individuals who are always observed regardless of treatment, and derive uniformly sharp bounds on this parameter under three increasingly restrictive sets of assumptions. The first result imposes standard MTE assumptions with an unrestricted sample selection mechanism. The second set of conditions imposes monotonicity of the sample selection variable with respect to treatment, considerably shrinking the identified set. Finally, we incorporate a stochastic dominance assumption which tightens the lower bound for the MTE. Our analysis extends to discrete instruments. The results rely on a mixture reformulation of the problem where the mixture weights are identified, extending Lee's (2009) trimming procedure to the MTE context. We propose estimators for the bounds derived and use data made available by Deb, Munking and Trivedi (2006) to empirically illustrate the usefulness of our approach.


title: Multi-Expert Human Action Recognition with Hierarchical Super-Class Learning

abstract: In still image human action recognition, existing studies have mainly leveraged extra bounding box information along with class labels to mitigate the lack of temporal information in still images; however, preparing extra data with manual annotation is time-consuming and also prone to human errors. Moreover, the existing studies have not addressed action recognition with long-tailed distribution. In this paper, we propose a two-phase multi-expert classification method for human action recognition to cope with long-tailed distribution by means of super-class learning and without any extra information. To choose the best configuration for each super-class and characterize inter-class dependency between different action classes, we propose a novel Graph-Based Class Selection (GCS) algorithm. In the proposed approach, a coarse-grained phase selects the most relevant fine-grained experts. Then, the fine-grained experts encode the intricate details within each super-class so that the inter-class variation increases. Extensive experimental evaluations are conducted on various public human action recognition datasets, including Stanford40, Pascal VOC 2012 Action, BU101+, and IHAR datasets. The experimental results demonstrate that the proposed method yields promising improvements. To be more specific, in IHAR, Sanford40, Pascal VOC 2012 Action, and BU101+ benchmarks, the proposed approach outperforms the state-of-the-art studies by 8.92%, 0.41%, 0.66%, and 2.11 % with much less computational cost and without any auxiliary annotation information. Besides, it is proven that in addressing action recognition with long-tailed distribution, the proposed method outperforms its counterparts by a significant margin.


title: Data-driven integration of norm-penalized mean-variance portfolios

abstract: Mean-variance optimization (MVO) is known to be sensitive to estimation error in its inputs. Norm penalization of MVO programs is a regularization technique that can mitigate the adverse effects of estimation error. We augment the standard MVO program with a convex combination of parameterized $L_1$ and $L_2$-norm penalty functions. The resulting program is a parameterized quadratic program (QP) whose dual is a box-constrained QP. We make use of recent advances in neural network architecture for differentiable QPs and present a data-driven framework for optimizing parameterized norm-penalties to minimize the downstream MVO objective. We present a novel technique for computing the derivative of the optimal primal solution with respect to the parameterized $L_1$-norm penalty by implicit differentiation of the dual program. The primal solution is then recovered from the optimal dual variables. Historical simulations using US stocks and global futures data demonstrate the benefit of the data-driven optimization approach.


title: Filamentation of Fast Radio Bursts in magnetar winds

abstract: Magnetars are the most promising progenitors of Fast Radio Bursts (FRBs). Strong radio waves propagating through the magnetar wind are subject to non-linear effects, including modulation/filamentation instabilities. We derive the dispersion relation for modulations of strong waves propagating in magnetically-dominated pair plasmas focusing on dimensionless strength parameters $a_0\lesssim 1$, and discuss implications for FRBs. As an effect of the instability, the FRB radiation intensity develops sheets perpendicular to the direction of the wind magnetic field. When the FRB front expands outside the radius where the instability ends, the radiation sheets are scattered due to diffraction. The FRB scattering timescale depends on the properties of the magnetar wind. In a cold wind, the typical scattering timescale is $\tau_{\rm sc}\sim{\rm\; \mu s-ms}$ at the frequency $\nu\sim 1{\rm\; GHz}$. The scattering timescale increases at low frequencies, with the scaling $\tau_{\rm sc}\propto\nu^{-2}$. The frequency-dependent broadening of the brightest pulse of FRB 181112 is consistent with this scaling. From the scattering timescale of the pulse, one can estimate that the wind Lorentz factor is larger than a few tens. In a warm wind, the scattering timescale can approach $\tau_{\rm sc}\sim{\rm\; ns}$. Then scattering produces a frequency modulation of the observed intensity with a large bandwidth, $\Delta\nu\sim 1/\tau_{\rm sc}\gtrsim 100{\rm\; MHz}$. Broadband frequency modulations observed in FRBs could be due to scattering in a warm magnetar wind.


title: On Approximate Sequencing Policies for Linear Storage Devices

abstract: This paper investigates sequencing policies for file reading requests in linear storage devices, such as magnetic tapes. Tapes are the technology of choice for long-term storage in data centers due to their low cost and reliability. However, their physical structure imposes challenges to data retrieval operations reflected in classic optimization and operations research problems. In this work, we provide a theoretical and numerical performance analysis of low-complexity algorithms under deterministic, stochastic, and online settings, which are key in practice due to their interpretability and the large scale of existing data services. In the deterministic setting, we show that traditional policies, such as first-in first-out (FIFO), have arbitrarily poor performance, and we develop and investigate new constant-factor approximations. For the stochastic setting, we present a fully polynomial-time approximation scheme that weighs files based on their access frequencies. Finally, we investigate an online extension and propose a new algorithm with constant competitive-factor guarantees. Our numerical analysis on synthetic and real-world data suggest that the proposed algorithms may significantly outperform policies currently adopted in practice with respect to average reading times.


title: Synapse Compression for Event-Based Convolutional-Neural-Network Accelerators

abstract: Manufacturing-viable neuromorphic chips require novel computer architectures to achieve the massively parallel and efficient information processing the brain supports so effortlessly. Emerging event-based architectures are making this dream a reality. However, the large memory requirements for synaptic connectivity are a showstopper for the execution of modern convolutional neural networks (CNNs) on massively parallel, event-based (spiking) architectures. This work overcomes this roadblock by contributing a lightweight hardware scheme to compress the synaptic memory requirements by several thousand times, enabling the execution of complex CNNs on a single chip of small form factor. A silicon implementation in a 12-nm technology shows that the technique increases the system's implementation cost by only 2%, despite achieving a total memory-footprint reduction of up to 374x compared to the best previously published technique.


title: Fermion and scalar two-component dark matter from a $Z_4$ symmetry

abstract: We revisit a two-component dark matter model in which the dark matter particles are a singlet fermion ($\psi$) and a singlet scalar ($S$), both stabilized by a single $Z_4$ symmetry. The model, which was proposed by Y. Cai and A. Spray, is remarkably simple, with its phenomenology determined by just five parameters: the two dark matter masses and three dimensionless couplings. In fact, $S$ interacts with the Standard Model particles via the usual Higgs-portal, whereas $\psi$ only interacts directly with $S$, via the Yukawa terms $\overline{\psi^c}(y_s+y_p\gamma^5)\psi\,S$. We consider the two possible mass hierarchies among the dark matter particles, $M_S<M_\psi$ and $M_\psi<M_S$, and numerically investigate the consistency of the model with current bounds. The main novelties of our analysis are the inclusion of the $y_p$ coupling, the update of the direct detection limits, and a more detailed characterization of the viable parameter space. For dark matter masses below $1.3$ TeV or so, we find that the model not only is compatible with all known constraints, but that it also gives rise to observable signals in future dark matter experiments. Our results show that both dark matter particles may be observed in direct detection experiments and that the most relevant indirect detection channel is due to the annihilation of $\psi$. We also argue that this setup can be extended to other $Z_N$ symmetries and additional dark matter particles.


title: Nonclassical correlations of radiation in relation to Bell nonlocality

abstract: We analyze nonclassical correlations between outcomes of measurements conducted on two spatial radiation modes. These correlations cannot be simulated with statistical mixtures of coherent states or, more generally, with non-negative phase-space functions of quantum states and measurements. We argue that nonclassical correlations are naturally related to Bell nonlocality, the former being a more general class of quantum correlations. Indeed, it is known that local realistic as well as noncontextual models correspond to non-negative solutions to a system of linear equations for the joint probability distributions of all observables. We demonstrate that nonclassical correlations correspond to a particular solution to this system, which may have negative values even if local realism is not violated. A modification of Bell inequalities enables us to test such correlations. At the same time, our approach leads to a formulation of Bell inequalities applicable also to continuous variables. The results are illustrated with two-mode squeezed vacuum states and with hybrid entangled states (Schr\"odinger-Cat states), one mode being analyzed by balanced and the other one by unbalanced homodyne detection.


title: Learning Body-Aware 3D Shape Generative Models

abstract: The shape of many objects in the built environment is dictated by their relationships to the human body: how will a person interact with this object? Existing data-driven generative models of 3D shapes produce plausible objects but do not reason about the relationship of those objects to the human body. In this paper, we learn body-aware generative models of 3D shapes. Specifically, we train generative models of chairs, an ubiquitous shape category, which can be conditioned on a given body shape or sitting pose. The body-shape-conditioned models produce chairs which will be comfortable for a person with the given body shape; the pose-conditioned models produce chairs which accommodate the given sitting pose. To train these models, we define a "sitting pose matching" metric and a novel "sitting comfort" metric. Calculating these metrics requires an expensive optimization to sit the body into the chair, which is too slow to be used as a loss function for training a generative model. Thus, we train neural networks to efficiently approximate these metrics. We use our approach to train three body-aware generative shape models: a structured part-based generator, a point cloud generator, and an implicit surface generator. In all cases, our approach produces models which adapt their output chair shapes to input human body specifications.