title: Magpy: A C++ accelerated Python package for simulating magnetic nanoparticle stochastic dynamics

abstract: Magpy is a C++ accelerated Python package for modelling and simulating the magnetic dynamics of nano-sized particles. Nanoparticles are modelled as a system of three-dimensional macrospins and simulated with a set of coupled stochastic differential equations (the Landau-Lifshitz-Gilbert equation), which are solved numerically using explicit or implicit methods. The results of the simulations may be used to compute equilibrium states, the dynamic response to external magnetic fields, and heat dissipation. Magpy is built on a C++ library, which is optimised for serial execution, and exposed through a Python interface utilising an embarrassingly parallel strategy. Magpy is free, open-source, and available on github under the 3-Clause BSD License.


title: H.E.S.S. highlights

abstract: The H.E.S.S. collaboration continues to run an array of five Imaging Atmospheric Cherenkov Telescopes to observe the Southern sky in very high energy gamma-rays. In this presentation, recent highlight results obtained with H.E.S.S. are briefly reviewed, with references to the relevant publications.


title: Detection of the Characteristic Pion-Decay Signature in the Molecular Clouds

abstract: The $\gamma$-ray emission from molecular clouds is widely believed to have a hadronic origin, while unequivocal evidence is still lacking. In this work, we analyze the Fermi-LAT Pass 8 publicly available data accumulated from 2008 August 4 to 2017 June 29 and report the significant detection of a characteristic $\pi^0$-decay feature in the $\gamma$-ray spectra of the molecular clouds Orion A and Orion B. This detection provides a direct evidence for the hadronic origin of their $\gamma$-ray emission.


title: Commutativity in Lagrangian and Hamiltonian Mechanics

abstract: The main result of this note is a characterization of the Poisson commutativity of Hamilton functions in terms of their principal action functions.


title: The QLBS Q-Learner Goes NuQLear: Fitted Q Iteration, Inverse RL, and Option Portfolios

abstract: The QLBS model is a discrete-time option hedging and pricing model that is based on Dynamic Programming (DP) and Reinforcement Learning (RL). It combines the famous Q-Learning method for RL with the Black-Scholes (-Merton) model's idea of reducing the problem of option pricing and hedging to the problem of optimal rebalancing of a dynamic replicating portfolio for the option, which is made of a stock and cash. Here we expand on several NuQLear (Numerical Q-Learning) topics with the QLBS model. First, we investigate the performance of Fitted Q Iteration for a RL (data-driven) solution to the model, and benchmark it versus a DP (model-based) solution, as well as versus the BSM model. Second, we develop an Inverse Reinforcement Learning (IRL) setting for the model, where we only observe prices and actions (re-hedges) taken by a trader, but not rewards. Third, we outline how the QLBS model can be used for pricing portfolios of options, rather than a single option in isolation, thus providing its own, data-driven and model independent solution to the (in)famous volatility smile problem of the Black-Scholes model.


title: Noncrossing partitions, Bruhat order and the cluster complex

abstract: We introduce two order relations on finite Coxeter groups which refine the absolute and the Bruhat order, and establish some of their main properties. In particular we study the restriction of these orders to noncrossing partitions and show that the intervals for these orders can be enumerated in terms of the cluster complex. The properties of our orders permit to revisit several results in Coxeter combinatorics, such as the Chapoton triangles and how they are related, the enumeration of reflections with full support, the bijections between clusters and noncrossing partitions.


title: Navigable maps of structural brain networks across species

abstract: Connectomes are spatially embedded networks whose architecture has been shaped by physical constraints and communication needs throughout evolution. Using a decentralized navigation protocol, we investigate the relationship between the structure of the connectomes of different species and their spatial layout. As a navigation strategy, we use greedy routing where nearest neighbors, in terms of geometric distance, are visited. We measure the fraction of successful greedy paths and their length as compared to shortest paths in the topology of connectomes. In Euclidean space, we find a striking difference between the navigability properties of mammalian and non-mammalian species, which implies the inability of Euclidean distances to fully explain the structural organization of their connectomes. In contrast, we find that hyperbolic space, the effective geometry of complex networks, provides almost perfectly navigable maps of connectomes for all species, meaning that hyperbolic distances are exceptionally congruent with the structure of connectomes. Hyperbolic maps therefore offer a quantitative meaningful representation of connectomes that suggests a new cartography of the brain based on the combination of its connectivity with its effective geometry rather than on its anatomy only. Hyperbolic maps also provide a universal framework to study decentralized communication processes in connectomes of different species and at different scales on an equal footing.


title: A Precise Distance to the Host Galaxy of the Binary Neutron Star Merger GW170817 Using Surface Brightness Fluctuations

abstract: The joint detection of gravitational waves and electromagnetic radiation from the binary neutron star (BNS) merger GW170817 has provided unprecedented insight into a wide range of physical processes: heavy element synthesis via the $r$-process; the production of relativistic ejecta; the equation of state of neutron stars and the nature of the merger remnant; the binary coalescence timescale; and a measurement of the Hubble constant via the "standard siren" technique. In detail, all of these results depend on the distance to the host galaxy of the merger event, NGC4993. In this paper we measure the surface brightness fluctuation (SBF) distance to NGC4993 in the F110W and F160W passbands of the Wide Field Camera 3 Infrared Channel on the Hubble Space Telescope (HST). For the preferred F110W passband we derive a distance modulus of $m{-}M=33.05\pm0.08\pm0.10$ mag, or a linear distance $d=40.7\pm1.4\pm1.9$ Mpc (random and systematic errors, respectively); a virtually identical result is obtained from the F160W data. This is the most precise distance to NGC4993 available to date. Combining our distance measurement with the corrected recession velocity of NGC4993 implies a Hubble constant $H_0=71.9\pm 7.1~km~s^{-1}~Mpc^{-1}$. A comparison of our result to the GW-inferred value of $H_0$ indicates a binary orbital inclination of $i\,{\gtrsim}\,137~\deg$. The SBF technique can be applied to early-type host galaxies of BNS mergers to ${\sim\,}100$ Mpc with HST and possibly as far as ${\sim\,}300$ Mpc with the James Webb Space Telescope, thereby helping to break the inherent distance-inclination degeneracy of the GW signals at distances where many future BNS mergers are likely to be detected.


title: Universes without the Weak Force: Astrophysical Processes with Stable Neutrons

abstract: We investigate a class of universes in which the weak interaction is not in operation. We consider how astrophysical processes are altered in the absence of weak forces, including Big Bang Nucleosynthesis (BBN), galaxy formation, molecular cloud assembly, star formation, and stellar evolution. Without weak interactions, neutrons no longer decay, and the universe emerges from its early epochs with a mixture of protons, neutrons, deuterium, and helium. The baryon-to-photon ratio must be smaller than the canonical value in our universe to allow free nucleons to survive the BBN epoch without being incorporated into heavier nuclei. At later times, the free neutrons readily combine with protons to make deuterium in sufficiently dense parts of the interstellar medium, and provide a power source before they are incorporated into stars. Almost all of the neutrons are incorporated into deuterium nuclei before stars are formed. As a result, stellar evolution proceeds primarily through strong interactions, with deuterium first burning into helium, and then helium fusing into carbon. Low-mass deuterium-burning stars can be long-lived, and higher mass stars can synthesize the heavier elements necessary for life. Although somewhat different from our own, such universes remain potentially habitable.


title: Toward Stronger Robustness of Network Controllability: A Snapback Network Model

abstract: A new complex network model, called q-snapback network, is introduced. Basic topological characteristics of the network, such as degree distribution, average path length, clustering coefficient and Pearson correlation coefficient, are evaluated. The typical 4-motifs of the network are simulated. The robustness of both state and structural controllabilities of the network against targeted and random node- and edge-removal attacks, with comparisons to the multiplex congruence network and the generic scale-free network, are presented. It is shown that the q-snapback network has the strongest robustness of controllabilities due to its advantageous inherent structure with many chain- and loop-motifs.


title: The $q$-Onsager Algebra and the Universal Askey-Wilson Algebra

abstract: Recently Pascal Baseilhac and Stefan Kolb obtained a PBW basis for the $q$-Onsager algebra $\mathcal O_q$. They defined the PBW basis elements recursively, and it is obscure how to express them in closed form. To mitigate the difficulty, we bring in the universal Askey-Wilson algebra $\Delta_q$. There is a natural algebra homomorphism $\natural \colon \mathcal O_q \to \Delta_q$. We apply $\natural $ to the above PBW basis, and express the images in closed form. Our results make heavy use of the Chebyshev polynomials of the second kind.


title: R*-operation and five-loop calculations

abstract: We sketch how the R*-operation can be used to compute the pole terms of Feynman diagrams. We identify computational difficulties when performing five-loop calculations, and provide four solutions that drastically reduce the number of terms that are generated. Using these methods, we have computed the beta function for Yang-Mills theory with fermions, the R-ratio in electron-positron annihilation, and Higgs decays to quarks and gluons at five-loop accuracy. The results for the beta function and Higgs decay width to gluons in the heavy-top limit are briefly discussed. There is no need for six-loop extensions of these calculations in the near future.


title: Four-loop results on anomalous dimensions and splitting functions in QCD

abstract: We report on recent progress on the flavour non-singlet splitting functions in perturbative QCD. The~exact four-loop (N^3LO) contribution to these functions has been obtained in the planar limit of a large number of colours. Phenomenologically sufficient approximate expressions have been obtained for the parts not exactly known so far. Both cases include results for the four-loop cusp and virtual anomalous dimensions which are relevant well beyond the evolution of non-singlet quark distributions, for which an accuracy of (well) below 1% has now been been reached.


title: Phononic Rogue Waves

abstract: We present a theoretical study of extreme events occurring in phononic lattices. In particular, we focus on the formation of rogue or freak waves, which are characterized by their localization in both spatial and temporal domains. We consider two examples. The first one is the prototypical nonlinear mass-spring system in the form of a homogeneous Fermi-Pasta-Ulam-Tsingou (FPUT) lattice with a polynomial potential. By deriving an approximation based on the nonlinear Schroedinger (NLS) equation, we are able to initialize the FPUT model using a suitably transformed Peregrine soliton solution of the NLS, obtaining dynamics that resembles a rogue wave on the FPUT lattice. We also show that Gaussian initial data can lead to dynamics featuring rogue wave for sufficiently wide Gaussians. The second example is a diatomic granular crystal exhibiting rogue wave like dynamics, which we also obtain through an NLS reduction and numerical simulations. The granular crystal (a chain of particles that interact elastically) is a widely studied system that lends itself to experimental studies. This study serves to illustrate the potential of such dynamical lattices towards the experimental observation of acoustic rogue waves.


title: Spectroscopic and photometric confirmation of chromospheric activity in four stars

abstract: We present analysis of medium resolution optical spectra and long term V band photometry of four cool stars, BD+13 5000, BD+11 3024, TYC 3557-919-1 and TYC 5163-1764-1. Our spectroscopic analysis reveals that the stars are giant or sub-giant from K0 or K1 spectral type, and all of them exhibit emission features in their Ca ii H& K lines. These features appear to be modulated with the rotation of the stars. Except BD+11 3024, we observe that the radial velocities of the target stars are not stable, which suggests that each of them might be a member of a binary system. Global analysis of photometric data indicates clear cyclic variation for BD+13 5000 and TYC 5163-1764-1 with a period of 8.0$\pm$0.3 and 5.04$\pm$0.04 year, respectively. Besides that, we observe a dramatic increase ($\sim$0\fm7) in the mean brightness of BD+11 3024, accompanied with a 2.87$\pm$0.12 cyclic variation, embedded into the global brightening trend, which indicates possible multiple cycles on this star.


title: On the vanishing contact structure for viscosity solutions of contact type Hamilton-Jacobi equations I: Cauchy problem

abstract: We study the representation formulae for the fundamental solutions and viscosity solutions of the Hamilton-Jacobi equations of contact type. We also obtain a vanishing contact structure result for relevant Cauchy problems which can be regarded as an extension to the vanishing discount problem.


title: Level Reciprocity in the twisted second moment of Rankin-Selberg L-functions

abstract: We prove an exact formula for the second moment of Rankin-Selberg $L$-functions $L(1/2,f \times g)$ twisted by $\lambda_f(p)$, where $g$ is a fixed holomorphic cusp form and $f$ is summed over automorphic forms of a given level $q$. The formula is a reciprocity relation that exchanges the twist parameter $p$ and the level $q$. The method involves the Bruggeman/Kuznetsov trace formula on both ends; finally the reciprocity relation is established by an identity of sums of Kloosterman sums.


title: Observable windows for the QCD axion through $N_\text{eff}$

abstract: We show that when the QCD axion is directly coupled to quarks with $c_q/f \, \partial_\mu a \, \bar{q} \gamma^\mu \gamma^5 q$, such as in DFSZ models, the dominant production mechanism in the early universe at temperatures $1 \, {\rm GeV}\lesssim T \lesssim 100 \,{\rm GeV}$ is obtained via $q \bar{q} \leftrightarrow g a$ and $q g \leftrightarrow q a$, where $g$ are gluons. Different heavy quarks $q_i$ can produce a thermal axion background that decouples at a temperature $T_i$: (1) top quark at $T_t \lesssim 100 \, {\rm GeV}$ for $f/c_t \lesssim 3\times 10^8 {\rm GeV}$; (2) bottom quark at $T_b \lesssim m_b$, for $f/c_b\lesssim 8\times 10^{7} {\rm GeV}$; (3) charm quark at $T_c \lesssim m_c $ for $f/c_c \lesssim 5\times 10^{7} {\rm GeV}$. Each of these cases corresponds to a contribution to the effective number of relativistic degrees of freedom, in the windows given by $0.027 \leq \Delta N_\text{eff}\leq 0.031$, $0.037 \leq \Delta N_\text{eff} \leq 0.039$ and $0.039 \leq \Delta N_\text{eff}$, respectively. These contributions are larger than the one obtained when thermalization happens only above the electroweak phase transition, $\Delta N_{\rm eff}\lesssim 0.027$, and are within reach of future CMB S4 experiments, thus opening an alternative window to detect the axion and to test the early universe at such temperatures.


title: Non-universal behaviour of helical two-dimensional three-component turbulence

abstract: The dynamics of two-dimensional three-component (2D3C) flows is relevant to describe the long-time evolution of strongly rotating flows and/or of conducting fluids with a strong mean magnetic field. We show that in the presence of a strong helical forcing, the out-of-plane component ceases to behave as a passive advected quantity and develops a nontrivial dynamics which deeply changes its large-scale properties. We show that a small-scale helicity injection correlates the input on the 2D component with the one on the out-of-plane component. As a result, the third component develops a non-trivial energy transfer. The latter is mediated by homochiral triads, confirming the strong 3D nature of the leading dynamical interactions. In conclusion, we show that the out-of-plane component in a 2D3C flow enjoys strong non-universal properties as a function of the degree of mirror symmetry of the small-scale forcing.


title: EUV Irradiance Inputs to Thermospheric Density Models: Open Issues and Path Forward

abstract: One of the objectives of the NASA LWS Institute on Nowcasting of Atmospheric Drag for LEO Spacecraft was to investigate whether and how to increase the accuracy of atmospheric drag models by improving the quality of the solar forcing inputs, namely Extreme Ultraviolet (EUV) irradiance information. In this focused review, we examine the status of and issues with EUV measurements and proxies, discuss recent promising developments, and suggest a number of ways to improve the reliability, availability, and forecast accuracy of EUV measurements in the next solar cycle.


title: Combined analysis of double Higgs production via gluon fusion at the HL-LHC in the effective field theory approach

abstract: We perform the combined analysis of the double Higgs production via gluon fusion in the $b\bar{b} \gamma\gamma$ and $b\bar{b}\tau^+\tau^-$ decay channels at the High-Luminosity LHC (HL-LHC). To validate our analysis, we reproduce the ATLAS result of the $b\bar{b} \gamma\gamma$ process including all contributions from fakes. For the $b\bar{b}\tau^+\tau^-$ decay channel, we perform the similar analysis to the CMS one. As an improvement, we also perform the multivariate analysis employing the boosted decision tree algorithm. Then, we derive 68% probability contours on anomalous Higgs couplings in the effective field theory (EFT) approach for various analyses. We find that the $b\bar{b}\tau^+\tau^-$ process outperforms the $b\bar{b}\gamma\gamma$ for the measurement of energy-growing operators, while adding the $b\bar{b}\tau^+\tau^-$ process is least beneficial for improving the precision of the Higgs self-coupling (mainly set by the $b\bar{b}\gamma\gamma$ process). We illustrate that the double Higgs production alone can be comparable to the single Higgs process in constraining the modification of the top Yukawa coupling in the positive direction. Focusing on the Higgs self-coupling as a special interest, we derive the precision as a function of various improvable parameters such as tag and mistag rates of tau leptons, heavy flavor jets, photon identification, diphoton mass resolution, and jet energy resolution to take into account future phenomenological studies. As an optimistic benchmark scenario, we illustrate that the 68% and 95% probability intervals of the Higgs self-coupling, $\lambda_3/\lambda_{3}^{SM}$, at the HL-LHC can reach $[0.2,\, 2.3]$ and $[-0.1,\, 3.5] \cup [4.0,\, 6.5]$, respectively, where the correlation among the EFT coefficients is taken into account.


title: Satellites Form Fast & Late: a Population Synthesis for the Galilean Moons

abstract: The satellites of Jupiter are thought to form in a circumplanetary disc. Here we address their formation and orbital evolution with a population synthesis approach, by varying the dust-to-gas ratio, the disc dispersal timescale and the dust refilling timescale. The circumplanetary disc initial conditions (density and temperature) are directly drawn from the results of 3D radiative hydrodynamical simulations. The disc evolution is taken into account within the population synthesis. The satellitesimals were assumed to grow via streaming instability. We find that the moons form fast, often within $10^4$ years, due to the short orbital timescales in the circumplanetary disc. They form in sequence, and many are lost into the planet due to fast type I migration, polluting Jupiter's envelope with typically 15 Earth-masses of metals. The last generation of moons can form very late in the evolution of the giant planet, when the disc has already lost more than the 99% of its mass. The late circumplanetary disc is cold enough to sustain water ice, hence not surprisingly the 85% of the moon population has icy composition. The distribution of the satellite-masses is peaking slightly above Galilean masses, up until a few Earth-masses, in a regime which is observable with the current instrumentation around Jupiter-analog exoplanets orbiting sufficiently close to their host stars. We also find that systems with Galilean-like masses occur in 20% of the cases and they are more likely when discs have long dispersion timescales and high dust-to-gas ratios.


title: A Gauss-Jacobi Kernel Compression Scheme for Fractional Differential Equations

abstract: A scheme for approximating the kernel $w$ of the fractional $\alpha$-integral by a linear combination of exponentials is proposed and studied. The scheme is based on the application of a composite Gauss-Jacobi quadrature rule to an integral representation of $w$. This results in an approximation of $w$ in an interval $[\delta,T]$, with $0<\delta$, which converges rapidly in the number $J$ of quadrature nodes associated with each interval of the composite rule. Using error analysis for Gauss-Jacobi quadratures for analytic functions, an estimate of the relative pointwise error is obtained. The estimate shows that the number of terms required for the approximation to satisfy a prescribed error tolerance is bounded for all $\alpha\in(0,1)$, and that $J$ is bounded for $\alpha\in(0,1)$, $T>0$, and $\delta\in(0,T)$.


title: Constraining Transfer Cross Sections Using Bayes' Theorem

abstract: Being able to rigorously quantify the uncertainties in reaction models is crucial to moving this field forward. Even though Bayesian methods are becoming increasingly popular in nuclear theory, they are yet to be implemented and applied in reaction theory problems. The purpose of this work is to investigate, using Bayesian methods, the uncertainties in the optical potentials generated from fits to elastic scattering data and the subsequent uncertainties in the transfer predictions. We also study the differences in two reaction models where the parameters are constrained in a similar manner, as well as the impact of reducing the experimental error bars on the data used to constrain the parameters. We use Bayes' Theorem combined with a Markov Chain Monte Carlo to determine posterior distributions for the parameters of the optical model, constrained by neutron-, proton-, and/or deuteron-target elastic scattering. These potentials are then used to predict transfer cross sections within the adiabatic wave approximation or the distorted-wave Born approximation. We study a number of reactions involving deuteron projectiles with energies in the range of $10-25$ MeV/u on targets with mass $A=48-208$. The case of $^{48}$Ca(d,p)$^{49}$Ca transfer to the ground state is described in detail. A comparative study of the effect of the size of experimental errors is also performed. Five transfer reactions are studied, and their results compiled in order to systematically identify trends. Uncertainties in transfer cross sections can vary significantly (25-100\%) depending on the reaction. While these uncertainties are reduced when smaller experimental error bars are used to constrain the potentials, this reduction is not trivially related to the error reduction. We also find smaller uncertainties when using the adiabatic formulation than when using distorted-wave Born approximation.


title: AGN Black hole mass estimates using polarization in broad emission Lines

abstract: The innermost regions in active galactic nuclei (AGNs) were not being spatially resolved so far but spectropolarimetry can provide us insight about their hidden physics and the geometry. From spectropolarimetric observations in broad emission lines and assuming equatorial scattering as a dominant polarization mechanism, it is possible to estimate the mass of supermassive black holes (SMBHs). We explore the possibilities and limits and to put constraints on the usage of the method for determining SMBH masses using polarization in broad emission lines by providing more in-depth theoretical modeling. Methods. We use the Monte Carlo radiative transfer code STOKES for exploring polarization of Type 1 AGNs. We model equatorial scattering using flared-disk geometry for a set of different SMBH masses assuming Thomson scattering. In addition to the Keplerian motion in the BLR, we also consider cases of additional radial inflows and vertical outflows. We model the profiles of polarization plane position angle, degree of polarization and total unpolarized line for different BLR geometries and different SMBH masses. Our modeling confirms that the method can be widely used for Type-1 AGNs when viewing inclinations are between 25 and 45 degrees. We show that the distance between the BLR and scattering region (SR) has a significant impact on the mass estimates and the best mass estimates are when the SR is situated at the distance 1.5-2.5 times larger than the outer BLR radius. Our models show that if Keplerian motion can be traced through the polarized line profile, then the direct estimation of the mass of the SMBH can be performed. When radial inflows or vertical outflows are present in the BLR, this method can be applied if velocities of the inflow/outflow are less than 500 km/s. We find that models for NGC4051, NGC4151, 3C273 and PG0844+349 are in good agreements with observations.


title: $p-$forms non-minimally coupled to gravity in Randall-Sundrum scenarios

abstract: In this paper we study the coupling of $p$-form fields with geometrical tensor fields, namely Ricci, Einstein, Horndeski and Riemann in Randall-Sundrum scenarios with co-dimension one. We consider delta-like and branes generated by a kink and a domain wall. We begin by a detailed study of the Kalb-Ramond (KR) field. The analysis of KR field is very rich since it is a tensorial object and more complex non-minimal couplings are possible. The generalization to $p$-forms can provide more information about the properties and structures that can possibly be universal in the geometrical localization mechanism. The zero mode is treated separately and conditions for localization of zero modes of $p-$forms are found for all the cases above and with this we arrive at the above conclusion about vector fields. Another property that can be tested is the absence of resonances found in the case of vector fields. For this we analyze the possible unstable massive modes for all the above cases via transmission coefficient. Our conclusion is that we have more probability to observe massive unstable modes in the Ricci and Riemann coupling.


title: Revealing compact structures of interstellar plasma in the Galaxy with RadioAstron

abstract: The aim of our work was to study the spatial structure of inhomogeneities of interstellar plasma in the directions of five pulsars: B0823+26, B0834+06, B1237+25, B1929+10, and B2016+28. Observations of these pulsars were made with RadioAstron space-ground radio interferometer at 324 MHz. We measured the angular size of the scattering disks to be in range between 0.63 and 3.2 mas. We determined the position of scattering screens on the line of sight. Independent estimates of the distances to the screens were made from the curvature of parabolic arcs revealed in the secondary spectra of four pulsars. The model of uniform distribution of inhomogeneities on the line of sight is not suitable. According to the results, we came to the conclusion that scattering is mainly produced by compact plasma layers and the uniform model of inhomogeneties distribution on the line of sight in not applicable.


title: Density matrix renormalization group for a highly degenerate quantum system: Sliding environment block approach

abstract: We present an infinite lattice DMRG sweeping procedure which can be used as a replacement for the standard infinite lattice blocking schemes. Although the scheme is generally applicable to any system, its main advantages are the correct representation of commensurability issues and the treatment of degenerate systems. As an example we apply the method to a spin chain featuring a highly degenerate ground state space where the new sweeping scheme provides an increase in performance as well as accuracy by many orders of magnitude compared to a recently published work.


title: The Foucault Pendulum Precession and the Additivity of Infinitesimal Rotations

abstract: The purpose of this paper is to give an intuitive explanation of the Foucault pendulum precession (Fpp) by exploiting the easily proved result that infinitesimal spatial rotations about different axes in three dimensions are additive. This allows to explain in a simpler way the Foucault precession by combining physical laws, intuition and a few elementary notions of astronomy. The paper is by no means intended to replace the material given in standard texts of classical mechanics; it is only an interdisciplinary description of the Foucault precession tailored for non specialists.


title: Anomalous Exponents in Strong Turbulence

abstract: To characterize fluctuations in a turbulent flow, one usually studies different moments of velocity increments and dissipation rate, $\overline{(v(x+r)-v(x))^{n}}\propto r^{\zeta_{n}}$ and $\overline{{\cal E}^{n}}\propto Re^{d_{n}}$, respectively. In high Reynolds number flows, the moments of different orders cannot be simply related to each other which is the signature of anomalous scaling, one of the most puzzling features of turbulent flows. High-order moments are related to extreme, rare events and our ability to quantitatively describe them is crucially important for meteorology, heat, mass transfer and other applications. In this work we present a solution to this problem in the particular case of the Navier-Stokes equations driven by a random force. A novel aspect of this work is that, unlike previous efforts which aimed at seeking solutions around the infinite Reynolds number limit, we concentrate on the vicinity of transitional Reynolds numbers $Re^{tr}$ where the first emergence of anomalous scaling is observed out of a low-$Re$ Gaussian background. The obtained closed expressions for anomalous scaling exponents $\zeta_{n}$ and $d_{n}$, which depend on the transition Reynolds number, agree well with experimental and numerical data in the literature and, when $n\gg 1$, $d_{n}\approx 0.19n \ln(n)$. The theory yields the energy spectrum $E(k)\propto k^{-\zeta_{2}-1}$ with $\zeta_{2}\approx 0.699$, different from the outcome of Kolmogorov's theory. It is also argued that fluctuations of dissipation rate and those of the transition point itself are responsible for both, deviation from Gaussian statistics and multiscaling of velocity field.


title: Cut Finite Elements for Convection in Fractured Domains

abstract: We develop a cut finite element method (CutFEM) for the convection problem in a so called fractured domain which is a union of manifolds of different dimensions such that a $d$ dimensional component always resides on the boundary of a $d+1$ dimensional component. This type of domain can for instance be used to model porous media with embedded fractures that may intersect. The convection problem can be formulated in a compact form suitable for analysis using natural abstract directional derivative and divergence operators. The cut finite element method is based on using a fixed background mesh that covers the domain and the manifolds are allowed to cut through a fixed background mesh in an arbitrary way. We consider a simple method based on continuous piecewise linear elements together with weak enforcement of the coupling conditions and stabilization. We prove a priori error estimates and present illustrating numerical examples.


title: Invariants of multidimensional time series based on their iterated-integral signature

abstract: We introduce a novel class of features for multidimensional time series, that are invariant with respect to transformations of the ambient space. The general linear group, the group of rotations and the group of permutations of the axes are considered. The starting point for their construction is Chen's iterated-integral signature.


title: Overcoming the vanishing gradient problem in plain recurrent networks

abstract: Plain recurrent networks greatly suffer from the vanishing gradient problem while Gated Neural Networks (GNNs) such as Long-short Term Memory (LSTM) and Gated Recurrent Unit (GRU) deliver promising results in many sequence learning tasks through sophisticated network designs. This paper shows how we can address this problem in a plain recurrent network by analyzing the gating mechanisms in GNNs. We propose a novel network called the Recurrent Identity Network (RIN) which allows a plain recurrent network to overcome the vanishing gradient problem while training very deep models without the use of gates. We compare this model with IRNNs and LSTMs on multiple sequence modeling benchmarks. The RINs demonstrate competitive performance and converge faster in all tasks. Notably, small RIN models produce 12%--67% higher accuracy on the Sequential and Permuted MNIST datasets and reach state-of-the-art performance on the bAbI question answering dataset.


title: The atmospheric parameters of FGK stars using wavelet analysis of CORALIE spectra

abstract: Atmospheric properties of F-,G- and K-type stars can be measured by spectral model fitting. These methods require data with good signal-to-noise ratio and reliable continuum normalisation. This is particularly challenging for the spectra we have obtained with the CORALIE \'{e}chelle spectrograph for FGK stars with transiting M-dwarf companions. The spectra tend to have low signal-to-noise ratios, which makes it difficult to analyse them using existing methods. Our aim is to create a reliable automated spectral analysis routine to determine $T_{\rm eff}$, [Fe/H], $V \sin i$ from the CORALIE spectra of FGK stars. We use wavelet decomposition to distinguish between noise, continuum trends, and stellar spectral features in the CORALIE spectra. A subset of wavelet coefficients from the target spectrum are compared to those from a grid of models in a Bayesian framework to determine the posterior probability distributions of the atmospheric parameters. By testing our method using synthetic spectra we found that our method converges on the best fitting atmospheric parameters. We test the wavelet method on 20 FGK exoplanet host stars for which higher quality data have been independently analysed using equivalent width measurements. We find that we can determine $T_{\rm eff}$ to a precision of $85$K, [Fe/H] to a precision of 0.06 dex and $V \sin i$ to a precision of 1.35 kms$^{-1}$ for stars with $V \sin i$ $\geq$ 5 kms$^{-1}$. We find an offset in metallicity $\approx -$0.18 dex relative to the equivelent width fitting method. We can determine $\log g$ to a precision of $0.13$ dex but measurements are only reliable to confirm dwarf-like surface gravity. The wavelet method can be used to determine $T_{\rm eff}$, [Fe/H] and $V \sin i$ for FGK stars from \'{e}chelle spectra. We find that our method is self consistent, and robust for spectra with SNR$\gtrapprox 40$.


title: Challenges of the Dynamic Detection of Functionally Similar Code Fragments

abstract: Classic clone detection approaches are hardly capable of finding redundant code that has been developed independently, i.e., is not the result of copy&paste. To automatically detect such functionally similar code of independent origin, we experimented with a dynamic detection approach that applies random testing to selected chunks of code similar to Jiang&Su's approach. We found that such an approach faces several limitations in its application to diverse Java systems. This paper details on our insights regarding these challenges of dynamic detection of functionally similar code fragments. Our findings support a substantiated discussion on detection approaches and serve as a starting point for future research.


title: The universality of electronic friction II: Equivalence of the quantum-classical Liouville equation approach with von Oppen's nonequilibrium Green's function methods out of equilibrium

abstract: In a recent publication [W. Dou, G. Miao, and J. E. Subotnik, Phys. Rev. Lett. 119, 046001 (2017)], using the quantum-classical Liouville equation (QCLE), we derived a very general form for the electronic friction felt by a molecule moving near one or many metal surfaces. Moreover, we have already proved the equivalence of the QCLE electronic friction with the Head-Gordon--Tully model as well as a generalized version of von Oppen's nonequilibrium Green's function (NEGF) method at equilibrium [W. Dou and J. E. Subotnik, Phys. Rev. B 96, 104305 (2017)]. In the present paper, we now further prove the equivalence between the QCLE friction and the NEGF friction for the case of multiple metal surfaces and an out-of-equilibrium electronic current. The present results conclude our recent claim that there is only one universal electronic friction tensor arising from the Born-Oppenheimer approximation.


title: Material Discovery and Design Principles for Stable, High Activity Perovskite Cathodes for Solid Oxide Fuel Cells

abstract: Critical to the development of improved solid oxide fuel cell (SOFC) technology are novel compounds with high oxygen reduction reaction (ORR) catalytic activity and robust stability under cathode operating conditions. Approximately 2145 distinct perovskite compositions are screened for potential use as high activity, stable SOFC cathodes, and it is verified that the screening methodology qualitatively reproduces the experimental activity, stability, and conduction properties of well-studied cathode materials. The calculated oxygen p-band center is used as a first principle-based descriptor of the surface exchange coefficient (k*), which in turn correlates with cathode ORR activity. Convex hull analysis is used under operating conditions in the presence of oxygen, hydrogen, and water vapor to determine thermodynamic stability. This search has yielded 52 potential cathode materials with good predicted stability in typical SOFC operating conditions and predicted k* on par with leading ORR perovskite catalysts. The established trends in predicted k* and stability are used to suggest methods of improving the performance of known promising compounds. The material design strategies and new materials discovered in the computational search help enable the development of high activity, stable compounds for use in future solid oxide fuel cells and related applications.


title: Least primitive root and simultaneous power-non residues

abstract: Let $p$ be a prime and let $g(p)$ be the least primitive root modulo $p$. We prove that for any $\epsilon>0$ and $p$ large enough the bound \begin{align} g(p)\ll p^{\frac{1}{4\sqrt{e}}+\epsilon} \nonumber \end{align} holds for most prime $p$ such that $p-1$ does not have small prime factors, but $2$. We also give an explicit description of the exceptional set.


title: SUPERSHARP - Segmented Unfolding Primary for Exoplanet Research via Spectroscopic High Angular Resolution Photography

abstract: We propose to search for biosignatures in the spectra of reflected light from about 100 Earth-sized planets that are already known to be orbiting in their habitable zones (HZ). For a sample of G and K type hosts, most of these planets will be between 25 and 50 milli-arcsec (mas) from their host star and 1 billion to 10 billion times fainter. To separate the planet's image from that of its host star at the wavelength (763nm) of the oxygen biosignature we need a telescope with an aperture of 16 metres. Furthermore, the intensity of the light from the host star at the position in the image of the exoplanet must be suppressed otherwise the exoplanet will be lost in the glare. This presents huge technical challenges. The Earth's atmosphere is turbulent which makes it impossible to achieve the required contrast from the ground at 763nm. The telescope therefore needs to be in space and to fit the telescope in the rocket fairing it must be a factor of 4 or more times smaller when folded than when operational. To obtain spectroscopy of the planet's biosignature at 763nm we need to use an integral field spectrometer (IFS) with a field of view (FOV) of 1000 x 1000 milli-arcsec (mas) and a spectral resolution of 100. This is a device that simultaneously takes many pictures of the exoplanet each at a slightly different wavelength which are then recorded as a data cube with two spatial dimensions and one wavelength dimension. In every data cube wavelength slice, the background light from the host star at the location of the planet image must be minimised. This is achieved via a coronagraph which blocks the light from the host star and active/adaptive optics techniques which continuously maintain very high accuracy optical alignment to make the images as sharp as possible. These are the technical challenges to be addressed in a design study.


title: Ideals modulo a prime

abstract: The main focus of this paper is on the problem of relating an ideal $I$ in the polynomial ring $\mathbb Q[x_1, \dots, x_n]$ to a corresponding ideal in $\mathbb F_p[x_1,\dots, x_n]$ where $p$ is a prime number; in other words, the \textit{reduction modulo $p$} of $I$. We first define a new notion of $\sigma$-good prime for $I$ which does depends on the term ordering $\sigma$, but not on the given generators of $I$. We relate our notion of $\sigma$-good primes to some other similar notions already in the literature. Then we introduce and describe a new invariant called the universal denominator which frees our definition of reduction modulo~$p$ from the term ordering, thus letting us show that all but finitely many primes are good for $I$. One characteristic of our approach is that it enables us to easily detect some bad primes, a distinct advantage when using modular methods.


title: Colour Rearrangement for Dipole Showers

abstract: We present an algorithm to rearrange the colour chains of dipole showers in the shower process according to the colour amplitudes of a simple matrix element. We implement the procedure in the dipole shower of Herwig and show comparisons to data.


title: Heavy Anionic Complex Creates a Unique Water Structure at a Soft Charged Interface

abstract: Ion hydration and interfacial water play crucial roles in numerous phenomena ranging from biological to industrial systems. Although biologically relevant (and mostly smaller) ions have been studied extensively in this context, very little experimental data exist about molecular scale behavior of heavy ions and their complexes at interfaces, especially under technologically significant conditions. It has recently been shown that PtCl62- complexes adsorb at positively charged interfaces in a two-step process that cannot fit into well-known empirical trends, such as Hofmeister series. Here, a combined vibrational sum frequency generation and molecular dynamics study reveals that a unique interfacial water structure is connected to this peculiar adsorption behavior. A novel sub-ensemble analysis of MD simulation results show that after adsorption, PtCl62- complexes partially retain their first and second hydration spheres, and it is possible to identify three different types of water molecules around them based on their orientational structures and hydrogen bonding strengths. These results have important implications for relating interfacial water structure and hydration enthalpy to the general understanding of specific ion effects. This in turn influences interpretation of heavy metal ion distribution across and reactivity within, liquid interfaces.


title: Fidelity deviation in quantum teleportation

abstract: We analyze the performance of quantum teleportation in terms of average fidelity and fidelity deviation. The average fidelity is defined as the average value of the fidelities over all possible input states and the fidelity deviation is their standard deviation, which is referred to as a concept of fluctuation or universality. In the analysis, we find the condition to optimize both measures under a noisy quantum channel---we here consider the so-called Werner channel. To characterize our results, we introduce a two-dimensional space defined by the aforementioned measures, in which the performance of the teleportation is represented as a point with the channel noise parameter. Through further analysis, we specify some regions drawn for different channel conditions, establishing the connection to the dissimilar contributions of the entanglement to the teleportation and the Bell inequality violation.


title: Long-Term Stability of Planets in the $\alpha$ Centauri System, II: Forced Eccentricities

abstract: We extend our study of the extent of the regions within the $\alpha$ Centauri AB star system where small planets are able to orbit for billion-year timescales (Quarles & Lissauer 2016, AJ 151, 111) to investigate the effects of minimizing the forced eccentricity of initial trajectories. We find that initially prograde, circumstellar orbits require a piecewise quadratic function to accurately approximate forced eccentricity as a function of semimajor axis, but retrograde orbits can be modeled using a linear function. Circumbinary orbits in the $\alpha$ Centauri AB system are less affected by the forced eccentricity. {Planets on circumstellar orbits that begin with eccentricity vectors near their forced values are generally stable, up to $\sim$$10^9$ yr, out to a larger semimajor axis than are planets beginning on circular orbits. The amount by which the region of stability expands is much larger for retrograde orbits than it is for prograde orbits. The location of the stability boundary for two planet systems on prograde, circular orbits is much more sensitive to the initial eccentricity state than it is for analogous single planet systems.


title: Origins of Hot Jupiters

abstract: Hot Jupiters were the first exoplanets to be discovered around main sequence stars and astonished us with their close-in orbits. They are a prime example of how exoplanets have challenged our textbook, solar-system inspired story of how planetary systems form and evolve. More than twenty years after the discovery of the first hot Jupiter, there is no consensus on their predominant origin channel. Three classes of hot Jupiter creation hypotheses have been proposed: in situ formation, disk migration, and high eccentricity tidal migration. Although no origin channel alone satisfactorily explains all the evidence, two major origins channels together plausibly account for properties of hot Jupiters themselves and their connections to other exoplanet populations.


title: Dynamical AdS/Yang-Mills Model

abstract: The AdS/CFT correspondence has provided new and useful insights into the nonperturbative regime of strongly coupled gauge theories. We construct a class of models meant to mimic Yang-Mills theory using the superpotential method. This method allows us to efficiently address the problem of solving all the equations of motion. The conformal symmetry is broken in the infrared by a dilaton field. Using a five-dimensional action we calculate the mass spectrum of scalar glueballs. This spectrum contains a tachyon, indicating an instability in the theory. We stabilize the theory by introducing a cosmological constant in the bulk and a pair of 3-branes, as in the Randall-Sundrum model. The scalar glueball masses computed by lattice gauge theory are then described very well by just a few parameters. Prospects for extending the model to other spins and parities and to finite temperature are considered.


title: Retardation effects in spectroscopic measurements of the Casimir-Polder interaction

abstract: Spectroscopy is a unique experimental tool for measuring the fundamental Casimir-Polder interaction between excited state atoms, or other polarizable quantum objects, and a macroscopic surface. Spectroscopic measurements probe atoms at nanometric distances away from the surface where QED retardation is usually negligible and the atom-surface interaction is proportional to the inverse cube of the separation distance, otherwise known as the van der Waals regime. Here we focus on selective reflection, one of the main spectroscopic probes of Casimir-Polder interactions. We calculate for the first time selective reflection spectra using the full, distance dependent, Casimir-Polder energy shift and linewidth. We demonstrate that retardation can have significant effects, in particular for experiments with low-lying energy states. We also show that the effective probing depth of selective reflection spectroscopy depends on the transition linewidth. Our analysis allows us to calculate selective reflection spectra with composite surfaces, such as metasurfaces, dielectric stacks, or even bi-dimensional materials.


title: Fluid and gyrofluid modeling of low-$\beta_e$ plasmas: phenomenology of kinetic Alfv\'en wave turbulence

abstract: Reduced fluid models including electron inertia and ion finite Larmor radius corrections are derived asymptotically, both from fluid basic equations and from a gyrofluid model. They apply to collisionless plasmas with small ion-to-electron equilibrium temperature ratio and low $\beta_e$, where $\beta_e$ indicates the ratio between the equilibrium electron pressure and the magnetic pressure exerted by a strong, constant and uniform magnetic guide field. The consistency between the fluid and gyrofluid approaches is ensured when choosing ion closure relations prescribed by the underlying ordering. A two-field reduction of the gyrofluid model valid for arbitrary equilibrium temperature ratio is also introduced, and is shown to have a noncanonical Hamiltonian structure. This model provides a convenient framework for studying kinetic Alfv\'en wave turbulence, from MHD to sub-$d_e$ scales (where $d_e$ holds for the electron skin depth). Magnetic energy spectra are phenomenologically determined within energy and generalized helicity cascades in the perpendicular spectral plane. Arguments based on absolute statistical equilibria are used to predict the direction of the transfers, pointing out that, within the sub-ion range associated with a $k_\perp^{-7/3}$ transverse magnetic spectrum, the generalized helicity could display an inverse cascade if injected at small scales, for example by reconnection processes.


title: Real Randomized Benchmarking

abstract: Randomized benchmarking provides a tool for obtaining precise quantitative estimates of the average error rate of a physical quantum channel. Here we define real randomized benchmarking, which enables a separate determination of the average error rate in the real and complex parts of the channel. This provides more fine-grained information about average error rates with approximately the same cost as the standard protocol. The protocol requires only averaging over the real Clifford group, a subgroup of the full complex Clifford group, and makes use of the fact that it forms an orthogonal 2-design. It therefore allows benchmarking of fault-tolerant gates for an encoding which does not contain the full Clifford group transversally. Furthermore, our results are especially useful when considering quantum computations on rebits (or real encodings of complex computations), in which case the real Clifford group now plays the role of the complex Clifford group when studying stabilizer circuits.


title: Anatomy of an online misinformation network

abstract: Massive amounts of fake news and conspiratorial content have spread over social media before and after the 2016 US Presidential Elections despite intense fact-checking efforts. How do the spread of misinformation and fact-checking compete? What are the structural and dynamic characteristics of the core of the misinformation diffusion network, and who are its main purveyors? How to reduce the overall amount of misinformation? To explore these questions we built Hoaxy, an open platform that enables large-scale, systematic studies of how misinformation and fact-checking spread and compete on Twitter. Hoaxy filters public tweets that include links to unverified claims or fact-checking articles. We perform k-core decomposition on a diffusion network obtained from two million retweets produced by several hundred thousand accounts over the six months before the election. As we move from the periphery to the core of the network, fact-checking nearly disappears, while social bots proliferate. The number of users in the main core reaches equilibrium around the time of the election, with limited churn and increasingly dense connections. We conclude by quantifying how effectively the network can be disrupted by penalizing the most central nodes. These findings provide a first look at the anatomy of a massive online misinformation diffusion network.


title: Luminous efficiency estimates of meteors -II. Application to Canadian Automated Meteor Observatory meteor events

abstract: Luminous efficiency is a necessary parameter for determining meteoroid mass from optical emission. Despite this importance, it is very poorly known, with previous results varying by up to two orders of magnitude for a given speed. We present the most recent study of luminous efficiency values determined with modern high-resolution instruments, by directly comparing dynamic and photometric meteoroid masses. Fifteen non-fragmenting meteoroids were used, with a further five clearly fragmenting events for comparison. Twelve of the fifteen non-fragmenting meteoroids had luminous efficiencies less than 1%, while the fragmenting meteoroids had upper limits of a few tens of percent. No clear trend with speed was seen, but there was a weak negative trend of luminous efficiency on meteoroid mass, implying that smaller meteoroids radiate more efficiently.


title: Prompt gamma-ray emission of GRB 170817A associated to GW 170817: A consistent picture

abstract: The short GRB 170817A associated to the first detection of gravitation waves from a Binary Neutron Star (BNS) merger was in many ways unusual. Possible explanations are emission from a cocoon or cocoon break out, off-axis view of a structured or uniform jet, and on-axis ultra-relativistic jet with reduced density and Lorentz factor. Here we use a phenomenological model of shock evolution and synchrotron/self-Compton emission to simulate the prompt emission of GRB 170817A and to test above proposals. We find that synchrotron emission from a mildly relativistic cocoon with a Lorentz factor of 2-3, as considered in the literature, generates a too soft, too long, and too bright prompt emission. Off-axis view of an structured jet with a Lorentz factor of about 10 can reproduce observations, but needs a very efficient transfer of kinetic energy to electrons in internal shocks, which is disfavored by particle in cell simulations. We also comment on cocoon breakout as a mechanism for generation of the prompt gamma-ray. A relativistic jet with a Lorentz factor of about 100 and a density lower than typical short GRBs seems to be the most plausible model and we conclude that GRB 170817A was intrinsically faint. Based on this result and findings of relativistic magnetohydrodynamics simulations of BNS merger in the literature we discuss physical and astronomical conditions, which may lead to such faint short GRBs. We identify small mass difference of progenitor neutron stars, their old age and reduced magnetic field, and anti-alignment of spin-orbit angular momentum induced by environmental gravitational disturbances during the lifetime of the BNS as causes for the faintness of GRB 170817A. We predict that BNS mergers at lower redshifts generate on average fainter GRBs.


title: Classifying braidings on fusion categories

abstract: We show that braidings on a fusion category $\mathcal{C}$ correspond to certain fusion subcategories of the center of $\mathcal{C}$ transversal to the canonical Lagrangian algebra. This allows to classify braidings on non-degenerate and group-theoretical fusion categories.


title: Non-Adversarial Unsupervised Word Translation

abstract: Unsupervised word translation from non-parallel inter-lingual corpora has attracted much research interest. Very recently, neural network methods trained with adversarial loss functions achieved high accuracy on this task. Despite the impressive success of the recent techniques, they suffer from the typical drawbacks of generative adversarial models: sensitivity to hyper-parameters, long training time and lack of interpretability. In this paper, we make the observation that two sufficiently similar distributions can be aligned correctly with iterative matching methods. We present a novel method that first aligns the second moment of the word distributions of the two languages and then iteratively refines the alignment. Extensive experiments on word translation of European and Non-European languages show that our method achieves better performance than recent state-of-the-art deep adversarial approaches and is competitive with the supervised baseline. It is also efficient, easy to parallelize on CPU and interpretable.


title: Reduced order models for fluid-structure interaction problems with applications in haemodynamics

abstract: This paper deals with fast simulations of the haemodynamics in large arteries by considering a reduced model of the associated fluid-structure interaction problem, which in turn allows an additional reduction in terms of the numerical discretisation. The resulting method is both accurate and computationally cheap. This goal is achieved by means of two levels of reduction: first, we describe the model equations with a reduced mathematical formulation which allows to write the fluid-structure interaction problem as a Navier-Stokes system with non-standard boundary conditions; second, we employ numerical reduction techniques to further and drastically lower the computational costs. The numerical reduction is obtained coupling two well-known techniques: the proper orthogonal decomposition and the reduced basis method, in particular the greedy algorithm. We start by reducing the numerical dimension of the problem at hand with a proper orthogonal decomposition and we measure the system energy with specific norms; this allows to take into account the different orders of magnitude of the state variables, the velocity and the pressure. Then, we introduce a strategy based on a greedy procedure which aims at enriching the reduced discretization space with low offline computational costs. As application, we consider a realistic haemodynamics problem with a perturbation in the boundary conditions and we show the good performances of the reduction techniques presented in the paper. The gains obtained in term of CPU time are of three orders of magnitude.


title: A methodology for calculating the latency of GPS-probe data

abstract: Crowdsourced GPS probe data has been gaining popularity in recent years as a source for real-time traffic information. Efforts have been made to evaluate the quality of such data from different perspectives. A quality indicator of any traffic data source is latency that describes the punctuality of data, which is critical for real-time operations, emergency response, and traveler information systems. This paper offers a methodology for measuring the probe data latency, with respect to a selected reference source. Although Bluetooth re-identification data is used as the reference source, the methodology can be applied to any other ground-truth data source of choice (i.e. Automatic License Plate Readers, Electronic Toll Tag). The core of the methodology is a maximum pattern matching algorithm that works with three different fitness objectives. To test the methodology, sample field reference data were collected on multiple freeways segments for a two-week period using portable Bluetooth sensors as ground-truth. Equivalent GPS probe data was obtained from a private vendor, and its latency was evaluated. Latency at different times of the day, the impact of road segmentation scheme on latency, and sensitivity of the latency to both speed slowdown, and recovery from slowdown episodes are also discussed.


title: Right sign of spin rotation operator

abstract: For the fermion transformation in the space all books of quantum mechanics propose to use the unitary operator $\widehat{U}_{\vec n}(\varphi)=\exp{(-i\frac\varphi2(\widehat\sigma\cdot\vec n))}$, where $\varphi$ is angle of rotation around the axis $\vec{n}$. But this operator turns the spin in inverse direction presenting the rotation to the left. The error of defining of $\widehat{U}_{\vec n}(\varphi)$ action is caused because the spin supposed as simple vector which is independent from $\widehat\sigma$-operator a priori. In this work it is shown that each fermion marked by number $i$ has own Pauli-vector $\widehat\sigma_i$ and both of them change together. If we suppose the global $\widehat\sigma$-operator and using the Bloch Sphere approach define for all fermions the common quantization axis $z$ the spin transformation will be the same: the right hand rotation around the axis $\vec{n}$ is performed by the operator $\widehat{U}^+_{\vec n}(\varphi)=\exp{(+i\frac\varphi2(\widehat\sigma\cdot\vec n))}$.


title: Genus zero Gopakumar-Vafa type invariants for Calabi-Yau 4-folds II: Fano 3-folds

abstract: In analogy with the Gopakumar-Vafa (GV) conjecture on Calabi-Yau (CY) 3-folds, Klemm and Pandharipande defined GV type invariants on Calabi-Yau 4-folds using Gromov-Witten theory and conjectured their integrality. In a joint work with Maulik and Toda, the author conjectured their genus zero invariants are $\mathrm{DT_4}$ invariants of one dimensional stable sheaves. In this paper, we study this conjecture on the total space of canonical bundle of a Fano 3-fold $Y$, which reduces to a relation between twisted GW and $\mathrm{DT_3}$ invariants on $Y$. Examples are computed for both compact and non-compact Fano 3-folds to support our conjecture.


title: Long-Term Stability of Tightly Packed Multi-Planet Systems in Prograde, Coplanar, Circumstellar Orbits within the $\alpha$ Centauri AB System

abstract: We perform long-term simulations, up to ten billion years, of closely-spaced configurations of 2 -- 6 planets, each as massive as the Earth, traveling on nested orbits about either stellar component in $\alpha$ Centauri AB. The innermost planet initially orbits at either the inner edge of its star's empirical habitable zone (HZ) or the inner edge of its star's conservative HZ. Although individual planets on low inclination, low eccentricity, orbits can survive throughout the habitable zones of both stars, perturbations from the companion star require that the minimum spacing of planets in multi-planet systems within the habitable zones of each star must be significantly larger than the spacing of similar multi-planet systems orbiting single stars in order to be long-lived. The binary companion induces a forced eccentricity upon the orbits of planets in orbit around either star. Planets on appropriately-phased circumstellar orbits with initial eccentricities equal to their forced eccentricities can survive on more closely spaced orbits than those with initially circular orbits, although the required spacing remains higher than for planets orbiting single stars. A total of up to nine planets on nested prograde orbits can survive for the current age of the system within the empirical HZs of the two stars, with five of these orbiting $\alpha$ Centauri B and four orbiting $\alpha$ Centauri A.


title: Scaling properties of multiscale equilibration

abstract: We investigate the lattice spacing dependence of the equilibration time for a recently proposed multiscale thermalization algorithm for Markov chain Monte Carlo simulations. The algorithm uses a renormalization-group matched coarse lattice action and prolongation operation to rapidly thermalize decorrelated initial configurations for evolution using a corresponding target lattice action defined at a finer scale. Focusing on non-topological long-distance observables in pure SU(3) gauge theory, we provide quantitative evidence that the slow modes of the Markov process, which provide the dominant contribution to the rethermalization time, have a suppressed contribution toward the continuum limit, despite their associated timescales increasing. Based on these numerical investigations, we conjecture that the prolongation operation used herein will produce ensembles that are indistinguishable from the target fine-action distribution for a sufficiently fine coupling at a given level of statistical precision, thereby eliminating the cost of rethermalization.


title: Regulation of Pom cluster dynamics in Myxococcus xanthus

abstract: Precise positioning of the cell division site is essential for the correct segregation of the genetic material into the two daughter cells. In the bacterium Myxococcus xanthus, the proteins PomX and PomY form a cluster on the chromosome that performs a biased random walk to midcell and positively regulates cell division there. PomZ, an ATPase, is necessary for tethering of the cluster to the nucleoid and regulates its movement towards midcell. It has remained unclear how the cluster dynamics change when the biochemical parameters, such as the attachment rates of PomZ to the nucleoid and the cluster, the ATP hydrolysis rate of PomZ or the mobility of PomZ dimers interacting with the nucleoid and cluster, are varied. To answer these questions, we investigate a one-dimensional model that includes the nucleoid, the Pom cluster and the PomZ protein. We find that a mechanism based on the diffusive PomZ fluxes on the nucleoid into the cluster can explain the latter's midnucleoid localization for a broad parameter range. Furthermore, there is an ATP hydrolysis rate that minimizes the time the cluster needs to reach midnucleoid. If the dynamics of PomZ dimers on the nucleoid is slow relative to the cluster's velocity, we observe oscillatory cluster movements around midnucleoid. To understand midnucleoid localization, we developed a semi-analytical approach that dissects the net movement of the cluster into its components: the difference in PomZ fluxes into the cluster from either side, the force exerted by a single PomZ dimer on the cluster and the effective friction coefficient of the cluster. Importantly, we predict that the Pom cluster oscillates around midnucleoid if the diffusivity of PomZ on the nucleoid is reduced. A similar approach to that applied here may also prove useful for cargo localization in ParABS systems.


title: An L+T Spectral Binary with Possible AB Doradus Kinematics

abstract: We present the identification of WISE~J135501.90$-$825838.9 as a spectral binary system with a slight possibility of planetary-mass components in the $130-200$\,Myr AB~Doradus moving group. Peculiarities in the near-infrared spectrum of this source suggest it to be a blended-light binary with L6.0$\pm$1.0 and T3.0$\pm$1.8 or L7.0$\pm$0.6 and T7.5$\pm$0.4 components. Its proper motion and radial velocity as a combined-light source yield a high membership probability for AB~Doradus. While the young L6+T3 case is underluminous in a color-magnitude diagram at the AB~Doradus kinematic distance, the young L7+T7.5 case could be viable. Gravity-sensitive indicators are more consistent with a field-age binary. If confirmed as a young object, member of AB~Doradus, we estimate masses of 11$\pm$1\,\mjup\ and 9$\pm$1\,\mjup\ with both component masses below the Deuterium burning mass limit. Otherwise, we find masses of $72^{+4}_{-5}$ and $61^{+6}_{-8}$ for the field L6+T3 case and $70^{+2}_{-4}$ and $42^{+5}_{-6}$ for the field L7+T7.5 case. Our identification of WISE~J135501.90$-$825838.9 as a candidate young spectral binary introduces a new technique for detecting and characterizing planetary-mass companions to young brown dwarfs.


title: Robust Foregrounds Removal for 21-cm Experiments

abstract: Direct detection of the Epoch of Reionization via the redshifted 21-cm line will have unprecedented implications on the study of structure formation in the early Universe. To fulfill this promise current and future 21-cm experiments will need to detect the weak 21-cm signal over foregrounds several order of magnitude greater. This requires accurate modeling of the galactic and extragalactic emission and of its contaminants due to instrument chromaticity, ionosphere and imperfect calibration. To solve for this complex modeling, we propose a new method based on Gaussian Process Regression (GPR) which is able to cleanly separate the cosmological signal from most of the foregrounds contaminants. We also propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere. Using this method, chromatic effects causing the so-called "wedge" are effectively eliminated (i.e. deconvolved) in the cylindrical ($k_{\perp}, k_{\parallel}$) power spectrum.


title: Latitude: A Model for Mixed Linear-Tropical Matrix Factorization

abstract: Nonnegative matrix factorization (NMF) is one of the most frequently-used matrix factorization models in data analysis. A significant reason to the popularity of NMF is its interpretability and the `parts of whole' interpretation of its components. Recently, max-times, or subtropical, matrix factorization (SMF) has been introduced as an alternative model with equally interpretable `winner takes it all' interpretation. In this paper we propose a new mixed linear--tropical model, and a new algorithm, called Latitude, that combines NMF and SMF, being able to smoothly alternate between the two. In our model, the data is modeled using the latent factors and latent parameters that control whether the factors are interpreted as NMF or SMF features, or their mixtures. We present an algorithm for our novel matrix factorization. Our experiments show that our algorithm improves over both baselines, and can yield interpretable results that reveal more of the latent structure than either NMF or SMF alone.


title: Ehrenfest statistical dynamics in chemistry: study of decoherence effects

abstract: In previous works, we introduced a geometric route to define our Ehrenfest Statistical Dynamics (ESD) and we proved that, for a simple toy-model, the resulting ESD does not preserve purity. We now take a step further: we investigate decoherence and pointer basis in the ESD model by considering some uncertainty in the degrees of freedom of a simple but realistic molecular model, consisting of two classical cores and one quantum electron. The Ehrenfest model is sometimes discarded as a valid approximation to non-adiabatic coupled quantum-classical dynamics because it does not describe the decoherence in the quantum subsystem. However, any rigorous statistical analysis of the Ehrenfest dynamics, such as the described ESD formalism, proves that decoherence exists. In this article, decoherence in ESD is studied by measuring the change in the quantum subsystem purity and by analysing the appearance of the pointer basis to which the system decoheres, which for our example is composed by the eigenstates of the electronic Hamiltonian.


title: Primordial black holes with an accurate QCD equation of state

abstract: Making use of definitive new lattice computations of the Standard Model thermodynamics during the quantum chromodynamic (QCD) phase transition, we calculate the enhancement in the mass distribution of primordial black holes (PBHs) due to the softening of the equation of state. We find that the enhancement peaks at approximately $0.7M_\odot$, with the formation rate increasing by at least two orders of magnitude due to the softening of the equation of state at this time, with a range of approximately $0.3M_\odot<M<1.4M_\odot$ at full width half-maximum. PBH formation is increased by a smaller amount for PBHs with masses spanning a large range, $10^{-3}M_\odot<M_{\rm PBH}<10^{3}M_\odot$, which includes the masses of the BHs that LIGO detected. The most significant source of uncertainty in the number of PBHs formed is now due to unknowns in the formation process, rather than from the phase transition. A near scale-invariant density power spectrum tuned to generate a population with mass and merger rate consistent with that detected by LIGO should also produce a much larger energy density of PBHs with solar mass. The existence of BHs below the Chandresekhar mass limit would be a smoking gun for a primordial origin and they could arguably constitute a significant fraction of the cold dark matter density. They also pose a challenge to inflationary model building which seek to produce the LIGO BHs without overproducing lighter PBHs.


title: Towards uncovering the structure of power fluctuations of wind farms

abstract: The structure of the turbulence-driven power fluctuations in a wind farm is fundamentally described from basic concepts. A derived tuning-free model, supported with experiments, reveals the underlying spectral content of the power fluctuations of a wind farm. It contains two power-law trends and oscillations in the relatively low- and high-frequency ranges. The former is mostly due to the turbulent interaction between the flow and the turbine properties; whereas the latter is due to the advection between turbine pairs. The spectral wind-farm scale power fluctuations $\Phi_P$ exhibits a power-law decay proportional to $f^{-5/3-2}$ in the region corresponding to the turbulence inertial subrange and at relatively large scales, $\Phi_P\sim f^{-2}$. Due to the advection and turbulent diffusion of large-scale structures, a spectral oscillation exists with the product of a sinusoidal behavior and an exponential decay in the frequency domain.


title: Enhanced polymer capture speed and extended translocation time in pressure-solvation traps

abstract: The efficiency of nanopore-based biosequencing techniques requires fast anionic polymer capture by like-charged pores followed by a prolonged translocation process. We show that this condition can be achieved by setting a pressure-solvation trap. Polyvalent cation addition to the KCl solution triggers the like-charge polymer-pore attraction. The attraction speeds-up the pressure-driven polymer capture but also traps the molecule at the pore exit, reducing the polymer capture time and extending the polymer escape time by several orders of magnitude. By direct comparison with translocation experiments [D. P. Hoogerheide et al., ACS Nano 8, 7384 (2014)], we characterize as well the electrohydrodynamics of polymers transport in pressure-voltage traps. We derive scaling laws that can accurately reproduce the pressure dependence of the experimentally measured polymer translocation velocity and time. We also find that during polymer capture, the electrostatic barrier on the translocating molecule slows down the liquid flow. This prediction identifies the streaming current measurement as a potential way to probe electrostatic polymer-pore interactions.


title: A First Option Calibration of the GARCH Diffusion Model by a PDE Method

abstract: Time-series calibrations often suggest that the GARCH diffusion model could also be a suitable candidate for option (risk-neutral) calibration. But unlike the popular Heston model, it lacks a fast, semi-analytic solution for the pricing of vanilla options, perhaps the main reason why it is not used in this way. In this paper we show how an efficient finite difference-based PDE solver can effectively replace analytical solutions, enabling accurate option calibrations in less than a minute. The proposed pricing engine is shown to be robust under a wide range of model parameters and combines smoothly with black-box optimizers. We use this approach to produce a first PDE calibration of the GARCH diffusion model to SPX options and present some benchmark results for future reference.


title: Neutron moderation spectrum considering inelastic scattering

abstract: For the first time an analytic expression was obtained for the inelastic neutron scattering law with an isotropic neutron source within the gas model, considering moderating medium temperature as a parameter. The inelastic scattering law is obtained, based on the solution of the kinematic problem of neutron inelastic scattering on a nucleus in laboratory coordinate system (L-system) in general case. I.e. in case not only a neutron but also a nucleus have arbitrary velocity vector in L-system. Analytic expressions are found for the neutron flux density and moderation spectrum in reactor fissile medium, both in case of the elastic scattering law, obtained earlier by the authors, and in case of the inelastic scattering law obtained in this paper. Both elastic and inelastic scattering laws are considered to be dependent on the medium temperature. The obtained expressions for neutron moderation spectra enable reinterpretation of physical nature of the processes that determine the shape of neutron spectrum in a wide energy range.


title: TASEP fluctuations with soft-shock initial data

abstract: We consider the totally asymmetric simple exclusion process with \emph{soft-shock} initial particle density, which is a step function increasing in the direction of flow and the step size chosen small to admit KPZ scaling. The initial configuration is deterministic and the dynamics create a shock. We prove that the fluctuations of a particle at the macroscopic position of the shock converge to the maximum of two independent GOE Tracy-Widom random variables, which establishes a conjecture of Ferrari and Nejjar. Furthermore, we show the joint fluctuations of particles near the shock are determined by the maximum of two lines described in terms of these two random variables. The microscopic position of the shock is then seen to be their difference. Our proofs rely on determinantal formulae and a novel factorization of the associated kernels.


title: A Formalization for Specifying and Implementing Correct Pull-Stream Modules

abstract: Pull-stream is a JavaScript demand-driven functional design pattern based on callback functions that enables the creation and easy composition of independent modules that are used to create streaming applications. It is used in popular open source projects and the community around it has created over a hundred compatible modules. While the description of the pull-stream design pattern may seem simple, it does exhibit complicated termination cases. Despite the popularity and large uptake of the pull-stream design pattern, there was no existing formal specification that could help programmers reason about the correctness of their implementations. Thus, the main contribution of this paper is to provide a formalization for specifying and implementing correct pull-stream modules based on the following: (1) we show the pull-stream design pattern is a form of declarative concurrent programming; (2) we present an event-based protocol language that supports our formalization, independently of JavaScript; (3) we provide the first precise and explicit definition of the expected sequences of events that happen at the interface of two modules, which we call the pull-stream protocol; (4) we specify reference modules that exhibit the full range of behaviors of the pull-stream protocol; (5) we validate our definitions against the community expectations by testing the existing core pull-stream modules against them and identify unspecified behaviors in existing modules. Our approach helps to better understand the pull-stream protocol, to ensure interoperability of community modules, and to concisely and precisely specify new pull-stream abstractions in papers and documentation.


title: Improved perturbative QCD formalism for $B_c$ meson decays

abstract: We derive the $k_T$ resummation for doubly heavy-flavored $B_c$ meson decays by including the charm quark mass effect into the known formula for a heavy-light system. The resultant Sudakov factor is employed in the perutrbative QCD study of the "golden channel" $B_c^+ \to J/\psi \pi^+$. With a reasonable model for the $B_c$ meson distribution amplitude, which maintains approximate on-shell conditions of both the partonic bottom and charm quarks, it is observed that the imaginary piece of the $B_c \to J/\psi$ transition form factor appears to be power suppressed, and the $B_c^+ \to J/\psi \pi^+$ branching ratio is not lower than $10^{-3}$. The above improved perutrbative QCD formalism is applicable to $B_c$ meson decays to other charmonia and charmed mesons.


title: Universal Language Model Fine-tuning for Text Classification

abstract: Inductive transfer learning has greatly impacted computer vision, but existing approaches in NLP still require task-specific modifications and training from scratch. We propose Universal Language Model Fine-tuning (ULMFiT), an effective transfer learning method that can be applied to any task in NLP, and introduce techniques that are key for fine-tuning a language model. Our method significantly outperforms the state-of-the-art on six text classification tasks, reducing the error by 18-24% on the majority of datasets. Furthermore, with only 100 labeled examples, it matches the performance of training from scratch on 100x more data. We open-source our pretrained models and code.


title: Upgrading from Gaussian Processes to Student's-T Processes

abstract: Gaussian process priors are commonly used in aerospace design for performing Bayesian optimization. Nonetheless, Gaussian processes suffer two significant drawbacks: outliers are a priori assumed unlikely, and the posterior variance conditioned on observed data depends only on the locations of those data, not the associated sample values. Student's-T processes are a generalization of Gaussian processes, founded on the Student's-T distribution instead of the Gaussian distribution. Student's-T processes maintain the primary advantages of Gaussian processes (kernel function, analytic update rule) with additional benefits beyond Gaussian processes. The Student's-T distribution has higher Kurtosis than a Gaussian distribution and so outliers are much more likely, and the posterior variance increases or decreases depending on the variance of observed data sample values. Here, we describe Student's-T processes, and discuss their advantages in the context of aerospace optimization. We show how to construct a Student's-T process using a kernel function and how to update the process given new samples. We provide a clear derivation of optimization-relevant quantities such as expected improvement, and contrast with the related computations for Gaussian processes. Finally, we compare the performance of Student's-T processes against Gaussian process on canonical test problems in Bayesian optimization, and apply the Student's-T process to the optimization of an aerostructural design problem.


title: Characterization of probability distribution convergence in Wasserstein distance by $L^{p}$-quantization error function

abstract: We establish conditions to characterize probability measures by their $L^{p}$-quantization error functions in both $\mathbb{R}^{d}$ and Hilbert settings. This characterization is two-fold: static (identity of two distributions) and dynamic (convergence for the $L^p$-Wasserstein distance). We first propose a criterion on the quantization level $N$, valid for any norm on $\mathbb{R}^{d}$ and any order $p$ based on a geometrical approach involving the Vorono\"i diagram. Then, we prove that in the $L^2$-case on a (separable) Hilbert space, the condition on the level $N$ can be reduced to $N=2$, which is optimal. More quantization based characterization cases on dimension 1 and a discussion of the completeness of a distance defined by the quantization error function can be found in the end of this paper.


title: Observers in Kerr spacetimes: the ergoregion on the equatorial plane

abstract: We perform a detailed analysis of the properties of stationary observers located on the equatorial plane of the ergosphere in a Kerr spacetime, including light-surfaces. This study highlights crucial differences between black hole and the super-spinner sources. In the case of Kerr naked singularities, the results allow us to distinguish between "weak" and "strong" singularities, corresponding to spin values close to or distant from the limiting case of extreme black holes, respectively. We derive important limiting angular frequencies for naked singularities. We especially study very weak singularities as resulting from the spin variation of black holes. We also explore the main properties of zero angular momentum observers for different classes of black hole and naked singularity spacetimes.


title: Jamming of Deformable Polygons

abstract: There are two main classes of physics-based models for two-dimensional cellular materials: packings of repulsive disks and the vertex model. These models have several disadvantages. For example, disk interactions are typically a function of particle overlap, yet the model assumes that the disks remain circular during overlap. The shapes of the cells can vary in the vertex model, however, the packing fraction is fixed at $\phi=1$. Here, we describe the deformable particle model (DPM), where each particle is a polygon composed of a large number of vertices. The total energy includes three terms: two quadratic terms to penalize deviations from the preferred particle area $a_0$ and perimeter $p_0$ and a repulsive interaction between DPM polygons that penalizes overlaps. We performed simulations to study the onset of jamming in packings of DPM polygons as a function of asphericity, ${\cal A} = p_0^2/4\pi a_0$. We show that the packing fraction at jamming onset $\phi_J({\cal A})$ grows with increasing ${\cal A}$, reaching confluence at ${\cal A} \approx 1.16$. ${\cal A}^*$ corresponds to the value at which DPM polygons completely fill the cells obtained from a surface-Voronoi tessellation. Further, we show that DPM polygons develop invaginations for ${\cal A} > {\cal A}^*$ with excess perimeter that grows linearly with ${\cal A}-{\cal A}^*$. We confirm that packings of DPM polygons are solid-like over the full range of ${\cal A}$ by showing that the shear modulus is nonzero.


title: An estimation of level sets for non local KPP equations with delay

abstract: We study the large time asymptotic behavior of the solutions of the linear parabolic equation with delay $(*)$: $u_{t}(t,x) = u_{xx}(t,x) - u(t,x) + \int_{\mathbb{R}} k(x-y) \, u (t-h, y)\, dy$, $x \in \R$, $\ t >0$, and $k(x) \in L^1(\R)$. As an application we get estimates on the measure of level sets of non local KPP type equations with delay. For this type of nonlinear equations we prove that, in contrast with the classical case, the solution to the initial value problem with data of compact support may not be persistent.


title: The Of?p stars of the Magellanic Clouds: Are they strongly magnetic?

abstract: All known Galactic Of?p stars have been shown to host strong, organized, magnetic fields. Recently, five Of?p stars have been discovered in the Magellanic Clouds. They posses photometric \citep{Naze} and spectroscopic \citep{Walborn} variability compatible with the Oblique Rotator Model (ORM). However, their magnetic fields have yet to be directly detected. We have developed an algorithm allowing for the synthesis of photometric observables based on the Analytic Dynamical Magnetosphere (ADM) model of \citet{Owocki}. We apply our model to OGLE photometry in order to constrain their magnetic geometries and surface dipole strengths. We predict that the field strengths for some of these candidate extra-Galactic magnetic stars may be within the detection limits of the FORS2 instrument.


title: LCD: Low Latency Command Dissemination for A Platoon of Vehicles

abstract: In a vehicular platoon, a lead vehicle that is responsible for managing the platoon's moving directions and velocity periodically disseminates control commands to following vehicles based on vehicle-to-vehicle communications. However, reducing command dissemination latency with multiple vehicles while ensuring successful message delivery to the tail vehicle is challenging. We propose a new linear dynamic programming algorithm using backward induction and interchange arguments to minimize the dissemination latency of the vehicles. Furthermore, a closed form of dissemination latency in vehicular platoon is obtained by utilizing Markov chain with M/M/1 queuing model. Simulation results confirm that the proposed dynamic programming algorithm improves the dissemination rate by at least 50.9%, compared to similar algorithms in the literature. Moreover, it also approximates the best performance with the maximum gap of up to 0.2 second in terms of latency.


title: New insights into the nature of transition disks from a complete disk survey of the Lupus star forming region

abstract: Transition disks with large dust cavities around young stars are promising targets for studying planet formation. Previous studies have revealed the presence of gas cavities inside the dust cavities hinting at recently formed, giant planets. However, many of these studies are biased towards the brightest disks in the nearby star forming regions, and it is not possible to derive reliable statistics that can be compared with exoplanet populations. We present the analysis of 11 transition disks with large cavities (>20 AU radius) from a complete disk survey of the Lupus star forming region, using ALMA Band 7 observations at 0.3" (22-30 AU radius) resolution of the 345 GHz continuum, 13CO and C18O 3-2 observations and the Spectral Energy Distribution of each source. Gas and dust surface density profiles are derived using the physical-chemical modeling code DALI. This is the first study of transition disks of large cavities within a complete disk survey within a star forming region. The dust cavity sizes range from 20-90 AU radius and in three cases, a gas cavity is resolved as well. The deep drops in gas density and large dust cavity sizes are consistent with clearing by giant planets. The fraction of transition disks with large cavities in Lupus is ~11%, which is inconsistent with exoplanet population studies of giant planets at wide orbits. Furthermore, we present a hypothesis of an evolutionary path for large massive disks evolving into transition disks with large cavities.


title: Constraints on the radiation temperature before inflation

abstract: We consider the short period of cosmic expansion ranging from the end of the Planck era to the beginning of inflation and set upper and lower limits on the temperature of the radiation at the commencement of the inflationary phase.


title: The Rayleigh-Taylor instability for the Verigin problem with and without phase transition

abstract: Isothermal compressible two-phase flows in a capillary are modeled with and without phase transition in the presence of gravity, employing Darcy's law for the velocity field. It is shown that the resulting systems are thermodynamically consistent in the sense that the available energy is a strict Lyapunov functional. In both cases, the equilibria with flat interface are identified. It is shown that the problems are well-posed in an $L_p$-setting and generate local semiflows in the proper state manifolds. The main result concerns the stability of equilibria with flat interface, i.e. the Rayleigh-Taylor instability.


title: Iron-peak elements Sc, V, Mn, Cu and Zn in Galactic bulge globular clusters

abstract: Globular clusters are tracers of the history of star formation and chemical enrichment in the early Galaxy. Their abundance pattern can help understanding their chemical enrichment processes. In particular, the iron-peak elements have been relatively little studied so far in the Galactic bulge. The abundances of iron-peak elements are derived for the sample clusters, and compared with bulge field, and thick disk stars. We derive abundances of the iron-peak elements Sc, V, Mn, Cu, and Zn in individual stars of five bulge globular clusters (NGC 6528, NGC 6553, NGC 6522, NGC 6558, HP1), and of the reference thick disk/inner halo cluster 47 Tucanae (NGC 104). High resolution spectra were obtained with the UVES spectrograph at the Very Large Telescope over the years. The sample globular clusters studied span metallicities in the range -1.2<Fe/H]<0.0. V and Sc appear to vary in lockstep with Fe, indicating that they are produced in the same supernovae as Fe. We find that Mn is deficient in metal-poor stars, confirming that it is underproduced in massive stars; Mn-over-Fe steadily increases at the higher metallicities due to a metallicity-dependent enrichment by supernovae of type Ia. Cu behaves as a secondary element, indicating its production in a weak-s process in massive stars. Zn has an alpha-like behaviour at low metallicities, which can be explained in terms of nucleosynthesis in hypernovae. At the metal-rich end, Zn decreases with increasing metallicity, similarly to the alpha-elements.


title: Local free energy approximations for the coarse-graining of adsorption phenomena

abstract: We investigate the coarse-graining of host-guest systems under the perspective of the local distribution of pore occupancies, along with the physical meaning and actual computability of the coarse-interaction terms. We show that the widely accepted approach, in which the contributions to the free energy given by the molecules located in two neighboring pores are estimated through Monte Carlo simulations where the two pores are kept separated from the rest of the system, leads to inaccurate results at high sorbate densities. In the coarse-graining strategy that we propose, which is based on the Bethe-Peierls approximation, density-independent interaction terms are instead computed according to local effective potentials that take into account the correlations between the pore pair and its surroundings by means of mean-field correction terms, without the need of simulating the pore pair separately. Use of the interaction parameters obtained this way allows the coarse-grained system to reproduce more closely the equilibrium properties of the original one. Results are shown for lattice-gases where the local free energy can be computed exactly, and for a system of Lennard-Jones particles under the effect of a static confining field.


title: When Does Stochastic Gradient Algorithm Work Well?

abstract: In this paper, we consider a general stochastic optimization problem which is often at the core of supervised learning, such as deep learning and linear classification. We consider a standard stochastic gradient descent (SGD) method with a fixed, large step size and propose a novel assumption on the objective function, under which this method has the improved convergence rates (to a neighborhood of the optimal solutions). We then empirically demonstrate that these assumptions hold for logistic regression and standard deep neural networks on classical data sets. Thus our analysis helps to explain when efficient behavior can be expected from the SGD method in training classification models and deep neural networks.


title: Indication of anisotropy in arrival directions of ultra-high-energy cosmic rays through comparison to the flux pattern of extragalactic gamma-ray sources

abstract: A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.