Showing 1–50 of 167 results for author: Peiris, H V

pop-cosmos: Scaleable inference of galaxy properties and redshifts with a data-driven population model

Authors: Stephen Thorp, Justin Alsing, Hiranya V. Peiris, Sinan Deger, Daniel J. Mortlock, Boris Leistedt, Joel Leja, Arthur Loureiro

Abstract: We present an efficient Bayesian method for estimating individual photometric redshifts and galaxy properties under a pre-trained population model (pop-cosmos) that was calibrated using purely photometric data. This model specifies a prior distribution over 16 stellar population synthesis (SPS) parameters using a score-based diffusion model, and includes a data model with detailed treatment of neb… ▽ More We present an efficient Bayesian method for estimating individual photometric redshifts and galaxy properties under a pre-trained population model (pop-cosmos) that was calibrated using purely photometric data. This model specifies a prior distribution over 16 stellar population synthesis (SPS) parameters using a score-based diffusion model, and includes a data model with detailed treatment of nebular emission. We use a GPU-accelerated affine invariant ensemble sampler to achieve fast posterior sampling under this model for 292,300 individual galaxies in the COSMOS2020 catalog, leveraging a neural network emulator (Speculator) to speed up the SPS calculations. We apply both the pop-cosmos population model and a baseline prior inspired by Prospector-$α$, and compare these results to published COSMOS2020 redshift estimates from the widely-used EAZY and LePhare codes. For the $\sim 12,000$ galaxies with spectroscopic redshifts, we find that pop-cosmos yields redshift estimates that have minimal bias ($\sim10^{-4}$), high accuracy ($σ_\text{MAD}=7\times10^{-3}$), and a low outlier rate ($1.6\%$). We show that the pop-cosmos population model generalizes well to galaxies fainter than its $r<25$ mag training set. The sample we have analyzed is $\gtrsim3\times$ larger than has previously been possible via posterior sampling with a full SPS model, with average throughput of 15 GPU-sec per galaxy under the pop-cosmos prior, and 0.6 GPU-sec per galaxy under the Prospector prior. This paves the way for principled modeling of the huge catalogs expected from upcoming Stage IV galaxy surveys. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 22 pages, 14 figures. For submission to ApJ

Deep learning insights into non-universality in the halo mass function

Authors: Ningyuan Guo, Luisa Lucie-Smith, Hiranya V. Peiris, Andrew Pontzen, Davide Piras

Abstract: The abundance of dark matter haloes is a key cosmological probe in forthcoming galaxy surveys. The theoretical understanding of the halo mass function (HMF) is limited by our incomplete knowledge of the origin of non-universality and its cosmological parameter dependence. We present a deep learning model which compresses the linear matter power spectrum into three independent factors which are nec… ▽ More The abundance of dark matter haloes is a key cosmological probe in forthcoming galaxy surveys. The theoretical understanding of the halo mass function (HMF) is limited by our incomplete knowledge of the origin of non-universality and its cosmological parameter dependence. We present a deep learning model which compresses the linear matter power spectrum into three independent factors which are necessary and sufficient to describe the $z=0$ HMF from the state-of-the-art AEMULUS emulator to sub-per cent accuracy in a $w$CDM$+N_\mathrm{eff}$ parameter space. Additional information about growth history does not improve the accuracy of HMF predictions if the matter power spectrum is already provided as input, because required aspects of the former can be inferred from the latter. The three factors carry information about the universal and non-universal aspects of the HMF, which we interrogate via the information-theoretic measure of mutual information. We find that non-universality is captured by recent growth history after matter-dark-energy equality and $N_\mathrm{eff}$ for $M\sim 10^{13} \, \mathrm{M_\odot}\, h^{-1}$ haloes, and by $Ω_{\rm m}$ for $M\sim 10^{15} \, \mathrm{M_\odot}\, h^{-1}$. The compact representation learnt by our model can inform the design of emulator training sets to achieve high emulator accuracy with fewer simulations. △ Less

Submitted 9 July, 2024; v1 submitted 24 May, 2024; originally announced May 2024.

Comments: 17 pages, 13 figures. Minor changes to match version accepted for publication in MNRAS

Data-Space Validation of High-Dimensional Models by Comparing Sample Quantiles

Authors: Stephen Thorp, Hiranya V. Peiris, Daniel J. Mortlock, Justin Alsing, Boris Leistedt, Sinan Deger

Abstract: We present a simple method for assessing the predictive performance of high-dimensional models directly in data space when only samples are available. Our approach is to compare the quantiles of observables predicted by a model to those of the observables themselves. In cases where the dimensionality of the observables is large (e.g. multiband galaxy photometry), we advocate that the comparison is… ▽ More We present a simple method for assessing the predictive performance of high-dimensional models directly in data space when only samples are available. Our approach is to compare the quantiles of observables predicted by a model to those of the observables themselves. In cases where the dimensionality of the observables is large (e.g. multiband galaxy photometry), we advocate that the comparison is made after projection onto a set of principal axes to reduce the dimensionality. We demonstrate our method on a series of two-dimensional examples. We then apply it to results from a state-of-the-art generative model for galaxy photometry (pop-cosmos; arXiv:2402.00935) that generates predictions of colors and magnitudes by forward simulating from a 16-dimensional distribution of physical parameters represented by a score-based diffusion model. We validate the predictive performance of this model directly in a space of nine broadband colors. Although motivated by this specific example, the techniques we present will be broadly useful for evaluating the performance of flexible, non-parametric population models of this kind, and can be readily applied to any setting where two sets of samples are to be compared. △ Less

Submitted 16 February, 2024; v1 submitted 1 February, 2024; originally announced February 2024.

Comments: 12 pages, 7 figures. Submitted to ApJS. Companion paper to arXiv:2402.00935

Detecting strongly-lensed type Ia supernovae with LSST

Authors: Nikki Arendse, Suhail Dhawan, Ana Sagués Carracedo, Hiranya V. Peiris, Ariel Goobar, Radek Wojtak, Catarina Alves, Rahul Biswas, Simon Huber, Simon Birrer, The LSST Dark Energy Science Collaboration

Abstract: Strongly-lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (… ▽ More Strongly-lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (SNIa) in LSST by quantifying the expected annual number of detections, the impact of stellar microlensing, follow-up feasibility, and how to best separate lensed and unlensed SNIa. We simulate SNIa lensed by galaxies, using the current LSST baseline v3.0 cadence, and find an expected number of 44 lensed SNIa detections per year. Microlensing effects by stars in the lens galaxy are predicted to lower the lensed SNIa detections by $\sim 8 \%$. The lensed events can be separated from the unlensed ones by jointly considering their colours and peak magnitudes. We define a `gold sample' of $\sim 10$ lensed SNIa per year with time delay $> 10$ days, $> 5$ detections before light-curve peak, and sufficiently bright ($m_i < 22.5$ mag) for follow-up observations. In three years of LSST operations, such a sample is expected to yield a $1.5\%$ measurement of the Hubble constant. △ Less

Submitted 16 May, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS; 15 pages, 13 figures

An Anti-halo Void Catalogue of the Local Super-Volume

Authors: Stephen Stopyra, Hiranya V. Peiris, Andrew Pontzen, Jens Jasche, Guilhem Lavaux

Abstract: We construct an anti-halo void catalogue of $150$ voids with radii $R > 10\,h^{-1}\mathrm{\,Mpc}$ in the Local Super-Volume ($<135\,h^{-1}\mathrm{\,Mpc}$ from the Milky Way), using posterior resimulation of initial conditions inferred by field-level inference with Bayesian Origin Reconstruction from Galaxies (\codefont{BORG}). We describe and make use of a new algorithm for creating a single, unif… ▽ More We construct an anti-halo void catalogue of $150$ voids with radii $R > 10\,h^{-1}\mathrm{\,Mpc}$ in the Local Super-Volume ($<135\,h^{-1}\mathrm{\,Mpc}$ from the Milky Way), using posterior resimulation of initial conditions inferred by field-level inference with Bayesian Origin Reconstruction from Galaxies (\codefont{BORG}). We describe and make use of a new algorithm for creating a single, unified void catalogue by combining different samples from the posterior. The catalogue is complete out to $135\,h^{-1}\mathrm{\,Mpc}$, with void abundances matching theoretical predictions. Finally, we compute stacked density profiles of those voids which are reliably identified across posterior samples, and show that these are compatible with $Λ$CDM expectations once environmental selection (e.g., the estimated $\sim 4\%$ under-density of the Local Super-Volume) is accounted for. △ Less

Submitted 16 May, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 11 pages, 4 figures. Updated to match version accepted by MNRAS. Minor clarifications and no changes to conclusion or catalogue

Measuring the nuclear equation of state with neutron star-black hole mergers

Authors: Nikhil Sarin, Hiranya V. Peiris, Daniel J. Mortlock, Justin Alsing, Samaya M. Nissanke, Stephen M. Feeney

Abstract: Gravitational-wave (GW) observations of neutron star-black hole (NSBH) mergers are sensitive to the nuclear equation of state (EOS). We present a new methodology for EOS inference with non-parametric Gaussian process (GP) priors, enabling direct constraints on the pressure at specific densities and the length-scale of correlations on the EOS. Using realistic simulations of NSBH mergers, incorporat… ▽ More Gravitational-wave (GW) observations of neutron star-black hole (NSBH) mergers are sensitive to the nuclear equation of state (EOS). We present a new methodology for EOS inference with non-parametric Gaussian process (GP) priors, enabling direct constraints on the pressure at specific densities and the length-scale of correlations on the EOS. Using realistic simulations of NSBH mergers, incorporating both GW and electromagnetic (EM) selection to ensure sample purity, we find that a GW detector network operating at O5-sensitivities will constrain the radius of a $\unit[1.4]{M_{\odot}}$ NS and the maximum NS mass with $1.6\%$ and $13\%$ precision, respectively. With the same sample, the projected constraint on the length-scale of correlations in the EOS is $\geq~\unit[3.2]{MeV~fm^{-3}}$. These results demonstrate strong potential for insights into the nuclear EOS from NSBH systems, provided they are robustly identified. △ Less

Submitted 9 July, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 10 pages, 4 Figures. Accepted in PRD

Generalized cold-atom simulators for vacuum decay

Authors: Alexander C. Jenkins, Ian G. Moss, Thomas P. Billam, Zoran Hadzibabic, Hiranya V. Peiris, Andrew Pontzen

Abstract: Cold-atom analogue experiments are a promising new tool for studying relativistic vacuum decay, allowing us to empirically probe early-Universe theories in the laboratory. However, existing analogue proposals place stringent requirements on the atomic scattering lengths that are challenging to realize experimentally. Here we generalize these proposals and show that any stable mixture between two s… ▽ More Cold-atom analogue experiments are a promising new tool for studying relativistic vacuum decay, allowing us to empirically probe early-Universe theories in the laboratory. However, existing analogue proposals place stringent requirements on the atomic scattering lengths that are challenging to realize experimentally. Here we generalize these proposals and show that any stable mixture between two states of a bosonic isotope can be used as a relativistic analogue. This greatly expands the range of suitable experimental setups, and will thus expedite efforts to study vacuum decay with cold atoms. △ Less

Submitted 17 November, 2023; v1 submitted 3 November, 2023; originally announced November 2023.

Comments: 11 pages, 4 figures, comments welcome; v2: reference added

Analog vacuum decay from vacuum initial conditions

Authors: Alexander C. Jenkins, Jonathan Braden, Hiranya V. Peiris, Andrew Pontzen, Matthew C. Johnson, Silke Weinfurtner

Abstract: Ultracold atomic gases can undergo phase transitions that mimic relativistic vacuum decay, allowing us to empirically test early-Universe physics in tabletop experiments. We investigate the physics of these analog systems, going beyond previous analyses of the classical equations of motion to study quantum fluctuations in the cold-atom false vacuum. We show that the fluctuation spectrum of this va… ▽ More Ultracold atomic gases can undergo phase transitions that mimic relativistic vacuum decay, allowing us to empirically test early-Universe physics in tabletop experiments. We investigate the physics of these analog systems, going beyond previous analyses of the classical equations of motion to study quantum fluctuations in the cold-atom false vacuum. We show that the fluctuation spectrum of this vacuum state agrees with the usual relativistic result in the regime where the classical analogy holds, providing further evidence for the suitability of these systems for studying vacuum decay. Using a suite of semiclassical lattice simulations, we simulate bubble nucleation from this analog vacuum state in a 1D homonuclear potassium-41 mixture, finding qualitative agreement with instanton predictions. We identify realistic parameters for this system that will allow us to study vacuum decay with current experimental capabilities, including a prescription for efficiently scanning over decay rates, and show that this setup will probe the quantum (rather than thermal) decay regime at temperatures $T\lesssim10\,\mathrm{nK}$. Our results help lay the groundwork for using upcoming cold-atom experiments as a new probe of nonperturbative early-Universe physics. △ Less

Submitted 4 January, 2024; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 17 pages, 11 figures, matches version published in PRD

Explaining dark matter halo density profiles with neural networks

Authors: Luisa Lucie-Smith, Hiranya V. Peiris, Andrew Pontzen

Abstract: We use explainable neural networks to connect the evolutionary history of dark matter halos with their density profiles. The network captures independent factors of variation in the density profiles within a low-dimensional representation, which we physically interpret using mutual information. Without any prior knowledge of the halos' evolution, the network recovers the known relation between the… ▽ More We use explainable neural networks to connect the evolutionary history of dark matter halos with their density profiles. The network captures independent factors of variation in the density profiles within a low-dimensional representation, which we physically interpret using mutual information. Without any prior knowledge of the halos' evolution, the network recovers the known relation between the early time assembly and the inner profile, and discovers that the profile beyond the virial radius is described by a single parameter capturing the most recent mass accretion rate. The results illustrate the potential for machine-assisted scientific discovery in complicated astrophysical datasets. △ Less

Submitted 19 January, 2024; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 7 pages, 5 figures. Minor changes to match version accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 132, 031001 (2024)

Towards Accurate Field-Level Inference of Massive Cosmic Structures

Authors: Stephen Stopyra, Hiranya V. Peiris, Andrew Pontzen, Jens Jasche, Guilhem Lavaux

Abstract: We investigate the accuracy requirements for field-level inference of cluster and void masses using data from galaxy surveys. We introduce a two-step framework that takes advantage of the fact that cluster masses are determined by flows on larger scales than the clusters themselves. First, we determine the integration accuracy required to perform field-level inference of cosmic initial conditions… ▽ More We investigate the accuracy requirements for field-level inference of cluster and void masses using data from galaxy surveys. We introduce a two-step framework that takes advantage of the fact that cluster masses are determined by flows on larger scales than the clusters themselves. First, we determine the integration accuracy required to perform field-level inference of cosmic initial conditions on these large scales, by fitting to late-time galaxy counts using the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. A 20-step COLA integrator is able to accurately describe the density field surrounding the most massive clusters in the Local Super-Volume ($<135\,h^{-1}\,\mathrm{Mpc}$), but does not by itself lead to converged virial mass estimates. Therefore we carry out `posterior resimulations', using full $N$-body dynamics while sampling from the inferred initial conditions, and thereby obtain estimates of masses for nearby massive clusters. We show that these are in broad agreement with existing estimates, and find that mass functions in the Local Super-Volume are compatible with $Λ$CDM. △ Less

Submitted 8 January, 2024; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 14 pages, 7 figures. Updated to match version accepted by MNRAS (published 2023/10/18). Minor changes and 1 new figure

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 1, January 2024, Pages 1244-1256

A robust estimator of mutual information for deep learning interpretability

Authors: Davide Piras, Hiranya V. Peiris, Andrew Pontzen, Luisa Lucie-Smith, Ningyuan Guo, Brian Nord

Abstract: We develop the use of mutual information (MI), a well-established metric in information theory, to interpret the inner workings of deep learning models. To accurately estimate MI from a finite number of samples, we present GMM-MI (pronounced $``$Jimmie$"$), an algorithm based on Gaussian mixture models that can be applied to both discrete and continuous settings. GMM-MI is computationally efficien… ▽ More We develop the use of mutual information (MI), a well-established metric in information theory, to interpret the inner workings of deep learning models. To accurately estimate MI from a finite number of samples, we present GMM-MI (pronounced $``$Jimmie$"$), an algorithm based on Gaussian mixture models that can be applied to both discrete and continuous settings. GMM-MI is computationally efficient, robust to the choice of hyperparameters and provides the uncertainty on the MI estimate due to the finite sample size. We extensively validate GMM-MI on toy data for which the ground truth MI is known, comparing its performance against established mutual information estimators. We then demonstrate the use of our MI estimator in the context of representation learning, working with synthetic data and physical datasets describing highly non-linear processes. We train deep learning models to encode high-dimensional data within a meaningful compressed (latent) representation, and use GMM-MI to quantify both the level of disentanglement between the latent variables, and their association with relevant physical quantities, thus unlocking the interpretability of the latent representation. We make GMM-MI publicly available. △ Less

Submitted 23 March, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: 30 pages, 8 figures. Minor changes to match version accepted for publication in Machine Learning: Science and Technology. GMM-MI available at https://github.com/dpiras/GMM-MI

Journal ref: Machine Learning: Science and Technology, Volume 4, Number 2, 025006, April 2023

Impact of Rubin Observatory cadence choices on supernovae photometric classification

Authors: Catarina S. Alves, Hiranya V. Peiris, Michelle Lochner, Jason D. McEwen, Richard Kessler

Abstract: The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will discover an unprecedented number of supernovae (SNe), making spectroscopic classification for all the events infeasible. LSST will thus rely on photometric classification, whose accuracy depends on the not-yet-finalized LSST observing strategy. In this work, we analyze the impact of cadence choices on classification perfor… ▽ More The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will discover an unprecedented number of supernovae (SNe), making spectroscopic classification for all the events infeasible. LSST will thus rely on photometric classification, whose accuracy depends on the not-yet-finalized LSST observing strategy. In this work, we analyze the impact of cadence choices on classification performance using simulated multi-band light curves. First, we simulate SNe with an LSST baseline cadence, a non-rolling cadence, and a presto-color cadence which observes each sky location three times per night instead of twice. Each simulated dataset includes a spectroscopically-confirmed training set, which we augment to be representative of the test set as part of the classification pipeline. Then, we use the photometric transient classification library snmachine to build classifiers. We find that the active region of the rolling cadence used in the baseline observing strategy yields a 25% improvement in classification performance relative to the background region. This improvement in performance in the actively-rolling region is also associated with an increase of up to a factor of 2.7 in the number of cosmologically-useful Type Ia supernovae relative to the background region. However, adding a third visit per night as implemented in presto-color degrades classification performance due to more irregularly sampled light curves. Overall, our results establish desiderata on the observing cadence related to classification of full SNe light curves, which in turn impacts photometric SNe cosmology with LSST. △ Less

Submitted 15 March, 2023; v1 submitted 27 October, 2022; originally announced October 2022.

Comments: 22 pages, 14 figures. Changed to match version accepted by the Astrophysical Journal Supplement Series (accepted 06/02/2023)

Searching For Dark Matter with Plasma Haloscopes

Authors: Alexander J. Millar, Steven M. Anlage, Rustam Balafendiev, Pavel Belov, Karl van Bibber, Jan Conrad, Marcel Demarteau, Alexander Droster, Katherine Dunne, Andrea Gallo Rosso, Jon E. Gudmundsson, Heather Jackson, Gagandeep Kaur, Tove Klaesson, Nolan Kowitt, Matthew Lawson, Alexander Leder, Akira Miyazaki, Sid Morampudi, Hiranya V. Peiris, Henrik S. Røising, Gaganpreet Singh, Dajie Sun, Jacob H. Thomas, Frank Wilczek , et al. (2 additional authors not shown)

Abstract: We summarise the recent progress of the Axion Longitudinal Plasma HAloscope (ALPHA) Consortium, a new experimental collaboration to build a plasma haloscope to search for axions and dark photons. The plasma haloscope is a novel method for the detection of the resonant conversion of light dark matter to photons. ALPHA will be sensitive to QCD axions over almost a decade of parameter space, potentia… ▽ More We summarise the recent progress of the Axion Longitudinal Plasma HAloscope (ALPHA) Consortium, a new experimental collaboration to build a plasma haloscope to search for axions and dark photons. The plasma haloscope is a novel method for the detection of the resonant conversion of light dark matter to photons. ALPHA will be sensitive to QCD axions over almost a decade of parameter space, potentially discovering dark matter and resolving the Strong CP problem. Unlike traditional cavity haloscopes, which are generally limited in volume by the Compton wavelength of the dark matter, plasma haloscopes use a wire metamaterial to create a tuneable artificial plasma frequency, decoupling the wavelength of light from the Compton wavelength and allowing for much stronger signals. We develop the theoretical foundations of plasma haloscopes and discuss recent experimental progress. Finally, we outline a baseline design for ALPHA and show that a full-scale experiment could discover QCD axions over almost a decade of parameter space. △ Less

Submitted 22 March, 2023; v1 submitted 30 September, 2022; originally announced October 2022.

Comments: Endorsers: Jens Dilling, Michael Febbraro, Stefan Knirck, and Claire Marvinney. 26 pages, 17 figures, version accepted in Physical Review D, typo corrected

Report number: FERMILAB-PUB-22-739-T

Journal ref: Physical Review D 107, 055013 (2023)

Stellar angular momentum can be controlled from cosmological initial conditions

Authors: Corentin Cadiou, Andrew Pontzen, Hiranya V. Peiris

Abstract: The angular momentum of galaxies controls the kinematics of their stars, which in turn drives observable quantities such as the apparent radius, the bulge fraction, and the alignment with other nearby structures. To show how angular momentum of galaxies is determined, we build high (${35}\,\mathrm{pc}$) resolution numerical experiments in which we increase or decrease the angular momentum of the L… ▽ More The angular momentum of galaxies controls the kinematics of their stars, which in turn drives observable quantities such as the apparent radius, the bulge fraction, and the alignment with other nearby structures. To show how angular momentum of galaxies is determined, we build high (${35}\,\mathrm{pc}$) resolution numerical experiments in which we increase or decrease the angular momentum of the Lagrangian patches in the early universe. We simulate three galaxies over their histories from $z=200$ to $z=2$, each with five different choices for the angular momentum (fifteen simulations in total). Our results show that altering early-universe angular momentum changes the timing and orbital parameters of mergers, which in turn changes the total stellar angular momentum within a galaxy's virial radius in a predictable manner. Of our three galaxies, one has no large satellite at $z=2$; in this case, the specific angular momentum is concentrated in the central galaxy. We modify its stellar angular momentum over $0.7\,\mathrm{dex}$ (from $61$ to $320\,\mathrm{kpc.km.s^{-1}}$) and show that this causes its effective radius to grow by $40\,\%$, its $v/σ$ parameter to grow by a factor $\times 2.6$ and its bulge fraction to decrease from $0.72$ to $0.57$. The ability to control angular momentum will allow future studies to probe the causal origin of scaling relations between galaxy mass, angular momentum and morphology, and to better understand the origin of galactic intrinsic alignments. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Comments: 8 pages, 3 figures. Submitted to MNRAS. Comments welcome

Modelling Populations of Kilonovae

Authors: Christian N. Setzer, Hiranya V. Peiris, Oleg Korobkin, Stephan Rosswog

Abstract: The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements, and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection curr… ▽ More The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements, and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models which are quick to evaluate, and are calibrated with more detailed multi-physics simulations, are needed to design observational strategies for kilonovae detection, and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from BNS mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of binary neutron star (BNS) mergers produce peak i-band absolute magnitudes $-20 \leq M_i \leq -11$. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution. △ Less

Submitted 6 March, 2023; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 14 pages, 12 figures, 2 tables, matches version accepted to MNRAS. Fixes a software bug from previous versions which incorrectly determined the kilonovae ejecta velocities, v4: added link to public code package https://github.com/cnsetzer/kilopop

Journal ref: MNRAS, Volume 520, Issue 2, April 2023, Pages 2829-2842

Mass Renormalization in Lattice Simulations of False Vacuum Decay

Authors: Jonathan Braden, Matthew C. Johnson, Hiranya V. Peiris, Andrew Pontzen, Silke Weinfurtner

Abstract: False vacuum decay, a quantum mechanical first-order phase transition in scalar field theories, is an important phenomenon in early universe cosmology. Recently, real-time semi-classical techniques based on ensembles of lattice simulations were applied to the problem of false vacuum decay. In this context, or any other lattice simulation, the effective potential experienced by long-wavelength mode… ▽ More False vacuum decay, a quantum mechanical first-order phase transition in scalar field theories, is an important phenomenon in early universe cosmology. Recently, real-time semi-classical techniques based on ensembles of lattice simulations were applied to the problem of false vacuum decay. In this context, or any other lattice simulation, the effective potential experienced by long-wavelength modes is not the same as the bare potential. To make quantitative predictions using the real-time semi-classical techniques, it is therefore necessary to understand the redefinition of model parameters and the corresponding deformation of the vacuum state, as well as stochastic contributions that require modeling of unresolved subgrid modes. In this work, we focus on the former corrections and compute the expected modification of the true and false vacuum effective mass, which manifests as a modified dispersion relationship for linear fluctuations about the vacuum. We compare these theoretical predictions to numerical simulations and find excellent agreement. Motivated by this, we use the effective masses to fix the shape of a parameterized effective potential, and explore the modeling uncertainty associated with non-linear corrections. We compute the decay rates in both the Euclidean and real-time formalisms, finding qualitative agreement in the dependence on the UV cutoff. These calculations further demonstrate that a quantitative understanding of the rates requires additional corrections. △ Less

Submitted 24 April, 2023; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 15 pages, 5 figures, accepted. Additional clarifications added. Main results unchanged

Journal ref: Phys. Rev. D 107 (2023) 8, 083509

Axion Dark Matter

Authors: C. B. Adams, N. Aggarwal, A. Agrawal, R. Balafendiev, C. Bartram, M. Baryakhtar, H. Bekker, P. Belov, K. K. Berggren, A. Berlin, C. Boutan, D. Bowring, D. Budker, A. Caldwell, P. Carenza, G. Carosi, R. Cervantes, S. S. Chakrabarty, S. Chaudhuri, T. Y. Chen, S. Cheong, A. Chou, R. T. Co, J. Conrad, D. Croon , et al. (130 additional authors not shown)

Abstract: Axions are well-motivated dark matter candidates with simple cosmological production mechanisms. They were originally introduced to solve the strong CP problem, but also arise in a wide range of extensions to the Standard Model. This Snowmass white paper summarizes axion phenomenology and outlines next-generation laboratory experiments proposed to detect axion dark matter. There are vibrant synerg… ▽ More Axions are well-motivated dark matter candidates with simple cosmological production mechanisms. They were originally introduced to solve the strong CP problem, but also arise in a wide range of extensions to the Standard Model. This Snowmass white paper summarizes axion phenomenology and outlines next-generation laboratory experiments proposed to detect axion dark matter. There are vibrant synergies with astrophysical searches and advances in instrumentation including quantum-enabled readout, high-Q resonators and cavities and large high-field magnets. This white paper outlines a clear roadmap to discovery, and shows that the US is well-positioned to be at the forefront of the search for axion dark matter in the coming decade. △ Less

Submitted 29 March, 2023; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: restore and expand author list

Discovering the building blocks of dark matter halo density profiles with neural networks

Authors: Luisa Lucie-Smith, Hiranya V. Peiris, Andrew Pontzen, Brian Nord, Jeyan Thiyagalingam, Davide Piras

Abstract: The density profiles of dark matter halos are typically modeled using empirical formulae fitted to the density profiles of relaxed halo populations. We present a neural network model that is trained to learn the mapping from the raw density field containing each halo to the dark matter density profile. We show that the model recovers the widely-used Navarro-Frenk-White (NFW) profile out to the vir… ▽ More The density profiles of dark matter halos are typically modeled using empirical formulae fitted to the density profiles of relaxed halo populations. We present a neural network model that is trained to learn the mapping from the raw density field containing each halo to the dark matter density profile. We show that the model recovers the widely-used Navarro-Frenk-White (NFW) profile out to the virial radius, and can additionally describe the variability in the outer profile of the halos. The neural network architecture consists of a supervised encoder-decoder framework, which first compresses the density inputs into a low-dimensional latent representation, and then outputs $ρ(r)$ for any desired value of radius $r$. The latent representation contains all the information used by the model to predict the density profiles. This allows us to interpret the latent representation by quantifying the mutual information between the representation and the halos' ground-truth density profiles. A two-dimensional representation is sufficient to accurately model the density profiles up to the virial radius; however, a three-dimensional representation is required to describe the outer profiles beyond the virial radius. The additional dimension in the representation contains information about the infalling material in the outer profiles of dark matter halos, thus discovering the splashback boundary of halos without prior knowledge of the halos' dynamical history. △ Less

Submitted 13 May, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 12 pages, 6 figures. Minor changes to match version accepted for publication in PRD

Next Generation Observatories -- Report from the Dawn VI Workshop; October 5-7 2021

Authors: D. H. Shoemaker, Stefan Ballmer, Matteo Barsuglia, E. Berger, Emanuele Berti, Duncan A. Brown, Poonam Chandra, Matthew Evans, Ke Fang, Wen-fai Fong, Andreas Freise, Peter Fritschel, Jenny Greene, C. J. Horowitz, Jeff Kissel, Brian Lantz, Paul D. Lasky, Harald Lueck, M. Coleman Miller, Alexander H. Nitz, David Ottaway, Hiranya V. Peiris, Michele Punturo, D. H. Reitze, Gary H. Sanders , et al. (11 additional authors not shown)

Abstract: The workshop Dawn VI: Next Generation Observatories took place online over three days, 5-7 October, 2021. More than 200 physicists and astronomers attended to contribute to, and learn from, a discussion of next-generation ground-based gravitational-wave detectors. The program was centered on the next generation of ground-based gravitational-wave observatories and their synergy with the greater lan… ▽ More The workshop Dawn VI: Next Generation Observatories took place online over three days, 5-7 October, 2021. More than 200 physicists and astronomers attended to contribute to, and learn from, a discussion of next-generation ground-based gravitational-wave detectors. The program was centered on the next generation of ground-based gravitational-wave observatories and their synergy with the greater landscape of scientific observatories of the 2030s. Cosmic Explorer (CE), a concept developed with US National Science Foundation support, was a particular focus; Einstein Telescope (ET), the European next generation concept, is an important complement and partner in forming a network. The concluding summary of the meeting expressed the sentiment that the observational science accessible to CE and ET, also in combination with data from other non-GW observatories, will stimulate a very broad community of analysts and yield insights which are exciting given the access to GWs from the entire universe. The need, and desire, for closer collaboration between ET and CE was expressed; a three-detector network is optimal for delivering much of the science. The science opportunities afforded by CE and ET are broad and compelling, impacting a wide range of disciplines in physics and high energy astrophysics. There was a consensus that CE is a concept that can deliver the promised science. A strong endorsement of Cosmic Explorer, as described in the CE Horizon Study, is a primary outcome of DAWN VI. △ Less

Submitted 20 February, 2022; v1 submitted 23 December, 2021; originally announced December 2021.

Comments: Proceedings of the Dawn VI workshop

New limits on light dark matter - proton cross section from the cosmic large-scale structure

Authors: Keir K. Rogers, Cora Dvorkin, Hiranya V. Peiris

Abstract: We set the strongest limits to-date on the velocity-independent dark matter (DM) - proton cross section $σ$ for DM masses $m = 10\,\mathrm{keV}$ to $100\,\mathrm{GeV}$, using large-scale structure traced by the Lyman-alpha forest: e.g., a 95% lower limit $σ< 6 \times 10^{-30}\,\mathrm{cm}^2$, for $m = 100\,\mathrm{keV}$. Our results complement direct detection, which has limited sensitivity to sub… ▽ More We set the strongest limits to-date on the velocity-independent dark matter (DM) - proton cross section $σ$ for DM masses $m = 10\,\mathrm{keV}$ to $100\,\mathrm{GeV}$, using large-scale structure traced by the Lyman-alpha forest: e.g., a 95% lower limit $σ< 6 \times 10^{-30}\,\mathrm{cm}^2$, for $m = 100\,\mathrm{keV}$. Our results complement direct detection, which has limited sensitivity to sub-GeV DM. We use an emulator of cosmological simulations, combined with data from the smallest cosmological scales used to-date, to model and search for the imprint of primordial DM-proton collisions. Cosmological bounds are improved by up to a factor of 25. △ Less

Submitted 4 April, 2022; v1 submitted 19 November, 2021; originally announced November 2021.

Comments: 12 pages, 5 figures. Minor changes to match version accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 128, 171301 (2022)

The large-scale environment of thermonuclear and core-collapse supernovae

Authors: Eleni Tsaprazi, Jens Jasche, Ariel Goobar, Hiranya V. Peiris, Igor Andreoni, Michael W. Coughlin, Christoffer U. Fremling, Matthew J. Graham, Mansi Kasliwal, Shri R. Kulkarni, Ashish A. Mahabal, Reed Riddle, Jesper Sollerman, Anastasios Tzanidakis

Abstract: The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of $\sim$ 3.8 Mpc, constrained by the 2M++ galaxy survey. We limit our analysis to redshift… ▽ More The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of $\sim$ 3.8 Mpc, constrained by the 2M++ galaxy survey. We limit our analysis to redshift $z<0.036$, using samples of 498 thermonuclear and 782 core-collapse SNe from the Zwicky Transient Facility's Bright Transient Survey and Census of the Local Universe catalogues. We detect clustering of SNe with high significance; the observed clustering of the two SNe populations is consistent with each other. Further, the clustering of SN hosts is consistent with that of the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) DR12 spectroscopic galaxy sample in the same redshift range. Using a tidal shear classifier, we classify the LSS into voids, sheets, filaments and knots. We find that both SNe and SDSS galaxies are predominantly found in sheets and filaments. SNe are significantly under-represented in voids and over-represented in knots compared to the volume fraction in these structures. This work opens the potential for using forthcoming wide-field deep SN surveys as a complementary LSS probe. △ Less

Submitted 1 December, 2021; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 8 pages, 4 figures, accepted by MNRAS

Considerations for optimizing photometric classification of supernovae from the Rubin Observatory

Authors: Catarina S. Alves, Hiranya V. Peiris, Michelle Lochner, Jason D. McEwen, Tarek Allam Jr, Rahul Biswas

Abstract: The Vera C. Rubin Observatory will increase the number of observed supernovae (SNe) by an order of magnitude; however, it is impossible to spectroscopically confirm the class for all the SNe discovered. Thus, photometric classification is crucial but its accuracy depends on the not-yet-finalized observing strategy of Rubin Observatory's Legacy Survey of Space and Time (LSST). We quantitatively ana… ▽ More The Vera C. Rubin Observatory will increase the number of observed supernovae (SNe) by an order of magnitude; however, it is impossible to spectroscopically confirm the class for all the SNe discovered. Thus, photometric classification is crucial but its accuracy depends on the not-yet-finalized observing strategy of Rubin Observatory's Legacy Survey of Space and Time (LSST). We quantitatively analyze the impact of the LSST observing strategy on SNe classification using simulated multi-band light curves from the Photometric LSST Astronomical Time-Series Classification Challenge (PLAsTiCC). First, we augment the simulated training set to be representative of the photometric redshift distribution per supernovae class, the cadence of observations, and the flux uncertainty distribution of the test set. Then we build a classifier using the photometric transient classification library snmachine, based on wavelet features obtained from Gaussian process fits, yielding similar performance to the winning PLAsTiCC entry. We study the classification performance for SNe with different properties within a single simulated observing strategy. We find that season length is important, with light curves of 150 days yielding the highest performance. Cadence also has an important impact on SNe classification; events with median inter-night gap <3.5 days yield higher classification performance. Interestingly, we find that large gaps (>10 days) in light curve observations do not impact performance if sufficient observations are available on either side, due to the effectiveness of the Gaussian process interpolation. This analysis is the first exploration of the impact of observing strategy on photometric supernova classification with LSST. △ Less

Submitted 29 October, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 18 pages, 13 figures. Changed to match version accepted by the Astrophysical Journal Supplement Series (accepted 28/10/2021). Software publicly available at https://github.com/LSSTDESC/snmachine

Journal ref: 2022 ApJS 258 23

Quantifying the Rarity of the Local Super-Volume

Authors: Stephen Stopyra, Hiranya V. Peiris, Andrew Pontzen, Jens Jasche, Priyamvada Natarajan

Abstract: We investigate the extent to which the number of clusters of mass exceeding $10^{15}\,M_{\odot}\,h^{-1}$ within the local super-volume ($<135\mathrm{\,Mpc}h^{-1}$) is compatible with the standard $Λ$CDM cosmological model. Depending on the mass estimator used, we find that the observed number $N$ of such massive structures can vary between $0$ and $5$. Adopting $N=5$ yields $Λ$CDM likelihoods as l… ▽ More We investigate the extent to which the number of clusters of mass exceeding $10^{15}\,M_{\odot}\,h^{-1}$ within the local super-volume ($<135\mathrm{\,Mpc}h^{-1}$) is compatible with the standard $Λ$CDM cosmological model. Depending on the mass estimator used, we find that the observed number $N$ of such massive structures can vary between $0$ and $5$. Adopting $N=5$ yields $Λ$CDM likelihoods as low as $2.4\times 10^{-3}$ (with $σ_8=0.81$) or $3.8\times 10^{-5}$ (with $σ_8=0.74$). However, at the other extreme ($N=0$), the likelihood is of order unity. Thus, while potentially very powerful, this method is currently limited by systematic uncertainties in cluster mass estimates. This motivates efforts to reduce these systematics with additional observations and improved modelling. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 7 pages, 3 figures, 1 table

The causal effect of environment on halo mass and concentration

Authors: Corentin Cadiou, Andrew Pontzen, Hiranya V. Peiris, Luisa Lucie-Smith

Abstract: Understanding the impact of environment on the formation and evolution of dark matter halos and galaxies is a crucial open problem. Studying statistical correlations in large simulated populations sheds some light on these impacts, but the causal effect of an environment on individual objects is harder to pinpoint. Addressing this, we present a new method for resimulating a single dark matter halo… ▽ More Understanding the impact of environment on the formation and evolution of dark matter halos and galaxies is a crucial open problem. Studying statistical correlations in large simulated populations sheds some light on these impacts, but the causal effect of an environment on individual objects is harder to pinpoint. Addressing this, we present a new method for resimulating a single dark matter halo in multiple large-scale environments. In the initial conditions, we 'splice' (i.e. insert) the Lagrangian region of a halo into different Gaussian random fields, while enforcing consistency with the statistical properties of $Λ$CDM. Applying this technique, we demonstrate that the mass of halos is primarily determined by the density structure inside their Lagrangian patches, while the halos' concentration is more strongly affected by environment. The splicing approach will also allow us to study, for example, the impact of the cosmic web on accretion processes and galaxy quenching. △ Less

Submitted 13 September, 2021; v1 submitted 7 July, 2021; originally announced July 2021.

Comments: 6 pages, 5 figures. Accepted 2021 September 10. Received 2021 September 9; in original form 2021 July 7

Results of the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC)

Authors: R. Hložek, K. A. Ponder, A. I. Malz, M. Dai, G. Narayan, E. E. O. Ishida, T. Allam Jr, A. Bahmanyar, R. Biswas, L. Galbany, S. W. Jha, D. O. Jones, R. Kessler, M. Lochner, A. A. Mahabal, K. S. Mandel, J. R. Martínez-Galarza, J. D. McEwen, D. Muthukrishna, H. V. Peiris, C. M. Peters, C. N. Setzer

Abstract: Next-generation surveys like the Legacy Survey of Space and Time (LSST) on the Vera C. Rubin Observatory will generate orders of magnitude more discoveries of transients and variable stars than previous surveys. To prepare for this data deluge, we developed the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC), a competition which aimed to catalyze the development of ro… ▽ More Next-generation surveys like the Legacy Survey of Space and Time (LSST) on the Vera C. Rubin Observatory will generate orders of magnitude more discoveries of transients and variable stars than previous surveys. To prepare for this data deluge, we developed the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC), a competition which aimed to catalyze the development of robust classifiers under LSST-like conditions of a non-representative training set for a large photometric test set of imbalanced classes. Over 1,000 teams participated in PLAsTiCC, which was hosted in the Kaggle data science competition platform between Sep 28, 2018 and Dec 17, 2018, ultimately identifying three winners in February 2019. Participants produced classifiers employing a diverse set of machine learning techniques including hybrid combinations and ensemble averages of a range of approaches, among them boosted decision trees, neural networks, and multi-layer perceptrons. The strong performance of the top three classifiers on Type Ia supernovae and kilonovae represent a major improvement over the current state-of-the-art within astronomy. This paper summarizes the most promising methods and evaluates their results in detail, highlighting future directions both for classifier development and simulation needs for a next generation PLAsTiCC data set. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: 20 pages, 14 figures

Prospects for Measuring the Hubble Constant with Neutron-Star-Black-Hole Mergers

Authors: Stephen M. Feeney, Hiranya V. Peiris, Samaya M. Nissanke, Daniel J. Mortlock

Abstract: Gravitational wave (GW) and electromagnetic (EM) observations of neutron-star-black-hole (NSBH) mergers can provide precise local measurements of the Hubble constant ($H_0$), ideal for resolving the current $H_0$ tension. We perform end-to-end analyses of realistic populations of simulated NSBHs, incorporating both GW and EM selection for the first time. We show that NSBHs could achieve unbiased 1… ▽ More Gravitational wave (GW) and electromagnetic (EM) observations of neutron-star-black-hole (NSBH) mergers can provide precise local measurements of the Hubble constant ($H_0$), ideal for resolving the current $H_0$ tension. We perform end-to-end analyses of realistic populations of simulated NSBHs, incorporating both GW and EM selection for the first time. We show that NSBHs could achieve unbiased 1.5-2.4% precision $H_0$ estimates by 2030. The achievable precision is strongly affected by the details of spin precession and tidal disruption, highlighting the need for improved modeling of NSBH mergers. △ Less

Submitted 29 April, 2021; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: Version published in Physical Review Letters. 10 pages, 3 figures. Code available at https://github.com/sfeeney/nsbh

Journal ref: Phys. Rev. Lett. 126, 171102 (2021)

Angular momentum evolution can be predicted from cosmological initial conditions

Authors: Corentin Cadiou, Andrew Pontzen, Hiranya V. Peiris

Abstract: The angular momentum of dark matter haloes controls their spin magnitude and orientation, which in turn influences the galaxies therein. However, the process by which dark matter haloes acquire angular momentum is not fully understood; in particular, it is unclear whether angular momentum growth is stochastic. To address this question, we extend the genetic modification technique to allow control… ▽ More The angular momentum of dark matter haloes controls their spin magnitude and orientation, which in turn influences the galaxies therein. However, the process by which dark matter haloes acquire angular momentum is not fully understood; in particular, it is unclear whether angular momentum growth is stochastic. To address this question, we extend the genetic modification technique to allow control over the angular momentum of any region in the initial conditions. Using this technique to produce a sequence of modified simulations, we can then investigate whether changes to the angular momentum of a specified region in the evolved universe can be accurately predicted from changes in the initial conditions alone. We find that the angular momentum in regions with modified initial conditions can be predicted between 2 and 4 times more accurately than expected from applying tidal torque theory. This result is masked when analysing the angular momentum of haloes, because particles in the outskirts of haloes dominate the angular momentum budget. We conclude that the angular momentum of Lagrangian patches is highly predictable from the initial conditions, with apparent chaotic behaviour being driven by stochastic changes to the arbitrary boundary defining the halo. △ Less

Submitted 9 February, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: 7 pages, 4 figures. Accepted for publication in MNRAS: 2021 February 3. Received 2021 January 30; in original form 2020 December 14

An emulator for the Lyman-$α$ forest in beyond-$Λ$CDM cosmologies

Authors: Christian Pedersen, Andreu Font-Ribera, Keir K. Rogers, Patrick McDonald, Hiranya V. Peiris, Andrew Pontzen, Anže Slosar

Abstract: Interpreting observations of the Lyman-$α$ forest flux power spectrum requires interpolation between a small number of expensive simulations. We present a Gaussian process emulator modelling the 1D flux power spectrum as a function of the amplitude and slope of the small-scale linear matter power spectrum, and the state of the intergalactic medium at the epoch of interest ($2 < z < 4$). This param… ▽ More Interpreting observations of the Lyman-$α$ forest flux power spectrum requires interpolation between a small number of expensive simulations. We present a Gaussian process emulator modelling the 1D flux power spectrum as a function of the amplitude and slope of the small-scale linear matter power spectrum, and the state of the intergalactic medium at the epoch of interest ($2 < z < 4$). This parameterisation enables the prediction of the flux power spectrum in extended cosmological models that are not explicitly included in the training set, eliminating the need to construct bespoke emulators for a number of extensions to $Λ$CDM. Our emulator is appropriate for cosmologies in which the linear matter power spectrum is described to percent level accuracy by just an amplitude and slope across the epoch of interest, and in the regime probed by eBOSS/DESI data. We demonstrate this for massive neutrino cosmologies, where the emulator is able to predict the flux power spectrum in a $Σm_ν=0.3$ eV neutrino cosmology to sub-percent accuracy, without including massive neutrinos in the training simulations. Further parameters would be required to describe models with sharp features in the linear power, such as warm or light axion dark matter. This work will facilitate the combination of upcoming DESI data with observations of the cosmic microwave background, to obtain constraints on neutrino mass and other extensions to $Λ$CDM cosmology. △ Less

Submitted 14 May, 2021; v1 submitted 30 November, 2020; originally announced November 2020.

Comments: 26 pages, 11 figures. Minor changes to match version published in JCAP

Journal ref: JCAP 05 (2021) 033

Deep learning insights into cosmological structure formation

Authors: Luisa Lucie-Smith, Hiranya V. Peiris, Andrew Pontzen, Brian Nord, Jeyan Thiyagalingam

Abstract: The evolution of linear initial conditions present in the early universe into extended halos of dark matter at late times can be computed using cosmological simulations. However, a theoretical understanding of this complex process remains elusive; in particular, the role of anisotropic information in the initial conditions in establishing the final mass of dark matter halos remains a long-standing… ▽ More The evolution of linear initial conditions present in the early universe into extended halos of dark matter at late times can be computed using cosmological simulations. However, a theoretical understanding of this complex process remains elusive; in particular, the role of anisotropic information in the initial conditions in establishing the final mass of dark matter halos remains a long-standing puzzle. Here, we build a deep learning framework to investigate this question. We train a three-dimensional convolutional neural network (CNN) to predict the mass of dark matter halos from the initial conditions, and quantify in full generality the amounts of information in the isotropic and anisotropic aspects of the initial density field about final halo masses. We find that anisotropies add a small, albeit statistically significant amount of information over that contained within spherical averages of the density field about final halo mass. However, the overall scatter in the final mass predictions does not change qualitatively with this additional information, only decreasing from 0.9 dex to 0.7 dex. Given such a small improvement, our results demonstrate that isotropic aspects of the initial density field essentially saturate the relevant information about final halo mass. Therefore, instead of searching for information directly encoded in initial conditions anisotropies, a more promising route to accurate, fast halo mass predictions is to add approximate dynamical information based e.g. on perturbation theory. More broadly, our results indicate that deep learning frameworks can provide a powerful tool for extracting physical insight into cosmological structure formation. △ Less

Submitted 1 March, 2024; v1 submitted 20 November, 2020; originally announced November 2020.

Comments: 17 pages, 10 figures. Accepted in PRD

Journal ref: Phys. Rev. D 109, 063524 (2024)

How to Build a Catalogue of Linearly-Evolving Cosmic Voids

Authors: Stephen Stopyra, Hiranya V. Peiris, Andrew Pontzen

Abstract: Cosmic voids provide a powerful probe of the origin and evolution of structures in the Universe because their dynamics can remain near-linear to the present day. As a result they have the potential to connect large scale structure at late times to early-Universe physics. Existing "watershed"-based algorithms, however, define voids in terms of their morphological properties at low redshift. The deg… ▽ More Cosmic voids provide a powerful probe of the origin and evolution of structures in the Universe because their dynamics can remain near-linear to the present day. As a result they have the potential to connect large scale structure at late times to early-Universe physics. Existing "watershed"-based algorithms, however, define voids in terms of their morphological properties at low redshift. The degree to which the resulting regions exhibit linear dynamics is consequently uncertain, and there is no direct connection to their evolution from the initial density field. A recent void definition addresses these issues by considering "anti-halos". This approach consists of inverting the initial conditions of an $N$-body simulation to swap overdensities and underdensities. After evolving the pair of initial conditions, anti-halos are defined by the particles within the inverted simulation that are inside halos in the original (uninverted) simulation. In this work, we quantify the degree of non-linearity of both anti-halos and watershed voids using the Zel'dovich approximation. We find that non-linearities are introduced by voids with radii less than $5\,\mathrm{Mpc}\,h^{-1}$, and that both anti-halos and watershed voids can be made into highly linear sets by removing these voids. △ Less

Submitted 24 November, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 8 pages, 5 figures. Changed to match version accepted by MNRAS (accepted 17/11/2020), minor clarifications added and no changes to the conclusions

General framework for cosmological dark matter bounds using $N$-body simulations

Authors: Keir K. Rogers, Hiranya V. Peiris

Abstract: We present a general framework for obtaining robust bounds on the nature of dark matter using cosmological $N$-body simulations and Lyman-alpha forest data. We construct an emulator of hydrodynamical simulations, which is a flexible, accurate and computationally-efficient model for predicting the response of the Lyman-alpha forest flux power spectrum to different dark matter models, the state of t… ▽ More We present a general framework for obtaining robust bounds on the nature of dark matter using cosmological $N$-body simulations and Lyman-alpha forest data. We construct an emulator of hydrodynamical simulations, which is a flexible, accurate and computationally-efficient model for predicting the response of the Lyman-alpha forest flux power spectrum to different dark matter models, the state of the intergalactic medium (IGM) and the primordial power spectrum. The emulator combines a flexible parameterization for the small-scale suppression in the matter power spectrum arising in "non-cold" dark matter models, with an improved IGM model. We then demonstrate how to optimize the emulator for the case of ultra-light axion dark matter, presenting tests of convergence. We also carry out cross-validation tests of the accuracy of flux power spectrum prediction. This framework can be optimized for the analysis of many other dark matter candidates, e.g., warm or interacting dark matter. Our work demonstrates that a combination of an optimized emulator and cosmological "effective theories," where many models are described by a single set of equations, is a powerful approach for robust and computationally-efficient inference from the cosmic large-scale structure. △ Less

Submitted 1 February, 2021; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 23 pages, 11 figures. Minor changes to match version accepted for publication in PRD. Companion letter: arXiv:2007.12705 "Strong bound on canonical ultra-light axion dark matter from the Lyman-alpha forest" (K. Rogers & H. Peiris)

Journal ref: Phys. Rev. D 103, 043526 (2021)

Strong bound on canonical ultra-light axion dark matter from the Lyman-alpha forest

Authors: Keir K. Rogers, Hiranya V. Peiris

Abstract: We present a new bound on the ultra-light axion (ULA) dark matter mass $m_\text{a}$, using the Lyman-alpha forest to look for suppressed cosmic structure growth: a 95% lower limit $m_\text{a} > 2 \times 10^{-20}\,\text{eV}$. This strongly disfavors ($> 99.7\%$ credibility) the canonical ULA with $10^{-22}\,\text{eV} < m_\text{a} < 10^{-21}\,\text{eV}$, motivated by the string axiverse and solution… ▽ More We present a new bound on the ultra-light axion (ULA) dark matter mass $m_\text{a}$, using the Lyman-alpha forest to look for suppressed cosmic structure growth: a 95% lower limit $m_\text{a} > 2 \times 10^{-20}\,\text{eV}$. This strongly disfavors ($> 99.7\%$ credibility) the canonical ULA with $10^{-22}\,\text{eV} < m_\text{a} < 10^{-21}\,\text{eV}$, motivated by the string axiverse and solutions to possible tensions in the cold dark matter model. We strengthen previous equivalent bounds by about an order of magnitude. We demonstrate the robustness of our results using an optimized emulator of improved hydrodynamical simulations. △ Less

Submitted 24 December, 2020; v1 submitted 24 July, 2020; originally announced July 2020.

Comments: 8 pages, 4 figures. Minor changes to match version accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 126, 071302 (2021)

Bayesian inflationary reconstructions from Planck 2018 data

Authors: Will J. Handley, Anthony N. Lasenby, Hiranya V. Peiris, Michael P. Hobson

Abstract: We present three non-parametric Bayesian primordial reconstructions using Planck 2018 polarization data: linear spline primordial power spectrum reconstructions, cubic spline inflationary potential reconstructions and sharp-featured primordial power spectrum reconstructions. All three methods conditionally show hints of an oscillatory feature in the primordial power spectrum in the multipole range… ▽ More We present three non-parametric Bayesian primordial reconstructions using Planck 2018 polarization data: linear spline primordial power spectrum reconstructions, cubic spline inflationary potential reconstructions and sharp-featured primordial power spectrum reconstructions. All three methods conditionally show hints of an oscillatory feature in the primordial power spectrum in the multipole range $\ell\sim20$ to $\ell\sim50$, which is to some extent preserved upon marginalization. We find no evidence for deviations from a pure power law across a broad observable window ($50\lesssim\ell\lesssim2000$), but find that parameterizations are preferred which are able to account for lack of resolution at large angular scales due to cosmic variance, and at small angular scales due to Planck instrument noise. Furthermore, the late-time cosmological parameters are unperturbed by these extensions to the primordial power spectrum. This work is intended to provide a background and give more details of the Bayesian primordial reconstruction work found in the Planck 2018 papers. △ Less

Submitted 11 November, 2019; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: 25 Pages, 15 Figures

Journal ref: Phys. Rev. D 100, 103511 (2019)

Planck 2018 results. V. CMB power spectra and likelihoods

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso , et al. (143 additional authors not shown)

Abstract: This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low… ▽ More This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $τ$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $τ$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $Λ$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$σ$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged) △ Less

Submitted 15 September, 2020; v1 submitted 30 July, 2019; originally announced July 2019.

Comments: Revised to match version published in Astronomy & Astrophysics

Journal ref: A&A 641, A5 (2020)

Target Neutrino Mass Precision for Determining the Neutrino Hierarchy

Authors: Constance Mahony, Boris Leistedt, Hiranya V. Peiris, Jonathan Braden, Benjamin Joachimi, Andreas Korn, Linda Cremonesi, Ryan Nichol

Abstract: Recent works combining neutrino oscillation and cosmological data to determine the neutrino hierarchy found a range of odds in favour of the normal hierarchy. These results arise from differing approaches to incorporating prior knowledge about neutrinos. We develop a hierarchy-agnostic prior and show that the hierarchy cannot be conclusively determined with current data. The determination of the h… ▽ More Recent works combining neutrino oscillation and cosmological data to determine the neutrino hierarchy found a range of odds in favour of the normal hierarchy. These results arise from differing approaches to incorporating prior knowledge about neutrinos. We develop a hierarchy-agnostic prior and show that the hierarchy cannot be conclusively determined with current data. The determination of the hierarchy is limited by the neutrino mass scale $Σ_ν$ measurement. We obtain a target precision of $σ(Σ_ν) = 0.014$ eV, necessary for conclusively establishing the normal hierarchy with future data. △ Less

Submitted 18 March, 2020; v1 submitted 9 July, 2019; originally announced July 2019.

Comments: v3: Version accepted in Phys Rev D. Minor clarifications and added robustness tests, conclusions unchanged

Journal ref: Phys. Rev. D 101, 083513 (2020)

An interpretable machine learning framework for dark matter halo formation

Authors: Luisa Lucie-Smith, Hiranya V. Peiris, Andrew Pontzen

Abstract: We present a generalization of our recently proposed machine learning framework, aiming to provide new physical insights into dark matter halo formation. We investigate the impact of the initial density and tidal shear fields on the formation of haloes over the mass range $11.4 \leq \log(M/M_{\odot}) \leq 13.4$. The algorithm is trained on an N-body simulation to infer the final mass of the halo t… ▽ More We present a generalization of our recently proposed machine learning framework, aiming to provide new physical insights into dark matter halo formation. We investigate the impact of the initial density and tidal shear fields on the formation of haloes over the mass range $11.4 \leq \log(M/M_{\odot}) \leq 13.4$. The algorithm is trained on an N-body simulation to infer the final mass of the halo to which each dark matter particle will later belong. We then quantify the difference in the predictive accuracy between machine learning models using a metric based on the Kullback-Leibler divergence. We first train the algorithm with information about the density contrast in the particles' local environment. The addition of tidal shear information does not yield an improved halo collapse model over one based on density information alone; the difference in their predictive performance is consistent with the statistical uncertainty of the density-only based model. This implies that our machine learning setup does not identify any significant role for the tidal shear in determining halo masses. This result is confirmed as we verify the ability of the initial conditions-to-halo mass mapping learnt from one simulation to generalize to independent simulations. Our work illustrates the broader potential of developing interpretable machine learning frameworks to gain physical understanding of non-linear large-scale structure formation. △ Less

Submitted 19 September, 2019; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: 13 pages, 8 figures. Minor changes to match version published in MNRAS. Accepted on 12/09/2019

Nonlinear Dynamics of the Cold Atom Analog False Vacuum

Authors: Jonathan Braden, Matthew C. Johnson, Hiranya V. Peiris, Andrew Pontzen, Silke Weinfurtner

Abstract: We investigate the nonlinear dynamics of cold atom systems that can in principle serve as quantum simulators of false vacuum decay. The analog false vacuum manifests as a metastable vacuum state for the relative phase in a two-species Bose-Einstein condensate (BEC), induced by a driven periodic coupling between the two species. In the appropriate low energy limit, the evolution of the relative pha… ▽ More We investigate the nonlinear dynamics of cold atom systems that can in principle serve as quantum simulators of false vacuum decay. The analog false vacuum manifests as a metastable vacuum state for the relative phase in a two-species Bose-Einstein condensate (BEC), induced by a driven periodic coupling between the two species. In the appropriate low energy limit, the evolution of the relative phase is approximately governed by a relativistic wave equation exhibiting true and false vacuum configurations. In previous work, a linear stability analysis identified exponentially growing short-wavelength modes driven by the time-dependent coupling. These modes threaten to destabilize the analog false vacuum. Here, we employ numerical simulations of the coupled Gross-Pitaevski equations (GPEs) to determine the non-linear evolution of these linearly unstable modes. We find that unless a physical mechanism modifies the GPE on short length scales, the analog false vacuum is indeed destabilized. We briefly discuss various physically expected corrections to the GPEs that may act to remove the exponentially unstable modes. To investigate the resulting dynamics in cases where such a removal mechanism exists, we implement a hard UV cutoff that excludes the unstable modes as a simple model for these corrections. We use this to study the range of phenomena arising from such a system. In particular, we show that by modulating the strength of the time-dependent coupling, it is possible to observe the crossover between a second and first order phase transition out of the false vacuum. △ Less

Submitted 9 October, 2019; v1 submitted 16 April, 2019; originally announced April 2019.

Comments: v2: Version accepted to JHEP. Explanatory material added. Results unchanged. v1: 21 pages+appendices, 15 figures

Journal ref: JHEP 10 (2019) 174

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

Primordial Non-Gaussianity

Authors: P. Daniel Meerburg, Daniel Green, Muntazir Abidi, Mustafa A. Amin, Peter Adshead, Zeeshan Ahmed, David Alonso, Behzad Ansarinejad, Robert Armstrong, Santiago Avila, Carlo Baccigalupi, Tobias Baldauf, Mario Ballardini, Kevin Bandura, Nicola Bartolo, Nicholas Battaglia, Daniel Baumann, Chetan Bavdhankar, José Luis Bernal, Florian Beutler, Matteo Biagetti, Colin Bischoff, Jonathan Blazek, J. Richard Bond, Julian Borrill , et al. (153 additional authors not shown)

Abstract: Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with… ▽ More Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale. △ Less

Submitted 14 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

First measurement of the Hubble constant from a dark standard siren using the Dark Energy Survey galaxies and the LIGO/Virgo binary-black-hole merger GW170814

Authors: The DES Collaboration, the LIGO Scientific Collaboration, the Virgo Collaboration, M. Soares-Santos, A. Palmese, W. Hartley, J. Annis, J. Garcia-Bellido, O. Lahav, Z. Doctor, M. Fishbach, D. E. Holz, H. Lin, M. E. S. Pereira, A. Garcia, K. Herner, R. Kessler, H. V. Peiris, M. Sako, S. Allam, D. Brout, A. Carnero Rosell, H. Y. Chen, C. Conselice, J. deRose , et al. (1181 additional authors not shown)

Abstract: We present a multi-messenger measurement of the Hubble constant H_0 using the binary-black-hole merger GW170814 as a standard siren, combined with a photometric redshift catalog from the Dark Energy Survey (DES). The luminosity distance is obtained from the gravitational wave signal detected by the LIGO/Virgo Collaboration (LVC) on 2017 August 14, and the redshift information is provided by the DE… ▽ More We present a multi-messenger measurement of the Hubble constant H_0 using the binary-black-hole merger GW170814 as a standard siren, combined with a photometric redshift catalog from the Dark Energy Survey (DES). The luminosity distance is obtained from the gravitational wave signal detected by the LIGO/Virgo Collaboration (LVC) on 2017 August 14, and the redshift information is provided by the DES Year 3 data. Black-hole mergers such as GW170814 are expected to lack bright electromagnetic emission to uniquely identify their host galaxies and build an object-by-object Hubble diagram. However, they are suitable for a statistical measurement, provided that a galaxy catalog of adequate depth and redshift completion is available. Here we present the first Hubble parameter measurement using a black-hole merger. Our analysis results in $H_0 = 75.2^{+39.5}_{-32.4}~{\rm km~s^{-1}~Mpc^{-1}}$, which is consistent with both SN Ia and CMB measurements of the Hubble constant. The quoted 68% credible region comprises 60% of the uniform prior range [20,140] ${\rm km~s^{-1}~Mpc^{-1}}$, and it depends on the assumed prior range. If we take a broader prior of [10,220] ${\rm km~s^{-1}~Mpc^{-1}}$, we find $H_0 = 78^{+ 96}_{-24}~{\rm km~s^{-1}~Mpc^{-1}}$ ($57\%$ of the prior range). Although a weak constraint on the Hubble constant from a single event is expected using the dark siren method, a multifold increase in the LVC event rate is anticipated in the coming years and combinations of many sirens will lead to improved constraints on $H_0$. △ Less

Submitted 22 March, 2019; v1 submitted 6 January, 2019; originally announced January 2019.

Comments: 18 pages, 2 Figures, submitted to ApJL

Report number: FERMILAB-PUB-18-629-AE

Journal ref: The Astrophysical Journal Letters, Volume 876, 1, L7 (2019)

Serendipitous Discoveries of Kilonovae in the LSST Main Survey: Maximising Detections of Sub-Threshold Gravitational Wave Events

Authors: Christian N. Setzer, Rahul Biswas, Hiranya V. Peiris, Stephan Rosswog, Oleg Korobkin, Ryan T. Wollaeger

Abstract: We investigate the ability of the Large Synoptic Survey Telescope (LSST) to discover kilonovae (kNe) from binary neutron star (BNS) and neutron star-black hole (NSBH) mergers, focusing on serendipitous detections in the Wide-Fast-Deep (WFD) survey. We simulate observations of kNe with proposed LSST survey strategies, paying particular attention to cadence choices that are compatible with the broad… ▽ More We investigate the ability of the Large Synoptic Survey Telescope (LSST) to discover kilonovae (kNe) from binary neutron star (BNS) and neutron star-black hole (NSBH) mergers, focusing on serendipitous detections in the Wide-Fast-Deep (WFD) survey. We simulate observations of kNe with proposed LSST survey strategies, paying particular attention to cadence choices that are compatible with the broader LSST cosmology programme. We find that if all kNe are identical to GW170817, the baseline survey strategy will yield 58 kNe over the survey lifetime. If we instead assume a representative population model of BNS kNe, we expect to detect only 27 kNe. However, we find the choice of survey strategy significantly impacts these numbers and can increase them to 254 kNe and 82 kNe over the survey lifetime, respectively. This improvement arises from an increased cadence of observations between different filters with respect to the baseline. We then consider the ability of the Advanced LIGO/Virgo (ALV) detector network to detect these BNS mergers. If the optimal survey strategy is adopted, 202 of the GW170817-like kNe and 56 of the BNS population model kNe are detected with LSST but are below the threshold for detection by the ALV network. This represents, for both models, an increase by a factor greater than 4.5 in the number of detected sub-threshold events over the baseline survey strategy. Such a population of sub-threshold events would provide an opportunity to conduct electromagnetic-triggered searches for signals in gravitational-wave detector data and assess selection effects in measurements of the Hubble constant from standard sirens, e.g., related to viewing angle effects. △ Less

Submitted 18 February, 2019; v1 submitted 26 December, 2018; originally announced December 2018.

Comments: 14 pages, 9 figures, 6 tables, to be submitted to MNRAS; v2: version accepted by MNRAS. Minor clarifications added

Report number: LA-UR-18-31553

An Emulator for the Lyman-alpha Forest

Authors: Simeon Bird, Keir K. Rogers, Hiranya V. Peiris, Licia Verde, Andreu Font-Ribera, Andrew Pontzen

Abstract: We present methods for interpolating between the 1-D flux power spectrum of the Lyman-$α$ forest, as output by cosmological hydrodynamic simulations. Interpolation is necessary for cosmological parameter estimation due to the limited number of simulations possible. We construct an emulator for the Lyman-$α$ forest flux power spectrum from $21$ small simulations using Latin hypercube sampling and G… ▽ More We present methods for interpolating between the 1-D flux power spectrum of the Lyman-$α$ forest, as output by cosmological hydrodynamic simulations. Interpolation is necessary for cosmological parameter estimation due to the limited number of simulations possible. We construct an emulator for the Lyman-$α$ forest flux power spectrum from $21$ small simulations using Latin hypercube sampling and Gaussian process interpolation. We show that this emulator has a typical accuracy of 1.5% and a worst-case accuracy of 4%, which compares well to the current statistical error of 3 - 5% at $z < 3$ from BOSS DR9. We compare to the previous state of the art, quadratic polynomial interpolation. The Latin hypercube samples the entire volume of parameter space, while quadratic polynomial emulation samples only lower-dimensional subspaces. The Gaussian process provides an estimate of the emulation error and we show using test simulations that this estimate is reasonable. We construct a likelihood function and use it to show that the posterior constraints generated using the emulator are unbiased. We show that our Gaussian process emulator has lower emulation error than quadratic polynomial interpolation and thus produces tighter posterior confidence intervals, which will be essential for future Lyman-$α$ surveys such as DESI. △ Less

Submitted 23 January, 2023; v1 submitted 11 December, 2018; originally announced December 2018.

Comments: 28 pages, 10 figures, accepted to JCAP with minor changes. Fixes a typographical error in eq. 2.3

Journal ref: JCAP 02 (2019) 050

Bayesian emulator optimisation for cosmology: application to the Lyman-alpha forest

Authors: Keir K. Rogers, Hiranya V. Peiris, Andrew Pontzen, Simeon Bird, Licia Verde, Andreu Font-Ribera

Abstract: The Lyman-alpha forest provides strong constraints on both cosmological parameters and intergalactic medium astrophysics, which are forecast to improve further with the next generation of surveys including eBOSS and DESI. As is generic in cosmological inference, extracting this information requires a likelihood to be computed throughout a high-dimensional parameter space. Evaluating the likelihood… ▽ More The Lyman-alpha forest provides strong constraints on both cosmological parameters and intergalactic medium astrophysics, which are forecast to improve further with the next generation of surveys including eBOSS and DESI. As is generic in cosmological inference, extracting this information requires a likelihood to be computed throughout a high-dimensional parameter space. Evaluating the likelihood requires a robust and accurate mapping between the parameters and observables, in this case the 1D flux power spectrum. Cosmological simulations enable such a mapping, but due to computational time constraints can only be evaluated at a handful of sample points; "emulators" are designed to interpolate between these. The problem then reduces to placing the sample points such that an accurate mapping is obtained while minimising the number of expensive simulations required. To address this, we introduce an emulation procedure that employs Bayesian optimisation of the training set for a Gaussian process interpolation scheme. Starting with a Latin hypercube sampling (other schemes with good space-filling properties can be used), we iteratively augment the training set with extra simulations at new parameter positions which balance the need to reduce interpolation error while focussing on regions of high likelihood. We show that smaller emulator error from the Bayesian optimisation propagates to smaller widths on the posterior distribution. Even with fewer simulations than a Latin hypercube, Bayesian optimisation shrinks the 95% credible volume by 90% and, e.g., the 1 sigma error on the amplitude of small-scale primordial fluctuations by 38%. This is the first demonstration of Bayesian optimisation applied to large-scale structure emulation, and we anticipate the technique will generalise to many other probes such as galaxy clustering, weak lensing and 21cm. △ Less

Submitted 5 February, 2019; v1 submitted 11 December, 2018; originally announced December 2018.

Comments: 23 pages, 4 figures. Minor changes to match version published in JCAP

Journal ref: JCAP02(2019)031

Optimizing the LSST Observing Strategy for Dark Energy Science: DESC Recommendations for the Wide-Fast-Deep Survey

Authors: Michelle Lochner, Daniel M. Scolnic, Humna Awan, Nicolas Regnault, Philippe Gris, Rachel Mandelbaum, Eric Gawiser, Husni Almoubayyed, Christian N. Setzer, Simon Huber, Melissa L. Graham, Renée Hložek, Rahul Biswas, Tim Eifler, Daniel Rothchild, Tarek Allam Jr, Jonathan Blazek, Chihway Chang, Thomas Collett, Ariel Goobar, Isobel M. Hook, Mike Jarvis, Saurabh W. Jha, Alex G. Kim, Phil Marshall , et al. (11 additional authors not shown)

Abstract: Cosmology is one of the four science pillars of LSST, which promises to be transformative for our understanding of dark energy and dark matter. The LSST Dark Energy Science Collaboration (DESC) has been tasked with deriving constraints on cosmological parameters from LSST data. Each of the cosmological probes for LSST is heavily impacted by the choice of observing strategy. This white paper is wri… ▽ More Cosmology is one of the four science pillars of LSST, which promises to be transformative for our understanding of dark energy and dark matter. The LSST Dark Energy Science Collaboration (DESC) has been tasked with deriving constraints on cosmological parameters from LSST data. Each of the cosmological probes for LSST is heavily impacted by the choice of observing strategy. This white paper is written by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy that will benefit all cosmological analyses with LSST. It is accompanied by the DESC DDF (Deep Drilling Fields) white paper (Scolnic et al.). We use a variety of metrics to understand the effects of the observing strategy on measurements of weak lensing, large-scale structure, clusters, photometric redshifts, supernovae, strong lensing and kilonovae. In order to reduce systematic uncertainties, we conclude that the current baseline observing strategy needs to be significantly modified to result in the best possible cosmological constraints. We provide some key recommendations: moving the WFD (Wide-Fast-Deep) footprint to avoid regions of high extinction, taking visit pairs in different filters, changing the 2x15s snaps to a single exposure to improve efficiency, focusing on strategies that reduce long gaps (>15 days) between observations, and prioritizing spatial uniformity at several intervals during the 10-year survey. △ Less

Submitted 14 December, 2018; v1 submitted 30 November, 2018; originally announced December 2018.

Comments: The LSST DESC response (WFD) to the Call for White Papers on LSST Cadence Optimization. Comments welcome

Unbiased Hubble constant estimation from binary neutron star mergers

Authors: Daniel J. Mortlock, Stephen M. Feeney, Hiranya V. Peiris, Andrew R. Williamson, Samaya M. Nissanke

Abstract: Gravitational wave (GW) observations of binary neutron star (BNS) mergers can be used to measure luminosity distances and hence, when coupled with estimates for the mergers' host redshifts, infer the Hubble constant, $H_0$. These observations are, however, affected by GW measurement noise, uncertainties in host redshifts and peculiar velocities, and are potentially biased by selection effects and… ▽ More Gravitational wave (GW) observations of binary neutron star (BNS) mergers can be used to measure luminosity distances and hence, when coupled with estimates for the mergers' host redshifts, infer the Hubble constant, $H_0$. These observations are, however, affected by GW measurement noise, uncertainties in host redshifts and peculiar velocities, and are potentially biased by selection effects and the mis-specification of the cosmological model or the BNS population. The estimation of $H_0$ from samples of BNS mergers with optical counterparts is tested here by using a phenomenological model for the GW strains that captures both the data-driven event selection and the distance-inclination degeneracy, while being simple enough to facilitate large numbers of simulations. A rigorous Bayesian approach to analyzing the data from such simulated BNS merger samples is shown to yield results that are unbiased, have the appropriate uncertainties, and are robust to model mis-specification. Applying such methods to a sample of $N \simeq 50$ BNS merger events, as LIGO+Virgo could produce in the next $\sim 5$ years, should yield robust and accurate Hubble constant estimates that are precise to a level of $\sim 2$ km s$^{-1}$ Mpc$^{-1}$, sufficient to reliably resolve the current tension between local and cosmological measurements of $H_0$. △ Less

Submitted 7 October, 2019; v1 submitted 28 November, 2018; originally announced November 2018.

Comments: 21 pages (including eight-page appendices); 9 figures. Matches version accepted by PRD

Journal ref: Phys. Rev. D 100, 103523 (2019)

Classification of Multiwavelength Transients with Machine Learning

Authors: K. Sooknunan, M. Lochner, Bruce A. Bassett, H. V. Peiris, R. Fender, A. J. Stewart, M. Pietka, P. A. Woudt, J. D. McEwen, O. Lahav

Abstract: With the advent of powerful telescopes such as the Square Kilometer Array and the Vera C. Rubin Observatory, we are entering an era of multiwavelength transient astronomy that will lead to a dramatic increase in data volume. Machine learning techniques are well suited to address this data challenge and rapidly classify newly detected transients. We present a multiwavelength classification algorith… ▽ More With the advent of powerful telescopes such as the Square Kilometer Array and the Vera C. Rubin Observatory, we are entering an era of multiwavelength transient astronomy that will lead to a dramatic increase in data volume. Machine learning techniques are well suited to address this data challenge and rapidly classify newly detected transients. We present a multiwavelength classification algorithm consisting of three steps: (1) interpolation and augmentation of the data using Gaussian processes; (2) feature extraction using wavelets; and (3) classification with random forests. Augmentation provides improved performance at test time by balancing the classes and adding diversity into the training set. In the first application of machine learning to the classification of real radio transient data, we apply our technique to the Green Bank Interferometer and other radio light curves. We find we are able to accurately classify most of the 11 classes of radio variables and transients after just eight hours of observations, achieving an overall test accuracy of 78 percent. We fully investigate the impact of the small sample size of 82 publicly available light curves and use data augmentation techniques to mitigate the effect. We also show that on a significantly larger simulated representative training set that the algorithm achieves an overall accuracy of 97 percent, illustrating that the method is likely to provide excellent performance on future surveys. Finally, we demonstrate the effectiveness of simultaneous multiwavelength observations by showing how incorporating just one optical data point into the analysis improves the accuracy of the worst performing class by 19 percent. △ Less

Submitted 8 March, 2021; v1 submitted 20 November, 2018; originally announced November 2018.

Comments: 16 pages, 12 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 502, Issue 1, March 2021, Pages 206-224

Cosmological Constraints from Multiple Probes in the Dark Energy Survey

Authors: DES Collaboration, T. M. C. Abbott, A. Alarcon, S. Allam, P. Andersen, F. Andrade-Oliveira, J. Annis, J. Asorey, A. Avelino, S. Avila, D. Bacon, N. Banik, B. A. Bassett, E. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, E. Bertin, J. Blazek, S. L. Bridle, D. Brooks, D. Brout, D. L. Burke, J. Calcino, H. Camacho , et al. (144 additional authors not shown)

Abstract: The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these pr… ▽ More The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $Ω_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade. △ Less

Submitted 6 May, 2019; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: 8 pages, 2 figures; v3 matches version accepted by PRL

Report number: FERMILAB-PUB-18-585-AE

Journal ref: Phys. Rev. Lett. 122, 171301 (2019)

Inferring properties of the local white dwarf population in astrometric and photometric surveys

Authors: Axel Widmark, Daniel J. Mortlock, Hiranya V. Peiris

Abstract: The Gaia mission will provide precise astrometry for an unprecedented number of white dwarfs (WDs), encoding information on stellar evolution, Type Ia supernovae progenitor scenarios, and the star formation and dynamical history of the Milky Way. With such a large data set, it is possible to infer properties of the WD population using only astrometric and photometric information. We demonstrate a… ▽ More The Gaia mission will provide precise astrometry for an unprecedented number of white dwarfs (WDs), encoding information on stellar evolution, Type Ia supernovae progenitor scenarios, and the star formation and dynamical history of the Milky Way. With such a large data set, it is possible to infer properties of the WD population using only astrometric and photometric information. We demonstrate a framework to accomplish this using a mock data set with SDSS ugriz photometry and Gaia astrometric information. Our technique utilises a Bayesian hierarchical model for inferring properties of a WD population while also taking into account all observational errors of individual objects, as well as selection and incompleteness effects. We demonstrate that photometry alone can constrain the WD population's distributions of temperature, surface gravity and phenomenological type, and that astrometric information significantly improves determination of the WD surface gravity distribution. We also discuss the possibility of identifying unresolved binary WDs using only photometric and astrometric information. △ Less

Submitted 1 August, 2018; originally announced August 2018.

Comments: 11 pages, 9 figures

Journal ref: MNRAS 485 (2019) 1

Planck 2018 results. X. Constraints on inflation

Authors: Planck Collaboration, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso , et al. (151 additional authors not shown)

Abstract: We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index… ▽ More We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (φ) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation. △ Less

Submitted 2 August, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: References added and minor improvements. BICEP2/Keck Array BK15 is used in the place of BICEP2/Keck Array BK14

Journal ref: A&A 641, A10 (2020)

Planck 2018 results. VI. Cosmological parameters

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, R. Battye, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese , et al. (157 additional authors not shown)

Abstract: We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $Λ$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $Λ$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis g… ▽ More We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $Λ$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $Λ$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $Ω_c h^2 = 0.120\pm 0.001$, baryon density $Ω_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $τ= 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100θ_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$Λ$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $Ω_m = 0.315\pm 0.007$; and matter fluctuation amplitude $σ_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$Λ$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_ν< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$Λ$CDM at over $2\,σ$, which pulls some parameters that affect the lensing amplitude away from the base-$Λ$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged) △ Less

Submitted 9 August, 2021; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 73 pages; Updated with published reionization result corrigendum on p59. Parameter tables and chains available at https://wiki.cosmos.esa.int/planck-legacy-archive/index.php/Cosmological_Parameters

Journal ref: A&A 641, A6 (2020)