title: CytonRL: an Efficient Reinforcement Learning Open-source Toolkit Implemented in C++

abstract: This paper presents an open-source enforcement learning toolkit named CytonRL (https://github.com/arthurxlw/cytonRL). The toolkit implements four recent advanced deep Q-learning algorithms from scratch using C++ and NVIDIA's GPU-accelerated libraries. The code is simple and elegant, owing to an open-source general-purpose neural network library named CytonLib. Benchmark shows that the toolkit achieves competitive performances on the popular Atari game of Breakout.


title: Two-field cosmological phase transitions and gravitational waves in the singlet Majoron model

abstract: In the singlet Majoron model, we study cosmological phase transitions (PTs) and their resulting gravitational waves (GWs), in the two-field phase space, without freezing any of the field directions. We first calculate the effective potential, at one loop and at finite temperature, of the Standard Model Higgs doublet together with one extra Higgs singlet. We make use of the public available Python package CosmoTransitions to simulate the two-dimensional cosmological PTs and evaluate the gravitational waves generated by first-order PTs. With the full 2D simulation, we are able not only to confirm the PTs' properties previously discussed in the literature, but also we find new patterns, such as strong first-order PTs tunneling from a vacuum located on one axis to another vacuum located on the second axis. The two-field phase space analysis presents a richer panel of cosmological PT patterns compared to analysis with a single-field approximation. The PTGW amplitudes turn out to be out of the reach for the space-borne gravitational wave interferometers such as LISA, DECIGO, BBO, TAIJI and TianQin when constraints from colliders physics are taken into account.


title: Explicit Construction of the Voronoi and Delaunay Cells of W(An) and W(Dn) Lattices and Their Facets

abstract: Voronoi and Delaunay (Delone) cells of the root and weight lattices of the Coxeter-Weyl groups W(an) and W(dn) are constructed. The face centered cubic (fcc) and body centered cubic (bcc)lattices are obtained in this context. Basic definitions are introduced such as parallelotope, fundamental simplex, contact polytope, root polytope, Voronoi cell, Delone cell, n-simplex, n-octahedron (cross polytope), n-cube and n-hemicube and their volumes are calculated. Voronoi cell of the root lattice is constructed as the dual of the root polytope which turns out to be the union of Delone cells. It is shown that the Delone cells centered at the origin of the root lattice An are the polytopes of the fundamental weights w1, w2, ..., wn and the Delone cells of the root lattice Dn are the polytopes obtained from the weights w1, wn-1, wn. A simple mechanism explains the tessellation of the root lattice by Delone cells. We prove that the (n-1)-facet of the Voronoi cell of the root lattice An is (n-1)-dimensional rhombohedron and similarly the (n-1) -facet of the Voronoi cell of the root lattice Dn is a dipyramid with a base of (n-2)-cube. Volume of the Voronoi cell is calculated via its (n-1) -facet which in turn can be obtained from the fundamental simplex. Tessellations of the root lattice with the Voronoi and Delone cells are explained by giving examples from lower dimensions. Similar considerations are also worked out for the weight lattices An* and Dn*.


title: Walk-Steered Convolution for Graph Classification

abstract: Graph classification is a fundamental but challenging issue for numerous real-world applications. Despite recent great progress in image/video classification, convolutional neural networks (CNNs) cannot yet cater to graphs well because of graphical non-Euclidean topology. In this work, we propose a walk-steered convolutional (WSC) network to assemble the essential success of standard convolutional neural networks as well as the powerful representation ability of random walk. Instead of deterministic neighbor searching used in previous graphical CNNs, we construct multi-scale walk fields (a.k.a. local receptive fields) with random walk paths to depict subgraph structures and advocate graph scalability. To express the internal variations of a walk field, Gaussian mixture models are introduced to encode principal components of walk paths therein. As an analogy to a standard convolution kernel on image, Gaussian models implicitly coordinate those unordered vertices/nodes and edges in a local receptive field after projecting to the gradient space of Gaussian parameters. We further stack graph coarsening upon Gaussian encoding by using dynamic clustering, such that high-level semantics of graph can be well learned like the conventional pooling on image. The experimental results on several public datasets demonstrate the superiority of our proposed WSC method over many state-of-the-arts for graph classification.


title: Demonstration of sensitivity increase in mercury free-spin-precession magnetometers due to laser-based readout for neutron electric dipole moment searches

abstract: We report on a laser based $^{199}$Hg co-magnetometer deployed in an experiment searching for a permanent electric dipole moment of the neutron. We demonstrate a more than five times increased signal to-noise-ratio in a direct comparison measurement with its $^{204}$Hg discharge bulb-based predecessor. An improved data model for the extraction of important system parameters such as the degrees of absorption and polarization is derived. Laser- and lamp-based data-sets can be consistently described by the improved model which permits to compare measurements using the two different light sources and to explain the increase in magnetometer performance. The laser-based magnetometer satisfies the magnetic field sensitivity requirements for the next generation nEDM experiments.

categories: physics.ins-det hep-ex nucl-ex physics.atom-ph


title: BigDL: A Distributed Deep Learning Framework for Big Data

abstract: This paper presents BigDL (a distributed deep learning framework for Apache Spark), which has been used by a variety of users in the industry for building deep learning applications on production big data platforms. It allows deep learning applications to run on the Apache Hadoop/Spark cluster so as to directly process the production data, and as a part of the end-to-end data analysis pipeline for deployment and management. Unlike existing deep learning frameworks, BigDL implements distributed, data parallel training directly on top of the functional compute model (with copy-on-write and coarse-grained operations) of Spark. We also share real-world experience and "war stories" of users that have adopted BigDL to address their challenges(i.e., how to easily build end-to-end data analysis and deep learning pipelines for their production data).


title: MOND from a brane-world picture

abstract: I describe a heuristic model where MOND dynamics emerge in a universe viewed as a nearly spherical brane embedded in a higher-dimensional flat space. The brane, described by $\xi(\Omega)$, is of density $\sigma$ ($\xi$ and $\Omega$ are the radial and angular coordinates in the embedding space). The brane and matter -- confined to the brane and of density $\rho(\Omega)\ll\sigma$ -- are coupled to a potential $\varepsilon(\xi)$. I restrict myself to shallow perturbations, $\xi(\Omega)=\ell_0+\zeta(\Omega)$, $|\zeta|\ll\ell_0$. A balanced brane implies $\hat a_0\equiv\varepsilon'(\ell_0)\sim T/\sigma\ell_0$, $T$ is the brane tension, yielding for the velocity of small brane perturbations $c^2\sim T/\sigma\sim \ell_0\hat a_0$. But, $\hat a_0$ plays the role of the MOND acceleration constant in local gravitational dynamics; so $\hat a_0\sim c^2/\ell_0$. What we, in the brane, perceive as the gravitational potential is $\phi\equiv\varepsilon[\xi(\Omega)]\approx \phi_0+\hat a_0\zeta$. Aspects of MOND that may emerge naturally as geometrical properties are: a. The special role of acceleration in MOND, and why it is an acceleration, $a_0$, that marks the transition from the standard dynamics much above $a_0$ to scale-invariant dynamics much below $a_0$. b. The intriguing connection of $a_0$ with cosmology. c. The Newtonian limit corresponds to local departure $|\zeta|\ll\ell_0$; i.e., $\phi-\phi_0\sim a_0\zeta\ll a_0\ell_0\sim c^2$ - whereas relativity enters when $|\zeta|\not\ll\ell_0$. The model also opens new vistas for extension, e.g., it points to possible dependence of $a_0$ on $\phi$, and to $a_0$ losing its status and meaning altogether in the relativistic regime. The required global balance of the brane might solve the `old' cosmological-constant problem. I discuss possible connections with the nearly-de-Sitter nature of our Universe. (Abridged.)


title: Spin transport in long-range interacting one-dimensional chain

abstract: We numerically study spin transport and nonequilibrium spin-density profiles in a clean one-dimensional spin-chain with long-range interactions, decaying as a power-law,$r^{-\alpha}$ with distance. We find two distinct regimes of transport: for $\alpha<1/2$, spin excitations relax instantaneously in the thermodynamic limit, and for $\alpha>1/2$, spin transport combines both diffusive and superdiffusive features. We show that while for $\alpha>3/2$ the spin diffusion coefficient is finite, transport in the system is never strictly diffusive, contrary to corresponding classical systems.


title: New physics vs new paradigms: distinguishing CPT violation from NSI

abstract: Our way of describing Nature is based on local relativistic quantum field theories, and then CPT symmetry, a natural consequence of Lorentz invariance, locality and hermiticity of the Hamiltonian, is one of the few if not the only prediction that all of them share. Therefore, testing CPT invariance does not test a particular model but the whole paradigm. Current and future long baseline experiments will assess the status of CPT in the neutrino sector at an unprecedented level and thus its distinction from similar experimental signatures arising from non-standard interactions is imperative. Whether the whole paradigm is at stake or just the standard model of neutrinos crucially depends on that.


title: Astrophysical radio background cannot explain the EDGES 21-cm signal: constraints from cooling of non-thermal electrons

abstract: Recently the EDGES experiment has claimed the detection of an absorption feature centered at 78 MHz. When interpreted as a signature of cosmic dawn, this feature appears at the correct wavelength (corresponding to a redshift range of $z\approx15-20$) but is larger by at least a factor of two in amplitude compared to the standard 21-cm models. One way to explain the excess radio absorption is by the enhancement of the diffuse radio background at $\nu = 1.42$ GHz ($\lambda=21$ cm) in the rest frame of the absorbing neutral hydrogen. Astrophysical scenarios, based on the acceleration of relativistic electrons by accretion on to supermassive black holes (SMBHs) and by supernovae (SN) from first stars, have been proposed to produce the enhanced radio background via synchrotron emission. In this Letter we show that either the synchrotron or the inverse-Compton (IC) cooling time for such electrons is at least three orders of magnitude shorter than the duration of the EDGES signal centered at $z \approx 17$, irrespective of the magnetic field strength. The synchrotron radio emission at 1.42 GHz due to rapidly cooling electrons is $\sim 10^3$ times smaller than the non-cooling estimate. Thus astrophysical scenarios for excess radio background proposed to explain the EDGES signal appear very unlikely.


title: $B$-decay anomalies and scalar leptoquarks in unified Pati-Salam models from noncommutative geometry

abstract: Motivated by possible scalar-leptoquark explanations of the recently reported $B$-decay anomalies, we investigate whether the required leptoquarks can be accommodated within models based on noncommutative geometry (NCG). The models considered have the gauge structure of Pati-Salam models, $SU(4)\times SU(2)_L\times SU(2)_R$, with gauge coupling unification at a single scale. In one of the models, we find a unique scalar leptoquark with quantum numbers $(3,1,-\frac{1}{3})_{321}$, originating from a complex multiplet $(6,1,1)_{422}$, which can potentially explain the $B$-decay anomalies if its mass is on the order of a few TeV. The unification of couplings can be realized with the inclusion of a single step of intermediate symmetry breaking. The scalar leptoquark under consideration does not contribute to proton decay due to the absence of diquark couplings, as dictated by the underlying noncommutative geometry.


title: Dihedral symmetry in $SU(N)$ Yang-Mills theory

abstract: We point out that charge conjugation and coordinate reflection symmetries do not commute with the center symmetry of $SU(N)$ YM theory when $N>2$. As a result, for generic values of the $\theta$ angle, the group of discrete zero-form symmetries of YM theory on e.g. the spacetime manifold $\mathbb{R}^3\times S^1$ includes the dihedral group $D_{2N}$ which is non-Abelian for $N>2$. At $\theta = \pi$, the non-Abelian factor in the symmetry group is enhanced to $D_{4N}$ due to discrete 't Hooft anomaly considerations. We illustrate these results in YM theory as well as in a simple quantum mechanical model, where we study representation theory as a function of $\theta$ angle.


title: On the existence of the Blandford-Znajek monopole for a slowly rotating Kerr black hole

abstract: The Blandford-Znajek monopole is a conjectured solution of force-free electrodynamics in the vicinity of a slowly rotating Kerr black hole, supposedly defined as a perturbation in small angular momentum. It is used to argue for the extraction of energy from rotating black holes by the Blandford-Znajek process. We set up a careful analysis of the perturbative definition of the Blandford-Znajek monopole, showing in particular that the regime in which it is defined allows to use the technique of matched asymptotic expansions. Our conclusion is that the Blandford-Znajek monopole, as it is defined, is not consistent with demanding physically reasonable boundary conditions far away from the event horizon. This puts into question the existence of the Blandford-Znajek monopole, at least in the limit of slow rotation of the Kerr black hole.


title: On the rates of type Ia supernovae originating from white dwarf collisions in quadruple star systems

abstract: We consider the evolution of stellar hierarchical quadruple systems in the 2+2 (two binaries orbiting each other's barycentre) and 3+1 (triple orbited by a fourth star) configurations. In our simulations, we take into account the effects of secular dynamical evolution, stellar evolution, tidal evolution and encounters with passing stars. We focus on type Ia supernovae (SNe Ia) driven by collisions of carbon-oxygen (CO) white dwarfs (WDs). Such collisions can arise from several channels: (1) collisions due to extremely high eccentricities induced by secular evolution, (2) collisions following a dynamical instability of the system, and (3) collisions driven by semisecular evolution. The systems considered here have initially wide inner orbits, with initial semilatus recti larger than 12 AU, implying no interaction if the orbits were isolated. However, taking into account dynamical evolution, we find that $\approx 0.4$ ($\approx 0.6$) of 2+2 (3+1) systems interact. In particular, Roche Lobe overflow can be triggered possibly in highly eccentric orbits, dynamical instability can ensue due to mass-loss-driven orbital expansion or secular evolution, or a semisecular regime can be entered. We compute the delay-time distributions (DTDs) of collision-induced SNe Ia, and find that they are flatter compared to the observed DTD. Moreover, our combined SNe Ia rates are $(3.7\pm0.7) \times 10^{-6} \,\mathrm{M}_\odot^{-1}$ and $(1.3\pm0.2) \times 10^{-6} \,\mathrm{M}_\odot^{-1}$ for 2+2 and 3+1 systems, respectively, three orders of magnitude lower compared to the observed rate, of order $10^{-3} \, \mathrm{M}_\odot^{-1}$. The low rates can be ascribed to interactions before the stars evolve to CO WDs. However, our results are lower limits given that we considered a subset of quadruple systems.


title: The Size-Luminosity Relationship of Quasar Narrow-Line Regions

abstract: The presence of an active galactic nucleus (AGN) can strongly affect its host. Due to the copious radiative power of the nucleus, the effects of radiative feedback can be detected over the entire host galaxy and sometimes well into the intergalactic space. In this paper we model the observed size-luminosity relationship of the narrow-line regions (NLRs) of AGN. We model the NLR as a collection of clouds in pressure equilibrium with the ionizing radiation, with each cloud producing line emission calculated by Cloudy. The sizes of the NLRs of powerful quasars are reproduced without any free parameters, as long as they contain massive ($10^5 - 10^7 M_\odot$) ionization-bounded clouds. At lower AGN luminosities the observed sizes are larger than the model sizes, likely due to additional unmodeled sources of ionization (e.g., star formation). We find that the observed saturation of sizes at $\sim 10$ kpc which is observed at high AGN luminosities ($L_\text{ion} \simeq 10^{46}$ erg/s) is naturally explained by optically thick clouds absorbing the ionizing radiation and preventing illumination beyond a critical distance. Using our models in combination with observations of the [O III]/IR ratio and the [O III] size -- IR luminosity relationship, we calculate the covering factor of the obscuring torus (and therefore the type 2 fraction within the quasar population) to be $f=0.5$, though this is likely an upper bound. Finally, because the gas behind the ionization front is invisible in ionized gas transitions, emission-based NLR mass calculations underestimate the mass of the NLR and therefore of the energetics of ionized-gas winds.


title: On the kinematics of a runaway Be star population

abstract: We explore the hypothesis that B type emission-line stars (Be stars) have their origin in mass-transfer binaries by measuring the fraction of runaway Be stars. We assemble the largest-to-date catalogue of 632 Be stars with 6D kinematics, exploiting the precise astrometry of the Tycho-Gaia Astrometric Solution (TGAS) from the first Gaia Data Release. Using binary stellar evolution simulations, we make predictions for the runaway and equatorial rotation velocities of a runaway Be star population. Accounting for observational biases, we calculate that if all classical Be stars originated through mass transfer in binaries, then $17.5\%$ of the Be stars in our catalogue should be runaways. The remaining $82.5\%$ should be in binaries with subdwarfs, white dwarfs or neutron stars, because those systems either remained bound post-supernova or avoided the supernova entirely. Using a Bayesian methodology, we compare the hypothesis that each Be star in our catalogue is a runaway to the null hypothesis that it is a member of the Milky Way disc. We find that $13.1^{+2.6}_{-2.4}\%$ of the Be stars in our catalogue are runaways, and identify a subset of 40 high-probability runaways. We argue that deficiencies in our understanding of binary stellar evolution, as well as the degeneracy between velocity dispersion and number of runaway stars, can explain the slightly lower runaway fraction. We thus conclude that all Be stars could be explained by an origin in mass-transfer binaries. This conclusion is testable with the second Gaia data release (DR2).


title: Dust Attenuation Curves in the Local Universe: Demographics and New Laws for Star-forming Galaxies and High-redshift Analogs

abstract: We study dust attenuation curves of 230,000 individual galaxies in the local universe, ranging from quiescent to intensely star-forming systems, using GALEX, SDSS, and WISE photometry calibrated on Herschel-ATLAS. We use a new method of constraining SED fits with infrared luminosity (SED+LIR fitting), and parameterized attenuation curves determined with the CIGALE SED fitting code. Attenuation curve slopes and UV bump strengths are reasonably well constrained independently from one another. We find that $A_{\lambda}/A_V$ attenuation curves exhibit a very wide range of slopes that are on average as steep as the SMC curve slope. The slope is a strong function of optical opacity. Opaque galaxies have shallower curves - in agreement with recent radiate transfer models. The dependence of slopes on the opacity produces an apparent dependence on stellar mass: more massive galaxies having shallower slopes. Attenuation curves exhibit a wide range of UV bump amplitudes, from none to MW-like; with an average strength 1/3 of the MW bump. Notably, local analogs of high-redshift galaxies have an average curve that is somewhat steeper than the SMC curve, with a modest UV bump that can be to first order ignored, as its effect on the near-UV magnitude is 0.1 mag. Neither the slopes nor the strengths of the UV bump depend on gas-phase metallicity. Functional forms for attenuation laws are presented for normal star-forming galaxies, high-z analogs and quiescent galaxies. We release the catalog of associated SFRs and stellar masses (GSWLC-2).


title: Holistic spectroscopy: Complete reconstruction of a wide-field, multi-object spectroscopic image using a photonic comb

abstract: The primary goal of Galactic archaeology is to learn about the origin of the Milky Way from the detailed chemistry and kinematics of millions of stars. Wide-field multi-fibre spectrographs are increasingly used to obtain spectral information for huge samples of stars. Some surveys (e.g. GALAH) are attempting to measure up to 30 separate elements per star. Stellar abundance spectroscopy is a subtle art that requires a very high degree of spectral uniformity across each of the fibres. However wide-field spectrographs are notoriously non-uniform due to the fast output optics necessary to image many fibre outputs onto the detector. We show that precise spectroscopy is possible with such instruments across all fibres by employing a photonic comb -- a device that produces uniformly spaced spots of light on the CCD to precisely map complex aberrations. Aberrations are parametrized by a set of orthogonal moments with $\sim100$ independent parameters. We then reproduce the observed image by convolving high resolution spectral templates with measured aberrations as opposed to extracting the spectra from the observed image. Such a forward modeling approach also trivializes some spectroscopic reduction problems like fibre cross-talk, and reliably extracts spectra with a resolution $\sim2.3$ times above the nominal resolution of the instrument. Our rigorous treatment of optical aberrations also encourages a less conservative spectrograph design in the future.


title: Analytic inclusion of the scale dependence of the anomalous dimension matrix in Standard Model Effective Theory

abstract: The renormalization group equations (RGEs) in Standard Model effective theory are usually either solved analytically, neglecting the scale dependence of gauge and Yukawa couplings, or numerically without such approximations. We present analytic solutions of RGEs that take into account the dominant scale dependence of the anomalous-dimension matrix due to the running of the QCD coupling and of the top-Yukawa coupling. We consider first the case for which a given operator is generated directly through mixing with the parent operator whose Wilson coefficient is non-vanishing at the new physics scale. Subsequently we consider the case of two-step running, in which two operators do not mix directly, but only via a third mediator operator. We generalize these solutions to an arbitrary number of operators and show how even in this case analytic solutions can be obtained.


title: The EDGE-CALIFA Survey: Molecular and Ionized Gas Kinematics in Nearby Galaxies

abstract: We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved $^{12}$CO(J=1-0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and H$\alpha$ rotation curves using the same geometric parameters out to $\gtrsim$1 $R_e$. We find that $\sim$75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence supported disk. Using observations of the H$\gamma$ transition also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and H$\alpha$ rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.


title: Quantum optics of spin waves through ac Stark modulation

abstract: We bring the set of linear quantum operations, important for many fundamental studies in photonic systems, to the material domain of collective excitations known as spin waves. Using the ac Stark effect we realize quantum operations on single excitations and demonstrate a spin-wave analogue of Hong-Ou-Mandel effect, realized via a beamsplitter implemented in the spin wave domain. Our scheme equips atomic-ensemble-based quantum repeaters with quantum information processing capability and can be readily brought to other physical systems, such as doped crystals or room-temperature atomic ensembles.


title: Entanglement Wedge Reconstruction and Entanglement of Purification

abstract: In the holographic correspondence, subregion duality posits that knowledge of the mixed state of a finite spacelike region of the boundary theory allows full reconstruction of a specific region of the bulk, known as the entanglement wedge. This statement has been proven for local bulk operators. In this paper, specializing first for simplicity to a Rindler wedge of AdS$_3$, we find that generic curves within the wedge are in fact not fully reconstructible with entanglement entropies in the corresponding boundary region, even after using the most general variant of hole-ography, which was recently shown to suffice for reconstruction of arbitrary spacelike curves in the Poincare patch. This limitation is an analog of the familiar phenomenon of entanglement shadows, which we call 'entanglement shade'. We overcome it by showing that the information about the nonreconstructible curve segments is encoded in a slight generalization of the concept of entanglement of purification, whose holographic dual has been discussed very recently. We introduce the notion of 'differential purification', and demonstrate that, in combination with differential entropy, it enables the complete reconstruction of all spacelike curves within an arbitrary entanglement wedge in any 3-dimensional bulk geometry.


title: Strategies for optimal single-shot discrimination of quantum measurements

abstract: In this work we study the problem of single-shot discrimination of von Neumann measurements, which we associate with measure-and-prepare channels. There are two possible approaches to this problem. The first one is simple and does not utilize entanglement. We focus only on the discrimination of classical probability distributions, which are outputs of the channels. We find necessary and sufficient criterion for perfect discrimination in this case. A more advanced approach requires the usage of entanglement. We quantify the distance between two measurements in terms of the diamond norm (called sometimes the completely bounded trace norm). We provide an exact expression for the optimal probability of correct distinction and relate it to the discrimination of unitary channels. We also state a necessary and sufficient condition for perfect discrimination and a semidefinite program which checks this condition. Our main result, however, is a cone program which calculates the distance between the measurements and hence provides an upper bound on the probability of their correct distinction. As a by-product, the program finds a strategy (input state) which achieves this bound. Finally, we provide a full description for the cases of Fourier matrices and mirror isometries.


title: Gravitational Bose-Einstein condensation in the kinetic regime

abstract: We study Bose-Einstein condensation and formation of Bose stars in the virialized dark matter halos/miniclusters by universal gravitational interactions. We prove that this phenomenon does occur and it is described by kinetic equation. We give expression for the condensation time. Our results suggest that Bose stars may form kinetically in the mainstream dark matter models such as invisible QCD axions and Fuzzy Dark Matter.

categories: astro-ph.CO astro-ph.GA cond-mat.stat-mech hep-ph


title: Fractionalized Metal in a Falicov-Kimball Model

abstract: Quantum Monte Carlo simulations reveal an exotic metallic phase with a single-particle gap but gapless spin and charge excitations and a nonsaturating resistivity in a two-dimensional SU(2) Falicov-Kimball model. An exact duality between this model and an unconstrained slave-spin theory leads to a classification of the phase as a fractionalized or orthogonal metal whose low-energy excitations have different quantum numbers than the original electrons. Whereas the fractionalized metal corresponds to the regime of disordered slave spins, the regime of ordered slave spins is a Fermi liquid. At a critical temperature, an Ising phase transition to a spontaneously generated constrained slave-spin theory of the Hubbard model is observed.


title: The average number of rational points on genus two curves is bounded

abstract: We prove that, when genus two curves $C/\mathbb{Q}$ with a marked Weierstass point are ordered by height, the average number of rational points $\#|C(\mathbb{Q})|$ is bounded. The argument follows the same ideas as the sphere-packing proof of boundedness of the average number of integral points on (quasiminimal Weierstrass models of) elliptic curves. That is, we bound the number of small-height points by hand, the number of medium-height points by establishing an explicit Mumford gap principle and using the theorem of Kabatiansky-Levenshtein on spherical codes (this technique goes back to work of Silverman, Helfgott, and Helfgott-Venkatesh), and the number of large-height points by using Bombieri-Vojta's proof of Faltings' theorem. Explicitly, in dealing with non-small-height points we prove that the number of rational points $(x,y)$ on $C_f: y^2 = f(x)$ satisfying $h(x) > 8 h(f)$ is $\ll 1.872^{\mathrm{rank}(\mathrm{Jac}(C)(\mathbb{Q}))}$, which has finite average by the theorem of Bhargava-Gross on the average size of $2$-Selmer groups of Jacobians over this family. We note that our arguments in the small-height and large-height cases extend to general genera $g\geq 2$, though for medium points we need to use Stoll's bounds on the non-Archimedean local height differences in genus $2$. For example, we prove that the number of rational points $P\in C(\mathbb{Q})$ with $h(P)\gg_g h(C)$ on $C/\mathbb{Q}$ smooth projective and of genus $g\geq 2$ is $\ll 1.872^{\mathrm{rank}(\mathrm{Jac}(C)(\mathbb{Q}))}$, and that in fact the base of the exponent can be reduced to $1.311$ once $g\gg 1$, though this is surely known to experts (the difference is the use of the Kabatiansky-Levenshtein bound in lieu of more elementary techniques).


title: A Complete Census of Luminous Stellar Variability on Day to Decade Timescales

abstract: Stellar photometric variability offers a novel probe of the interior structure and evolutionary state of stars. Here we present a census of stellar variability on day to decade timescales across the color-magnitude diagram for 73,000 stars brighter than $M_I$=-5 in the Whirlpool Galaxy (M51). Our Cycle 24 HST program acquired V and I-band images over 34 epochs spanning one year with pseudo-random cadences enabling sensitivity to periods from days to months. We supplement these data with archival V and I-band HST data obtained in 1995 and 2005 providing sensitivity to variability on decade timescales. At least 50% of stars brighter than $M_I$=-7 show strong evidence for variability within our Cycle 24 data; among stars with V-I>2 the variability fraction rises to ~100%. Large amplitude variability (>0.3 mag) on decade timescales is restricted to red supergiants and very luminous blue stars. Both populations display fairly smooth variability on month-year timescales. The Cepheid instability strip is clearly visible in our data, although the variability fraction within this region never exceeds ~10%. The location of variable stars across the color magnitude diagram broadly agrees with theoretical sources of variability, including the instability strip, red supergiant pulsational instabilities, long-period fundamental mode pulsations, and radiation-dominated envelopes in massive stars. Our data can be used to place stringent constraints on the precise onset of these various instabilities and their lifetimes and growth rates.


title: Differentiating Disk and Black Hole Driven Jets with EHT Images of Variability in M87

abstract: Millimetre-wavelength very long baseline interferometric (mm-VLBI) observations of M87 by the Event Horizon Telescope (EHT) should provide a unique opportunity to observe and characterize the origins of jet variability already seen at longer wavelengths. Synchrotron spot models have been used to model variability near black holes; this work extends these by allowing spots to shear and deform in the jet velocity field. Depending on the position of the spot, shearing forces can significantly alter the structure of the spot, producing distinct signals in reconstructed images and light curves. The maximum intensity of the shearing spot can vary by as much as a factor of five depending on the spot azimuthal launch position, but the intensity decay time depends most significantly on the spot radial launch position. Spots launched by a black hole driven jet exhibit distinct arc structures in reconstructed images, and exhibit brighter and shorter-lived enhancements of the light curve. Spots launched by a wind-driven jet have exhibit much simpler structures in the image, and longer-lived light curve enhancements than spots launched by a black hole driven jet.


title: Non-Vacuous Generalization Bounds at the ImageNet Scale: A PAC-Bayesian Compression Approach

abstract: Modern neural networks are highly overparameterized, with capacity to substantially overfit to training data. Nevertheless, these networks often generalize well in practice. It has also been observed that trained networks can often be "compressed" to much smaller representations. The purpose of this paper is to connect these two empirical observations. Our main technical result is a generalization bound for compressed networks based on the compressed size. Combined with off-the-shelf compression algorithms, the bound leads to state of the art generalization guarantees; in particular, we provide the first non-vacuous generalization guarantees for realistic architectures applied to the ImageNet classification problem. As additional evidence connecting compression and generalization, we show that compressibility of models that tend to overfit is limited: We establish an absolute limit on expected compressibility as a function of expected generalization error, where the expectations are over the random choice of training examples. The bounds are complemented by empirical results that show an increase in overfitting implies an increase in the number of bits required to describe a trained network.


title: Introduction to Effective Field Theories

abstract: Lecture notes from the 2017 Les Houches Summer School on Effective Field Theories. The lectures covered introductory material on EFTs as used in high energy physics to compute experimentally observable quantities. Other lectures at the school covered a wide range of applications in greater depth.


title: Setting the Stage for Cosmic Chronometers. I. Assessing the Impact of Young Stellar Populations on Hubble Parameter Measurements

abstract: The expansion history of the Universe can be constrained in a cosmology-independent way by measuring the differential age evolution of cosmic chronometers. This yields a measurement of the Hubble parameter $H(z)$ as a function of redshift. The most reliable cosmic chronometers known so far are extremely massive and passively evolving galaxies. Age-dating these galaxies is, however, a difficult task, and even a small contribution of an underlying young stellar population could, in principle, affect the age estimate and its cosmological interpretation. We present several spectral indicators to detect, quantify and constrain such contamination in old galaxies, and study how their combination can be used to maximize the purity of cosmic chronometers selection. In particular, we analyze the CaII H/K ratio, the presence (or absence) of H$\alpha$ and [OII] emission lines, higher order Balmer absorption lines, and UV flux; each indicator is especially sensitive to a particular age range, allowing us to detect young components ranging between 10 Myr and 1 Gyr. The combination of these indicators minimizes the contamination to a level below 1% in the case of ideal data. More importantly, it offers a way to control the systematic error on $H(z)$ as a function of the contamination by young stellar populations. We show that for our previous measurements of the Hubble parameter, the possible bias induced by the presence of a younger component is well below the current errors. We envision that these indicators will be instrumental in paving the road for a robust and reliable dating of the old population and its cosmological interpretation.


title: The Aemulus Project I: Numerical Simulations for Precision Cosmology

abstract: The rapidly growing statistical precision of galaxy surveys has lead to a need for ever-more precise predictions of the observables used to constrain cosmological and galaxy formation models. The primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper we present one such suite of simulations, forming the basis for the AEMULUS Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h^-1 Gpc)^3 simulations with mass resolution and force softening of 3.51\times 10^10 (Omega_m / 0.3) ~ h^-1 M_sun and 20 ~ h^-1 kpc respectively in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) for halos with more than 500 (200) particles respectively and scales of r>200 ~ h^-1 kpc in real space or k ~ 3 h Mpc^-1 in harmonic space for z\le 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k ~ 4 h^-1 Mpc. We provide the halo catalogs and snapshots detailed in this work to the community at https://AemulusProject.github.io.


title: The Aemulus Project II: Emulating the Halo Mass Function

abstract: Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a halo mass function emulator and demonstrate improved accuracy relative to state-of-the-art analytic models. In this work, mass is defined using an overdensity criteria of 200 relative to the mean background density. Our emulator is constructed from the AEMULUS simulations, a suite of 40 N-body simulations with snapshots from z=3 to z=0. These simulations cover the flat wCDM parameter space allowed by recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results, varying the parameters w, Omega_m, Omega_b, sigma_8, N_{eff}, n_s, and H_0. We validate our emulator using five realizations of seven different cosmologies, for a total of 35 test simulations. These test simulations were not used in constructing the emulator, and were run with fully independent initial conditions. We use our test simulations to characterize the modeling uncertainty of the emulator, and introduce a novel way of marginalizing over the associated systematic uncertainty. We confirm non-universality in our halo mass function emulator as a function of both cosmological parameters and redshift. Our emulator achieves better than 1% precision over much of the relevant parameter space, and we demonstrate that the systematic uncertainty in our emulator will remain a negligible source of error for cluster abundance studies through at least the LSST Year 1 data set.


title: The Aemulus Project III: Emulation of the Galaxy Correlation Function

abstract: Using the N-body simulations of the AEMULUS Project, we construct an emulator for the non-linear clustering of galaxies in real and redshift space. We construct our model of galaxy bias using the halo occupation framework, accounting for possible velocity bias. The model includes 15 parameters, including both cosmological and galaxy bias parameters. We demonstrate that our emulator achieves ~ 1% precision at the scales of interest, 0.1<r<10 h^{-1} Mpc, and recovers the true cosmology when tested against independent simulations. Our primary parameters of interest are related to the growth rate of structure, f, and its degenerate combination fsigma_8. Using this emulator, we show that the constraining power on these parameters monotonically increases as smaller scales are included in the analysis, all the way down to 0.1 h^{-1} Mpc. For a BOSS-like survey, the constraints on fsigma_8 from r<30 h^{-1} Mpc scales alone are more than a factor of two tighter than those from the fiducial BOSS analysis of redshift-space clustering using perturbation theory at larger scales. The combination of real- and redshift-space clustering allows us to break the degeneracy between f and sigma_8, yielding a 9% constraint on f alone for a BOSS-like analysis. The current AEMULUS simulations limit this model to surveys of massive galaxies. Future simulations will allow this framework to be extended to all galaxy target types, including emission-line galaxies.


title: Universal Dependency Parsing for Hindi-English Code-switching

abstract: Code-switching is a phenomenon of mixing grammatical structures of two or more languages under varied social constraints. The code-switching data differ so radically from the benchmark corpora used in NLP community that the application of standard technologies to these data degrades their performance sharply. Unlike standard corpora, these data often need to go through additional processes such as language identification, normalization and/or back-transliteration for their efficient processing. In this paper, we investigate these indispensable processes and other problems associated with syntactic parsing of code-switching data and propose methods to mitigate their effects. In particular, we study dependency parsing of code-switching data of Hindi and English multilingual speakers from Twitter. We present a treebank of Hindi-English code-switching tweets under Universal Dependencies scheme and propose a neural stacking model for parsing that efficiently leverages part-of-speech tag and syntactic tree annotations in the code-switching treebank and the preexisting Hindi and English treebanks. We also present normalization and back-transliteration models with a decoding process tailored for code-switching data. Results show that our neural stacking parser is 1.5% LAS points better than the augmented parsing model and our decoding process improves results by 3.8% LAS points over the first-best normalization and/or back-transliteration.


title: The GALAH survey: An abundance, age, and kinematic inventory of the solar neighbourhood made with TGAS

abstract: The overlap between the spectroscopic Galactic Archaeology with HERMES (GALAH) survey & $Gaia$ provides a high-dimensional chemodynamical space of unprecedented size. We present a first analysis of a subset of this overlap, of 7066 dwarf, turn-off, & sub-giant stars. [...] We investigate correlations between chemical compositions, ages, & kinematics for this sample. Stellar parameters & elemental abundances are derived from the GALAH spectra with the spectral synthesis code SME. [...] We report Li, C, O, Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, as well as Ba & we note that we employ non-LTE calculations for Li, O, Al, & Fe. We show that the use of astrometric & photometric data improves the accuracy of the derived spectroscopic parameters, especially $\log g$. [...] we recover the result that stars of the high-$\alpha$ sequence are typically older than stars in the low-$\alpha$ sequence, the latter spanning $-0.7<$[Fe/H]$<+0.5$. While these two sequences become indistinguishable in [$\alpha$/Fe] vs. [Fe/H] at the metal-rich regime, we find that age can be used to separate stars from the extended high-$\alpha$ & the low-$\alpha$ sequence even in this regime. [...] we find that the old stars ($>8$ Gyr have lower angular momenta $L_z$ than the Sun, which implies that they are on eccentric orbits & originate from the inner disk. Contrary to some previous smaller scale studies we find a continuous evolution in the high-$\alpha$-sequence up to super-solar [Fe/H] rather than a gap, which has been interpreted as a separate "high-$\alpha$ metal-rich" population. Stars in our sample that are younger than 10 Gyr, are mainly found on the low $\alpha$-sequence & show a gradient in $L_z$ from low [Fe/H] ($L_z>L_{z,\odot}$) towards higher [Fe/H] ($L_z<L_{z,\odot}$), which implies that the stars at the ends of this sequence are likely not originating from the close solar vicinity.


title: Egocentric 6-DoF Tracking of Small Handheld Objects

abstract: Virtual and augmented reality technologies have seen significant growth in the past few years. A key component of such systems is the ability to track the pose of head mounted displays and controllers in 3D space. We tackle the problem of efficient 6-DoF tracking of a handheld controller from egocentric camera perspectives. We collected the HMD Controller dataset which consist of over 540,000 stereo image pairs labelled with the full 6-DoF pose of the handheld controller. Our proposed SSD-AF-Stereo3D model achieves a mean average error of 33.5 millimeters in 3D keypoint prediction and is used in conjunction with an IMU sensor on the controller to enable 6-DoF tracking. We also present results on approaches for model based full 6-DoF tracking. All our models operate under the strict constraints of real time mobile CPU inference.


title: Limits of multiplicative inhomogeneous random graphs and L\'evy trees: The continuum graphs

abstract: Motivated by limits of critical inhomogeneous random graphs, we construct a family of sequences of measured metric spaces that we call continuous multiplicative graphs, that are expected to be the universal limit of graphs related to the multiplicative coalescent (the Erd\H{o}s--R\'enyi random graph, more generally the so-called rank-one inhomogeneous random graphs of various types, and the configuration model). At the discrete level, the construction relies on a new point of view on (discrete) inhomogeneous random graphs that involves an embedding into a Galton--Watson forest. The new representation allows us to demonstrate that a processus that was already present in the pionnering work of Aldous [Ann. Probab., vol.~25, pp.~812--854, 1997] and Aldous and Limic [Electron. J. Probab., vol.~3, pp.~1--59, 1998] about the multiplicative coalescent actually also (essentially) encodes the limiting metric: The discrete embedding of random graphs into a Galton--Watson forest is paralleled by an embedding of the encoding process into a L\'evy process which is crucial in proving the very existence of the local time functionals on which the metric is based; it also yields a transparent approach to compactness and fractal dimensions of the continuous objects. In a companion paper, we show that the continuous L\'evy graphs are indeed the scaling limit of inhomogeneous random graphs.


title: Chiral supercurrent through a quantum Hall weak link

abstract: We use a microscopic model to calculate properties of the supercurrent carried by chiral edge states of a quantum Hall weak link. This "chiral" supercurrent is qualitatively distinct from the usual Josephson supercurrent in that it cannot be mediated by a single edge alone, i.e., both right and left going edges are needed. Moreover, chiral supercurrent was previously shown to obey an unusual current-phase relation with period $2 \phi_0=h/e$, which is twice as large as the period of conventional Josephson junctions. We show that the "chiral" nature of this supercurrent is sharply defined, and is robust to interactions to infinite order in perturbation theory. We compare our results with recent experimental findings of Amet et al [Science, 352(6288)] and find that quantitative agreement in magnitude of the supercurrent can be attained by making reasonable but critical assumptions about the superconductor quantum Hall interface.


title: Weak-Lensing Mass Calibration of ACTPol Sunyaev-Zel'dovich Clusters with the Hyper Suprime-Cam Survey

abstract: We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.


title: Probing the Top-Higgs Yukawa CP Structure in dileptonic $t \bar t h$ with $M_2$-Assisted Reconstruction

abstract: Constraining the Higgs boson properties is a cornerstone of the LHC program. We study the potential to directly probe the Higgs-top CP-structure via the $t\bar{t}h$ channel at the LHC with the Higgs boson decaying to a bottom pair and top-quarks in the dileptonic mode. We show that a combination of laboratory and $t\bar{t}$ rest frame observables display large CP-sensitivity, exploring the spin correlations in the top decays. To efficiently reconstruct our final state, we present a method based on simple mass minimization and prove its robustness to shower, hadronization and detector effects. In addition, the mass reconstruction works as an extra relevant handle for background suppression. Based on our results, we demonstrate that the Higgs-top CP-phase $(\alpha)$ can be probed up to $\cos\alpha< 0.7$ at the high luminosity LHC.


title: The Dirichlet problem for semi-linear equations

abstract: We study the Dirichlet problem for the semi--linear partial differential equations ${\rm div}\,(A\nabla u)=f(u)$ in simply connected domains $D$ of the complex plane $\mathbb C$ with continuous boundary data. We prove the existence of the weak solutions $u$ in the class $C\cap W^{1,2}_{\rm loc}(D)$ if a Jordan domain $D$ satisfies the quasihyperbolic boundary condition by Gehring--Martio. An example of such a domain that fails to satisfy the standard (A)--condition by Ladyzhenskaya--Ural'tseva and the known outer cone condition is given. We also extend our results to simply connected non-Jordan domains formulated in terms of the prime ends by Caratheodory. Our approach is based on the theory of the logarithmic potential, singular integrals, the Leray--Schauder technique and a factorization theorem in \cite{GNR2017}. This theorem allows us to represent $u$ in the form $u=U\circ\omega,$ where $\omega(z)$ stands for a quasiconformal mapping of $D$ onto the unit disk ${\mathbb D}$, generated by the measurable matrix function $A(z),$ and $U$ is a solution of the corresponding quasilinear Poisson equation in the unit disk ${\mathbb D}$. In the end, we give some applications of these results to various processes of diffusion and absorption in anisotropic and inhomogeneous media.


title: On the End-Fire (Super) Directivity of an Array of Two Elementary Dipoles

abstract: The concept of source currents of a radiating source can be employed to express directivity in some particular cases analytically. For an antenna array, this concept can be combined with concepts of mutual radiation intensity and the mutual power of array elements. Using these approaches, we treat the evaluation of the directivity for an elementary dipole, and its end-fire array, which is of special attention due to superdirective properties. A closed form expression for the directivity of this array with out-of-phase excitation is derived. It is observed that end-fire directivity can be further enhanced by optimizing the excitation currents of the array. Their optimal relative phase and corresponding increased directivity are also found analytically. The results are validated by a full-wave simulator.


title: Constraining cosmic ray acceleration in young star clusters using multi-wavelength observations

abstract: We use $1$D and $3$D two-fluid cosmic ray (CR) hydrodynamic simulations to investigate the role of CRs in the vicinity of a compact young star cluster. We model a self-gravitating cloud (density profile $\rho \propto r^{-1}$), include important thermal and non-thermal processes, and explore two different CR injection scenarios. We show that if internal shocks in the wind-driving region are the main site for CR acceleration, then the resulting $\gamma$-ray luminosity ($L_{\rm \gamma}$) can reach $\approx 5\%$ of the mechanical luminosity ($L_{\rm w}$), independent of the fraction of wind energy ($\sim 1-20\%$) injected into CRs. In contrast, if the forward/reverse shock of a bubble is the injection site then $L_{\rm \gamma}$ increases linearly with the CR injection fraction, as expected analytically. We find that the X-ray luminosity ($L_{\rm x}$) in the forward/reverse shock injection scenario is $\gtrsim 10^{-3} L_{\rm w}$, which is $\sim 10$ times larger than in the central wind-driving injection case. We predict the corresponding range of the synchrotron radio luminosity. We show how multi-wavelength observations can constrain the CR parameters. Comparing the predicted multi-wavelength luminosities with those of 30 Doradus we identify the reverse shock as the most probable CR injection site, and that thermal conduction is important. We do not find significant dynamical impact of CRs in our models.


title: Tips for Deciphering and Quick Calculation of Radiation Spectra

abstract: Radiation spectra from ultra-relativistic electrons in thin [$T\ll l_f(\omega)$] and thick [$T\gg l_f(\omega)$] targets are discussed. The method of simplified averaging is described by examples of Landau-Pomeranchuk-Migdal effect and radiation at doughnut scattering. General infrared and ultraviolet asymptotic properties of radiation spectra are discussed.


title: M-PACT: An Open Source Platform for Repeatable Activity Classification Research

abstract: There are many hurdles that prevent the replication of existing work which hinders the development of new activity classification models. These hurdles include switching between multiple deep learning libraries and the development of boilerplate experimental pipelines. We present M-PACT to overcome existing issues by removing the need to develop boilerplate code which allows users to quickly prototype action classification models while leveraging existing state-of-the-art (SOTA) models available in the platform. M-PACT is the first to offer four SOTA activity classification models, I3D, C3D, ResNet50+LSTM, and TSN, under a single platform with reproducible competitive results. This platform allows for the generation of models and results over activity recognition datasets through the use of modular code, various preprocessing and neural network layers, and seamless data flow. In this paper, we present the system architecture, detail the functions of various modules, and describe the basic tools to develop a new model in M-PACT.


title: Uniform Substitution for Differential Game Logic

abstract: This paper presents a uniform substitution calculus for differential game logic (dGL). Church's uniform substitutions substitute a term or formula for a function or predicate symbol everywhere. After generalizing them to differential game logic and allowing for the substitution of hybrid games for game symbols, uniform substitutions make it possible to only use axioms instead of axiom schemata, thereby substantially simplifying implementations. Instead of subtle schema variables and soundness-critical side conditions on the occurrence patterns of logical variables to restrict infinitely many axiom schema instances to sound ones, the resulting axiomatization adopts only a finite number of ordinary dGL formulas as axioms, which uniform substitutions instantiate soundly. This paper proves soundness and completeness of uniform substitutions for the monotone modal logic dGL. The resulting axiomatization admits a straightforward modular implementation of dGL in theorem provers.


title: Evaluating Massive MIMO Precoding based on 3D-Channel Measurements with a Spider Antenna

abstract: Massive Multiple-Input Multiple-Output (MIMO)communications uses a large number of antennas at the base station to increase the data rate and user density in future wireless systems. For simulation, it has become common practice to use i.i.d. complex Gaussian matrix entries to obtain an average MIMO channel behavior. More refined models have been devised and proposed to standardization bodies; yet, channel modeling remains an active area of research, as current models tend to be, still, quite limited, e.g., when it comes to evaluating clustering algorithms, with regions of spatial orthogonality for concurrent scheduling of users, which is an essential concept in massive MIMO precoding. For this, spatial correlations need to be included. To further refine channel modeling, we have built a "spider antenna" prototype that allows spatially continuous measurements in three dimensions, enabling a high-resolution sampling over, initially, a volume of 2m x 2m x 2m for indoor measurements. Several experiments have been conducted to illustrate the new insights to be gained when studying user orthogonality, clustering and precoding in a massive MIMO context. Furthermore, the influence of antenna array geometry and user spacing on the achievable rate over actually measured channels is studied.


title: Confidence intervals for the area under the receiver operating characteristic curve in the presence of ignorable missing data

abstract: Receiver operating characteristic (ROC) curves are widely used as a measure of accuracy of diagnostic tests and can be summarized using the area under the ROC curve (AUC). Often, it is useful to construct a confidence intervals for the AUC, however, since there are a number of different proposed methods to measure variance of the AUC, there are thus many different resulting methods for constructing these intervals. In this manuscript, we compare different methods of constructing Wald-type confidence interval in the presence of missing data where the missingness mechanism is ignorable. We find that constructing confidence intervals using multiple imputation (MI) based on logistic regression (LR) gives the most robust coverage probability and the choice of CI method is less important. However, when missingness rate is less severe (e.g. less than 70%), we recommend using Newcombe's Wald method for constructing confidence intervals along with multiple imputation using predictive mean matching (PMM).


title: Periodic orbits around Kerr Sen black holes

abstract: We investigate periodic orbits and zoom-whirl behaviors around a Kerr Sen black hole with a rational number $q$ in terms of three integers $(z,w,v)$, from which one can immediately read off the number of leaves(or zooms), the ordering of the leaves, and the number of whirls. The characteristic of zoom-whirl periodic orbits is the precession of multi-leaf orbits in the strong field regime. This feature is analogous to the counterpart in the Kerr space-time. Finally, we analyze the impact of the charge parameter $b$ on the zoom-whirl periodic orbits. Compared to the periodic orbits around the Kerr black hole, it is found that typically lower energies are required for the same orbits in the Kerr Sen black hole.


title: Automorphic Forms on Feit's Hermitian Lattices

abstract: We consider the genus of $20$ classes of unimodular Hermitian lattices of rank $12$ over the Eisenstein integers. This set is the domain for a certain space of algebraic modular forms. We find a basis of Hecke eigenforms, and guess global Arthur parameters for the associated automorphic representations, which recover the computed Hecke eigenvalues. Congruences between Hecke eigenspaces, combined with the assumed parameters, recover known congruences for classical modular forms, and support new instances of conjectured Eisenstein congruences for $\mathbb{U}(2,2)$ automorphic forms.


title: Area Rate Evaluation based on Spatial Clustering of massive MIMO Channel Measurements

abstract: Channel models for massive MIMO are typically based on matrices with complex Gaussian entries, extended by the Kronecker and Weichselberger model. One reason for observing a gap between modeled and actual channel behavior is the absence of spatial consistency in many such models, that is, spatial correlations over an area in the x, y-dimensions are not accounted for, making it difficult to study, e.g., area-throughput measures. In this paper, we propose an algorithm that can distinguish between regions of non-line-of-sight (NLoS) and line-of-sight (LoS) via a rank-metric criterion combined with a spiral search. With a k-means clustering algorithm a throughput per region (i.e., cluster) can be calculated, leading to what we refer to as "area-throughput". For evaluating the proposed orthogonality clustering scheme we use a simple filtered MIMO channel model which is spatially consistent, with known degrees of freedom. Moreover, we employ actual (spatially consistent) area channel measurements based on spatial sampling using a spider antenna and show that the proposed algorithm can be used to estimate the degrees of freedom, and, subsequently, the number of users that maximizes the throughput per square meter.


title: Subcarrier-Interlaced FDD for Faster-than-TDD Channel Tracking in Massive MIMO Systems

abstract: Canonical Massive MIMO uses time division duplex (TDD) to exploit channel reciprocity within the coherence time, avoiding feedback of channel state information (CSI), as is required for precoding at the base station. We extend the idea of exploiting reciprocity to the coherence bandwidth, allocating subcarriers of a multicarrier (OFDM)-based system in an interlaced fashion to up- and downlink, respectively, referred to as (OFDM-)subcarrier interlaced FDD (IFDD). Exploiting this "two-dimensional" channel reciprocity within both time and frequency coherence does not require any CSI feedback and, even more so, offers faster-than-TDD channel tracking. To address imperfections of actual hardware, we conducted measurements, verifying the practical viability of IFDD. As it turns out, the scheme incurs similar transmit/receive isolation issues (self-interference) as well-known from the full-duplex (FDup) literature. As will be shown, such hardware challenges like self-interference or frequency offsets are not critical for massive MIMO operation, but can be neglected as the number of antennas grows large. Finally, we quantify how IFDD allows to track the channel variations much faster than TDD over a wide range of commonly used pilot symbol rates.


title: Multiple solutions and their asymptotics for laminar flows through a porous channel with different permeabilities

abstract: The existence and multiplicity of similarity solutions for the steady, incompressible and fully developed laminar flows in a uniformly porous channel with two permeable walls are investigated. We shall focus on the so-called asymmetric case where the upper wall is with an amount of flow injection and the lower wall with a different amount of suction. We show that there exist three solutions designated as type $I$, type $II$ and type $III$ for the asymmetric case. The numerical results suggest that a unique solution exists for the Reynolds number $0\leq R<14.10$ and two additional solutions appear for $R>14.10$. The corresponding asymptotic solution for each of the multiple solutions is constructed by the method of boundary layer correction or matched asymptotic expansion for the most difficult high Reynolds number case. Asymptotic solutions are all verified by their corresponding numerical solutions.


title: Oversampling and undersampling in de Branges spaces arising from regular Schr\"odinger operators

abstract: The classical results on oversampling and undersampling (or aliasing), of functions in Paley-Wiener spaces, are generalized to the case of functions in de Branges spaces arising from regular Schr\"odinger operators with a wide range of potentials.


title: Venus as an Analog for Hot Earths

abstract: Here, we evaluate our nearest planetary neighbor, Venus, as an exemplar of the runaway greenhouse state that bounds the inner edge of the habitable zone. Despite its current hellish surface environment, Venus may once have been habitable with oceans of surface liquid water. Over time, it lost its potentially clement environment as the sun brightened. Today, it represents the end-state of habitable planet evolution, and it therefore provides valuable lessons on habitability as a planetary process. Beyond the solar system, exo-Venus analogs are likely common types of planets, and we likely have already discovered many of Venus' sisters orbiting other stars. Furthermore, our near-future exoplanet detection and characterization methods are biased towards observing them. Therefore, it is instructive to consider what Venus can teach us about exo-Venus analogs. By observing exo-Venus planets of differing ages in in different astrophysical contexts, these distant hot terrestrial worlds may likewise allow us to witness processes that occurred on Venus in the past.


title: Chronos: A Unifying Optimization Framework for Speculative Execution of Deadline-critical MapReduce Jobs

abstract: Meeting desired application deadlines in cloud processing systems such as MapReduce is crucial as the nature of cloud applications is becoming increasingly mission-critical and deadline-sensitive. It has been shown that the execution times of MapReduce jobs are often adversely impacted by a few slow tasks, known as stragglers, which result in high latency and deadline violations. While a number of strategies have been developed in existing work to mitigate stragglers by launching speculative or clone task attempts, none of them provides a quantitative framework that optimizes the speculative execution for offering guaranteed Service Level Agreements (SLAs) to meet application deadlines. In this paper, we bring several speculative scheduling strategies together under a unifying optimization framework, called Chronos, which defines a new metric, Probability of Completion before Deadlines (PoCD), to measure the probability that MapReduce jobs meet their desired deadlines. We systematically analyze PoCD for popular strategies including Clone, Speculative-Restart, and Speculative-Resume, and quantify their PoCD in closed-form. The result illuminates an important tradeoff between PoCD and the cost of speculative execution, measured by the total (virtual) machine time required under different strategies. We propose an optimization problem to jointly optimize PoCD and execution cost in different strategies, and develop an algorithmic solution that is guaranteed to be optimal. Chronos is prototyped on Hadoop MapReduce and evaluated against three baseline strategies using both experiments and trace-driven simulations, achieving 50% net utility increase with up to 80% PoCD and 88% cost improvements.


title: Social Behavior of Drivers and Traffic

abstract: Drivers are agents, they are members of the group of drivers. Human groupings have hierarchical structure. The civilization consists of societies, societies consist of groups, and groups consist of individuals. We will consider the group of drivers. Every social group has power to organize individuals and use them for its own purposes. This is true also for the group of drivers. An agent, a driver, we call a member of a social group of drivers, we assume that his-her special properties are defined: they are agents driving a car. A driver-s behavior is influenced to some degree by the need to associate with other drivers and to obtain the approval of other drivers in the group, this is a general property of agents in a group. A driver equates his-her needs with those of the other drivers from the group. From this we explain how three empirically observed dependencies of personal driver radius dependence on some factors enabled us to identify quantity which characterizes verification of information by a driver. The smaller personal radius the larger process of verification of verbal and non-verbal information about the other person in other to accept the risk that this person will be closer to us as a person (driver). We expect that this conclusion is a general conclusion.


title: Accelerating Human-in-the-loop Machine Learning: Challenges and Opportunities

abstract: Development of machine learning (ML) workflows is a tedious process of iterative experimentation: developers repeatedly make changes to workflows until the desired accuracy is attained. We describe our vision for a "human-in-the-loop" ML system that accelerates this process: by intelligently tracking changes and intermediate results over time, such a system can enable rapid iteration, quick responsive feedback, introspection and debugging, and background execution and automation. We finally describe Helix, our preliminary attempt at such a system that has already led to speedups of up to 10x on typical iterative workflows against competing systems.


title: Ideal polyhedral surfaces in Fuchsian manifolds

abstract: Let $S_{g,n}$ be a surface of genus $g > 1$ with $n>0$ punctures equipped with a complete hyperbolic cusp metric. Then it can be uniquely realized as the boundary metric of an ideal Fuchsian polyhedron. In the present paper we give a new variational proof of this result. We also give an alternative proof of the existence and uniqueness of a hyperbolic polyhedral metric with prescribed curvature in a given conformal class.


title: The GALAH survey: Co-orbiting stars and chemical tagging

abstract: We present a study using the second data release of the GALAH survey of stellar parameters and elemental abundances of 15 pairs of stars identified by Oh et al 2017. They identified these pairs as potentially co-moving pairs using proper motions and parallaxes from Gaia DR1. We find that 11 very wide (>1.7 pc) pairs of stars do in fact have similar Galactic orbits, while a further four claimed co-moving pairs are not truly co-orbiting. Eight of the 11 co-orbiting pairs have reliable stellar parameters and abundances, and we find that three of those are quite similar in their abundance patterns, while five have significant [Fe/H] differences. For the latter, this indicates that they could be co-orbiting because of the general dynamical coldness of the thin disc, or perhaps resonances induced by the Galaxy, rather than a shared formation site. Stars such as these, wide binaries, debris of past star formation episodes, and coincidental co-orbiters, are crucial for exploring the limits of chemical tagging in the Milky Way.


title: PCMSolver: an Open-Source Library for Solvation Modeling

abstract: PCMSolver is an open-source library for continuum electrostatic solvation. It can be combined with any quantum chemistry code and requires a minimal interface with the host program, greatly reducing programming effort. As input, PCMSolver needs only the molecular geometry to generate the cavity and the expectation value of the molecular electrostatic potential on the cavity surface. It then returns the solvent polarization back to the host program. The design is powerful and versatile: minimal loss of performance is expected, and a standard single point self-consistent field implementation requires no more than 2 days of work. We provide a brief theoretical overview, followed by two tutorials: one aimed at quantum chemistry program developers wanting to interface their code with PCMSolver, the other aimed at contributors to the library. We finally illustrate past and ongoing work, showing the library's features, combined with several quantum chemistry programs.


title: High-throughput high-resolution multi-scale three dimensional laser printing

abstract: Throughput and resolution are both of critical importance for modern manufacturing, however, they are usually contradictory with each other. Here, a high-throughput high-resolution three dimensional (3D) printing is demonstrated by incorporating a simultaneous spatiotemporal focusing (SSTF) scheme in two-photon polymerization (TPP). Remarkably, the SSTF can ensure generation of a spherical focal spot of which the size depends linearly on the laser power. Thus, the resolution can be continuously adjusted from sub-10 {\mu}m to ~40 {\mu}m by only increasing the laser power. A multi-scale 3D structure is fabricated using this technique, for which the regions of coarse and fine feature sizes are produced at high and low resolutions, respectively.


title: The Nucleus of Active Asteroid 311P/(2013 P5) PANSTARRS

abstract: The unique inner-belt asteroid 311P/PANSTARRS (formerly P/2013 P5) is notable for its sporadic, comet-like ejection of dust in nine distinct epochs spread over $\sim$250 days in 2013. This curious behavior has been interpreted as the product of localized, equator-ward landsliding from the surface of an asteroid rotating at the brink of instability. We obtained new Hubble Space Telescope observations to directly measure the nucleus and to search for evidence of its rapid rotation. However, instead of providing photometric evidence for rapid nucleus rotation, our data set a lower limit to the lightcurve period, $P \ge$ 5.4 hour. The dominant feature of the lightcurve is a V-shaped minimum, $\sim$0.3 magnitudes deep, that is suggestive of an eclipsing binary. Under this interpretation, the time-series data are consistent with a secondary/primary mass ratio, $m_s/m_p \sim$ 1:6, a ratio of separation/primary radius, $r/r_p \sim$ 4 and an orbit period $\sim$0.8 days. These properties lie within the range of other asteroid binaries that are thought to be formed by rotational breakup. While the lightcurve period is long, centripetal dust ejection is still possible if one or both components rotates rapidly ($\lesssim$ 2 hour) and has a small lightcurve variation because of azimuthal symmetry. Indeed, radar observations of asteroids in critical rotation reveal "muffin-shaped" morphologies which are closely azimuthally symmetric and which show minimal lightcurves. Our data are consistent with 311P being a close binary in which one or both components rotates near the centripetal limit. The mass loss in 2013 suggests that breakup occurred recently and could even be on-going. A search for fragments that might have been recently ejected beyond the Hill sphere reveals none larger than effective radius $r_e \sim$ 10 m.


title: Characterization of OCS-HCCCCH and N2O-HCCCCH dimers: Theory and experiment

abstract: The infrared spectra of the weakly-bound dimers OCS-HCCCCH, in the region of the nu1 fundamental band of OCS (2050 cm-1), and N2O-HCCCCH, in the region of the nu1 fundamental band of N2O (2200 cm-1), were observed in a pulsed supersonic slit jet expansion probed with tunable diode/QCL lasers. Both OCS-HCCCCH and N2O-HCCCCH were found to have planar structure with side-by-side monomer units having nearly parallel axes. These bands have hybrid rotational structure which allows for estimates of the orientation of OCS and N2O in the plane of their respective dimers. Analogous bands for OCS-DCCCCD and N2O-DCCCCD were also observed and found to be consistent with the normal isotopologues. Various levels of ab initio calculations were performed to find stationary points on the potential energy surface, optimized structures and interaction energies. Four stable geometries were found for OCS-HCCCCH and three for N2O-HCCCCH. The rotational parameters at CCSD(T*)-F12c level of theory give results in very good agreement with those obtained from the observed spectra. In both dimers, the experimental structure corresponds to the lowest energy isomer.


title: Symplectic Geometry and Connectivity of Spaces of Frames

abstract: Frames provide redundant, stable representations of data which have important applications in signal processing. We introduce a connection between symplectic geometry and frame theory and show that many important classes of frames have natural symplectic descriptions. Symplectic tools seem well-adapted to addressing a number of important questions about frames; in this paper we focus on the frame homotopy conjecture posed in 2002 and recently proved by Cahill, Mixon, and Strawn, which says that the space of finite unit norm tight frames is connected. We give a simple symplectic proof of a double generalization of the frame homotopy conjecture, showing that spaces of complex frames with arbitrary prescribed norms and frame operators are connected. To spark further investigation, we also suggest a number of fundamental questions in frame theory which seem amenable to a symplectic approach.


title: The GALAH and TESS-HERMES surveys: high-resolution spectroscopy of luminous supergiants in the Magellanic Clouds and Bridge

abstract: We report the serendipitous observations of 571 luminous supergiants in the Magellanic Clouds by the spectroscopic GALAH and TESS-HERMES surveys: 434 stars in the Large Magellanic Cloud and 137 in the Small Magellanic Cloud. We also find one star that appears associated with structured star formation in the Magellanic Bridge. Both of these surveys are aimed at the local volume of the Galaxy but have simple, magnitude-limited selection functions that mean they include some observations of luminous extra-Galactic stars. The surveys determine stellar parameter and abundances using The Cannon, a data-driven generative modelling approach. In this work, we explore the results from The Cannon when it is fed the spectra of these intrinsically luminous supergiants in the Magellanic Clouds, which are well outside the normal bounds of The Cannon's training set. We find that, although the parameters are astrophysically incorrect, the $v\sin i$ and the abundances of lithium, barium, and magnesium are excellent discriminants of these stars. It shows that in the future, with an expanded training set, it should be possible to determine accurate values for these types of stars.


title: Investigating Cybersecurity Issues In Active Traffic Management Systems

abstract: Active Traffic Management (ATM) systems have been introduced by transportation agencies to manage recurrent and non-recurrent congestion. ATM systems rely on the interconnectivity of components made possible by wired and/or wireless networks. Unfortunately, this connectivity that supports ATM systems also provides potential system access points that results in vulnerability to cyberattacks. This is becoming more pronounced as ATM systems begin to integrate internet of things (IoT) devices. Hence, there is a need to rigorously evaluate ATM systems for cyberattack vulnerabilities, and explore design concepts that provide stability and graceful degradation in the face of cyberattacks. In this research, a prototype ATM system along with a real-time cyberattack monitoring system were developed for a 1.5-mile section of I-66 in Northern Virginia. The monitoring system detects deviation from expected operation of an ATM system by comparing lane control states generated by the ATM system with lane control states deemed most likely by the monitoring system. In case of any deviation between two sets of states, the monitoring system displays the lane control states generated by the back-up data source. In a simulation experiment, the prototype ATM system and cyberattack monitoring system were subject to emulated cyberattacks. The evaluation results showed that when subject to cyberattack, the mean speed reduced by 15% compared to the case with the ATM system and was similar to the baseline case. This illustrates that the effectiveness of the ATM system was negated by cyberattacks. The monitoring system however, allowed the ATM system to revert to an expected safe state and reduced the negative impact of cyberattacks. These results illustrate the need to revisit ATM system design concepts as a means to protect against cyberattacks in addition to traditional system intrusion prevention approaches.


title: Densely Connected High Order Residual Network for Single Frame Image Super Resolution

abstract: Deep convolutional neural networks (DCNN) have been widely adopted for research on super resolution recently, however previous work focused mainly on stacking as many layers as possible in their model, in this paper, we present a new perspective regarding to image restoration problems that we can construct the neural network model reflecting the physical significance of the image restoration process, that is, embedding the a priori knowledge of image restoration directly into the structure of our neural network model, we employed a symmetric non-linear colorspace, the sigmoidal transfer, to replace traditional transfers such as, sRGB, Rec.709, which are asymmetric non-linear colorspaces, we also propose a "reuse plus patch" method to deal with super resolution of different scaling factors, our proposed methods and model show generally superior performance over previous work even though our model was only roughly trained and could still be underfitting the training set.


title: Principal Component Analysis as a Tool for Characterizing Black Hole Images and Variability

abstract: We explore the use of principal component analysis (PCA) to characterize high-fidelity simulations and interferometric observations of the millimeter emission that originates near the horizons of accreting black holes. We show mathematically that the Fourier transforms of eigenimages derived from PCA applied to an ensemble of images in the spatial-domain are identical to the eigenvectors of PCA applied to the ensemble of the Fourier transforms of the images, which suggests that this approach may be applied to modeling the sparse interferometric Fourier-visibilities produced by an array such as the Event Horizon Telescope (EHT). We also show that the simulations in the spatial domain themselves can be compactly represented with a PCA-derived basis of eigenimages allowing for detailed comparisons between variable observations and time-dependent models, as well as for detection of outliers or rare events within a time series of images. Furthermore, we demonstrate that the spectrum of PCA eigenvalues is a diagnostic of the power spectrum of the structure and, hence, of the underlying physical processes in the simulated and observed images.


title: Dynamo Processes Constrained by Solar and Stellar Observations

abstract: Our understanding of stellar dynamos has largely been driven by the phenomena we have observed of our own Sun. Yet, as we amass longer-term datasets for an increasing number of stars, it is clear that there is a wide variety of stellar behavior. Here we briefly review observed trends that place key constraints on the fundamental dynamo operation of solar-type stars to fully convective M dwarfs, including: starspot and sunspot patterns, various magnetism-rotation correlations, and mean field flows such as differential rotation and meridional circulation. We also comment on the current insight that simulations of dynamo action and flux emergence lend to our working knowledge of stellar dynamo theory. While the growing landscape of both observations and simulations of stellar magnetic activity work in tandem to decipher dynamo action, there are still many puzzles that we have yet to fully understand.


title: Tree Morphology for Phenotyping from Semantics-Based Mapping in Orchard Environments

abstract: Measuring tree morphology for phenotyping is an essential but labor-intensive activity in horticulture. Researchers often rely on manual measurements which may not be accurate for example when measuring tree volume. Recent approaches on automating the measurement process rely on LIDAR measurements coupled with high-accuracy GPS. Usually each side of a row is reconstructed independently and then merged using GPS information. Such approaches have two disadvantages: (1) they rely on specialized and expensive equipment, and (2) since the reconstruction process does not simultaneously use information from both sides, side reconstructions may not be accurate. We also show that standard loop closure methods do not necessarily align tree trunks well. In this paper, we present a novel vision system that employs only an RGB-D camera to estimate morphological parameters. A semantics-based mapping algorithm merges the two-sides 3D models of tree rows, where integrated semantic information is obtained and refined by robust fitting algorithms. We focus on measuring tree height, canopy volume and trunk diameter from the optimized 3D model. Experiments conducted in real orchards quantitatively demonstrate the accuracy of our method.


title: An information-theoretic on-line update principle for perception-action coupling

abstract: Inspired by findings of sensorimotor coupling in humans and animals, there has recently been a growing interest in the interaction between action and perception in robotic systems [Bogh et al., 2016]. Here we consider perception and action as two serial information channels with limited information-processing capacity. We follow [Genewein et al., 2015] and formulate a constrained optimization problem that maximizes utility under limited information-processing capacity in the two channels. As a solution we obtain an optimal perceptual channel and an optimal action channel that are coupled such that perceptual information is optimized with respect to downstream processing in the action module. The main novelty of this study is that we propose an online optimization procedure to find bounded-optimal perception and action channels in parameterized serial perception-action systems. In particular, we implement the perceptual channel as a multi-layer neural network and the action channel as a multinomial distribution. We illustrate our method in a NAO robot simulator with a simplified cup lifting task.


title: Automatic Detection of Indoor and Outdoor Scenarios using NMEA Message Data from GPS Receivers

abstract: Detection of indoor and outdoor scenarios is an important resource for many types of activities such as multisensor navigation and location-based services. This research presents the use of NMEA data provided by GPS receivers to characterize different types of scenarios automatically. A set of static tests was performed to evaluate metrics such as number of satellites, positioning solution geometry and carrier-to-receiver noise-density ratio values to detect possible patterns to determine indoor and outdoor scenarios. Subsequently, validation tests are applied to verify that parameters obtained are adequate.


title: Persistence probability of a random polynomial arising from evolutionary game theory

abstract: In this paper, we obtain an asymptotic formula for the persistence probability in the positive real line of a random polynomial arising from evolutionary game theory. It corresponds to the probability that a multi-player two-strategy random evolutionary game has no internal equilibria. The key ingredient is to approximate the sequence of random polynomials indexed by their degrees by an appropriate centered stationary Gaussian process.

categories: math.AP math-ph math.DS math.MP math.PR q-bio.PE


title: Powell moves and the Goeritz group

abstract: In 1980 J. Powell proposed that five specific elements sufficed to generate the Goeritz group of any Heegaard splitting of $S^3$, extending work of Goeritz on genus $2$ splittings. Here we prove that Powell's conjecture was correct for splittings of genus $3$ as well, and discuss a framework for deciding the truth of the conjecture for higher genus splittings.


title: White-light hyperbolic Airy beams

abstract: Ultra-broadband hyperbolic Airy beams are experimentally realized and spatio-temporally characterized. Transmissive (refractive) and reflective cubic phase elements were used to impose a cubic phase on the input beam. Nondispersing beams are produced in reflective geometry, while the main lobe of the hyperbolic Airy beam created with a transmissive refractive phase element suffered from lateral dispersion.


title: Long-term optical, UV, and X-ray continuum variations in the changing-look AGN HE 1136-2304

abstract: A strong outburst in the X-ray continuum and a change of its Seyfert spectral type was detected in HE 1136-2304 in 2014. The spectral type changed from nearly Seyfert 2 type (1.95) to Seyfert 1.5 type in comparison to previous observations taken ten to twenty years before. In a subsequent variability campaign we wanted to investigate whether this outburst was a single event or whether the variability pattern following the outburst was similar to those seen in other variable Seyfert galaxies. In addition to a SALT spectral variability campaign, we carried out optical continuum as well as X-ray and UV (Swift) monitoring studies from 2014 to 2017. HE 1136-2304 strongly varied on timescales of days to months from 2014 to 2017. No systematic trends were found in the variability behavior following the outburst in 2014. A general decrease in flux would have been expected for a tidal disruption event. This could not be confirmed. More likely the flux variations are connected to irregular fluctuations in the accretion rate. The strongest variability amplitudes have been found in the X-ray regime: HE 1136-2304 varied by a factor of eight during 2015. The amplitudes of the continuum variability (from the UV to the optical) systematically decreased with wavelength following a power law F_var = a ${\times}$ {\lambda}^-c with c = 0.84. There is a trend that the B-band continuum shows a delay of three light days with respect to the variable X-ray flux. The Seyfert type 1.5 did not change despite the strong continuum variations for the period between 2014 and 2017.


title: Efimov Physics in Quenched Unitary Bose Gases

abstract: We study the impact of three-body physics in quenched unitary Bose gases, focusing on the role of the Efimov effect. Using a local density model, we solve the three-body problem and determine three-body decay rates at unitary, finding density-dependent, log-periodic Efimov oscillations, violating the expected continuous scale-invariance in the system. We find that the breakdown of continuous scale-invariance, due to Efimov physics, manifests also in the earliest stages of evolution after the interaction quench to unitarity, where we find the growth of a substantial population of Efimov states for densities in which the interparticle distance is comparable to the size of an Efimov state. This agrees with the early-time dynamical growth of three-body correlations at unitarity [Colussi {\em et al}., Phys. Rev. Lett. 120, 100401 (2018)]. By varying the sweep rate away from unitarity, we also find a departure from the usual Landau-Zener analysis for state transfer when the system is allowed to evolve at unitarity and develop correlations.


title: Weak* and entropy approximation of nonhyperbolic measures: a geometrical approach

abstract: We study $C^1$-robustly transitive and nonhyperbolic diffeomorphisms having a partially hyperbolic splitting with one-dimensional central bundle whose strong un-/stable foliations are both minimal. {In dimension $3$, an important class of examples of such systems is given by those with a simple closed periodic curve tangent to the central bundle.} We prove that there is a $C^1$-open and dense subset of such diffeomorphisms such that every nonhyperbolic ergodic measure (i.e. with zero central exponent) can be approximated in the weak$\ast$ topology and in entropy by measures supported in basic sets with positive (negative) central Lyapunov exponent. Our method also allows to show how entropy changes across measures with central Lyapunov exponent close to zero. We also prove that any nonhyperbolic ergodic measure is in the intersection of the convex hulls of the measures with positive central exponent and with negative central exponent.


title: On the radial acceleration relation for galaxies in a wide range of redshifts

abstract: The thermal history of the Universe is introduced within the Debye Gravitational Theory (DGT), a thermodynamic theory of induced gravity, and allows to obtain the evolution of systems with the redshift. DGT reproduce the ESO VLT observations, showing falling outer rotation curves for galaxies at redshift above 0.77. A scaling law is observed in the radial acceleration relation (RAR) of galaxies. For accelerations smaller than $\sim 10^{-10}ms^{-1}$ the observed acceleration $g_{obs}$ decreases more slowly than acceleration generated by the baryonic mass $g_{bar}$, following always the relation $g_{obs} \sim \sqrt{g_{bar}}$. The RAR does not care about the specific properties of the galaxy, the relation exists in nearby high-mass elliptical and low-mass spheroidal galaxies. In this paper, through a straightforward analysis, we show that according to DGT, the RAR scaling law observed in nearby galaxies is broken, when are considered distant galaxies. The extreme case happen for galaxies at redshift above 0.77, because according to DGT, the galaxies have declining rotation curves, and the rotation velocity falling faster than the Keplerian-law curve. Then in this case $g_{obs}$ decreases more faster than $g_{bar}$, this means that $g_{obs}$ is systematically lower than $g_{bar}$. We show that DGT prediction for the RAR of galaxies at high redshift is in agreement from those obtained from the falling rotation curves observed by VLT telescope.


title: Trace class Markov chains for the Normal-Gamma Bayesian shrinkage model

abstract: High-dimensional data, where the number of variables exceeds or is comparable to the sample size, is now pervasive in many scientific applications. In recent years, Bayesian shrinkage models have been developed as effective and computationally feasible tools to analyze such data, especially in the context of linear regression. In this paper, we focus on the Normal-Gamma shrinkage model developed by Griffin and Brown. This model subsumes the popular Bayesian lasso model, and a three-block Gibbs sampling algorithm to sample from the resulting intractable posterior distribution has been developed by Griffin and Brown. We consider an alternative two-block Gibbs sampling algorithm and rigorously demonstrate its advantage over the three-block sampler by comparing specific spectral properties. In particular, we show that the Markov operator corresponding to the two-block sampler is trace class (and hence Hilbert-Schmidt), whereas the operator corresponding to the three-block sampler is not even Hilbert-Schmidt. The trace class property for the two-block sampler implies geometric convergence for the associated Markov chain, which justifies the use of Markov chain CLT's to obtain practical error bounds for MCMC based estimates. Additionally, it facilitates theoretical comparisons of the two-block sampler with sandwich algorithms which aim to improve performance by inserting inexpensive extra steps in between the two conditional draws of the two-block sampler.