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

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

pop-cosmos: A comprehensive picture of the galaxy population from COSMOS data

Authors: Justin Alsing, Stephen Thorp, Sinan Deger, Hiranya Peiris, Boris Leistedt, Daniel Mortlock, Joel Leja

Abstract: We present pop-cosmos: a comprehensive model characterizing the galaxy population, calibrated to $140,938$ galaxies from the Cosmic Evolution Survey (COSMOS) with photometry in $26$ bands from the ultra-violet to the infra-red. We construct a detailed forward model for the COSMOS data, comprising: a population model describing the joint distribution of galaxy characteristics and its evolution (par… ▽ More We present pop-cosmos: a comprehensive model characterizing the galaxy population, calibrated to $140,938$ galaxies from the Cosmic Evolution Survey (COSMOS) with photometry in $26$ bands from the ultra-violet to the infra-red. We construct a detailed forward model for the COSMOS data, comprising: a population model describing the joint distribution of galaxy characteristics and its evolution (parameterized by a flexible score-based diffusion model); a state-of-the-art stellar population synthesis (SPS) model connecting galaxies' instrinsic properties to their photometry; and a data-model for the observation, calibration and selection processes. By minimizing the optimal transport distance between synthetic and real data we are able to jointly fit the population- and data-models, leading to robustly calibrated population-level inferences that account for parameter degeneracies, photometric noise and calibration, and selection effects. We present a number of key predictions from our model of interest for cosmology and galaxy evolution, including the mass function and redshift distribution; the mass-metallicity-redshift and fundamental metallicity relations; the star-forming sequence; the relation between dust attenuation and stellar mass, star formation rate and attenuation-law index; and the relation between gas-ionization and star formation. Our model encodes a comprehensive picture of galaxy evolution that faithfully predicts galaxy colors across a broad redshift ($z<4$) and wavelength range. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 27 pages, 13 figures. See also companion paper by Thorp et. al. in the same arXiv mailing

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)

Hierarchical Bayesian inference of photometric redshifts with stellar population synthesis models

Authors: Boris Leistedt, Justin Alsing, Hiranya Peiris, Daniel Mortlock, Joel Leja

Abstract: We present a Bayesian hierarchical framework to analyze photometric galaxy survey data with stellar population synthesis (SPS) models. Our method couples robust modeling of spectral energy distributions with a population model and a noise model to characterize the statistical properties of the galaxy populations and real observations, respectively. By self-consistently inferring all model paramete… ▽ More We present a Bayesian hierarchical framework to analyze photometric galaxy survey data with stellar population synthesis (SPS) models. Our method couples robust modeling of spectral energy distributions with a population model and a noise model to characterize the statistical properties of the galaxy populations and real observations, respectively. By self-consistently inferring all model parameters, from high-level hyper-parameters to SPS parameters of individual galaxies, one can separate sources of bias and uncertainty in the data.We demonstrate the strengths and flexibility of this approach by deriving accurate photometric redshifts for a sample of spectroscopically-confirmed galaxies in the COSMOS field, all with 26-band photometry and spectroscopic redshifts. We achieve a performance competitive with publicly-released photometric redshift catalogs based on the same data. Prior to this work, this approach was computationally intractable in practice due to the heavy computational load of SPS model calls; we overcome this challenge using with neural emulators. We find that the largest photometric residuals are associated with poor calibration for emission line luminosities and thus build a framework to mitigate these effects. This combination of physics-based modeling accelerated with machine learning paves the path towards meeting the stringent requirements on the accuracy of photometric redshift estimation imposed by upcoming cosmological surveys. The approach also has the potential to create new links between cosmology and galaxy evolution through the analysis of photometric datasets. △ Less

Submitted 4 November, 2022; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 16 pages, 6 figures. Version accepted in APJS

Forward modeling of galaxy populations for cosmological redshift distribution inference

Authors: Justin Alsing, Hiranya Peiris, Daniel Mortlock, Joel Leja, Boris Leistedt

Abstract: We present a forward modeling framework for estimating galaxy redshift distributions from photometric surveys. Our forward model is composed of: a detailed population model describing the intrinsic distribution of physical characteristics of galaxies, encoding galaxy evolution physics; a stellar population synthesis model connecting the physical properties of galaxies to their photometry; a data-m… ▽ More We present a forward modeling framework for estimating galaxy redshift distributions from photometric surveys. Our forward model is composed of: a detailed population model describing the intrinsic distribution of physical characteristics of galaxies, encoding galaxy evolution physics; a stellar population synthesis model connecting the physical properties of galaxies to their photometry; a data-model characterizing the observation and calibration processes for a given survey; and, explicit treatment of selection cuts, both into the main analysis sample and subsequent sorting into tomographic redshift bins. This approach has the appeal that it does not rely on spectroscopic calibration data, provides explicit control over modeling assumptions, and builds a direct bridge between photo-$z$ inference and galaxy evolution physics. In addition to redshift distributions, forward modeling provides a framework for drawing robust inferences about the statistical properties of the galaxy population more generally. We demonstrate the utility of forward modeling by estimating the redshift distributions for the Galaxy And Mass Assembly (GAMA) and Vimos VLT Deep (VVDS) surveys, validating against their spectroscopic redshifts. Our baseline model is able to predict tomographic redshift distributions for GAMA and VVDS with a bias of $Δでるたz \lesssim 0.003$ and $Δでるたz \simeq 0.01$ on the mean redshift respectively -- comfortably accurate enough for Stage III cosmological surveys -- without any hyper-parameter tuning (i.e., prior to doing any fitting to those data). We anticipate that with additional hyper-parameter fitting and modeling improvements, forward modeling can provide a path to accurate redshift distribution inference for Stage IV surveys. △ Less

Submitted 17 November, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: Accepted to ApJS 18 September 2022

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

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: a pioneering process of community-focused experimental design

Authors: Federica B. Bianco, Željko Ivezić, R. Lynne Jones, Melissa L. Graham, Phil Marshall, Abhijit Saha, Michael A. Strauss, Peter Yoachim, Tiago Ribeiro, Timo Anguita, Franz E. Bauer, Eric C. Bellm, Robert D. Blum, William N. Brandt, Sarah Brough, Màrcio Catelan, William I. Clarkson, Andrew J. Connolly, Eric Gawiser, John Gizis, Renee Hlozek, Sugata Kaviraj, Charles T. Liu, Michelle Lochner, Ashish A. Mahabal , et al. (21 additional authors not shown)

Abstract: Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core scienc… ▽ More Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue. △ Less

Submitted 1 September, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: Submitted as the opening paper of the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

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

The Impact of Observing Strategy on Cosmological Constraints with LSST

Authors: Michelle Lochner, Dan Scolnic, Husni Almoubayyed, Timo Anguita, Humna Awan, Eric Gawiser, Satya Gontcho A Gontcho, Philippe Gris, Simon Huber, Saurabh W. Jha, R. Lynne Jones, Alex G. Kim, Rachel Mandelbaum, Phil Marshall, Tanja Petrushevska, Nicolas Regnault, Christian N. Setzer, Sherry H. Suyu, Peter Yoachim, Rahul Biswas, Tristan Blaineau, Isobel Hook, Marc Moniez, Eric Neilsen, Hiranya Peiris , et al. (2 additional authors not shown)

Abstract: The generation-defining Vera C. Rubin Observatory will make state-of-the-art measurements of both the static and transient universe through its Legacy Survey for Space and Time (LSST). With such capabilities, it is immensely challenging to optimize the LSST observing strategy across the survey's wide range of science drivers. Many aspects of the LSST observing strategy relevant to the LSST Dark En… ▽ More The generation-defining Vera C. Rubin Observatory will make state-of-the-art measurements of both the static and transient universe through its Legacy Survey for Space and Time (LSST). With such capabilities, it is immensely challenging to optimize the LSST observing strategy across the survey's wide range of science drivers. Many aspects of the LSST observing strategy relevant to the LSST Dark Energy Science Collaboration, such as survey footprint definition, single visit exposure time and the cadence of repeat visits in different filters, are yet to be finalized. Here, we present metrics used to assess the impact of observing strategy on the cosmological probes considered most sensitive to survey design; these are large-scale structure, weak lensing, type Ia supernovae, kilonovae and strong lens systems (as well as photometric redshifts, which enable many of these probes). We evaluate these metrics for over 100 different simulated potential survey designs. Our results show that multiple observing strategy decisions can profoundly impact cosmological constraints with LSST; these include adjusting the survey footprint, ensuring repeat nightly visits are taken in different filters and enforcing regular cadence. We provide public code for our metrics, which makes them readily available for evaluating further modifications to the survey design. We conclude with a set of recommendations and highlight observing strategy factors that require further research. △ Less

Submitted 12 April, 2021; originally announced April 2021.

Comments: 49 pages, 17 figures

Journal ref: ApJS 259 58, 2022

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)

GenetIC -- a new initial conditions generator to support genetically modified zoom simulations

Authors: Stephen Stopyra, Andrew Pontzen, Hiranya Peiris, Nina Roth, Martin Rey

Abstract: We present genetIC, a new code for generating initial conditions for cosmological N-body simulations. The code allows precise, user-specified alterations to be made to arbitrary regions of the simulation (while maintaining consistency with the statistical ensemble). These "genetic modifications" allow, for example, the history, mass, or environment of a target halo to be altered in order to study… ▽ More We present genetIC, a new code for generating initial conditions for cosmological N-body simulations. The code allows precise, user-specified alterations to be made to arbitrary regions of the simulation (while maintaining consistency with the statistical ensemble). These "genetic modifications" allow, for example, the history, mass, or environment of a target halo to be altered in order to study the effect on their evolution. The code natively supports initial conditions with nested zoom regions at progressively increasing resolution. Modifications in the high-resolution region must propagate self-consistently onto the lower resolution grids; to enable this while maintaining a small memory footprint, we introduce a Fourier-space filtering approach to generating fields at variable resolution. Due to a close correspondence with modifications, constrained initial conditions can also be produced by genetIC (for example with the aim of matching structures in the local Universe). We test the accuracy of modifications performed within zoom initial conditions. The code achieves sub-percent precision, which is easily sufficient for current applications in galaxy formation. △ Less

Submitted 1 December, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: 23 pages, 8 figures, github repository: https://github.com/pynbody/genetIC, Accepted by ApJS 23/11/20, changed to match accepted version

SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry

Authors: Justin Alsing, Hiranya Peiris, Joel Leja, ChangHoon Hahn, Rita Tojeiro, Daniel Mortlock, Boris Leistedt, Benjamin D. Johnson, Charlie Conroy

Abstract: We present SPECULATOR - a fast, accurate, and flexible framework for emulating stellar population synthesis (SPS) models for predicting galaxy spectra and photometry. For emulating spectra, we use principal component analysis to construct a set of basis functions, and neural networks to learn the basis coefficients as a function of the SPS model parameters. For photometry, we parameterize the magn… ▽ More We present SPECULATOR - a fast, accurate, and flexible framework for emulating stellar population synthesis (SPS) models for predicting galaxy spectra and photometry. For emulating spectra, we use principal component analysis to construct a set of basis functions, and neural networks to learn the basis coefficients as a function of the SPS model parameters. For photometry, we parameterize the magnitudes (for the filters of interest) as a function of SPS parameters by a neural network. The resulting emulators are able to predict spectra and photometry under both simple and complicated SPS model parameterizations to percent-level accuracy, giving a factor of $10^3$-$10^4$ speed up over direct SPS computation. They have readily-computable derivatives, making them amenable to gradient-based inference and optimization methods. The emulators are also straightforward to call from a GPU, giving an additional order-of-magnitude speed-up. Rapid SPS computations delivered by emulation offers a massive reduction in the computational resources required to infer the physical properties of galaxies from observed spectra or photometry and simulate galaxy populations under SPS models, whilst maintaining the accuracy required for a range of applications. △ Less

Submitted 15 April, 2020; v1 submitted 26 November, 2019; originally announced November 2019.

Comments: 15 pages, 9 figures, accepted by ApJS April 2020

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