title: Spectral Gap of Replica Exchange Langevin Diffusion on Mixture Distributions

abstract: Langevin diffusion (LD) is one of the main workhorses for sampling problems. However, its convergence rate can be significantly reduced if the target distribution is a mixture of multiple densities, especially when each component concentrates around a different mode. Replica exchange Langevin diffusion (ReLD) is a sampling method that can circumvent this issue. In particular, ReLD adds another LD sampling a high-temperature version of the target density, and exchange the locations of two LDs according to a Metropolis-Hasting type of law. This approach can be further extended to multiple replica exchange Langevin diffusion (mReLD), where $K$ additional LDs are added to sample distributions at different temperatures and exchanges take place between neighboring-temperature processes. While ReLD and mReLD have been used extensively in statistical physics, molecular dynamics, and other applications, there is little existing analysis on its convergence rate and choices of temperatures. This paper closes these gaps assuming the target distribution is a mixture of log-concave densities. We show ReLD can obtain constant or even better convergence rates even when the density components of the mixture concentrate around isolated modes. We also show using mReLD with $K$ additional LDs can achieve the same result while the exchange frequency only needs to be $(1/K)$-th power of the one in ReLD.


title: Learning Hamiltonian Monte Carlo in R

abstract: Hamiltonian Monte Carlo (HMC) is a powerful tool for Bayesian computation. In comparison with the traditional Metropolis-Hastings algorithm, HMC offers greater computational efficiency, especially in higher dimensional or more complex modeling situations. To most statisticians, however, the idea of HMC comes from a less familiar origin, one that is based on the theory of classical mechanics. Its implementation, either through Stan or one of its derivative programs, can appear opaque to beginners. A lack of understanding of the inner working of HMC, in our opinion, has hindered its application to a broader range of statistical problems. In this article, we review the basic concepts of HMC in a language that is more familiar to statisticians, and we describe an HMC implementation in R, one of the most frequently used statistical software environments. We also present hmclearn, an R package for learning HMC. This package contains a general-purpose HMC function for data analysis. We illustrate the use of this package in common statistical models. In doing so, we hope to promote this powerful computational tool for wider use. Example code for common statistical models is presented as supplementary material for online publication.


title: Cranks for Ramanujan-type congruences of $k$-colored partitions

abstract: Dyson famously provided combinatorial explanations for Ramanujan's partition congruences modulo $5$ and $7$ via his rank function, and postulated that an invariant explaining all of Ramanujan's congruences modulo $5$, $7$, and $11$ should exist. Garvan and Andrews-Garvan later discovered such an invariant called the crank, fulfilling Dyson's goal. Many further examples of congruences of partition functions are known in the literature. In this paper, we provide a framework for discovering and proving the existence of such invariants for families of congruences and partition functions. As a first example, we find a family of crank functions that simultaneously explains most known congruences for colored partition functions. The key insight is to utilize a powerful recent theory of theta blocks due to Gritsenko, Skoruppa, and Zagier. The method used here should be useful in the study of other combinatorial functions.


title: A bound on the degree of singular vectors for the exceptional Lie superalgebra $E(5,10)$

abstract: We use the language of Lie pseudoalgebras to gain information about the representation theory of the simple infinite-dimensional linearly compact Lie superalgebra of exceptional type $E(5,10)$. This technology allows us to prove that the degree of singular vectors in minimal Verma modules is $\leq 14$. A few technical adjustments allow us to refine the bound, proving that the degree must always be $\leq 12$ and it is actually, except for a finite number of cases, $\leq 10$.


title: Supremum, infimum and hyperlimits in the non-Archimedean ring of Colombeau generalized numbers

abstract: It is well-known that the notion of limit in the sharp topology of sequences of Colombeau generalized numbers $\widetilde{\mathbb{R}}$ does not generalize classical results. E.g.~the sequence $\frac{1}{n}\not\to0$ and a sequence $(x_{n})_{n\in\mathbb{N}}$ converges \emph{if} and only if $x_{n+1}-x_{n}\to0$. This has several deep consequences, e.g.~in the study of series, analytic generalized functions, or sigma-additivity and classical limit theorems in integration of generalized functions. The lacking of these results is also connected to the fact that $\widetilde{\mathbb{R}}$ is necessarily not a complete ordered set, e.g.~the set of all the infinitesimals does not have neither supremum nor infimum. We present a solution of these problems with the introduction of the notions of hypernatural number, hypersequence, close supremum and infimum. In this way, we can generalize all the classical theorems for the hyperlimit of a hypersequence. The paper explores ideas that can be applied to other non-Archimedean settings.


title: A closed-form expression for the kth power of semicirculant and r-circulant matrices

abstract: We derive a closed-form expression for the kth power of semicirculant matrices by using the determinant of certain matrices. As an application, a closed-form expression for the kth power of r-circulant matrices is also povided.


title: Carleman Estimate for Ultrahyperbolic Operators and Improved Interior Control for Wave Equations

abstract: In this article, we present a novel Carleman estimate for ultrahyperbolic operators, in $ \mathbb{R}^m_t \times \mathbb{R}^n_x $. Then, we use a special case of this estimate to obtain improved observability results for wave equations with time-dependent lower order terms. The key improvements are: (1) we obtain smaller observation regions compared to standard Carleman estimate results, and (2) we also prove observability when the observation point lies inside the domain. Finally, as a corollary of the observability result, we obtain improved interior controllability for the wave equation.


title: Near-Optimal SQ Lower Bounds for Agnostically Learning Halfspaces and ReLUs under Gaussian Marginals

abstract: We study the fundamental problems of agnostically learning halfspaces and ReLUs under Gaussian marginals. In the former problem, given labeled examples $(\mathbf{x}, y)$ from an unknown distribution on $\mathbb{R}^d \times \{ \pm 1\}$, whose marginal distribution on $\mathbf{x}$ is the standard Gaussian and the labels $y$ can be arbitrary, the goal is to output a hypothesis with 0-1 loss $\mathrm{OPT}+\epsilon$, where $\mathrm{OPT}$ is the 0-1 loss of the best-fitting halfspace. In the latter problem, given labeled examples $(\mathbf{x}, y)$ from an unknown distribution on $\mathbb{R}^d \times \mathbb{R}$, whose marginal distribution on $\mathbf{x}$ is the standard Gaussian and the labels $y$ can be arbitrary, the goal is to output a hypothesis with square loss $\mathrm{OPT}+\epsilon$, where $\mathrm{OPT}$ is the square loss of the best-fitting ReLU. We prove Statistical Query (SQ) lower bounds of $d^{\mathrm{poly}(1/\epsilon)}$ for both of these problems. Our SQ lower bounds provide strong evidence that current upper bounds for these tasks are essentially best possible.


title: Graph-theoretic optimization for edge consensus

abstract: We consider network structures that optimize the $\mathcal{H}_2$ norm of weighted, time scaled consensus networks, under a minimal representation of such consensus networks described by the edge Laplacian. We show that a greedy algorithm can be used to find the minimum-$\mathcal{H}_2$ norm spanning tree, as well as how to choose edges to optimize the $\mathcal{H}_2$ norm when edges are added back to a spanning tree. In the case of edge consensus with a measurement model considering all edges in the graph, we show that adding edges between slow nodes in the graph provides the smallest increase in the $\mathcal{H}_2$ norm.


title: Partitioned Least Squares

abstract: In this paper we propose a variant of the linear least squares model allowing practitioners to partition the input features into groups of variables that they require to contribute similarly to the final result. The output allows practitioners to assess the importance of each group and of each variable in the group. We formally show that the new formulation is not convex and provide two alternative methods to deal with the problem: one non-exact method based on an alternating least squares approach; and one exact method based on a reformulation of the problem using an exponential number of sub-problems whose minimum is guaranteed to be the optimal solution. We formally show the correctness of the exact method and also compare the two solutions showing that the exact solution provides better results in a fraction of the time required by the alternating least squares solution (assuming that the number of partitions is small). For the sake of completeness, we also provide an alternative branch and bound algorithm that can be used in place of the exact method when the number of partitions is too large, and a proof of NP-completeness of the optimization problem introduced in this paper.


title: How does the extracellular matrix affect the rigidity of an embedded spheroid?

abstract: Cellularized tissue and polymer networks can both transition from floppy to rigid as a function of their control parameters, and, yet, the two systems often mechanically interact, which may affect their respective rigidities. To study this interaction, we consider a vertex model with interfacial tension (a spheroid) embedded in a spring network in two dimensions. We identify two regimes with different global spheroid shapes, governed by the pressure resulting from competition between interfacial tension and tension in the network. In the first regime, the tissue remains compact, while in the second, a cavitation-like instability leads to the emergence of gaps at the tissue-network interface. Intriguingly, compression of the tissue promotes fluidization, while tension promotes cellular alignment and rigidification, with the mechanisms driving rigidification differing on either side of the instability.


title: Coloured noise time series as appropriate models for environmental variation in artificial evolutionary systems

abstract: Ecological, environmental and geophysical time series consistently exhibit the characteristics of coloured (1/f^\b{eta}) noise. Here we briefly survey the literature on coloured noise, population persistence and related evolutionary dynamics, before introducing coloured noise as an appropriate model for environmental variation in artificial evolutionary systems. To illustrate and explore the effects of different noise colours, a simple evolutionary model that examines the trade-off between specialism and generalism in fluctuating environments is applied. The results of the model clearly demonstrate a need for greater generalism as environmental variability becomes `whiter', whilst specialisation is favoured as environmental variability becomes `redder'. Pink noise, sitting midway between white and red noise, is shown to be the point at which the pressures for generalism and specialism balance, providing some insight in to why `pinker' noise is increasingly being seen as an appropriate model of typical environmental variability. We go on to discuss how the results presented here feed in to a wider discussion on evolutionary responses to fluctuating environments. Ultimately we argue that Artificial Life as a field should embrace the use of coloured noise to produce models of environmental variability.


title: Composed Fine-Tuning: Freezing Pre-Trained Denoising Autoencoders for Improved Generalization

abstract: We focus on prediction problems with structured outputs that are subject to output validity constraints, e.g. pseudocode-to-code translation where the code must compile. While labeled input-output pairs are expensive to obtain, "unlabeled" outputs, i.e. outputs without corresponding inputs, are freely available (e.g. code on GitHub) and provide information about output validity. We can capture the output structure by pre-training a denoiser to denoise corrupted versions of unlabeled outputs. We first show that standard fine-tuning after pre-training destroys some of this structure. We then propose composed fine-tuning, which fine-tunes a predictor composed with the pre-trained denoiser, which is frozen to preserve output structure. For two-layer ReLU networks, we prove that composed fine-tuning significantly reduces the complexity of the predictor, thus improving generalization. Empirically, we show that composed fine-tuning improves over standard fine-tuning on two pseudocode-to-code translation datasets (3% and 6% relative). The improvement from composed fine-tuning is magnified on out-of-distribution (OOD) examples (4% and 25% relative).


title: Reputation for Playing Mixed Actions: A Characterization Theorem

abstract: A patient player privately observes a persistent state that directly affects his myopic opponents' payoffs, and can be one of the several commitment types that plays the same mixed action in every period. I characterize the set of environments under which the patient player obtains at least his commitment payoff in all equilibria regardless of his stage-game payoff function. Due to interdependent values, the patient player cannot guarantee his mixed commitment payoff by imitating the mixed-strategy commitment type, and small perturbations to a pure commitment action can significantly reduce the patient player's guaranteed equilibrium payoff.


title: Enhanced space charge limited current for curved electron emitters

abstract: The maximum current that can be transported across a vacuum diode is limited by forces arising due to space charge. In a planar diode configuration, the space charge limited (SCL) current density from a planar emitting patch is given by the Child-Langmuir (CL) law $J_{CL} \sim V_g^{3/2}/D^2$ where $V_g$ is the potential difference across the diode and $D$ is the separation between the anode and cathode. We show here analytically using the nonlinear line charge model that for a curved emitter in a planar diode configuration, the limiting current obeys the scaling relationship $J_{SCL} \sim \gamma_a V_g^{3/2}/D^2$ where $\gamma_a$ is the apex field enhancement factor of the curved emitter. For an emitter with large height ($h$) to apex radius of curvature ($R_a$) ratio, the limiting current far exceeds the planar value. The result is verified using the particle-in-cell code PASUPAT for two curved emitters shapes.

categories: physics.plasm-ph cond-mat.mes-hall physics.acc-ph physics.app-ph


title: Sensitivity of stellar physics to the equation of state

abstract: The formation and evolution of stars depends on various physical aspects of stellar matter, including the equation of state (EOS) and transport properties. Although often dismissed as `ideal gas-like' and therefore simple, states occurring in stellar matter are dense plasmas, and the EOS has not been established precisely. EOS constructed using multi-physics approaches found necessary for laboratory studies of warm dense matter give significant variations in stellar regimes, and vary from the EOS commonly used in simulations of the formation and evolution of stars. We have investigated the sensitivity of such simulations to variations in the EOS, for sun-like and low-mass stars, and found a strong sensitivity of the lifetime of the Sun and of the lower luminosity limit for red dwarfs. We also found a significant sensitivity in the lower mass limit for red dwarfs. Simulations of this type are also used for other purposes in astrophysics, including the interpretation of absolute magnitude as mass, the conversion of inferred mass distribution to the initial mass function using predicted lifetimes, simulations of star formation from nebulae, simulations of galactic evolution, and the baryon census used to bound the exotic contribution to dark matter. Although many of the sensitivities of stellar physics to the EOS are large, some of the inferred astrophysical quantities are also constrained by independent measurements, although the constraints may be indirect and non-trivial. However, it may be possible to use such measurements to constrain the EOS more than presently possible by established plasma theory.


title: Synchronisation phase as an indicator of persistent quantum correlations between subsystems

abstract: Spontaneous synchronisation is a collective phenomenon that can occur in both dynamical classical and quantum systems. Here, we analyse the spontaneous synchronisation dynamics of vibrations assisting energy transfer in a bio-inspired system. We find the emergence of a constant non-zero `synchronisation phase' between synchronised vibrational displacements as the natural frequencies of the oscillators are detuned. This phase difference arises from the asymmetric participation of local modes in the long-lived synchronised state. Furthermore, we investigate the relationships between the synchronisation phase, detuning and the degree of quantum correlations between the synchronising subsystems and find that the synchronisation phase captures how quantum correlations persistently exceed classical correlations during the dynamics. We show that our analysis applies to a variety of spontaneously synchronising open quantum systems. Our work therefore opens up a promising avenue to investigate non-trivial quantum phenomena in complex biomolecular and nano-scale chemical systems.


title: Model-based Reinforcement Learning for Semi-Markov Decision Processes with Neural ODEs

abstract: We present two elegant solutions for modeling continuous-time dynamics, in a novel model-based reinforcement learning (RL) framework for semi-Markov decision processes (SMDPs), using neural ordinary differential equations (ODEs). Our models accurately characterize continuous-time dynamics and enable us to develop high-performing policies using a small amount of data. We also develop a model-based approach for optimizing time schedules to reduce interaction rates with the environment while maintaining the near-optimal performance, which is not possible for model-free methods. We experimentally demonstrate the efficacy of our methods across various continuous-time domains.


title: Enhanced photocatalytic activity of plasmonic Au nanoparticles incorporated MoS$_2$ nanosheets for degradation of organic dyes

abstract: In the present paper, we investigate the effect of plasmonic Au nanoparticles (NPs) decoration on the photocatalytic efficiency of MoS$_2$ nanosheets. The Au NPs are grown on the surface of chemically exfoliated MoS$_2$ nanosheets by chemical reduction method. Au-MoS$_2$ nanostructures (NSs) are characterized by X-ray diffractometer, Raman spectrometer, absorption spectrophotometer, and transmission electron microscopy. Exfoliated MoS$_2$ and Au-MoS$_2$ NSs are used to study the photocatalytic degradation of organic dyes methyl red (MR) and methylene blue (MB). Under UV-Visible light irradiation, pristine MoS$_2$ shows photo degradation efficiencies between 30% to 46.9% for MR and 23.3% to 44% for MB, with varying exposure times from 30 to 120 min, respectively. However, Au-MoS$_2$ NSs with maximum Au NPs concentration show enhanced degradation efficiency from 70.2 to 96.7% for MR, and from 65.2 to 94.3% for MB. The manifold enhancement of degradation efficiency for both the dyes with Au-MoS$_2$ NSs may be attributed to the presence of Au NPs acting as charge trapping sites in the NSs. We believe this study would have potential application in battling the ill-effects of environmental degradation, which poses a major threat to humans as well as biodiversity.


title: Towards the Study of Morphological Processing of the Tangkhul Language

abstract: There is no or little work on natural language processing of Tangkhul language. The current work is a humble beginning of morphological processing of this language using an unsupervised approach. We use a small corpus collected from different sources of text books, short stories and articles of other topics. Based on the experiments carried out, the morpheme identification task using morphessor gives reasonable and interesting output despite using a small corpus.


title: Totally positive kernels, Polya frequency functions, and their transforms

abstract: The composition operators preserving total non-negativity and total positivity for various classes of kernels are classified, following three themes. Letting a function act by post composition on kernels with arbitrary domains, it is shown that such a composition operator maps the set of totally non-negative kernels to itself if and only if the function is constant or linear, or just linear if it preserves total positivity. Symmetric kernels are also discussed, with a similar outcome. These classification results are a byproduct of two matrix-completion results and the second theme: an extension of A.M. Whitney's density theorem from finite domains to subsets of the real line. This extension is derived via a discrete convolution with modulated Gaussian kernels. The third theme consists of analyzing, with tools from harmonic analysis, the preservers of several families of totally non-negative and totally positive kernels with additional structure: continuous Hankel kernels on an interval, P\'olya frequency functions, and P\'olya frequency sequences. The rigid structure of post-composition transforms of totally positive kernels acting on infinite sets is obtained by combining several specialized situations settled in our present and earlier works.


title: Estimation of Covid-19 Prevalence from Serology Tests: A Partial Identification Approach

abstract: We propose a partial identification method for estimating disease prevalence from serology studies. Our data are results from antibody tests in some population sample, where the test parameters, such as the true/false positive rates, are unknown. Our method scans the entire parameter space, and rejects parameter values using the joint data density as the test statistic. The proposed method is conservative for marginal inference, in general, but its key advantage over more standard approaches is that it is valid in finite samples even when the underlying model is not point identified. Moreover, our method requires only independence of serology test results, and does not rely on asymptotic arguments, normality assumptions, or other approximations. We use recent Covid-19 serology studies in the US, and show that the parameter confidence set is generally wide, and cannot support definite conclusions. Specifically, recent serology studies from California suggest a prevalence anywhere in the range 0%-2% (at the time of study), and are therefore inconclusive. However, this range could be narrowed down to 0.7%-1.5% if the actual false positive rate of the antibody test was indeed near its empirical estimate (~0.5%). In another study from New York state, Covid-19 prevalence is confidently estimated in the range 13%-17% in mid-April of 2020, which also suggests significant geographic variation in Covid-19 exposure across the US. Combining all datasets yields a 5%-8% prevalence range. Our results overall suggest that serology testing on a massive scale can give crucial information for future policy design, even when such tests are imperfect and their parameters unknown.


title: Tuning Chern Number in Quantum Anomalous Hall Insulators

abstract: The quantum anomalous Hall (QAH) state is a two-dimensional topological insulating state that has quantized Hall resistance of h/Ce2 and vanishing longitudinal resistance under zero magnetic field, where C is called the Chern number. The QAH effect has been realized in magnetic topological insulators (TIs) and magic-angle twisted bilayer graphene. Despite considerable experimental efforts, the zero magnetic field QAH effect has so far been realized only for C = 1. Here we used molecular beam epitaxy to fabricate magnetic TI multilayers and realized the QAH effect with tunable Chern number C up to 5. The Chern number of these QAH insulators is tuned by varying the magnetic doping concentration or the thickness of the interior magnetic TI layers in the multilayer samples. A theoretical model is developed to understand our experimental observations and establish phase diagrams for QAH insulators with tunable Chern numbers. The realization of QAH insulators with high tunable Chern numbers facilitates the potential applications of dissipationless chiral edge currents in energy-efficient electronic devices and opens opportunities for developing multi-channel quantum computing and higher-capacity chiral circuit interconnects.


title: Radiography in high mass X-ray binaries -- Micro-structure of the stellar wind through variability of the column density

abstract: In high mass X-ray binaries (HMXBs), an accreting compact object orbits a high mass star which loses mass through a dense and inhomogeneous wind. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass loss rates for high mass stars. We explore the impact of clumpiness on the variability of the column density with a simplified wind model. In particular, we focus on the standard deviation of the column density and the characteristic duration of enhanced absorption episodes, and compare them with analytical predictions based on the porosity length. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a clump in front of the compact object. We provide a recipe to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. The coherence time scale grants direct access to the size of the clumps while their mass can be deduced separately from the amplitude of the variability. If it is due to unaccreted passing-by clumps, the high column density variations in some HMXBs requires high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are hardly compatible with the ones derived from first principles. Alternatively, other components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism still to be identified, or a dense environment in the immediate vicinity of the accretor such as an accretion disk, an outflow or a spherical shell around the magnetosphere of the accreting neutron star.


title: Extrapolated Sequential Constraint Method for Variational Inequality over the Intersection of Fixed-Point Sets

abstract: This paper deals with the solving of variational inequality problem where the constrained set is given as the intersection of a number of fixed-point sets. To this end, we present an extrapolated sequential constraint method. At each iteration, the proposed method is updated based on the ideas of a hybrid conjugate gradient method used to accelerate the well-known hybrid steepest descent method, and an extrapolated cyclic cutter method for solving a common fixed point problem. We prove strong convergence of the method under some suitable assumptions of step-size sequences. We finally show the numerical efficiency of the proposed method compared to some existing methods.


title: On the Iteration Complexity of Hypergradient Computation

abstract: We study a general class of bilevel problems, consisting in the minimization of an upper-level objective which depends on the solution to a parametric fixed-point equation. Important instances arising in machine learning include hyperparameter optimization, meta-learning, and certain graph and recurrent neural networks. Typically the gradient of the upper-level objective (hypergradient) is hard or even impossible to compute exactly, which has raised the interest in approximation methods. We investigate some popular approaches to compute the hypergradient, based on reverse mode iterative differentiation and approximate implicit differentiation. Under the hypothesis that the fixed point equation is defined by a contraction mapping, we present a unified analysis which allows for the first time to quantitatively compare these methods, providing explicit bounds for their iteration complexity. This analysis suggests a hierarchy in terms of computational efficiency among the above methods, with approximate implicit differentiation based on conjugate gradient performing best. We present an extensive experimental comparison among the methods which confirm the theoretical findings.


title: Griffiths-McCoy singularity on the diluted Chimera graph: Monte Carlo simulations and experiments on the quantum hardware

abstract: The Griffiths-McCoy singularity is a phenomenon characteristic of low-dimensional disordered quantum spin systems, in which the magnetic susceptibility shows singular behavior as a function of the external field even within the paramagnetic phase. We study whether this phenomenon is observed in the transverse-field Ising model with disordered ferromagnetic interactions on the quasi-two-dimensional diluted Chimera graph both by quantum Monte Carlo simulations and by extensive experiments on the D-Wave quantum annealer used as a quantum simulator. From quantum Monte Carlo simulations, evidence is found for the existence of the Griffiths-McCoy singularity in the paramagnetic phase. The experimental approach on the quantum hardware produces results that are less clear-cut due to the intrinsic noise and errors in the analog quantum device but can nonetheless be interpreted to be consistent with the existence of the Griffiths-McCoy singularity as in the Monte Carlo case. This is the first experimental approach based on an analog quantum simulator to study the subtle phenomenon of Griffiths-McCoy singularities in a disordered quantum spin system, through which we have clarified the capabilities and limitations of the D-Wave quantum annealer as a quantum simulator.


title: Systematics in the XENON1T data: the 15-keV anti-axion

abstract: The XENON1T collaboration has found an excess of electron recoil events in their Science Run 1 data below ~7 keV with a spectral shape consistent with that expected from a solar-axion-induced signal. The claimed statistical significance of the solar-axion model over the null hypothesis is 3.5$\sigma$. In this work we provide evidence for mismodeling in the electron recoil data that may decrease the local significance of the axion model to as low as $p \approx 0.1$. To reach this conclusion, we search for a signal with the spectral template of the solar axion model, but shifted to higher (unphysical) energies above ~7 keV. We find that the distribution of significances found from this side-band analysis does not follow the expected chi-square distribution, which allows us to quantify the extent to which mismodeling may be affecting the interpretation of the data at energies below ~7 keV. For example, we find that there is an even higher-significance fit to the data when the solar axion model is shifted upwards in energy by ~15 keV and allowed to have a negative normalization.


title: Distributed Link Removal Strategy for Networked Meta-Population Epidemics and its Application to the Control of the COVID-19 Pandemic

abstract: In this paper, we investigate the distributed link removal strategy for networked meta-population epidemics. In particular, a deterministic networked susceptible-infected-recovered (SIR) model is considered to describe the epidemic evolving process. In order to curb the spread of epidemics, we present the spectrum-based optimization problem involving the Perron-Frobenius eigenvalue of the matrix constructed by the network topology and transition rates. A modified distributed link removal strategy is developed such that it can be applied to the SIR model with heterogeneous transition rates on weighted digraphs. The proposed approach is implemented to control the COVID-19 pandemic by using the reported infected and recovered data in each state of Germany. The numerical experiment shows that the infected percentage can be significantly reduced by using the distributed link removal strategy.


title: Efficient Enumerations for Minimal Multicuts and Multiway Cuts

abstract: Let $G = (V, E)$ be an undirected graph and let $B \subseteq V \times V$ be a set of terminal pairs. A node/edge multicut is a subset of vertices/edges of $G$ whose removal destroys all the paths between every terminal pair in $B$. The problem of computing a {\em minimum} node/edge multicut is NP-hard and extensively studied from several viewpoints. In this paper, we study the problem of enumerating all {\em minimal} node multicuts. We give an incremental polynomial delay enumeration algorithm for minimal node multicuts, which extends an enumeration algorithm due to Khachiyan et al. (Algorithmica, 2008) for minimal edge multicuts. Important special cases of node/edge multicuts are node/edge {\em multiway cuts}, where the set of terminal pairs contains every pair of vertices in some subset $T \subseteq V$, that is, $B = T \times T$. We improve the running time bound for this special case: We devise a polynomial delay and exponential space enumeration algorithm for minimal node multiway cuts and a polynomial delay and space enumeration algorithm for minimal edge multiway cuts.


title: Amplitude estimation via maximum likelihood on noisy quantum computer

abstract: Recently we find several candidates of quantum algorithms that may be implementable in near-term devices for estimating the amplitude of a given quantum state, which is a core sub- routine in various computing tasks such as the Monte Carlo methods. One of those algorithms is based on the maximum likelihood estimate with parallelized quantum circuits. In this paper, we extend this method so that it incorporates the realistic noise effect, and then give an experimental demonstration on a superconducting IBM Quantum device. The maximum likelihood estimator is constructed based on the model assuming the depolarization noise. We then formulate the problem as a two-parameters estimation problem with respect to the target amplitude parameter and the noise parameter. In particular we show that there exist anomalous target values, where the Fisher information matrix becomes degenerate and consequently the estimation error cannot be improved even by increasing the number of amplitude amplifications. The experimental demonstration shows that the proposed maximum likelihood estimator achieves quantum speedup in the number of queries, though the estimation error saturates due to the noise. This saturated value of estimation error is consistent to the theory, which implies the validity of the depolarization noise model and thereby enables us to predict the basic requirement on the hardware components (particularly the gate error) in quantum computers to realize the quantum speedup in the amplitude estimation task.


title: Hadronic light-by-light contribution to $(g-2)_\mu$ from lattice QCD with SU(3) flavor symmetry

abstract: We perform a lattice QCD calculation of the hadronic light-by-light contribution to $(g-2)_\mu$ at the SU(3) flavor-symmetric point $m_\pi=m_K\simeq 420\,$MeV. The representation used is based on coordinate-space perturbation theory, with all QED elements of the relevant Feynman diagrams implemented in continuum, infinite Euclidean space. As a consequence, the effect of using finite lattices to evaluate the QCD four-point function of the electromagnetic current is exponentially suppressed. Thanks to the SU(3)-flavor symmetry, only two topologies of diagrams contribute, the fully connected and the leading disconnected. We show the equivalence in the continuum limit of two methods of computing the connected contribution, and introduce a sparse-grid technique for computing the disconnected contribution. Thanks to our previous calculation of the pion transition form factor, we are able to correct for the residual finite-size effects and extend the tail of the integrand. We test our understanding of finite-size effects by using gauge ensembles differing only by their volume. After a continuum extrapolation based on four lattice spacings, we obtain $a_\mu^{\rm hlbl} = (65.4\pm 4.9 \pm 6.6)\times 10^{-11}$, where the first error results from the uncertainties on the individual gauge ensembles and the second is the systematic error of the continuum extrapolation. Finally, we estimate how this value will change as the light-quark masses are lowered to their physical values.


title: Object Files and Schemata: Factorizing Declarative and Procedural Knowledge in Dynamical Systems

abstract: Modeling a structured, dynamic environment like a video game requires keeping track of the objects and their states declarative knowledge) as well as predicting how objects behave (procedural knowledge). Black-box models with a monolithic hidden state often fail to apply procedural knowledge consistently and uniformly, i.e., they lack systematicity. For example, in a video game, correct prediction of one enemy's trajectory does not ensure correct prediction of another's. We address this issue via an architecture that factorizes declarative and procedural knowledge and that imposes modularity within each form of knowledge. The architecture consists of active modules called object files that maintain the state of a single object and invoke passive external knowledge sources called schemata that prescribe state updates. To use a video game as an illustration, two enemies of the same type will share schemata but will have separate object files to encode their distinct state (e.g., health, position). We propose to use attention to determine which object files to update, the selection of schemata, and the propagation of information between object files. The resulting architecture is a drop-in replacement conforming to the same input-output interface as normal recurrent networks (e.g., LSTM, GRU) yet achieves substantially better generalization on environments that have multiple object tokens of the same type, including a challenging intuitive physics benchmark.


title: On Matrix Consequence (Extended Abstract)

abstract: These results are a contribution to the model theory of matrix consequence. We give a semantic characterization of uniform and couniform consequence relations. These properties have never been treated individually, at least in a semantic manner. We consider these notions from a purely semantic point of view and separately, introducing the notion of a uniform bundle/atlas and that of a couniform class of logical matrices. Then, we show that any uniform bundle defines a uniform consequence; and if a structural consequence is uniform, then its Lindenbaum atlas is uniform. Thus, any structural consequence is uniform if, and only if, it is determined by a uniform bundle/atlas. On the other hand, any couniform set of matrices defines a couniform structural consequence. Also, the Lindenbaum atlas of a couniform structural consequence is couniform. Thus, any structural consequence is couniform if, and only if, it is determined by a couniform bundle/atlas. We then apply these observations to compare structural consequence relations that are defined in different languages when one language is a primitive extension of another. We obtain that for any structural consequence defined in a language having (at least) a denumerable set of sentential variables, if this consequence is uniform and couniform, then it and the \emph{ W\'{o}jcicki's consequence} corresponding to it, which is defined in any primitive extension of the given language, are determined by one and the same atlas which is both uniform and couniform.


title: A Local Singularity Analysis for the Ricci Flow and its Applications to Ricci Flows with Bounded Scalar Curvature

abstract: We develop a refined singularity analysis for the Ricci flow by investigating curvature blow-up rates locally. We first introduce general definitions of Type I and Type II singular points and show that these are indeed the only possible types of singular points. In particular, near any singular point the Riemannian curvature tensor has to blow up at least at a Type I rate, generalising a result of Enders, Topping and the first author that relied on a global Type I assumption. We also prove analogous results for the Ricci tensor, as well as a localised version of Sesum's result, namely that the Ricci curvature must blow up near every singular point of a Ricci flow, again at least at a Type I rate. Finally, we show some applications of the theory to Ricci flows with bounded scalar curvature.


title: Self-Supervised MultiModal Versatile Networks

abstract: Videos are a rich source of multi-modal supervision. In this work, we learn representations using self-supervision by leveraging three modalities naturally present in videos: visual, audio and language streams. To this end, we introduce the notion of a multimodal versatile network -- a network that can ingest multiple modalities and whose representations enable downstream tasks in multiple modalities. In particular, we explore how best to combine the modalities, such that fine-grained representations of the visual and audio modalities can be maintained, whilst also integrating text into a common embedding. Driven by versatility, we also introduce a novel process of deflation, so that the networks can be effortlessly applied to the visual data in the form of video or a static image. We demonstrate how such networks trained on large collections of unlabelled video data can be applied on video, video-text, image and audio tasks. Equipped with these representations, we obtain state-of-the-art performance on multiple challenging benchmarks including UCF101, HMDB51, Kinetics600, AudioSet and ESC-50 when compared to previous self-supervised work. Our models are publicly available.


title: A probabilistic mechanism for quark confinement

abstract: The confinement of quarks is one of the enduring mysteries of modern physics. There is a longstanding physics heuristic that confinement is a consequence of `unbroken center symmetry'. This article gives mathematical confirmation of this heuristic, by rigorously defining of center symmetry in lattice gauge theories and proving that a theory is confining when center symmetry is unbroken. Furthermore, a sufficient condition for unbroken center symmetry is given: It is shown that if the center of the gauge group is nontrivial, and correlations decay exponentially under arbitrary boundary conditions, then center symmetry does not break.


title: Constraint Algebra in Bigravity

abstract: The constraint algebra is derived in the second order tetrad Hamiltonian formalism of the bigravity. This is done by a straightforward calculation without involving any insights, implicit functions, and Dirac brackets. The tetrad approach is the only way to present the bigravity action as a linear functional of lapses and shifts, and the Hassan-Rosen transform (characterized as "a complicated redefinition of the shift variable" according to the authors) appears here not as an ansatz but as fixing of a Lagrange multiplier. A comparison of this approach with the other ones is provided.


title: Fast Radio Burst Trains from Magnetar Oscillations

abstract: Quasi-periodic oscillations inferred during rare magnetar giant flare tails were initially interpreted as torsional oscillations of the neutron star (NS) crust, and have been more recently described as global core+crust perturbations. Similar frequencies are also present in high signal-to-noise magnetar short bursts. In magnetars, disturbances of the field are strongly coupled to the NS crust regardless of the triggering mechanism of short bursts. For low-altitude magnetospheric magnetar models of fast radio bursts (FRBs) associated with magnetar short bursts, such as the low-twist model, crustal oscillations may be associated with additional radio bursts in the encompassing short burst event (as recently suggested for SGR 1935+2154). Given the large extragalactic volume probed by wide-field radio transient facilities, this offers the prospect of studying NS crusts leveraging samples far more numerous than galactic high-energy magnetar bursts by studying statistics of sub-burst structure or clustered trains of FRBs. We explore the prospects for distinguishing NS equation of state models with increasingly larger future sets of FRB observations. Lower $l$-number eigenmodes (corresponding to FRB time intervals of $\sim5-50$ ms) are likely less susceptible than high-$l$ modes to confusion by systematic effects associated with the NS crust physics, magnetic field, and damping. They may be more promising in their utility, and also may corroborate models where FRBs arise from mature magnetars. Future observational characterization of such signals can also determine whether they can be employed as cosmological "standard oscillators" to constrain redshift, or can be used to constrain the mass of FRB-producing magnetars when reliable redshifts are available.


title: Learning Robot Skills with Temporal Variational Inference

abstract: In this paper, we address the discovery of robotic options from demonstrations in an unsupervised manner. Specifically, we present a framework to jointly learn low-level control policies and higher-level policies of how to use them from demonstrations of a robot performing various tasks. By representing options as continuous latent variables, we frame the problem of learning these options as latent variable inference. We then present a temporal formulation of variational inference based on a temporal factorization of trajectory likelihoods,that allows us to infer options in an unsupervised manner. We demonstrate the ability of our framework to learn such options across three robotic demonstration datasets.


title: Liquid Resource Types

abstract: This article presents liquid resource types, a technique for automatically verifying the resource consumption of functional programs. Existing resource analysis techniques trade automation for flexibility -- automated techniques are restricted to relatively constrained families of resource bounds, while more expressive proof techniques admitting value-dependent bounds rely on handwritten proofs. Liquid resource types combine the best of these approaches, using logical refinements to automatically prove precise bounds on a program's resource consumption. The type system augments refinement types with potential annotations to conduct an amortized resource analysis. Importantly, users can annotate data structure declarations to indicate how potential is allocated within the type, allowing the system to express bounds with polynomials and exponentials, as well as more precise expressions depending on program values. We prove the soundness of the type system, provide a library of flexible and reusable data structures for conducting resource analysis, and use our prototype implementation to automatically verify resource bounds that previously required a manual proof.


title: True to the Model or True to the Data?

abstract: A variety of recent papers discuss the application of Shapley values, a concept for explaining coalitional games, for feature attribution in machine learning. However, the correct way to connect a machine learning model to a coalitional game has been a source of controversy. The two main approaches that have been proposed differ in the way that they condition on known features, using either (1) an interventional or (2) an observational conditional expectation. While previous work has argued that one of the two approaches is preferable in general, we argue that the choice is application dependent. Furthermore, we argue that the choice comes down to whether it is desirable to be true to the model or true to the data. We use linear models to investigate this choice. After deriving an efficient method for calculating observational conditional expectation Shapley values for linear models, we investigate how correlation in simulated data impacts the convergence of observational conditional expectation Shapley values. Finally, we present two real data examples that we consider to be representative of possible use cases for feature attribution -- (1) credit risk modeling and (2) biological discovery. We show how a different choice of value function performs better in each scenario, and how possible attributions are impacted by modeling choices.


title: Vision-Based Goal-Conditioned Policies for Underwater Navigation in the Presence of Obstacles

abstract: We present Nav2Goal, a data-efficient and end-to-end learning method for goal-conditioned visual navigation. Our technique is used to train a navigation policy that enables a robot to navigate close to sparse geographic waypoints provided by a user without any prior map, all while avoiding obstacles and choosing paths that cover user-informed regions of interest. Our approach is based on recent advances in conditional imitation learning. General-purpose, safe and informative actions are demonstrated by a human expert. The learned policy is subsequently extended to be goal-conditioned by training with hindsight relabelling, guided by the robot's relative localization system, which requires no additional manual annotation. We deployed our method on an underwater vehicle in the open ocean to collect scientifically relevant data of coral reefs, which allowed our robot to operate safely and autonomously, even at very close proximity to the coral. Our field deployments have demonstrated over a kilometer of autonomous visual navigation, where the robot reaches on the order of 40 waypoints, while collecting scientifically relevant data. This is done while travelling within 0.5 m altitude from sensitive corals and exhibiting significant learned agility to overcome turbulent ocean conditions and to actively avoid collisions.


title: Transformers are RNNs: Fast Autoregressive Transformers with Linear Attention

abstract: Transformers achieve remarkable performance in several tasks but due to their quadratic complexity, with respect to the input's length, they are prohibitively slow for very long sequences. To address this limitation, we express the self-attention as a linear dot-product of kernel feature maps and make use of the associativity property of matrix products to reduce the complexity from $\mathcal{O}\left(N^2\right)$ to $\mathcal{O}\left(N\right)$, where $N$ is the sequence length. We show that this formulation permits an iterative implementation that dramatically accelerates autoregressive transformers and reveals their relationship to recurrent neural networks. Our linear transformers achieve similar performance to vanilla transformers and they are up to 4000x faster on autoregressive prediction of very long sequences.


title: Webster Sequences, Apportionment Problems, and Just-In-Time Sequencing

abstract: Given a real number $\alpha \in (0,1)$, we define the Webster sequence of density $\alpha$ to be $W_\alpha = (\lceil(n-1/2) / \alpha\rceil)_{n\in\mathbb{N}}$, where $\lceil x \rceil$ is the ceiling function. It is known that if $\alpha$ and $\beta$ are irrational with $\alpha + \beta = 1$, then $W_\alpha$ and $W_\beta$ partition $\mathbb{N}$. On the other hand, an analogous result for three-part partitions does not hold: There does not exist a partition of $\mathbb{N}$ into sequences $W_\alpha, W_\beta, W_\gamma$ with $\alpha, \beta, \gamma $ irrational. In this paper, we consider the question of how close one can come to a three-part partition of $\mathbb{N}$ into Webster sequences with prescribed densities $\alpha, \beta, \gamma $. We show that if $\alpha, \beta, \gamma $ are irrational with $\alpha + \beta +\gamma = 1$, there exists a partition of $\mathbb{N}$ into sequences of densities $\alpha, \beta, \gamma$, in which one of the sequences is a Webster sequence and the other two are "almost" Webster sequences that are obtained from Webster sequences by perturbing some elements by at most 1. We also provide two efficient algorithms to construct such partitions. The first algorithm outputs the $n$th element of each sequence in $O(1)$ time and the second algorithm gives the assignment of the $m$th positive integer to a sequence in $O(1)$ time. We show that the results are best-possible in several respects. Moreover, we describe applications of these results to apportionment and optimal scheduling problems.


title: On dynamics of nonmagnetic accretion disks

abstract: Axisymmetric accretion disks in vicinity of a central compact body are studied. In the case of non-viscous disk it is proven that all solutions for the midplane circular velocity are unstable. Hence, the pure hydrodynamic turbulence in accretion disks is possible. It is disproved the well-known arguments that an inviscid accretion disk must be sub-Keplerian. It is also demonstrated that the regular asymptotic solutions, often used in astrophysics, can lead to erroneous conclusions. It is proven that a laminar viscous disk can be approximated with a great precision by the vortex motion. Assuming that a turbulent gas tends to flow with minimal losses, we have shown that a turbulent disk tends to be Keplerian.


title: An Imitation Learning Approach for Cache Replacement

abstract: Program execution speed critically depends on increasing cache hits, as cache hits are orders of magnitude faster than misses. To increase cache hits, we focus on the problem of cache replacement: choosing which cache line to evict upon inserting a new line. This is challenging because it requires planning far ahead and currently there is no known practical solution. As a result, current replacement policies typically resort to heuristics designed for specific common access patterns, which fail on more diverse and complex access patterns. In contrast, we propose an imitation learning approach to automatically learn cache access patterns by leveraging Belady's, an oracle policy that computes the optimal eviction decision given the future cache accesses. While directly applying Belady's is infeasible since the future is unknown, we train a policy conditioned only on past accesses that accurately approximates Belady's even on diverse and complex access patterns, and call this approach Parrot. When evaluated on 13 of the most memory-intensive SPEC applications, Parrot increases cache miss rates by 20% over the current state of the art. In addition, on a large-scale web search benchmark, Parrot increases cache hit rates by 61% over a conventional LRU policy. We release a Gym environment to facilitate research in this area, as data is plentiful, and further advancements can have significant real-world impact.


title: Network broadcast analysis and control of turbulent flows

abstract: We present a network-based modal analysis technique that identifies key dynamical paths along which perturbations amplify over a time-varying base flow. This analysis is built upon the Katz centrality, which reveals the flow structures that can effectively spread perturbations over a time-evolving network of vortical interactions on the base flow. Motivated by the resolvent form of the Katz function, we take the singular value decomposition of the resulting communicability matrix, complementing the resolvent analysis for fluid flows. The right-singular vectors, referred to as the broadcast modes, give insights into the sensitive regions where introduced perturbations can be effectively spread and amplified over the entire fluid-flow network that evolves in time. We apply this analysis to a two-dimensional decaying isotropic turbulence. The broadcast mode reveals that vortex dipoles are important structures in spreading perturbations. By perturbing the flow with the principal broadcast mode, we demonstrate the utility of the insights gained from the present analysis to effectively modify the evolution of turbulent flows. The current network-inspired work presents a novel use of network analysis to guide flow control efforts, in particular for time-varying base flows.


title: The Many Faces of Robustness: A Critical Analysis of Out-of-Distribution Generalization

abstract: We introduce four new real-world distribution shift datasets consisting of changes in image style, image blurriness, geographic location, camera operation, and more. With our new datasets, we take stock of previously proposed methods for improving out-of-distribution robustness and put them to the test. We find that using larger models and artificial data augmentations can improve robustness on real-world distribution shifts, contrary to claims in prior work. We find improvements in artificial robustness benchmarks can transfer to real-world distribution shifts, contrary to claims in prior work. Motivated by our observation that data augmentations can help with real-world distribution shifts, we also introduce a new data augmentation method which advances the state-of-the-art and outperforms models pretrained with 1000 times more labeled data. Overall we find that some methods consistently help with distribution shifts in texture and local image statistics, but these methods do not help with some other distribution shifts like geographic changes. Our results show that future research must study multiple distribution shifts simultaneously, as we demonstrate that no evaluated method consistently improves robustness.


title: The Heterogeneity Hypothesis: Finding Layer-Wise Differentiated Network Architectures

abstract: In this paper, we tackle the problem of convolutional neural network design. Instead of focusing on the design of the overall architecture, we investigate a design space that is usually overlooked, i.e. adjusting the channel configurations of predefined networks. We find that this adjustment can be achieved by shrinking widened baseline networks and leads to superior performance. Based on that, we articulate the heterogeneity hypothesis: with the same training protocol, there exists a layer-wise differentiated network architecture (LW-DNA) that can outperform the original network with regular channel configurations but with a lower level of model complexity. The LW-DNA models are identified without extra computational cost or training time compared with the original network. This constraint leads to controlled experiments which direct the focus to the importance of layer-wise specific channel configurations. LW-DNA models come with advantages related to overfitting, i.e. the relative relationship between model complexity and dataset size. Experiments are conducted on various networks and datasets for image classification, visual tracking and image restoration. The resultant LW-DNA models consistently outperform the baseline models. Code is available at https://github.com/ofsoundof/Heterogeneity_Hypothesis.


title: Mutual Information in Coupled Double Quantum Dots: A Simple Analytic Model for Potential Artificial Consciousness

abstract: The integrated information theory is thought to be a key clue towards the theoretical understanding of consciousness. In this study, we propose a simple numerical model comprising a set of coupled double quantum dots, where the disconnection of the elements is represented by the removal of Coulomb interaction between the quantum dots, for the quantitative investigation of integrated information. As a measure of integrated information, we calculate the mutual information in the model system, as the Kullback-Leibler divergence between the connected and disconnected status, through the probability distribution of the electronic states from the master transition-rate equations. We reasonably demonstrate that the increase in the strength of interaction between the quantum dots leads to higher mutual information, owing to the larger divergence in the probability distributions of the electronic states. Our model setup could be a useful basic tool for numerical analyses in the field of integrated information theory.


title: Filtering of stationary Gaussian statistical experiments

abstract: This article proposes a new filtering model for stationary Gaussian Markov statistical experiments, given by diffusion-type difference stochastic equations.


title: Any Three Longest Paths In A Connected Graph Has A Common Vertex

abstract: A question was raised in 1995 at the British Combinatorial Conference: Do any three longest paths in a connected graph have a vertex in common? In this paper, it is shown that the answer to that question is yes.


title: Mathematical modelling of charge transport in graphene heterojunctions

abstract: A typical graphene heterojunction device can be divided into two classical zones, where the transport is basically diffusive, separated by a "quantum active region" (e.g., a locally gated region), where the charge carriers are scattered according to the laws of quantum mechanics.In this paper we derive a mathematical model of such a device, where the classical regions are described by drift-diffusion equations and the quantum zone is seen as an interface where suitable transmission conditions are imposed that take into account the quantum scattering process. Numerical simulations show good agreement with experimental data.


title: Long Working Distance Apodized Grating Coupler

abstract: We design a focusing grating coupler by a simultaneous apodization of the filling factor and the period. In addition to in plane focusing to the input waveguide providing a total length of less than 70{\mu}m, a further apodization of the curvature allows out of plane focusing into a fiber set 150{\mu}m away from the grating surface. The design is proposed for novel semiconductor on insulator waveguides such as GaInP. The coupling efficiency is calculated from two dimensional simulation which is about 45%.


title: Faster Digital Quantum Simulation by Symmetry Protection

abstract: Simulating the dynamics of quantum systems is an important application of quantum computers and has seen a variety of implementations on current hardware. We show that by introducing quantum gates implementing unitary transformations generated by the symmetries of the system, one can induce destructive interference between the errors from different steps of the simulation, effectively giving faster quantum simulation by symmetry protection. We derive rigorous bounds on the error of a symmetry-protected simulation algorithm and identify conditions for optimal symmetry protection. In particular, when the symmetry transformations are chosen as powers of a unitary, the error of the algorithm is approximately projected to the so-called quantum Zeno subspaces. We prove a bound on this approximation error, exponentially improving a recent result of Burgarth, Facchi, Gramegna, and Pascazio. We apply the symmetry protection technique to the simulations of the XXZ Heisenberg interactions with local disorder and the Schwinger model in quantum field theory. For both systems, the technique can reduce the simulation error by several orders of magnitude over the unprotected simulation. Finally, we provide numerical evidence suggesting that the technique can also protect simulation against other types of coherent, temporally correlated errors, such as the $1/f$ noise commonly found in solid-state experiments.


title: Primordial black holes from cusp collapse on cosmic strings

abstract: Primordial black holes (PBHs) are of fundamental interest in cosmology and astrophysics, and have received much attention as a dark matter candidate and as a potential source of gravitational waves. One possible PBH formation mechanism is the gravitational collapse of cosmic strings. Thus far, the entirety of the literature on PBH production from cosmic strings has focused on the collapse of (quasi)circular cosmic string loops, which make up only a tiny fraction of the cosmic loop population. We demonstrate here a novel PBH formation mechanism: the collapse of a small segment of cosmic string in the neighbourhood of a cusp. Using the hoop conjecture, we show that collapse is inevitable whenever a cusp appears on a macroscopically-large loop, forming a PBH whose rest mass is smaller than the mass of the loop by a factor of the dimensionless string tension squared, $(G\mu)^2$. Since cusps are generic features of cosmic string loops, and do not rely on finely-tuned loop configurations like circular collapse, this implies that cosmic strings produce PBHs in far greater numbers than has previously been recognised. The resulting PBHs are highly spinning and boosted to ultrarelativistic velocities; they populate a unique region of the BH mass-spin parameter space, and are therefore a "smoking gun" observational signature of cosmic strings. We derive new constraints on $G\mu$ from the evaporation of cusp-collapse PBHs, and update existing constraints on $G\mu$ from gravitational-wave searches.


title: Vector Leptoquarks Beyond Tree Level II: $\mathcal{O}(\alpha_s)$ Corrections and Radial Modes

abstract: We clarify open issues in relating low- and high-energy observables, at next-to-leading order accuracy, in models with a massive leptoquark embedded in a flavor non-universal $SU(4) \times SU(3) \times SU(2) \times U(1)$ gauge group. Extending previous work on this subject, we present a complete analysis of the $\mathcal{O}(\alpha_s)$ corrections to the matching conditions of semileptonic operators at the high scale. These corrections are not negligible, but they do not exceed the 10% level and are subleading compared to the $\mathcal{O}(\alpha_4)$ corrections proportional to the leading leptoquark coupling, which is expected to be much larger than the QCD coupling in the parameter space region of phenomenological interest. We further analyze the impact of radial modes, both at $\mathcal{O}(\alpha_4)$ and at $\mathcal{O}(\alpha_s)$ accuracy, highlighting their role in the renormalization of the theory.


title: Dark Matter Annihilation Can Produce a Detectable Antihelium Flux through $\bar{\Lambda}_b$ Decays

abstract: Recent observations by the Alpha Magnetic Spectrometer (AMS-02) have tentatively detected a handful of cosmic-ray antihelium events. Such events have long been considered as smoking-gun evidence for new physics, because astrophysical antihelium production is expected to be negligible. However, the dark-matter-induced antihelium flux is also expected to fall below current sensitivities, particularly in light of existing antiproton constraints. Here, we demonstrate that a previously neglected standard model process -- the production of antihelium through the displaced-vertex decay of $\bar{\Lambda}_b$-baryons -- can significantly boost the dark matter induced antihelium flux. This process can triple the standard prompt-production of antihelium, and more importantly, entirely dominate the production of the high-energy antihelium nuclei reported by AMS-02.


title: Implications of the XENON1T Excess on the Dark Matter Interpretation

abstract: The dark matter interpretation for a recent observation of excessive electron recoil events at the XENON1T detector seems challenging because its velocity is not large enough to give rise to recoiling electrons of $\mathcal{O}({\rm keV})$. Fast-moving or boosted dark matter scenarios are receiving attention as a remedy for this issue, rendering the dark matter interpretation a possibility to explain the anomaly. We investigate various scenarios where such dark matter of spin 0 and 1/2 interacts with electrons via an exchange of vector, pseudo-scalar, or scalar mediators. We find parameter values not only to reproduce the excess but to be consistent with existing bounds. Our study suggests that the scales of mass and coupling parameters preferred by the excess can be mostly affected by the type of mediator, and that significantly boosted dark matter can explain the excess depending on the mediator type and its mass choice. The method proposed in this work is general, and hence readily applicable to the interpretation of observed data in the dark matter direct detection experiment.


title: Strong Coupling Limit of the Holstein-Hubbard Model

abstract: We analyze the quantum phase diagram of the Holstein-Hubbard model using an asymptotically exact strong-coupling expansion. We find all sorts of interesting phases including a pair-density wave (PDW), a charge 4e (and even a charge 6e) superconductor, regimes of phase separation, and a variety of distinct charge-density-wave (CDW), spin-density-wave (SDW) and superconducting regimes. We chart the crossovers that occur as a function of the degree of retardation, i.e. the ratio of characteristic phonon frequency to the strength of interactions.


title: Ellipticity of photon emission from strongly magnetized hot QCD plasma

abstract: By making use of an explicit representation for the imaginary part of the photon polarization tensor in terms of transitions between the Landau levels of light quarks, we study the angular dependence of direct photon emission from a strongly magnetized quark-gluon plasma. Because of the magnetic field, the leading order photon rate comes from the three processes of the zeroth order in the coupling constant $\alpha_s$: (i) the quark splitting ($q\rightarrow q+\gamma $), (ii) the antiquark splitting ($\bar{q} \rightarrow \bar{q}+\gamma $), and (iii) the quark-antiquark annihilation ($q + \bar{q}\rightarrow \gamma$). In a wide range of moderately high temperatures, $T\gtrsim m_\pi$, and strong magnetic fields, $|eB|\gtrsim m_\pi^2$, the direct photon production is dominated by the two splitting processes. We show that the Landau-level quantization of quark states plays an important role in the energy and angular dependence of the photon emission. Among other things, it leads to a nontrivial momentum dependence of the photon ellipticity coefficient $v_2$, which takes negative values at small transverse momenta and positive values at large transverse momenta. The crossover between the two regimes occurs around $k_T\simeq \sqrt{|eB|}$. In application to heavy-ion collisions, this suggests that a large value of $v_2$ for the direct photons could be explained in part by the magnetic field in the quark-gluon plasma.


title: Towards the classification of rank-$r$ $\mathcal{N}=2$ SCFTs. Part I: twisted partition function and central charge formulae

abstract: We derive explicit formulae to compute the $a$ and $c$ central charges of four dimensional $\mathcal{N}=2$ superconformal field theories (SCFTs) directly from Coulomb branch related quantities. The formulae apply at arbitrary rank. We also discover general properties of the low-energy limit behavior of the flavor symmetry of $\mathcal{N}=2$ SCFTs which culminate with our $\mathcal{N}=2$ UV-IR simple flavor condition. This is done by determining precisely the relation between the integrand of the partition function of the topologically twisted version of the 4d $\mathcal{N}=2$ SCFTs and the singular locus of their Coulomb branches. The techniques developed here are extensively applied to many rank-2 SCFTs, including new ones, in a companion paper.


title: Unification of parton and coupled-wire approaches to quantum magnetism in two dimensions

abstract: The fractionalization of microscopic degrees of freedom is a remarkable manifestation of strong interactions in quantum many-body systems. Analytical studies of this phenomenon are primarily based on two distinct frameworks: field theories of partons and emergent gauge fields, or coupled arrays of one-dimensional quantum wires. We unify these approaches for two-dimensional spin systems. Via exact manipulations, we demonstrate how parton gauge theories arise in microscopic wire arrays and explicitly relate spin operators to emergent quasiparticles and gauge-field monopoles. This correspondence allows us to compute physical correlation functions within both formulations and leads to a straightforward algorithm for constructing parent Hamiltonians for a wide range of exotic phases. We exemplify this technique for several chiral and non-chiral quantum spin liquids.


title: Dust Populations in the Iconic Vega Planetary System Resolved by ALMA

abstract: The Vega planetary system hosts the archetype of extrasolar Kuiper belts, and is rich in dust from the sub-au region out to 100's of au, suggesting intense dynamical activity. We present ALMA mm observations that detect and resolve the outer dust belt from the star for the first time. The interferometric visibilities show that the belt can be fit by a Gaussian model or by power-law models with a steep inner edge (at 60-80 au). The belt is very broad, extending out to at least 150-200 au. We strongly detect the star and set a stringent upper limit to warm dust emission previously detected in the infrared. We discuss three scenarios that could explain the architecture of Vega's planetary system, including the new {ALMA} constraints: no outer planets, a chain of low-mass planets, and a single giant planet. The planet-less scenario is only feasible if the outer belt was born with the observed sharp inner edge. If instead the inner edge is currently being truncated by a planet, then the planet must be $\gtrsim$6 M$_{\oplus}$ and at $\lesssim71$ au to have cleared its chaotic zone within the system age. In the planet chain scenario, outward planet migration and inward scattering of planetesimals could produce the hot and warm dust observed in the inner regions of the system. In the single giant planet scenario, an asteroid belt could be responsible for the warm dust, and mean motion resonances with the planet could put asteroids on star-grazing orbits, producing the hot dust.


title: Data-Driven Stellar Models

abstract: We develop a data-driven model to map stellar parameters (effective temperature, surface gravity and metallicity) accurately and precisely to broad-band stellar photometry. This model must, and does, simultaneously constrain the passband-specific dust reddening vector in the Milky Way. The model uses a neural network to learn the (de-reddened) absolute magnitude in one band and colors across many bands, given stellar parameters from spectroscopic surveys and parallax constraints from Gaia. To demonstrate the effectiveness of this approach, we train our model on a dataset with spectroscopic parameters from LAMOST, APOGEE and GALAH, Gaia parallaxes, and optical and near-infrared photometry from Gaia, Pan-STARRS~1, 2MASS and WISE. Testing the model on these datasets leads to an excellent fit and a precise - and by construction accurate - prediction of the color-magnitude diagrams in many bands. This flexible approach rigorously links spectroscopic and photometric surveys, and also results in an improved, stellar-temperature-dependent reddening vector. As such, it provides a simple and accurate method for predicting photometry in stellar evolutionary models. Our model will form a basis to infer stellar properties, distances and dust extinction from photometric data, which should be of great use in 3D mapping of the Milky Way. Our trained model may be obtained at https://doi.org/10.5281/zenodo.3902382.


title: The Ophiuchus DIsc Survey Employing ALMA (ODISEA). II. The effect of stellar multiplicity on disc properties

abstract: We present Adaptive Optics (AO) near infrared (NIR) observations using VLT/NACO and Keck/NIRC2 of ODISEA targets. ODISEA is an ALMA survey of the entire population of circumstellar discs in the Ophiuchus molecular cloud. From the whole sample of ODISEA we select all the discs that are not already observed in the NIR with AO and that are observable with NACO or NIRC2. The NIR-ODISEA survey consists of 147 stars observed in NIR AO imaging for the first time, as well as revisiting almost all the binary systems of Ophiuchus present in the literature (20 out of 21). In total, we detect 20 new binary systems and one triple system. For each of them we calculate the projected separation and position angle of the companion, as well as their NIR and millimeter flux ratios. From the NIR contrast we derived the masses of the secondaries, finding that 9 of them are in the sub-stellar regime (30-50 \MJup). Discs in multiple systems reach a maximum total dust mass of $\sim$ 50 M$_{\oplus}$, while discs in single stars can reach a dust mass of 200 M$_{\oplus}$. Discs with masses above 10 M$_{\oplus}$ are found only around binaries with projected separations larger than $\sim$ 110 au. The maximum disc size is also larger around single star than binaries. However, since most discs in Ophiuchus are very small and low-mass, the effect of visual binaries is relatively weak in the general disc population.


title: The Fermionic Universal One-Loop Effective Action

abstract: Recent development of path integral matching techniques based on the covariant derivative expansion has made manifest a universal structure of one-loop effective Lagrangians. The universal terms can be computed once and for all to serve as a reference for one-loop matching calculations and to ease their automation. Here we present the fermionic universal one-loop effective action (UOLEA), resulting from integrating out heavy fermions with scalar, pseudo-scalar, vector and axial-vector couplings. We also clarify the relation of the new terms computed here to terms previously computed in the literature and those that remain to complete the UOLEA. Our results can be readily used to efficiently obtain analytical expressions for effective operators arising from heavy fermion loops.


title: Scaling K2. III. Comparable Planet Occurrence in the FGK Samples of Campaign 5 and Kepler

abstract: Using our K2 Campaign 5 fully automated planet detection data set (43 planets), which has corresponding measures of completeness and reliability, we infer an underlying planet population model for the FGK dwarfs sample (9,257 stars). Implementing a broken power-law for both the period and radius distribution, we find an overall planet occurrence of $1.00^{+1.07}_{-0.51}$ planets per star within a period range of 0.5-38 days. Making similar cuts and running a comparable analysis on the Kepler sample (2,318 planets; 94,222 stars), we find an overall occurrence of $1.10\pm0.05$ planets per star. Since the Campaign 5 field is nearly 120 angular degrees away from the Kepler field, this occurrence similarity offers evidence that the Kepler sample may provide a good baseline for Galactic inferences. Furthermore, the Kepler stellar sample is metal-rich compared to the K2 Campaign 5 sample, thus a finding of occurrence parity may reduce the role of metallicity in planet formation. However, a weak ($1.5\sigma$) difference, in agreement with metal-driven formation, is found when assuming the Kepler model power-laws for the K2 Campaign 5 sample and optimizing only the planet occurrence factor. This weak trend indicates further investigation of metallicity dependent occurrence is warranted once a larger sample of uniformly vetted K2 planet candidates are made available.


title: High-precision distance measurements with classical pulsating stars

abstract: Classical Cepheid and RR Lyrae variables are radially pulsating stars that trace young and old-age stellar populations, respectively. These classical pulsating stars are the most sensitive probes for the precision stellar astrophysics and the extragalactic distance measurements. Despite their extensive use as standard candles thanks to their well-defined Period-Luminosity relations, distance measurements based on these objects suffer from their absolute primary calibrations, metallicity effects, and other systematic uncertainties. Here, I present a review of classical Cepheid, RR Lyrae, and Type II Cepheid variables starting with a historical introduction and describing their basic evolutionary and pulsational properties. I will focus on recent theoretical and observational efforts to establish absolute scale for these standard candles at multiple wavelengths. The application of these classical pulsating stars to high-precision cosmic distance scale will be discussed along with observational systematics. I will summarize with an outlook for further improvements in our understanding of these classical pulsators in the upcoming era of extremely large telescopes.


title: Plasma lensing in comparison to gravitational lensing -- Formalism and degeneracies

abstract: Gravitational and plasma lensing share the same mathematical formalism in the limit of geometrical optics. Both phenomena can be effectively described by a projected, two-dimensional deflection potential whose gradient causes an instantaneous light deflection in a single, thin lens plane. We highlight the differences in the time-delay and lensing equations that occur because plasma lensing is caused by a potential directly proportional to the deflecting electron number density and gravitational lensing is caused by a potential related to the deflecting mass density by a Poisson equation. Since we treat plasma and gravitational lensing as thin-screen effective theories, their degeneracies are both caused by the unknown distribution of deflecting objects. Deriving the formalism-intrinsic degeneracies for plasma lensing, we find that they are analogous to those occurring in gravitational lensing. To break the degeneracies, galaxies and galaxy-cluster scale strong gravitational lenses must rely on additional assumptions or complementary observations. Physically realistic assumptions to arrive at self-consistent lens and source reconstructions can be provided by simulations and analytical effective theories. In plasma lensing, a deeper understanding of the deflecting electron density distributions is still under development, so that a model-based comprehensive lens reconstruction is not yet possible. However, we show that transient lenses and multi-wavelength observations help to break the arising degeneracies. We conclude that the development of an observation-based inference of local lens properties seems currently the best way to further probe the morphologies of plasma electron densities. Due to the simpler evidence-based breaking of the lensing degeneracies, we expect to obtain tighter constraints on the local plasma electron densities than on the gravitationally deflecting masses.


title: CoMaLit -- VI. Intrinsic scatter in stacked relations. The weak lensing AMICO galaxy clusters in KiDS-DR3

abstract: Unbiased and precise mass calibration of galaxy clusters is crucial to fully exploit galaxy clusters as cosmological probes. Stacking of weak lensing signal allows us to measure observable-mass relations down to less massive halos halos without extrapolation. We propose a Bayesian inference method to constrain the intrinsic scatter of the mass proxy in stacked analyses. The scatter of the stacked data is rescaled with respect to the individual scatter based on the number of binned clusters. We apply this method to the galaxy clusters detected with the AMICO (Adaptive Matched Identifier of Clustered Objects) algorithm in the third data release of the Kilo-Degree Survey. The results confirm the optical richness as a low scatter mass proxy. Based on the optical richness and the calibrated weak lensing mass-richness relation, mass of individual objects down to ~10^13 solar masses can be estimated with a precision of ~20 per cent.


title: The Standard Model, The Exceptional Jordan Algebra, and Triality

abstract: Jordan, Wigner and von Neumann classified the possible algebras of quantum mechanical observables, and found they fell into 4 "ordinary" families, plus one remarkable outlier: the exceptional Jordan algebra. We point out an intriguing relationship between the complexification of this algebra and the standard model of particle physics, its minimal left-right-symmetric $SU(3)\times SU(2)_{L}\times SU(2)_{R}\times U(1)$ extension, and $Spin(10)$ unification. This suggests a geometric interpretation, where a single generation of standard model fermions is described by the tangent space $(\mathbb{C}\otimes\mathbb{O})^{2}$ of the complex octonionic projective plane, and the existence of three generations is related to $SO(8)$ triality.


title: HFQPOs and discoseismic mode excitation in eccentric, relativistic discs. II. Magnetohydrodynamic simulations

abstract: Trapped inertial oscillations (r-modes) provide a promising explanation for high-frequency quasi-periodic oscillations (HFQPOs) observed in the emission from black hole X-ray binary systems. An eccentricity (or warp) can excite r-modes to large amplitudes, but concurrently the oscillations are likely damped by magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability (MRI). We force eccentricity in global, unstratified, zero-net flux MHD simulations of relativistic accretion discs, and find that a sufficiently strong disc distortion generates trapped inertial waves despite this damping. In our simulations, eccentricities above ~ 0.03 in the inner disc excite trapped waves. In addition to the competition between r-mode damping and driving, we observe that larger amplitude eccentric structures modify and in some cases suppress MRI turbulence. Given the variety of distortions (warps as well as eccentricities) capable of amplifying r-modes, the robustness of trapped inertial wave excitation in the face of MRI turbulence in our simulations provides support for a discoseismic explanation for HFQPOs.


title: XENON1T solar axion and the Higgs boson emerging from the dark

abstract: In a recent letter we proposed a new non-thermal mechanism of Dark Matter production based on vacuum misalignment, where both the Higgs boson and a very light pseudo-scalar $\eta$ emerge from the Dark sector. In this letter, we identify the parameter space in a composite scenario where the light pseudo-scalar can be produced in the sun and explain the XENON1T excess in electron recoil data. The model's Dark Matter candidate has a mass around $50$ TeV and out of range for Direct Detection. Testable predictions include Gravitational waves at frequencies in the Hz range from a cosmological phase transition, an exotic decay $Z \to \gamma + \mbox{inv.}$ with rates $4 \div 16 \cdot 10^{-12}$ testable at a future Tera-Z collider, and an enhancement by $17\div 40$ % of the branching ratio $K_L \to \pi^0 + \mbox{inv.}$, not enough to explain the KOTO anomaly. All these predictions may be confirmed by future experiments.


title: No Redshift Evolution in the Broad Line Region Metallicity up to z=7.54: Deep NIR Spectroscopy of ULAS J1342+0928

abstract: We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.


title: Reinforcement Learning for Digital Quantum Simulation

abstract: Digital quantum simulation is a promising application for quantum computers. Their free programmability provides the potential to simulate the unitary evolution of any many-body Hamiltonian with bounded spectrum by discretizing the time evolution operator through a sequence of elementary quantum gates, typically achieved using Trotterization. A fundamental challenge in this context originates from experimental imperfections for the involved quantum gates, which critically limits the number of attainable gates within a reasonable accuracy and therefore the achievable system sizes and simulation times. In this work, we introduce a reinforcement learning algorithm to systematically build optimized quantum circuits for digital quantum simulation upon imposing a strong constraint on the number of allowed quantum gates. With this we consistently obtain quantum circuits that reproduce physical observables with as little as three entangling gates for long times and large system sizes. As concrete examples we apply our formalism to a long range Ising chain and the lattice Schwinger model. Our method makes larger scale digital quantum simulation possible within the scope of current experimental technology.


title: Classical algebraic structures in string theory effective actions

abstract: We study generic properties of string theory effective actions obtained by classically integrating out massive excitations from string field theories based on cyclic homotopy algebras of $A_\infty$ or $L_\infty$ type. We construct observables in the UV theory and we discuss their fate after integration-out. Furthermore, we discuss how to compose two subsequent integrations of degrees of freedom (horizontal composition) and how to integrate out degrees of freedom after deforming the UV theory with a new consistent interaction (vertical decomposition). We then apply our general results to the open bosonic string using Witten's open string field theory. There we show how the horizontal composition can be used to systematically integrate out the Nakanishi-Lautrup field from the set of massless excitations, ending with a non-abelian $A_\infty$-gauge theory for just the open string gluon. Moreover we show how the vertical decomposition can be used to construct effective open-closed couplings by deforming Witten OSFT with a tadpole given by the Ellwood invariant. Also, we discuss how the effective theory controls the possibility of removing the tadpole in the microscopic theory, giving a new framework for studying D-branes deformations induced by changes in the closed string background.


title: Entanglement Entropy of Fermions from Wigner Functions: Excited States and Open Quantum Systems

abstract: We formulate a new ``Wigner characteristics'' based method to calculate entanglement entropies of subsystems of Fermions using Keldysh field theory. This bypasses the requirements of working with complicated manifolds for calculating R\'{e}nyi entropies for many body systems. We provide an exact analytic formula for R\'{e}nyi and von-Neumann entanglement entropies of non-interacting open quantum systems, which are initialised in arbitrary Fock states. We use this formalism to look at entanglement entropies of momentum Fock states of one-dimensional Fermions. We show that the entanglement entropy of a Fock state can scale either logarithmically or linearly with subsystem size, depending on whether the number of discontinuities in the momentum distribution is smaller or larger than the subsystem size. This classification of states in terms number of blocks of occupied momenta allows us to analytically estimate the number of critical and non-critical Fock states for a particular subsystem size. We also use this formalism to describe entanglement dynamics of an open quantum system starting with a single domain wall at the center of the system. Using entanglement entropy and mutual information, we understand the dynamics in terms of coherent motion of the domain wall wavefronts, creation and annihilation of domain walls and incoherent exchange of particles with the bath.


title: Straight Lightning as a Signature of Macroscopic Dark Matter

abstract: Macroscopic dark matter (macros) is a broad class of alternative candidates to particle dark matter. These candidates would transfer energy to matter primarily through elastic scattering. A sufficiently large macro passing through the atmosphere would produce a straight channel of ionized plasma. If the cross-section of the macro is $\sigma_x \gtrapprox 6 \times 10^{-9} cm^2$, then under atmospheric conditions conducive to lightning (eg. a thunderstorm) the plasma channel would be sufficient to seed a lightning strike with a single leader. This is entirely unlike ordinary bolt lightning in which a long sequence of hundreds or thousands of few-meter-long leaders are strung together. This macro-induced lightning would be extremely straight, and thus highly distinctive. Neither wind shear nor magnetohydrodynamic instabilities would markedly spoil its straightness. The only photographically documented case of a straight lightning bolt is probably not straight enough to have been macro-induced. We estimate the region of macro parameter space that could be probed by a search for straight lightning from the number of thunderstorms happening on Earth at any time. We also estimate the parameter space that can be probed by carefully monitoring Jupiter, e.g. using the Hubble Space Telescope. All code and data is available at https://github.com/cwru-pat/macro_lightning.


title: $\textsf{Xsec}$: the cross-section evaluation code

abstract: The evaluation of higher-order cross-sections is an important component in the search for new physics, both at hadron colliders and elsewhere. For most new physics processes of interest, total cross-sections are known at next-to-leading order (NLO) in the strong coupling $\alpha_s$, and often beyond, via either higher-order terms at fixed powers of $\alpha_s$, or multi-emission resummation. However, the computation time for such higher-order cross-sections is prohibitively expensive, and precludes efficient evaluation in parameter-space scans beyond two dimensions. Here we describe the software tool $\textsf{xsec}$, which allows for fast evaluation of cross-sections based on the use of machine-learning regression, using distributed Gaussian processes trained on a pre-generated sample of parameter points. This first version of the code provides all NLO Minimal Supersymmetric Standard Model strong-production cross-sections at the LHC, for individual flavour final states, evaluated in a fraction of a second. Moreover, it calculates regression errors, as well as estimates of errors from higher-order contributions, from uncertainties in the parton distribution functions, and from the value of $\alpha_s$. While we focus on a specific phenomenological model of supersymmetry, the method readily generalises to any process where it is possible to generate a sufficient training sample.


title: Octonions, trace dynamics and non-commutative geometry: a case for unification in spontaneous quantum gravity

abstract: We have recently proposed a new matrix dynamics at the Planck scale, building on the theory of trace dynamics. This is a Lagrangian dynamics in which the matrix degrees of freedom are made from Grassmann numbers, and the Lagrangian is trace of a matrix polynomial. Matrices made from even grade elements of the Grassmann algebra are called bosonic, and those made from odd grade elements are called fermionic: together they describe an `aikyon'. In the present article we provide a basic definition of spin angular momentum in this matrix dynamics, and introduce a bosonic (fermionic) configuration variable conjugate to the spin of a boson (fermion). We then show that at energies below Planck scale, where the matrix dynamics reduces to quantum theory, fermions have half-integer spin (in multiples of Planck's constant), and bosons have integral spin. We also show that this definition of spin agrees with the conventional understanding of spin in relativistic quantum mechanics. Consequently, we obtain an elementary proof for the spin-statistics connection. We then motivate why an octonionic space is the natural space in which an aikyon evolves. The group of automorphisms in this space is the exceptional Lie group $G_2$ which has fourteen generators [could they stand for the twelve vector bosons and two degrees of freedom of the graviton? ]. The aikyon also resembles a closed string, and it has been suggested in the literature that 10-D string theory can be represented as a 2-D string in the 8-D octonionic space. From the work of Cohl Furey and others it is known that the Dixon algebra made from the four division algebras [real numbers, complex numbers, quaternions and octonions] can possibly describe the symmetries of the standard model. In the present paper we outline how in our work the Dixon algebra arises naturally, and could lead to a unification of gravity with the standard model.