title: Deep learning velocity signals allows to quantify turbulence intensity

abstract: Turbulence, the ubiquitous and chaotic state of fluid motions, is characterized by strong and statistically non-trivial fluctuations of the velocity field, over a wide range of length- and time-scales, and it can be quantitatively described only in terms of statistical averages. Strong non-stationarities hinder the possibility to achieve statistical convergence, making it impossible to define the turbulence intensity and, in particular, its basic dimensionless estimator, the Reynolds number. Here we show that by employing Deep Neural Networks (DNN) we can accurately estimate the Reynolds number within $15\%$ accuracy, from a statistical sample as small as two large-scale eddy-turnover times. In contrast, physics-based statistical estimators are limited by the rate of convergence of the central limit theorem, and provide, for the same statistical sample, an error at least $100$ times larger. Our findings open up new perspectives in the possibility to quantitatively define and, therefore, study highly non-stationary turbulent flows as ordinarily found in nature as well as in industrial processes.

categories: physics.flu-dyn cond-mat.stat-mech cs.AI cs.LG


title: Atomic Services: sustainable ecosystem of smart city services through pan-European collaboration

abstract: In a world with an ever increasing urbanization, governance is investigating innovative solutions to sustain the society evolution. Internet-of-Things promises huge benefits for cities and the proliferation of smart city deployments demonstrates the common acceptance of IoT as basis for many solutions. The city pilots developments occurred in parallel and with different designs thus creating fragmentation of IoT. The European project SynchroniCity aims to synchronize 8 smart cities to establish a shared environment fostering a self-sustained business growth. In this article we present the collaborative methodology and shared efforts spent towards the creation of a common ecosystem for the development of smart city services. Our design evolves around the concept of "atomic services" that implements a single functional block to be composed for full-fledged smart city services. This creates opportunities for diverse stakeholders to participate to a global smart cities market. The methodology and outcome of our efforts will be followed by 10 new cities globally, thus expanding the market range for IoT stakeholders


title: Improved Concentration Bounds for Gaussian Quadratic Forms

abstract: For a wide class of monotonic functions $f$, we develop a Chernoff-style concentration inequality for quadratic forms $Q_f \sim \sum\limits_{i=1}^n f(\eta_i) (Z_i + \delta_i)^2$, where $Z_i \sim N(0,1)$. The inequality is expressed in terms of traces that are rapid to compute, making it useful for bounding p-values in high-dimensional screening applications. The bounds we obtain are significantly tighter than those that have been previously developed, which we illustrate with numerical examples.


title: A Relativized Alon Second Eigenvalue Conjecture for Regular Base Graphs IV: An Improved Sidestepping Theorem

abstract: This is the fourth in a series of articles devoted to showing that a typical covering map of large degree to a fixed, regular graph has its new adjacency eigenvalues within the bound conjectured by Alon for random regular graphs. In this paper we prove a {\em Sidestepping Theorem} that is more general and easier to use than earlier theorems of this kind. Such theorems concerns a family probability spaces $\{{\mathcal{M}}_n\}$ of $n\times n$ matrices, where $n$ varies over some infinite set, $N$, of natural numbers. Many trace methods use simple "Markov bounds" to bound the expected spectral radius of elements of ${\mathcal{M}}_n$: this consists of choosing one value, $k=k(n)$, for each $n\in N$, and proving expected spectral radius bounds based on the expected value of the trace of the $k=k(n)$-power of elements of ${\mathcal{M}}_n$. {\em Sidestepping} refers to bypassing such simple Markov bounds, obtaining improved results using a number of values of $k$ for each fixed $n\in N$. In more detail, if the $M\in {\mathcal{M}}_n$ expected value of ${\rm Trace}(M^k)$ has an asymptotic expansion in powers of $1/n$, whose coefficients are "well behaved" functions of $k$, then one can get improved bounds on the spectral radius of elements of ${\mathcal{M}}_n$ that hold with high probability. Such asymptotic expansions are shown to exist in the third article in this series for the families of matrices that interest us; in the fifth and sixth article in this series we will apply the Sidestepping Theorem in this article to prove the main results in this series of articles. This article is independent of all other articles in this series; it can be viewed as a theorem purely in probability theory, concerning random matrices or, equivalently, the $n$ random variables that are the eigenvalues of the elements of ${\mathcal{M}}_n$.


title: Momentum Contrast for Unsupervised Visual Representation Learning

abstract: We present Momentum Contrast (MoCo) for unsupervised visual representation learning. From a perspective on contrastive learning as dictionary look-up, we build a dynamic dictionary with a queue and a moving-averaged encoder. This enables building a large and consistent dictionary on-the-fly that facilitates contrastive unsupervised learning. MoCo provides competitive results under the common linear protocol on ImageNet classification. More importantly, the representations learned by MoCo transfer well to downstream tasks. MoCo can outperform its supervised pre-training counterpart in 7 detection/segmentation tasks on PASCAL VOC, COCO, and other datasets, sometimes surpassing it by large margins. This suggests that the gap between unsupervised and supervised representation learning has been largely closed in many vision tasks.


title: Integer Programming Approaches to Balanced Connected $k$-Partition

abstract: We address the problem of partitioning a vertex-weighted connected graph into $k$ connected subgraphs that have similar weights, for a fixed integer $k\geq 2$. This problem, known as the \emph{balanced connected $k$-partition problem} ($BCP_k$), is defined as follows. Given a connected graph $G$ with nonnegative weights on the vertices, find a partition $\{V_i\}_{i=1}^k$ of $V(G)$ such that each class $V_i$ induces a connected subgraph of $G$, and the weight of a class with the minimum weight is as large as possible. It is known that $BCP_k$ is $NP$-hard even on bipartite graphs and on interval graphs. It has been largely investigated under different approaches and perspectives. On the practical side, $BCP_k$ is used to model many applications arising in police patrolling, image processing, cluster analysis, operating systems and robotics. We propose three integer linear programming formulations for the balanced connected $k$-partition problem. The first one contains only binary variables and a potentially large number of constraints that are separable in polynomial time. Some polyhedral results on this formulation, when all vertices have unit weight, are also presented. The other formulations are based on flows and have a polynomial number of constraints and variables. Preliminary computational experiments have shown that the proposed formulations outperform the other formulations presented in the literature.


title: Symmetry and Singularity Properties of Steen-Ermakov-Milne-Pinney Equations

abstract: We examine the general element of the class of ordinary differential equations, $yy^{(n+1)}+\alpha y'y^{(n)}=0$, for its Lie point symmetries. We observe that the algebraic properties of this class of equations display an attractive set of patterns, the general member of the class can have three type of Algebra, $(n+1)A_1 \oplus_s\{A_1 \oplus sl(2,R)\}$, $A_1 \oplus sl(2,R)$ or $A_2 \oplus A_1$, for different values of $\alpha$. We look at the singularity properties of these equations for various values of $\alpha$.


title: Recombination: A family of Markov chains for redistricting

abstract: Redistricting is the problem of partitioning a set of geographical units into a fixed number of districts, subject to a list of often-vague rules and priorities. In recent years, the use of randomized methods to sample from the vast space of districting plans has been gaining traction in courts of law for identifying partisan gerrymanders, and it is now emerging as a possible analytical tool for legislatures and independent commissions. In this paper, we set up redistricting as a graph partition problem and introduce a new family of Markov chains called Recombination (or ReCom) on the space of graph partitions. The main point of comparison will be the commonly used Flip walk, which randomly changes the assignment label of a single node at a time. We present evidence that ReCom mixes efficiently, especially in contrast to the slow-mixing Flip, and provide experiments that demonstrate its qualitative behavior. We demonstrate the advantages of ReCom on real-world data and explain both the challenges of the Markov chain approach and the analytical tools that it enables. We close with a short case study involving the Virginia House of Delegates.


title: Cyber Risk at the Edge: Current and future trends on Cyber Risk Analytics and Artificial Intelligence in the Industrial Internet of Things and Industry 4.0 Supply Chains

abstract: Digital technologies have changed the way supply chain operations are structured. In this article, we conduct systematic syntheses of literature on the impact of new technologies on supply chains and the related cyber risks. A taxonomic/cladistic approach is used for the evaluations of progress in the area of supply chain integration in the Industrial Internet of Things and Industry 4.0, with a specific focus on the mitigation of cyber risks. An analytical framework is presented, based on a critical assessment with respect to issues related to new types of cyber risk and the integration of supply chains with new technologies. This paper identifies a dynamic and self-adapting supply chain system supported with Artificial Intelligence and Machine Learning (AI/ML) and real-time intelligence for predictive cyber risk analytics. The system is integrated into a cognition engine that enables predictive cyber risk analytics with real-time intelligence from IoT networks at the edge. This enhances capacities and assist in the creation of a comprehensive understanding of the opportunities and threats that arise when edge computing nodes are deployed, and when AI/ML technologies are migrated to the periphery of IoT networks.


title: Artificial Intelligence Strategies for National Security and Safety Standards

abstract: Recent advances in artificial intelligence (AI) have lead to an explosion of multimedia applications (e.g., computer vision (CV) and natural language processing (NLP)) for different domains such as commercial, industrial, and intelligence. In particular, the use of AI applications in a national security environment is often problematic because the opaque nature of the systems leads to an inability for a human to understand how the results came about. A reliance on 'black boxes' to generate predictions and inform decisions is potentially disastrous. This paper explores how the application of standards during each stage of the development of an AI system deployed and used in a national security environment would help enable trust. Specifically, we focus on the standards outlined in Intelligence Community Directive 203 (Analytic Standards) to subject machine outputs to the same rigorous standards as analysis performed by humans.


title: Balanced Policy Evaluation and Learning for Right Censored Data

abstract: Individualized treatment rules can lead to better health outcomes when patients have heterogeneous responses to treatment. Very few individualized treatment rule estimation methods are compatible with a multi-treatment observational study with right censored survival outcomes. In this paper we extend policy evaluation methods to the right censored data setting. Existing approaches either make restrictive assumptions about the structure of the data, or use inverse weighting methods that increase the variance of the estimator resulting in decreased performance. We propose a method which uses balanced policy evaluation combined with an imputation approach to remove right censoring. We show that the proposed imputation approach is compatible with a large number of existing survival models and can be used to extend any individualized treatment rule estimation method to the right censored data setting. We establish the rate at which the imputed values converge to the conditional expected survival times, as well as consistency guarantees and regret bounds for the combined balanced policy with imputation approach. In simulation studies, we demonstrate the improved performance of our approach compared to existing methods. We also apply our method to data from the University of North Carolina Center for AIDS Research HIV Clinical Cohort.


title: Entanglement of Propagating Optical Modes via a Mechanical Interface

abstract: Many applications of quantum information processing (QIP) require distribution of quantum states in networks, both within and between distant nodes. Optical quantum states are uniquely suited for this purpose, as they propagate with ultralow attenuation and are resilient to ubiquitous thermal noise. Mechanical systems are then envisioned as versatile interfaces between photons and a variety of solid-state QIP platforms. Here, we demonstrate a key step towards this vision, and generate entanglement between two propagating optical modes, by coupling them to the same, cryogenic mechanical system. The entanglement persists at room temperature, where we verify the inseparability of the bipartite state and fully characterize its logarithmic negativity by homodyne tomography. We detect, without any corrections, correlations corresponding to a logarithmic negativity of $E_\mathrm{N}=0.35$. Combined with quantum interfaces between mechanical systems and solid-state qubit processors already available or under development, this paves the way for mechanical systems enabling long-distance quantum information networking over optical fiber networks.


title: Update on splitting of lifetimes of $c$ and $b$ hyperons within the heavy quark expansion and decays $\Xi_Q \to \Lambda_Q \pi$

abstract: It is argued that the results of recent precision measurements by LHCb of the lifetimes of charmed and bottom hyperons are very well consistent with the description within heavy quark expansion and allow to accurately determine matrix elements of light-flavor nonsinglet four-quark operators over the hyperons and to accurately reproduce the difference of lifemes of $\Lambda_b$ and $\Xi_b^-$. When combined with the recent LHCb results on the decay $\Xi_b^- \to \Lambda_b \pi^-$ this leads to prediction of a lower bound on the rate of the decays $\Xi_c \to \Lambda_c \pi$.


title: Reporting on Decision-Making Algorithms and some Related Ethical Questions

abstract: Companies report on their financial performance for decades. More recently they have also started to report on their environmental impact and their social responsibility. The latest trend is now to deliver one single integrated report where all stakeholders of the company can easily connect all facets of the business with their impact considered in a broad sense. The main purpose of this integrated approach is to avoid delivering data related to disconnected silos, which consequently makes it very difficult to globally assess the overall performance of an entity or a business line. In this paper, we focus on how companies report on risks and ethical issues related to the increasing use of Artificial Intelligence (AI). We explain some of these risks and potential issues. Next, we identify some recent initiatives by various stakeholders to define a global ethical framework for AI. Finally, we illustrate with four cases that companies are very shy to report on these facets of AI.


title: Antithetic integral feedback for the robust control of monostable and oscillatory biomolecular circuits

abstract: Biomolecular feedback systems are now a central application area of interest within control theory. While classical control techniques provide invaluable insight into the function and design of both natural and synthetic biomolecular systems, there are certain aspects of biological control that have proven difficult to analyze with traditional methods. To this end, we describe here how the recently developed tools of dominance analysis can be used to gain insight into the nonlinear behavior of the antithetic integral feedback circuit, a recently discovered control architecture which implements integral control of arbitrary biomolecular processes using a simple feedback mechanism. We show that dominance theory can predict both monostability and periodic oscillations in the circuit, depending on the corresponding parameters and architecture. We then use the theory to characterize the robustness of the asymptotic behavior of the circuit in a nonlinear setting.


title: The phonetic bases of vocal expressed emotion: natural versus acted

abstract: Can vocal emotions be emulated? This question has been a recurrent concern of the speech community, and has also been vigorously investigated. It has been fueled further by its link to the issue of validity of acted emotion databases. Much of the speech and vocal emotion research has relied on acted emotion databases as valid proxies for studying natural emotions. To create models that generalize to natural settings, it is crucial to work with valid prototypes -- ones that can be assumed to reliably represent natural emotions. More concretely, it is important to study emulated emotions against natural emotions in terms of their physiological, and psychological concomitants. In this paper, we present an on-scale systematic study of the differences between natural and acted vocal emotions. We use a self-attention based emotion classification model to understand the phonetic bases of emotions by discovering the most 'attended' phonemes for each class of emotions. We then compare these attended-phonemes in their importance and distribution across acted and natural classes. Our tests show significant differences in the manner and choice of phonemes in acted and natural speech, concluding moderate to low validity and value in using acted speech databases for emotion classification tasks.


title: Analysis of minima for geodesic and chordal cost for a minimal 2D pose-graph SLAM problem

abstract: In this paper, we show that for a minimal pose-graph problem, even in the ideal case of perfect measurements and spherical covariance, using the so-called "wrap function" when comparing angles results in multiple suboptimal local minima. We numerically estimate regions of attraction to these local minima for some numerical examples, and give evidence to show that they are of nonzero measure. In contrast, under the same assumptions, we show that the \textit{chordal distance} representation of angle error has a unique minimum up to periodicity. For chordal cost, we also search for initial conditions that fail to converge to the global minimum, and find that this occurs with far fewer points than with geodesic cost.


title: Real time observation of a stationary magneton

abstract: The magnetic dipole field geometry of subatomic elementary particles like the electron differs from the classical macroscopic field imprint of a bar magnet. It resembles more like an eight figure or else joint double quantum-dots instead of the classical, spherical more uniform field of a bar magnet. This actual subatomic quantum magnetic field of an electron at rest, is called Quantum Magnet or else a Magneton. Normally, a macroscale bar magnet should behave like a relative giant Quantum Magnet with identical magnetic dipole field imprint since all of its individual magnetons collectively inside the material, dipole moments are uniformly aligned forming the total net field of the magnet. However due to Quantum Decoherence (QDE) phenomenon at the macroscale and macroscopic magnetic field imaging sensors limitations which cannot pickup these rapid quantum magnetization fluctuations, this field is masked and not visible at the macroscale. By using the relative inexpensive submicron resolution Ferrolens quantum magnetic optical sensor and method, we can actually make this net magneton field visible on macroscale magnets. We call this net total field herein, Quantum Field of Magnet (QFM) differentiating it therefore from the field of the single subatomic magneton thus quantum magnet. Additionally, the unique potential of the Ferrolens device to display also the magnetic flux lines of this macroscopically projected giant Magenton gives us the unique opportunity to study the individual magnetic flux lines geometrical pattern that of a single subatomic magneton. We describe this particular magnetic flux of the magneton observed, quantum magnetic flux. Therefore a novel observation has been made that the QFM of the Magnet-Magneton consists of a dipole vortex shaped magnetic flux geometrical pattern responsible for creating the classical macroscopic N-S field of the magnet.


title: Constraining the abundance of primordial black holes with gravitational lensing of gravitational waves at LIGO frequencies

abstract: Gravitational waves from binary black holes that are gravitationally lensed can be distorted by small microlenses along the line of sight. Microlenses with masses of a few tens of solar masses, and that are close to a critical curve in the lens plane, can introduce a time delay of a few millisecond. Such time delay would result in distinctive interference patterns in the gravitational wave that can be measured with current experiments such as LIGO/Virgo. We consider the particular case of primordial black holes with masses between 5 and 50 solar masses acting as microlenses. We study the effect of a population of primordial black holes constituting a fraction of the dark matter, and contained in a macrolens (galaxy or cluster), over gravitational waves that are being lensed by the combined effect of the macrolens plus microlenses. We find that at the typical magnifications expected for observed GW events, the fraction of dark matter in the form of compact microlenses, such as primordial black holes, can be constrained to percent level. Similarly, if a small percentage of the dark matter is in the form of microlenses with a few tens of solar masses, at sufficiently large magnification factors, all gravitational waves will show interference effects. These effects could have an impact on the inferred parameters. The effect is more important for macroimages with negative parity, which usually arrive after the macroimages with positive parity.


title: Dynamical axion misalignment with small instantons

abstract: We present a new mechanism to relax the initial misalignment angle of the QCD axion and raise the cosmological bound on the axion decay constant. The QCD axion receives a contribution from small UV instantons during inflation, which raises its mass to the inflationary Hubble scale. This makes the axion start rolling down its potential early on. In the scenario, the standard model Yukawa couplings of quarks are dynamical, being of order one during the inflationary era and reducing to their standard model values once it ends. This means that after inflation the contribution of the small instantons is suppressed, and the axion potential reduces to the standard one from the usual IR instantons. As a result, when the axion starts to oscillate again after inflation, the initial misalignment angle is suppressed due to the dynamics during inflation. While the general idea of dynamical axion misalignment has been discussed in the literature before, we present in detail the major bottleneck on the mismatching between the minima of the axion potentials during and after inflation, and how it is circumvented in our scenario via the Froggatt-Nielsen mechanism. Taking into account of all the constraints, we find that the axion decay constant could be raised to the GUT scale, $10^{15}$ GeV, in our scenario.


title: Axion Dark Matter, Proton Decay and Unification

abstract: We discuss the possibility to predict the QCD axion mass in the context of grand unified theories. We investigate the implementation of the DFSZ mechanism in the context of renormalizable SU(5) theories. In the simplest theory, the axion mass can be predicted with good precision in the range $m_a = (2-16)$ neV, and there is a strong correlation between the predictions for the axion mass and proton decay rates. In this context, we predict an upper bound for the proton decay channels with antineutrinos, $\tau(p\to K^+ \bar{\nu}) \lesssim 4 \times 10^{37} \text{ yr}$ and $\tau(p \to \pi^+ \bar{\nu}) \lesssim 2 \times 10^{36}\text{ yr}$. This theory can be considered as the minimal realistic grand unified theory with the DFSZ mechanism and it can be fully tested by proton decay and axion experiments.


title: Flat space holography and complex SYK

abstract: We provide the first steps towards a flat space holographic correspondence in two bulk spacetime dimensions. The gravity side is described by a conformally transformed version of the matterless Callan-Giddings-Harvey-Strominger model. The field theory side follows from the complex Sachdev-Ye-Kitaev model in the limit of large specific heat and vanishing compressibility. We derive the boundary action analogous to the Schwarzian as the key link between gravity and field theory sides and show that it coincides with a geometric action discovered recently by one of us, see 1908.08089.


title: Hierarchy of relaxation timescales in local random Liouvillians

abstract: To characterize the generic behavior of open quantum systems, we consider random, purely dissipative Liouvillians with a notion of locality. We find that the positivity of the map implies a sharp separation of the relaxation timescales according to the locality of observables. Specifically, we analyze a spin-1/2 system of size $\ell$ with up to $n$-body Lindblad operators, which are $n$-local in the complexity-theory sense. Without locality ($n=\ell$), the complex Liouvillian spectrum densely covers a "lemon"-shaped support, in agreement with recent findings [Phys. Rev. Lett. 123, 140403;arXiv:1905.02155]. However, for local Liouvillians ($n<\ell$), we find that the spectrum is composed of several dense clusters with random matrix spacing statistics, each featuring a lemon-shaped support wherein all eigenvectors correspond to $n$-body decay modes. This implies a hierarchy of relaxation timescales of $n$-body observables, which we verify to be robust in the thermodynamic limit.


title: From VOAs to short star products in SCFT

abstract: We build a bridge between two algebraic structures in SCFT: a VOA in the Schur sector of 4d $\mathcal{N}=2$ theories and an associative algebra in the Higgs sector of 3d $\mathcal{N}=4$. The natural setting is a 4d $\mathcal{N}=2$ SCFT placed on $S^3\times S^1$: by sending the radius of $S^1$ to zero, we recover the 3d $\mathcal{N}=4$ theory, and the corresponding VOA on the torus degenerates to the associative algebra on the circle. We prove that: 1) the Higgs branch operators remain in the cohomology; 2) all the Schur operators of the non-Higgs type are lifted by line operators wrapped on the $S^1$; 3) no new cohomology classes are added. We show that the algebra in 3d is given by the quotient $\mathcal{A}_H = {\rm Zhu}_{s}(V)/N$, where ${\rm Zhu}_{s}(V)$ is the non-commutative Zhu algebra of the VOA $V$ (for ${s}\in{\rm Aut}(V)$), and $N$ is a certain ideal. This ideal is the null space of the (${s}$-twisted) trace map $T_{s}: {\rm Zhu}_{s}(V) \to \mathbb{C}$ determined by the torus 1-point function in the high temperature (or small complex structure) limit. It therefore equips $\mathcal{A}_H$ with a non-degenerate (twisted) trace, leading to a short star-product according to the recent results of Etingof and Stryker. The map $T_{s}$ is easy to determine for unitary VOAs, but has a much subtler structure for non-unitary and non-$C_2$-cofinite VOAs of our interest. We comment on relation to the Beem-Rastelli conjecture on the Higgs branch and the associated variety. A companion paper will explore further details, examples, and some applications of these ideas.


title: Quantum percolation of monopole paths and the response of quantum spin ice

abstract: We consider quantum spin ice in a temperature regime in which its response is dominated by the coherent motion of a dilute gas of monopoles. The hopping amplitude of a monopole is sensitive to the configuration of its surrounding spins, taken to be quasi-static on the relevant timescales. This leads to well-known blocked directions in the monopole motion; we find that these are sufficient to reduce the coherent propagation of monopoles to quantum diffusion. This result is robust against disorder, as a direct consequence of the ground-state degeneracy, which disrupts the quantum interference processes needed for weak localization. Moreover, recent work [Tomasello et al., Phys. Rev. Lett. 123, 067204 (2019)] has shown that the monopole hopping amplitudes are roughly bimodal: for $\approx 1/3$ of the flippable spins surrounding a monopole, these amplitudes are extremely small. We exploit this structure to construct a theory of quantum monopole motion in spin ice. In the limit where the slow hopping terms are set to zero, the monopole wavefunctions appear to be fractal; we explain this observation via a mapping to quantum percolation on trees. The fractal, non-ergodic nature of monopole wavefunctions manifests itself in the low-frequency behavior of monopole spectral functions, and is consistent with experimental observations.


title: High Resolution Imaging in the Visible with Faint Reference Stars on Large Ground-Based Telescopes

abstract: Astronomers working with faint targets will benefit greatly from improved image quality on current and planned ground-based telescopes. At present, most adaptive optic systems are targeted at the highest resolution with bright guide stars. We demonstrate a significantly new approach to measuring low-order wavefront errors by using a pupil-plane curvature wavefront sensor design. By making low order wavefront corrections we can deliver significant improvements in image resolution in the visible on telescopes in the 2.5m to 8.2m range on good astronomical sites. As a minimum the angular resolution will be improved by a factor of 2.5 to 3 under any reasonable conditions and, with further correction and image selection, even sharper images may be obtained routinely. We re-examine many of the assumptions about what may be achieved with faint reference stars to achieve this performance. We show how our new design of curvature wavefront sensor combined with wavefront fitting routines based on radon transforms allow this performance to be achieved routinely. Simulations over a wide range of conditions match the performance already achieved in runs with earlier versions of the hardware described. Reference stars significantly fainter than I 17m may be used routinely to produce images with a near diffraction limited core and halo much smaller than that delivered by natural seeing.


title: Zodiacal Exoplanets in Time (ZEIT) IX: a flat transmission spectrum and a highly eccentric orbit for the young Neptune K2-25b as revealed by Spitzer

abstract: Transiting planets in nearby young clusters offer the opportunity to study the atmospheres and dynamics of planets during their formative years. To this end, we focused on K2-25b -- a close-in ($P$=3.48 days), Neptune-sized exoplanet orbiting a M4.5 dwarf in the 650Myr Hyades cluster. We combined photometric observations of K2-25 covering a total of 44 transits and spanning >2 yr, drawn from a mix of space-based telescopes (Spitzer Space Telescope and K2) and ground-based facilities (Las Cumbres Observatory Global Telescope network and MEarth). The transit photometry spanned 0.6--4.5$\mu$m, which enabled our study of K2-25b's transmission spectrum. We combined and fit each dataset at a common wavelength within a Markov Chain Monte Carlo framework, yielding consistent planet parameters. The resulting transit depths ruled out a solar-composition atmosphere for K2-25b for the range of expected planetary masses and equilibrium temperature at a $>4\sigma$ confidence level, and are consistent with a flat transmission spectrum. Mass constraints and transit observations at a finer grid of wavelengths (e.g., from the Hubble Space Telescope) are needed to make more definitive statements about the presence of clouds or an atmosphere of high mean molecular weight. Our precise measurements of K2-25b's transit duration also enabled new constraints on the eccentricity of K2-25's orbit. We find K2-25b's orbit to be eccentric ($e>0.20$) for all reasonable stellar densities and independent of the observation wavelength or instrument. The high eccentricity is suggestive of a complex dynamical history and motivates future searches for additional planets or stellar companions.


title: A closer look at the spur, blob, wiggle, and gaps in GD-1

abstract: The GD-1 stream is one of the longest and coldest stellar streams discovered to date, and one of the best objects for constraining the dark matter properties of the Milky Way. Using data from {\it Gaia} DR2 we study the proper motions, distance, morphology and density of the stream to uncover small scale perturbations. The proper motion cleaned data shows a clear distance gradient across the stream, ranging from 7 to 12 kpc. However, unlike earlier studies that found a continuous gradient, we uncover a distance minimum at $\varphi_{1}\approx$-50 deg, after which the distance increases again. We can reliably trace the stream between -85$<\varphi_{1}<$15 deg, showing an even further extent to GD-1 beyond the earlier extension of \citet{Price-Whelan18a}. We constrain the stream track and density using a Boolean matched filter approach and find three large under densities and find significant residuals in the stream track lining up with these gaps. In particular, a gap is visible at $\varphi_{1}$=-3 deg, surrounded by a clear sinusoidal wiggle. We argue that this wiggle is due to a perturbation since it has the wrong orientation to come from a progenitor. We compute a total initial stellar mass of the stream segment of 1.58$\pm$0.07$\times$10$^{4}$ M$_{\odot}$. With the extended view of the spur in this work, we argue that the spur may be unrelated to the adjacent gap in the stream. Finally, we show that an interaction with the Sagittarius dwarf can create features similar to the spur.


title: Galaxy Formation with BECDM -- II. Cosmic Filaments and First Galaxies

abstract: Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.


title: Mixed Scalar-Current bootstrap in three dimensions

abstract: We study the mixed system of correlation functions involving a scalar field charged under a global $U(1)$ symmetry and the associated conserved spin-1 current $J_\mu$. Using numerical bootstrap techniques we obtain bounds on new observables not accessible in the usual scalar bootstrap. We then specialize to the $O(2)$ model and extract rigorous bounds on the three-point function coefficient of two currents and the unique relevant scalar singlet, as well as those of two currents and the stress tensor. Using these results, and comparing with a quantum Monte Carlo simulation of the $O(2)$ model conductivity, we give estimates of the thermal one-point function of the relevant singlet and the stress tensor. We also obtain new bounds on operators in various sectors.


title: Stellar population models based on the SDSS-IV MaStar library of stellar spectra. I. Intermediate-age/old models

abstract: We use the first release of the SDSS/MaStar stellar library comprising ~9000, high S/N spectra, to calculate integrated spectra of stellar population models. The models extend over the wavelength range 0.36-1.03 micron and share the same spectral resolution (R~1800) and flux calibration as the SDSS-IV/MaNGA galaxy data. The parameter space covered by the stellar spectra collected thus far allows the calculation of models with ages and chemical composition in the range t>200 Myr, -2 <=[Z/H]<= + 0.35, which will be extended as MaStar proceeds. Notably, the models include spectra for dwarf Main Sequence stars close to the core H-burning limit, as well as spectra for cold, metal-rich giants. Both stellar types are crucial for modelling lambda>0.7 micron absorption spectra. Moreover, a better parameter coverage at low metallicity allows the calculation of models as young as 500 Myr and the full account of the Blue Horizontal Branch phase of old populations. We present models adopting two independent sets of stellar parameters (T_eff, logg, [Z/H]). In a novel approach, their reliability is tested 'on the fly' using the stellar population models themselves. We perform tests with Milky Way and Magellanic Clouds globular clusters, finding that the new models recover their ages and metallicities remarkably well, with systematics as low as a few per cent for homogeneous calibration sets. We also fit a MaNGA galaxy spectrum, finding residuals of the order of a few per cent comparable to the state-of-art models, but now over a wider wavelength range.


title: Three-body losses of a polarized Fermi gas near a p-wave Feshbach resonance in effective field theory

abstract: We study three-body recombination of fully spin-polarized ${}^6$Li atoms that are interacting resonantly in relative p-waves. Motivated by a recent experiment, we focus on negative scattering volumes where three atoms recombine into a deep dimer and another atom. We calculate the three-body recombination rate using a Faddeev equation derived from effective field theory. In particular, we study the magnetic field and temperature dependences of the loss rate and use the recombination data to determine the effective range of the p-wave atom-atom interaction. We also predict the existence of a shallow three-body bound state that manifests itself as a prominent feature in the energy-dependent three-body recombination rate.


title: N\"other's Second Theorem as an Obstruction to Charge Quantization

abstract: While it is a standard result in field theory that the scaling dimension of conserved currents and their associated gauge fields are determined strictly by dimensional analysis and hence cannot change under any amount of renormalization, it is also the case that the standard conservation laws for currents, $dJ=0$, remain unchanged in form if any differential operator that commutes with the total exterior derivative, $[d,\hat Y]=0$, multiplies the current. Such an operator, effectively changing the dimension of the current, increases the allowable gauge transformations in electromagnetism and is at the heart of N\"other's second theorem. We review here our recent work on one particular instance of this theorem, namely fractional electromagnetism and highlight the holographic dilaton models that exhibit such behavior and the physical consequences this theory has for charge quantization. Namely, the standard electromagnetic gauge and the fractional counterpart cannot both yield integer values of Planck's constant when they are integrated around a closed loop, thereby leading to a breakdown of charge quantization.


title: Cosmology dependence of galaxy cluster scaling relations

abstract: The abundance of galaxy clusters as a function of mass and redshift is a well known powerful cosmological probe, which relies on underlying modelling assumptions on the mass-observable relations (MOR). Some of the MOR parameters can be constrained directly from multi-wavelength observations, as the normalization at some reference cosmology, the mass-slope, the redshift evolution and the intrinsic scatter. However, the cosmology dependence of MORs cannot be tested with multi-wavelength observations alone. We use {\tt Magneticum} simulations to explore the cosmology dependence of galaxy cluster scaling relations. We run fifteen hydro-dynamical cosmological simulations varying $\Omega_m$, $\Omega_b$, $h_0$ and $\sigma_8$ (around a reference cosmological model). The MORs considered are gas mass, baryonic mass, gas temperature, $Y$ and velocity dispersion as a function of virial mass. We verify that the mass and redshift slopes and the intrinsic scatter of the MORs are nearly independent of cosmology with variations significantly smaller than current observational uncertainties. We show that the gas mass and baryonic mass sensitively depends only on the baryon fraction, velocity dispersion and gas temperature on $h_0$, and $Y$ on both baryon fraction and $h_0$. We investigate the cosmological implications of our MOR parameterization on a mock catalog created for an idealized eROSITA-like experiment. We show that our parametrization introduces a strong degeneracy between the cosmological parameters and the normalization of the MOR. Finally, the parameter constraints derived at different overdensity ($\Delta_{500c}$), for X-ray bolometric gas luminosity, and for different subgrid physics prescriptions are shown in the appendix.


title: Adaptive filtering of projective quantum measurements using discrete stochastic methods

abstract: Adaptive filtering is a powerful class of control theoretic concepts useful in extracting information from noisy data sets or performing forward prediction in time for a dynamic system. The broad utilization of the associated algorithms makes them attractive targets for similar problems in the quantum domain. To date, however, the construction of adaptive filters for quantum systems has typically been carried out in terms of stochastic differential equations for weak, continuous quantum measurements, as used in linear quantum systems such as optical cavities. Discretized measurement models are not as easily treated in this framework, but are frequently employed in quantum information systems leveraging projective measurements. This paper presents a detailed analysis of several technical innovations that enable classical filtering of discrete projective measurements, useful for adaptively learning system-dynamics, noise properties, or hardware performance variations in classically correlated measurement data from quantum devices. In previous work we studied a specific case of this framework, in which noise and calibration errors on qubit arrays could be efficiently characterized in space; here, we present a generalized analysis of filtering in quantum systems and demonstrate that the traditional convergence properties of nonlinear classical filtering hold using single-shot projective measurements. These results are important early demonstrations indicating that a range of concepts and techniques from classical nonlinear filtering theory may be applied to the characterization of quantum systems involving discretized projective measurements, paving the way for broader adoption of control theoretic techniques in quantum technology.


title: The Cloud Factory I: Generating resolved filamentary molecular clouds from galactic-scale forces

abstract: We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitational potential for computational efficiency. The simulations have nested levels of resolution, with the lowest layer tied to tracer particles injected into individual cloud complexes. These tracer refinement regions are embedded in the larger simulation so continue to experience forces from outside the cloud. This allows the simulations to act as a laboratory for testing the effect of galactic environment on star formation. Here we introduce our method and investigate the effect of galactic environment on filamentary clouds. We find that cloud complexes formed after a clustered burst of feedback, have shorter lengths and are less likely to fragment compared to quiescent clouds (e.g. the Musca filament) or those dominated by the galactic potential (e.g. Nessie). Spiral arms and differential rotation preferentially align filaments, but strong feedback randomises them. Long filaments formed within the cloud complexes are necessarily coherent with low internal velocity gradients, which has implications for the formation of filamentary star-clusters. Cloud complexes formed in regions dominated by supernova feedback have fewer star-forming cores, and these are more widely distributed. These differences show galactic-scale forces can have a significant impact on star formation within molecular clouds.


title: Implicit Hamiltonian Monte Carlo for Sampling Multiscale Distributions

abstract: Hamiltonian Monte Carlo (HMC) has been widely adopted in the statistics community because of its ability to sample high-dimensional distributions much more efficiently than other Metropolis-based methods. Despite this, HMC often performs sub-optimally on distributions with high correlations or marginal variances on multiple scales because the resulting stiffness forces the leapfrog integrator in HMC to take an unreasonably small stepsize. We provide intuition as well as a formal analysis showing how these multiscale distributions limit the stepsize of leapfrog and we show how the implicit midpoint method can be used, together with Newton-Krylov iteration, to circumvent this limitation and achieve major efficiency gains. Furthermore, we offer practical guidelines for when to choose between implicit midpoint and leapfrog and what stepsize to use for each method, depending on the distribution being sampled. Unlike previous modifications to HMC, our method is generally applicable to highly non-Gaussian distributions exhibiting multiple scales. We illustrate how our method can provide a dramatic speedup over leapfrog in the context of the No-U-Turn sampler (NUTS) applied to several examples.


title: An Introduction to Artificial Intelligence and Solutions to the Problems of Algorithmic Discrimination

abstract: There is substantial evidence that Artificial Intelligence (AI) and Machine Learning (ML) algorithms can generate bias against minorities, women, and other protected classes. Federal and state laws have been enacted to protect consumers from discrimination in credit, housing, and employment, where regulators and agencies are tasked with enforcing these laws. Additionally, there are laws in place to ensure that consumers understand why they are denied access to services and products, such as consumer loans. In this article, we provide an overview of the potential benefits and risks associated with the use of algorithms and data, and focus specifically on fairness. While our observations generalize to many contexts, we focus on the fairness concerns raised in consumer credit and the legal requirements of the Equal Credit and Opportunity Act. We propose a methodology for evaluating algorithmic fairness and minimizing algorithmic bias that aligns with the provisions of federal and state anti-discrimination statutes that outlaw overt, disparate treatment, and, specifically, disparate impact discrimination. We argue that while the use of AI and ML algorithms heighten potential discrimination risks, these risks can be evaluated and mitigated, but doing so requires a deep understanding of these algorithms and the contexts and domains in which they are being used.


title: Can astrophysical neutrinos trace the origin of the detected ultra-high energy cosmic rays?

abstract: Since astrophysical neutrinos are produced in the interactions of cosmic rays, identifying the origin of cosmic rays using directional correlations with neutrinos is one of the most interesting possibilities of the field. For that purpose, especially the Ultra-High Energy Cosmic Rays (UHECRs) are promising, as they are deflected less by extragalactic and Galactic magnetic fields than cosmic rays at lower energies. However, photo-hadronic interactions of the UHECRs limit their horizon, while neutrinos do not interact over cosmological distances. We study the possibility to search for anisotropies by investigating neutrino-UHECR correlations from the theoretical perspective, taking into account the UHECR horizon, magnetic-field deflections, and the cosmological source evolution. Under the assumption that the neutrinos and UHECRs all come from the same source class, we demonstrate that the non-observation of neutrino multiplets strongly constrains the possibility to find neutrino-UHECR correlations.


title: Mechanistic insight into the limiting factors of graphene-based environmental sensors

abstract: Graphene has demonstrated great promise for technological use, yet control over material growth and understanding of how material imperfections affect the performance of devices are challenges that hamper the development of applications. In this work we reveal new insight into the connections between the performance of the graphene devices as environmental sensors and the microscopic details of the interactions at the sensing surface. Specifically, we monitor changes in the resistance of the chemical-vapour deposition grown graphene devices as exposed to different concentrations of ethanol. We perform thermal surface treatments after the devices are fabricated, use scanning probe microscopy to visualize their effects on the graphene sensing surface down to nanometer scale and correlate them with the measured performance of the device as an ethanol sensor. Our observations are compared to theoretical calculations of charge transfers between molecules and the graphene surface. We find that, although often overlooked, the surface cleanliness after device fabrication is responsible for the device performance and reliability. These results further our understanding of the mechanisms of sensing in graphene-based environmental sensors and pave the way to optimizing such devices, especially for their miniaturization, as with decreasing size of the active zone the potential role of contaminants will rise.


title: What do you mean, BERT? Assessing BERT as a Distributional Semantics Model

abstract: Contextualized word embeddings, i.e. vector representations for words in context, are naturally seen as an extension of previous noncontextual distributional semantic models. In this work, we focus on BERT, a deep neural network that produces contextualized embeddings and has set the state-of-the-art in several semantic tasks, and study the semantic coherence of its embedding space. While showing a tendency towards coherence, BERT does not fully live up to the natural expectations for a semantic vector space. In particular, we find that the position of the sentence in which a word occurs, while having no meaning correlates, leaves a noticeable trace on the word embeddings and disturbs similarity relationships.


title: Variational circuit compiler for quantum error correction

abstract: Quantum error correction is vital for implementing universal quantum computing. A key component is the encoding circuit that maps a product state of physical qubits into the encoded multipartite entangled logical state. Known methods are typically not 'optimal' either in terms of the circuit depth (and therefore the error burden) or the specifics of the target platform, i.e. the native gates and topology of a system. This work introduces a variational compiler for efficiently finding the encoding circuit of general quantum error correcting codes with given quantum hardware. Focusing on the noisy intermediate scale quantum regime, we show how to systematically compile the circuit following an optimising process seeking to minimise the number of noisy operations that are allowed by the noisy quantum hardware or to obtain the highest fidelity of the encoded state with noisy gates. We demonstrate our method by deriving novel encoders for logic states of the five qubit code and the seven qubit Steane code. We describe ways to augment the discovered circuits with error detection. Our method is applicable quite generally for compiling the encoding circuits of quantum error correcting codes.


title: Flyby-induced misalignments in planet-hosting discs

abstract: We now have several observational examples of misaligned broken protoplanetary discs, where the disc inner regions are strongly misaligned with respect to the outer disc. Current models suggest that this disc structure can be generated with an internal misaligned companion (stellar or planetary), but the occurrence rate of these currently unobserved companions remains unknown. Here we explore whether a strong misalignment between the inner and outer disc can be formed without such a companion. We consider a disc that has an existing gap --- essentially separating the disc into two regions --- and use a flyby to disturb the discs, leading to a misalignment. Despite considering the most optimistic parameters for this scenario, we find maximum misalignments between the inner and outer disc of $\sim$45$^{\circ}$ and that these misalignments are short-lived. We thus conclude that the currently observed misaligned discs must harbour internal, misaligned companions.


title: Predicting Unobserved Space For Planning via Depth Map Augmentation

abstract: Safe and efficient path planning is crucial for autonomous mobile robots. A prerequisite for path planning is to have a comprehensive understanding of the 3D structure of the robot's environment. On MAVs this is commonly achieved using low-cost sensors, such as stereo or RGB-D cameras. These sensors may fail to provide depth measurements in textureless or IR-absorbing areas and have limited effective range. In path planning, this results in inefficient trajectories or failure to recognize a feasible path to the goal, hence significantly impairing the robot's mobility. Recent advances in deep learning enables us to exploit prior experience about the shape of the world and hence to infer complete depth maps from color images and additional sparse depth measurements. In this work, we present an augmented planning system and investigate the effects of employing state-of-the-art depth completion techniques, specifically trained to augment sparse depth maps originating from RGB-D sensors, semi-dense methods and stereo matchers. We extensively evaluate our approach in online path planning experiments based on simulated data, as well as global path planning experiments based on real world MAV data. We show that our augmented system, provided with only sparse depth perception, can reach on-par performance to ground truth depth input in simulated online planning experiments. On real world MAV data the augmented system demonstrates superior performance compared to a planner based on very dense RGB-D depth maps.


title: Interval matrices: realization of ranks by rational matrices

abstract: Let $\alpha $ be a $p \times q$ interval matrix with $p \geq q$ and with the endpoints of all its entries in the set of the rational numbers. We prove that, if $\alpha $ contains a rank-$r$ real matrix with $r \in \{2, q-2,q-1,q\}$, then it contains a rank-$r$ rational matrix.


title: Capturing Non-Gaussianity of the Large-Scale Structure with Weighted Skew-Spectra

abstract: The forthcoming generation of wide-field galaxy surveys will probe larger volumes and galaxy densities, thus allowing for a much larger signal-to-noise ratio for higher-order clustering statistics, in particular the galaxy bispectrum. Extracting this information, however, is more challenging than using the power spectrum due to more complex theoretical modeling, as well as significant computational cost of evaluating the bispectrum signal and the error budget. To overcome these challenges, several proxy statistics have been proposed in the literature, which partially or fully capture the information in the bispectrum, while being computationally less expensive than the bispectrum. One such statistics are {\it weighted skew-spectra}, which are cross-spectra of the density field and appropriately weighted quadratic fields. Using Fisher forecasts, we show that the information in these skew-spectra is equivalent to that in the bispectrum for parameters that appear as amplitudes in the bispectrum model, such as galaxy bias parameters or the amplitude of primordial non-Gaussianity. We consider three shapes of the primordial bispectrum: local, equilateral and that due to massive particles with spin two during inflation. To obtain constraints that match those from a measurement of the full bispectrum, we find that it is crucial to account for the full covariance matrix of the skew-spectra.


title: Portable and wireless signal transducer for field testing of environmental sensors based on 2D materials

abstract: In this paper we present the design and fabrication of a portable device for environmental monitoring applications. This novel hand-held apparatus monitors the changes in the resistance of a sensing surface with a high accuracy and resolution and transmits the recorded data wirelessly to a cellphone. Such a design offers a solution for field testing of environmental sensors. The tested sensing surface in this study is based on an ultrathin material: graphene, which is placed on the surface of a Si/SiO2 wafer. This signal transducer and wireless communication system form together an ideal platform to harvest the sensitivity and selectivity of 2D materials for gas sensing applications.


title: Constraints on the Engines of Fast Radio Bursts

abstract: We model the sample of fast radio bursts (FRB), including the newly discovered CHIME repeaters, using the synchrotron blast wave model of Metzger, Margalit & Sironi (2019). This model postulates that FRBs are precursor radiation from ultra-relativistic magnetized shocks generated as flare ejecta from a central engine collide with an effectively stationary external medium. Downward drifting of the burst frequency structure naturally arises from deceleration of the blast-wave. The data are consistent with FRBs being produced by flares of energy $E_{\rm flare} \sim 10^{43}-10^{46}(f_{\xi}/10^{-3})^{-4/5}$ erg, where $f_{\xi}$ is the maser efficiency, and minimum bulk Lorentz factors $\Gamma \approx 10^2-10^3$, which generate the observed FRBs at shock radii $r_{\rm sh} \sim 10^{12}-10^{13}$ cm. We infer upstream densities $n_{\rm ext}(r_{\rm sh}) \sim 10^{2}-10^{4}$ cm$^{-3}$ and radial profiles $n_{\rm ext} \propto r^{-k}$ showing a range of slopes $k \approx [-2,1]$ (which are seen to evolve between bursts), broadly consistent with the upstream medium being the inner edge of an ion-loaded shell released by a recent energetic flare. The burst timescales, energetics, rates, and external medium properties are consistent with repeating FRBs arising from young, hyper-active flaring magnetars, but the methodology presented is generally applicable to any central engine which injects energy impulsively into a dense magnetized medium. Uncertainties and variations of the model are discussed, including the effects of the strong electric field of the FRB wave (strength parameter $a \gg 1$) on the upstream medium. One-dimensional particle-in-cell simulations of magnetized shocks into a pair plasma are presented which demonstrate that high maser efficiency can be preserved, even in the limit $a \gg 1$ in which the FRB wave accelerates the upstream electrons to ultra-relativistic speeds.


title: Generalized and multiplexed $q$-plates: experimental implementation

abstract: In this paper we generalize the concept of $q$-plate, allowing arbitrary functions of both the radial and the azimuthal variables, and study their effect on uniformly polarized beams in the near and far-field regime. This gives a tool for achieving beams with hybrid states of polarization (SoPs), and alternative phase and intensity distributions. We also implement an experimental device based on a liquid crystal on silicon (LCoS) display for emulating these generalized $q$-plates. Moreover, we propose an application that takes advantage of the pixelated nature of this kind of devices for creating arbitrary superpositions of vector and vortex beams by representing onto the LCoS randomized combinations of two different $q$-plates, i.e. multiplexed $q$-plates. Great agreement is found between theoretical and experimental results.


title: Maximizing the Partial Decode-and-Forward Rate in the Gaussian MIMO Relay Channel

abstract: It is known that circularly symmetric Gaussian signals are the optimal input signals for the partial decode-and-forward (PDF) coding scheme in the Gaussian multiple-input multiple-output (MIMO) relay channel, but there is currently no method to find the optimal covariance matrices nor to compute the optimal achievable PDF rate since the optimization is a non-convex problem in its original formulation. In this paper, we show that it is possible to find a convex reformulation of the problem by means of an approximation and a primal decomposition. We derive an explicit solution for the inner problems as well as an explicit gradient for the outer problem, so that the efficient cutting-plane method can be applied for solving the outer problem. As the accuracy of this provably convergent algorithm might be impaired by the previous approximation, we additionally propose a modified algorithm, for which we cannot give a converge guarantee, but which provides rigorous upper and lower bounds to the optimum of the original rate maximization. In numerical simulations, these bounds become tight in all considered instances of the problem, showing that the proposed method manages to find the global optimum in all these instances.


title: Dynamic mode decomposition based analysis of flow over a sinusoidally pitching airfoil

abstract: Dynamic mode decomposition (DMD) has proven to be a valuable tool for the analysis of complex flow-fields but the application of this technique to flows with moving boundaries is not straightforward. This is due to the difficulty in accounting in the DMD formulation, for a body of non-zero thickness moving through the field of interest. This work presents a method for decomposing the flow on or near a moving boundary by a change of reference frame, followed by a correction to the computed modes that is determined by the frequency spectrum of the motion. The correction serves to recover the modes of the underlying flow dynamics, while removing the effect of change in reference frame. This method is applied to flow over sinusoidally pitching airfoils, and the DMD analysis is used to derive useful insights regarding flow-induced pitch oscillations of these airfoils.


title: A scalable weakly-synchronous algorithm for solving partial differential equations

abstract: Synchronization overheads pose a major challenge as applications advance towards extreme scales. In current large-scale algorithms, synchronization as well as data communication delay the parallel computations at each time step in a time-dependent partial differential equation (PDE) solver. This creates a new scaling wall when moving towards exascale. We present a weakly-synchronous algorithm based on novel asynchrony-tolerant (AT) finite-difference schemes that relax synchronization at a mathematical level. We utilize remote memory access programming schemes that have been shown to provide significant speedup on modern supercomputers, to efficiently implement communications suitable for AT schemes, and compare to two-sided communications that are state-of-practice. We present results from simulations of Burgers' equation as a model of multi-scale strongly non-linear dynamical systems. Our algorithm demonstrate excellent scalability of the new AT schemes for large-scale computing, with a speedup of up to $3.3$x in communication time and $2.19$x in total runtime. We expect that such schemes can form the basis for exascale PDE algorithms.


title: Constrained Bayesian ICA for Brain Connectome Inference

abstract: Brain connectomics is a developing field in neurosciences which strives to understand cognitive processes and psychiatric diseases through the analysis of interactions between brain regions. However, in the high-dimensional, low-sample, and noisy regimes that typically characterize fMRI data, the recovery of such interactions remains an ongoing challenge: how can we discover patterns of co-activity between brain regions that could then be associated to cognitive processes or psychiatric disorders? In this paper, we investigate a constrained Bayesian ICA approach which, in comparison to current methods, simultaneously allows (a) the flexible integration of multiple sources of information (fMRI, DWI, anatomical, etc.), (b) an automatic and parameter-free selection of the appropriate sparsity level and number of connected submodules and (c) the provision of estimates on the uncertainty of the recovered interactions. Our experiments, both on synthetic and real-life data, validate the flexibility of our method and highlight the benefits of integrating anatomical information for connectome inference.


title: Machine Learning Based Network Vulnerability Analysis of Industrial Internet of Things

abstract: It is critical to secure the Industrial Internet of Things (IIoT) devices because of potentially devastating consequences in case of an attack. Machine learning and big data analytics are the two powerful leverages for analyzing and securing the Internet of Things (IoT) technology. By extension, these techniques can help improve the security of the IIoT systems as well. In this paper, we first present common IIoT protocols and their associated vulnerabilities. Then, we run a cyber-vulnerability assessment and discuss the utilization of machine learning in countering these susceptibilities. Following that, a literature review of the available intrusion detection solutions using machine learning models is presented. Finally, we discuss our case study, which includes details of a real-world testbed that we have built to conduct cyber-attacks and to design an intrusion detection system (IDS). We deploy backdoor, command injection, and Structured Query Language (SQL) injection attacks against the system and demonstrate how a machine learning based anomaly detection system can perform well in detecting these attacks. We have evaluated the performance through representative metrics to have a fair point of view on the effectiveness of the methods.


title: ADMX SLIC: Results from a Superconducting LC Circuit Investigating Cold Axions

abstract: Axions are a promising cold dark matter candidate. Haloscopes, which use the conversion of axions to photons in the presence of a magnetic field to detect axions, are the basis of microwave cavity searches such as the Axion Dark Matter eXperiment (ADMX). To search for lighter, low frequency axions in the sub $2\times10^{-7}$ eV (50 MHz) range, a tunable lumped-element LC circuit has been proposed. For the first time, through ADMX SLIC (Superconducting Lc circuit Investigating Cold axions), a resonant LC circuit was used to probe this region of axion mass-coupling space. The detector used a superconducting LC circuit with piezoelectric driven capacitive tuning. The axion mass and corresponding frequency range $1.7498 -1.7519 \times10^{-7}$ eV (42.31 -- 42.36 MHz), $1.7734 - 1.7738 \times10^{-7}$ eV (42.88 -- 42.89 MHz), and $1.8007 - 1.8015 \times10^{-7}$ eV (43.54 -- 43.56 MHz) was covered at magnetic fields of 4.5 T, 5.0 T, and 7.0 T respectively. Exclusion results from the search data, for coupling below $10^{-12} \text{GeV}^{-1}$ are presented.


title: On the construction of the linearization of a nonlinear connection

abstract: The construction of a linear connection on a pullback bundle from a connection on a vector bundle is explained in terms of fiberwise linear approximation. This procedure clarifies the geometric meaning of the linearized connection as well as the associated parallel transport and curvature.


title: Factor Group-Sparse Regularization for Efficient Low-Rank Matrix Recovery

abstract: This paper develops a new class of nonconvex regularizers for low-rank matrix recovery. Many regularizers are motivated as convex relaxations of the matrix rank function. Our new factor group-sparse regularizers are motivated as a relaxation of the number of nonzero columns in a factorization of the matrix. These nonconvex regularizers are sharper than the nuclear norm; indeed, we show they are related to Schatten-$p$ norms with arbitrarily small $0 < p \leq 1$. Moreover, these factor group-sparse regularizers can be written in a factored form that enables efficient and effective nonconvex optimization; notably, the method does not use singular value decomposition. We provide generalization error bounds for low-rank matrix completion which show improved upper bounds for Schatten-$p$ norm reglarization as $p$ decreases. Compared to the max norm and the factored formulation of the nuclear norm, factor group-sparse regularizers are more efficient, accurate, and robust to the initial guess of rank. Experiments show promising performance of factor group-sparse regularization for low-rank matrix completion and robust principal component analysis.


title: On the Relativized Alon Second Eigenvalue Conjecture VI: Sharp Bounds for Ramanujan Base Graphs

abstract: This is the sixth in a series of articles devoted to showing that a typical covering map of large degree to a fixed, regular graph has its new adjacency eigenvalues within the bound conjectured by Alon for random regular graphs. In this article we show that if the fixed graph is regular Ramanujan, then the {\em algebraic power} of the model of random covering graphs is $+\infty$. This implies a number of interesting results, such as (1) one obtains the upper and lower bounds---matching to within a multiplicative constant---for the probability that a random covering map has some new adjacency eigenvalue outside the Alon bound, and (2) with probability smaller than any negative power of the degree of the covering map, some new eigenvalue fails to be within the Alon bound without the covering map containing one of finitely many "tangles" as a subgraph (and this tangle containment event has low probability).


title: On the Upsilon invariant of fibered knots and right-veering open books

abstract: We give a sufficient condition using the Ozsv\'ath-Stipsicz-Szab\'o concordance invariant Upsilon for the monodromy of the open book decomposition of a fibered knot to be right-veering. As an application, we generalize a result of Baker on ribbon concordances between fibered knots. Following Baker, we conclude that either fibered knots $K$ in $S^{3}$ satisfying that $\Upsilon'(t) = -g(K)$ for some $t \in [0,1)$ are unique in their smooth concordance classes or there exists a counterexample to the Slice-Ribbon Conjecture.


title: Optimization Models for Estimating Transit Network Origin-Destination Flows with AVL/APC Data

abstract: With the increasing adoption of Automatic Vehicle Location (AVL) and Automatic Passenger Count (APC) technologies by transit agencies, a massive amount of time-stamped and location-based passenger boarding and alighting count data can be collected on a continuous basis. The availability of such large-scale transit data offers new opportunities to produce estimates for Origin-Destination (O-D) flows, helping inform transportation planning and transit management. However, the state-of-the-art methodologies for AVL/APC data analysis mostly tackle the O-D flow estimation problem within routes and barely infer the transfer activities across the entire transit network. This paper proposes three optimization models to identify transfers and approximate network-level O-D flows by minimizing the deviations between estimated and observed proportions or counts of transferring passengers: A Quadratic Integer Program (QIP), a feasible rounding procedure for the Quadratic Convex Programming (QCP) relaxation of the QIP, and an Integer Program (IP). The inputs of the models are readily available by applying the various route-level flow estimation algorithms to the automatically collected AVL/APC data and the output of the models is a network O-D estimation at varying geographical resolutions. The optimization models were evaluated on a case study for Ann Arbor-Ypsilanti area in Michigan. The IP model outperforms the QCP approach in terms of accuracy and remains tractable from an efficiency standpoint, contrary to the QIP. Its estimated O-D matrix achieves an R-Squared metric of 95.57% at the Traffic Analysis Zone level and 92.39% at the stop level, compared to the ground-truth estimates inferred from the state-of-practice trip-chaining methods.


title: Comments on the CKN Bound

abstract: Cohen, Kaplan, and Nelson (CKN) conjectured that the UV and IR cutoffs of effective quantum field theories coupled to gravity are not independent, but are connected by the physics of black holes. We interpret the CKN bound as a scale-dependent depletion of the QFT density of states and discuss various aspects of the bound on small and large scales. For laboratory experiments, we argue that the bound provides small corrections to ordinary quantum field theory, which we estimate to be of order $m_e/M_p$ for $g-2$ of the electron. On large scales, we suggest a modification of the CKN bound due to the presence of cosmological horizons and discuss the connection with entropy bounds.


title: The Number of Threshold Words on $n$ Letters Grows Exponentially for Every $n\geq 27$

abstract: For every $n\geq 27$, we show that the number of $n/(n-1)^+$-free words (i.e., threshold words) of length $k$ on $n$ letters grows exponentially in $k$. This settles all but finitely many cases of a conjecture of Ochem.


title: Kolmogorovian active turbulence of a sparse assembly of interacting swimmers

abstract: Active matter, composed of self-propelled entities, forms a wide class of out-of-equilibrium systems that display striking collective behaviors among which the so-called active turbulence where spatially and time disordered flow patterns spontaneously arise in a variety of {active systems}. De facto, the active turbulence naming suggests a connection with a second seminal class of out-of-equilibrium systems, fluid turbulence, and yet of very different nature with energy injected at global system scale rather than at the elementary scale of single constituents. Indeed the existence of a possible strong-tie between active and canonical turbulence remains an open question and a field of profuse research. Using an assembly of self-propelled interfacial particles, we show experimentally that this active system shares remarkable quantitative similarities with canonical fluid turbulence, as described by the celebrated 1941 phenomenology of Kolmogorov. Making active matter entering into the universality class of fluid turbulence not only benefits to its future development but may also provide new insights for the longstanding description of turbulent flows, arguably one of the biggest remaining mysteries in classical physics.


title: Learning Internal Representations (COLT 1995)

abstract: Probably the most important problem in machine learning is the preliminary biasing of a learner's hypothesis space so that it is small enough to ensure good generalisation from reasonable training sets, yet large enough that it contains a good solution to the problem being learnt. In this paper a mechanism for {\em automatically} learning or biasing the learner's hypothesis space is introduced. It works by first learning an appropriate {\em internal representation} for a learning environment and then using that representation to bias the learner's hypothesis space for the learning of future tasks drawn from the same environment. An internal representation must be learnt by sampling from {\em many similar tasks}, not just a single task as occurs in ordinary machine learning. It is proved that the number of examples $m$ {\em per task} required to ensure good generalisation from a representation learner obeys $m = O(a+b/n)$ where $n$ is the number of tasks being learnt and $a$ and $b$ are constants. If the tasks are learnt independently ({\em i.e.} without a common representation) then $m=O(a+b)$. It is argued that for learning environments such as speech and character recognition $b\gg a$ and hence representation learning in these environments can potentially yield a drastic reduction in the number of examples required per task. It is also proved that if $n = O(b)$ (with $m=O(a+b/n)$) then the representation learnt will be good for learning novel tasks from the same environment, and that the number of examples required to generalise well on a novel task will be reduced to $O(a)$ (as opposed to $O(a+b)$ if no representation is used). It is shown that gradient descent can be used to train neural network representations and experiment results are reported providing strong qualitative support for the theoretical results.


title: Prime stars multiplexes

abstract: This work investigates the class of prime star multiplexes, in which each of its layers $i$, $i=1,2, \ldots M$, consists of a regular cycle graph where any node has $2J_i$ neighbors. In a process that does not affect the cyclic topology, it is assumed that, before the multiplex is assembled, the nodes are labeled differently in each individual layer. As the setup requires that all representations of the same node in the different $M$ layers must be linked by inter-layers connections, the resulting multiplex pattern can be highly complex. This can be better visualized if one assumes that in one layer the nodes are labeled in the sequentially ascending order and that the nodes with the same label are drawn on the top of the other, so that all inter-layer connections are represented by vertical lines. In such cases, the other $M-1$ layers are characterized by long distance shortcuts. As a consequence, in spite of sharing the same internal topological structure, the multiplex ends up with very dissimilar layers. For prime number of nodes, a regular star geometry arises by requiring that the neighbor labels of the $M-1$ layers differ by a constant value $p_i>1$. For $M=2$, we use analytical and numerical approaches to provide a thorough characterization of the multiplex topological properties, of the inter layer dissimilarity, and of the diffusive dynamical processes taking place on them. For the sake of definitiveness, it is considered that each node in the sequentially labeled layer is characterized by $J_1\geq1$. In the other layer, we fix $J_2\equiv1$, while $p>1$ becomes a proxy of layer dissimilarity.


title: AdS Q-Soliton and Inhomogeneously mass-deformed ABJM Model

abstract: We study dual geometries to a deformed ABJM model with spatially dependent source functions at finite temperature. These source functions are proportional to the mass function $m(x)= m_0 \sin k x$ and its derivative $m'(x)$. As dual geometries, we find hairy black branes and AdS solitons corresponding to deconfinement phase and confining phase of the dual field theory, respectively. It turns out that the hairy AdS solitons have lower free energy than the black branes when the Hawking temperature is smaller than the confining scale. Therefore the dual system undergoes the first order phase transition. Even though our study is limited to the so-called Q-lattice ansatz, the solution space contains a set of solutions dual to a supersymmetric mass deformation. As a physical quantity to probe the confining phase, we investigate the holographic entanglement entropy and discuss its behavior in terms of modulation effect.


title: Pricing Multi-Interval Dispatch under Uncertainty Part I: Dispatch-Following Incentives

abstract: Pricing multi-interval economic dispatch of electric power under operational uncertainty is considered in this two-part paper. Part I investigates dispatch-following incentives of profit-maximizing generators and shows that, under mild conditions, no uniform-pricing scheme for the rolling-window economic dispatch provides dispatch-following incentives that avoid discriminative out-of-the-market uplifts. A nonuniform pricing mechanism, referred to as the temporal locational marginal pricing (TLMP), is proposed. As an extension of the standard locational marginal pricing (LMP), TLMP takes into account both generation and ramping-induced opportunity costs. It eliminates the need for the out-of-the-market uplifts and guarantees full dispatch-following incentives regardless of the accuracy of the demand forecasts used in the dispatch. It is also shown that, under TLMP, a price-taking market participant has incentives to bid truthfully with its marginal cost of generation. Part II of the paper extends the theoretical results developed in Part I to more general network settings. It investigates a broader set of performance measures, including the incentives of the truthful revelation of ramping limits, revenue adequacy of the operator, consumer payments, generator profits, and price volatility under the rolling-window dispatch model with demand forecast errors.


title: Regular orbits of quasisimple linear groups I

abstract: Let $G \leq \mathrm{GL}(V)$ be a group with a unique subnormal quasisimple subgroup $E(G)$ that acts absolutely irreducibly on $V$. A base for $G$ acting on $V$ is a set of vectors with trivial pointwise stabiliser in $G$. In this paper we determine the minimal base size of $G$ when $E(G)/Z(E(G))$ is a finite simple group of Lie type in cross-characteristic. We show that $G$ has a regular orbit on $V$, with specific exceptions, for which we find the base size.


title: Spectroscopic tests for short-range modifications of Newtonian and post-Newtonian potentials

abstract: There are theoretical frameworks, such as the large extra dimension models, which predict the strengthening of the gravitational field in short distances. Here we obtain new empiric constraints for deviations of standard gravity in the atomic length scale from analyses of recent and accurate data of hydrogen spectroscopy. The new bounds, extracted from 1S-3S transition, are compared with previous limits given by antiprotonic Helium spectroscopy. Independent constraints are also determined by investigating the effects of gravitational spin-orbit coupling on the atomic spectrum. We show that the analysis of the influence of that interaction, which is responsible for the spin precession phenomena, on the fine structure of the states can be employed as a test of a post-Newtonian potential in the atomic domain. The constraints obtained here from 2P_{1/2}-2P_{3/2} transition in hydrogen are tighter than previous bounds determined from measurements of the spin precession in an electron-nucleus scattering.


title: Visual-Inertial Localization for Skid-Steering Robots with Kinematic Constraints

abstract: While visual localization or SLAM has witnessed great progress in past decades, when deploying it on a mobile robot in practice, few works have explicitly considered the kinematic (or dynamic) constraints of the real robotic system when designing state estimators. To promote the practical deployment of current state-of-the-art visual-inertial localization algorithms, in this work we propose a low-cost kinematics-constrained localization system particularly for a skid-steering mobile robot. In particular, we derive in a principle way the robot's kinematic constraints based on the instantaneous centers of rotation (ICR) model and integrate them in a tightly-coupled manner into the sliding-window bundle adjustment (BA)-based visual-inertial estimator. Because the ICR model parameters are time-varying due to, for example, track-to-terrain interaction and terrain roughness, we estimate these kinematic parameters online along with the navigation state. To this end, we perform in-depth the observability analysis and identify motion conditions under which the state/parameter estimation is viable. The proposed kinematics-constrained visual-inertial localization system has been validated extensively in different terrain scenarios.


title: Computing Equilibria in Binary Networked Public Goods Games

abstract: Public goods games study the incentives of individuals to contribute to a public good and their behaviors in equilibria. In this paper, we examine a specific type of public goods game where players are networked and each has binary actions, and focus on the algorithmic aspects of such games. First, we show that checking the existence of a pure-strategy Nash equilibrium is NP-complete. We then identify tractable instances based on restrictions of either utility functions or of the underlying graphical structure. In certain cases, we also show that we can efficiently compute a socially optimal Nash equilibrium. Finally, we propose a heuristic approach for computing approximate equilibria in general binary networked public goods games, and experimentally demonstrate its effectiveness.


title: Generation and detection of coherent longitudinal acoustic waves in ultrathin 1T'-MoTe2

abstract: Layered transition metal dichalcogenides have attracted substantial attention owing to their versatile functionalities and compatibility with current nanofabrication technologies. Thus, noninvasive means to determine the mechanical properties of nanometer (nm) thick specimens are of increasing importance. Here, we report on the detection of coherent longitudinal acoustic phonon modes generated by impulsive femtosecond (fs) optical excitation. Broadband fs-transient absorption experiments in 1T'-MoTe2 flakes as a function of thickness (7 nm - 30 nm) yield a longitudinal sound speed of 2990 m/s. In addition, temperature dependent measurements unveil a linear decrease of the normalized Young's modulus with a slope of -0.002 per K and no noticeable change caused by the Td - 1T' structural phase transition or variations in film thickness.


title: Evolving LMXBs: CARB Magnetic Braking

abstract: The formation of low-mass X-ray binaries (LMXBs) is an ongoing challenge in stellar evolution. The important subset of LMXBs are the binary systems with a neutron star (NS) accretor. In NS LMXBs with non-degenerate donors, the mass transfer is mainly driven by magnetic braking. The discrepancies between the observed mass transfer (MT) rates and the theoretical models were known for a while. Theory predictions of the MT rates are too weak and differ by an order of magnitude or more. Recently, we showed that with the standard magnetic braking, it is not possible to find progenitor binary systems such that they could reproduce -- at any time of their evolution -- most of the observed persistent NS LMXBs. In this ${\it Letter}$ we present a modified magnetic braking prescription, CARB (Convection And Rotation Boosted). CARB magnetic braking combines two recent improvements in understanding stellar magnetic fields and magnetized winds -- the dependence of the magnetic field strength on the outer convective zone and the dependence of the Alfv\`en radius on the donor's rotation. Using this new magnetic braking prescription, we can reproduce the observed mass transfer rates at the detected mass ratio and orbital period for all well-observed to-the-date Galactic persistent NS LMXBs. For the systems where the effective temperature of the donor stars is known, theory agrees with observations as well.


title: The Assembly of the First Massive Black Holes

abstract: The existence of $\approx$10^9 Msun supermassive black holes (SMBHs) within the first billion year of the universe has stimulated numerous ideas for the prompt formation and rapid growth of BHs in the early universe. Here we review ways in which the seeds of massive BHs may have first assembled, how they may have subsequently grown as massive as $\approx$10^9 Msun, and how multi-messenger observations could distinguish between different SMBH assembly scenarios. We conclude the following: (1) The ultra-rare $\approx$10^9 Msun SMBHs represent only the tip of the iceberg. Early BHs likely fill a continuum from stellar-mass (approx. 10 Msun) to the super-massive ($\approx$10^9 Msun) regime, reflecting a range of initial masses and growth histories. (2) Stellar-mass BHs were likely left behind by the first generation of stars at redshifts as high as z=30, but their initial growth was typically stunted due to the shallow potential wells of their host galaxies. (3) Conditions in some larger, metal-poor galaxies soon became conducive to the rapid formation and growth of massive `seed' holes, via gas accretion and by mergers in dense stellar clusters. (4) BH masses depend on the environment (such as the number and properties of nearby radiation sources and the local baryonic streaming velocity), and on the metal enrichment and assembly history of the host galaxy. (5) Distinguishing between assembly mechanisms will be difficult, but a combination of observations by LISA (probing massive BH growth via mergers) and by deep multi-wavelength electromagnetic observations (probing growth via gas accretion) is particularly promising.


title: On the Expansion, Age, and Origin of the Puzzling Shell/Pulsar Wind Nebula G310.6-1.6

abstract: We present a 142-ks Chandra observation of the enigmatic combination supernova remnant G310.6-1.6 consisting of a bright pulsar-wind nebula driven by an energetic pulsar, surrounded by a highly circular, very faint shell with a featureless, probably synchrotron, spectrum. Comparison with an observation 6 years earlier shows no measurable expansion of the shell, though some features in the pulsar-wind nebula have moved. We find an expansion age of at least 2500 yr, implying a current shock velocity less than about 1000 km/s. We place severe upper limits on thermal emission from the shell; if the shell locates the blast wave, a Sedov interpretation would require the remnant to be very young, about 1000 yr, and to have resulted from a dramatically sub-energetic supernova, ejecting << 0.02 M_sun with energy E < 3 x 10^47 erg. Even a merger-induced collapse of a white dwarf to a neutron star, with a low-energy explosion, is unlikely to produce such an event. Other explanations seem equally unlikely.


title: Extended FLRW Models: dynamical cancellation of cosmological anisotropies

abstract: We investigate a corner of the Bianchi models that has not received much attention: "extended FLRW models" (eFLRW) defined as a cosmological model with underlying anisotropic Bianchi geometry that nevertheless expands isotropically and can be mapped onto a reference FLRW model with the same expansion history. In order to investigate the stability and naturalness of such models in a dynamical systems context, we consider spatially homogeneous models that contain a massless scalar field $\varphi$ and a non-tilted perfect fluid obeying an equation of state $p=w\rho$. Remarkably, we find that matter anisotropies and geometrical anisotropies tend to cancel out dynamically. Hence, the expansion is asymptotically isotropic under rather general conditions. Although extended FLRW models require a special matter sector with anisotropies that are 'fine-tuned" relative to geometrical anisotropies, our analysis shows that such solutions are dynamically preferred attractors in general relativity. Specifically, we prove that all locally rotationally symmetric Bianchi type III universes with space-like $\nabla_\mu\varphi$ are asymptotically shear-free, for all $w\in[-1,1]$. Moreover, all shear-free equilibrium sets with anisotropic spatial curvature are proved to be stable with respect to all homogeneous perturbations for $w\geq -1/3$.

categories: gr-qc astro-ph.CO hep-th math-ph math.MP


title: On the Mean Subtree Order of Graphs Under Edge Addition

abstract: For a graph $G$, the mean subtree order of $G$ is the average order of a subtree of $G$. In this note, we provide counterexamples to a recent conjecture of Chin, Gordon, MacPhee, and Vincent, that for every connected graph $G$ and every pair of distinct vertices $u$ and $v$ of $G$, the addition of the edge between $u$ and $v$ increases the mean subtree order. In fact, we show that the addition of a single edge between a pair of nonadjacent vertices in a graph of order $n$ can decrease the mean subtree order by as much as $n/3$ asymptotically. We propose the weaker conjecture that for every connected graph $G$ which is not complete, there exists a pair of nonadjacent vertices $u$ and $v$, such that the addition of the edge between $u$ and $v$ increases the mean subtree order. We prove this conjecture in the special case that $G$ is a tree.


title: Unitary Subharmonic Response and Floquet Majorana Modes

abstract: Detection and manipulation of excitations with non-Abelian statistics, such as Majorana fermions, are essential for creating topological quantum computers. To this end, we show the connection between the existence of such localized particles and the phenomenon of unitary subharmonic response (SR) in periodically driven systems. In particular, starting from highly nonequilibrium initial states, the unpaired Majorana modes exhibit spin oscillations with twice the driving period, are localized, and can have exponentially long lifetimes in clean systems. While the lifetime of SR is limited in translationally invariant systems, we show that disorder can be engineered to stabilize the subharmonic response of Majorana modes. A viable observation of this phenomenon can be achieved using modern multiqubit hardware, such as superconducting circuits and cold atomic systems.

categories: cond-mat.str-el cond-mat.stat-mech math-ph math.MP quant-ph


title: Response to NITRD, NCO, NSF Request for Information on "Update to the 2016 National Artificial Intelligence Research and Development Strategic Plan"

abstract: We present a response to the 2018 Request for Information (RFI) from the NITRD, NCO, NSF regarding the "Update to the 2016 National Artificial Intelligence Research and Development Strategic Plan." Through this document, we provide a response to the question of whether and how the National Artificial Intelligence Research and Development Strategic Plan (NAIRDSP) should be updated from the perspective of Fermilab, America's premier national laboratory for High Energy Physics (HEP). We believe the NAIRDSP should be extended in light of the rapid pace of development and innovation in the field of Artificial Intelligence (AI) since 2016, and present our recommendations below. AI has profoundly impacted many areas of human life, promising to dramatically reshape society --- e.g., economy, education, science --- in the coming years. We are still early in this process. It is critical to invest now in this technology to ensure it is safe and deployed ethically. Science and society both have a strong need for accuracy, efficiency, transparency, and accountability in algorithms, making investments in scientific AI particularly valuable. Thus far the US has been a leader in AI technologies, and we believe as a national Laboratory it is crucial to help maintain and extend this leadership. Moreover, investments in AI will be important for maintaining US leadership in the physical sciences.


title: AI-optimized detector design for the future Electron-Ion Collider: the dual-radiator RICH case

abstract: Advanced detector R&D requires performing computationally intensive and detailed simulations as part of the detector-design optimization process. We propose a general approach to this process based on Bayesian optimization and machine learning that encodes detector requirements. As a case study, we focus on the design of the dual-radiator Ring Imaging Cherenkov (dRICH) detector under development as part of the particle-identification system at the future Electron-Ion Collider (EIC). The EIC is a US-led frontier accelerator project for nuclear physics, which has been proposed to further explore the structure and interactions of nuclear matter at the scale of sea quarks and gluons. We show that the detector design obtained with our automated and highly parallelized framework outperforms the baseline dRICH design within the assumptions of the current model. Our approach can be applied to any detector R&D, provided that realistic simulations are available.