title: On Human Predictions with Explanations and Predictions of Machine Learning Models: A Case Study on Deception Detection

abstract: Humans are the final decision makers in critical tasks that involve ethical and legal concerns, ranging from recidivism prediction, to medical diagnosis, to fighting against fake news. Although machine learning models can sometimes achieve impressive performance in these tasks, these tasks are not amenable to full automation. To realize the potential of machine learning for improving human decisions, it is important to understand how assistance from machine learning models affects human performance and human agency. In this paper, we use deception detection as a testbed and investigate how we can harness explanations and predictions of machine learning models to improve human performance while retaining human agency. We propose a spectrum between full human agency and full automation, and develop varying levels of machine assistance along the spectrum that gradually increase the influence of machine predictions. We find that without showing predicted labels, explanations alone slightly improve human performance in the end task. In comparison, human performance is greatly improved by showing predicted labels (>20% relative improvement) and can be further improved by explicitly suggesting strong machine performance. Interestingly, when predicted labels are shown, explanations of machine predictions induce a similar level of accuracy as an explicit statement of strong machine performance. Our results demonstrate a tradeoff between human performance and human agency and show that explanations of machine predictions can moderate this tradeoff.

categories: cs.AI cs.CL cs.CY physics.soc-ph stat.ML


title: Ratio of black hole to galaxy mass of an extremely red dust-obscured galaxy at z = 2.52

abstract: We present a near-infrared (NIR) spectrum of WISE J1042+1641, an extremely red dust-obscured galaxy (DOG), which has been observed with the LIRIS on the 4.2m William Hershel Telescope. This object was selected as a hyper-luminous DOG candidate at z ~ 2 by combining the optical and IR photometric data based on the SDSS and WISE, although its redshift had not yet been confirmed. Based on the LIRIS observation, we confirmed its redshift of 2.521 and total IR luminosity of log(L_IR/L_sun) = 14.57, which satisfies the criterion for an extremely luminous IR galaxy (ELIRG). Moreover, we indicate that this object seems to have an extremely massive black hole with M_BH = 10^10.92 M_sun based on the broad Halpha line: the host stellar mass is derived as M_star = 10^13.55 M_sun by a fit of the spectral energy distribution. Very recently, it has been reported that this object is an anomalous gravitationally lensed quasar based on near-IR high-resolution imaging data obtained with the Hubble Space Telescope. Its magnification factor has also been estimated with some uncertainty (i.e., mu = 53-122). We investigate the ratio of the black hole to galaxy mass, which is less strongly affected by a lensing magnification factor, instead of the absolute values of the luminosities and masses. We find that the M_BH/M_star ratio (i.e., 0.0140-0.0204) is significantly higher than the local relation, following a sequence of unobscured quasars instead of obscured objects (e.g., submillimeter galaxies) at the same redshift. Moreover, the LIRIS spectrum shows strongly blueshifted oxygen lines with an outflowing velocity of ~ 1100 km/s, and our Swift X-ray observation also supports that this source is an absorbed AGN with an intrinsic column density of N_H = 4.9 x 10^23 cm^-2. These results imply that WISE J1042+1641 is in a blow-out phase at the end of the buried rapid black hole growth.


title: Exact localized and ballistic eigenstates in disordered chaotic spin ladders and the Fermi-Hubbard model

abstract: We demonstrate the existence of exact atypical many-body eigenstates in a class of disordered, interacting one-dimensional quantum systems that includes the Fermi-Hubbard model as a special case. These atypical eigenstates, which generically have finite energy density and are exponentially many in number, are populated by noninteracting excitations. They can exhibit Anderson localization with area-law eigenstate entanglement or, surprisingly, ballistic transport at any disorder strength. These properties differ strikingly from those of typical eigenstates nearby in energy, which we show give rise to diffusive transport as expected in a chaotic quantum system. We discuss how to observe these atypical eigenstates in cold-atom experiments realizing the Fermi-Hubbard model, and comment on the robustness of their properties.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.quant-gas nlin.CD quant-ph


title: Cosmic Reionization History and Dark Matter Scenarios

abstract: We perform an analysis of the cosmic reionization in the standard cold dark matter (CDM) paradigm and in alternative dark matter scenarios. Building upon the work of Corasaniti et al. (2017), we predict the reionization history for CDM, for warm dark matter (WDM), late-forming dark matter (LFDM) and ultra-light axion dark matter (ADM) models which reproduce state-of-art measurements of the galaxy luminosity function at very high-redshifts $6\le z\le 10$. To this purpose we adopt a reionization model parametrized in terms of the limiting UV-magnitude of galaxies contributing to the reionization $M_{\rm lim}$ and the average effective escape fraction of UV photons reaching the intergalactic medium $\tilde{f}$. For each DM model we compute the redshift evolution of the Thomson scattering optical depth $\tau_e(z)$ and the comoving ionization fraction $Q_{\rm HII}(z)$. We find the different DM models to have similar reionization histories. Differences with respect to the CDM case increase at fainter limiting UV-magnitudes and are degenerate with the effect of varying the reionization model parameters. Using Planck's determination of the integrated optical depth in combination with measurements of the neutral hydrogen fraction at different redshifts, we infer constraints on $\tilde{f}$ and $M_{\rm lim}$. The results are largely independent of the assumed DM scenario, in particular for $M_{\rm lim}\gtrsim -13$ we obtain that the effective escape fraction lies in the range $0.07\lesssim \tilde{f}\lesssim 0.15$ at $2\sigma$.


title: Inclusive heavy flavor jet production with semi-inclusive jet functions: from proton to heavy-ion collisions

abstract: The past several years have witnessed important developments in the QCD theory of jet production and jet substructure in hadronic collisions. In the framework of soft-collinear effective theory, semi-inclusive jet functions and semi-inclusive fragmenting jet functions have allowed us to combine higher order calculations with resummation of potentially large logarithms of the jet radius, $\ln R$. Very recently, the semi-inclusive jet functions for partons fragmenting into heavy flavor jets were computed by Dai, Kim and Leibovich. In this paper we show how the formalism can be extended to c-jet and b-jet production in heavy ion collisions. The semi-inclusive jet functions for heavy flavor jets in a QCD medium are evaluated up to the next-to-leading order in $\alpha_s$ and first order in opacity. For phenomenological applications, we also consider the inclusion of the cold nuclear matter effects and the jet energy dissipation due to collisional interactions in matter. We present the numerical predictions for the cross sections and the corresponding nuclear modification factors in proton-nucleus and nucleus-nucleus collisions and compare our results to data from the Large Hadron Collider.


title: Vector-Boson Fusion Higgs Pair Production at N$^3$LO

abstract: We calculate the next-to-next-to-next-to-leading order (N$^3$LO) QCD corrections to vector-boson fusion (VBF) Higgs pair production in the limit in which there is no partonic exchange between the two protons. We show that the inclusive cross section receives negligible corrections at this order, while the scale variation uncertainties are reduced by a factor four. We present differential distributions for the transverse momentum and rapidity of the final state Higgs bosons, and show that there is almost no kinematic dependence to the third order corrections. Finally we study the impact of deviations from the Standard Model in the trilinear Higgs coupling, and show that the structure of the higher order corrections does not depend on the self-coupling. These results are implemented in the latest release of the proVBFH-incl program.


title: Mergers, tidal interactions, and mass exchange in a population of disc globular clusters: II. Long-term evolution

abstract: Globular clusters (GCs), the oldest stellar systems observed in the Milky Way, have for long been considered single stellar populations. As such, they provided an ideal laboratory to understand stellar dynamics and primordial star formation processes. However, during the last two decades, observations unveiled their real, complex nature. Beside their pristine stars, GCs host one or more helium enriched and possibly younger stellar populations whose formation mechanism is still unknown. Even more puzzling is the existence of GCs showing star by star iron spreads. Using detailed N-body simulations we explore the hypothesis that these anomalies in metallicity could be the result of mutual stripping and mergers between a primordial population of disc GCs. In the first paper of this series we proved, both with analytical arguments and short-term N-body simulations, that disc GCs have larger fly-by and close encounter rates with respect to halo clusters. These interactions lead to mass exchange and even mergers that form new GCs, possibly showing metallicity spreads. Here, by means of long-term direct N-body simulations, we provide predictions on the dynamical properties of GCs that underwent these processes. The comparison of our predictions with available and future observational data could provide insights on the origin of GCs and on the Milky Way build-up history as a whole.


title: A Spin-2 Conjecture on the Swampland

abstract: We consider effective theories with massive fields that have spins larger than or equal to two. We conjecture a universal cutoff scale on any such theory that depends on the lightest mass of such fields. This cutoff corresponds to the mass scale of an infinite tower of states, signalling the breakdown of the effective theory. The cutoff can be understood as the Weak Gravity Conjecture applied to the St\"uckelberg gauge field in the mass term of the high spin fields. A strong version of our conjecture applies even if the graviton itself is massive, so to massive gravity. We provide further evidence for the conjecture from string theory.


title: Synthetic nebular emission from massive galaxies II: ultraviolet-line diagnostics of dominant ionizing sources

abstract: We compute synthetic optical and ultraviolet (UV) emission-line properties of galaxies in a full cosmological framework by coupling, in post-processing, new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies. Our self-consistent modelling accounts for nebular emission from young stars and accreting black holes (BHs). We investigate which optical- and UV-line diagnostic diagrams can best help to discern between the main ionizing sources, as traced by the ratio of BH accretion to star formation rates in model galaxies, over a wide range of redshifts. At low redshift, simulated star-forming galaxies, galaxies dominated by active galactic nuclei and composite galaxies are appropriately differentiated by standard selection criteria in the classical [OIII]$\lambda$5007/H$\beta$ versus [NII]$\lambda$6584/H$\alpha$ diagram. At redshifts $z \gt 1$, however, this optical diagram fails to discriminate between active and inactive galaxies at metallicities below $0.5\ Z_\odot$. To robustly classify the ionizing radiation of such metal-poor galaxies, which dominate in the early Universe, we confirm 3 previous, and propose 11 novel diagnostic diagrams based on equivalent widths and luminosity ratios of UV emission lines, such as EW(OIII]$\lambda$1663) versus OIII]$\lambda$1663/HeII$\lambda$1640, CIII]$\lambda$1908/HeII$\lambda$1640 versus OIII]$\lambda$1663/HeII$\lambda$1640, and CIV$\lambda$1550/CIII]$\lambda$1908 versus CIII]$\lambda$1908/CII$\lambda$2326. We formulate associated UV selection criteria and discuss some caveats of our results (e.g., uncertainties in the modelling of the HeII$\lambda$1640 line). These UV diagnostic diagrams are potentially important for the interpretation of high-quality spectra of very distant galaxies to be gathered by next-generation telescopes, such as the James Webb Space Telescope.


title: SO(10) paths to dark matter

abstract: The grand-unification gauge group SO(10) contains matter parity as a discrete subgroup. This symmetry could be at the origin of dark matter stability. The properties of the dark matter candidates depend on the path along which SO(10) is broken, in particular through Pati-Salam or left-right symmetric subgroups. We systematically determine the non-supersymmetric dark matter scenarios that can be realized along the various paths. We emphasize that the dark matter candidates may have colored or electrically charged partners at low scale that belong to the same SO(10) multiplet. These states, which in many cases are important for co-annihilation, could be observed more easily than the dark matter particle. We determine the structure of the tree-level and loop-induced mass splittings between the dark matter candidate and their partners and discuss the possible phenomenological implications.


title: Measuring the local matter density using Gaia DR2

abstract: We determine the total dynamical matter density in the solar neighbourhood using the second Gaia data release (DR2). The dynamical matter density distribution is inferred in a framework of a Bayesian hierarchical model, which accounts for position and velocity of all individual stars, as well as the full error covariance matrix of astrometric observables, in a joint fit of the vertical velocity distribution and stellar number density distribution. This was done for eight separate data samples, with different cuts in observed absolute magnitude, each containing about 25,000 stars. The model for the total matter density does not rely on any underlying baryonic model, although we assumed that it is symmetrical, smooth, and monotonically decreasing with distance from the mid-plane. We infer a density distribution which is strongly peaked in the region close to the Galactic plane ($\lesssim 60$ pc), for all eight stellar samples. Assuming a baryonic model and a dark matter halo of constant density, this corresponds to a surplus surface density of approximately 5--9 M$_\odot$/pc$^2$. For the Sun's position and vertical velocity with respect to the Galactic plane, we infer $Z_\odot = 4.76\pm 2.27$ pc and $W_\odot = 7.24\pm 0.19$ km/s. These results suggest a surplus of matter close to the Galactic plane, possibly explained by an underestimated density of cold gas. We discuss possible systematic effects that could bias our result, for example unmodelled non-equilibrium effects, and how to account for such effects in future extensions of this work.


title: Developing the MeV potential of DUNE: Detailed considerations of muon-induced spallation and other backgrounds

abstract: The Deep Underground Neutrino Experiment (DUNE) could be revolutionary for MeV neutrino astrophysics, because of its huge detector volume, unique event reconstruction capabilities, and excellent sensitivity to the $\nu_e$ flavor. However, its backgrounds are not yet known. A major background is expected due to muon spallation of argon, which produces unstable isotopes that later beta decay. We present the first comprehensive study of MeV spallation backgrounds in argon, detailing isotope production mechanisms and decay properties, analyzing beta energy and time distributions, and proposing experimental cuts. We show that above a nominal detection threshold of 5-MeV electron energy, the most important backgrounds are --- surprisingly --- due to low-A isotopes, such as Li, Be, and B, even though high-A isotopes near argon are abundantly produced. We show that spallation backgrounds can be powerfully rejected by simple cuts, with clear paths for improvements. We compare these background rates to rates of possible MeV astrophysical neutrino signals in DUNE, including solar neutrinos (detailed in a companion paper [Capozzi et al. arXiv:1808.08232 [hep-ph]]), supernova burst neutrinos, and the diffuse supernova neutrino background. Further, to aid trigger strategies, in the Appendixes we quantify the rates of single and multiple MeV events due to spallation, radiogenic neutron capture, and other backgrounds, including through pileup. Our overall conclusion is that DUNE has high potential for MeV neutrino astrophysics, but reaching this potential requires new experimental initiatives.


title: Non-perturbative probability distribution function for cosmological counts in cells

abstract: We present a non-perturbative calculation of the 1-point probability distribution function (PDF) for the spherically-averaged matter density field. The PDF is represented as a path integral and is evaluated using the saddle-point method. It factorizes into an exponent given by a spherically symmetric saddle-point solution and a prefactor produced by fluctuations. The exponent encodes the leading sensitivity of the PDF to the dynamics of gravitational clustering and statistics of the initial conditions. In contrast, the prefactor has only a weak dependence on cosmology. It splits into a monopole contribution which is evaluated exactly, and a factor corresponding to aspherical fluctuations. The latter is crucial for the consistency of the calculation: neglecting it would make the PDF incompatible with translational invariance. We compute the aspherical prefactor using a combination of analytic and numerical techniques. We demonstrate the factorization of spurious enhanced contributions of large bulk flows and their cancellation due the equivalence principle. We also identify the sensitivity to the short-scale physics and argue that it must be properly renormalized. The uncertainty associated with the renormalization procedure gives an estimate of the theoretical error. For zero redshift, the precision varies from sub percent for moderate density contrasts to tens of percent at the tails of the distribution. It improves at higher redshifts. We compare our results with N-body simulation data and find an excellent agreement.


title: The topological objects near the chiral crossover transition in QCD

abstract: We study the underlying topology of gauge fields in 2+1 flavor QCD with domain wall fermions on lattices of size $32^3\times 8$, at and immediately above the chiral crossover transition. Using valence overlap fermions with exact index theorem, we focus on its zero modes for different choices of periodicity phases along the temporal direction. Our studies show that the zero modes are due to fractionally charged topological objects, the instanton-dyons. We further provide qualitative study of the interactions between those and compare with the available semi-classical results, finding remarkably accurate agreement in all cases.


title: Determining Metastable Ion Lifetime and History through Wave-Particle Interaction

abstract: Laser-induced fluorescence (LIF) performed on metastable ions is frequently used to probe the dynamics of ground-state ion motions in many laboratory plasmas. However, these measurements place restrictions on the metastable ion lifetime. Metastable states are produced from direct ionization of neutral atoms as well as ions in other electronic states, of which the former will only faithfully represent processes that act on the ion dynamics in a time shorter than the metastable lifetime. We present here the first experimental study of this type of systematic effect using wave-particle interaction in an argon multidipole plasma. The metastable lifetime and relative fraction of metastables produced from preexisting ions, necessary for correcting the LIF measurement errors, can be determined by fitting the experimental results with the theory we propose.


title: Excitation and depletion of the asteroid belt in the early instability scenario

abstract: Containing only a few percent the mass of the moon, the current asteroid belt is around three to four orders of magnitude smaller that its primordial mass inferred from disk models. Yet dynamical studies have shown that the asteroid belt could not have been depleted by more than about an order of magnitude over the past ~4 Gyr. The remainder of the mass loss must have taken place during an earlier phase of the solar system's evolution. An orbital instability in the outer solar system occurring during the process of terrestrial planet formation can reproduce the broad characteristics of the inner solar system. Here, we test the viability of this model within the constraints of the main belt's low present-day mass and orbital structure. While previous studies modeled asteroids as massless test particles because of limited computing power, our work uses GPU (Graphics Processing Unit) acceleration to model a fully self-gravitating asteroid belt. We find that depletion in the main belt is related to the giant planets' exact evolution within the orbital instability. Simulations that produce the closest matches to the giant planets' current orbits deplete the main belt by two to three orders of magnitude. These simulated asteroid belts are also good matches to the actual asteroid belt in terms of their radial mixing and broad orbital structure.


title: Anisotropic stars as ultracompact objects in General Relativity

abstract: Anisotropic stresses are ubiquitous in nature, but their modeling in General Relativity is poorly understood and frame dependent. We introduce the first study on the dynamical properties of anisotropic self-gravitating fluids in a covariant framework. Our description is particularly useful in the context of tests of the black hole paradigm, wherein ultracompact objects are used as black hole mimickers but otherwise lack a proper theoretical framework. We show that: (i) anisotropic stars can be as compact and as massive as black holes, even for very small anisotropy parameters; (ii) the nonlinear dynamics of the 1+1 system is in good agreement with linearized calculations, and shows that configurations below the maximum mass are nonlinearly stable; (iii) strongly anisotropic stars have vanishing tidal Love numbers in the black-hole limit; (iv) their formation will usually be accompanied by gravitational-wave echoes at late times.


title: Fully differential Vector-Boson Fusion Higgs Pair Production at Next-to-Next-to-Leading Order

abstract: We calculate the fully differential next-to-next-to-leading order (NNLO) QCD corrections to vector-boson fusion (VBF) Higgs pair production. This calculation is achieved in the limit in which there is no colored cross-talk between the colliding protons, using the projection-to-Born method. We present differential cross sections of key observables, showing corrections of up to 3-4% at this order after typical VBF cuts, with the total cross section receiving contributions of about 2%. In contrast to single Higgs VBF production, we find that the NNLO corrections are for the most part within the next-to-leading order scale uncertainty bands.


title: Two Terrestrial Planet Families With Different Origins

abstract: The potentially important role of stellar irradiation in envelope removal for planets with diameters of $\lessapprox$ 2 R$_{\Earth}$ has been inferred both through theoretical work and the observed bimodal distribution of small planet occurrence as a function of radius. We examined the trends for small planets in the three-dimensional radius-insolation-density space and find that the terrestrial planets divide into two distinct families based on insolation. The lower insolation family merges with terrestrial planets and small bodies in the solar system and is thus Earth-like. The higher insolation terrestrial planet family forms a bulk-density continuum with the sub-Neptunes, and is thus likely to be composed of remnant cores produced by photoevaporation. Based on the density-radius relationships, we suggest that both terrestrial families show evidence of density enhancement through collisions. Our findings highlight the important role that both photoevaporation and collisions have in determining the density of small planets.


title: The Mono-Tau Menace: From $B$ Decays to High-$p_T$ Tails

abstract: We investigate the crossing-symmetry relation between $b\to c\tau^-\bar\nu$ decay and $b\bar c\to \tau^-\bar\nu$ scattering to derive direct correlations of New Physics in semi-tauonic $B$-meson decays and the mono-tau signature at the LHC ($pp\to\tau_h X$ + MET). Using an exhaustive set of effective operators and heavy mediators we find that the current ATLAS and CMS data constrain scenarios addressing anomalies in $B$-decays. Pure tensor solutions, completed by leptoquark, and right-handed solutions, completed by $W^\prime_R$ or leptoquark, are challenged by our analysis. Furthermore, the sensitivity that will be achieved in the high-luminosity phase of the LHC will probe $all$ the possible scenarios that explain the anomalies. Finally, we note that the LHC is also competitive in the $b\to u$ transitions and bounds in some cases are currently better than those from $B$ decays.


title: Control and Coherence Time Enhancement of the 0-$\pi $ Qubit

abstract: Kitaev's 0-$\pi $ qubit encodes quantum information in two protected, near-degenerate states of a superconducting quantum circuit. In a recent work, we have shown that the coherence times of a realistic 0-$\pi $ device can surpass that of today's best superconducting qubits [Groszkowski et al., New Journal of Physics 20 043053 (2018)]. Here we address controllability of the 0-$\pi $ qubit. Specifically, we investigate the potential for dispersive control and readout, and introduce a new, fast and high-fidelity single-qubit gate that can interpolate smoothly between logical X and Z. We characterize the action of this gate using a multi-level treatment of the device, and analyze the impact of circuit element disorder and deviations in control and circuit parameters from their optimal values. Furthermore, we propose a cooling scheme to decrease the photon shot-noise dephasing rate, which we previously found to limit the coherence times of 0-$\pi $ devices within reach of current experiments. Using this approach, we predict coherence time enhancements between one and three orders of magnitude, depending on parameter regime.


title: Integrability-protected adiabatic reversibility in quantum spin chains

abstract: We consider the out-of-equilibrium dynamics of the Heisenberg anisotropic quantum spin--$1/2$ chain threaded by a time-dependent magnetic flux. In the spirit of the recently developed generalized hydrodynamics (GHD), we exploit the integrability of the model for any flux values to derive an exact description of the dynamics in the limit of slowly varying flux: the state of the system is described at any time by a time-dependent generalized Gibbs ensemble. Two dynamical regimes emerge according to the value of the anisotropy $\Delta$. For $|\Delta| > 1$, reversibility is preserved: the initial state is always recovered whenever the flux is brought back to zero. On the contrary, for $|\Delta| < 1$, instabilities of quasiparticles produce irreversible dynamics as confirmed by the dramatic growth of entanglement entropy. In this regime, the standard GHD description becomes incomplete and we complement it via a maximum entropy production principle. We test our predictions against numerical simulations finding excellent agreement.


title: Testing quantum decoherence at DUNE

abstract: We address some theoretical issues of the quantum decoherence phenomenon within the neutrino oscillation framework and carry out various tests under DUNE simulated experimental environment. On the theoretical side, we provide a general expression for an invariant decoherence matrix under a quantum basis rotation. On the simulated experimental side, considering a rotation invariant and a non-invariant decoherence matrix, we study the impact on the fitting of the standard oscillation parameters, the sensitivity in the mass hierarchy and CP violation, combining the neutrino and antineutrino mode and all available neutrino oscillation probabilities channels. Furthermore, a sensitivity for the decoherence parameter of order $10^{-24}\text{GeV}$ at 3$\sigma$, is obtained for our best case.


title: The Navier-Stokes-Vlasov-Fokker-Planck system as a scaling limit of particles in a fluid

abstract: Convergence of a system of particles, interacting with a fluid, to Navier-Stokes-Vlasov-Fokker-Planck system is studied. The interaction between particles and fluid is described by Stokes drag force. The empirical measure of particles is proved to converge to the Vlasov-Fokker-Planck component of the system and the velocity of the fluid coupled with the particles converges in the uniform topology to the the Navier-Stokes component. A new uniqueness result for the PDE system is added.


title: Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D

abstract: We present an asymptotically and unconditionally stable numerical method to account for the momentum transfer between multiple species. Momentum is conserved to machine precision. This implies that the asymptotic equilibrium corresponds to the velocity of the center of mass. Aimed at studying dust dynamics, we implement this numerical method in the publicly available code FARGO3D. To validate our implementation, we develop a test suite for an arbitrary number of species, based on analytical or exact solutions of problems related to perfect damping, damped sound waves, shocks, local and global gas-dust radial drift in a disk and, linear streaming instability. In particular, we obtain first-order, steady-state solutions for the radial drift of multiple dust species in protoplanetary disks, in which the pressure gradient is not necessarily small. We additionally present non-linear shearing-box simulations of the streaming instability and compare them with previous results obtained with Lagrangian particles. We successfully validate our implementation by recovering the solutions from the test suite to second- and first-order accuracy in space and time, respectively. From this, we conclude that our scheme is suitable, and very robust, to study the self-consistent dynamics of several fluids. In particular, it can be used for solving the collisions between gas and dust in protoplanetary disks, with any degree of coupling.


title: The LOFAR Two-metre Sky Survey - II. First data release

abstract: The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

categories: astro-ph.GA astro-ph.CO astro-ph.HE astro-ph.IM


title: The LOFAR Two-metre Sky Survey (LoTSS) III. First Data Release: optical/IR identifications and value-added catalogue

abstract: The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release.


title: The LOFAR Two-metre Sky Survey IV. First Data Release: Photometric redshifts and rest-frame magnitudes

abstract: The LOFAR Two-metre Sky Survey (LoTSS) is a sensitive, high-resolution 120-168 MHz survey of the Northern sky. The LoTSS First Data Release (DR1) presents 424 square degrees of radio continuum observations over the HETDEX Spring Field (10h45m00s $<$ right ascension $<$ 15h30m00s and 45$^\circ$00$'$00$'$ $<$ declination $<$ 57$^\circ$00$'$00$''$) with a median sensitivity of 71$\mu$Jy/beam and a resolution of 6$''$. In this paper we present photometric redshifts (photo-$z$) for 94.4% of optical sources over this region that are detected in the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) 3$\pi$ steradian survey. Combining the Pan-STARRS optical data with mid-infrared photometry from the Wide-field Infrared Survey Explorer, we estimate photo-$z$s using a novel hybrid photometric redshift methodology optimised to produce the best possible performance for the diverse sample of radio continuum selected sources. For the radio-continuum detected population, we find an overall scatter in the photo-$z$ of 3.9% and an outlier fraction ($\left | z_{\rm{phot}} - z_{\rm{spec}} \right | / (1+z_{\rm{spec}}) > 0.15$) of 7.9%. We also find that, at a given redshift, there is no strong trend in photo-$z$ quality as a function of radio luminosity. However there are strong trends as a function of redshift for a given radio luminosity, a result of selection effects in the spectroscopic sample and/or intrinsic evolution within the radio source population. Additionally, for the sample of sources in the LoTSS First Data Release with optical counterparts, we present rest-frame optical and mid-infrared magnitudes based on template fits to the consensus photometric (or spectroscopic when available) redshift.


title: Evolutionary phases of merging clusters as seen by LOFAR

abstract: Massive, merging galaxy clusters often host giant, diffuse radio sources that arise from shocks and turbulence; hence, radio observations can be useful for determining the merger state of a cluster. In preparation for a larger study, we selected three clusters -- Abell 1319, Abell 1314, and RXC J1501.3+4220 (Z7215) -- making use of the new LOFAR Two-Metre Sky Survey (LoTSS) at 120-168 MHz, and together with archival data, show that these clusters appear to be in pre-merging, merging, and post-merging states, respectively. We argue that Abell 1319 is likely in its pre-merging phase, where three separate cluster components are about to merge. There are no radio halos nor radio relics detected in this system. Abell 1314 is a highly-disturbed, low-mass cluster which is likely in the process of merging. This low-mass system does not show a radio halo, however, we argue that the merger activates mechanisms that cause electron re-acceleration in the large 800 kpc radio tail associated with IC~711. In the cluster Z7215 we discover diffuse radio emission at the cluster center, and we classify this emission as a radio halo, although it is dimmer and smaller than expected by the radio halo power versus cluster mass correlation. We suggest that the disturbed cluster Z7215 is in its post-merging phase. Systematic studies of this kind over a larger sample of clusters observed with LoTSS will help constrain the time scales involved in turbulent re-acceleration and the subsequent energy losses of the underlying electrons.


title: The spectacular cluster chain Abell 781 as observed with LOFAR, GMRT, and XMM-Newton

abstract: A number of merging galaxy clusters show the presence of large-scale radio emission associated with the intra-cluster medium (ICM). These synchrotron sources are generally classified as radio haloes and radio relics. Whilst it is commonly accepted that mergers play a crucial role in the formation of radio haloes and relics, not all the merging clusters show the presence of giant diffuse radio sources and this provides important information concerning current models. The Abell 781 complex is a spectacular system composed of an apparent chain of clusters on the sky. Its main component is undergoing a merger and hosts peripheral emission that is classified as a candidate radio relic and a disputed radio halo. We used new LOw Frequency ARay (LOFAR) observations at 143 MHz and archival Giant Metrewave Radio Telescope (GMRT) observations at 325 and 610 MHz to study radio emission from non-thermal components in the ICM of Abell 781. Complementary information came from XMM-Newton data, which allowed us to investigate the connection with the thermal emission and its complex morphology. The origin of the peripheral emission is still uncertain. We speculate that it is related to the interaction between a head tail radio galaxy and shock. However, the current data allow us only to set an upper limit of $\mathcal{M} < 1.4$ on the Mach number of this putative shock. Instead, we successfully characterise the surface brightness and temperature jumps of a shock and two cold fronts in the main cluster component of Abell 781. Their positions suggest that the merger is involving three substructures. We do not find any evidence for a radio halo either at the centre of this system or in the other clusters of the chain. We place an upper limit to the diffuse radio emission in the main cluster of Abell 781 that is a factor of 2 below the current radio power-mass relation for giant radio haloes.


title: The origin of radio emission in broad absorption line quasars: Results from the LOFAR Two-metre Sky Survey

abstract: We present a study of the low-frequency radio properties of broad absorption line quasars (BALQSOs) from the LOFAR Two-metre Sky-Survey Data Release 1 (LDR1). The value-added LDR1 catalogue contains Pan-STARRS counterparts, which we match with the Sloan Digital Sky Survey (SDSS) DR7 and DR12 quasar catalogues. We find that BALQSOs are twice as likely to be detected at 144$\,$MHz than their non-BAL counterparts, and BALQSOs with low-ionisation species present in their spectra are three times more likely to be detected than those with only high-ionisation species. The BALQSO fraction at 144$\,$MHz is constant with increasing radio luminosity, which is inconsistent with previous results at 1.4$\,$GHz, indicating that observations at the different frequencies may be tracing different sources of radio emission. We cross-match radio sources between the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey and LDR1, which provides a bridge via the LDR1 Pan-STARRS counterparts to identify BALQSOs in SDSS. Consequently we expand the sample of BALQSOs detected in FIRST by a factor of three. The LDR1-detected BALQSOs in our sample are almost exclusively radio-quiet (\logr $\,<2$), with radio sizes at 144$\,$MHz typically less than $200\,$kpc; these radio sizes tend to be larger than those at 1.4$\,$GHz, suggesting more extended radio emission at low frequencies. We find that although the radio detection fraction increases with increasing balnicity index (BI), there is no correlation between BI and either low-frequency radio power or radio-loudness. This suggests that both radio emission and BI may be linked to the same underlying process, but are spatially distinct phenomena.


title: LoTSS/HETDEX: Disentangling star formation and AGN activity in gravitationally-lensed radio-quiet quasars

abstract: Determining the star-forming properties of radio-quiet quasars is important for understanding the co-evolution of star formation and black hole accretion. Here, we present the detection of the gravitationally-lensed radio-quiet quasars SDSS J1055+4628, SDSS J1313+5151 and SBS 1520+530 at 144 MHz that fall in the HETDEX Spring Field targeted in the LOFAR Two-metre Sky Survey (LoTSS) first full data release. We compare their radio and far-infrared luminosities relative to the radio-infrared correlation and find that their radio luminosities can be explained by star formation. The implied star formation rates derived from their radio and infrared luminosities are between 20 and 300 $\rm{M_{\odot}~yr^{-1}}$. These detections represent the first study of gravitationally lensed sources with LOFAR, opening a new frequency window for investigating the star-forming properties of high-redshift quasars at radio wavelengths. We consider the implications for future data releases and estimate that many of the objects in our parent sample will be detected during LoTSS, significantly increasing the fraction of gravitationally lensed radio-quiet quasars with radio detections.


title: LoTSS/HETDEX: Optical quasars I. Low-frequency radio properties of optically selected quasars

abstract: The radio-loud/radio-quiet (RL/RQ) dichotomy in quasars is still an open question. Although it is thought that accretion onto supermassive black holes in the centre the host galaxies of quasars is responsible for some radio continuum emission, there is still a debate as to whether star formation or active galactic nuclei (AGN) activity dominate the radio continuum luminosity. To date, radio emission in quasars has been investigated almost exclusively using high-frequency observations in which the Doppler boosting might have an important effect on the measured radio luminosity, whereas extended structures, best observed at low radio frequencies, are not affected by the Doppler enhancement. We used a sample of quasars selected by their optical spectra in conjunction with sensitive and high-resolution low-frequency radio data provided by the LOw Frequency ARray (LOFAR) as part of the LOFAR Two-Metre Sky Survey (LoTSS) to investigate their radio properties using the radio loudness parameter ($\mathcal{R} = \frac{L_{\mathrm{144-MHz}}}{L_{\mathrm{i\,band}}}$). The examination of the radio continuum emission and RL/RQ dichotomy in quasars exhibits that quasars show a wide continuum of radio properties (i.e. no clear bimodality in the distribution of $\mathcal{R}$). Radio continuum emission at low frequencies in low-luminosity quasars is consistent with being dominated by star formation. We see a significant albeit weak dependency of $\mathcal{R}$ on the source nuclear parameters. For the first time, we are able to resolve radio morphologies of a considerable number of quasars. All these crucial results highlight the impact of the deep and high-resolution low-frequency radio surveys that foreshadow the compelling science cases for the Square Kilometre Array (SKA).


title: The intergalactic magnetic field probed by a giant radio galaxy

abstract: Cosmological simulations predict that an intergalactic magnetic field (IGMF) pervades the large scale structure (LSS) of the Universe. Measuring the IGMF is important to determine its origin (i.e. primordial or otherwise). Using data from the LOFAR Two Metre Sky Survey (LoTSS), we present the Faraday rotation measure (RM) and depolarisation properties of the giant radio galaxy J1235+5317, at a redshift of $z = 0.34$ and 3.38 Mpc in size. We find a mean RM difference between the lobes of $2.5\pm0.1$ rad/m$^2$ , in addition to small scale RM variations of ~0.1 rad/m$^2$ . From a catalogue of LSS filaments based on optical spectroscopic observations in the local universe, we find an excess of filaments intersecting the line of sight to only one of the lobes. Associating the entire RM difference to these LSS filaments leads to a gas density-weighted IGMF strength of ~0.3 {\mu}G. However, direct comparison with cosmological simulations of the RM contribution from LSS filaments gives a low probability (~5%) for an RM contribution as large as 2.5 rad/m$^2$ , for the case of IGMF strengths of 10 to 50 nG. It is likely that variations in the RM from the Milky Way (on 11' scales) contribute significantly to the mean RM difference, and a denser RM grid is required to better constrain this contribution. In general, this work demonstrates the potential of the LOFAR telescope to probe the weak signature of the IGMF. Future studies, with thousands of sources with high accuracy RMs from LoTSS, will enable more stringent constraints on the nature of the IGMF.


title: The MAGNUM survey: different gas properties in the outflowing and disk components in nearby active galaxies with MUSE

abstract: We investigated the interstellar medium properties of the disc and outflowing gas in the central regions of nine nearby Seyfert galaxies, all characterised by prominent outflows. These objects are part of the Measuring Active Galactic Nuclei Under MUSE Microscope survey, which aims to probe their physical conditions and ionisation mechanism by exploiting MUSE unprecedented sensitivity. We studied the different properties of the gas in the disc and outflow with spatially and kinematically resolved maps by dividing the strongest emission lines in velocity bins. We associated the core of the lines with the disc, consistent with the stellar velocity, and the redshifted and the blueshifted wings with the outflow. We find that the outflowing gas is characterised by higher values of density and ionisation parameter than the disc, which presents a higher dust extinction. Moreover, we distinguish high- and low-ionisation regions across the portion of spatially resolved narrow-line region traced by the outflowing gas. The high-ionisation regions characterised by the lowest [NII]/H{\alpha} and [SII]/H{\alpha} line ratios generally trace the innermost parts along the axis of the emitting cones where the [SIII]/[SII] line ratio is enhanced, while the low-ionisation regions follow the cone edges and/or the regions perpendicular to the axis of the outflows, also characterised by a higher [OIII] velocity dispersion. A possible scenario to explain these features relies on the presence of two distinct populations of line emitting clouds: one is optically thin to the radiation and is characterised by the highest excitation, while the other is optically thick and is impinged by a filtered, and thus harder, radiation field which generates strong low-excitation lines. The highest values of [NII]/H{\alpha} and [SII]/H{\alpha} line ratios may be due to shocks and/or a hard filtered radiation field from the AGN.


title: Influence of lattice dynamics on lithium-ion conductivity: A first-principles study

abstract: In the context of novel solid electrolytes for solid-state batteries, first-principles calculations are becoming increasingly more popular due to their ability to reproduce and predict accurately the energy, structural, and dynamical properties of fast-ion conductors. In order to accelerate the discovery of new superionic conductors is convenient to establish meaningful relations between ionic transport and simple materials descriptors. Recently, several experimental studies on lithium fast-ion conductors have suggested a correlation between lattice softness and enhanced ionic conductivity due to a concomitant decrease in the activation energy for ion migration, $E_{a}$. In this article, we employ extensive \emph{ab initio} molecular dynamics simulations based on density functional theory to substantiate the links between ionic transport and lattice dynamics in a number of structurally and chemically distinct lithium superionic conductors. Our first-principles results show no evidence for a direct and general correlation between $E_{a}$, or the hopping attempt frequency, and lattice softness. However, we find that, in agreement with recent observations, the pre-exponential factor of lithium diffusivity, $D_{0}$, follows the Meyer-Neldel rule $\propto \exp{\left(E_{a}/\langle \omega \rangle\right)}$, where $\langle \omega \rangle$ represents an average phonon frequency. Hence, lattice softness can be identified with enhanced lithium diffusivity but only within families of superionic materials presenting very similar migration activation energies, due to larger $D_{0}$. On the technical side, we show that neglection of temperature effects in first-principles estimation of $E_{a}$ may lead to huge inaccuracies of $\sim 10$\%. The limitations of zero-temperature harmonic approaches in modeling of lithium-ion conductors are also illustrated.


title: The homotopy of C-motivic modular forms

abstract: A C-motivic modular forms spectrum mmf has recently been constructed. This article presents detailed computational information on the Adams spectral sequence for mmf. This information is essential for computing with the C-motivic and classical Adams spectral sequences that compute the C-motivic and classical stable homotopy groups of spheres.


title: The ATESP 5 GHz radio survey IV. 19, 38, and 94 GHz observations and radio spectral energy distributions

abstract: It is now established that the faint radio population is a mixture of star-forming galaxies and faint active galactic nuclei (AGNs), with the former dominating below S(1.4GHz) \sim 100 muJy and the latter at larger flux densities. The faint radio AGN component can itself be separated into two main classes, mainly based on the host-galaxy properties: sources associated with red/early-type galaxies (like radio galaxies) are the dominant class down to \sim 100 muJy; quasar/Seyfert--like sources contribute an additional 10-20\%. One of the major open questions regarding faint radio AGNs is the physical process responsible for their radio emission. This work aims at investigating this issue, with particular respect to the AGN component associated with red/early-type galaxies. Such AGNs show, on average, flatter radio spectra than radio galaxies and are mostly compact (<= 30 kpc in size). Various scenarios have been proposed to explain their radio emission. For instance they could be core/core-jet dominated radio galaxies, low-power BL LACs, or advection-dominated accretion flow (ADAF) systems. We used the Australia Telescope Compact Array (ATCA) to extend a previous follow-up multi-frequency campaign to 38 and 94 GHz. (abridged)


title: EFSIS: Ensemble Feature Selection Integrating Stability

abstract: Ensemble learning that can be used to combine the predictions from multiple learners has been widely applied in pattern recognition, and has been reported to be more robust and accurate than the individual learners. This ensemble logic has recently also been more applied in feature selection. There are basically two strategies for ensemble feature selection, namely data perturbation and function perturbation. Data perturbation performs feature selection on data subsets sampled from the original dataset and then selects the features consistently ranked highly across those data subsets. This has been found to improve both the stability of the selector and the prediction accuracy for a classifier. Function perturbation frees the user from having to decide on the most appropriate selector for any given situation and works by aggregating multiple selectors. This has been found to maintain or improve classification performance. Here we propose a framework, EFSIS, combining these two strategies. Empirical results indicate that EFSIS gives both high prediction accuracy and stability.


title: Do Spectroscopic Dense Gas Fractions Track Molecular Cloud Surface Densities?

abstract: We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.


title: Theory of coherent phase modes in insulating Josephson junction arrays

abstract: Recent microwave reflection measurements of Josephson junction ladders have suggested the presence of nearly coherent collective charge oscillations deep in the insulating phase. Here we develop a qualitative understanding of such coherent charge modes by studying the local dynamical conductivity of the insulating phase of a finite length sine-Gordon model. By considering parameters near the non-interacting Fermion limit where the charge operator dominantly couples to soliton-antisoliton pairs of the sine-Gordon model, we find that the local dynamical resistance shows an array of sharp peaks in frequency representing coherent phase oscillations on top of an incoherent background. The strength of the coherent peaks relative to the incoherent background increases as a powerlaw in frequency as well as exponentially as the Luttinger parameter approaches a critical value. The dynamical conductivity also clearly shows the insulating gap. We then compare the results in the high frequency limit to a perturbative estimate of phase-slip-induced decay of plasmons in the Josephson junction ladder.


title: LOFAR observations of the XMM-LSS field

abstract: We present observations of the XMM Large-Scale Structure (XMM-LSS) field observed with the LOw Frequency ARray (LOFAR) at 120-168 MHz. Centred at a J2000 declination of $-4.5^{\circ}$, this is a challenging field to observe with LOFAR because of its low elevation with respect to the array. The low elevation of this field reduces the effective collecting area of the telescope, thereby reducing sensitivity. This low elevation also causes the primary beam to be elongated in the north-south direction, which can introduce side lobes in the synthesised beam in this direction. However the XMM-LSS field is a key field to study because of the wealth of ancillary information, encompassing most of the electromagnetic spectrum. The field was observed for a total of 12 hours from three four-hour LOFAR tracks using the Dutch array. The final image presented encompasses $\sim 27$ deg$^2$, which is the region of the observations with a $>$50\% primary beam response. Once combined, the observations reach a central rms of 280 $\mu$Jy beam$^{-1}$ at 144 MHz and have an angular resolution of $7.5 \times \ 8.5$". We present our catalogue of detected sources and investigate how our observations compare to previous radio observations. This includes investigating the flux scale calibration of these observations compared to previous measurements, the implied spectral indices of the sources, the observed source counts and corrections to obtain the true source counts, and finally the clustering of the observed radio sources.


title: Radio-loud AGN in the first LoTSS data release: The lifetimes and environmental impact of jet-driven sources

abstract: We constructed a sample of 23,344 radio-loud active galactic nuclei (RLAGN) from the catalogue derived from the LOFAR Two-Metre Sky Survey (LoTSS) survey of the HETDEX Spring field. Although separating AGN from star-forming galaxies remains challenging, the combination of spectroscopic and photometric techniques we used gives us one of the largest available samples of candidate RLAGN. We used the sample, combined with recently developed analytical models, to investigate the lifetime distribution of RLAGN. We show that large or giant powerful RLAGN are probably the old tail of the general RLAGN population, but that the low-luminosity RLAGN candidates in our sample, many of which have sizes $<100$ kpc, either require a very different lifetime distribution or have different jet physics from the more powerful objects. We then used analytical models to develop a method of estimating jet kinetic powers for our candidate objects and constructed a jet kinetic luminosity function based on these estimates. These values can be compared to observational quantities, such as the integrated radiative luminosity of groups and clusters, and to the predictions from models of RLAGN feedback in galaxy formation and evolution. In particular, we show that RLAGN in the local Universe are able to supply all the energy required per comoving unit volume to counterbalance X-ray radiative losses from groups and clusters and thus prevent the hot gas from cooling. Our computation of the kinetic luminosity density of local RLAGN is in good agreement with other recent observational estimates and with models of galaxy formation.


title: Close-up view of an ongoing merger between the NGC 4839 group and the Coma cluster - a post-merger scenario

abstract: We study a merger of the NGC 4839 group with the Coma cluster using X-ray observations from the XMM-Newton and Chandra telescopes. X-ray data show two prominent features: (i) a long (~600 kpc in projection) and bent tail of cool gas trailing (towards south-west) the optical center of NGC 4839, and ii) a "sheath" region of enhanced X-ray surface brightness enveloping the group, which is due to hotter gas. While at first glance the X-ray images suggest that we are witnessing the first infall of NGC 4839 into the Coma cluster core, we argue that a post-merger scenario provides a better explanation of the observed features and illustrate this with a series of numerical simulations. In this scenario, the tail is formed when the group, initially moving to the south-west, reverses its radial velocity after crossing the apocenter, the ram pressure ceases and the ram-pressure-displaced gas falls back toward the center of the group and overshoots it. Shortly after the apocenter passage, the optical galaxy, dark matter and gaseous core move in a north-east direction, while the displaced gas continues moving to the south-west. The "sheath" is explained as being due to interaction of the re-infalling group with its own tail of stripped gas mixed with the Coma gas. In this scenario, the shock, driven by the group before reaching the apocenter, has already detached from the group and would be located close to the famous relic to the south-west of the Coma cluster.


title: Learning to synthesize: splitting and recombining low and high spatial frequencies for image recovery

abstract: Deep Neural Network (DNN)-based image reconstruction, despite many successes, often exhibits uneven fidelity between high and low spatial frequency bands. In this paper we propose the Learning Synthesis by DNN (LS-DNN) approach where two DNNs process the low and high spatial frequencies, respectively, and, improving over [30], the two DNNs are trained separately and a third DNN combines them into an image with high fidelity at all bands. We demonstrate LS-DNN in two canonical inverse problems: super-resolution (SR) in diffraction-limited imaging (DLI), and quantitative phase retrieval (QPR). Our results also show comparable or improved performance over perceptual-loss based SR [21], and can be generalized to a wider range of image recovery problems.


title: An unusual continuous paramagnetic-limited superconducting phase transition in 2D NbSe$_2$

abstract: Time reversal and spatial inversion are two key symmetries for conventional Bardeen-Cooper-Schrieffer (BCS) superconductivity. Breaking inversion symmetry can lead to mixed-parity Cooper pairing and unconventional superconducting properties. Two-dimensional (2D) NbSe$_2$ has emerged as a new non-centrosymmetric superconductor with the unique out-of-plane or Ising spin-orbit coupling (SOC). Here, we report the observation of an unusual continuous paramagnetic-limited superconductor-normal metal transition in 2D NbSe$_2$. Using tunneling spectroscopy under high in-plane magnetic fields, we observe a continuous closing of the superconducting gap at the upper critical field at low temperatures, in stark contrast to the abrupt first-order transition observed in BCS thin film superconductors. The paramagnetic-limited continuous transition arises from a large spin susceptibility of the superconducting phase due to the Ising SOC. The result is further supported by self-consistent mean-field calculations based on the ab initio band structure of 2D NbSe$_2$. Our findings establish 2D NbSe$_2$ as a promising platform for exploring novel spin-dependent superconducting phenomena and device concepts, such as equal-spin Andreev reflection and topological superconductivity.


title: On nontrivially knotted closed trajectories and positive lower bound for the energy

abstract: In this paper we suggest a soluiton to this question Does the existence of any nontrivially knotted closed trajectory for a field presents a positive lower bound for the energy and hence blocks the relaxation of the field into arbitrarily small energies


title: Photovoltaic effect in multi-domain ferroelectric perovskite oxides

abstract: We propose a device model that elucidates the role of domain walls in the photovoltaic effect in multi-domain ferroelectric perovskites. The model accounts for the intricate interplay between ferroelectric polarization, space charges, photo-generation and electronic transport. When applied to bismuth ferrite, results show a significant electric potential step across both 71-degree and 109-degree domain walls, which in turn contributes to the photovoltaic (PV) effect. We also find a strong correlation between polarization and oxygen octahedra tilts, which indicates the nontrivial role of the latter in the PV effect. The domain wall-based PV effect is further shown to be additive in nature, allowing for the possibility of generating above-bandgap voltage


title: The environments of radio-loud AGN from the LOFAR Two-Metre Sky Survey (LoTSS)

abstract: An understanding of the relationship between radio-loud active galaxies and their large-scale environments is essential for realistic modelling of radio-galaxy evolution and environmental impact, for understanding AGN triggering and life cycles, and for calibrating galaxy feedback in cosmological models. We use the LOFAR Two-Metre Sky Survey (LoTSS) Data Release 1 catalogues to investigate this relationship. We cross-matched 8,745 radio-loud AGN with 0.08<z<0.4, selected from LoTSS, with two Sloan Digital Sky Survey cluster catalogues, and find that only 10 percent of LoTSS AGN in this redshift range have an association, so that the majority of low-redshift AGN (including a substantial fraction of the most radio-luminous objects) must inhabit haloes with M < 10^14 M_sun. We find that the probability of a cluster association, and the richness of the associated cluster, is correlated with AGN radio luminosity, and for the cluster population, the number of associated AGN and the radio luminosity of the brightest associated AGN is richness-dependent. We demonstrate that these relations are not driven solely by host-galaxy stellar mass, supporting models in which large-scale environment is influential in driving AGN jet activity. At the lowest radio luminosities we find that the minority of objects with a cluster association are located at larger mean cluster-centre distances than more luminous AGN, which appears to be driven primarily by host-galaxy mass. Finally, we also find that FRI radio galaxies inhabit systematically richer environments than FRIIs. The work presented here demonstrates the potential of LoTSS for AGN environmental studies. In future, the full northern-sky LoTSS catalogue, together with the use of deeper optical/IR imaging data and spectroscopic follow-up with WEAVE-LOFAR, will provide opportunities to extend this type of work to much larger samples and higher redshifts.


title: Optimal Transport Classifier: Defending Against Adversarial Attacks by Regularized Deep Embedding

abstract: Recent studies have demonstrated the vulnerability of deep convolutional neural networks against adversarial examples. Inspired by the observation that the intrinsic dimension of image data is much smaller than its pixel space dimension and the vulnerability of neural networks grows with the input dimension, we propose to embed high-dimensional input images into a low-dimensional space to perform classification. However, arbitrarily projecting the input images to a low-dimensional space without regularization will not improve the robustness of deep neural networks. Leveraging optimal transport theory, we propose a new framework, Optimal Transport Classifier (OT-Classifier), and derive an objective that minimizes the discrepancy between the distribution of the true label and the distribution of the OT-Classifier output. Experimental results on several benchmark datasets show that, our proposed framework achieves state-of-the-art performance against strong adversarial attack methods.


title: On the maximal multiplicity of block sizes in a random set partition

abstract: We study the asymptotic behavior of the maximal multiplicity $M_n=M_n(\sigma)$ of the blocks in a set partition of $[n]=\{1,2,...,n\}$, assuming that $\sigma$ is chosen uniformly at random from the set of all such partitions. Let $W=W(n)$ be the unique positive root of the equation $We^W=n$ and let $f_n$ be the fractional part of $W(n)$. Furthermore, let $R_n=W^{\lfloor W\rfloor}/\lfloor W\rfloor !$ and let $\vartheta_n=\min{\{f_n,1-f_n\}}$. We show that, over a subsequence $\{n_k\}_{k\ge 1}$, $(M_{n_k}-R_{n_k})/\sqrt{R_{n_k}}$ converges weakly, as $k\to\infty$, to $\max{\{Z_1,Z_2-u\}}$, where $Z_1$ and $Z_2$ are two independent copies of a standard normal random variable and either $u=\left(\frac{1}{2\pi}\right)^{1/4}\lim_{k\to\infty}\vartheta_{n_k}\frac{\sqrt{n_k}}{\log^{7/4}{n_k}}\in [0,\infty)$ or $u=\infty$. The proof uses the saddle point method. A comparison with the similar statistic for random integer partitions of $n$ is also given.


title: Diagnostics of a Nuclear Starburst: Water and Methanol Masers

abstract: We test models of starburst driven outflows using observations of the 22.2 GHz H$_2$O and 36.2 GHz class I CH$_3$OH maser lines. We have observed the starburst galaxy NGC 253 using the Karl G. Jansky Very Large Array. We present evidence for entrainment of star-forming dense-molecular gas in the outflow of NGC 253. We also show that H$_2$O masers are associated with forming super star clusters and not with supernova remnants. We detect four new 36 GHz CH$_3$OH masers in the central kpc and show possible evidence for a star-formation origin of two class I CH$_3$OH masers. Such high resolution observations are essential for understanding the origin of these masers.


title: Atomki Anomaly in Family-Dependent $U(1)'$ Extension of the Standard Model

abstract: In the context of a gauge invariant, non-anomalous and family-dependent (non-universal) $U(1)'$ extension of the Standard Model, wherein a new high scale mechanism generates masses and couplings for the first two fermion generations and the standard Higgs mechanism does so for the third one, we find solutions to the anomaly observed by the Atomki collaboration in the decay of excited states of Beryllium, in the form of a very light $Z'$ state, stemming from the $U(1)'$ symmetry breaking, with significant axial couplings so as to evade a variety of low scale experimental constraints.


title: Systematic effects in LOFAR data: A unified calibration strategy

abstract: Context: New generation low-frequency telescopes are exploring a new parameter space in terms of depth and resolution. The data taken with these interferometers, for example with the LOw Frequency ARray (LOFAR), are often calibrated in a low signal-to-noise ratio regime and the removal of critical systematic effects is challenging. The process requires an understanding of their origin and properties. Aim: In this paper we describe the major systematic effects inherent to next generation low-frequency telescopes, such as LOFAR. With this knowledge, we introduce a data processing pipeline that is able to isolate and correct these systematic effects. The pipeline will be used to calibrate calibrator observations as the first step of a full data reduction process. Methods: We processed two LOFAR observations of the calibrator 3C196: the first using the Low Band Antenna (LBA) system at 42-66 MHz and the second using the High Band Antenna (HBA) system at 115-189 MHz. Results: We were able to isolate and correct for the effects of clock drift, polarisation misalignment, ionospheric delay, Faraday rotation, ionospheric scintillation, beam shape, and bandpass. The designed calibration strategy produced the deepest image to date at 54 MHz. The image has been used to confirm that the spectral energy distribution of the average radio source population tends to flatten at low frequencies. Conclusions: We prove that LOFAR systematic effects can be described by a relatively small number of parameters. Furthermore, the identification of these parameters is fundamental to reducing the degrees of freedom when the calibration is carried out on fields that are not dominated by a strong calibrator.


title: Asymptotics of Operator Semigroups via the Semigroup at Infinity

abstract: We systematize and generalize recent results of Gerlach and Gl\"uck on the strong convergence and spectral theory of bounded (positive) operator semigroups $(T_s)_{s\in S}$ on Banach spaces (lattices). (Here, $S$ can be an arbitrary commutative semigroup, and no topological assumptions neither on $S$ nor on its representation are required.) To this aim, we introduce the "semigroup at infinity" and give useful criteria ensuring that the well-known Jacobs--de Leeuw--Glicksberg splitting theory can be applied to it. Next, we confine these abstract results to positive semigroups on Banach lattices with a quasi-interior point. In that situation, the said criteria are intimately linked to so-called AM-compact operators (which entail kernel operators and compact operators); and they imply that the original semigroup asymptotically embeds into a compact group of positive invertible operators on an atomic Banach lattice. By means of a structure theorem for such group representations (reminiscent of the Peter--Weyl theorem and its consequences for Banach space representations of compact groups) we are able to establish quite general conditions implying the strong convergence of the original semigroup. Finally, we show how some classical results of Greiner (1982), Davies (2005), Keicher (2006) and Arendt (2008) and more recent ones by Gerlach and Gl\"uck (2017) are covered and extended through our approach.


title: The First Swift Intensive AGN Accretion Disk Reverberation Mapping Survey

abstract: Swift intensive accretion disk reverberation mapping of four AGN yielded light curves sampled $\sim$200-350 times in 0.3-10 keV X-ray and six UV/optical bands. Uniform reduction and cross-correlation analysis of these datasets yields three main results: 1) The X-ray/UV correlations are much weaker than those within the UV/optical, posing severe problems for the lamp-post reprocessing model in which variations in a central X-ray corona drive and power those in the surrounding accretion disk. 2) The UV/optical interband lags are generally consistent with $ \tau \propto \lambda^{4/3} $ as predicted by the centrally illuminated thin accretion disk model. While the average interband lags are somewhat larger than predicted, these results alone are not inconsistent with the thin disk model given the large systematic uncertainties involved. 3) The one exception is the U band lags, which are on average a factor of $\sim$2.2 larger than predicted from the surrounding band data and fits. This excess appears due to diffuse continuum emission from the broad-line region (BLR). The precise mixing of disk and BLR components cannot be determined from these data alone. The lags in different AGN appear to scale with mass or luminosity. We also find that there are systematic differences between the uncertainties derived by javelin vs. more standard lag measurement techniques, with javelin reporting smaller uncertainties by a factor of 2.5 on average. In order to be conservative only standard techniques were used in the analyses reported herein.


title: Model change detection with application to machine learning

abstract: Model change detection is studied, in which there are two sets of samples that are independently and identically distributed (i.i.d.) according to a pre-change probabilistic model with parameter $\theta$, and a post-change model with parameter $\theta'$, respectively. The goal is to detect whether the change in the model is significant, i.e., whether the difference between the pre-change parameter and the post-change parameter $\|\theta-\theta'\|_2$ is larger than a pre-determined threshold $\rho$. The problem is considered in a Neyman-Pearson setting, where the goal is to maximize the probability of detection under a false alarm constraint. Since the generalized likelihood ratio test (GLRT) is difficult to compute in this problem, we construct an empirical difference test (EDT), which approximates the GLRT and has low computational complexity. Moreover, we provide an approximation method to set the threshold of the EDT to meet the false alarm constraint. Experiments with linear regression and logistic regression are conducted to validate the proposed algorithms.


title: Tukey-Inspired Video Object Segmentation

abstract: We investigate the problem of strictly unsupervised video object segmentation, i.e., the separation of a primary object from background in video without a user-provided object mask or any training on an annotated dataset. We find foreground objects in low-level vision data using a John Tukey-inspired measure of "outlierness". This Tukey-inspired measure also estimates the reliability of each data source as video characteristics change (e.g., a camera starts moving). The proposed method achieves state-of-the-art results for strictly unsupervised video object segmentation on the challenging DAVIS dataset. Finally, we use a variant of the Tukey-inspired measure to combine the output of multiple segmentation methods, including those using supervision during training, runtime, or both. This collectively more robust method of segmentation improves the Jaccard measure of its constituent methods by as much as 28%.


title: The Aharoni--Korman conjecture for $N$-free posets with no infinite antichain

abstract: We give a necessary and sufficient condition for a $P_4$-free graph to be a cograph. This allows us to obtain a simple proof of the fact that finite $P_4$-free graphs are finite cographs. We also prove that $N$-free chain complete posets and $N$-free posets with no infinite antichains are series-parallel. As a consequence, we obtain that every $N$-free poset with no infinite antichain has a chain and a partition into antichains so that each part intersects the chain. This answers a conjecture of Aharoni and Korman (Order \textbf{9} (1992) 245--253) in this case.


title: The slice spectral sequence of a $C_4$-equivariant height-4 Lubin-Tate theory

abstract: We completely compute the slice spectral sequence of the $C_4$-spectrum $BP^{((C_4))}\langle 2 \rangle$. After periodization and $K(4)$-localization, this spectrum is equivalent to a height-4 Lubin-Tate theory $E_4$ with $C_4$-action induced from the Goerss-Hopkins-Miller theorem. In particular, our computation shows that $E_4^{hC_{12}}$ is 384-periodic.


title: Blazars in the LOFAR Two-Metre Sky Survey First Data Release

abstract: Historically, the blazar population has been poorly understood at low frequencies because survey sensitivity and angular resolution limitations have made it difficult to identify megahertz counterparts. We used the LOFAR Two-Metre Sky Survey (LoTSS) first data release value-added catalogue (LDR1) to study blazars in the low-frequency regime with unprecedented sensitivity and resolution. We identified radio counterparts to all $98$ known sources from the Third \textit{Fermi}-LAT Point Source Catalogue (3FGL) or Roma-BZCAT Multi-frequency Catalogue of Blazars ($5^{\mathrm{th}}$ edition) that fall within the LDR1 footprint. Only the 3FGL unidentified $\gamma$-ray sources (UGS) could not be firmly associated with an LDR1 source; this was due to source confusion. We examined the redshift and radio luminosity distributions of our sample, finding flat-spectrum radio quasars (FSRQs) to be more distant and more luminous than BL Lacertae objects (BL Lacs) on average. Blazars are known to have flat spectra in the gigahertz regime but we found this to extend down to $144$ MHz, where the radio spectral index, $\alpha$, of our sample is $-0.17 \pm 0.14$. For BL Lacs, $\alpha = -0.13 \pm 0.16$ and for FSRQs, $\alpha = -0.15 \pm 0.17$. We also investigated the radio-to-$\gamma$-ray connection for the $30$ $\gamma$-ray-detected sources in our sample. We find Pearson's correlation coefficient is $0.45$ ($p = 0.069$). This tentative correlation and the flatness of the spectral index suggest that the beamed core emission contributes to the low-frequency flux density. We compare our sample distribution with that of the full LDR1 on colour-colour diagrams, and we use this information to identify possible radio counterparts to two of the four UGS within the LDR1 field. We will refine our results as LoTSS continues.


title: A conditional space-time POD formalism for intermittent and rare events: example of acoustic bursts in turbulent jets

abstract: We present a conditional space-time proper orthogonal decomposition (POD) formulation that is tailored to the eduction of the average, rare or intermittent event from an ensemble of realizations of a fluid process. By construction, the resulting spatio-temporal modes are coherent in space and over a pre-defined finite time horizon and optimally capture the variance, or energy of the ensemble. For the example of intermittent acoustic radiation from a turbulent jet, we introduce a conditional expectation operator that focuses on the loudest events, as measured by a pressure probe in the far-field and contained in the tail of the pressure signal's probability distribution. Applied to high-fidelity simulation data, the method identifies a statistically significant `prototype', or average acoustic burst event that is tracked over time. Most notably, the burst event can be traced back to its precursor, which opens up the possibility of prediction of an imminent burst. We furthermore investigate the mechanism underlying the prototypical burst event using linear stability theory and find that its structure and evolution is accurately predicted by optimal transient growth theory. The jet-noise problem demonstrates that the conditional space-time POD formulation applies even for systems with probability distributions that are not heavy-tailed, i.e. for systems in which events overlap and occur in rapid succession.


title: Long range dipole-dipole interaction in atomic vapors probed by double-quantum two-dimensional coherent spectroscopy

abstract: Optical double-quantum two-dimensional coherent spectroscopy (2DCS) was implemented to probe interatomic dipole-dipole interactions in both potassium and rubidium atomic vapors. The dipole-dipole interaction was detected at densities of $4.81 \times 10^8$ cm$^{-3}$ and $8.40 \times 10^9$ cm$^{-3}$ for potassium and rubidium, respectively, corresponding to a mean interatomic separation of 15.8 $\mu$m or $3.0\times 10^5a_0$ for potassium and 6.1 $\mu$m or $1.2\times 10^5a_0$ for rubidium, where $a_0$ is the Bohr radius. We report the lowest atomic density at which dipole-dipole interactions are detected. The experimental results confirm the long range nature of the dipole-dipole interaction which is critical for understanding many-body physics in atoms/molecules. The long range interaction also has implications in atom-based applications involving many-body interactions. Additionally, we demonstrated that double-quantum 2DCS is sufficiently sensitive to probe dipole-dipole interaction at densities that can be achieved with cold atom in a magneto-optical trap, paving the way for double-quantum 2DCS studies of cold atoms and molecules. The method can also open a new avenue to study long-range interactions in solid states systems such as quantum dots and color centers in diamonds.


title: Chip-based Resonance Raman Spectroscopy Using Tantalum Pentoxide Waveguides

abstract: Blood analysis is an important diagnostic tool, as it provides a wealth of information about the patient health. Raman spectroscopy is a promising tool in blood analysis, but widespread clinical application is limited by its low signal strength, as well as complex and costly instrumentation. The growing field of waveguide-based Raman spectroscopy tries to solve these challenges by working towards fully integrated Raman sensors with increased interaction areas. In this work, we demonstrate on-chip resonance Raman measurements of hemoglobin, a crucial component of blood, at 532 nm excitation using a tantalum pentoxide (Ta2O5) waveguide platform. We have also characterized the background signal from Ta2O5 waveguide material when excited at 532 nm. In addition, we demonstrate spontaneous Raman measurements of iso-propanol and methanol using the same platform. Our results suggest that Ta2O5 is a promising waveguide platform for resonance Raman spectroscopy at 532 nm, and in particular for blood analysis.


title: dS/dS and $T\bar T$

abstract: The $T\bar T$ deformation of a conformal field theory has a dual description as a cutoff $AdS_3$ spacetime, at least at the level of pure 3d gravity. We generalize this deformation in such a way that it builds up a patch of bulk $dS_3$ spacetime instead. At each step along the trajectory in the space of $2d$ theories, the theory is deformed by a specific combination of $T\bar T$ and the two-dimensional cosmological constant. This provides a concrete holographic dual for the warped throat on the gravity side of the dS/dS duality, at leading order in large central charge. We also analyze a sequence of excitations of this throat on both sides of the duality, as well as the entanglement entropy. Our results point toward a mechanism for obtaining de Sitter solutions starting from seed conformal field theories with AdS duals.


title: Mitigating Architectural Mismatch During the Evolutionary Synthesis of Deep Neural Networks

abstract: Evolutionary deep intelligence has recently shown great promise for producing small, powerful deep neural network models via the organic synthesis of increasingly efficient architectures over successive generations. Existing evolutionary synthesis processes, however, have allowed the mating of parent networks independent of architectural alignment, resulting in a mismatch of network structures. We present a preliminary study into the effects of architectural alignment during evolutionary synthesis using a gene tagging system. Surprisingly, the network architectures synthesized using the gene tagging approach resulted in slower decreases in performance accuracy and storage size; however, the resultant networks were comparable in size and performance accuracy to the non-gene tagging networks. Furthermore, we speculate that there is a noticeable decrease in network variability for networks synthesized with gene tagging, indicating that enforcing a like-with-like mating policy potentially restricts the exploration of the search space of possible network architectures.


title: Hypergeometric function and Modular Curvature II. Connes-Moscovici functional relation after Lesch's work

abstract: As the second part of the sequel, we investigate the variation of rearrangement operators (more precisely, the spectral functions behind) arising in the study of modular geometry on noncommutative (two) tori. We initiate a systematic approach by introducing transformations corresponding to basic operations in calculus, like differentiation and integration by parts. As for applications, we extend, in a uniform way, the Connes-Moscovici's functional relations on noncommutative two tori attached to the variation of second heat coefficients to noncommutative tori of arbitrary dimension. Moreover, those transformations lead to more internal relations among the hypergeometric family obtained in part I of the sequel, which allows us to obtain, the first time, a computer-aid free verification of those Connes-Moscovici type functional relations.


title: Calibrating the relation of low-frequency radio continuum to star formation rate at 1 kpc scale with LOFAR

abstract: Radio continuum (RC) emission in galaxies allows us to measure star formation rates (SFRs) unaffected by extinction due to dust, of which the low-frequency part is uncontaminated from thermal (free-free) emission. We calibrate the conversion from the spatially resolved 140 MHz RC emission to the SFR surface density ($\Sigma_{\rm SFR}$) at 1 kpc scale. We used recent observations of three galaxies (NGC 3184, 4736, and 5055) from the LOFAR Two-metre Sky Survey (LoTSS), and archival LOw-Frequency ARray (LOFAR) data of NGC 5194. Maps were created with the facet calibration technique and converted to radio $\Sigma_{\rm SFR}$ maps using the Condon relation. We compared these maps with hybrid $\Sigma_{\rm SFR}$ maps from a combination of GALEX far-ultraviolet and Spitzer 24 $\mu\rm m$ data using plots tracing the relation at $1.2\times 1.2$-kpc$^2$ resolution. The RC emission is smoothed with respect to the hybrid $\Sigma_{\rm SFR}$ owing to the transport of cosmic-ray electrons (CREs). This results in a sublinear relation $(\Sigma_{\rm SFR})_{\rm RC} \propto [(\Sigma_{\rm SFR})_{\rm hyb}]^{a}$, where $a=0.59\pm 0.13$ (140 MHz) and $a=0.75\pm 0.10$ (1365 MHz). Both relations have a scatter of $\sigma = 0.3~\rm dex$. If we restrict ourselves to areas of young CREs ($\alpha > -0.65$; $I_\nu \propto \nu^\alpha$), the relation becomes almost linear at both frequencies with $a\approx 0.9$ and a reduced scatter of $\sigma = 0.2~\rm dex$. We then simulate the effect of CRE transport by convolving the hybrid $\Sigma_{\rm SFR}$ maps with a Gaussian kernel until the RC-SFR relation is linearised; CRE transport lengths are $l=1$-5 kpc. Solving the CRE diffusion equation, we find diffusion coefficients of $D=(0.13$-$1.5) \times 10^{28} \rm cm^2\,s^{-1}$ at 1 GeV. A RC-SFR relation at $1.4$ GHz can be exploited to measure SFRs at redshift $z \approx 10$ using $140$ MHz observations.


title: Informed MCMC with Bayesian Neural Networks for Facial Image Analysis

abstract: Computer vision tasks are difficult because of the large variability in the data that is induced by changes in light, background, partial occlusion as well as the varying pose, texture, and shape of objects. Generative approaches to computer vision allow us to overcome this difficulty by explicitly modeling the physical image formation process. Using generative object models, the analysis of an observed image is performed via Bayesian inference of the posterior distribution. This conceptually simple approach tends to fail in practice because of several difficulties stemming from sampling the posterior distribution: high-dimensionality and multi-modality of the posterior distribution as well as expensive simulation of the rendering process. The main difficulty of sampling approaches in a computer vision context is choosing the proposal distribution accurately so that maxima of the posterior are explored early and the algorithm quickly converges to a valid image interpretation. In this work, we propose to use a Bayesian Neural Network for estimating an image dependent proposal distribution. Compared to a standard Gaussian random walk proposal, this accelerates the sampler in finding regions of the posterior with high value. In this way, we can significantly reduce the number of samples needed to perform facial image analysis.


title: Topological Sieving of Rings According to their Rigidity

abstract: We present a novel mechanism for resolving the mechanical rigidity of nanoscopic circular polymers that flow in a complex environment. The emergence of a regime of negative differential mobility induced by topological interactions between the rings and the substrate is the key mechanism for selective sieving of circular polymers with distinct flexibilities. A simple model accurately describes the sieving process observed in molecular dynamics simulations and yields experimentally verifiable analytical predictions, which can be used as a reference guide for improving filtration procedures of circular filaments. The topological sieving mechanism we propose ought to be relevant also in probing the microscopic details of complex substrates.


title: Private Selection from Private Candidates

abstract: Differentially Private algorithms often need to select the best amongst many candidate options. Classical works on this selection problem require that the candidates' goodness, measured as a real-valued score function, does not change by much when one person's data changes. In many applications such as hyperparameter optimization, this stability assumption is much too strong. In this work, we consider the selection problem under a much weaker stability assumption on the candidates, namely that the score functions are differentially private. Under this assumption, we present algorithms that are near-optimal along the three relevant dimensions: privacy, utility and computational efficiency. Our result can be seen as a generalization of the exponential mechanism and its existing generalizations. We also develop an online version of our algorithm, that can be seen as a generalization of the sparse vector technique to this weaker stability assumption. We show how our results imply better algorithms for hyperparameter selection in differentially private machine learning, as well as for adaptive data analysis.


title: Dynamics of Dirac observables in canonical cosmological perturbation theory

abstract: The relational formalism based on geometrical clocks and Dirac observables in linearized canonical cosmological perturbation theory is used to introduce an efficient method to find evolution equations for gauge invariant variables. Our method generalizes an existing technique by Pons, Salisbury and Sundermeyer [1, 2] to relate the evolution of gauge invariant observables with the one of gauge variant quantities, and is applied as a demonstration for the longitudinal and spatially flat gauges. Gauge invariant evolution equations for the Bardeen potential and the Mukhanov-Sasaki variable are derived in the extended ADM phase space. Our method establishes a full agreement at the dynamical level between the canonical and conventional cosmological perturbation theory at the linear order using Dirac observables.


title: Ultrafast Structural Dynamics of Photo-Reactions Revealed by Model-Independent X-ray Cross-Correlation Analysis

abstract: We applied angular X-ray Cross-Correlation analysis (XCCA) to scattering images from a femtosecond resolution LCLS X-ray free-electron laser (XFEL) pump-probe experiment with solvated PtPOP ([Pt$_2$(P$_2$O$_5$H$_2$)$_4$]$^{4-}$) metal complex molecules. The molecules were pumped with linear polarized laser pulses creating an excited state population with a preferred orientational (alignment) direction. Two time scales of $1.9\pm1.5$ ps and $46\pm10$ ps were revealed by model-independent XCCA, associated with an internal structural changes and rotational dephasing, respectively. Our studies illustrate the potential of XCCA to reveal hidden structural information in a model independent analysis of time evolution of solvated metal complex molecules.


title: A Map of Knowledge

abstract: Knowledge representation has gained in relevance as data from the ubiquitous digitization of behaviors amass and academia and industry seek methods to understand and reason about the information they encode. Success in this pursuit has emerged with data from natural language, where skip-grams and other linear connectionist models of distributed representation have surfaced scrutable relational structures which have also served as artifacts of anthropological interest. Natural language is, however, only a fraction of the big data deluge. Here we show that latent semantic structure, comprised of elements from digital records of our interactions, can be informed by behavioral data and that domain knowledge can be extracted from this structure through visualization and a novel mapping of the literal descriptions of elements onto this behaviorally informed representation. We use the course enrollment behaviors of 124,000 students at a public university to learn vector representations of its courses. From these behaviorally informed representations, a notable 88% of course attribute information were recovered (e.g., department and division), as well as 40% of course relationships constructed from prior domain knowledge and evaluated by analogy (e.g., Math 1B is to Math H1B as Physics 7B is to Physics H7B). To aid in interpretation of the learned structure, we create a semantic interpolation, translating course vectors to a bag-of-words of their respective catalog descriptions. We find that the representations learned from enrollments resolved course vectors to a level of semantic fidelity exceeding that of their catalog descriptions, depicting a vector space of high conceptual rationality. We end with a discussion of the possible mechanisms by which this knowledge structure may be informed and its implications for data science.


title: On polycyclic codes over a finite chain ring

abstract: Galois images of polycyclic codes over a finite chain ring $S$ and their annihilator dual are investigated. The case when a polycyclic codes is Galois-disjoint over the ring $S,$ is characterized and, the trace codes and restrictions of free polycyclic codes over $S$ are also determined givind an analogue of Delsarte theorem among trace map, any S -linear code and its annihilator dual.


title: D\'ecomposition au-dessus des param\`etres de Langlands elliptiques

abstract: We prove that, over any elliptic global Langlands parameter $\sigma$, the cuspidal cohomology groups of moduli stacks of shtukas are given by a formula involving a finite dimensional representation of the centralizer of $\sigma$. It is a first step in the direction of Arthur-Kottwitz conjectures.


title: ShapeSearch: A Flexible and Efficient System for Shape-based Exploration of Trendlines

abstract: Identifying trendline visualizations with desired patterns is a common and fundamental data exploration task. Existing visual analytics tools offer limited flexibility and expressiveness for such tasks, especially when the pattern of interest is under-specified and approximate, and do not scale well when the pattern searching needs are ad-hoc, as is often the case. We propose ShapeSearch, an efficient and flexible pattern-searching tool, that enables the search for desired patterns via multiple mechanisms: sketch, natural-language, and visual regular expressions. We develop a novel shape querying algebra, with a minimal set of primitives and operators that can express a large number of ShapeSearch queries, and design a natural-language and regex-based parser to automatically parse and translate user queries to the algebra representation. To execute these queries within interactive response times, ShapeSearch uses a fast shape algebra-based execution engine with query-aware optimizations, and perceptually-aware scoring methodologies. We present a thorough evaluation of the system, including a general-purpose user study, a case study involving genomic data analysis, as well as performance experiments, comparing against state-of-the-art time series shape matching approaches---that together demonstrate the usability and scalability of ShapeSearch.