title: From the master equation to mean field game limit theory: A central limit theorem

abstract: Mean field games (MFGs) describe the limit, as $n$ tends to infinity, of stochastic differential games with $n$ players interacting with one another through their common empirical distribution. Under suitable smoothness assumptions that guarantee uniqueness of the MFG equilibrium, a form of law of large of numbers (LLN), also known as propagation of chaos, has been established to show that the MFG equilibrium arises as the limit of the sequence of empirical measures of the $n$-player game Nash equilibria, including the case when player dynamics are driven by both idiosyncratic and common sources of noise. The proof of convergence relies on the so-called master equation for the value function of the MFG, a partial differential equation on the space of probability measures. In this work, under additional assumptions, we establish a functional central limit theorem (CLT) that characterizes the limiting fluctuations around the LLN limit as the unique solution of a linear stochastic PDE. The key idea is to use the solution to the master equation to construct an associated McKean-Vlasov interacting $n$-particle system that is sufficiently close to the Nash equilibrium dynamics of the $n$-player game for large $n$. We then derive the CLT for the latter from the CLT for the former. Along the way, we obtain a new multidimensional CLT for McKean-Vlasov systems. We also illustrate the broader applicability of our methodology by applying it to establish a CLT for a specific linear-quadratic example that does not satisfy our main assumptions, and we explicitly solve the resulting stochastic PDE in this case.


title: Polynomial Heisenberg algebras, multiphoton coherent states and geometric phases

abstract: In this paper we will realize the polynomial Heisenberg algebras through the harmonic oscillator. We are going to construct then the Barut-Girardello coherent states, which coincide with the so-called multiphoton coherent states, and we will analyze the corresponding Heisenberg uncertainty relation and Wigner distribution function for some particular cases. We will show that these states are intrinsically quantum and cyclic, with a period being a fraction of the oscillator period. The associated geometric phases will be as well evaluated.


title: A Conditional Gradient Framework for Composite Convex Minimization with Applications to Semidefinite Programming

abstract: We propose a conditional gradient framework for a composite convex minimization template with broad applications. Our approach combines smoothing and homotopy techniques under the CGM framework, and provably achieves the optimal $\mathcal{O}(1/\sqrt{k})$ convergence rate. We demonstrate that the same rate holds if the linear subproblems are solved approximately with additive or multiplicative error. In contrast with the relevant work, we are able to characterize the convergence when the non-smooth term is an indicator function. Specific applications of our framework include the non-smooth minimization, semidefinite programming, and minimization with linear inclusion constraints over a compact domain. Numerical evidence demonstrates the benefits of our framework.


title: On the Impact of Fixed Point Hardware for Optical Fiber Nonlinearity Compensation Algorithms

abstract: Nonlinearity mitigation using digital signal processing has been shown to increase the achievable data rates of optical fiber transmission links. One especially effective technique is digital back propagation (DBP), an algorithm capable of simultaneously compensating for linear and nonlinear channel distortions. The most significant barrier to implementing this technique, however, is its high computational complexity. In recent years, there have been several proposed alternatives to DBP with reduced computational complexity, although such techniques have not demonstrated performance benefits commensurate with the complexity of implementation. In order to fully characterize the computational requirements of DBP, there is a need to model the algorithm behavior when constrained to the logic used in a digital coherent receiver. Such a model allows for the analysis of any signal recovery algorithm in terms of true hardware complexity which, crucially, includes the bit-depth of the multiplication operation. With a limited bit depth, there is quantization noise, introduced with each arithmetic operation, and it can no longer be assumed that the conventional DBP algorithm will outperform its low complexity alternatives. In this work, DBP and a single nonlinear step DBP implementation, the \textit{Enhanced Split Step Fourier} method (ESSFM), were compared with linear equalization using a generic software model of fixed point hardware. The requirements of bit depth and fast Fourier transform (FFT) size are discussed to examine the optimal operating regimes for these two schemes of digital nonlinearity compensation. For a 1000 km transmission system, it was found that (assuming an optimized FFT size), in terms of SNR, the ESSFM algorithm outperformed the conventional DBP for all hardware resolutions up to 13 bits.


title: Energy-, momentum-, density-, and positivity-preserving spatio-temporal discretizations for the nonlinear Landau collision operator with exact H-theorems

abstract: This paper explores energy-, momentum-, density-, and positivity-preserving spatio-temporal discretizations for the nonlinear Landau collision operator. We discuss two approaches, namely direct Galerkin formulations and discretizations of the underlying infinite-dimensional metriplectic structure of the collision integral. The spatial discretizations are chosen to reproduce the time-continuous conservation laws that correspond to Casimir invariants and to guarantee the positivity of the distribution function. Both the direct and the metriplectic discretization are demonstrated to have exact H-theorems and unique, physically exact equilibrium states. Most importantly, the two approaches are shown to coincide, given the chosen Galerkin method. A temporal discretization, preserving all of the mentioned properties, is achieved with so-called discrete gradients. Hence the proposed algorithm successfully translates all properties of the infinite-dimensional time-continuous Landau collision operator to time- and space-discrete sparse-matrix equations suitable for numerical simulation.


title: Entropy bounds for grammar compression

abstract: Grammar compression represents a string as a context free grammar. Achieving compression requires encoding such grammar as a binary string; there are a few commonly used encodings. We bound the size of practically used encodings for several heuristical compression methods, including \RePair and \Greedy algorithms: the standard encoding of \RePair, which combines entropy coding and special encoding of a grammar, achieves $1.5|S|H_k(S)$, where $H_k(S)$ is $k$-th order entropy of $S$. We also show that by stopping after some iteration we can achieve $|S|H_k(S)$. This is particularly interesting, as it explains a phenomenon observed in practice: introducing too many nonterminals causes the bit-size to grow. We generalize our approach to other compression methods like \Greedy and a wide class of irreducible grammars as well as to other practically used bit encodings (including naive, which uses fixed-length codes). Our approach not only proves the bounds but also partially explains why \Greedy and \RePair are much better in practice than other grammar based methods. In some cases we argue that our estimations are optimal. The tools used in our analysis are of independent interest: we prove the new, optimal, bounds on the zeroth order entropy of parsing of a string.


title: Eigenvector Computation and Community Detection in Asynchronous Gossip Models

abstract: We give a simple distributed algorithm for computing adjacency matrix eigenvectors for the communication graph in an asynchronous gossip model. We show how to use this algorithm to give state-of-the-art asynchronous community detection algorithms when the communication graph is drawn from the well-studied stochastic block model. Our methods also apply to a natural alternative model of randomized communication, where nodes within a community communicate more frequently than nodes in different communities. Our analysis simplifies and generalizes prior work by forging a connection between asynchronous eigenvector computation and Oja's algorithm for streaming principal component analysis. We hope that our work serves as a starting point for building further connections between the analysis of stochastic iterative methods, like Oja's algorithm, and work on asynchronous and gossip-type algorithms for distributed computation.


title: Aurigaia: mock Gaia DR2 stellar catalogues from the Auriga cosmological simulations

abstract: We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $\Lambda$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. These publicly available catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that: i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A- and B-dwarf stars up to radii of ~15 kpc; and ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to: i) develop and test analysis methods for the Gaia DR2 data; ii) gauge the limitations and biases of the data and iii) interpret the data in the light of theoretical predictions from realistic $ab$ $initio$ simulations of galaxy formation in the $\Lambda$CDM cosmological model.


title: From the master equation to mean field game limit theory: Large deviations and concentration of measure

abstract: We study a sequence of symmetric $n$-player stochastic differential games driven by both idiosyncratic and common sources of noise, in which players interact with each other through their empirical distribution. The unique Nash equilibrium empirical measure of the $n$-player game is known to converge, as $n$ goes to infinity, to the unique equilibrium of an associated mean field game. Under suitable regularity conditions, in the absence of common noise, we complement this law of large numbers result with non-asymptotic concentration bounds for the Wasserstein distance between the $n$-player Nash equilibrium empirical measure and the mean field equilibrium. We also show that the sequence of Nash equilibrium empirical measures satisfies a weak large deviation principle, which can be strengthened to a full large deviation principle only in the absence of common noise. For both sets of results, we first use the master equation, an infinite-dimensional partial differential equation that characterizes the value function of the mean field game, to construct an associated McKean-Vlasov interacting $n$-particle system that is exponentially close to the Nash equilibrium dynamics of the $n$-player game for large $n$, by refining estimates obtained in our companion paper. Then we establish a weak large deviation principle for McKean-Vlasov systems in the presence of common noise. In the absence of common noise, we upgrade this to a full large deviation principle and obtain new concentration estimates for McKean-Vlasov systems. Finally, in two specific examples that do not satisfy the assumptions of our main theorems, we show how to adapt our methodology to establish large deviations and concentration results.


title: Beam Transport Systems for Particle Therapy

abstract: The beam transport system between accelerator and patient treatment location in a particle therapy facility is described. After some general layout aspects the major beam handling tasks of this system are discussed. These are energy selection, an optimal transport of the particle beam to the beam delivery device and the gantry, a device that is able to rotate a beam delivery system around the patient, so that the tumour can be irradiated from almost any direction. Also the method of pencil beam scanning is described and how this is implemented within a gantry. Using this method the particle dose is spread over the tumour volume to the prescribed dose distribution.


title: Measurement Errors in R

abstract: This paper presents an R package to handle and represent measurements with errors in a very simple way. We briefly introduce the main concepts of metrology and propagation of uncertainty, and discuss related R packages. Building upon this, we introduce the 'errors' package, which provides a class for associating uncertainty metadata, automated propagation and reporting. Working with 'errors' enables transparent, lightweight, less error-prone handling and convenient representation of measurements with errors. Finally, we discuss the advantages, limitations and future work of computing with errors.


title: On the circular correlation coefficients for bivariate von Mises distributions on a torus

abstract: This paper studies circular correlations for the bivariate von Mises sine and cosine distributions. These are two simple and appealing models for bivariate angular data with five parameters each that have interpretations comparable to those in the ordinary bivariate normal model. However, the variability and association of the angle pairs cannot be easily deduced from the model parameters unlike the bivariate normal. Thus to compute such summary measures, tools from circular statistics are needed. We derive analytic expressions and study the properties of the Jammalamadaka-Sarma and Fisher-Lee circular correlation coefficients for the von Mises sine and cosine models. Likelihood-based inference of these coefficients from sample data is then presented. The correlation coefficients are illustrated with numerical and visual examples, and the maximum likelihood estimators are assessed on simulated and real data, with comparisons to their non-parametric counterparts. Implementations of these computations for practical use are provided in our R package BAMBI.


title: Approximate Abstractions of Markov Chains with Interval Decision Processes (Extended Version)

abstract: This work introduces a new abstraction technique for reducing the state space of large, discrete-time labelled Markov chains. The abstraction leverages the semantics of interval Markov decision processes and the existing notion of approximate probabilistic bisimulation. Whilst standard abstractions make use of abstract points that are taken from the state space of the concrete model and which serve as representatives for sets of concrete states, in this work the abstract structure is constructed considering abstract points that are not necessarily selected from the states of the concrete model, rather they are a function of these states. The resulting model presents a smaller one-step bisimulation error, when compared to a like-sized, standard Markov chain abstraction. We outline a method to perform probabilistic model checking, and show that the computational complexity of the new method is comparable to that of standard abstractions based on approximate probabilistic bisimulations.


title: Reachability and Distances under Multiple Changes

abstract: Recently it was shown that the transitive closure of a directed graph can be updated using first-order formulas after insertions and deletions of single edges in the dynamic descriptive complexity framework by Dong, Su, and Topor, and Patnaik and Immerman. In other words, Reachability is in DynFO. In this article we extend the framework to changes of multiple edges at a time, and study the Reachability and Distance queries under these changes. We show that the former problem can be maintained in DynFO$(+, \times)$ under changes affecting O($\frac{\log n}{\log \log n}$) nodes, for graphs with $n$ nodes. If the update formulas may use a majority quantifier then both Reachability and Distance can be maintained under changes that affect O($\log^c n$) nodes, for fixed $c \in \mathbb{N}$. Some preliminary results towards showing that distances are in DynFO are discussed.


title: Furstenberg systems of bounded multiplicative functions and applications

abstract: We prove a structural result for measure preserving systems naturally associated with any finite collection of multiplicative functions that take values on the complex unit disc. We show that these systems have no irrational spectrum and their building blocks are Bernoulli systems and infinite-step nilsystems. One consequence of our structural result is that strongly aperiodic multiplicative functions satisfy the logarithmically averaged variant of the disjointness conjecture of Sarnak for a wide class of zero entropy topological dynamical systems, which includes all uniquely ergodic ones. We deduce that aperiodic multiplicative functions with values plus or minus one have super-linear block growth. Another consequence of our structural result is that products of shifts of arbitrary multiplicative functions with values on the unit disc do not correlate with any totally ergodic deterministic sequence of zero mean. Our methodology is based primarily on techniques developed in a previous article of the authors where analogous results were proved for the M\"obius and the Liouville function. A new ingredient needed is a result obtained recently by Tao and Ter\"av\"ainen related to the odd order cases of the Chowla conjecture.


title: Probing hidden spin order with interpretable machine learning

abstract: The search of unconventional magnetic and nonmagnetic states is a major topic in the study of frustrated magnetism. Canonical examples of those states include various spin liquids and spin nematics. However, discerning their existence and the correct characterization is usually challenging. Here we introduce a machine-learning protocol that can identify general nematic order and their order parameter from seemingly featureless spin configurations, thus providing comprehensive insight on the presence or absence of hidden orders. We demonstrate the capabilities of our method by extracting the analytical form of nematic order parameter tensors up to rank 6. This may prove useful in the search for novel spin states and for ruling out spurious spin liquid candidates.

categories: cond-mat.stat-mech cond-mat.str-el physics.comp-ph stat.ML


title: Cosmological constraints on an exponential interaction in the dark sector

abstract: Cosmological models where dark matter (DM) and dark energy (DE) interact with each other are the general scenarios in compared to the non-interacting models. The interaction is usually motivated from the phenomenological ground and thus there is no such rule to prefer a particular interaction between DM and DE. Being motivated, in this work, allowing an exponential interaction between DM and DE in a spatially flat homogeneous and isotropic universe, we explore the dynamics of the universe through the constraints of the free parameters where the strength of the interaction is characterized by the dimensionless coupling parameter $\xi$ and the equation of state (EoS) for DE, $w_x$, is supposed to be a constant. The interaction scenario is fitted using the latest available observational data. Our analyses report that the observational data permit a non-zero value of $\xi$ but it is very small and consistent with $\xi =0$. From the constraints on $w_x$, we find that both phantom ($w_x< -1$) and quintessence ($w_x> -1$) regimes are equally allowed but $w_x$ is very close to `$-1$'. The overall results indicate that at the background level, the interaction model cannot be distinguished from the base $\Lambda $-cold dark matter model while from the perturbative analyses, the interaction model mildly deviates from the base model. We highlight that, even if we allow DM and DE to interact in an exponential manner, but according to the observational data, the evidence for a non-zero coupling is very small.


title: False Information on Web and Social Media: A Survey

abstract: False information can be created and spread easily through the web and social media platforms, resulting in widespread real-world impact. Characterizing how false information proliferates on social platforms and why it succeeds in deceiving readers are critical to develop efficient detection algorithms and tools for early detection. A recent surge of research in this area has aimed to address the key issues using methods based on feature engineering, graph mining, and information modeling. Majority of the research has primarily focused on two broad categories of false information: opinion-based (e.g., fake reviews), and fact-based (e.g., false news and hoaxes). Therefore, in this work, we present a comprehensive survey spanning diverse aspects of false information, namely (i) the actors involved in spreading false information, (ii) rationale behind successfully deceiving readers, (iii) quantifying the impact of false information, (iv) measuring its characteristics across different dimensions, and finally, (iv) algorithms developed to detect false information. In doing so, we create a unified framework to describe these recent methods and highlight a number of important directions for future research.


title: Covariant quark model for the electromagnetic structure of light nucleon resonances

abstract: We present estimates from the covariant spectator quark model for the electromagnetic transition form factors for the resonances $N(1440)\frac{1}{2}^+$, $N(1535)\frac{1}{2}^-$, $N(1520)\frac{3}{2}^-$ and $\Delta(1232)\frac{3}{2}^+$, at intermediate and large square momentum transfer ($Q^2$). The calculations associated to the $N(1440)\frac{1}{2}^+$, $N(1535)\frac{1}{2}^-$ and $N(1520)\frac{3}{2}^-$ states are based exclusively on the parametrizations derived for the nucleon. In the case of the $\Delta(1232)\frac{3}{2}^+$ (isospin 3/2), we use lattice QCD data to estimate the radial structure, and take into account the well known effects associated with the pion cloud dressing of the baryons. Our estimates are based mainly on the valence quark degrees of freedom and are in good agreement with the data for $Q^2 > 2$ GeV$^2$, with a few exceptions. The present predictions can be tested in a near future for large $Q^2$ at the Jefferson Lab -- 12 GeV upgrade.


title: The Runge Example for Interpolation and Wilkinson's Examples for Rootfinding

abstract: We look at two classical examples in the theory of numerical analysis, namely the Runge example for interpolation and Wilkinson's example (actually two examples) for rootfinding. We use the modern theory of backward error analysis and conditioning, as instigated and popularized by Wilkinson, but refined by Farouki and Rajan. By this means, we arrive at a satisfactory explanation of the puzzling phenomena encountered by students when they try to fit polynomials to numerical data, or when they try to use numerical rootfinding to find polynomial zeros. Computer algebra, with its controlled, arbitrary precision, plays an important didactic role.


title: Spatio-Temporal Neural Networks for Space-Time Series Forecasting and Relations Discovery

abstract: We introduce a dynamical spatio-temporal model formalized as a recurrent neural network for forecasting time series of spatial processes, i.e. series of observations sharing temporal and spatial dependencies. The model learns these dependencies through a structured latent dynamical component, while a decoder predicts the observations from the latent representations. We consider several variants of this model, corresponding to different prior hypothesis about the spatial relations between the series. The model is evaluated and compared to state-of-the-art baselines, on a variety of forecasting problems representative of different application areas: epidemiology, geo-spatial statistics and car-traffic prediction. Besides these evaluations, we also describe experiments showing the ability of this approach to extract relevant spatial relations.


title: A Theory of Transfers: Duality and convolution

abstract: We introduce and study the permanence properties of the class of linear transfers between probability measures. This class contains all cost minimizing mass transports, but also martingale mass transports, the Schrodinger bridge associated to a reversible Markov process, and the weak mass transports of Tala- grand, Marton, Gozlan and others. The class also includes various stochastic mass transports to which Monge-Kantorovich theory does not apply. We also introduce the cone of convex transfers, which include any p-power (p > 1) of a linear transfer, but also the logarithmic entropy, the Donsker-Varadhan infor- mation and certain free energy functionals. This first paper is mostly about exhibiting examples that point to the pervasiveness of the concept in the important work on correlating probability distributions. Duality formulae for general transfer inequalities follow in a very natural way. We also study the infinite self-convolution of a linear transfer in order to establish the existence of generalized weak KAM solutions that could be applied to the stochastic counterpart of Fathi-Mather theory.


title: Summer Research Report: Towards Incremental Lazard Cylindrical Algebraic Decomposition

abstract: Cylindrical Algebraic Decomposition (CAD) is an important tool within computational real algebraic geometry, capable of solving many problems to do with polynomial systems over the reals, but known to have worst-case computational complexity doubly exponential in the number of variables. It has long been studied by the Symbolic Computation community and is implemented in a variety of computer algebra systems, however, it has also found recent interest in the Satisfiability Checking community for use with SMT-solvers. The SCSC Project seeks to build bridges between these communities. The present report describes progress made during a Research Internship in Summer 2017 funded by the EU H2020 SCSC CSA. We describe a proof of concept implementation of an Incremental CAD algorithm in Maple, where CADs are built and refined incrementally by polynomial constraint, in contrast to the usual approach of a single computation from a single input. This advance would make CAD of use to SMT-solvers who search for solutions by constantly reformulating logical formula and querying solvers like CAD for whether a logical solution is admissible. We describe experiments for the proof of concept, which clearly display the computational advantages when compared to iterated re-computation. In addition, the project implemented this work under the recently verified Lazard projection scheme (with corresponding Lazard evaluation). That is the minimal complete CAD method in theory, and this is the first documented implementation.


title: Upgraded Photometric System of The 85-cm Telescope at Xinglong Station

abstract: The 85-cm telescope at the Xinglong station is a well-operated prime focus system with high science outputs. The telescope has been upgraded since 2014 with new corrector, filters and camera, which are provided by Beijing Normal University (BNU). The filter set is Johnson-Cousins UBVRI system. We report the test results of the new system including the bias, dark current, linearity, gain and readout noise of the CCD camera . Then we derive accurate instrumental calibration coefficients in UBVRI bands with Landolt standard stars in the photometric nights. Finally, we give the limiting magnitudes with various exposure time and signal-to-noise ratio for observers as references.


title: Lorentz violation, M\"oller scattering and finite temperature

abstract: Lorentz and CPT symmetries may be violated in new physics that emerges at very high energy scale, i.e., at the Planck scale. The differential cross section of the M\"oller scattering, due to Lorentz violation at finite temperature is calculated. Lorentz-violating effects emerge from an interaction vertex due to a CPT-odd nonminimal coupling in the covariant derivative. The finite temperature effects are determined using the Thermo Field Dynamics (TFD) formalism.


title: Gromov-Hausdorff limits of K\"ahler manifolds with Ricci curvature bounded below

abstract: We show that non-collapsed Gromov-Hausdorff limits of polarized Kahler manifolds, with Ricci curvature bounded below, are normal projective varieties, and the metric singularities of the limit space are precisely given by a countable union of analytic subvarieties. This extends a fundamental result of Donaldson-Sun, in which 2-sided Ricci curvature bounds were assumed. As a basic ingredient we show that, under lower Ricci curvature bounds, almost Euclidean balls in Kahler manifolds admit good holomorphic coordinates. Further applications are integral bounds for the scalar curvature on balls, and a rigidity theorem for Kahler manifolds with almost Euclidean volume growth.


title: Multiple Plasma Stream Electromechanical Model for Pulsed Plasma Thruster

abstract: Conventional snowplow/1-D models assume the existence of single completely ionized plasma sheet in the discharge chamber of pulsed plasma thrusters (PPTs). However, extensive experimental analysis and observations suggest the presence of multiple plasma streams and partial ionization of the propellant gas. A multi-layered plasma-neutral mass sheet one-dimensional numerical model is developed for PPTs in the present paper. The model accounts for the secondary plasma sheets generated due to the crowbar breakdowns as well as the neutral gas generated after the main discharge. The model accurately calculates the position of the plasma and neutral mass sheets, mean exit velocity, thrust delivered and the mass of individual plasma sheet developed in PPTs. The model is validated and discussed for two different ablative pulsed plasma thrusters. In the last section, the plasma sheet model and electrical circuit model are also validated for a liquid propellant PPT developed at the University of Tokyo. Results show that the multi-sheet assumption provides an accurate prediction for the discharge current and performance parameters. The multiple sheet assumption model accounts for most of the experimentally observed phenomena's which the conventional slug models tend to ignore.


title: Keys in the Clouds: Auditable Multi-device Access to Cryptographic Credentials

abstract: Personal cryptographic keys are the foundation of many secure services, but storing these keys securely is a challenge, especially if they are used from multiple devices. Storing keys in a centralized location, like an Internet-accessible server, raises serious security concerns (e.g. server compromise). Hardware-based Trusted Execution Environments (TEEs) are a well-known solution for protecting sensitive data in untrusted environments, and are now becoming available on commodity server platforms. Although the idea of protecting keys using a server-side TEE is straight-forward, in this paper we validate this approach and show that it enables new desirable functionality. We describe the design, implementation, and evaluation of a TEE-based Cloud Key Store (CKS), an online service for securely generating, storing, and using personal cryptographic keys. Using remote attestation, users receive strong assurance about the behaviour of the CKS, and can authenticate themselves using passwords while avoiding typical risks of password-based authentication like password theft or phishing. In addition, this design allows users to i) define policy-based access controls for keys; ii) delegate keys to other CKS users for a specified time and/or a limited number of uses; and iii) audit all key usages via a secure audit log. We have implemented a proof of concept CKS using Intel SGX and integrated this into GnuPG on Linux and OpenKeychain on Android. Our CKS implementation performs approximately 6,000 signature operations per second on a single desktop PC. The latency is in the same order of magnitude as using locally-stored keys, and 20x faster than smart cards.


title: Measuring Within and Between Group Inequality in Early-Life Mortality Over Time: A Bayesian Approach with Application to India

abstract: Most studies on inequality in infant and child mortality compare average mortality rates between large groups of births, for example, comparing births from different countries, income groups, ethnicities, or different times. These studies do not measure within-group disparities. The few studies that have measured within-group variability in infant and child mortality have used tools from the income inequality literature, such as Gini indices. We show that the latter are inappropriate for infant and child mortality. We develop novel tools that are appropriate for analyzing infant and child mortality inequality, including inequality measures, covariate adjustments, and ANOVA methods. We illustrate how to handle uncertainty about complex inference targets, including ensembles of probabilities and kernel density estimates. We illustrate our methodology using a large data set from India, where we estimate infant and child mortality risk for over 400,000 births using a Bayesian hierarchical model. We show that most of the variance in mortality risk exists within groups of births, not between them, and thus that within-group mortality needs to be taken into account when assessing inequality in infant and child mortality. Our approach has broad applicability to many health indicators.


title: Solving generalized Abel's integral equations of the first and second kinds via Taylor-collocation method

abstract: The aim of this paper is to present an efficient numerical procedure to approximate the generalized Abel's integral equations of the first and second kinds. For this reason, the Taylor polynomials and the collocation method are applied. Also, the error analysis of presented method is illustrated. Several examples are approximated and the numerical results show the accuracy and efficiency of this method.


title: Light-weight Head Pose Invariant Gaze Tracking

abstract: Unconstrained remote gaze tracking using off-the-shelf cameras is a challenging problem. Recently, promising algorithms for appearance-based gaze estimation using convolutional neural networks (CNN) have been proposed. Improving their robustness to various confounding factors including variable head pose, subject identity, illumination and image quality remain open problems. In this work, we study the effect of variable head pose on machine learning regressors trained to estimate gaze direction. We propose a novel branched CNN architecture that improves the robustness of gaze classifiers to variable head pose, without increasing computational cost. We also present various procedures to effectively train our gaze network including transfer learning from the more closely related task of object viewpoint estimation and from a large high-fidelity synthetic gaze dataset, which enable our ten times faster gaze network to achieve competitive accuracy to its current state-of-the-art direct competitor.


abstract: We study the interior Bernoulli free boundary problem for the infinity Laplacian. Our results cover existence, uniqueness, and characterization of solutions (above a threshold representing the "infinity Bernoulli constant"), their regularity, and their relationship with the solutions to the interior Bernoulli problem for the $p$-laplacian.


title: Isothermal Fragmentation: Is there a low-mass cut-off?

abstract: The evolution of self-gravitating clouds of isothermal gas forms the basis of many star formation theories. Therefore it is important to know under what conditions such a cloud will undergo homologous collapse into a single, massive object, or will fragment into a spectrum of smaller ones. And if it fragments, do initial conditions (e.g. Jeans mass, sonic mass) influence the mass function of the fragments, as predicted by many theories of star formation? In this paper we show that the relevant parameter separating homologous collapse from fragmentation is not the Mach number of the initial turbulence (as suspected by many), but the infall Mach number $\mathcal{M}_{\rm infall}\sim\sqrt{G M/(R c_s^2)}$, equivalent to the number of Jeans masses in the initial cloud $N_J$. We also show that fragmenting clouds produce a power-law mass function with slopes close to the expected -2 (i.e. equal mass in all logarithmic mass intervals). However, the low-mass cut-off of this mass function is entirely numerical; the initial properties of the cloud have no effect on it. In other words, if $\mathcal{M}_{\rm infall}\gg 1$, fragmentation proceeds without limit to masses much smaller than the initial Jeans mass.


title: A note on continuous-stage Runge-Kutta methods

abstract: We provide a note on continuous-stage Runge-Kutta methods (csRK) for solving initial value problems of first-order ordinary differential equations. Such methods, as an interesting and creative extension of traditional Runge-Kutta (RK) methods, can give us a new perspective on RK discretization and it may enlarge the application of RK approximation theory in modern mathematics and engineering fields. A highlighted advantage of investigation of csRK methods is that we do not need to study the tedious solution of multi-variable nonlinear algebraic equations stemming from order conditions. In this note, we will discuss and promote the recently-developed csRK theory. In particular, we will place emphasis on structure-preserving algorithms including symplectic methods, symmetric methods and energy-preserving methods which play a central role in the field of geometric numerical integration.


title: Metashooting: A Novel Tool for Free Energy Reconstruction from Polymorphic Phase Transition Mechanisms

abstract: We introduce a novel scheme for the mechanistic investigation of solid-solid phase transitions, which we dub \textit{metashooting}. Combining transition path sampling molecular dynamics and metadynamics, this scheme allows for both a complete mechanistic analysis and a detailed mapping of the free energy surface. This is illustrated by performing \textit{metashooting} calculations on the pressure-induced B4/B3 $\rightarrow$ B1 phase transition in ZnO. The resulting free energy map helps to clarify the role of intermediate configurations along this activated process and the competition between different mechanistic regimes with superior accuracy. We argue that \textit{metashooting} can be efficiently applied to a broader class of activated processes


title: Semianalytic Calculation of Gravitational Wave Spectrum Nonlinearly Induced from Primordial Curvature Perturbations

abstract: Whether or not the primordial gravitational wave (GW) produced during inflation is sufficiently strong to be observable, GWs are necessarily produced from the primordial curvature perturbations in the second order of perturbation. The induced GWs can be enhanced by curvature perturbations enhanced at small scales or by the presence of matter-dominated stages of the cosmological history. We analytically calculate the integral in the expression of the power spectrum of the induced GWs which is a universal part independent of the spectrum of the primordial curvature perturbations. This makes the subsequent numerical integrals significantly easy. In simple cases, we derive fully analytic formulas for the induced GW spectrum.


title: Sparsely-Observed Pulsating Red Giants in the AAVSO Observing Program

abstract: This paper reports on time-series analysis of 156 pulsating red giants (21 SRa, 52 SRb, 33 SR, 50 Lb) in the AAVSO (American Association of Variable Star Observers) observing program for which there are no more than 150-250 observations in total. Some results were obtained for 68 of these stars: 17 SRa, 14 SRb, 20 SR, and 17 Lb. These results generally include only an average period and amplitude. Many, if not most of the stars are undoubtedly more complex; pulsating red giants are known to have wandering periods, variable amplitudes, and often multiple periods including "long secondary periods" of unknown origin. These results (or lack thereof) raise the question of how the AAVSO should best manage the observation of these and other sparsely-observed pulsating red giants.


title: Electronic properties of ultra-thin YCrO3 films

abstract: We report on the heteroepitaxial stabilization of YCrO3 ultra-thin films on LSAT (001) substrate. Using a combination of resonant X-ray absorption spectroscopy (XAS) and atomic multiplet cluster calculation, the electronic structure of YCrO3 thin film was investigated. Polarization dependent Cr L3,2 edge XAS measurement reveal the presence of an anomalous orbital polarization uncharacteristic of a 3d3 electronic system. Atomic multiplet calculations demonstrate the critical importance of charge transfer energy, Coulomb correlation strength and hopping interaction in stabilizing this unusual orbital polarized state likely connected to the bulk multiferroicity.


title: Electron nematic fluid in a strained Sr3Ru2O7 film

abstract: Sr3Ru2O7 belongs to the family of layered strontium ruthenates and exhibits a range of unusual emergent properties, such as electron nematic behavior and metamagnetism. Here, we show that epitaxial film strain significantly modifies these phenomena. In particular, we observe enhanced magnetic interactions and an electron nematic phase that extends to much higher temperatures and over a larger magnetic field range than in bulk single crystals. Furthermore, the films show an unusual anisotropic non-Fermi liquid behavior that is controlled by the direction of the applied magnetic field. At high magnetic fields the metamagnetic transition to a ferromagnetic phase recovers isotropic Fermi-liquid behavior. The results support the interpretation that these phenomena are linked to the special features of the Fermi surface, which can be tuned by both film strain and an applied magnetic field.


title: The Radio Background Below 100 MHz

abstract: The recent detection of the "cosmic dawn" redshifted 21 cm signal at 78 MHz by the EDGES experiment differs significantly from theoretical predictions. In particular, the absorption trough is roughly a factor of two stronger than the most optimistic theoretical models. The early interpretations of the origin of this discrepancy fall into two categories. The first is that there is increased cooling of the gas due to interactions with dark matter, while the second is that the background radiation field includes a contribution from a component in addition to the cosmic microwave background. In this paper we examine the feasibility of the second idea using new data from the first station of the Long Wavelength Array. The data span 40 to 80 MHz and provide important constraints on the present-day background in a frequency range where there are few surveys with absolute temperature calibration suitable for measuring the strength of the radio monopole. We find support for a strong, diffuse radio background that was suggested by the ARCARDE 2 results in the 3 to 10 GHz range. We find that this background is well modeled by a power law with a spectral index of $-$2.58$\pm$0.05 and a temperature at the rest frame 21 cm frequency of 603$^{+102}_{-92}$ mK.


title: Perron-Frobenius Theorem for Rectangular Tensors and Directed Hypergraphs

abstract: For any positive integers $r$, $s$, $m$, $n$, an $(r,s)$-order $(n,m)$-dimensional rectangular tensor ${\cal A}=(a_{i_1\cdots i_r}^{j_1\cdots j_s}) \in ({\mathbb R}^n)^r\times ({\mathbb R}^m)^s$ is called partially symmetric if it is invariant under any permutation on the lower $r$ indexes and any permutation on the upper $s$ indexes. Such partially symmetric rectangular tensor arises naturally in studying directed hypergraphs. Ling and Qi [Front. Math. China, 2013] first studied the $(p,q)$-spectral radius (or singular values) and proved a Perron-Fronbenius theorem for such tensors when both $p,q \geq r+s$. We improved their results by extending to all $(p,q)$ satisfying $\frac{r}{p} +\frac{s}{q}\leq 1$. We also proved the Perron-Fronbenius theorem for general nonnegative $(r,s)$-order $(n,m)$-dimensional rectangular tensors when $\frac{r}{p}+\frac{s}{q}>1$. We essentially showed that this is best possible without additional conditions on $\cal A$. Finally, we applied these results to study the $(p,q)$-spectral radius of $(r,s)$-uniform directed hypergraphs.


title: Tidal Deformabilities and Radii of Neutron Stars from the Observation of GW170817

abstract: We use gravitational-wave observations of the binary neutron star merger GW170817 to explore the tidal deformabilities and radii of neutron stars. We perform Bayesian parameter estimation with the source location and distance informed by electromagnetic observations. We also assume that the two stars have the same equation of state; we demonstrate that for stars with masses comparable to the component masses of GW170817, this is effectively implemented by assuming that the stars' dimensionless tidal deformabilities are determined by the binary's mass ratio $q$ by $\Lambda_1/\Lambda_2 = q^6$. We investigate different choices of prior on the component masses of the neutron stars. We find that the tidal deformability and 90$\%$ credible interval is $\tilde{\Lambda}=222^{+420}_{-138}$ for a uniform component mass prior, $\tilde{\Lambda}=245^{+453}_{-151}$ for a component mass prior informed by radio observations of Galactic double neutron stars, and $\tilde{\Lambda}=233^{+448}_{-144}$ for a component mass prior informed by radio pulsars. We find a robust measurement of the common areal radius of the neutron stars across all mass priors of $8.9 \le \hat{R} \le 13.2$ km, with a mean value of $\langle \hat{R} \rangle = 10.8$ km. Our results are the first measurement of tidal deformability with a physical constraint on the star's equation of state and place the first lower bounds on the deformability and areal radii of neutron stars using gravitational waves.


title: Leveraging Friendship Networks for Dynamic Link Prediction in Social Interaction Networks

abstract: On-line social networks (OSNs) often contain many different types of relationships between users. When studying the structure of OSNs such as Facebook, two of the most commonly studied networks are friendship and interaction networks. The link prediction problem in friendship networks has been heavily studied. There has also been prior work on link prediction in interaction networks, independent of friendship networks. In this paper, we study the predictive power of combining friendship and interaction networks. We hypothesize that, by leveraging friendship networks, we can improve the accuracy of link prediction in interaction networks. We augment several interaction link prediction algorithms to incorporate friendships and predicted friendships. From experiments on Facebook data, we find that incorporating friendships into interaction link prediction algorithms results in higher accuracy, but incorporating predicted friendships does not when compared to incorporating current friendships.


title: On the local cartesian closure of exact completions

abstract: This paper presents a necessary and sufficient condition on a category with weak finite limits for its exact completion to be (locally) cartesian closed. A paper by Carboni and Rosolini already claimed such a characterisation using a different property on the base category, but we shall show that weak finite limits are not enough for their proof to go through. We shall also indicate how to strengthen the hypothesis for that proof to work. It will become clear that, in the case of ex/lex completions, their characterisation is still valid and it coincides with the one presented here.


title: On group rings: a translation of "Sur les anneaux de groupes" by Guy Renault

abstract: (Translation from the original French:) We characterize the rings $A$ and groups $G$ for which the group rings $A [G]$ are local, semi-local, or left perfect. The recent work of M. P. Malliavin and J. L. Pascaud permits the completion of results of [an earlier publication by the author] on self-injective group rings.


title: The random normal matrix model: insertion of a point charge

abstract: In this article, we study microscopic properties of a two-dimensional eigenvalue ensemble near a conical singularity arising from insertion of a point charge in the bulk of the support of eigenvalues. In particular, we characterize all rotationally symmetric scaling limits ('Mittag-Leffler fields') and obtain universality of them when the underlying potential is algebraic. Applications include a result on the asymptotic distribution of $\log|p_n(\zeta)|$ where $p_n$ is the characteristic polynomial of an $n$:th order random normal matrix.


title: Large Scale Scene Text Verification with Guided Attention

abstract: Many tasks are related to determining if a particular text string exists in an image. In this work, we propose a new framework that learns this task in an end-to-end way. The framework takes an image and a text string as input and then outputs the probability of the text string being present in the image. This is the first end-to-end framework that learns such relationships between text and images in scene text area. The framework does not require explicit scene text detection or recognition and thus no bounding box annotations are needed for it. It is also the first work in scene text area that tackles suh a weakly labeled problem. Based on this framework, we developed a model called Guided Attention. Our designed model achieves much better results than several state-of-the-art scene text reading based solutions for a challenging Street View Business Matching task. The task tries to find correct business names for storefront images and the dataset we collected for it is substantially larger, and more challenging than existing scene text dataset. This new real-world task provides a new perspective for studying scene text related problems. We also demonstrate the uniqueness of our task via a comparison between our problem and a typical Visual Question Answering problem.


title: Anomalous orbital structure in two-dimensional titanium dichalcogenides

abstract: Generally, lattice distortions play a key role in determining the ground states of materials. Although it is well known that trigonal distortions are generic to most two-dimensional transition metal dichalcogenides, the impact of this structural distortion on the electronic structure has not been understood conclusively. Here, by using a combination of polarization dependent X-ray absorption spectroscopy (XAS), X-ray photoelectron spectroscopy (XPS) and atomic multiplet cluster calculations, we have investigated the electronic structure of titanium dichalcogenides TiX2 (X=S, Se, Te), where the magnitude of the trigonal distortion increase monotonically from S to Se and Te. Our results reveal the presence of an anomalous and large crystal filed splitting. This unusual kind of crystal field splitting is likely responsible for the unconventional electronic structure of TiX2 compounds. Our results also indicate the drawback of the distorted crystal field picture in explaining the observed electronic ground state of these materials and emphasize the key importance of metal-ligand hybridization and electronic correlation in defining the electronic structures near Fermi energy.


title: Constraining Solar position and velocity with a Nearby Hypervelocity Star

abstract: Gravitational 3-body interaction among binary stars and the supermassive black hole (SMBH) at the center of the Milky Way occasionally ejects a hypervelocity star (HVS) with a velocity of ~1000 km/s. Due to the ejection location, such a HVS initially has negligible azimuthal angular momentum Lz ~ 0 kpc km/s. Even if the halo is mildly triaxial, Lz of a recently ejected nearby HVS remains negligible, since its flight time from the Galactic Center is too short to accumulate noticeable torque. However, if we make a wrong assumption about the Solar position and velocity, such a HVS would apparently have noticeable non-zero azimuthal angular momentum, due to the wrong reflex motion of the Sun. Conversely, with precise astrometric data for a nearby HVS, we can measure the Solar position and velocity by assuming that the HVS has zero azimuthal angular momentum. Based on this idea, here we propose a method to estimate the Galactocentric distance of the Sun R0 and the Galactocentric Solar azimuthal velocity Vsun by using a HVS. We demonstrate with mock data for a nearby HVS candidate that the Gaia astrometric data, along with the currently available constraint on Vsun/R0 from the proper motion measurement of Sgr A*, can constrain R0 and Vsun with uncertainties of ~0.27 kpc and ~7.8 km/s (or fractional uncertainties of 3 percent), respectively. Our method will be a promising tool to constrain (R0, Vsun), given that Gaia is expected to discover many nearby HVSs in the near future.


title: Hydrodynamics and two-dimensional dark lump solitons for polariton superfluids

abstract: We study a two-dimensional incoherently pumped exciton-polariton condensate described by an open-dissipative Gross-Pitaevskii equation for the polariton dynamics coupled to a rate equation for the exciton density. Adopting a hydrodynamic approach, we use multiscale expansion methods to derive several models appearing in the context of shallow water waves with viscosity. In particular, we derive a Boussinesq/Benney-Luke type equation and its far-field expansion in terms of Kadomtsev-Petviashvili-I (KP-I) equations for right- and left-going waves. From the KP-I model, we predict the existence of vorticity-free, weakly (algebraically) localized two-dimensional dark-lump solitons. We find that, in the presence of dissipation, dark lumps exhibit a lifetime three times larger than that of planar dark solitons. Direct numerical simulations show that dark lumps do exist, and their dissipative dynamics is well captured by our analytical approximation. It is also shown that lump-like and vortex-like structures can spontaneously be formed as a result of the transverse "snaking" instability of dark soliton stripes.


title: Stationary-Complete Spacetimes with non-standard splittings and pre-Randers metrics

abstract: Using the relativistic Fermat's principle, we establish a bridge between stationary-complete manifolds which satisfy the observer-manifold condition and pre-Randers metrics, namely, Randers metrics without any restriction on the one-form. As a consequence, we give a description of the causal ladder of such spacetimes in terms of the elements associated with the pre-Randers metric: its geodesics and the associated distance. We obtain, as applications of this interplay, the description of conformal maps of Killing submersions, and existence and multiplicity results for geodesics of pre-Randers metrics and magnetic geodesics.


title: Solving DWF Dirac Equation Using Multisplitting Preconditioned Conjugate Gradient

abstract: We show that using the multisplitting algorithm as a preconditioner for conjugate gradient inversion of the domain wall Dirac operator could effectively reduce the inter- node communication cost, at the expense of performing more on-node floating point operations. This method could be useful for supercomputers with far more on-node flops than inter-node communication bandwidth.


title: Transverse instability of electron phase-space holes in multi-dimensional Maxwellian plasmas

abstract: The stability of an initially one-dimensional electron hole to perturbations varying sinusoidally transverse to its trapping direction is analysed in detail. It is shown that the expected low-frequency eigenmode of the linearized Vlasov-Poisson system consists of a shift-mode, proportional to the gradient of the equilibrium potential. The resulting dispersion relation is that the total jetting force exerted by a perturbed hole on the particles balances the electric restoring tension of the hole. The tension is quantitatively small and can often be ignored. The particle force is expressed as integrals of equilibrium parameters over the hole and is shown at low frequency to be exactly equal to what has recently been found (by different analysis) to express `kinematic' hole momentum conservation. The mechanism of instability has nothing to do with the previously hypothesized transverse electron focusing. The unmagnetized growth rate $\gamma(k)$ is found numerically and is in excellent agreement with recent kinematic estimates. Magnetic field stabilization of the transverse mode is also evaluated. The resulting stability boundary for Maxwellian holes is in reasonable agreement with previously published criteria based on particle simulation. It arises from a change of trapped force sign across the resonance between bounce and cyclotron frequencies.


title: Single-User mmWave Massive MIMO: SVD-based ADC Bit Allocation and Combiner Design

abstract: In this paper, we propose a Singular-Value-Decomposition-based variable-resolution Analog to Digital Converter (ADC) bit allocation design for a single-user Millimeter wave massive Multiple-Input Multiple-Output receiver. We derive the optimality condition for bit allocation under a power constraint. This condition ensures optimal receiver performance in the Mean Squared Error (MSE) sense. We derive the MSE expression and show that it approaches the Cramer-Rao Lower Bound (CRLB). The CRLB is seen to be a function of the analog combiner, the digital combiner, and the bit allocation matrix. We attempt to minimize the CRLB with respect to the bit allocation matrix by making suitable assumptions regarding the structure of the combiners. In doing so, the bit allocation design reduces to a set of simple inequalities consisting of ADC bits, channel singular values and covariance of the quantization noise along each RF path. This results in a simple and computationally efficient bit allocation algorithm. Using simulations, we show that the MSE performance of our proposed bit allocation is very close to that of the Full Search (FS) bit allocation. We also show that the computational complexity of our proposed method has an order of magnitude improvement compared to FS and Genetic Algorithm based bit allocation of $\cite{Zakir1}$


title: Geometric pinning and antimixing in scaffolded lipid vesicles

abstract: Unravelling the physical mechanisms behind the organisation of lipid domains is a central goal in cell biology and membrane biophysics. Previous studies on cells and model lipid bilayers featuring phase-separated domains found an intricate interplay between the membrane geometry and its chemical composition. However, the lack of a model system with simultaneous control over the membrane shape and conservation of its composition precluded a fundamental understanding of curvature-induced effects. Here, we present a new class of multicomponent vesicles supported by colloidal scaffolds of designed shape. We find that the domain composition adapts to the geometry, giving rise to a novel "antimixed" state. Theoretical modelling allowed us to link the pinning of domains by regions of high curvature to the material parameters of the membrane. Our results provide key insights into the phase separation of cellular membranes and on curved surfaces in general.


title: Towards Symbolic Reinforcement Learning with Common Sense

abstract: Deep Reinforcement Learning (deep RL) has made several breakthroughs in recent years in applications ranging from complex control tasks in unmanned vehicles to game playing. Despite their success, deep RL still lacks several important capacities of human intelligence, such as transfer learning, abstraction and interpretability. Deep Symbolic Reinforcement Learning (DSRL) seeks to incorporate such capacities to deep Q-networks (DQN) by learning a relevant symbolic representation prior to using Q-learning. In this paper, we propose a novel extension of DSRL, which we call Symbolic Reinforcement Learning with Common Sense (SRL+CS), offering a better balance between generalization and specialization, inspired by principles of common sense when assigning rewards and aggregating Q-values. Experiments reported in this paper show that SRL+CS learns consistently faster than Q-learning and DSRL, achieving also a higher accuracy. In the hardest case, where agents were trained in a deterministic environment and tested in a random environment, SRL+CS achieves nearly 100% average accuracy compared to DSRL's 70% and DQN's 50% accuracy. To the best of our knowledge, this is the first case of near perfect zero-shot transfer learning using Reinforcement Learning.


title: Black-box Adversarial Attacks with Limited Queries and Information

abstract: Current neural network-based classifiers are susceptible to adversarial examples even in the black-box setting, where the attacker only has query access to the model. In practice, the threat model for real-world systems is often more restrictive than the typical black-box model where the adversary can observe the full output of the network on arbitrarily many chosen inputs. We define three realistic threat models that more accurately characterize many real-world classifiers: the query-limited setting, the partial-information setting, and the label-only setting. We develop new attacks that fool classifiers under these more restrictive threat models, where previous methods would be impractical or ineffective. We demonstrate that our methods are effective against an ImageNet classifier under our proposed threat models. We also demonstrate a targeted black-box attack against a commercial classifier, overcoming the challenges of limited query access, partial information, and other practical issues to break the Google Cloud Vision API.


title: On capacities of the two-user union channel with complete feedback

abstract: The exact values of the optimal symmetric rate point in the Cover--Leung capacity region of the two-user union channel with complete feedback were determined by Willems when the size of the input alphabet is 2, and by Vinck, Hoeks and Post when the size is at least 6. We complete this line of research when the size of the input alphabet is 3, 4 or 5. The proof hinges on the technical lemma that concerns the maximal joint entropy of two independent random variables in terms of their probability of equality. For the zero-error capacity region, using superposition coding, we provide a practical near-optimal communication scheme which improves all the previous explicit constructions.


title: Revisiting nonrelativistic limits of field theory: Antiparticles

abstract: Different subtleties and problems associated with a nonrelativistic limit of the field theory to the Schroedinger theory are discussed. In this paper, we revisit different cases of the nonrelativistic limit of a real and complex scalar field in the level of the Lagrangian and the equation of motion. We develop the nonrelativistic limit of the Dirac equation and action in the way that the nonrelativistic limit of spin-$\frac{1}{2}$ wave functions of particles and antiparticles appear simultaneously. We study the effect of a potential like $U(\phi)\propto \phi^4$ which can be attributed to axion dark matter field in this limit. We develop a formalism for studying the nonrelativistic limit of antiparticles in the quantum mechanics. We discussed the non-local approach for the nonrelativistic limit and its problems.


title: CMB Spectral Distortions from Cooling Macroscopic Dark Matter

abstract: We propose a new mechanism by which dark matter (DM) can affect the early universe. The hot interior of a macroscopic DM, or macro, can behave as a heat reservoir so that energetic photons are emitted from its surface. This results in spectral distortions (SDs) of the cosmic microwave background. The SDs depend on the density and the cooling processes of the interior, and the surface composition of the Macros. We use neutron stars as a model for nuclear-density Macros and find that the spectral distortions are mass-independent for fixed density. In our work, we find that, for Macros of this type that constitute 100$\%$ of the dark matter, the $\mu$ and $y$ distortions can be above detection threshold for typical proposed next-generation experiments such as PIXIE.


title: Imaging the zigzag Wigner crystal in confinement-tunable quantum wires

abstract: The existence of Wigner crystallization, one of the most significant hallmarks of strong electron correlations, has to date only been definitively observed in two-dimensional systems. In one-dimensional (1D) quantum wires Wigner crystals correspond to regularly spaced electrons; however, weakening the confinement and allowing the electrons to relax in a second dimension is predicted to lead to the formation of a new ground state constituting a zigzag chain with nontrivial spin phases and properties. Here we report the observation of such zigzag Wigner crystals by use of on-chip charge and spin detectors employing electron focusing to image the charge density distribution and probe their spin properties. This experiment demonstrates both the structural and spin phase diagrams of the 1D Wigner crystallization. The existence of zigzag spin chains and phases which can be electrically controlled in semiconductor systems may open avenues for experimental studies of Wigner crystals and their technological applications in spintronics and quantum information.


title: Towards Learning Sparsely Used Dictionaries with Arbitrary Supports

abstract: Dictionary learning is a popular approach for inferring a hidden basis or dictionary in which data has a sparse representation. Data generated from the dictionary A (an n by m matrix, with m > n in the over-complete setting) is given by Y = AX where X is a matrix whose columns have supports chosen from a distribution over k-sparse vectors, and the non-zero values chosen from a symmetric distribution. Given Y, the goal is to recover A and X in polynomial time. Existing algorithms give polytime guarantees for recovering incoherent dictionaries, under strong distributional assumptions both on the supports of the columns of X, and on the values of the non-zero entries. In this work, we study the following question: Can we design efficient algorithms for recovering dictionaries when the supports of the columns of X are arbitrary? To address this question while circumventing the issue of non-identifiability, we study a natural semirandom model for dictionary learning where there are a large number of samples $y=Ax$ with arbitrary k-sparse supports for x, along with a few samples where the sparse supports are chosen uniformly at random. While the few samples with random supports ensures identifiability, the support distribution can look almost arbitrary in aggregate. Hence existing algorithmic techniques seem to break down as they make strong assumptions on the supports. Our main contribution is a new polynomial time algorithm for learning incoherent over-complete dictionaries that works under the semirandom model. Additionally the same algorithm provides polynomial time guarantees in new parameter regimes when the supports are fully random. Finally using these techniques, we also identify a minimal set of conditions on the supports under which the dictionary can be (information theoretically) recovered from polynomial samples for almost linear sparsity, i.e., $k=\tilde{O}(n)$.


title: The discrete twofold Ellis-Gohberg inverse problem

abstract: In this paper a twofold inverse problem for orthogonal matrix functions in the Wiener class is considered. The scalar-valued version of this problem was solved by Ellis and Gohberg in 1992. Under reasonable conditions, the problem is reduced to an invertibility condition on an operator that is defined using the Hankel and Toeplitz operators associated to the Wiener class functions that comprise the data set of the inverse problem. It is also shown that in this case the solution is unique. Special attention is given to the case that the Hankel operator of the solution is a strict contraction and the case where the functions are matrix polynomials.


title: Frustrated Magnetism in Mott Insulating (V$_{1-x}$Cr$_x$)$_2$O$_3$

abstract: V2O3 famously features all four combinations of paramagnetic vs antiferromagnetic, and metallic vs insulating states of matter in response to %-level doping, pressure in the GPa range, and temperature below 300 K. Using time-of-flight neutron spectroscopy combined with density functional theory calculations of magnetic interactions, we have mapped and analyzed the inelastic magnetic neutron scattering cross section over a wide range of energy and momentum transfer in the chromium stabilized antiferromagnetic and paramagnetic insulating phases (AFI & PI). Our results reveal an important magnetic frustration and degeneracy of the PI phase which is relieved by the rhombohedral to monoclinic transition at $T_N=185$ K due to a significant magneto-elastic coupling. This leads to the recognition that magnetic frustration is an inherent property of the paramagnetic phase in $\rm (V_{1-x}Cr_x)_2O_3$ and plays a key role in suppressing the magnetic long range ordering temperature and exposing a large phase space for the paramagnetic Mott metal-insulator transition to occur.


title: Zero-Shot Visual Imitation

abstract: The current dominant paradigm for imitation learning relies on strong supervision of expert actions to learn both 'what' and 'how' to imitate. We pursue an alternative paradigm wherein an agent first explores the world without any expert supervision and then distills its experience into a goal-conditioned skill policy with a novel forward consistency loss. In our framework, the role of the expert is only to communicate the goals (i.e., what to imitate) during inference. The learned policy is then employed to mimic the expert (i.e., how to imitate) after seeing just a sequence of images demonstrating the desired task. Our method is 'zero-shot' in the sense that the agent never has access to expert actions during training or for the task demonstration at inference. We evaluate our zero-shot imitator in two real-world settings: complex rope manipulation with a Baxter robot and navigation in previously unseen office environments with a TurtleBot. Through further experiments in VizDoom simulation, we provide evidence that better mechanisms for exploration lead to learning a more capable policy which in turn improves end task performance. Videos, models, and more details are available at https://pathak22.github.io/zeroshot-imitation/


title: Benchmarking projective simulation in navigation problems

abstract: Projective simulation (PS) is a model for intelligent agents with a deliberation capacity that is based on episodic memory. The model has been shown to provide a flexible framework for constructing reinforcement-learning agents, and it allows for quantum mechanical generalization, which leads to a speed-up in deliberation time. PS agents have been applied successfully in the context of complex skill learning in robotics, and in the design of state-of-the-art quantum experiments. In this paper, we study the performance of projective simulation in two benchmarking problems in navigation, namely the grid world and the mountain car problem. The performance of PS is compared to standard tabular reinforcement learning approaches, Q-learning and SARSA. Our comparison demonstrates that the performance of PS and standard learning approaches are qualitatively and quantitatively similar, while it is much easier to choose optimal model parameters in case of projective simulation, with a reduced computational effort of one to two orders of magnitude. Our results show that the projective simulation model stands out for its simplicity in terms of the number of model parameters, which makes it simple to set up the learning agent in unknown task environments.


title: Constrained optimal design of automotive radar arrays using the Weiss-Weinstein Bound

abstract: We propose a design strategy for optimizing antenna positions in linear arrays for far-field Direction of Arrival (DoA) estimation of narrow-band sources in collocated MIMO radar. Our methodology allows to consider any spatial constraints and number of antennas, using as optimization function the Weiss- Weinstein bound formulated for an observation model with random target phase and known SNR, over a pre-determined Field-of-View (FoV). Optimized arrays are calculated for the typical case of a 77GHz MIMO radar of 3Tx and 4Rx channels. Simulations demonstrate a performance improvement of the proposed arrays compared to the corresponding uniform and minimum redundancy arrays for a wide regime of SNR values.


title: A data-driven framework for sparsity-enhanced surrogates with arbitrary mutually dependent randomness

abstract: The challenge of quantifying uncertainty propagation in real-world systems is rooted in the high-dimensionality of the stochastic input and the frequent lack of explicit knowledge of its probability distribution. Traditional approaches show limitations for such problems. To address these difficulties, we have developed a general framework of constructing surrogate models on spaces of stochastic input with arbitrary probability measure irrespective of the mutual dependencies between individual components and the analytical form. The present Data-driven Sparsity-enhancing Rotation for Arbitrary Randomness (DSRAR) framework includes a data-driven construction of multivariate polynomial basis for arbitrary mutually dependent probability measure and a sparsity enhancement rotation procedure. This sparsity-enhancing rotation method was initially proposed in our previous work [1] for Gaussian distributions, which may not be feasible for non-Gaussian distributions due to the loss of orthogonality after the rotation. To remedy such difficulties, we developed the new approach to construct orthonormal polynomials for arbitrary mutually dependent (amdP) randomness, ensuring the constructed basis maintains the orthogonality with respect to the density of the rotated random vector, where directly applying the regular polynomial chaos including arbitrary polynomial chaos (aPC) [2] shows limitations due to the assumption of the mutual independence between the components of the random inputs. The developed DSRAR framework leads to accurate recovery of a sparse representation of the target functions. The effectiveness of our method is demonstrated in challenging problems such as PDEs and realistic molecular systems where the underlying density is implicitly represented by a large collection of sample data, as well as systems with explicitly given non-Gaussian probabilistic measures.


title: Room temperature fiber laser at 3.92 {\mu}m

abstract: Rare-earth-doped fiber lasers are promising contenders in the development of spectroscopy, free-space communications and countermeasure applications in the 3 - 5 {\mu}m spectral region. However, given the limited transparency of the commonly used fluorozirconate glass fiber, these systems have only achieved wavelength coverage up to 3.8 {\mu}m, hence fueling the development of more suitable fiber glass compositions. To this extent, we propose in this Letter a novel heavily holmium-doped fluoroindate fiber, providing extended transparency up to 5 {\mu}m, to demonstrate the longest wavelength room temperature fiber laser at 3.92 {\mu}m. Achieving $\approx$ 200 mW of output power when cladding pumped by a commercial 888 nm laser diode, this demonstration paves the way for powerful mid-infrared fiber lasers emitting at and beyond 4 {\mu}m.


title: Faster Response in Bounded-Update-Rate, Discrete-time Networks using Delayed Self-Reinforcement

abstract: The response speed of a network impacts the efficacy of its actions to external stimuli. However, for a given bound on the update rate, the network-response speed is limited by the need to maintain stability. This work increases the network-response speed without having to increase the update rate by using delayed self-reinforcement (DSR), where each agent uses its already available information from the network to strengthen its individual update law. Example simulation results are presented that show more than an order of magnitude improvement in the response speed (quantified using the settling time) with the proposed DSR approach.


title: A symmetric formula for hypergeometric series

abstract: In terms of Dougall's $_2H_2$ series identity and the series rearrangement method, we establish an interesting symmetric formula for hypergeometric series. Then it is utilized to derive a known nonterminating form of Saalsch\"{u}tz's theorem. Similarly, we also show that Bailey's $_6\psi_6$ series identity implies the nonterminating form of Jackson's $_8\phi_7$ summation formula. Considering the reversibility of the proofs, it is routine to show that Dougall's $_2H_2$ series identity is equivalent to a known nonterminating form of Saalsch\"{u}tz's theorem and Bailey's $_6\psi_6$ series identity is equivalent to the nonterminating form of Jackson's $_8\phi_7$ summation formula.


title: Parameter Transfer Unit for Deep Neural Networks

abstract: Parameters in deep neural networks which are trained on large-scale databases can generalize across multiple domains, which is referred as "transferability". Unfortunately, the transferability is usually defined as discrete states and it differs with domains and network architectures. Existing works usually heuristically apply parameter-sharing or fine-tuning, and there is no principled approach to learn a parameter transfer strategy. To address the gap, a parameter transfer unit (PTU) is proposed in this paper. The PTU learns a fine-grained nonlinear combination of activations from both the source and the target domain networks, and subsumes hand-crafted discrete transfer states. In the PTU, the transferability is controlled by two gates which are artificial neurons and can be learned from data. The PTU is a general and flexible module which can be used in both CNNs and RNNs. Experiments are conducted with various network architectures and multiple transfer domain pairs. Results demonstrate the effectiveness of the PTU as it outperforms heuristic parameter-sharing and fine-tuning in most settings.


title: Sunspot Drawings by Japanese Official Astronomers in 1749-1750

abstract: Sunspot observations with telescopes in 18th century were carried out in Japan as well. One of these sunspot observations is recorded in an account called Sansaizusetsu narabini Kansei irai Jissoku Zusetsu (Charts of Three Worlds and Diagrams of Actual Observations since Kansei Era). We analyze manuscripts of this account to show in total 15 sunspot drawings in 1749-1750. These observations were carried out by contemporary official astronomers in Japan, with telescopes covered by zongurasus (< zonglas in Dutch, corresponding to "sunglass" in English). We count their group number of sunspots to locate them in long-term solar activity and show that their observations were situated around the solar maximum in 1749 or 1750. We also computed their locations and areas, while we have to admit the difference of variant manuscripts with one another. These observational records show the spread of sunspot observations not only in Europe but also in Japan and hence may contribute to crosscheck or possibly improve the known sunspot indices.


title: QSAR Classification Modeling for Bioactivity of Molecular Structure via SPL-Logsum

abstract: Quantitative structure-activity relationship (QSAR) modelling is effective 'bridge' to search the reliable relationship related bioactivity to molecular structure. A QSAR classification model contains a lager number of redundant, noisy and irrelevant descriptors. To address this problem, various of methods have been proposed for descriptor selection. Generally, they can be grouped into three categories: filters, wrappers, and embedded methods. Regularization method is an important embedded technology, which can be used for continuous shrinkage and automatic descriptors selection. In recent years, the interest of researchers in the application of regularization techniques is increasing in descriptors selection , such as, logistic regression(LR) with $L_1$ penalty. In this paper, we proposed a novel descriptor selection method based on self-paced learning(SPL) with Logsum penalized LR for predicting the bioactivity of molecular structure. SPL inspired by the learning process of humans and animals that gradually learns from easy samples(smaller losses) to hard samples(bigger losses) samples into training and Logsum regularization has capacity to select few meaningful and significant molecular descriptors, respectively. Experimental results on simulation and three public QSAR datasets show that our proposed SPL-Logsum method outperforms other commonly used sparse methods in terms of classification performance and model interpretation.


title: Solid deuterium surface degradation at ultracold neutron sources

abstract: Solid deuterium (sD_2) is used as an efficient converter to produce ultracold neutrons (UCN). It is known that the sD_2 must be sufficiently cold, of high purity and mostly in its ortho-state in order to guarantee long lifetimes of UCN in the solid from which they are extracted into vacuum. Also the UCN transparency of the bulk sD_2 material must be high because crystal inhomogeneities limit the mean free path for elastic scattering and reduce the extraction efficiency. Observations at the UCN sources at Paul Scherrer Institute and at Los Alamos National Laboratory consistently show a decrease of the UCN yield with time of operation after initial preparation or later treatment (`conditioning') of the sD_2. We show that, in addition to the quality of the bulk sD_2, the quality of its surface is essential. Our observations and simulations support the view that the surface is deteriorating due to a build-up of D_2 frost-layers under pulsed operation which leads to strong albedo reflections of UCN and subsequent loss. We report results of UCN yield measurements, temperature and pressure behavior of deuterium during source operation and conditioning, and UCN transport simulations. This, together with optical observations of sD_2 frost formation on initially transparent sD_2 in offline studies with pulsed heat input at the North Carolina State University UCN source results in a consistent description of the UCN yield decrease.


title: Distributed Distributional Deterministic Policy Gradients

abstract: This work adopts the very successful distributional perspective on reinforcement learning and adapts it to the continuous control setting. We combine this within a distributed framework for off-policy learning in order to develop what we call the Distributed Distributional Deep Deterministic Policy Gradient algorithm, D4PG. We also combine this technique with a number of additional, simple improvements such as the use of $N$-step returns and prioritized experience replay. Experimentally we examine the contribution of each of these individual components, and show how they interact, as well as their combined contributions. Our results show that across a wide variety of simple control tasks, difficult manipulation tasks, and a set of hard obstacle-based locomotion tasks the D4PG algorithm achieves state of the art performance.


title: MIKiS: the Multi-Instrument Kinematic Survey of Galactic Globular Clusters. I. Velocity dispersion profiles and rotation signals of 11 globular clusters

abstract: We present the first results of the Multi-Instrument Kinematic Survey of Galactic Globular Clusters, a project aimed at exploring the internal kinematics of a representative sample of Galactic globular clusters from the radial velocity of individual stars, covering the entire radial extension of each system. This is achieved by exploiting the formidable combination of multi-object and integral field unit spectroscopic facilities of the ESO Very Large Telescope. As a first step, here we discuss the results obtained for 11 clusters from high and medium resolution spectra acquired through a combination of FLAMES and KMOS observations. We provide the first kinematical characterization of NGC 1261 and NGC 6496. In all the surveyed systems, the velocity dispersion profile declines at increasing radii, in agreement with the expectation from the King model that best fits the density/luminosity profile. In the majority of the surveyed systems we find evidence of rotation within a few half-mass radii from the center. These results are in general overall agreement with the predictions of recent theoretical studies, suggesting that the detected signals could be the relic of significant internal rotation set at the epoch of the cluster's formation.


title: State Distribution-aware Sampling for Deep Q-learning

abstract: A critical and challenging problem in reinforcement learning is how to learn the state-action value function from the experience replay buffer and simultaneously keep sample efficiency and faster convergence to a high quality solution. In prior works, transitions are uniformly sampled at random from the replay buffer or sampled based on their priority measured by temporal-difference (TD) error. However, these approaches do not fully take into consideration the intrinsic characteristics of transition distribution in the state space and could result in redundant and unnecessary TD updates, slowing down the convergence of the learning procedure. To overcome this problem, we propose a novel state distribution-aware sampling method to balance the replay times for transitions with skew distribution, which takes into account both the occurrence frequencies of transitions and the uncertainty of state-action values. Consequently, our approach could reduce the unnecessary TD updates and increase the TD updates for state-action value with more uncertainty, making the experience replay more effective and efficient. Extensive experiments are conducted on both classic control tasks and Atari 2600 games based on OpenAI gym platform and the experimental results demonstrate the effectiveness of our approach in comparison with the standard DQN approach.


title: Progressive instability in circular masonry columns

abstract: The instability behaviour of eccentrically loaded circular masonry columns is investigated. Two approaches are considered for the analysis. One is based on a semi-analytical formulation of the relevant boundary-value problem for a no-tension material response; the other employs a plastic-damage-contact constitutive model, the CraftS model, to capture the complex microstructural behaviour of the material. The latter has been implemented in the finite element program LUSAS and has been already successfully employed to describe progressive instability in eccentrically loaded brickwork wallettes of rectangular cross section. Equilibrium paths and limit load estimations are computed for both analysis approaches for a range of column aspect ratios and load eccentricities. It is shown that the type of material response becomes less important for specimens with height-to-diameter aspect ratios greater than 7.5 and for loads applied to points in the kernel of the cross section, while for higher eccentricities the presence of a tensile strength increases considerably the limit load. The damage evolution predicted by the models is also investigated for selected cases, showing that the formulation based on the no-tension material is able to capture with good agreement the damaged zone of the column for loads with low eccentricities. For the same type of loading, a useful design formula is provided.


title: Integrability and non-Integrability in N=2 SCFTs and their Holographic Backgrounds

abstract: We show that the string worldsheet theory of Gaiotto-Maldacena holographic duals to N=2 superconformal field theories generically fails to be classically integrable. We demonstrate numerically that the dynamics of a winding string configuration possesses a non-vanishing Lyapunov exponent. Furthermore an analytic study of the Normal Variational Equation fails to yield a Liouvillian solution. An exception to the generic non-integrability of such backgrounds is provided by the non-Abelian T-dual of $AdS_5 \times S^5$. Here by virtue of the canonical transformation nature of the T-duality classical integrability is known to be present.


title: The Origin of the 300 km s$^{-1}$ Stream Near Segue 1

abstract: We present a search for new members of the 300 km s$^{-1}$ stream (300S) near the dwarf galaxy Segue 1 using wide-field survey data. We identify 11 previously unknown bright stream members in the APOGEE-2 and SEGUE-1 and 2 spectroscopic surveys. Based on the spatial distribution of the high-velocity stars, we confirm for the first time that this kinematic structure is associated with a 24$^{\circ}$-long stream seen in SDSS and Pan-STARRS imaging data. The 300S stars display a metallicity range of $-2.17 < {\rm [Fe/H]} < -1.24$, with an intrinsic dispersion of 0.21$_{-0.09}^{+0.12}$ dex. They also have chemical abundance patterns similar to those of Local Group dwarf galaxies, as well as that of the Milky Way halo. Using the open-source code galpy to model the orbit of the stream, we find that the progenitor of the stream passed perigalacticon about 70 Myr ago, with a closest approach to the Galactic Center of about 4.1 kpc. Using Pan-STARRS DR1 data, we obtain an integrated stream luminosity of $4 \times 10^3$ L$_{\odot}$. We conclude that the progenitor of the stream was a dwarf galaxy that is probably similar to the satellites that were accreted to build the present-day Milky Way halo.


title: De Sitter Holography and Entanglement Entropy

abstract: We propose a new example of entanglement knitting spacetime together, satisfying a series of checks of the corresponding von Neumann and Renyi entropies. The conjectured dual of de Sitter in d+1 dimensions involves two coupled CFT sectors constrained by residual d-dimensional gravity. In the d=2 case, the gravitational constraints and the CFT spectrum are relatively tractable. We identify a finite portion of each CFT Hilbert space relevant for de Sitter. Its maximum energy level coincides with the transition to the universal Cardy behavior for theories with a large central charge and a sparse light spectrum, derived by Hartman, Keller, and Stoica. Significant interactions between the two CFTs, derived previously for other reasons, suggest a maximally mixed state upon tracing out one of the two sectors; we derive this by determining the holographic Renyi entropies. The resulting entanglement entropy matches the Gibbons-Hawking formula for de Sitter entropy, including the numerical coefficient. Finally, we interpret the Gibbons-Hawking horizon entropy in terms of the Ryu-Takayanagi entropy, and explore the time evolution of the entanglement entropy.


title: Quasars probing quasars X: The quasar pair spectral database

abstract: The rare close projection of two quasars on the sky provides the opportunity to study the host galaxy environment of a foreground quasar in absorption against the continuum emission of a background quasar. For over a decade the "Quasars probing quasars" series has utilized this technique to further the understanding of galaxy formation and evolution in the presence of a quasar at z>2, resolving scales as small as a galactic disc and from bound gas in the circumgalactic medium to the diffuse environs of intergalactic space. Presented here, is the public release of the quasar pair spectral database utilized in these studies. In addition to projected pairs at z>2, the database also includes quasar pair members at z<2, gravitational lens candidates and quasars closely separated in redshift that are useful for small-scale clustering studies. In total the database catalogs 5627 distinct objects, with 4083 lying within 5' of at least one other source. A spectral library contains 3582 optical and near-infrared spectra for 3028 of the cataloged sources. As well as reporting on 54 newly discovered quasar pairs, we outline the key contributions made by this series over the last ten years, summarize the imaging and spectroscopic data used for target selection, discuss the target selection methodologies, describe the database content and explore some avenues for future work. Full documentation for the spectral database, including download instructions, are supplied at the following URL http://specdb.readthedocs.io/en/latest/


title: Gravitational free energy in topological AdS/CFT

abstract: We define and study a holographic dual to the topological twist of $\mathcal{N}=4$ gauge theories on Riemannian three-manifolds. The gravity duals are solutions to four-dimensional $\mathcal{N}=4$ gauged supergravity, where the three-manifold arises as a conformal boundary. Following our previous work, we show that the renormalized gravitational free energy of such solutions is independent of the boundary three-metric, as required for a topological theory. We then go further, analyzing the geometry of supersymmetric bulk solutions. Remarkably, we are able to show that the gravitational free energy of any smooth four-manifold filling of any three-manifold is always zero. Aided by this analysis, we prove a similar result for topological AdS$_5$/CFT$_4$. We comment on the implications of these results for the large $N$ limits of topologically twisted gauge theories in three and four dimensions, including the ABJM theory and $\mathcal{N}=4$ $SU(N)$ super-Yang-Mills, respectively.


title: Bound on Eigenstate Thermalization from Transport

abstract: We show that macroscopic thermalization and transport impose constraints on matrix elements entering the Eigenstate Thermalization Hypothesis (ETH) ansatz and require them to be correlated. It is often assumed that the ETH reduces to Random Matrix Theory (RMT) below the Thouless energy scale. We show this conventional picture is not self-consistent. We prove that energy scale at which the RMT behavior emerges has to be parametrically smaller than the inverse timescale of the slowest thermalization mode coupled to the operator of interest. We argue that the timescale marking the onset of the RMT behavior is the same timescale at which hydrodynamic description of transport breaks down.

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


title: A Deeper Look at the New Milky Way Satellites: Sagittarius II, Reticulum II, Phoenix II, and Tucana III

abstract: We present deep Magellan/Megacam stellar photometry of four recently discovered faint Milky Way satellites: Sagittarius II (Sgr II), Reticulum II (Ret II), Phoenix II (Phe II), and Tucana III (Tuc III). Our photometry reaches ~2-3 magnitudes deeper than the discovery data, allowing us to revisit the properties of these new objects (e.g., distance, structural properties, luminosity measurements, and signs of tidal disturbance). The satellite color-magnitude diagrams show that they are all old (~13.5 Gyr) and metal-poor ([Fe/H]$\lesssim-2.2$). Sgr II is particularly interesting as it sits in an intermediate position between the loci of dwarf galaxies and globular clusters in the size-luminosity plane. The ensemble of its structural parameters is more consistent with a globular cluster classification, indicating that Sgr II is the most extended globular cluster in its luminosity range. The other three satellites land directly on the locus defined by Milky Way ultra-faint dwarf galaxies of similar luminosity. Ret II is the most elongated nearby dwarf galaxy currently known for its luminosity range. Our structural parameters for Phe II and Tuc III suggest that they are both dwarf galaxies. Tuc III is known to be associated with a stellar stream, which is clearly visible in our matched-filter stellar density map. The other satellites do not show any clear evidence of tidal stripping in the form of extensions or distortions. Finally, we also use archival HI data to place limits on the gas content of each object.


title: Surface Topological Order and a new 't Hooft Anomaly of Interaction Enabled 3+1D Fermion SPTs

abstract: Symmetry protected topological (SPT) phases are well understood in the context of free fermions and in the context of interacting but essentially bosonic models. Recently it has been realized that intrinsically fermionic SPTs exist which only appear in interacting models. Here we show that the 3+1 dimensional realizations of these phases have surface states characterized by a new 't Hooft anomaly, captured by a $H^3(G, Z_2)$ class. This is encoded in the anomalous action of symmetry on the surface states with topological order, which must necessarily permute the anyons. We discuss in detail an example with symmetry group $G = Z_2 \times Z_4$. Using a network model of the surface we derive a candidate surface topological order given by a $Z_4$ gauge theory. We relate our findings to anomalies valued in $H^3$ with various coefficients introduced previously in both bosonic and fermionic settings, and describe a general framework that unifies these various anomalies.