title: Blind Modulo Analog-to-Digital Conversion of Vector Processes

abstract: In a growing number of applications, there is a need to digitize a (possibly high) number of correlated signals whose spectral characteristics are challenging for traditional analog-to-digital converters (ADCs). Examples, among others, include multiple-input multiple-output systems where the ADCs must acquire at once several signals at a very wide but sparsely and dynamically occupied bandwidth supporting diverse services. In such scenarios, the resolution requirements can be prohibitively high. As an alternative, the recently proposed modulo-ADC architecture can in principle require dramatically fewer bits in the conversion to obtain the target fidelity, but requires that spatiotemporal information be known and explicitly taken into account by the analog and digital processing in the converter, which is frequently impractical. Building on our recent work, we address this limitation and develop a blind version of the architecture that requires no such knowledge in the converter. In particular, it features an automatic modulo-level adjustment and a fully adaptive modulo-decoding mechanism, allowing it to asymptotically match the characteristics of the unknown input signal. Simulation results demonstrate the successful operation of the proposed algorithm.


title: GLADE+: An Extended Galaxy Catalogue for Multimessenger Searches with Advanced Gravitational-wave Detectors

abstract: We present GLADE+, an extended version of the GLADE galaxy catalogue introduced in our previous paper for multimessenger searches with advanced gravitational-wave detectors. GLADE+ combines data from six separate but not independent astronomical catalogues: the GWGC, 2MPZ, 2MASS XSC, HyperLEDA, and WISExSCOSPZ galaxy catalogues, and the SDSS-DR16Q quasar catalogue. To allow corrections of CMB-frame redshifts for peculiar motions, we calculated peculiar velocities along with their standard deviations of all galaxies having $B$-band magnitude data within redshift $z=0.05$ using the "Bayesian Origin Reconstruction from Galaxies" formalism. GLADE+ is complete up to luminosity distance $d_L=47^{+4}_{-2}$ Mpc in terms of the total expected $B$-band luminosity of galaxies, and contains all of the brightest galaxies giving 90\% of the total $B$-band and $K$-band luminosity up to $d_L\simeq 130$ Mpc. We include estimations of stellar masses and individual binary neutron star merger rates for galaxies with $W1$ magnitudes. These parameters can help in ranking galaxies in a given gravitational wave localization volume in terms of their likelihood of being hosts, thereby possibly reducing the number of pointings and total integration time needed to find the electromagnetic counterpart.

categories: astro-ph.CO astro-ph.GA astro-ph.HE gr-qc


title: Thermal runaway of silicon-based laser sails

abstract: Laser sail-based spacecraft -- where a powerful earth-based laser propels a lightweight outer-space vehicle -- have been recently proposed by the Breakthrough Starshot Initiative as a means of reaching relativistic speeds for interstellar spacetravel. The laser intensity at the sail required for this task is at least 1 GW$\cdot$m$^{-2}$ and, at such high intensities, thermal management of the sail becomes a significant challenge even when using materials with low absorption coefficients. Silicon has been proposed as one leading candidate material for the sail due to its low sub-bandgap absorption and high index of refraction, which allows for low-mass-density designs. However, here we show that the temperature-dependent bandgap of silicon combined with two-photon absorption processes can lead to thermal runaway for even the most optimistic viable assumptions of the material quality. From our calculations, we set bounds on the maximum laser intensities that can be used for a thermally stable, Si-based laser sail.


title: Parameter Tuning Strategies for Metaheuristic Methods Applied to Discrete Optimization of Structural Design

abstract: This paper presents several strategies to tune the parameters of metaheuristic methods for (discrete) design optimization of reinforced concrete (RC) structures. A novel utility metric is proposed, based on the area under the average performance curve. The process of modelling, analysis and design of realistic RC structures leads to objective functions for which the evaluation is computationally very expensive. To avoid costly simulations, two types of surrogate models are used. The first one consists of the creation of a database containing all possible solutions. The second one uses benchmark functions to create a discrete sub-space of them, simulating the main features of realistic problems. Parameter tuning of four metaheuristics is performed based on two strategies. The main difference between them is the parameter control established to perform partial assessments. The simplest strategy is suitable to tune good `generalist' methods, i.e., methods with good performance regardless the parameter configuration. The other one is more expensive, but is well suited to assess any method. Tuning results prove that Biogeography-Based Optimization, a relatively new evolutionary algorithm, outperforms other methods such as GA or PSO for such optimization problems, due to its particular approach of applying recombination and mutation operators.


title: Fast, high precision dynamics in quantum optimal control theory

abstract: Quantum optimal control theory is becoming increasingly crucial as quantum devices become more precise, but the need to quickly optimize these systems classically remains a significant bottleneck in their operation. Here we present a new theoretical quantum control framework for much faster optimization than the state of the art by replacing standard time propagation with a product of short-time propagators, each calculated using the Magnus expansion. The derived formulas for exact series terms and their gradients, based on earlier approximate integrals in a simulation setting, allow us to subsume the high cost of calculating commutators and integrals as an initial overhead. This provides an order of magnitude speedup for quantum control optimization.


title: Constraints on planetary and asteroid-mass primordial black holes from continuous gravitational-wave searches

abstract: We present new constraints on the merging rates of planetary-mass and asteroid-mass primordial black hole binaries using limits on continuous waves(quasi-monochromatic, quasi-infinite duration signals) derived from an all-sky search for isolated compact objects in the first six months of the third observing run (O3a) of LIGO/Virgo. We calculate the merging rates of these binaries in a model-independent way, and convert them to constraints on the primordial black hole abundance with minimal modelling assumptions. Our results show that we are sensitive to sources at most $\mathcal{O}(10$ pc) away for systems with chirp masses of $\mathcal{O}(10^{-5}M_\odot)$ at gravitational-wave frequencies around 30-40 Hz. These results also show that continuous-wave searches could in the future directly probe the existence of planetary-mass and asteroid-mass primordial black holes, especially those in binaries with asymmetric mass ratios. Furthermore, they demonstrate that new methods accounting for the full nonlinear gravitational-wave frequency evolution are needed to improve constraints on primordial black holes.


title: Non-Uniqueness of Non-Linear Optical Response

abstract: In recent years, non-linear optical phenomena have attracted much attention, with a particular focus on the engineering and exploitation of non-linear responses. Comparatively little study has however been devoted to the driving fields that generate these responses. In this work, we demonstrate that the relationship between a driving field and the optical response it induces is non-unique. Using a generic model for a strongly interacting system, we show that multiple candidate driving field exists, which will all generate the same response. Consequently, it is possible show that the optical response is not sufficient to determine the internal dynamics of the system, and that different solutions for the driving field will do different amounts of work on a system. This non-uniqueness phenomenon may in future be utilised to engineer internal system states without modifying its optical response.


title: Exploring the Endogenous Nature of Meme Stocks Using the Log-Periodic Power Law Model and Confidence Indicator

abstract: This study examined the endogenous nature of negative bubbles forming in meme stocks with the Log-Periodic Power Law (LPPL) Confidence Indicator (CI). A meme stock is a stock that has gained a significant amount of attention on a large social media platform such as Yahoo! or Reddit. This study examined four meme stocks including Tesla, Inc. (TSLA), GameStop Corp. (GME), Koss Corporation (KOSS), and AMC Entertainment Holdings Inc (AMC). The CI was able to detect numerous bubbles forming in meme stocks, but had difficulty in significantly predicting social media-induced exogenous rallies. This may have been due to price movements affected by external causes such as short squeezes. However, the model did provide proof for the formation of previous bubbles that could have been a catalyst for the meme stocks rallies. This study outlines the real unpredictability of many black-swan events, and further studies could be done examining exogenous bubbles.


title: Kempe Equivalent List Edge-Colorings of Planar Graphs

abstract: For a list assignment $L$ and an $L$-coloring $\varphi$, a Kempe swap in $\varphi$ is \emph{$L$-valid} if it yields another $L$-coloring. Two $L$-colorings are \emph{$L$-equivalent} if we can form one from another by a sequence of $L$-valid Kempe swaps. And a graph $G$ is \emph{$L$-swappable} if every two of its $L$-colorings are $L$-equivalent. We consider $L$-swappability of line graphs of planar graphs with large maximum degree. Let $G$ be a planar graph with $\Delta(G)\ge 9$ and let $H$ be the line graph of $G$. If $L$ is a $(\Delta(G)+1)$-assignment to $H$, then $H$ is $L$-swappable. Let $G$ be a planar graph with $\Delta(G)\ge 15$ and let $H$ be the line graph of $G$. If $L$ is a $\Delta(G)$-assignment to $H$, then $H$ is $L$-swappable. The first result is analogous to one for $L$-choosability by Borodin, which was later strengthened by Bonamy. The second result is analogous to another for $L$-choosability by Borodin, which was later strengthened by Borodin, Kostochka, and Woodall.


title: Beyond Pick-and-Place: Tackling Robotic Stacking of Diverse Shapes

abstract: We study the problem of robotic stacking with objects of complex geometry. We propose a challenging and diverse set of such objects that was carefully designed to require strategies beyond a simple "pick-and-place" solution. Our method is a reinforcement learning (RL) approach combined with vision-based interactive policy distillation and simulation-to-reality transfer. Our learned policies can efficiently handle multiple object combinations in the real world and exhibit a large variety of stacking skills. In a large experimental study, we investigate what choices matter for learning such general vision-based agents in simulation, and what affects optimal transfer to the real robot. We then leverage data collected by such policies and improve upon them with offline RL. A video and a blog post of our work are provided as supplementary material.


title: EpiBeds: Data informed modelling of the COVID-19 hospital burden in England

abstract: The first year of the COVID-19 pandemic put considerable strain on the national healthcare system in England. In order to predict the effect of the local epidemic on hospital capacity in England, we used a variety of data streams to inform the construction and parameterisation of a hospital progression model, which was coupled to a model of the generalised epidemic. We named this model EpiBeds. Data from a partially complete patient-pathway line-list was used to provide initial estimates of the mean duration that individuals spend in the different hospital compartments. We then fitted EpiBeds using complete data on hospital occupancy and hospital deaths, enabling estimation of the proportion of individuals that follow different clinical pathways, and the reproduction number of the generalised epidemic. The construction of EpiBeds makes it straightforward to adapt to different patient pathways and settings beyond England. As part of the UK response to the pandemic, EpiBeds has provided weekly forecasts to the NHS for hospital bed occupancy and admissions in England, Wales, Scotland, and Northern Ireland.


title: Time-crystalline long-range order in squeezed ground state

abstract: It is widely believed that ground-state time crystals are not realizable in realistic macroscopic systems. In particular, Watanabe and Oshikawa proved a theorem that implies the absence of the time-dependent long-range order (TDLRO) in the ground states of short-range many-body systems. However, this theorem does not forbid the presence of the ground-state TDLRO for macroscopic quantities. In this work, we investigate a simple bosonic model with a squeezed ground state and point out that the time-dependence of the ground-state TDLRO for the number operator is proportional to the square of the average number in the infinite-squeezing limit, or equivalently, the infinite average-number limit. This result implies the presence of the TDLRO for macroscopic boson number. We also discuss the physical implementations in optical, spin, and tight-binding systems, including the variants. We find an example with the macroscopic TDLRO whose essence is low-lying-state physics and another with $marginal$ TDLRO at the quantum critical point. In addition, we reconsider the definition of the ground-state time crystal in terms of the Floquet picture.

categories: cond-mat.mes-hall cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el quant-ph


title: Planning Sensing Sequences for Subsurface 3D Tumor Mapping

abstract: Surgical automation has the potential to enable increased precision and reduce the per-patient workload of overburdened human surgeons. An effective automation system must be able to sense and map subsurface anatomy, such as tumors, efficiently and accurately. In this work, we present a method that plans a sequence of sensing actions to map the 3D geometry of subsurface tumors. We leverage a sequential Bayesian Hilbert map to create a 3D probabilistic occupancy model that represents the likelihood that any given point in the anatomy is occupied by a tumor, conditioned on sensor readings. We iteratively update the map, utilizing Bayesian optimization to determine sensing poses that explore unsensed regions of anatomy and exploit the knowledge gained by previous sensing actions. We demonstrate our method's efficiency and accuracy in three anatomical scenarios including a liver tumor scenario generated from a real patient's CT scan. The results show that our proposed method significantly outperforms comparison methods in terms of efficiency while detecting subsurface tumors with high accuracy.


title: GRAPE for Fast and Scalable Graph Processing and random walk-based Embedding

abstract: Graph Representation Learning (GRL) methods opened new avenues for addressing complex, real-world problems represented by graphs. However, many graphs used in these applications comprise millions of nodes and billions of edges and are beyond the capabilities of current methods and software implementations. We present GRAPE, a software resource for graph processing and embedding that can scale with big graphs by using specialized and smart data structures, algorithms, and a fast parallel implementation of random walk-based methods. Compared with state-of-the-art software resources, GRAPE shows an improvement of orders of magnitude in empirical space and time complexity, as well as a competitive edge and node label prediction performance. GRAPE comprises about 1.7 million well-documented lines of Python and Rust code and provides 69 node embedding methods, 25 inference models, a collection of efficient graph processing utilities and over 80,000 graphs from the literature and other sources. Standardized interfaces allow seamless integration of third-party libraries, while ready-to-use and modular pipelines permit an easy-to-use evaluation of GRL methods, therefore also positioning GRAPE as a software resource to perform a fair comparison between methods and libraries for graph processing and embedding.


title: Crystal Diffusion Variational Autoencoder for Periodic Material Generation

abstract: Generating the periodic structure of stable materials is a long-standing challenge for the material design community. This task is difficult because stable materials only exist in a low-dimensional subspace of all possible periodic arrangements of atoms: 1) the coordinates must lie in the local energy minimum defined by quantum mechanics, and 2) global stability also requires the structure to follow the complex, yet specific bonding preferences between different atom types. Existing methods fail to incorporate these factors and often lack proper invariances. We propose a Crystal Diffusion Variational Autoencoder (CDVAE) that captures the physical inductive bias of material stability. By learning from the data distribution of stable materials, the decoder generates materials in a diffusion process that moves atomic coordinates towards a lower energy state and updates atom types to satisfy bonding preferences between neighbors. Our model also explicitly encodes interactions across periodic boundaries and respects permutation, translation, rotation, and periodic invariances. We significantly outperform past methods in three tasks: 1) reconstructing the input structure, 2) generating valid, diverse, and realistic materials, and 3) generating materials that optimize a specific property. We also provide several standard datasets and evaluation metrics for the broader machine learning community.


title: Last Iterate Risk Bounds of SGD with Decaying Stepsize for Overparameterized Linear Regression

abstract: Stochastic gradient descent (SGD) has been shown to generalize well in many deep learning applications. In practice, one often runs SGD with a geometrically decaying stepsize, i.e., a constant initial stepsize followed by multiple geometric stepsize decay, and uses the last iterate as the output. This kind of SGD is known to be nearly minimax optimal for classical finite-dimensional linear regression problems (Ge et al., 2019). However, a sharp analysis for the last iterate of SGD in the overparameterized setting is still open. In this paper, we provide a problem-dependent analysis on the last iterate risk bounds of SGD with decaying stepsize, for (overparameterized) linear regression problems. In particular, for last iterate SGD with (tail) geometrically decaying stepsize, we prove nearly matching upper and lower bounds on the excess risk. Moreover, we provide an excess risk lower bound for last iterate SGD with polynomially decaying stepsize and demonstrate the advantage of geometrically decaying stepsize in an instance-wise manner, which complements the minimax rate comparison made in prior works.


title: ABO: Dataset and Benchmarks for Real-World 3D Object Understanding

abstract: We introduce Amazon Berkeley Objects (ABO), a new large-scale dataset designed to help bridge the gap between real and virtual 3D worlds. ABO contains product catalog images, metadata, and artist-created 3D models with complex geometries and physically-based materials that correspond to real, household objects. We derive challenging benchmarks that exploit the unique properties of ABO and measure the current limits of the state-of-the-art on three open problems for real-world 3D object understanding: single-view 3D reconstruction, material estimation, and cross-domain multi-view object retrieval.


title: Single-component superconducting state in UTe2 at 2 K

abstract: UTe2 is a newly-discovered unconventional superconductor wherein multicomponent topological superconductivity is anticipated based on the presence of two superconducting transitions and time-reversal symmetry breaking in the superconducting state. The observation of two superconducting transitions, however, remains controversial. Here we demonstrate that UTe2 single crystals displaying an optimal superconducting transition temperature at 2 K exhibit a single transition and remarkably high quality supported by their small residual heat capacity in the superconducting state and large residual resistance ratio. Our results shed light on the intrinsic superconducting properties of UTe2 and bring into question whether UTe2 is a multicomponent superconductor at ambient pressure.


title: Stabilizing two-qubit entanglement by mimicking a squeezed environment

abstract: It is well known that qubits immersed in a squeezed vacuum environment exhibit many exotic phenomena, including dissipative entanglement stabilization. Here, we show that these effects only require interference between excitation and decay processes, and can be faithfully mimicked without non-classical light using simple classical temporal modulation. We present schemes that harnesses this idea to stabilize entanglement between two remote qubits coupled via a transmission line or waveguide, where either the qubit-waveguide coupling is modulated, or the qubits are directly driven. We analyze the resilience of these approaches against various imperfections, and also characterize the trade-off between the speed and quality of entanglement stabilization. Our protocols are compatible with state of the art cavity QED systems.


title: A note on Benli and Pollack's paper on small prime power residues

abstract: In this note, we revisit a result of Benli's and Pollack's on the number of small prime $k^{th}$ power residues. The proof is based on their idea of using reciprocity laws, but the argument is simplified and we prove a slightly stronger bound.


title: Subjective Difficulty in a Verbal Recognition-based Memory Task: Exploring Brain-behaviour Relationships at the Individual Level in Healthy Young Adults

abstract: The vast majority of fMRI studies of task-related brain activity utilize common levels of task demands and analyses that rely on the central tendencies of the data. This approach does not take into account perceived difficulty nor regional variations in brain activity between people. The results are findings of brain-behavior relationships that weaken as sample sizes increase. Participants of the current study included twenty-six healthy young adults evenly split between the sexes. The current work utilizes five parametrically modulated levels of memory load centered around each individualâs predetermined working memory cognitive capacity. Principal components analyses (PCA) identified the group-level central tendency of the data. After removing the group effect from the data, PCA identified individual-level patterns of brain activity across the five levels of task demands. Expression of the group effect significantly differed between the sexes across all load levels. Expression of the individual level patterns demonstrated a significant load by sex interaction. Furthermore, expressions of the individual maps make better predictors of response time behavior than group-derived maps.


title: Manipulating Ferrimagnets by Fields and Currents

abstract: Ferrimagnets (FIMs) can function as high-frequency antiferromagnets while being easy to detect as ferromagnets, offering unique opportunities for ultrafast device applications. While the physical behavior of FIMs near the compensation point has been widely studied, there lacks a generic understanding of FIMs where the ratio of sublattice spins can vary freely between the ferromagnetic and antiferromagnetic limits. Here we investigate the physical properties of a model two-sublattice FIM manipulated by static magnetic fields and current-induced torques. By continuously varying the ratio of sublattice spins, we clarify how the dynamical chiral modes in an FIM are intrinsically connected to their ferro- and antiferromagnetic counterparts, which reveals unique features not visible near the compensation point. In particular, we find that current-induced torques can trigger spontaneous oscillation of the terahertz exchange mode. Compared with its realization in antiferromagnets, a spin-torque oscillator using FIMs not only has a reduced threshold current density but also can be self-stabilized, obviating the need for dynamic feedback.


title: New Multielement Isotopic Compositions of Presolar SiC Grains: Implications for Their Stellar Origins

abstract: We report NanoSIMS Si and Mg-Al isotopic data (and C, N, and Ti isotopic data when available) for 85 submicron- to micron-sized presolar SiC grains from the CM2 Murchison meteorite, including 60 MS, 8 AB1, 8 X, 7 AB2, and 2 Y grains. The MS and Y grain data demonstrate that (1) C and N contamination mainly appears as surface contamination, and sufficient presputtering is needed to expose a clean grain surface for obtaining intrinsic C and N signals, and (2) Mg and Al contamination appears as adjacent grains and rims, and high-resolution imaging and the choice of small regions of interest during data reduction together are effective in suppressing the contamination. Our results strongly indicate that previous studies on presolar SiC grains could have sampled differing degrees of contamination for C, N, Mg, and Al. Compared to the literature data, our new MS and Y grains are in better agreement with carbon star observations for both the C and N isotopic ratios. By comparing our new, tighter distributions of 12C/13C, 14N/15N, and initial 26Al/27Al ratios for MS and Y grains with FRUITY AGB stellar models, we provide more stringent constraints on the occurrence of cool bottom processing and the production of 26Al in N-type carbon stars, classical asymptotic giant branch stars.


title: StARformer: Transformer with State-Action-Reward Representations for Visual Reinforcement Learning

abstract: Reinforcement Learning (RL) can be considered as a sequence modeling task: given a sequence of past state-action-reward experiences, an agent predicts a sequence of next actions. In this work, we propose State-Action-Reward Transformer (StARformer) for visual RL, which explicitly models short-term state-action-reward representations (StAR-representations), essentially introducing a Markovian-like inductive bias to improve long-term modeling. Our approach first extracts StAR-representations by self-attending image state patches, action, and reward tokens within a short temporal window. These are then combined with pure image state representations -- extracted as convolutional features, to perform self-attention over the whole sequence. Our experiments show that StARformer outperforms the state-of-the-art Transformer-based method on image-based Atari and DeepMind Control Suite benchmarks, in both offline-RL and imitation learning settings. StARformer is also more compliant with longer sequences of inputs. Our code is available at https://github.com/elicassion/StARformer.


title: Open-Set Recognition: a Good Closed-Set Classifier is All You Need?

abstract: The ability to identify whether or not a test sample belongs to one of the semantic classes in a classifier's training set is critical to practical deployment of the model. This task is termed open-set recognition (OSR) and has received significant attention in recent years. In this paper, we first demonstrate that the ability of a classifier to make the 'none-of-above' decision is highly correlated with its accuracy on the closed-set classes. We find that this relationship holds across loss objectives and architectures, and further demonstrate the trend both on the standard OSR benchmarks as well as on a large-scale ImageNet evaluation. Second, we use this correlation to boost the performance of a maximum logit score OSR 'baseline' by improving its closed-set accuracy, and with this strong baseline achieve state-of-the-art on a number of OSR benchmarks. Similarly, we boost the performance of the existing state-of-the-art method by improving its closed-set accuracy, but the resulting discrepancy with the strong baseline is marginal. Our third contribution is to present the 'Semantic Shift Benchmark' (SSB), which better respects the task of detecting semantic novelty, in contrast to other forms of distribution shift also considered in related sub-fields, such as out-of-distribution detection. On this new evaluation, we again demonstrate that there is negligible difference between the strong baseline and the existing state-of-the-art. Project Page: https://www.robots.ox.ac.uk/~vgg/research/osr/


title: Towards formalization and monitoring of microscopic traffic parameters using temporal logic

abstract: Smart cities are revolutionizing the transportation infrastructure by the integration of technology. However, ensuring that various transportation system components are operating as expected and in a safe manner is a great challenge. In this work, we propose the use of formal methods as a means to specify and reason about the traffic network's complex properties. Formal methods provide a flexible tool to define the safe operation of the traffic network by capturing non-conforming behavior, exploring various possible states of the traffic scene, and detecting any inconsistencies within it. Hence, we develop specification-based monitoring for the analysis of traffic networks using the formal language, Signal Temporal Logic. We develop monitors that identify safety-related behavior such as conforming to speed limits and maintaining appropriate headway. The framework is tested using a calibrated micro-simulated highway scenario and offline specification-based monitoring is applied to individual vehicle trajectories to understand whether they violate or satisfy the defined safety specifications. Statistical analysis of the outputs show that our approach can differentiate violating from conforming vehicle trajectories based on the defined specifications. This work can be utilized by traffic management centers to study the traffic stream properties, identify possible hazards, and provide valuable feedback for automating the traffic monitoring systems.


title: A Brief Introduction to Automatic Differentiation for Machine Learning

abstract: Machine learning and neural network models in particular have been improving the state of the art performance on many artificial intelligence related tasks. Neural network models are typically implemented using frameworks that perform gradient based optimization methods to fit a model to a dataset. These frameworks use a technique of calculating derivatives called automatic differentiation (AD) which removes the burden of performing derivative calculations from the model designer. In this report we describe AD, its motivations, and different implementation approaches. We briefly describe dataflow programming as it relates to AD. Lastly, we present example programs that are implemented with Tensorflow and PyTorch, which are two commonly used AD frameworks.


title: Micro-droplet nucleation through solvent exchange in a turbulent buoyant jet

abstract: Solvent exchange is a process involving mixing between a good solvent with dissolved solute and a poor solvent. The process creates local oversaturation which causes the nucleation of minute solute droplets. Such ternary systems on a macro-scale have remained unexplored in the turbulent regime. We experimentally study the solvent exchange process by injecting mixtures of ethanol and trans-anethole into water, forming a turbulent buoyant jet in upward direction. Locally, turbulent mixing causes oversaturation of the trans-anethole following turbulent entrainment. We optically measure the concentration of the nucleated droplets using a light attenuation technique and find that the radial concentration profile is sub-Gaussian. In contrast to the entrainment-based models, the spatial evolution of the oversaturation reveals continuous droplet nucleation downstream and radially across the jet, which we attribute to the limited mixing capacity of the jet. Though we are far from a full quantitative understanding, this work extends the knowledge on solvent exchange into the turbulent regime, and brings in a novel type of flow, broadening the scope of multicomponent, multiphase turbulent jets with phase transition.


title: Diagonalization of Polynomial-Time Deterministic Turing Machines via Nondeterministic Turing Machines

abstract: The diagonalization technique was invented by Georg Cantor to show that there are more real numbers than algebraic numbers and is very important in theoretical computer science. In this work, we enumerate all of the polynomial-time deterministic Turing machines and diagonalize against all of them by a universal nondeterministic Turing machine. As a result, we obtain that there is a language $L_d$ not accepted by any polynomial-time deterministic Turing machines but accepted by a nondeterministic Turing machine running within time $O(n^k)$ for any $k\in\mathbb{N}_1$. Based on these, we further show that $L_d\in\mathcal{NP}$ . That is, we present a proof that $\mathcal{P}$ and $\mathcal{NP}$ differ. Meanwhile, we show that there exists a language $L_s$ in $\mathcal{P}$ but the machine accepting it also runs within time $O(n^k)$ for all $k\in\mathbb{N}_1$. Lastly, we show that if $\mathcal{P}^O=\mathcal{NP}^O$ and on some rational base assumptions then the set $P^O$ of all polynomial-time deterministic oracle Turing machines with oracle $O$ is not enumerable, thus demonstrating that diagonalization technique (via a universal nondeterministic oracle Turing machine) will generally not apply to the relativized versions of the $\mathcal{P}$ versus $\mathcal{NP}$ problem.


title: Global Convergence of Triangularized Orthogonalization-free Method

abstract: This paper proves the global convergence of a triangularized orthogonalization-free method (TriOFM). TriOFM, in general, applies a triangularization idea to the gradient of an objective function and removes the rotation invariance in minimizers. More precisely, in this paper, the TriOFM works as an eigensolver for sizeable sparse matrices and obtains eigenvectors without any orthogonalization step. Due to the triangularization, the iteration is a discrete-time flow in a non-conservative vector field. The global convergence relies on the stable manifold theorem, whereas the convergence to stationary points is proved in detail in this paper. We provide two proofs inspired by the noisy power method and the noisy optimization method, respectively.


title: A Candidate for a Supergravity Anomaly

abstract: Using elementary BRS cohomology theory, this paper describes a supergravity anomaly analogous to, but very different from, the well known gauge and gravitational anomalies. It closely resembles the known gauge anomalies, but it results from a triangle diagram with two gravitinos and a gauge vector boson, rather than three gauge vector bosons, or two gravitons and a vector boson. A model that is likely to generate this supergravity anomaly is described. The coefficient of this anomaly, in perturbation theory, in a theory with unbroken supersymmetry, appears to be zero, because no relevant diagrams are linearly divergent. However, when, and only when, the theory has spontaneously broken supergravity, there are counterterms in the action which contribute to linearly divergent diagrams that can generate the anomaly. From the relevant Feynman diagrams in the theory, the general form of the anomaly can be conjectured. It is proportional to the VEV $\left< D^a \right>$ of the auxiliary field for the vector boson. So removing the anomaly generates a requirement that the effective spontaneous breaking of supergravity needs to be of the purely chiral type with $\left< F^i \right> \not = 0$ and with $\left< D^a \right>=0$.


title: Eukaryotic gene regulation at equilibrium, or non?

abstract: Models of transcriptional regulation that assume equilibrium binding of transcription factors have been very successful at predicting gene expression from sequence in bacteria. However, analogous equilibrium models do not perform as well in eukaryotes, most likely due to the largely out-of-equilibrium nature of eukaryotic regulatory processes. These processes come with unavoidable energy expenditure at the molecular level to support precise, reliable, or fast gene expression responses that could correspond to evolutionarily optimized regulatory strategies. Unfortunately, the space of possible non-equilibrium mechanisms is vast and predominantly uninteresting. The key question is therefore how this space can be navigated efficiently, to focus on mechanisms and models that are biologically relevant. In this review, we advocate for the normative role of theory - theory that prescribes rather than just describes - in providing such a focus. Theory should expand its remit beyond inferring models from data (by fitting), towards identifying non-equilibrium gene regulatory schemes (by optimizing) that may have been evolutionarily selected due to their favorable functional characteristics which outperform regulation at equilibrium. This approach can help us navigate the expanding complexity of regulatory architectures, and refocus the questions about the observed biological mechanisms from how they work towards why they have evolved in the first place. We illustrate our reasoning by toy examples for which we provide simulation code.


title: A Cross-Platform Benchmark for Interval Computation Libraries

abstract: Interval computation is widely used to certify computations that use floating point operations to avoid pitfalls related to rounding error introduced by inaccurate operations. Despite its popularity and practical benefits, support for interval arithmetic is not standardized nor available in mainstream programming languages. We propose the first benchmark for interval computations, coupled with reference solutions computed with exact arithmetic, and compare popular C and C++ libraries over different architectures, operating systems, and compilers. The benchmark allows identifying limitations in existing implementations, and provides a reliable guide on which library to use on each system. We believe that our benchmark will be useful for developers of future interval libraries, as a way to test the correctness and performance of their algorithms.


title: Mechanism for particle fractionalization and universal edge physics in quantum Hall fluids

abstract: Advancing a microscopic framework that rigorously unveils the underlying topological hallmarks of fractional quantum Hall (FQH) fluids is a prerequisite for making progress in the classification of strongly-coupled topological matter. Here we advance a second-quantization framework that helps reveal an exact fusion mechanism for particle fractionalization in FQH fluids, and uncover the fundamental structure behind the condensation of non-local operators characterizing topological order in the lowest-Landau-level (LLL). We show the first exact analytic computation of the quasielectron Berry connections leading to its fractional charge and exchange statistics, and perform Monte Carlo simulations that numerically confirm the fusion mechanism for quasiparticles. Thus, for instance, two quasiholes plus one electron of charge $e$ lead to an exact quasielectron of fractional charge $e/3$, and exchange statistics $1/3$, in a $\nu=1/3$ Laughlin fluid. We express, in a compact manner, the sequence of (both bosonic and fermionic) Laughlin second-quantized states highlighting the lack of local condensation. Furthermore, we present a rigorous constructive subspace bosonization dictionary for the bulk fluid and establish universal long-distance behavior of edge excitations by formulating a conjecture based on the DNA, or root state, of the FQH fluid.

categories: cond-mat.str-el cond-mat.mes-hall math-ph math.MP quant-ph


title: Low-Energy Signals from the Formation of Dark Matter-Nuclear Bound States

abstract: Dark matter particles may bind with nuclei if there exists an attractive force of sufficient strength. We show that a dark photon mediator of mass $\sim (10 - 100)$ MeV that kinetically mixes with Standard Model electromagnetism at the level of $\sim 10^{-3}$ generates keV-scale binding energies between dark matter and heavy elements, while forbidding the ability to bind with light elements. In underground direct detection experiments, the formation of such bound states liberates keV-scale energy in the form of electrons and photons, giving rise to mono-energetic electronic signals with a time-structure that may contain daily and seasonal modulations. We show that data from liquid-xenon detectors provides exquisite sensitivity to this scenario, constraining the galactic abundance of such dark particles to be at most $\sim 10^{-18} - 10^{-12}$ of the galactic dark matter density for masses spanning $\sim (1 - 10^5)$ GeV. However, an exponentially small fractional abundance of these dark particles is enough to explain the observed electron recoil excess at XENON1T.


title: Coalescence of black hole--neutron star binaries

abstract: We review the current status of general relativistic studies for coalescences of black hole--neutron star binaries. First, high-precision computations of black hole--neutron star binaries in quasiequilibrium circular orbits are summarized, focusing on the quasiequilibrium sequences and the mass-shedding limit. Next, the current status of numerical-relativity simulations for the merger of black hole--neutron star binaries is described. We summarize our understanding for the merger process, tidal disruption and its criterion, properties of the merger remnant and ejected material, gravitational waveforms, and gravitational-wave spectra. We also discuss expected electromagnetic counterparts to black hole--neutron star coalescences.


title: Extracting work from correlated many-body quantum systems

abstract: The presence of correlations in the input state of a non-interacting many-body quantum system can lead to an increase in the amount of work we can extract from it under global unitary processes (ergotropy). The present work explore such effect on translationally invariant systems relaying on the Matrix Product Operator formalism to define a measure of how much they are correlated. We observe that in the thermodynamic limit of large number of sites, complete work extraction can be attained for relatively small correlation strength (a reduction of a 2 factor in dB unit). Most importantly such an effect appears not to be associated with the presence of quantum correlations (e.g. entanglement) in the input state (classical correlation sources), and to be attainable by only using incoherent ergotropy. As a byproduct of our analysis we also present a rigorous formulation of the heuristic typicality argument first formulated in [Alicki and Fannes, 2013], which gives the maximum work extractable for a set of many identical quantum systems in the asymptotic limit.


title: Rossby numbers and stiffness values inferred from gravity-mode asteroseismology of rotating F- and B-type dwarfs: consequences for mixing, transport, magnetism, and convective penetration

abstract: Multi-dimensional (magneto-)hydrodynamical simulations of physical processes in stellar interiors depend on a multitude of uncalibrated free parameters, which set the spatial and time scales of their computations. We aim to provide an asteroseismic calibration of the wave and convective Rossby numbers, and of the stiffness at the interface between the convective core and radiative envelope of intermediate-mass stars. We deduce these quantities for rotating dwarfs from the observed properties of their identified gravity and gravito-inertial modes. We rely on near-core rotation rates and asteroseismic models of 26 B- and 37 F-type dwarf pulsators derived from 4-year Kepler space photometry, high-resolution spectroscopy and Gaia astrometry in the literature to deduce their convective and wave Rossby numbers. We compute the stiffness at the convection/radiation interface from the inferred maximum buoyancy frequency at the interface and the convective turnover frequency in the core. We use those asteroseismically inferred quantities to make predictions of convective penetration levels, local flux levels of gravito-inertial waves triggered by the convective core, and of the cores' potential rotational and magnetic states. Our sample of 63 gravito-inertial mode pulsators covers near-core rotation rates from almost zero up to the critical rate. The frequencies of their identified modes lead to models with stiffness values between $10^{2.69}$ and $10^{3.60}$ for the B-type pulsators, while those of F-type stars cover the range from $10^{3.47}$ to $10^{4.52}$. The convective Rossby numbers derived from the maximum convective diffusion coefficient in the convective core, based on mixing length theory and a value of the mixing length coefficient relevant for these pulsators, vary between $10^{-2.3}$ and $10^{-0.8}$ for B-type stars and $10^{-3}$ and $10^{-1.5}$ for F-type stars. (abridged)


title: Factorization and complex couplings in SYK and in Matrix Models

abstract: We consider the factorization problem in toy models of holography, in SYK and in Matrix Models. In a theory with fixed couplings, we introduce a fictitious ensemble averaging by inserting a projector onto fixed couplings. We compute the squared partition function and find that at large $N$ for a typical choice of the fixed couplings it can be approximated by two terms: a "wormhole" plus a "pair of linked half-wormholes". This resolves the factorization problem. We find that the second, half-wormhole, term can be thought of as averaging over the imaginary part of the couplings. In SYK, this reproduces known results from a different perspective. In a matrix model with an arbitrary potential, we propose the form of the "pair of linked half-wormholes" contribution. In GUE, we check that errors are indeed small for a typical choice of the hamiltonian. Our computation relies on a result by Brezin and Zee for a correlator of resolvents in a "deterministic plus random" ensemble of matrices.


title: Low frequency radio properties of the $z>5$ quasar population

abstract: Optically luminous quasars at $z > 5$ are important probes of super-massive black hole (SMBH) formation. With new and future radio facilities, the discovery of the brightest low-frequency radio sources in this epoch would be an important new probe of cosmic reionization through 21-cm absorption experiments. In this work, we systematically study the low-frequency radio properties of a sample of 115 known spectroscopically confirmed $z>5$ quasars using the second data release of the Low Frequency Array (LOFAR) Two Metre Sky survey (LoTSS-DR2), reaching noise levels of $\sim$80 $\mu$Jy beam$^{-1}$ (at 144 MHz) over an area of $\sim5720$ deg$^2$. We find that 41 sources (36%) are detected in LoTSS-DR2 at $>2 \sigma$ significance and we explore the evolution of their radio properties (power, spectral index, and radio loudness) as a function of redshift and rest-frame ultra-violet properties. We obtain a median spectral index of $-0.29^{+0.10}_{-0.09}$ by stacking 93 quasars using LoTSS-DR2 and Faint Images of the Radio Sky at Twenty Centimetres (FIRST) data at 1.4 GHz, in line with observations of quasars at $z<3$. We compare the radio loudness of the high-$z$ quasar sample to a lower-$z$ quasar sample at $z\sim2$ and find that the two radio loudness distributions are consistent with no evolution, although the low number of high-z quasars means that we cannot rule out weak evolution. Furthermore, we make a first order empirical estimate of the $z=6$ quasar radio luminosity function, which is used to derive the expected number of high-$z$ sources that will be detected in the completed LoTSS survey. This work highlights the fact that new deep radio observations can be a valuable tool in selecting high-$z$ quasar candidates for follow-up spectroscopic observations by decreasing contamination of stellar dwarfs and reducing possible selection biases introduced by strict colour cuts.


title: Investigating the Effect of Natural Language Explanations on Out-of-Distribution Generalization in Few-shot NLI

abstract: Although neural models have shown strong performance in datasets such as SNLI, they lack the ability to generalize out-of-distribution (OOD). In this work, we formulate a few-shot learning setup and examine the effects of natural language explanations on OOD generalization. We leverage the templates in the HANS dataset and construct templated natural language explanations for each template. Although generated explanations show competitive BLEU scores against groundtruth explanations, they fail to improve prediction performance. We further show that generated explanations often hallucinate information and miss key elements that indicate the label.


title: A detailed study of the bridge of excess X-ray emission between the galaxy clusters Abell 2029 and Abell 2033

abstract: We examine Suzaku, XMM-Newton, and Chandra observations of the Abell 2029/2033 system to investigate the nature of a bridge of X-ray emission joining the two galaxy clusters. By modelling the contributions from the outskirts of the two clusters, and excluding the emission from the southern infalling group and the background group LOS9, we find a significant excess of X-ray emission between the two clusters at the level of 6.5-7.0$\sigma$, depending on the choice of model, that cannot be explained by the overlap of the clusters. This excess component to the surface brightness is consistent with being emission from a filament with roughly 1.0 Mpc wide. The derived emission measure for the gas associated with the filament yields an average gas density of $3.7^{+1.0}_{-0.7} \times 10^{-5}$ cm$^{-3}$, corresponding roughly to 160 times the mean baryon density of the Universe. The Suzaku X-ray spectrum of the excess emission indicates that it is significantly colder ($1.4_{-0.5}^{+0.7}$ keV) than the cluster outskirts emission from the two clusters ($\sim$ 5 keV), statistically consistent with the temperature expected from the hottest and densest parts of the warm-hot intergalactic medium (WHIM). The geometry, density, and temperature are similar to those found from X-ray studies of the Abell 222/223 filament.


title: Inside an Asymptotically Flat Hairy Black Hole

abstract: We study the interior of a recently constructed family of asymptotically flat, charged black holes that develop (charged) scalar hair as one increases their charge at fixed mass. Inside the horizon, these black holes resemble the interior of a holographic superconductor. There are analogs of the Josephson oscillations of the scalar field, and the final Kasner singularity depends very sensitively on the black hole parameters near the onset of the instability. In an Appendix, we give a general argument that Cauchy horizons cannot exist in a large class of stationary black holes with scalar hair.


title: The Spitzer/IRAC Legacy over the GOODS Fields: Full-Depth 3.6, 4.5, 5.8 and 8.0um Mosaics and Photometry for > 9000 Galaxies at z~3.5-10 from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS)

abstract: We present the deepest Spitzer/IRAC $3.6$, $4.5$, $5.8$ and $8.0\mu$m wide-area mosaics yet over the GOODS-N and GOODS-S fields as part of the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) project. We reduced and mosaicked in a self-consistent way observations taken by the 11 different Spitzer/IRAC programs over the two GOODS fields from 12 years of Spitzer cryogenic and warm mission data. The cumulative depth in the $3.6\mu$m and $4.5\mu$m bands amounts to $\sim 4260$ hr, $\sim 1220$ hr of which are new very deep observations from the GREATS program itself. In the deepest area, the full-depth mosaics reach $\gtrsim200$ hr over an area of $\sim100$ arcmin$^2$, corresponding to a sensitivity of $\sim29$ AB magnitude at $3.6\mu$m ($1\sigma$ for point sources). Archival cryogenic $5.8\mu$m and $8.0\mu$m band data (a cumulative 976 hr) are also included in the release. The mosaics are projected onto the tangential plane of CANDELS/GOODS at a $0.3''$ pixel$^{-1}$ scale. This paper describes the methodology enabling, and the characteristics of, the public release of the mosaic science images, the corresponding coverage maps in the four IRAC bands, and the empirical Point-Spread Functions (PSFs). These PSFs enable mitigation of the source blending effects by taking into account the complex position-dependent variation in the IRAC images. The GREATS data products are in the Infrared Science Archive (IRSA). We also release the deblended $3.6$-to-$8.0\mu$m photometry for $9192$ Lyman-Break galaxies at $z\sim3.5-10$. GREATS will be the deepest mid-infrared imaging until JWST and, as such, constitutes a major resource for characterizing early galaxy assembly.


title: Warping Away Gravitational Instabilities in Protoplanetary Discs

abstract: We perform 3D SPH simulations of warped, non-coplanar gravitationally unstable discs to show that as the warp propagates through the self-gravitating disc, it heats up the disc rendering it gravitationally stable. Thus losing their spiral structure and appearing completely axisymmetric. In their youth, protoplanetary discs are expected to be massive and self-gravitating, which results in non-axisymmetric spiral structures. However recent observations of young protoplanetary discs with ALMA have revealed that discs with large-scale spiral structure are rarely observed in the midplane. Instead, axisymmetric discs with some also having ring & gap structures are more commonly observed. Our work invloving warps, non-coplanar disc structures that are expected to commonly occur in young discs, potentially resolves this discrepancy between observations and theoretical predictions. We demonstrate that they are able to suppress the large-scale spiral structure of self-gravitating protoplanetary discs.


title: The Neuroscience of Advanced Scientific Concepts

abstract: Cognitive neuroscience methods can identify the fMRI-measured neural representation of familiar individual concepts, such as apple, and decompose them into meaningful neural and semantic components. This approach was applied here to determine the neural representations and underlying dimensions of representation of far more abstract physics concepts related to matter and energy, such as fermion and dark matter, in the brains of 10 Carnegie Mellon physics faculty members who thought about the main properties of each of the concepts. One novel dimension coded the measurability vs. immeasurability of a concept. Another novel dimension of representation evoked particularly by post-classical concepts was associated with four types of cognitive processes, each linked to particular brain regions: (1) Reasoning about intangibles, taking into account their separation from direct experience and observability; (2) Assessing consilience with other, firmer knowledge; (3) Causal reasoning about relations that are not apparent or observable; and (4) Knowledge management of a large knowledge organization consisting of a multi-level structure of other concepts. Two other underlying dimensions, previously found in physics students, periodicity, and mathematical formulation, were also present in this faculty sample. The data were analyzed using factor analysis of stably responding voxels, a Gaussian-na\"ive Bayes machine-learning classification of the activation patterns associated with each concept, and a regression model that predicted activation patterns associated with each concept based on independent ratings of the dimensions of the concepts. The findings indicate that the human brain systematically organizes novel scientific concepts in terms of new dimensions of neural representation.


title: Thermal melting of a quantum electron solid in the presence of strong disorder: Anderson versus Wigner

abstract: We consider temperature-induced melting of a Wigner solid in one dimensional (1D) and two dimensional (2D) lattices of electrons interacting via the long-range Coulomb interaction in the presence of strong disorder arising from charged impurities in the system. The system simulates semiconductor-based 2D electron layers where Wigner crystallization is often claimed to be observed experimentally. Using exact diagonalization and utilizing the inverse participation ratio as well as conductance to distinguish between the localized insulating solid phase and the extended metallic liquid phase, we find that the effective melting temperature may be strongly enhanced by disorder since the disordered crystal typically could be in a localized glassy state incorporating the combined nonperturbative physics of both Anderson localization and Wigner crystallization. This disorder-induced enhancement of the melting temperature may explain why experiments often manage to observe insulating disorder-pinned Wigner solids in spite of the experimental temperature being decisively far above the theoretical melting temperature of the pristine Wigner crystal phase in many cases.


title: Super stellar abundances of alkali metals suggest significant migration for Hot Jupiters

abstract: We investigate the origin of the measured over-abundance of alkali metals in the atmospheres of hot gas giants, relative to both their host stars and their atmospheric water abundances. We show that formation exterior to the water snow line followed by inward disc-driven migration results in excess accretion of oxygen-poor, refractory-rich material from within the snow-line. This naturally leads to enrichment of alkali metals in the planetary atmosphere relative to the bulk composition of its host star but relative abundances of water that are similar to the stellar host. These relative abundances cannot be explained by in situ formation which places the refractory elements in the planetary deep interior rather than the atmosphere. We therefore suggest that the measured compositions of the atmospheres of hot Jupiters are consistent with significant migration for at least a subset of hot gas giants. Our model makes robust predictions about atmospheric composition that can be confirmed with future data from JWST and Ariel.


title: The impact of cosmic rays on dynamical balance and disk-halo interaction in Lstar disk galaxies

abstract: Cosmic rays (CRs) are an important component in the interstellar medium (ISM), but their effect on the dynamics of the disk-halo interface (< 10 kpc from the disk) is still unclear. We study the influence of CRs on the gas above the disk with high-resolution FIRE-2 cosmological simulations of late-type Lstar galaxies at redshift around zero. We compare runs with and without CR feedback (with constant anisotropic diffusion around 3e29 cm^2/s and streaming). Our simulations capture the relevant disk halo interactions, including outflows, inflows, and galactic fountains. Extra-planar gas in all of the runs satisfies dynamical balance, where total pressure balances the weight of the overlying gas. While the kinetic pressure from non-uniform motion (>1-kpc scale) dominates in the midplane, thermal and bulk pressures (or CR pressure if included) take over at large heights. We find that with CR feedback, (1) the warm (1e4K) gas is slowly accelerated by CRs; (2) the hot (> 5e5K) gas scale height is suppressed; (3) the warm-hot (2e4-5e5K) medium becomes the most volume-filling phase in the disk-halo interface. We develop a novel conceptual model of the near-disk gas dynamics in low-redshift Lstar galaxies: With CRs, the disk-halo interface is filled with CR-driven warm winds and hot super-bubbles that are propagating into the CGM with a small fraction falling back to the disk. Without CRs, most outflows from hot superbubbles are trapped by the existing hot halo and gravity, so typically they form galactic fountains.


title: Galaxy cluster strong lensing cosmography: cosmological constraints from a sample of regular galaxy clusters

abstract: Cluster strong lensing cosmography is a promising probe of the background geometry of the Universe and several studies have emerged, thanks to the increased quality of observations using space and ground-based telescopes. For the first time, we use a sample of five cluster strong lenses to measure the values of cosmological parameters and combine them with those from classical probes. In order to assess the degeneracies and the effectiveness of strong-lensing cosmography in constraining the background geometry of the Universe, we adopt four cosmological scenarios. We find good constraining power on the total matter density of the Universe ($\Omega_{\rm m}$) and the equation of state of the dark energy parameter $w$. For a flat $w$CDM cosmology, we find $\Omega_{\rm m} = 0.30_{-0.11}^{+0.09}$ and $w=-1.12_{-0.32}^{+0.17}$ from strong lensing only. Interestingly, we show that the constraints from the Cosmic Microwave Background (CMB) are improved by factors of 2.5 and 4.0 on $\Omega_{\rm m}$ and $w$, respectively, when combined with our posterior distributions in this cosmological model. In a scenario where the equation of state of dark energy evolves with redshift, the strong lensing constraints are compatible with a cosmological constant (i.e. $w=-1$). In a curved cosmology, our strong lensing analyses can accommodate a large range of values for the curvature of the Universe of $\Omega_{\rm k}=0.28_{-0.21}^{+0.16}$. In all cosmological scenarios, we show that our strong lensing constraints are complementary and in good agreement with measurements from the CMB, baryon acoustic oscillations and Type Ia supernovae. Our results show that cluster strong lensing cosmography is a potentially powerful probe to be included in the cosmological analyses of future surveys.


title: LYRA II: Cosmological dwarf galaxy formation with inhomogeneous Population III enrichment

abstract: We present the simulation of a $2\times10^9 M_\odot$ halo mass cosmological dwarf galaxy run to $z=0$ at 4 solar mass gas resolution with resolved supernova feedback. We compare three simple subgrid implementations for the inhomogeneous chemical enrichment from Population III stars and compare them to constraints from Local Group dwarf galaxies. The employed model, LYRA, is a novel high resolution galaxy formation model built for the moving mesh code AREPO, which is marked by a resolved multi-phase interstellar medium, single stars and individual supernova events. The resulting reionization relic is characterized by a short ($<1.5$ Gyr) star formation history that is repeatedly brought to a standstill by violent bursts of feedback. Star formation is reignited for a short duration due to a merger at $z\approx4$ and then again at $z\approx0.2-0$ after sustained gas accretion. Our model $z=0$ galaxy matches the stellar mass, size, stellar kinematics and metallicity relations of Local Group dwarf galaxies well. The dark matter profile does not exhibit a core in any version of the model. We show that the host halo masses of Population III stars affect the assembly history of dwarf galaxies. This manifests itself through the initial gaseous collapse in the progenitor halos, affecting the central density of the stellar component and through the accretion of luminous substructure.


title: Candidate isolated neutron stars in the 4XMM-DR10 catalog of X-ray sources

abstract: Most isolated neutron stars have been discovered thanks to the detection of their pulsed non-thermal emission, at wavelengths spanning from radio to gamma-rays. However, if the beamed non-thermal radiation does not intercept our line of sight or it is too faint or absent, isolated neutron stars can also be detected through their thermal emission, which peaks in the soft X-ray band and is emitted nearly isotropically. In the past thirty years, several thermally-emitting isolated neutron stars have been discovered thanks to X-ray all-sky surveys, observations targeted at the center of supernova remnants, or as serendipitous X-ray sources. Distinctive properties of these relatively rare X-ray sources are very soft spectra and high ratios of X-ray to optical flux. The recently released 4XMM-DR10 catalog contains more than half a million X-ray sources detected with the XMM-Newton telescope in the 0.2-10 keV range in observations carried out from 2000 to 2019. Based on a study of the spectral properties of these sources and on cross-correlations with catalogs of possible counterparts, we have carried out a search of isolated neutron stars, finding four potential candidates. The spectral and long-term variability analysis of these candidates, using also Chandra and Swift-XRT data, allowed us to point out the most interesting sources deserving further multiwavelength investigations.


title: Length and area generating functions for height-restricted Motzkin meanders

abstract: We derive the length and area generating function of planar height-restricted forward-moving discrete paths of increments +1, 0, or -1 with arbitrary starting and ending points, the so-called Motzkin meanders, and the more general length-area generating functions for Motzkin paths with markers monitoring the number of passages from the two height boundaries ('floor' and 'ceiling') and the time spent there. The results are obtained by embedding Motzkin paths in a two-step anisotropic Dyck path process and using propagator, exclusion statistics and bosonization techniques. We also present a cluster expansion of the logarithm of the generating functions that makes their polynomial structure explicit. These results are relevant to the derivation of statistical mechanical properties of physical systems such as polymers, vesicles, and solid-on-solid interfaces.

categories: math-ph cond-mat.stat-mech hep-th math.CO math.MP


title: The discovery of rest-frame UV colour gradients and a diversity of dust morphologies in bright z ~ 7 Lyman-break galaxies

abstract: We present deep ALMA dust continuum observations for a sample of luminous ($M_{\rm UV} < -22$) star-forming galaxies at $z \simeq 7$. We detect five of the six sources in the far-infrared (FIR), providing key constraints on the obscured star-formation rate (SFR) and the infrared-excess-$\beta$ (IRX-$\beta$) relation without the need for stacking. Despite the galaxies showing blue rest-frame UV slopes ($\beta \simeq -2$) we find that 35-75 percent of the total SFR is obscured. We find the IRX-$\beta$ relation derived for these $z \simeq 7$ sources is consistent with that found for local star-burst galaxies. Using our relatively high-resolution (FWHM $\simeq 0.7\,{\rm arcsec}$) observations we identify a diversity of dust morphologies in the sample. We find both compact emission that appears offset relative to the unobscured components and extended dust emission that is co-spatial with the rest-frame UV light. In the majority of the sources we detect strong rest-frame UV colour gradients (with up to $\Delta \beta \simeq 0.7$-$1.4$) as probed by the multi-band UltraVISTA ground-based data. The observed redder colours are spatially correlated with the location of the FIR detection. Our results show that even in bright Lyman-break galaxies at $z \simeq 7$ the peak of the star-formation is typically hosted by the fainter, redder, regions in the rest-frame UV, which have an obscured fraction of $f_{\rm obs} \ge 0.8$. As well as demonstrating the importance of dust obscured star-formation within the Epoch of Reionization, these observations provide an exciting taster of the rich spatially resolved datasets that will be obtained from JWST and high-resolution ALMA follow-up at these redshifts.


title: Integrating three-loop modular graph functions and transcendentality of string amplitudes

abstract: Modular graph functions (MGFs) are $\mathrm{SL}(2,\mathbb{Z})$-invariant functions on the Poincar\'e upper half-plane associated with Feynman graphs of a conformal scalar field on a torus. The low-energy expansion of genus-one superstring amplitudes involves suitably regularized integrals of MGFs over the fundamental domain for $\mathrm{SL}(2,\mathbb{Z})$. In earlier work, these integrals were evaluated for all MGFs up to two loops and for higher loops up to weight six. These results led to the conjectured uniform transcendentality of the genus-one four-graviton amplitude in Type II superstring theory. In this paper, we explicitly evaluate the integrals of several infinite families of three-loop MGFs and investigate their transcendental structure. Up to weight seven, the structure of the integral of each individual MGF is consistent with the uniform transcendentality of string amplitudes. Starting at weight eight, the transcendental weights obtained for the integrals of individual MGFs are no longer consistent with the uniform transcendentality of string amplitudes. However, in all the cases we examine, the violations of uniform transcendentality take on a special form given by the integrals of triple products of non-holomorphic Eisenstein series. If Type II superstring amplitudes do exhibit uniform transcendentality, then the special combinations of MGFs which enter the amplitudes must be such that these integrals of triple products of Eisenstein series precisely cancel one another. Whether this indeed is the case poses a novel challenge to the conjectured uniform transcendentality of genus-one string amplitudes.


title: Low Latency Detection of Massive Black Hole Binaries

abstract: The next decade is expected to see the launch of one or more space based gravitational wave detectors: the European lead Laser Interferometer Space Antenna (LISA); and one or more Chinese mission concepts, Taiji and TianQin. One of the primary scientific targets for these missions are the mergers of black holes with masses between $10^3 M_\odot$ and $10^8 M_\odot$. These systems may produce detectable electromagnetic signatures in additional to gravitational waves due to the presence of gas in mini-disks around each black hole, and a circumbinary disk surrounding the system. The electromagnetic emission may occur before, during and after the merger. In order to have the best chance of capturing all phases of the emission it is imperative that the gravitational wave signals can be detected in low latency, and used to produce reliable estimates for the sky location and distance to help guide the search for counterparts. The low latency detection also provides a starting point for the ``global fit'' of the myriad signals that are simultaneously present in the data. Here a low latency analysis pipeline is presented that is capable of analyzing months of data in just a few hours using a laptop from the last decade. The problem of performing a global fit is avoided by whitening out the bright foreground produced by nearby galactic binaries. The performance of the pipeline is illustrated using simulated data from the LISA Data Challenge.


title: The non-chiral intermediate Heisenberg ferromagnet equation

abstract: We present and solve a soliton equation which we call the non-chiral intermediate Heisenberg ferromagnet (ncIHF) equation. This equation, which depends on a parameter $\delta >0$, describes the time evolution of two coupled spin densities propagating on the real line, and in the limit $\delta \to \infty$ it reduces to two decoupled half-wave maps (HWM) equations of opposite chirality. We show that the ncIHF equation is related to the A-type hyperbolic spin Calogero-Moser (CM) system in two distinct ways: (i) it is obtained as a particular continuum limit of a Inozemtsev-type spin chain related to this CM system, (ii) it has multi-soliton solutions obtained by a spin-pole ansatz and with parameters satisfying the equations of motion of a complexified version of this CM system. The integrability of the ncIHF equation is shown by constructing a Lax pair. We also propose a periodic variant of the ncIHF equation related to the A-type elliptic spin CM system.


title: Residual Abundances in GALAH DR3: Implications for Nucleosynthesis and Identification of Unique Stellar Populations

abstract: We investigate the [X/Mg] abundances of 16 elements for 82,910 Galactic disk stars from GALAH+ DR3. We fit the median trends of low-Ia and high-Ia populations with a two-process model, which describes stellar abundances in terms of a prompt core-collapse and delayed Type-Ia supernova component. For each sample star, we fit the amplitudes of these two components and compute the residual $\Delta$[X/H] abundances from this two-parameter fit. We find RMS residuals $\lesssim 0.07$ dex for well-measured elements and correlated residuals among some elements (such as Ba, Y, and Zn) that indicate common enrichment sources. From a detailed investigation of stars with large residuals, we infer that roughly $40\%$ of the large deviations are physical and $60\%$ are caused by problematic data such as unflagged binarity, poor wavelength solutions, and poor telluric subtraction. As one example of a population with distinctive abundance patterns, we identify 15 stars that have 0.3-0.6 dex enhancements of Na but normal abundances of other elements from O to Ni and positive average residuals of Cu, Zn, Y, and Ba. We measure the median elemental residuals of 14 open clusters, finding systematic $\sim0.1-0.4$ dex enhancements of O, Ca, K, Y, and Ba and $\sim0.2$ dex depletion of Cu in young clusters. Finally, we present a restricted three-process model where we add an asymptotic giant branch star (AGB) component to better fit Ba and Y. With the addition of the third process, we identify a population of stars, preferentially young, that have much higher AGB enrichment than expected from their SNIa enrichment.


title: Molecular Graph Generation via Geometric Scattering

abstract: Graph neural networks (GNNs) have been used extensively for addressing problems in drug design and discovery. Both ligand and target molecules are represented as graphs with node and edge features encoding information about atomic elements and bonds respectively. Although existing deep learning models perform remarkably well at predicting physicochemical properties and binding affinities, the generation of new molecules with optimized properties remains challenging. Inherently, most GNNs perform poorly in whole-graph representation due to the limitations of the message-passing paradigm. Furthermore, step-by-step graph generation frameworks that use reinforcement learning or other sequential processing can be slow and result in a high proportion of invalid molecules with substantial post-processing needed in order to satisfy the principles of stoichiometry. To address these issues, we propose a representation-first approach to molecular graph generation. We guide the latent representation of an autoencoder by capturing graph structure information with the geometric scattering transform and apply penalties that structure the representation also by molecular properties. We show that this highly structured latent space can be directly used for molecular graph generation by the use of a GAN. We demonstrate that our architecture learns meaningful representations of drug datasets and provides a platform for goal-directed drug synthesis.


title: Widely-tunable, doubly-resonant Raman scattering on diamond in an open microcavity

abstract: Raman lasers based on bulk diamond are a valuable resource for generating coherent light in wavelength regimes where no common laser diodes are available. Nevertheless, the widespread use of such lasers is limited by their high threshold power requirements on the order of several Watts. Using on-chip microresonators, a significant reduction of the lasing threshold by more than two orders of magnitude has been shown. However, these resonators lack a continuous tuning mechanism and, mainly due to fabrication limitations, their implementation in the visible remains elusive. Here, we propose a platform for a diamond Raman laser in the visible. The device is based on a miniaturized, open-access Fabry-Perot cavity. Our microcavity provides widely-tunable doubly-resonant enhancement of Raman scattering from high quality single-crystalline diamond. We demonstrate a $>$THz continuous tuning range of doubly-resonant Raman scattering, a range limited only by the reflective stopband of the mirrors. Based on the experimentally determined quality factors exceeding $300\,000$, our theoretical analysis suggests that, with realistic improvements, a sub-mW threshold is readily within reach. Our findings pave the way to the creation of a universal low-power frequency shifter, a potentially valuable addition to the nonlinear optics toolbox.


title: Witnessing Objectivity on a Quantum Computer

abstract: Understanding the emergence of objectivity from the quantum realm has been a long standing issue strongly related to the quantum to classical crossover. Quantum Darwinism provides an answer, interpreting objectivity as consensus between independent observers. Quantum computers provide an interesting platform for such experimental investigation of quantum Darwinism, fulfilling their initial intended purpose as quantum simulators. Here we assess to what degree current NISQ devices can be used as experimental platforms in the field of quantum Darwinism. We do this by simulating an exactly solvable stochastic collision model, taking advantage of the analytical solution to benchmark the experimental results.


title: Congruence properties of combinatorial sequences via Walnut and the Rowland-Yassawi-Zeilberger automaton

abstract: Certain famous combinatorial sequences, such as the Catalan numbers and the Motzkin numbers, when taken modulo a prime power, can be computed by finite automata. Many theorems about such sequences can therefore be proved using Walnut, which is an implementation of a decision procedure for proving various properties of automatic sequences. In this paper we explore some results (old and new) that can be proved using this method.


title: Exploring The Holographic Swampland

abstract: We extend studies of holographic aspects of moduli stabilisation scenarios to both fibred versions of LVS and the type IIA DGKT flux vacua. We study the holographic properties of the low-energy moduli Lagrangian that describes both the AdS vacuum and also small perturbations about it. For type IIA vacua in the large-volume regime, the CFT data (operator dimensions and higher-point interactions) take a universal form independent of the many arbitrary flux choices, as was previously found for LVS stabilisation. For these IIA vacua the conformal dimensions of the dual operators are also, surprisingly, all integers, although we do not understand a deeper reason why this is so. In contrast to behaviour previously found for LVS and KKLT, the fibred models also admit cases of mixed double-trace operators (for two different axion fields) where the anomalous dimension is positive.


title: The Early-type Stars from LAMOST survey: Atmospheric parameters

abstract: Massive stars play key roles in many astrophysical processes. Deriving atmospheric parameters of massive stars is important to understand their physical properties and thus are key inputs to trace their evolution. Here we report our work on adopting the data-driven technique Stellar LAbel Machine ({\tt SLAM}) with the non-LTE TLUSTY synthetic spectra as the training dataset to estimate the stellar parameters of LAMOST optical spectra for early-type stars. We apply two consistency tests to verify this machine learning method and compare stellar labels given by {\tt SLAM} with that in literature for several objects having high-resolution spectra. We provide the stellar labels of effective temperature ($T_\mathrm{eff}$), surface gravity ($\log{g}$), metallicity ([M/H]), and projected rotational velocity ($v\sin{i}$) for 3,931 and 578 early-type stars from LAMOST Low-Resolution Survey (LAMOST-LRS) and Medium-Resolution Survey (LAMOST-MRS), respectively. To estimate the average statistical uncertainties of our results, we calculated the standard deviation between the predicted stellar label and the pre-labeled published values from the high-resolution spectra. The uncertainties of the four parameters are $\sigma(T_\mathrm{eff}) = 2,185 $K, $\sigma(\log{g}) = 0.29$ dex, and $\sigma(v\sin{i}) = 11\, \rm km\,s^{-1}$ for MRS, and $\sigma(T_\mathrm{eff}) = 1,642 $K, $\sigma(\log{g}) = 0.25$ dex, and $\sigma(v\sin{i}) = 42\, \rm km\,s^{-1}$ for LRS spectra, respectively. We notice that parameters of $T_\mathrm{eff}$, $\log{g}$ and [M/H] can be better constrained using LRS spectra rather than using MRS spectra, most likely due to their broad wavelength coverage, while $v\sin{i}$ is constrained better by MRS spectra than by LRS spectra, probably due to the relatively accurate line profiles of MRS spectra.


title: Infrared signatures of quantum bounce in a minisuperspace analysis of Lema\^{\i}tre-Tolman-Bondi dust collapse

abstract: In the quantum treatment of collapsing matter, a viable mode analysis is supposed to give information regarding emission during the collapse. Nevertheless, partly owing to operator ordering ambiguities involved in a typical quantum gravity analysis, the availability of such well-defined modes is not guaranteed. We study the mode decomposition of the unitarily evolving wave packet constructed for the quantum model of spherically symmetric dust collapsing in a marginally bound Lema\^{\i}tre-Tolman-Bondi (LTB) model. We identify the observable depicting mode decomposition, and using the freedom of operator ordering ambiguity, we obtain the Hermitian extension of this operator alongside the Hermitian Hamiltonian. After identifying incoming and outgoing modes with this operator's eigenstates, we estimate their contributions to the radiation profile. The infrared sector of this process demonstrates some characteristic features which turn out to be highly sensitive to the near-bounce dynamics of the dust cloud. Near the epoch of classical singularity, there is a significant contribution from incoming/outgoing modes of small wavenumber in the expanding/collapsing phase of the dust cloud, which keeps on decreasing as one moves away from the singularity. The information of the bounce is carried over to the infrared modes through a flip from largely incoming to largely outgoing radiation as the evolution progresses from collapsing to expanding phase, much before the information of bounce comes about to any observer. In the infrared sector, the saturation value of the amplitude marks the bounce radius. Thus, we argue that the information of the short scale physics is essentially carried over to the longest wavelength in this quantum gravity model, which we argue is rather more prominent for low energy processes.


title: Streaming Instability with Multiple Dust Species: II. Turbulence and Dust-Gas Dynamics at Nonlinear Saturation

abstract: The streaming instability is a fundamental process that can drive dust-gas dynamics and ultimately planetesimal formation in protoplanetary discs. As a linear instability, it has been shown that its growth with a distribution of dust sizes can be classified into two distinct regimes, fast- and slow-growth, depending on the dust-size distribution and the total dust-to-gas density ratio $\epsilon$. Using numerical simulations of an unstratified disc, we bring three cases in different regimes into nonlinear saturation. We find that the saturation states of the two fast-growth cases are similar to its single-species counterparts. The one with maximum dimensionless stopping time $\tau_\mathrm{s,max}=0.1$ and $\epsilon=2$ drives turbulent vertical dust-gas vortices, while the other with $\tau_\mathrm{s,max}=2$ and $\epsilon=0.2$ leads to radial traffic jams and filamentary structures of dust particles. The dust density distribution for the former is flat in low densities, while the one for the latter has a low-end cutoff. By contrast, the one slow-growth case results in a virtually quiescent state. Moreover, we find that in the fast-growth regime, significant dust segregation by size occurs, with large particles moving towards dense regions while small particles remain in the diffuse regions, and the mean radial drift of each dust species is appreciably altered from the (initial) drag-force equilibrium. The former effect may skew the spectral index derived from multi-wavelength observations and change the initial size distribution of a pebble cloud for planetesimal formation. The latter along with turbulent diffusion may influence the radial transport and mixing of solid materials in young protoplanetary discs.


title: The Grand Canonical Multiverse and the Small Cosmological Constant

abstract: We consider the Multiverse as an ensemble of universes. Using standard statistical physics analysis we get that the Cosmological Constant (CC) is exponentially small. The small and finite CC is achieved without any anthropic reasoning. We then quantize the CC. The quantization allows a precise summation of the possible contributions and using the measured value of the CC yields a prediction on the temperature of the Multiverse that we define. Furthermore, quantization allows the interpretation of a single Universe as a superposition of different eigenstates with different energy levels rather than the existence of an actual Multiverse.


title: On Permutation Invariant Problems in Large-Scale Inference

abstract: Simultaneous statistical inference problems are at the basis of almost any scientific discovery process. We consider a class of simultaneous inference problems that are invariant under permutations, meaning that all components of the problem are oblivious to the labelling of the multiple instances under consideration. For any such problem we identify the optimal solution which is itself permutation invariant, the most natural condition one could impose on the set of candidate solutions. Interpreted differently, for any possible value of the parameter we find a tight (non-asymptotic) lower bound on the statistical performance of any procedure that obeys the aforementioned condition. By generalizing the standard decision theoretic notions of permutation invariance, we show that the results apply to a myriad of popular problems in simultaneous inference, so that the ultimate benchmark for each of these problems is identified. The connection to the nonparametric empirical Bayes approach of Robbins is discussed in the context of asymptotic attainability of the bound uniformly in the parameter value.


title: Fifty Years of 3He-rich Events

abstract: The early 1970s saw a new and surprising feature in the composition of solar energetic particles (SEPs), resonant enhancements up to 10,000-fold in the ratio 3He/4He that could even make 3He dominant over H in rare events. It was soon learned that these events also had enhancements in the abundances of heavier elements, such as a factor of ~10 enhancements in Fe/O, which was later seen to be part of a smooth increase in enhancements vs. mass-to-charge ratio A/Q from H to Pb, rising by a factor of ~1000. These events were also associated with streaming 10 - 100 keV electrons that produce type III radio bursts. In recent years we have found these "impulsive" SEP events to be accelerated in islands of magnetic reconnection from plasma temperatures of 2 - 3 MK on open field lines in solar jets. Similar reconnection on closed loops traps the energy of the particles to produce hot (>10 MK), bright flares. Sometimes impulsive SEP intensities are boosted by shock waves when the jets launch fast coronal mass ejections. No single theory yet explains both the sharp resonance in 3He and the smooth increase up to heavier elements; two processes seem to occur. Sometimes the efficient acceleration even exhausts the rare 3He in the source region, limiting its fluence.


title: Two-dimensional gravity from vanishing metrical dimensions

abstract: We obtain a dynamical formulation of two-dimensional gravity from a non-Einsteinian phase in higher dimensions $(D=3+2n)$. The formalism is associated with (at least) one extra dimension of vanishing proper length, thus being inequivalent to either a Kaluza-Klein compactification or the Mann-Ross dimensional reduction defined upon a singular limit. The emergent solutions admit any arbitrary curvature in contrast with Jackiw-Teitelboim constant curvature gravity. We present the static and homogeneous solutions as explicit examples. The effective field equations are shown to remain unaffected by the inclusion of higher Lovelock terms beyond Einstein.


title: StateAFL: Greybox Fuzzing for Stateful Network Servers

abstract: Fuzzing network servers is a technical challenge, since the behavior of the target server depends on its state over a sequence of multiple messages. Existing solutions are costly and difficult to use, as they rely on manually-customized artifacts such as protocol models, protocol parsers, and learning frameworks. The aim of this work is to develop a greybox fuzzer (StateaAFL) for network servers that only relies on lightweight analysis of the target program, with no manual customization, in a similar way to what the AFL fuzzer achieved for stateless programs. The proposed fuzzer instruments the target server at compile-time, to insert probes on memory allocations and network I/O operations. At run-time, it infers the current protocol state of the target server by taking snapshots of long-lived memory areas, and by applying a fuzzy hashing algorithm (Locality-Sensitive Hashing) to map memory contents to a unique state identifier. The fuzzer incrementally builds a protocol state machine for guiding fuzzing. We implemented and released StateaAFL as open-source software. As a basis for reproducible experimentation, we integrated StateaAFL with a large set of network servers for popular protocols, with no manual customization to accomodate for the protocol. The experimental results show that the fuzzer can be applied with no manual customization on a large set of network servers for popular protocols, and that it can achieve comparable, or even better code coverage and bug detection than customized fuzzing. Moreover, our qualitative analysis shows that states inferred from memory better reflect the server behavior than only using response codes from messages.


title: Linear Relations of Siegel Poincare Series and Nonvanshing of the Central Value of Spinor L-functions

abstract: In this paper, we will first investigate the linear relations of a one parameter family of Siegel Poincar\'e series. Then we give the applications to the non-vanishing of Fourier coefficients of Siegel cusp eigenforms and the central values.


title: Not all noise is accounted equally: How differentially private learning benefits from large sampling rates

abstract: Learning often involves sensitive data and as such, privacy preserving extensions to Stochastic Gradient Descent (SGD) and other machine learning algorithms have been developed using the definitions of Differential Privacy (DP). In differentially private SGD, the gradients computed at each training iteration are subject to two different types of noise. Firstly, inherent sampling noise arising from the use of minibatches. Secondly, additive Gaussian noise from the underlying mechanisms that introduce privacy. In this study, we show that these two types of noise are equivalent in their effect on the utility of private neural networks, however they are not accounted for equally in the privacy budget. Given this observation, we propose a training paradigm that shifts the proportions of noise towards less inherent and more additive noise, such that more of the overall noise can be accounted for in the privacy budget. With this paradigm, we are able to improve on the state-of-the-art in the privacy/utility tradeoff of private end-to-end CNNs.


title: Neural Network Weights Do Not Converge to Stationary Points: An Invariant Measure Perspective

abstract: This work examines the deep disconnect between existing theoretical analyses of gradient-based algorithms and the practice of training deep neural networks. Specifically, we provide numerical evidence that in large-scale neural network training (e.g., ImageNet + ResNet101, and WT103 + TransformerXL models), the neural network's weights do not converge to stationary points where the gradient of the loss is zero. Remarkably, however, we observe that even though the weights do not converge to stationary points, the progress in minimizing the loss function halts and training loss stabilizes. Inspired by this observation, we propose a new perspective based on ergodic theory of dynamical systems to explain it. Rather than studying the evolution of weights, we study the evolution of the distribution of weights. We prove convergence of the distribution of weights to an approximate invariant measure, thereby explaining how the training loss can stabilize without weights necessarily converging to stationary points. We further discuss how this perspective can better align optimization theory with empirical observations in machine learning practice.


title: Causal Discovery from Conditionally Stationary Time Series

abstract: Causal discovery, i.e., inferring underlying causal relationships from observational data, has been shown to be highly challenging for AI systems. In time series modeling context, traditional causal discovery methods mainly consider constrained scenarios with fully observed variables and/or data from stationary time-series. We develop a causal discovery approach to handle a wide class of non-stationary time-series that are conditionally stationary, where the non-stationary behaviour is modeled as stationarity conditioned on a set of (possibly hidden) state variables. Named State-Dependent Causal Inference (SDCI), our approach is able to recover the underlying causal dependencies, provably with fully-observed states and empirically with hidden states. The latter is confirmed by experiments on synthetic linear system and nonlinear particle interaction data, where SDCI achieves superior performance over baseline causal discovery methods. Improved results over non-causal RNNs on modeling NBA player movements demonstrate the potential of our method and motivate the use of causality-driven methods for forecasting.


title: Cancelling the vacuum energy and Weyl anomaly in the standard model with dimension-zero scalar fields

abstract: The standard model is a remarkably consistent and complete quantum field theory but its coupling to gravity and the Higgs field remain problematic, as reflected in the cosmological constant problem, the Weyl anomaly, and the hierarchy puzzle. We point out that 36 conformally-coupled dimension-zero scalar fields can simultaneously cancel the vacuum energy and both terms in the Weyl anomaly, if the Higgs field is emergent. The cancellation is highly non-trivial: given the standard model gauge group $SU(3)\times SU(2)\times U(1)$, it requires precisely $48$ Weyl fermions, {\it i.e.}, three generations of standard model fermions, including right-handed neutrinos. Due to a large additional gauge symmetry, the new scalars contribute no new local degrees of freedom or particle states. Their only physical state is their vacuum state, in which they possess a scale invariant power spectrum extending to long wavelengths. This suggests a new explanation for the primordial scalar perturbations in cosmology, not requiring inflation. We also discuss how the Higgs field might emerge as a composite object.


title: Long GMC Lifetimes: Using the Method and Data of Meidt et al. 2015 with A Correction

abstract: Meidt et al. 2015 derived short lifetimes of 20-30Myr for giant molecular clouds (GMCs) in M51. Their novel approach utilizes a decline of the GMC population during their inter-arm passage from one spiral arm to the next. Using the inter-arm travel time "t_travel" as a fiducial clock, they converted the decline rate to GMC lifetimes. They implicitly adopted zero as the pattern speed of the spiral arms, resulting in a very short t_travel. However, it is well established that the spiral arm pattern in M51 is rotating (Meidt et al. 2008, 2013), and that t_travel should be measured with respect to the rotating spiral pattern. Here we use the same method and data of Meidt et al. 2015 and re-evaluate GMC lifetimes by accounting for the pattern speed given by Meidt et al. 2013. This correction gives a longer t_travel and consequently longer GMC lifetimes of 60-500Myr.


title: On Kitaoka's conjecture and lifting problem for universal quadratic forms

abstract: For a totally positive definite quadratic form over the ring of integers of a totally real number field $K$, we show that there are only finitely many totally real field extensions of $K$ of a fixed degree over which the form is universal (namely, those that have a short basis in a suitable sense). Along the way we give a general construction of a universal form of rank bounded by $D(\log D)^{d-1}$, where $d$ is the degree of $K$ over $\mathbb Q$ and $D$ is its discriminant. Furthermore, for any fixed degree we prove (weak) Kitaoka's conjecture that there are only finitely many totally real number fields with a universal ternary quadratic form.


title: On Deligne's conjecture for symmetric sixth $L$-functions of Hilbert modular forms

abstract: In this paper, we prove Deligne's conjecture for symmetric sixth $L$-functions of Hilbert modular forms. We extend the result of Morimoto based on a different approach. We define automorphic periods associated to globally generic $C$-algebraic cuspidal automorphic representations of ${\rm GSp}_4$ over totally real number fields whose archimedean components are (limits of) discrete series representations. We show that the algebraicity of critical $L$-values for ${\rm GSp}_4 \times {\rm GL}_2$ can be expressed in terms of these periods. In the case of Kim-Ramakrishnan-Shahidi lifts of ${\rm GL}_2$, we establish period relations between the automorphic periods and powers of Petersson norm of Hilbert modular forms. The conjecture for symmetric sixth $L$-functions then follows from these period relations and our previous work on the algebraicity of critical values for the adjoint $L$-functions for ${\rm GSp}_4$.