title: Retrieval-Augmented Convolutional Neural Networks for Improved Robustness against Adversarial Examples

abstract: We propose a retrieval-augmented convolutional network and propose to train it with local mixup, a novel variant of the recently proposed mixup algorithm. The proposed hybrid architecture combining a convolutional network and an off-the-shelf retrieval engine was designed to mitigate the adverse effect of off-manifold adversarial examples, while the proposed local mixup addresses on-manifold ones by explicitly encouraging the classifier to locally behave linearly on the data manifold. Our evaluation of the proposed approach against five readily-available adversarial attacks on three datasets--CIFAR-10, SVHN and ImageNet--demonstrate the improved robustness compared to the vanilla convolutional network.


title: Online Coloring of Short Intervals

abstract: We study the online graph coloring problem restricted to the intersection graphs of intervals with lengths in $[1,\sigma]$. For $\sigma=1$ it is the class of unit interval graphs, and for $\sigma=\infty$ the class of all interval graphs. Our focus is on intermediary classes. We present a $(1+\sigma)$-competitive algorithm, which beats the state of the art for $1 < \sigma < 2$, and proves that the problem we study can be strictly easier than online coloring of general interval graphs. On the lower bound side, we prove that no algorithm is better than $5/3$-competitive for any $\sigma>1$, nor better than $7/4$-competitive for any $\sigma>2$, and that no algorithm beats the $5/2$ asymptotic competitive ratio for all, arbitrarily large, values of $\sigma$. That last result shows that the problem we study can be strictly harder than unit interval coloring. Our main technical contribution is a recursive composition of strategies, which seems essential to prove any lower bound higher than $2$.


title: Fast and accurate circularization of a Rydberg atom

abstract: Preparation of a so-called circular state in a Rydberg atom where the projection of the electron angular momentum takes its maximum value is challenging due to the required amount of angular momentum transfer. Currently available protocols for circular state preparation are either accurate but slow or fast but error-prone. Here, we show how to use quantum optimal control theory to derive pulse shapes that realize fast and accurate circularization of a Rydberg atom. In particular, we present a theoretical proposal for optimized radio-frequency pulses that achieve high fidelity in the shortest possible time, given current experimental limitations on peak amplitudes and spectral bandwidth. We also discuss the fundamental quantum speed limit for circularization of a Rydberg atom, when lifting these constraints.


title: Faster Algorithms for some Optimization Problems on Collinear Points

abstract: We propose faster algorithms for the following three optimization problems on $n$ collinear points, i.e., points in dimension one. The first two problems are known to be NP-hard in higher dimensions. 1- Maximizing total area of disjoint disks: In this problem the goal is to maximize the total area of nonoverlapping disks centered at the points. Acharyya, De, and Nandy (2017) presented an $O(n^2)$-time algorithm for this problem. We present an optimal $\Theta(n)$-time algorithm. 2- Minimizing sum of the radii of client-server coverage: The $n$ points are partitioned into two sets, namely clients and servers. The goal is to minimize the sum of the radii of disks centered at servers such that every client is in some disk, i.e., in the coverage range of some server. Lev-Tov and Peleg (2005) presented an $O(n^3)$-time algorithm for this problem. We present an $O(n^2)$-time algorithm, thereby improving the running time by a factor of $\Theta(n)$. 3- Minimizing total area of point-interval coverage: The $n$ input points belong to an interval $I$. The goal is to find a set of $n$ disks of minimum total area, covering $I$, such that every disk contains at least one input point. We present an algorithm that solves this problem in $O(n^2)$ time.


title: Relative volume comparison of Ricci Flow and its applications

abstract: In this paper, we derive a relative volume comparison estimate along Ricci flow and apply it to studying the Gromov-Hausdorff convergence of K\"ahler-Ricci flow on a minimal manifold. This new estimate generalizes Perelman's no local collapsing estimate and can be regarded as an analogue of the Bishop-Gromov volume comparison for Ricci flow.


title: Conjugacy growth of commutators

abstract: For the free group $F_r$ on $r>1$ generators (respectively, the free product $G_1 * G_2$ of two nontrivial finite groups $G_1$ and $G_2$), we obtain the asymptotic for the number of conjugacy classes of commutators in $F_r$ (respectively, $G_1 * G_2$) with a given word length in a fixed set of free generators (respectively, the set of generators given by the nontrivial elements of $G_1$ and $G_2$). Our result is proven by using the classification of commutators in free groups and in free products by Wicks, and builds on the works of Rivin and Sharp, who asymptotically counted the conjugacy classes of commutator-subgroup elements in $F_r$ with a given word length.


title: Image Reconstruction Techniques in Neutron and Gamma-Ray Spectroscopy: Improving Lunar Prospector Data

abstract: We present improved resolution maps of the Lunar Prospector Neutron Spectrometer thermal, epithermal and fast neutron data and Gamma-Ray Spectrometer Th-line fluxes via global application of pixon image reconstruction techniques. With the use of mock data sets, we show that the pixon image reconstruction method compares favorably with other methods that have been used in planetary neutron and gamma-ray spectroscopy. The improved thermal neutron maps are able to clearly distinguish variations in composition across the lunar surface, including within the lunar basins of Hertzsprung and Schrodinger. The improvement in resolution reveals a correlation between albedo and thermal neutron flux within the craters. The consequent increase in dynamic range confirms that Hertzsprung basin contains one of the most anorthositic parts of the lunar crust, including nearly pure anorthite over a region tens of km in diameter. At Orientale, the improvement in spatial resolution of the epithermal neutron data show that there is a mismatch between measures of regolith maturity that sample the surface and those that probe the near-subsurface, which suggests a complex layering scenario.


title: Estimation of Local Degree Distributions via Local Weighted Averaging and Monte Carlo Cross-Validation

abstract: Owing to their capability of summarising interactions between elements of a system, networks have become a common type of data in many fields. As networks can be inhomogeneous, in that different regions of the network may exhibit different topologies, an important topic concerns their local properties. This paper focuses on the estimation of the local degree distribution of a vertex in an inhomogeneous network. The contributions are twofold: we propose an estimator based on local weighted averaging, and we set up a Monte Carlo cross-validation procedure to pick the parameters of this estimator. Under a specific modelling assumption we derive an oracle inequality that shows how the model parameters affect the precision of the estimator. We illustrate our method by several numerical experiments, on both real and synthetic data, showing in particular that the approach considerably improves upon the natural, empirical estimator.


title: Bell measurement rules out supraquantum correlations

abstract: The so called bipartite non-signaling boxes are systems whose statistics is constrained solely by the principle of no instantaneous signaling between distant locations. Such systems can exhibit much stronger correlations than those admitted by quantum mechanics. Inspired by quantum logic approach of Tylec and Ku\'{s}, J. Phys. A: Math. Theor. 48 (2015) 505303, we consider non-signaling boxes with three inputs per party, and extend the set of measurements with just a single {\it global} measurement - one that mimics quantum two-party Bell measurement. We then show that this seemingly mild extension completely rules out supraquantum correlations: the resulting system admits precisely quantum mechanical correlations of two qubits. We also consider non-maximally entangled measurements, obtaining interpolation between quantum and full no-signaling theory. Our study paves a way to a general programme of amending no-signaling theories with some measurements inherited from quantum mechanics, leading to various interpolations between non-signaling boxes and quantum mechanics.


title: Missing Data in Sparse Transition Matrix Estimation for Sub-Gaussian Vector Autoregressive Processes

abstract: High-dimensional time series data exist in numerous areas such as finance, genomics, healthcare, and neuroscience. An unavoidable aspect of all such datasets is missing data, and dealing with this issue has been an important focus in statistics, control, and machine learning. In this work, we consider a high-dimensional estimation problem where a dynamical system, governed by a stable vector autoregressive model, is randomly and only partially observed at each time point. Our task amounts to estimating the transition matrix, which is assumed to be sparse. In such a scenario, where covariates are highly interdependent and partially missing, new theoretical challenges arise. While transition matrix estimation in vector autoregressive models has been studied previously, the missing data scenario requires separate efforts. Moreover, while transition matrix estimation can be studied from a high-dimensional sparse linear regression perspective, the covariates are highly dependent and existing results on regularized estimation with missing data from i.i.d.~covariates are not applicable. At the heart of our analysis lies 1) a novel concentration result when the innovation noise satisfies the convex concentration property, as well as 2) a new quantity for characterizing the interactions of the time-varying observation process with the underlying dynamical system.


title: Goldstone modes and photonization for higher form symmetries

abstract: We discuss generalized global symmetries and their breaking. We extend Goldstone's theorem to higher form symmetries by showing that a perimeter law for an extended $p$-dimensional defect operator charged under a continuous $p$-form generalized global symmetry necessarily results in a gapless mode in the spectrum. We also show that a $p$-form symmetry in a conformal theory in $2(p+1)$ dimensions has a free realization. In four dimensions this means any 1-form symmetry in a $CFT_4$ can be realized by free Maxwell electrodynamics, i.e. the current can be photonized. The photonized theory has infinitely many conserved 0-form charges that are constructed by integrating the symmetry currents against suitable 1-forms. We study these charges by developing a twistor-based formalism that is a 4d analogue of the usual holomorphic complex analysis familiar in $CFT_2$. The charges are shown to obey an algebra with central extension, which is an analogue of the 2d Abelian Kac-Moody algebra for higher form symmetries.


title: Typical and Generic Ranks in Matrix Completion

abstract: We consider the problem of exact low-rank matrix completion from a geometric viewpoint: given a partially filled matrix M, we keep the positions of specified and unspecified entries fixed, and study how the minimal completion rank depends on the values of the known entries. If the entries of the matrix are complex numbers, then for a fixed pattern of locations of specified and unspecified entries there is a unique completion rank which occurs with positive probability. We call this rank the generic completion rank. Over the real numbers there can be multiple ranks that occur with positive probability; we call them typical completion ranks. We introduce these notions formally, and provide a number of inequalities and exact results on typical and generic ranks for different families of patterns of known and unknown entries.

categories: math.ST cs.NA math.AG math.NA math.RA stat.TH


title: Best Arm Identification for Contaminated Bandits

abstract: This paper studies active learning in the context of robust statistics. Specifically, we propose a variant of the Best Arm Identification problem for \emph{contaminated bandits}, where each arm pull has probability $\varepsilon$ of generating a sample from an arbitrary contamination distribution instead of the true underlying distribution. The goal is to identify the best (or approximately best) true distribution with high probability, with a secondary goal of providing guarantees on the quality of this distribution. The primary challenge of the contaminated bandit setting is that the true distributions are only partially identifiable, even with infinite samples. To address this, we develop tight, non-asymptotic sample complexity bounds for high-probability estimation of the first two robust moments (median and median absolute deviation) from contaminated samples. These concentration inequalities are the main technical contributions of the paper and may be of independent interest. Using these results, we adapt several classical Best Arm Identification algorithms to the contaminated bandit setting and derive sample complexity upper bounds for our problem. Finally, we provide matching information-theoretic lower bounds on the sample complexity (up to a small logarithmic factor).


title: Representations of Sparse Distributed Networks: A Locality-Sensitive Approach

abstract: In 1999, Brodal and Fagerberg (BF) gave an algorithm for maintaining a low outdegree orientation of a dynamic uniformly sparse graph. Specifically, for a dynamic graph on $n$-vertices, with arboricity bounded by $\alpha$ at all times, the BF algorithm supports edge updates in $O(\log n)$ amortized update time, while keeping the maximum outdegree in the graph bounded by $O(\alpha)$. Such an orientation provides a basic data structure for uniformly sparse graphs, which found applications to a plethora of dynamic graph algorithms. A significant weakness of the BF algorithm is the possible \emph{temporary} blowup of the maximum outdegree, following edge insertions. Although BF eventually reduces all outdegrees to $O(\alpha)$, local memory usage at the vertices, which is an important quality measure in distributed systems, cannot be bounded. We show how to modify the BF algorithm to guarantee that the outdegrees of all vertices are bounded by $O(\alpha)$ at all times, without hurting any of its other properties, and present an efficient distributed implementation of the modified algorithm. This provides the \emph{first} representation of distributed networks in which the local memory usage at all vertices is bounded by the arboricity (which is essentially the average degree of the densest subgraph) rather than the maximum degree. For settings where there are no local memory constraints, we take the temporary outdegree blowup to the extreme and allow a permanent outdegree blowup. This allows us to address the second significant weakness of the BF algorithm -- its inherently \emph{global} nature: An insertion of an edge $(u,v)$ may trigger changes in the orientations of edges that are far away from $u$ and $v$. We suggest an alternative \emph{local} scheme, which does not guarantee any outdegree bound on the vertices, yet is just as efficient as the BF scheme for various applications.


title: Bell correlations depth in many-body systems

abstract: While the interest in multipartite nonlocality has grown in recent years, its existence in large quantum systems is difficult to confirm experimentally. This is mostly due to the inadequacy of standard multipartite Bell inequalities to many-body systems: such inequalities usually rely on expectation values involving many parties and require an individual addressing of each party. In a recent work [J. Tura et al. Science 344, 6189 (2014)] some of us proposed simpler Bell inequalities overcoming such difficulties, opening the way for the detection of Bell correlations with trusted collective measurements through Bell correlation witnesses [R. Schmied et al. Science 352, 441 (2016)], hence demonstrating the presence of Bell correlations with assumptions on the statistics. Here, we address the question of assessing the number of particles sharing genuinely nonlocal correlations in a multipartite system. This endeavour is a priori challenging, as known Bell inequalities for genuine nonlocality suffer from the above shortcomings, plus a number of measurement settings scaling exponentially with the system size. We first show that most of these constraints drop once the witnesses corresponding to these inequalities are expressed: in systems where multipartite expectation values can be evaluated, these witnesses can reveal genuine nonlocality for an arbitrary number of particles with just two collective measurements. We then introduce a general framework focused on two-body Bell-like inequalities. We show that they also provide information about the number of particles that are genuinely nonlocal. Then, we characterize all such inequalities for a finite system size. We provide witnesses of Bell correlation depth $k\leq6$ for any number of parties, within experimental reach. A violation for depth $6$ is achieved with existing data from an ensemble of 480 atoms.


title: Memory Tagging and how it improves C/C++ memory safety

abstract: Memory safety in C and C++ remains largely unresolved. A technique usually called "memory tagging" may dramatically improve the situation if implemented in hardware with reasonable overhead. This paper describes two existing implementations of memory tagging: one is the full hardware implementation in SPARC; the other is a partially hardware-assisted compiler-based tool for AArch64. We describe the basic idea, evaluate the two implementations, and explain how they improve memory safety. This paper is intended to initiate a wider discussion of memory tagging and to motivate the CPU and OS vendors to add support for it in the near future.


title: Orthogonal Basis Function Over the Unit Circle with the Minimax Property

abstract: We construct an orthogonal basis of functions defined over the unit circle as the product of the common sinusoidal functions of the azimuth angle by radial functions which are essentially sines of a polynomials of the radial distance to the origin. The main impetus of this approach is to generate basis functions where the minima and maxima along both coordinates, the azimuth and the distance r to the center, have the same amplitude, akin to the Chebyshev polynomial basis of the one-dimensional unit interval. The construction is based on numerical evaluation of the overlap integrals, which have the format of generalized Fresnel integrals.


title: Leading-order Stokes flows near a corner

abstract: Singular solutions of the Stokes equations play important roles in a variety of fluid dynamics problems. They allow the calculation of exact flows, are the basis of the boundary integral methods used in numerical computations, and can be exploited to derive asymptotic flows in a wide range of physical problems. The most fundamental singular solution is the flow's Green function due to a point force, termed the Stokeslet. Its expression is classical both in free space and near a flat surface. Motivated by problems in biological physics occurring near corners, we derive in this paper the asymptotic behaviour for the Stokeslet both near and far from a corner geometry by using complex analysis on a known double integral solution for corner flows. We investigate all possible orientations of the point force relative to the corner and all corner geometries from acute to obtuse. The case of salient corners is also addressed for point forces aligned with both walls. We use experiments on beads sedimenting in corn syrup to qualitatively test the applicability of our results. The final results and scaling laws will allow to address the role of hydrodynamic interactions in problems from colloidal science to microfluidics and biological physics.


title: A stable flat universe with variable Cosmological constant in f(R,T) Gravity

abstract: In this paper, a general FRW cosmological model has been constructed in $f(R,T)$ gravity reconstruction with variable cosmological constant. A number of solutions to the field equations has been generated by utilizing a form for the Hubble parameter that leads to Berman's law of constant deceleration parameter $q=m-1$. The possible decelerating and accelerating solutions have been investigated. For ($q>0$) we get a stable flat decelerating radiation-dominated universe at $q=1$. For ($q<0$) we get a stable accelerating solution describing a flat universe with positive energy density and negative cosmological constant. Nonconventional mechanisms that are expected to address the late-time acceleration with negative cosmological constant have been discussed.


title: Moving mesh finite difference solution of non-equilibrium radiation diffusion equations

abstract: A moving mesh finite difference method based on the moving mesh partial differential equation is proposed for the numerical solution of the 2T model for multi-material, non-equilibrium radiation diffusion equations. The model involves nonlinear diffusion coefficients and its solutions stay positive for all time when they are positive initially. Nonlinear diffusion and preservation of solution positivity pose challenges in the numerical solution of the model. A coefficient-freezing predictor-corrector method is used for nonlinear diffusion while a cutoff strategy with a positive threshold is used to keep the solutions positive. Furthermore, a two-level moving mesh strategy and a sparse matrix solver are used to improve the efficiency of the computation. Numerical results for a selection of examples of multi-material non-equilibrium radiation diffusion show that the method is capable of capturing the profiles and local structures of Marshak waves with adequate mesh concentration. The obtained numerical solutions are in good agreement with those in the existing literature. Comparison studies are also made between uniform and adaptive moving meshes and between one-level and two-level moving meshes.


title: The Hiperwall Visualization Platform for Big Data Research

abstract: In the era of Big Data, with the increasing use of large-scale data-driven applications, the visualization of very large high-resolution images and extracting useful information (searching for specific targets or rare signal events) from these images can pose challenges to the current video-wall display technologies. At Bellarmine University, we have set up an Advanced Visualization and Computational Lab (AVCL) using a state-of-the-art next generation video-wall technology, called Hiperwall (Highly Interactive Parallelized Display Wall). The 16 feet wide by 4.5 feet high Hiperwall visualization system consists of eight display tiles that are arranged in a 4x2 tile format and has an effective resolution of 16.5 Megapixels. Using Hiperwall, we can perform interactive visual data analytics of large images by conducting comparative views of multiple large images in Astronomy and multiple data events in experimental High Energy Physics (HEP). Users can display a single large image across all the display tiles, or view many different images simultaneously on multiple display tiles. Hiperwall enables simultaneous visualization of multiple high resolution images and its contents on the entire display wall without loss of clarity. Hiperwall's middleware also allows researchers in geographically diverse locations to collaborate on large scientific experiments. In this paper we will provide a description of a new generation of display wall setup at Bellarmine University that is based on the Hiperwall technology, which is a robust visualization system for Big Data research.


title: Interaction of orthogonal-polarized waves in 1D metamaterial with Kerr nonlinearity

abstract: A theoretical study of wave propagation in 1D metamaterial is presented. A system of nonlinear evolution equation for electromagnetic waves with both polarizations account is derived by means of projection operators method for general nonlinearity and dispersion. The system describes interaction of opposite directed waves with a given polarization. The particular case of Kerr nonlinearity and Drude dispersion is considered. In such approximation it results in the correspondent systems of nonlinear equations that generalizes the Sch\"{a}fer-Wayne one. Particular solutions in case of slow-varying envelopes are found, plotted and analyzed in gigahertz range. Travelling wave solution for the system of equation of interaction of orthogonal-polarized waves is also obtained and the correspondent nonlinear dispersion relations are written in explicit form.


title: Magnon-induced superconductivity in a topological insulator coupled to ferro- and antiferromagnetic insulators

abstract: We study the effective interactions between Dirac fermions on the surface of a three-dimensional topological insulator due to the proximity coupling to the magnetic fluctuations in a ferromagnetic or antiferromagnetic insulator. Our results show that the magnetic fluctuations can mediate attractive interactions between Dirac fermions of both Amperean and BCS type. In the ferromagnetic case, we find pairing between fermions with parallel momenta, so-called Amperean pairing, whenever the effective Lagrangian for the magnetic fluctuations does not contain a quadratic term. The pairing interaction also increases with increasing Fermi momentum, and is in agreement with previous studies in the limit of high chemical potential. If a quadratic term is present, the pairing is instead of BCS type above a certain chemical potential. In the antiferromagnetic case, BCS pairing occurs when the ferromagnetic coupling between magnons on the same sublattice exceeds the antiferromagnetic coupling between magnons on different sublattices. Outside this region in parameter space, we again find that Amperean pairing is realized.


title: Associating Fast Radio Bursts with Extragalactic Radio Sources: General Methodology and a Search for a Counterpart to FRB 170107

abstract: The discovery of a repeating fast radio burst has led to the first precise localization, an association with a dwarf galaxy, and the identification of a coincident persistent radio source. However, further localizations are required to determine the nature of FRBs, the sources powering them, and the possibility of multiple populations. Here we investigate the use of associated persistent radio sources to establish FRB counterparts, taking into account the localization area and the persistent source flux density. Due to the lower areal number density of radio sources compared to faint optical sources, robust associations can be achieved for less precise localizations as compared to direct optical host galaxy associations. For generally larger localizations which preclude robust associations, the number of candidate hosts can be reduced based on the ratio of radio-to-optical brightness. We find that confident associations with $\sim 0.01-$1 mJy sources, comparable to the luminosity of the persistent source associated with FRB 121102 over the redshift range $z \approx 0.1 - 1$, require FRB localizations of $\lesssim 20''$. In the absence of a robust association, constraints can be placed on the luminosity of an associated radio source as a function of localization and DM. For DM $\approx 1000 \rm \ pc \ cm^{-3}$, an upper limit comparable to the luminosity of the FRB 121102 persistent source can be placed if the localization is $\lesssim 10''$. We apply our analysis to the case of the ASKAP FRB 170107, using optical and radio observations of the localization region. We identify two candidate hosts based on a ratio of radio-to-optical brightness of $\gtrsim 100$. We find that if one of these is associated with FRB 170107, the resulting radio luminosity ($1 \times 10^{29} - 4 \times 10^{30} \ \rm erg \ s^{-1} \ Hz^{-1}$) is comparable to the luminosity of the FRB 121102 persistent source.


title: About 30% of Sun-like Stars Have Kepler-like Planetary Systems: A Study of their Intrinsic Architecture

abstract: We constrain the intrinsic architecture of Kepler planetary systems by modeling the observed multiplicities of the transiting planets (tranets) and their transit timing variations (TTVs). We robustly determine that the fraction of Sun-like stars with Kepler-like planets, $\eta_{\rm Kepler}$, is $30\pm3\%$. Here Kepler-like planets are planets that have radii $R_{\rm p} \gtrsim R_\oplus$ and orbital periods $P<400$~days. Our result thus significantly revises previous claims that more than 50\% of Sun-like stars have such planets. Combining with the average number of Kepler planets per star ($\sim0.9$), we obtain that on average each planetary system has $3.0\pm0.3$ planets within 400 days. We also find that the dispersion in orbital inclinations of planets within a given planetary system, $\sigma_{i,k}$, is a steep function of its number of planets, $k$. This can be parameterized as $\sigma_{i,k}\propto k^\alpha$ and we find that $-4<\alpha<-2$ at 2-$\sigma$ level. Such a distribution well describes the observed multiplicities of both tranets and TTVs with no excess of single tranets. Therefore we do not find evidence supporting the so-called "Kepler dichotomy." Together with a previous study on orbital eccentricities, we now have a consistent picture: the fewer planets in a system, the hotter it is dynamically. We discuss briefly possible scenarios that lead to such a trend. Despite our Solar system not belonging to the Kepler club, it is interesting to notice that the Solar system also has three planets within 400 days and that the inclination dispersion is similar to Kepler systems of the same multiplicity.


title: Quantum simulation of zero temperature quantum phases and incompressible states of light via non-Markovian reservoir engineering techniques

abstract: We review recent theoretical developments on the stabilization of strongly correlated quantum fluids of light in driven-dissipative photonic devices through novel non-Markovian reservoir engineering techniques. This approach allows to compensate losses and refill selectively the photonic population so to sustain a desired steady-state. It relies in particular on the use of a frequency-dependent incoherent pump which can be implemented, e.g., via embedded two-level systems maintained at a strong inversion of population. As specific applications of these methods, we discuss the generation of Mott Insulator (MI) and Fractional Quantum Hall (FQH) states of light. As a first step, we present the case of a narrowband emission spectrum and show how this allows for the stabilization of MI and FQH states under the condition that the photonic states are relatively flat in energy. As soon as the photonic bandbwidth becomes comparable to the emission linewidth, important non-equilibrium signatures and entropy generation appear. As a second step, we review a more advanced configuration based on reservoirs with a broadband frequency distribution, and we highlight the potential of this configuration for the quantum simulation of equilibrium quantum phases at zero temperature with tunable chemical potential. As a proof of principle we establish the applicability of our scheme to the Bose-Hubbard model by confirming the presence of a perfect agreement with the ground-state predictions both in the Mott Insulating and superfluid regions, and more generally in all parts of the parameter space. Future prospects towards the quantum simulation of more complex configurations are finally outlined, along with a discussion of our scheme as a concrete realization of quantum annealing.


title: The Monte Carlo photoionization and moving-mesh radiation hydrodynamics code CMacIonize

abstract: We present the public Monte Carlo photoionization and moving-mesh radiation hydrodynamics code CMacIonize, which can be used to simulate the self-consistent evolution of HII regions surrounding young O and B stars, or other sources of ionizing radiation. The code combines a Monte Carlo photoionization algorithm that uses a complex mix of hydrogen, helium and several coolants in order to self-consistently solve for the ionization and temperature balance at any given type, with a standard first order hydrodynamics scheme. The code can be run as a post-processing tool to get the line emission from an existing simulation snapshot, but can also be used to run full radiation hydrodynamical simulations. Both the radiation transfer and the hydrodynamics are implemented in a general way that is independent of the grid structure that is used to discretize the system, allowing it to be run both as a standard fixed grid code, but also as a moving-mesh code.


title: Identifying Young $Kepler$ Planet Host Stars from Keck-HIRES Spectra of Lithium

abstract: The lithium doublet at 6708 \AA\ provides an age diagnostic for main sequence FGK dwarfs. We measured the abundance of lithium in 1305 stars with detected transiting planets from the Kepler Mission using high-resolution spectroscopy. Our catalog of lithium measurements from this sample have a range of abundance from A(Li) = 3.11 $\pm$ 0.07 to an upper limit of $-$0.84 dex. For a magnitude-limited sample that comprises 960 of the 1305 stars, our Keck-HIRES spectra have a median S/N = 45 per pixel at $\sim$6700 \AA\ with spectral resolution $\frac{\lambda}{\Delta \lambda}$ = $R$ = 55,000. We identify 80 young stars that have A(Li) values greater than the Hyades at their respective effective temperatures; these stars are younger than $\sim$650 Myr old, the approximate age of the Hyades. We then compare the distribution of A(Li) with planet size, multiplicity, orbital period, and insolation flux. We find larger planets preferentially in younger systems, with an A-D two-sided test p-value = 0.002, a $>3\sigma$ confidence that the older and younger planet samples do not come from the same parent distribution. This is consistent with planet inflation/photoevaporation at early ages. The other planet parameters ($Kepler$ planet multiplicity, orbital period, and insolation flux) are uncorrelated with age.


title: The frequency of very young galaxies in the local Universe: I. A test for galaxy formation and cosmological models

abstract: In the local Universe, the existence of very young galaxies (VYGs), having formed at least half their stellar mass in the last 1 Gyr, is debated. We predict the present-day fraction of VYGs among central galaxies as a function of galaxy stellar mass. For this, we apply to high mass resolution Monte-Carlo halo merger trees (MCHMTs) three (one) analytical models of galaxy formation, where the ratio of stellar to halo mass (mass growth rate) is a function of halo mass and redshift. Galaxy merging is delayed until orbital decay by dynamical friction. With starbursts associated with halo mergers, our models predict typically one percent of VYGs up to galaxy masses of $10^{10}$ M$_\odot$, falling rapidly at higher masses, and VYGs are usually associated with recent major mergers of their haloes. Without these starbursts, two of the models have VYG fractions reduced by 1 or 2 dex at low or intermediate stellar masses, and VYGs are rarely associated with major halo mergers. In comparison, the state-of-the-art semi-analytical model (SAM) of Henriques et al. produces only 0.01% of VYGs at intermediate masses. Finally, the Menci et al. SAM run on MCMHTs with Warm Dark Matter cosmology generates 10 times more VYGs at masses below $10^8$ M$_\odot$ than when run with Cold Dark Matter. The wide range in these VYG fractions illustrates the usefulness of VYGs to constrain both galaxy formation and cosmological models.


title: Phase Space Entanglement Spectrum

abstract: We generalize the position- and momentum-space entanglement cuts to a family of cuts corresponding to regions in the classical phase space. We explicitly compute the entanglement spectra of free fermionic many-body wavefunctions for a family of phase space entanglement cuts that continuously interpolates between position- and momentum-space cuts. For inversion symmetric wavefunctions, the phase space entanglement spectrum possess a chiral symmetry, to which a topological index can be associated.

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


title: Large velocity dispersion of molecular gas in bars of strongly barred galaxies, NGC 1300 and NGC 5383

abstract: We carried out $^{12}$CO(J = 1 -0) observations toward bar and arm regions of strongly barred galaxies, NGC 1300 and NGC 5383, with the Nobeyama 45-m radio telescope (beamsize of 1-2 kpc in the galaxies). The aim of the observations is to qualitatively examine a new scenario for the suppression of star formation in bars based on recent high-resolution numerical simulations: higher speed collisions between molecular clouds in the bar region compared with the arm region suppress the massive star formation. CO emissions were detected from all the regions, indicating the presence of the molecular gases in the strong bars without associating clear HII regions. In both galaxies, the velocity width of the CO line profile tends to be larger in the bar region than in the arm region, which is qualitatively consistent with the new scenario.


title: The Sagittarius dwarf galaxy: Where did all the gas go?

abstract: The remarkable 1994 discovery of the Sagittarius dwarf galaxy (Sgr) revealed that, together with the Magellanic Clouds (MCs), there are at least three major dwarf galaxies, each with a total mass of order 10^10 - 10^11 solar masses, falling onto the Galaxy in the present epoch. Beyond a Galactic radius of 300 kpc, dwarfs tend to retain their gas. At roughly 50 kpc, the MCs have experienced substantial gas stripping as evidenced by the Magellanic Stream which extends from them. Since Sgr experienced star formation long after it fell into the Galaxy, it is interesting to explore just how and when this dwarf lost its gas. To date, there has been no definitive detection of an associated gas component. We revisit recent simulations of the stellar and dark matter components of Sgr but, for the first time, include gas that is initially bound to the infalling galaxy. We find that the gas stripping was 30 - 50 % complete at its first disc crossing ~2.7 Gyr ago, then entirely stripped at its last disc crossing ~1 Gyr ago. Our timeline is consistent with the last substantial burst of star formation in Sgr which occurred about the time of the last disc crossing. We discuss the consequences of gas stripping and conclude that the vast majority of the stripped gas was fully settled onto the Galaxy by ~300 Myr ago. It is highly unlikely that any of the high- or intermediate-velocity clouds have a direct association with the Sgr dwarf.


title: Chemical Abundances of Globular Clusters in NGC 5128 (Centaurus A)

abstract: We perform a detailed abundance analysis on integrated-light spectra of 20 globular clusters (GCs) in the early-type galaxy NGC 5128 (Centaurus A). The GCs were observed with X-Shooter on the VLT. The cluster sample spans a metallicity range of $-1.92 < $ [Fe/H] $< -0.13$ dex. Using theoretical isochrones we compute synthetic integrated-light spectra and iterate the individual abundances until the best fit to the observations is obtained. We measured abundances of Mg, Ca, and Ti, and find a slightly higher enhancement in NGC 5128 GCs with metallicities [Fe/H] < $-$0.75 dex, of the order of $\sim$0.1 dex, than in the average values observed in the MW for GCs of the same metallicity. If this $\alpha$-enhancement in the metal-poor GCs in NGC 5128 is genuine, it could hint at a chemical enrichment history different than that experienced by the MW. We also measure Na abundances in 9 out of 20 GCs. We find evidence for intra-cluster abundance variations in 6 of these clusters where we see enhanced [Na/Fe] > $+$0.25 dex. We obtain the first abundance measurements of Cr, Mn, and Ni for a sample of the GC population in NGC 5128 and find consistency with the overall trends observed in the MW, with a slight enhancement ($<$0.1 dex) in the Fe-peak abundances measured in the NGC 5128.


title: QCD Shear Viscosity at (almost) NLO

abstract: We compute the shear viscosity of QCD with matter, including almost all next-to-leading order corrections -- that is, corrections suppressed by one power of $g$ relative to leading order. We argue that the still missing terms are small. The next-to-leading order corrections are large and bring $\eta/s$ down by more than a factor of 3 at physically relevant couplings. The perturbative expansion is problematic even at $T \simeq 100$ GeV. The largest next-to-leading order correction to $\eta/s$ arises from modifications to the qhat parameter, which determines the rate of transverse momentum diffusion. We also explore quark number diffusion, and shear viscosity in pure-glue QCD and in QED.


title: Astrophysics with New Horizons: Making the Most of a Generational Opportunity

abstract: The outer solar system provides a unique, quiet vantage point from which to observe the universe around us, where measurements could enable several niche astrophysical science cases that are too difficult to perform near Earth. NASA's New Horizons mission comprises an instrument package that provides imaging capability from ultraviolet (UV) to near-infrared (near-IR) wavelengths with moderate spectral resolution located beyond the orbit of Pluto. A carefully designed survey with New Horizons can optimize the use of expendable propellant and the limited data telemetry bandwidth to allow several measurements, including a detailed understanding of the cosmic extragalactic background light; studies of the local and extragalactic UV background; measurements of the properties of dust and ice in the outer solar system; confirmation and characterization of transiting exoplanets; determinations of the mass of dark objects using gravitational microlensing; and rapid follow-up of transient events. New Horizons is currently in an extended mission designed to focus on Kuiper Belt science that will conclude in 2021. The astrophysics community has a unique, generational opportunity to use this mission for astronomical observation at heliocentric distances beyond 50 au in the next decade. In this paper, we discuss the potential science cases for such an extended mission, and provide an initial assessment of the most important operational requirements and observation strategies it would require. We conclude that New Horizons is capable of transformative science, and that it would make a valuable and unique asset for astrophysical science that is unlikely to be replicated in the near future.

categories: astro-ph.IM astro-ph.CO astro-ph.EP astro-ph.GA astro-ph.SR


title: Tidal disruption of dwarf spheroidal galaxies: the strange case of Crater II

abstract: Dwarf spheroidal galaxies of the Local Group obey a relationship between the line-of-sight velocity dispersion and half-light radius, although there are a number of dwarfs that lie beneath this relation with suppressed velocity dispersion. The most discrepant of these (in the Milky Way) is the `feeble giant' Crater II. Using analytic arguments supported by controlled numerical simulations of tidally-stripped flattened two-component dwarf galaxies, we investigate interpretations of Crater II within standard galaxy formation theory. Heavy tidal disruption is necessary to explain the velocity-dispersion suppression which is plausible if the proper motion of Crater II is $(\mu_{\alpha*},\mu_\delta)=(-0.21\pm0.09,-0.24\pm0.09)\mathrm{mas\,yr}^{-1}$. Furthermore, we demonstrate that the velocity dispersion of tidally-disrupted systems is solely a function of the total mass loss even for weakly-embedded and flattened systems. The half-light radius evolution depends more sensitively on orbital phase and the properties of the dark matter profile. The half-light radius of weakly-embedded cusped systems rapidly decreases producing some tension with the Crater II observations. This tension is alleviated by cored dark matter profiles, in which the half-light radius can grow after tidal disruption. The evolution of flattened galaxies is characterised by two competing effects: tidal shocking makes the central regions rounder whilst tidal distortion produces a prolate tidally-locked outer envelope. After $\sim70$ per cent of the central mass is lost, tidal distortion becomes the dominant effect and the shape of the central regions of the galaxy tends to a universal prolate shape irrespective of the initial shape.


title: Emergent Multi-flavor QED3 at the Plateau Transition between Fractional Chern Insulators: Applications to graphene heterostructures

abstract: Recent experiments in graphene heterostructures have observed Chern insulators - integer and fractional Quantum Hall states made possible by a periodic substrate potential. Here we study theoretically the competition between different Chern insulators, which can be tuned by the amplitude of the periodic potential, leads to a new family of quantum critical points described by QED$_3$-Chern-Simons theory. At these critical points, $N_f$ flavors of Dirac fermions interact through an emergent U$(1)$ gauge theory at Chern-Simons level $K$, and remarkably, the entire family (with any $N_f$ or $K$) can be realized at special values of the external magnetic field. Transitions between particle-hole conjugate Jain states realize "pure" QED$_3$ in which multiple flavors of Dirac fermion interact with a Maxwell U$(1)$ gauge field. The multi-flavor nature of the critical point leads to an emergent SU$(N_f)$ symmetry. Specifically, at the transition from a $\nu=$1/3 to 2/3 quantum Hall state, the emergent SU(3) symmetry predicts an octet of charge density waves with enhanced susceptibilities, which is verified by DMRG numerical simulations on microscopic models applicable to graphene heterostructures. We propose experiments on Chern insulators that could resolve open questions in the study of 2+1 dimensional conformal field theories and test recent duality inspired conjectures.


title: Precision Holography for $\mathcal{N}=2^{*}$ on $S^4$ from type IIB Supergravity

abstract: We find a new supersymmetric solution of type IIB supergravity which is holographically dual to the planar limit of the four-dimensional $\mathcal{N}=2^*$ supersymmetric Yang-Mills theory on $S^4$. We study a probe fundamental string in this background which is dual to a supersymmetric Wilson loop in the $\mathcal{N}=2^*$ theory. Using holography we calculate the expectation value of this line operator to leading order in the 't Hooft coupling. The result is a non-trivial function of the mass parameter of the $\mathcal{N}=2^*$ theory that precisely matches the result from supersymmetric localization.


title: The hELENa project - II. Abundance distribution trends of early-type galaxies: from dwarfs to giants

abstract: In this second paper of T$h$e role of $E$nvironment in shaping $L$ow-mass $E$arly-type $N$earby g$a$laxies (hELENa) series we study [Mg/Fe] abundance distribution trends of early-type galaxies observed with the SAURON integral field unit, spanning a wide range in mass and local environment densities: 20 low-mass early-types (dEs) of Sybilska et al. (2017) and 258 massive early types (ETGs) of the $ATLAS^{3D}$ project, all homogeneously reduced and analyzed. We show that the [Mg/Fe] ratios scale with velocity dispersion ($\sigma$) at fixed [Fe/H] and that they evolve with [Fe/H] along similar paths for all early-types, grouped in bins of increasing local and global {$\sigma$}, as well as the second velocity moment $V_{rms}$, indicating a common inside-out formation pattern. We then place our dEs on the [Mg/Fe] $vs.$ [Fe/H] diagram of Local Group galaxies and show that dEs occupy the same region and show a similar trend line slope in the diagram as the high-metallicity stars of the Milky Way and the Large Magellanic Cloud. This finding extends the similar trend found for dwarf spheroidal $vs.$ dwarf irregular galaxies and supports the notion that dEs have evolved from late-type galaxies that have lost their gas at a point of their evolution, which likely coincided with them entering denser environments.


title: The Dark Matter equation of state through cosmic history

abstract: Cold Dark Matter (CDM) is a crucial constituent of the current concordance cosmological model. Having a vanishing equation of state (EoS), its energy density scales with the inverse cosmic volume and is thus uniquely described by a single number, its present abundance. We test the inverse cosmic volume law for Dark Matter (DM) by allowing its EoS to vary independently in eight redshift bins in the range $z=10^5$ and $z=0$. We use the latest measurements of the Cosmic Microwave Background radiation from the Planck satellite and supplement them with Baryon Acoustic Oscillation (BAO) data from the 6dF and SDSS-III BOSS surveys, and with the Hubble Space Telescope (HST) key project data. We find no evidence for nonzero EoS in any of the eight redshift bins. With Planck data alone, the DM abundance is most strongly constrained around matter-radiation equality $\omega^{\rm eq}_g = 0.1193^{+0.0036}_{-0.0035}$ (95% c.l.), whereas its present day value is more weakly constrained $\omega^{(0)}_g = 0.16^{+0.12}_{-0.10}$ (95% c.l.). Adding BAO or HST data does not significantly change the $\omega^{\rm eq}_g$ constraint, while $\omega^{(0)}_g$ tightens to $0.160^{+0.069}_{-0.065} $ (95% c.l.) and $0.124^{+0.081}_{-0.067}$ (95% c.l.) respectively. Our results constrain for the first time the level of "coldness" required of the DM across various cosmological epochs and show that the DM abundance is strictly positive at all times.


title: Non-ergodic delocalized states for efficient population transfer within a narrow band of the energy landscape

abstract: We analyze the role of coherent tunneling that gives rise to bands of delocalized quantum states providing a coherent pathway for population transfer (PT) between computational states with similar energies. Given an energy function ${\cal E}(z)$ of a binary optimization problem and a bit-string $z_i$ with atypically low energy, our goal is to find other bit-strings with energies within a narrow window around ${\cal E}(z_i)$. We study PT due to quantum evolution under a transverse field $B_\perp$ of an n-qubit system that encodes ${\cal E}(z)$. We focus on a simple yet nontrivial model: $M$ randomly chosen "marked" bit-strings ($2^n \gg M$) are assigned energies in the interval ${\cal E}(z)\in[-n -W/2, n + W/2]$ with $W << B_\perp$, while the rest of the states are assigned energy $0$. The PT starts at a marked state $z_i$ and ends up in a superposition of $\sim \Omega$ marked states inside the PT window. The scaling of a typical runtime for PT with $n$ and $\Omega$ is the same as in the multi-target Grover's algorithm, except for a factor that is equal to $\exp(n \,B_{\perp}^{-2}/2)$ for $n \gg B_{\perp}^{2} \gg 1$. Unlike the Hamiltonians used in analog quantum search algorithms, the model we consider is non-integrable, and the transverse field delocalizes the marked states. PT protocol is not sensitive to the value of B and may be initialized at a marked state. We develop microscopic theory of PT. Under certain conditions, the band of the system eigenstates splits into mini-bands of non-ergodic delocalized states, whose width obeys a heavy-tailed distribution directly related to that of PT runtimes. We find analytical form of this distribution by solving nonlinear cavity equations for the random matrix ensemble. We argue that our approach can be applied to study the PT protocol in other transverse field spin glass models, with a potential quantum advantage over classical algorithms.


title: Theory of thermionic emission from a two-dimensional conductor and its application to a graphene-semiconductor Schottky junction

abstract: The standard theory of thermionic emission developed for three-dimensional semiconductors does not apply to two-dimensional materials even for making qualitative predictions because of the vanishing out-of-plane quasiparticle velocity. This study reveals the fundamental origin of the out-of-plane charge carrier motion in a two-dimensional conductor due to the finite quasiparticle lifetime and huge uncertainty of the out-of-plane momentum. The theory is applied to a Schottky junction between graphene and a bulk semiconductor to derive a thermionic constant, which, in contrast to the conventional Richardson constant, is determined by the Schottky barrier height and Fermi level in graphene.


title: Lagrangian formulation, generalizations and quantization of null Maxwell's knots

abstract: Knotted solutions to electromagnetism are investigated as an independent subsector of the theory. We write down a Lagrangian and a Hamiltonian formulation of Bateman's construction for the knotted electromagnetic solutions. We introduce a general definition of the null condition and generalize the construction of Maxwell's theory to massless free complex scalar, its dual two form field, and to a massless DBI scalar. We set up the framework for quantizing the theory both in a path integral approach, as well as the canonical Dirac method for a constrained system. We make several observations about the semi-classical quantization of systems of null configurations.


title: Entanglement of Purification in Free Scalar Field Theories

abstract: We compute the entanglement of purification (EoP) in a 2d free scalar field theory with various masses. This quantity measures correlations between two subsystems and is reduced to the entanglement entropy when the total system is pure. We obtain explicit numerical values by assuming minimal gaussian wave functionals for the purified states. We find that when the distance between the subsystems is large, the EoP behaves like the mutual information. However, when the distance is small, the EoP shows a characteristic behavior which qualitatively agrees with the conjectured holographic computation and which is different from that of the mutual information. We also study behaviors of mutual information in purified spaces and violations of monogamy/strong superadditivity.


title: Quantum and Classical Phases of the Pyrochlore Heisenberg Model with Competing Interactions

abstract: We investigate the quantum Heisenberg model on the pyrochlore lattice for a generic spin $S$ in the presence of nearest-neighbor $J_{1}$ and second-nearest-neighbor $J_{2}$ exchange interactions. By employing the pseudofermion functional renormalization group method, we find, for $S=1/2$ and $S=1$, an extended quantum-spin-liquid phase centered around $J_{2}=0$, which is shown to be robust against the introduction of breathing anisotropy. The effects of temperature, quantum fluctuations, breathing anisotropies, and a $J_{2}$ coupling on the nature of the scattering profile, and the pinch points, in particular, are studied. For the magnetic phases of the $J_{1}$-$J_{2}$ model, quantum fluctuations are shown to renormalize phase boundaries compared to the classical model and to modify the ordering wave vectors of spiral magnetic states, while no new magnetic orders are stabilized.


title: On the Dynamics of Near-Extremal Black Holes

abstract: We analyse the dynamics of near-extremal Reissner-Nordstr\"om black holes in asymptotically four-dimensional Anti-de Sitter space (AdS$_4$). We work in the spherically symmetric approximation and study the thermodynamics and the response to a probe scalar field. We find that the behaviour of the system, at low energies and to leading order in our approximations, is well described by the Jackiw-Teitelboim (JT) model of gravity. In fact, this behaviour can be understood from symmetry considerations and arises due to the breaking of time reparametrisation invariance. The JT model has been analysed in considerable detail recently and related to the behaviour of the SYK model. Our results indicate that features in these models which arise from symmetry considerations alone are more general and present quite universally in near-extremal black holes.


title: Human Perceptions of Fairness in Algorithmic Decision Making: A Case Study of Criminal Risk Prediction

abstract: As algorithms are increasingly used to make important decisions that affect human lives, ranging from social benefit assignment to predicting risk of criminal recidivism, concerns have been raised about the fairness of algorithmic decision making. Most prior works on algorithmic fairness normatively prescribe how fair decisions ought to be made. In contrast, here, we descriptively survey users for how they perceive and reason about fairness in algorithmic decision making. A key contribution of this work is the framework we propose to understand why people perceive certain features as fair or unfair to be used in algorithms. Our framework identifies eight properties of features, such as relevance, volitionality and reliability, as latent considerations that inform people's moral judgments about the fairness of feature use in decision-making algorithms. We validate our framework through a series of scenario-based surveys with 576 people. We find that, based on a person's assessment of the eight latent properties of a feature in our exemplar scenario, we can accurately (> 85%) predict if the person will judge the use of the feature as fair. Our findings have important implications. At a high-level, we show that people's unfairness concerns are multi-dimensional and argue that future studies need to address unfairness concerns beyond discrimination. At a low-level, we find considerable disagreements in people's fairness judgments. We identify root causes of the disagreements, and note possible pathways to resolve them.


title: Hyperbolic Lattices in Circuit Quantum Electrodynamics

abstract: After close to two decades of research and development, superconducting circuits have emerged as a rich platform for both quantum computation and quantum simulation. Lattices of superconducting coplanar waveguide (CPW) resonators have been shown to produce artificial materials for microwave photons, where weak interactions can be introduced either via non-linear resonator materials or strong interactions via qubit-resonator coupling. Here, we highlight the previously-overlooked property that these lattice sites are deformable and allow the realization of tight-binding lattices which are unattainable, even in conventional solid-state systems. In particular, we show that networks of CPW resonators can create a new class of materials which constitute regular lattices in an effective hyperbolic space with constant negative curvature. We present numerical simulations of a series of hyperbolic analogs of the kagome lattice which show unusual densities of states with a spectrally-isolated degenerate flat band. We also present a proof-of-principle experimental realization of one of these lattices. This paper represents the first step towards on-chip quantum simulation of materials science and interacting particles in curved space.


title: Inclusive prompt photon production in electron-nucleus scattering at small x

abstract: We compute the differential cross-section for inclusive prompt photon production in deeply inelastic scattering (DIS) of electrons on nuclei at small $x$ in the framework of the Color Glass Condensate (CGC) effective theory. The leading order (LO) computation in this framework resums leading logarithms in $x$ as well as power corrections to all orders in $Q_{s,A}^2/Q^2$, where $Q_{s,A}(x)$ is the nuclear saturation scale. This LO result is proportional to universal dipole and quadrupole Wilson line correlators in the nucleus. In the soft photon limit, the Low-Burnett-Kroll theorem allows us to recover existing results on inclusive DIS dijet production. The $k_{\perp}$ and collinearly factorized expressions for prompt photon production in DIS are also recovered in a leading twist approximation to our result. In the latter case, our result corresponds to the dominant next-to-leading order (NLO) perturbative QCD contribution at small $x$. We next discuss the computation of the NLO corrections to inclusive prompt photon production in the CGC framework. In particular, we emphasize the advantages for higher order computations in inclusive photon production, and for fully inclusive DIS, arising from the simple momentum space structure of the dressed quark and gluon "shock wave" propagators in the "wrong" light cone gauge $A^-=0$ for a nucleus moving with $P^{+}_{N} \rightarrow \infty$.


title: The Degrees of Freedom of Area Regge Calculus: Dynamics, Non-metricity, and Broken Diffeomorphisms

abstract: Discretization of general relativity is a promising route towards quantum gravity. Discrete geometries have a finite number of degrees of freedom and can mimic aspects of quantum geometry. However, selection of the correct discrete freedoms and description of their dynamics has remained a challenging problem. We explore classical area Regge calculus, an alternative to standard Regge calculus where instead of lengths, the areas of a simplicial discretization are fundamental. There are a number of surprises: though the equations of motion impose flatness we show that diffeomorphism symmetry is broken for a large class of area Regge geometries. This is due to degrees of freedom not available in the length calculus. In particular, an area discretization only imposes that the areas of glued simplicial faces agrees; their shapes need not be the same. We enumerate and characterize these non-metric, or `twisted', degrees of freedom and provide tools for understanding their dynamics. The non-metric degrees of freedom also lead to fewer invariances of the area Regge action---in comparison to the length action---under local changes of the triangulation (Pachner moves). This means that invariance properties can be used to classify the dynamics of spin foam models. Our results lay a promising foundation for understanding the dynamics of the non-metric degrees of freedom in loop quantum gravity and spin foams.


title: Pion and kaon valence-quark quasiparton distributions

abstract: Algebraic Ansaetze for the Poincar\'e-covariant Bethe-Salpeter wave functions of the pion and kaon are used to calculate their light-front wave functions (LFWFs), parton distribution amplitudes (PDAs), quasi-PDAs (qPDAs), valence parton distribution functions (PDFs), and quasi-PDFs (qPDFs). The LFWFs are broad, concave functions; and the scale of flavour-symmetry violation in the kaon is roughly 15%, being set by the ratio of emergent masses in the $s$-and $u$-quark sectors. qPDAs computed with longitudinal momentum $P_z =1.75\,$GeV provide a semiquantitatively accurate representation of the objective PDA; but even with $P_z=3\,$GeV, they cannot provide information about this amplitude's endpoint behaviour. On the valence-quark domain, similar outcomes characterise qPDFs. In this connection, however, the ratio of kaon-to-pion $u$-quark qPDFs is found to provide a good approximation to the true PDF ratio on $0.3\lesssim x \lesssim 0.8$, suggesting that with existing resources computations of ratios of quasi-parton-distributions can yield results that support empirical comparison.


title: How black holes shape globular clusters: Modeling NGC 3201

abstract: Numerical simulations have shown that black holes (BHs) can strongly influence the evolution and present-day observational properties of globular clusters (GCs). Using a Monte Carlo code, we construct GC models that match the Milky Way (MW) cluster NGC 3201, the first cluster in which a stellar-mass BH was identified through radial-velocity measurements. We predict that NGC 3201 contains $\gtrsim 200$ stellar-mass BHs. Furthermore, we explore the dynamical formation of main sequence-BH binaries and demonstrate that systems similar to the observed BH binary in NGC 3201 are produced naturally. Additionally, our models predict the existence of bright blue-straggler-BH binaries unique to core-collapsed clusters, which otherwise retain few BHs.


title: Quantum Cosmic No-Hair Theorem and Inflation

abstract: We consider implications of the quantum extension of the inflationary no hair theorem. We show that when the quantum state of inflation is picked to ensure the validity of the EFT of fluctuations, it takes only ${\cal O}(10)$ efolds of inflation to erase the effects of the initial distortions on the inflationary observables. Thus the Bunch-Davies vacuum is a very strong quantum attractor during inflation. We also consider bouncing universes, where the initial conditions seem to linger much longer and the quantum `balding' by evolution appears to be less efficient.


title: Involutive categories, colored $\ast$-operads and quantum field theory

abstract: Involutive category theory provides a flexible framework to describe involutive structures on algebraic objects, such as anti-linear involutions on complex vector spaces. Motivated by the prominent role of involutions in quantum (field) theory, we develop the involutive analogs of colored operads and their algebras, named colored $\ast$-operads and $\ast$-algebras. Central to the definition of colored $\ast$-operads is the involutive monoidal category of symmetric sequences, which we obtain from a general product-exponential $2$-adjunction whose right adjoint forms involutive functor categories. For $\ast$-algebras over $\ast$-operads we obtain involutive analogs of the usual change of color and operad adjunctions. As an application, we turn the colored operads for algebraic quantum field theory into colored $\ast$-operads. The simplest instance is the associative $\ast$-operad, whose $\ast$-algebras are unital and associative $\ast$-algebras.


title: Formal moduli problems and formal derived stacks

abstract: This paper presents a survey on formal moduli problems. It starts with an introduction to pointed formal moduli problems and a sketch of proof of a Theorem (independently proven by Lurie and Pridham) which gives a precise mathematical formulation for Drinfeld's derived deformation theory philosophy, which gives a correspondence between formal moduli problems and differential graded Lie algebras. The second part deals with Lurie's general theory of deformation contexts, which we present in a slightly different way than the original paper, emphasising the (more symmetric) notion of Koszul duality contexts and morphisms thereof. In the third part, we explain how to apply this machinery to the case of non-split formal moduli problems under a given derived affine scheme; this situation has been dealt with recently by Joost Nuiten, and requires to replace differential graded Lie algebras with differential graded Lie algebroids. In the last part, we globalize this to the more general setting of formal thickenings of derived stacks, and suggest an alternative approach to results of Gaitsgory and Rozenblyum.


title: Mass determination of the 1:3:5 near-resonant planets transiting GJ 9827 (K2-135)

abstract: Aims. GJ 9827 (K2-135) has recently been found to host a tightly packed system consisting of three transiting small planets whose orbital periods of 1.2, 3.6, and 6.2 days are near the 1:3:5 ratio. GJ 9827 hosts the nearest planetary system (d = $30.32\pm1.62$ pc) detected by Kepler and K2 . Its brightness (V = 10.35 mag) makes the star an ideal target for detailed studies of the properties of its planets. Results. We find that GJ 9827 b has a mass of $M_\mathrm{b}=3.74^{+0.50}_{-0.48}$ $M_\oplus$ and a radius of $R_\mathrm{b}=1.62^{+0.17}_{-0.16}$ $R_\oplus$, yielding a mean density of $\rho_\mathrm{b} = 4.81^{+1.97}_{-1.33}$ g cm$^{-3}$. GJ 9827 c has a mass of $M_\mathrm{c}=1.47^{+0.59}_{-0.58}$ $M_\oplus$, radius of $R_\mathrm{c}=1.27^{+0.13}_{-0.13}$ $R_\oplus$, and a mean density of $\rho_\mathrm{c}= 3.87^{+2.38}_{-1.71}$ g cm$^{-3}$. For GJ 9827 d we derive $M_\mathrm{d}=2.38^{+0.71}_{-0.69}$ $M_\oplus$, $R_\mathrm{d}=2.09^{+0.22}_{-0.21}$ $R_\oplus$, and $\rho_\mathrm{d}= 1.42^{+0.75}_{-0.52}$ g cm$^{-3}$. Conclusions. GJ 9827 is one of the few known transiting planetary systems for which the masses of all planets have been determined with a precision better than 30%. This system is particularly interesting because all three planets are close to the limit between super-Earths and mini-Neptunes. We also find that the planetary bulk compositions are compatible with a scenario where all three planets formed with similar core/atmosphere compositions, and we speculate that while GJ 9827 b and GJ 9827 c lost their atmospheric envelopes, GJ 9827 d maintained its atmosphere, owing to the much lower stellar irradiation. This makes GJ 9827 one of the very few systems where the dynamical evolution and the atmospheric escape can be studied in detail for all planets, helping us to understand how compact systems form and evolve.


title: Constructing exact representations of quantum many-body systems with deep neural networks

abstract: We develop a constructive approach to generate artificial neural networks representing the exact ground states of a large class of many-body lattice Hamiltonians. It is based on the deep Boltzmann machine architecture, in which two layers of hidden neurons mediate quantum correlations among physical degrees of freedom in the visible layer. The approach reproduces the exact imaginary-time Hamiltonian evolution, and is completely deterministic. In turn, compact and exact network representations for the ground states are obtained without stochastic optimization of the network parameters. The number of neurons grows linearly with the system size and total imaginary time, respectively. Physical quantities can be measured by sampling configurations of both physical and neuron degrees of freedom. We provide specific examples for the transverse-field Ising and Heisenberg models by implementing efficient sampling. As a compact, classical representation for many-body quantum systems, our approach is an alternative to the standard path integral, and it is potentially useful also to systematically improve on numerical approaches based on the restricted Boltzmann machine architecture.

categories: cond-mat.dis-nn cond-mat.stat-mech cond-mat.str-el physics.comp-ph quant-ph


title: Biorthogonal vectors, sesquilinear forms and some physical operators

abstract: Continuing the analysis undertaken in previous articles, we discuss some features of non-self-adjoint operators and sesquilinear forms which are defined starting from two biorthogonal families of vectors, like the so-called generalized Riesz systems, enjoying certain properties. In particular we discuss what happens when they forms two $\D$-quasi bases.


title: Evolution of the Stellar Mass--Metallicity Relation - I: Galaxies in the z~0.4 Cluster Cl0024

abstract: We present the stellar mass-stellar metallicity relationship (MZR) in the Cl0024+1654 galaxy cluster at z~0.4 using full spectrum stellar population synthesis modeling of individual quiescent galaxies. The lower limit of our stellar mass range is $M_*=10^{9.7}M_\odot$, the lowest galaxy mass at which individual stellar metallicity has been measured beyond the local universe. We report a detection of an evolution of the stellar MZR with observed redshift at $0.037\pm0.007$ dex per Gyr, consistent with the predictions from hydrodynamical simulations. Additionally, we find that the evolution of the stellar MZR with observed redshift can be explained by an evolution of the stellar MZR with their formation time, i.e., when the single stellar population (SSP)-equivalent ages of galaxies are taken into account. This behavior is consistent with stars forming out of gas that also has an MZR with a normalization that decreases with redshift. Lastly, we find that over the observed mass range, the MZR can be described by a linear function with a shallow slope, ($[Fe/H] \propto (0.16 \pm 0.03) \log M_*$). The slope suggests that galaxy feedback, in terms of mass-loading factor, might be mass-independent over the observed mass and redshift range.


title: Constructing an Explicit AdS/CFT Correspondence with Cartan Geometry

abstract: An explicit AdS/CFT correspondence is shown for the Lie group $SO(4,2)$. The Lie symmetry structures allow for the construction of two physical theories through the tools of Cartan geometry. One is a gravitational theory that has anti-de Sitter symmetry. The other is also a gravitational theory but is conformally symmetric and lives on 8-dimensional biconformal space. These "extra" four dimensions have the degrees of freedom used to construct a Yang-Mills theory. The two theories, based on AdS or conformal symmetry, have a natural correspondence in the context of their Lie algebras alone where neither SUSY, nor holography, is necessary.


title: Radial velocities of RR Lyrae stars in and around NGC 6441

abstract: Detailed elemental abundance patterns of metal-poor ([Fe/H] ~ -1~dex) stars in the Galactic bulge indicate that a number of them are consistent with globular cluster (GC) stars and may be former members of dissolved GCs. This would indicate that a few per cent of the Galactic bulge was built up from destruction and/or evaporation of globular clusters. Here an attempt is made to identify such presumptive stripped stars originating from the massive, inner Galaxy globular cluster NGC~6441 using its rich RR Lyrae variable star (RRL) population. We present radial velocities of forty RRLs centered on the globular cluster NGC~6441. All of the 13 RRLs observed within the cluster tidal radius have velocities consistent with cluster membership, with an average radial velocity of 24 +- 5~km/s and a star-to-star scatter of 11~km/s. This includes two new RRLs that were previously not associated with the cluster. Eight RRLs with radial velocities consistent with cluster membership but up to three time the distance from the tidal radius are also reported. These potential extra-tidal RRLs also have exceptionally long periods, which is a curious characteristic of the NGC~6441 RRL population that hosts RRLs with periods longer than seen anywhere else in the Milky Way. As expected of stripped cluster stars, most are inline with the cluster's orbit. Therefore, either the tidal radius of NGC~6441 is underestimated and/or we are seeing dissolving cluster stars stemming from NGC~6441 that are building up the old spheroidal bulge.


title: OGLE-2016-BLG-1266: A Probable Brown-Dwarf/Planet Binary at the Deuterium Fusion Limit

abstract: We report the discovery, via the microlensing method, of a new very-low-mass binary system. By combining measurements from Earth and from the Spitzer telescope in Earth-trailing orbit, we are able to measure the microlensing parallax of the event, and find that the lens likely consists of an $(12.0 \pm 0.6) M_{\rm J}$ + $(15.7 \pm 1.5) M_{\rm J}$ super-Jupiter / brown-dwarf pair. The binary is located at a distance of $(3.08 \pm 0.18)$ kpc in the Galactic Plane, and the components have a projected separation of $(0.43 \pm 0.03)$ AU. Two alternative solutions with much lower likelihoods are also discussed, an 8- and 6-$M_{\rm J}$ model and a 90- and 70-$M_{\rm J}$ model. Although disfavored at the 3-$\sigma$ and 5-$\sigma$ levels, these alternatives cannot be rejected entirely. We show how the more-massive of these models could be tested with future direct imaging.


title: Reinforcement and Imitation Learning for Diverse Visuomotor Skills

abstract: We propose a model-free deep reinforcement learning method that leverages a small amount of demonstration data to assist a reinforcement learning agent. We apply this approach to robotic manipulation tasks and train end-to-end visuomotor policies that map directly from RGB camera inputs to joint velocities. We demonstrate that our approach can solve a wide variety of visuomotor tasks, for which engineering a scripted controller would be laborious. In experiments, our reinforcement and imitation agent achieves significantly better performances than agents trained with reinforcement learning or imitation learning alone. We also illustrate that these policies, trained with large visual and dynamics variations, can achieve preliminary successes in zero-shot sim2real transfer. A brief visual description of this work can be viewed in https://youtu.be/EDl8SQUNjj0


title: Conjugate Bayes for probit regression via unified skew-normal distributions

abstract: Regression models for dichotomous data are ubiquitous in statistics. Besides being useful for inference on binary responses, these methods serve also as building blocks in more complex formulations, such as density regression, nonparametric classification and graphical models. Within the Bayesian framework, inference proceeds by updating the priors for the coefficients, typically set to be Gaussians, with the likelihood induced by probit or logit regressions for the responses. In this updating, the apparent absence of a tractable posterior has motivated a variety of computational methods, including Markov Chain Monte Carlo routines and algorithms which approximate the posterior. Despite being routinely implemented, Markov Chain Monte Carlo strategies face mixing or time-inefficiency issues in large p and small n studies, whereas approximate routines fail to capture the skewness typically observed in the posterior. This article proves that the posterior distribution for the probit coefficients has a unified skew-normal kernel, under Gaussian priors. Such a novel result allows efficient Bayesian inference for a wide class of applications, especially in large p and small-to-moderate n studies where state-of-the-art computational methods face notable issues. These advances are outlined in a genetic study, and further motivate the development of a wider class of conjugate priors for probit models along with methods to obtain independent and identically distributed samples from the unified skew-normal posterior.


title: Design of iMacros-based Data Crawler and the Behavioral Analysis of Facebook Users

abstract: Obtaining the desired dataset is still a prime challenge faced by researchers while analyzing Online Social Network (OSN) sites. Application Programming Interfaces (APIs) provided by OSN service providers for retrieving data impose several unavoidable restrictions which make it difficult to get a desirable dataset. In this paper, we present an iMacros technology-based data crawler called IMcrawler, capable of collecting every piece of information which is accessible through a browser from the Facebook website within the legal framework which permits access to publicly shared user content on OSNs. The proposed crawler addresses most of the challenges allied with web data extraction approaches and most of the APIs provided by OSN service providers. Two broad sections have been extracted from Facebook user profiles, namely, Personal Information and Wall Activities. The present work is the first attempt towards providing the detailed description of crawler design for the Facebook website.


title: A Robust Real-Time Automatic License Plate Recognition Based on the YOLO Detector

abstract: Automatic License Plate Recognition (ALPR) has been a frequent topic of research due to many practical applications. However, many of the current solutions are still not robust in real-world situations, commonly depending on many constraints. This paper presents a robust and efficient ALPR system based on the state-of-the-art YOLO object detector. The Convolutional Neural Networks (CNNs) are trained and fine-tuned for each ALPR stage so that they are robust under different conditions (e.g., variations in camera, lighting, and background). Specially for character segmentation and recognition, we design a two-stage approach employing simple data augmentation tricks such as inverted License Plates (LPs) and flipped characters. The resulting ALPR approach achieved impressive results in two datasets. First, in the SSIG dataset, composed of 2,000 frames from 101 vehicle videos, our system achieved a recognition rate of 93.53% and 47 Frames Per Second (FPS), performing better than both Sighthound and OpenALPR commercial systems (89.80% and 93.03%, respectively) and considerably outperforming previous results (81.80%). Second, targeting a more realistic scenario, we introduce a larger public dataset, called UFPR-ALPR dataset, designed to ALPR. This dataset contains 150 videos and 4,500 frames captured when both camera and vehicles are moving and also contains different types of vehicles (cars, motorcycles, buses and trucks). In our proposed dataset, the trial versions of commercial systems achieved recognition rates below 70%. On the other hand, our system performed better, with recognition rate of 78.33% and 35 FPS.


title: Shampoo: Preconditioned Stochastic Tensor Optimization

abstract: Preconditioned gradient methods are among the most general and powerful tools in optimization. However, preconditioning requires storing and manipulating prohibitively large matrices. We describe and analyze a new structure-aware preconditioning algorithm, called Shampoo, for stochastic optimization over tensor spaces. Shampoo maintains a set of preconditioning matrices, each of which operates on a single dimension, contracting over the remaining dimensions. We establish convergence guarantees in the stochastic convex setting, the proof of which builds upon matrix trace inequalities. Our experiments with state-of-the-art deep learning models show that Shampoo is capable of converging considerably faster than commonly used optimizers. Although it involves a more complex update rule, Shampoo's runtime per step is comparable to that of simple gradient methods such as SGD, AdaGrad, and Adam.


title: Transport Properties of Near Surface InAs Two-dimensional Heterostructures

abstract: Two-dimensional electron systems (2DESs) confined to the surface of narrowband semiconductors have attracted great interest since they can easily integrate with superconductivity (or ferromagnetism) enabling new possibilities in hybrid device architectures and study of exotic states in proximity of superconductors. In this work, we study indium arsenide heterostructures where combination of clean interface with superconductivity, high mobility and spin-orbit coupling can be achieved. The weak antilocalization measurements indicate presence of strong spin-orbit coupling at high densities. We study the magnetotransport as a function of top barrier and density and report clear observation of integer quantum Hall states. We report improved electron mobility reaching up to 44,000 cm$^{2}$/Vs in undoped heterstructures and well developed integer quantum Hall states starting as low as 2.5~T.


title: Porous medium equation with a blow-up nonlinearity and a non-decreasing constraint

abstract: The final goal of this paper is to prove existence of local (strong) solutions to a (fully nonlinear) porous medium equation with blow-up term and nondecreasing constraint. To this end, the equation, arising in the context of Damage Mechanics, is reformulated as a mixed form of two different types of doubly nonlinear evolution equations. Global (in time) solutions to some approximate problems are constructed by performing a time discretization argument and by taking advantage of energy techniques based on specific structures of the equation. Moreover, a variational comparison principle for (possibly non-unique) approximate solutions is established and it also enables us to obtain a local solution as a limit of approximate ones.


title: Finite volume mass gap and free energy of the SU(N)xSU(N) chiral sigma model

abstract: We compute the free energy in the presence of a chemical potential coupled to a conserved charge in the effective SU(N)xSU(N) scalar field theory to third order for asymmetric volumes in general d-dimensions, using dimensional regularization. We also compute the mass gap in a finite box with periodic boundary conditions.


title: General volumes in the Orlicz-Brunn-Minkowski theory and a related Minkowski Problem I

abstract: The general volume of a star body, a notion that includes the usual volume, the $q$th dual volumes, and many previous types of dual mixed volumes, is introduced. A corresponding new general dual Orlicz curvature measure is defined that specializes to the $(p,q)$-dual curvature measures introduced recently by Lutwak, Yang, and Zhang. General variational formulas are established for the general volume of two types of Orlicz linear combinations. One of these is applied to the Minkowski problem for the new general dual Orlicz curvature measure, giving in particular a solution to the Minkowski problem posed by Lutwak, Yang, and Zhang for the $(p,q)$-dual curvature measures when $p>0$ and $q<0$. A dual Orlicz-Brunn-Minkowski inequality for general volumes is obtained, as well as dual Orlicz-Minkowski-type inequalities and uniqueness results for star bodies. Finally, a very general Minkowski-type inequality, involving two Orlicz functions, two convex bodies, and a star body, is proved, that includes as special cases several others in the literature, in particular one due to Lutwak, Yang, and Zhang for the $(p,q)$-mixed volume.


title: The gravitational field of static p-branes in linearized ghost-free gravity

abstract: We study the gravitational field of static p-branes in D-dimensional Minkowski space in the framework of linearized ghost-free (GF) gravity. The concrete models of GF gravity we consider are parametrized by the non-local form factors $\exp(-\Box/\mu^2)$ and $\exp(\Box^2/\mu^4)$, where $\mu^{-1}$ is the scale of non-locality. We show that the singular behavior of the gravitational field of p-branes in General Relativity is cured by short-range modifications introduced by the non-localities, and we derive exact expressions of the regularized gravitational fields, whose geometry can be written as a warped metric. For large distances compared to the scale of non-locality, $\mu r\rightarrow\infty$, our solutions approach those found in linearized General Relativity.


title: Optimal Investment Decision Under Switching regimes of Subsidy Support

abstract: We address the problem of making a managerial decision when the investment project is subsidized, which results in the resolution of an infinite-horizon optimal stopping problem of a switching diffusion driven by either an homogeneous or an inhomogeneous continuous-time Markov chain. We provide a characterization of the value function (and optimal strategy) of the optimal stopping problem. On the one hand, broadly, we can prove that the value function is the unique viscosity solution to a system of HJB equations. On the other hand, when the Markov chain is homogeneous and the switching diffusion is one-dimensional, we obtain stronger results: the value function is the difference between two convex functions.


title: i3PosNet: Instrument Pose Estimation from X-Ray in temporal bone surgery

abstract: Purpose: Accurate estimation of the position and orientation (pose) of surgical instruments is crucial for delicate minimally invasive temporal bone surgery. Current techniques lack in accuracy and/or line-of-sight constraints (conventional tracking systems) or expose the patient to prohibitive ionizing radiation (intra-operative CT). A possible solution is to capture the instrument with a c-arm at irregular intervals and recover the pose from the image. Methods: i3PosNet infers the position and orientation of instruments from images using a pose estimation network. Said framework considers localized patches and outputs pseudo-landmarks. The pose is reconstructed from pseudo-landmarks by geometric considerations. Results: We show i3PosNet reaches errors less than 0.05mm. It outperforms conventional image registration-based approaches reducing average and maximum errors by at least two thirds. i3PosNet trained on synthetic images generalizes to real x-rays without any further adaptation. Conclusion: The translation of Deep Learning based methods to surgical applications is difficult, because large representative datasets for training and testing are not available. This work empirically shows sub-millimeter pose estimation trained solely based on synthetic training data.


title: Phases of quantum dimers from ensembles of classical stochastic trajectories

abstract: We study the connection between the phase behaviour of quantum dimers and the dynamics of classical stochastic dimers. At the so-called Rokhsar-Kivelson (RK) point a quantum dimer Hamiltonian is equivalent to the Markov generator of the dynamics of classical dimers. A less well understood fact is that away from the RK point the quantum-classical connection persists: in this case the Hamiltonian corresponds to a non-stochastic "tilted" operator that encodes the statistics of time-integrated observables of the classical stochastic problem. This implies a direct relation between the phase behaviour of quantum dimers and properties of ensembles of stochastic trajectories of classical dimers. We make these ideas concrete by studying fully packed dimers on the square lattice. Using transition path sampling - supplemented by trajectory umbrella sampling - we obtain the large deviation statistics of dynamical activity in the classical problem, and show the correspondence between the phase behaviour of the classical and quantum systems. The transition at the RK point between quantum phases of distinct order corresponds, in the classical case, to a trajectory phase transition between active and inactive dynamical phases. Furthermore, from the structure of stochastic trajectories in the active dynamical phase we infer that the ground state of quantum dimers has columnar order to one side of the RK point. We discuss how these results relate to those from quantum Monte Carlo, and how our approach may generalise to other problems.


title: Stochastic inflation with quantum and thermal noise

abstract: We add a thermal noise to Starobinsky equation of slow roll inflation. We calculate the number of e-folds of the stochastic system. The power spectrum and the spectral index are evaluated from the fluctuations of e-folds using an expansion in the quantum and thermal noise terms.


title: Near-Linear Time Local Polynomial Nonparametric Estimation with Box Kernels

abstract: Local polynomial regression (Fan and Gijbels 1996) is an important class of methods for nonparametric density estimation and regression problems. However, straightforward implementation of local polynomial regression has quadratic time complexity which hinders its applicability in large-scale data analysis. In this paper, we significantly accelerate the computation of local polynomial estimates by novel applications of multi-dimensional binary indexed trees (Fenwick 1994). Both time and space complexity of our proposed algorithm is nearly linear in the number of input data points. Simulation results confirm the efficiency and effectiveness of our proposed approach.


title: An adaptive timestepping methodology for particle advance in coupled CFD-DEM simulations

abstract: An adpative integration technique for time advancement of particle motion in the context of coupled computational fluid dynamics (CFD) - discrete element method (DEM) simulations is presented in this work. CFD-DEM models provide an accurate description of multiphase physical systems where a granular phase exists in an underlying continuous medium. The time integration of the granular phase in these simulations present unique computational challenges due to large variations in time scales associated with particle collisions. The algorithm presented in this work uses a local time stepping approach to resolve collisional time scales for only a subset of particles that are in close proximity to potential collision partners, thereby resulting in substantial reduction of computational cost. This approach is observed to be 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.


title: Lossy Compression of Decimated Gaussian Random Walks

abstract: We consider the problem of estimating a Gaussian random walk from a lossy compression of its decimated version. Hence, the encoder operates on the decimated random walk, and the decoder estimates the original random walk from its encoded version under a mean squared error (MSE) criterion. It is well-known that the minimal distortion in this problem is attained by an estimate-and-compress (EC) source coding strategy, in which the encoder first estimates the original random walk and then compresses this estimate subject to the bit constraint. In this work, we derive a closed-form expression for this minimal distortion as a function of the bitrate and the decimation factor. Next, we consider a compress-and-estimate (CE) source coding scheme, in which the encoder first compresses the decimated sequence subject to an MSE criterion (with respect to the decimated sequence), and the original random walk is estimated only at the decoder. We evaluate the distortion under CE in a closed form and show that there exists a nonzero gap between the distortion under the two schemes. This difference in performance illustrates the importance of having the decimation factor at the encoder.


title: The overlooked role of stellar variability in LMC intermediate-age clusters

abstract: Broad, extended main sequence turnoffs seen in the majority of the intermediate-age (1-3 Gyr) LMC star clusters, have been interpreted as the result of an extended star formation history and/or the effect of extreme stellar rotation. A more fundamental explanation may be given by stellar variability. For clusters in these age range, the instability strip crosses the upper main sequence producing a number of variable stars (known as Delta Scuti) which, if nor properly taken into account, could appear as an extended turnoff. First results of a variability program in the LMC cluster NGC 1846 reveals a sizeable number of this type of variables, although still too low to produce a meaningful broadening, with the caveat that the true variable content of the center of this and other clusters in the LMC will only be revealed with a dedicated HST program.


title: A graph-theoretic framework for algorithmic design of experiments

abstract: In this paper, we demonstrate that considering experiments in a graph-theoretic manner allows us to exploit automorphisms of the graph to reduce the number of evaluations of candidate designs for those experiments, and thus find optimal designs faster. We show that the use of automorphisms for reducing the number of evaluations required of an optimality criterion function is effective on designs where experimental units have a network structure. Moreover, we show that we can take block designs with no apparent network structure, such as one-way blocked experiments, row-column experiments, and crossover designs, and add block nodes to induce a network structure. Considering automorphisms can thus reduce the amount of time it takes to find optimal designs for a wide class of experiments.


title: Data-dependent PAC-Bayes priors via differential privacy

abstract: The Probably Approximately Correct (PAC) Bayes framework (McAllester, 1999) can incorporate knowledge about the learning algorithm and (data) distribution through the use of distribution-dependent priors, yielding tighter generalization bounds on data-dependent posteriors. Using this flexibility, however, is difficult, especially when the data distribution is presumed to be unknown. We show how an {\epsilon}-differentially private data-dependent prior yields a valid PAC-Bayes bound, and then show how non-private mechanisms for choosing priors can also yield generalization bounds. As an application of this result, we show that a Gaussian prior mean chosen via stochastic gradient Langevin dynamics (SGLD; Welling and Teh, 2011) leads to a valid PAC-Bayes bound given control of the 2-Wasserstein distance to an {\epsilon}-differentially private stationary distribution. We study our data-dependent bounds empirically, and show that they can be nonvacuous even when other distribution-dependent bounds are vacuous.


title: Magnetic field dependence of electronic properties of MoS$_2$ quantum dots with different edges

abstract: Using the tight-binding approach, we investigate the energy spectrum of square, triangular and hexagonal MoS$_2$ quantum dots (QDs) in the presence of a perpendicular magnetic field. Novel edge states emerge in MoS$_2$ QDs, which are distributed over the whole edge which we call ring states. The ring states are robust in the presence of spin-orbit coupling (SOC). The corresponding energy levels of the ring states oscillate as function of the perpendicular magnetic field which are related to Aharonov-Bohm oscillations. Oscillations in the magnetic field dependence of the energy levels and the peaks in the magneto-optical spectrum emerge (disappear) as the ring states are formed (collapsed). The period and the amplitude of the oscillation decreases with the size of the MoS$_2$ QDs.


title: A partial correlation vine based approach for modeling and forecasting multivariate volatility time-series

abstract: A novel approach for dynamic modeling and forecasting of realized covariance matrices is proposed. Realized variances and realized correlation matrices are jointly estimated. The one-to-one relationship between a positive definite correlation matrix and its associated set of partial correlations corresponding to any vine specification is used for data transformation. The model components therefore are realized variances as well as realized standard and partial correlations corresponding to a daily log-return series. As such, they have a clear practical interpretation. A method to select a regular vine structure, which allows for parsimonious time-series and dependence modeling of the model components, is introduced. Being algebraically independent the latter do not underlie any algebraic constraint. The proposed model approach is outlined in detail and motivated along with a real data example on six highly liquid stocks. The forecasting performance is evaluated both with respect to statistical precision and in the context of portfolio optimization. Comparisons with Cholesky decomposition based benchmark models support the excellent prediction ability of the proposed model approach.