-
Baryonification extended to thermal Sunyaev Zeldovich
Authors:
Giovanni Aricò,
Raul E. Angulo
Abstract:
Baryonification algorithms model the impact of galaxy formation and feedback on the matter field in gravity-only simulations by adopting physically motivated parametric prescriptions. In this paper, we extend these models to describe gas temperature and pressure, allowing for a self-consistent modelling of the thermal Sunyaev-Zeldovich effect, weak gravitational lensing, and their cross-correlatio…
▽ More
Baryonification algorithms model the impact of galaxy formation and feedback on the matter field in gravity-only simulations by adopting physically motivated parametric prescriptions. In this paper, we extend these models to describe gas temperature and pressure, allowing for a self-consistent modelling of the thermal Sunyaev-Zeldovich effect, weak gravitational lensing, and their cross-correlation, down to small scales. We validate our approach by showing that it can simultaneously reproduce the electron pressure, gas, stellar, and dark matter power spectra as measured in all BAHAMAS hydrodynamical simulations. Specifically, with only two additional free parameters, we can fit the electron pressure auto- and cross-power spectra at 10% while reproducing the suppression in the matter power spectrum induced by baryons at the per cent level, for different AGN feedback strengths in BAHAMAS. Furthermore, we reproduce BAHAMAS convergence and thermal Sunyaev Zeldovich angular power spectra within 1% and 10% accuracy, respectively, down to l = 5000. When used jointly with cosmological rescaling algorithms, the baryonification presented here will allow a fast and accurate exploration of cosmological and astrophysical scenarios. Therefore, it can be employed to create mock catalogues, lightcones, and large training sets for emulators aimed at interpreting forthcoming multi-wavelength observations of the large-scale structure of the Universe.
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
Euclid. I. Overview of the Euclid mission
Authors:
Euclid Collaboration,
Y. Mellier,
Abdurro'uf,
J. A. Acevedo Barroso,
A. Achúcarro,
J. Adamek,
R. Adam,
G. E. Addison,
N. Aghanim,
M. Aguena,
V. Ajani,
Y. Akrami,
A. Al-Bahlawan,
A. Alavi,
I. S. Albuquerque,
G. Alestas,
G. Alguero,
A. Allaoui,
S. W. Allen,
V. Allevato,
A. V. Alonso-Tetilla,
B. Altieri,
A. Alvarez-Candal,
A. Amara,
L. Amendola
, et al. (1086 additional authors not shown)
Abstract:
The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14…
▽ More
The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.
△ Less
Submitted 22 May, 2024;
originally announced May 2024.
-
Cosmic shear with small scales: DES-Y3, KiDS-1000 and HSC-DR1
Authors:
Carlos García-García,
Matteo Zennaro,
Giovanni Aricò,
David Alonso,
Raul E. Angulo
Abstract:
We present a cosmological analysis of the combination of the DES-Y3, KiDS-1000 and HSC-DR1 weak lensing samples under a joint harmonic-space pipeline making use of angular scales down to $\ell_{\rm max}=4500$, corresponding to significantly smaller scales ($δθ\sim2.4'$) than those commonly used in cosmological weak lensing studies. We are able to do so by accurately modelling non-linearities and t…
▽ More
We present a cosmological analysis of the combination of the DES-Y3, KiDS-1000 and HSC-DR1 weak lensing samples under a joint harmonic-space pipeline making use of angular scales down to $\ell_{\rm max}=4500$, corresponding to significantly smaller scales ($δθ\sim2.4'$) than those commonly used in cosmological weak lensing studies. We are able to do so by accurately modelling non-linearities and the impact of baryonic effects using Baccoemu. We find $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3}=0.795^{+0.015}_{-0.017}$, in relatively good agreement with CMB constraints from Planck (less than $\sim1.8σ$ tension), although we obtain a low value of $Ω_{\rm m}=0.212^{+0.017}_{-0.032}$, in tension with Planck at the $\sim3σ$ level. We show that this can be recast as an $H_0$ tension if one parametrises the amplitude of fluctuations and matter abundance in terms of variables without hidden dependence on $H_0$. Furthermore, we find that this tension reduces significantly after including a prior on the distance-redshift relationship from BAO data, without worsening the fit. In terms of baryonic effects, we show that failing to model and marginalise over them on scales $\ell\lesssim2000$ does not significantly affect the posterior constraints for DES-Y3 and KiDS-1000, but has a mild effect on deeper samples, such as HSC-DR1. This is in agreement with our ability to only mildly constrain the parameters of the Baryon Correction Model with these data
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
DES Y3 cosmic shear down to small scales: constraints on cosmology and baryons
Authors:
Giovanni Aricò,
Raul E. Angulo,
Matteo Zennaro,
Sergio Contreras,
Angela Chen,
Carlos Hernández-Monteagudo
Abstract:
We present the first analysis of cosmic shear measured in DES Y3 that employs the entire range of angular scales in the data. To achieve this, we build upon recent advances in the theoretical modelling of weak lensing provided by a combination of $N$-body simulations, physical models of baryonic processes, and neural networks. Specifically, we use BACCOemu to model the linear and nonlinear matter…
▽ More
We present the first analysis of cosmic shear measured in DES Y3 that employs the entire range of angular scales in the data. To achieve this, we build upon recent advances in the theoretical modelling of weak lensing provided by a combination of $N$-body simulations, physical models of baryonic processes, and neural networks. Specifically, we use BACCOemu to model the linear and nonlinear matter power spectrum including baryonic physics, allowing us to robustly exploit scales smaller than those used by the DES Collaboration. We show that the additional data produce cosmological parameters that are tighter but consistent with those obtained from larger scales, while also constraining the distribution of baryons. In particular, we measure the mass scale at which haloes have lost half of their gas, $\log\,M_{\rm c}=14.38^{+0.60}_{-0.56}\log(h^{-1}{\rm M_{ \odot}})$, and a parameter that quantifies the weighted amplitudes of the present-day matter inhomogeneities, $S_8=0.799^{+0.023}_{-0.015}$. Our constraint on $S_8$ is statistically compatible with that inferred from the Planck satellite's data at the $0.9σ$ level. We find instead a $1.4σ$ shift in comparison to that from the official DES Y3 cosmic shear, because of different choices in the modelling of intrinsic alignment, non-linearities, baryons, and lensing shear ratios. We conclude that small scales in cosmic shear data contain valuable astrophysical and cosmological information and thus should be included in standard analyses.
△ Less
Submitted 7 August, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
Consistent and simultaneous modelling of galaxy clustering and galaxy-galaxy lensing with Subhalo Abundance Matching
Authors:
Sergio Contreras,
Raul E. Angulo,
Jonás Chaves-Montero,
Simon D. M. White,
Giovanni Aricò
Abstract:
The spatial distribution of galaxies and their gravitational lensing signal offer complementary tests of galaxy formation physics and cosmology. However, their synergy can only be fully exploited if both probes are modelled accurately and consistently. In this paper, we demonstrate that this can be achieved using an extension of Subhalo Abundance Matching, dubbed SHAMe. Specifically, we use mock c…
▽ More
The spatial distribution of galaxies and their gravitational lensing signal offer complementary tests of galaxy formation physics and cosmology. However, their synergy can only be fully exploited if both probes are modelled accurately and consistently. In this paper, we demonstrate that this can be achieved using an extension of Subhalo Abundance Matching, dubbed SHAMe. Specifically, we use mock catalogues built from the TNG300 hydrodynamical simulation to show that SHAMe can simultaneously model the multipoles of the redshift-space galaxy correlation function and galaxy-galaxy lensing, without noticeable bias within the statistical sampling uncertainties of a SDSS volume and on scales r = [0.6-30] Mpc/h. Modelling the baryonic processes in galaxy-galaxy lensing with a baryonification scheme allows SHAMe's range of validity to be extended to r = [0.1-30] Mpc/h. Remarkably, our model achieves this level of precision with just five free parameters beyond those describing the baryonification model. At fixed cosmology, we find that galaxy-galaxy lensing provides a general consistency test but little additional information on galaxy modelling parameters beyond that encoded in the redshift-space multipoles. It does, however, improve constraints if only the projected correlation function is available, as in surveys with only photometric redshifts. We expect SHAMe to have a higher fidelity across a wider range of scales than more traditional methods such as Halo Occupation Distribution modelling. Thus it should provide a significantly more powerful and more robust tool for analysing next-generation large-scale surveys.
△ Less
Submitted 10 January, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.
-
The cosmology dependence of the concentration-mass-redshift relation
Authors:
Daniel López-Cano,
Raúl E. Angulo,
Aaron D. Ludlow,
M. Zennaro,
S. Contreras,
Jonás Chaves-Montero,
G. Aricò
Abstract:
The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift -- the so-called concentration-mass-redshift relation, denoted $c(M,z)$. We present here an extensive study of the cosmology-dependence of $c(M,z)$ that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological paramet…
▽ More
The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift -- the so-called concentration-mass-redshift relation, denoted $c(M,z)$. We present here an extensive study of the cosmology-dependence of $c(M,z)$ that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological parameters were varied: $σ_{8}$, $Ω_{\mathrm{M}}$, $Ω_{\mathrm{b}}$, $n_{\mathrm{s}}$, $h$, $M_ν$, $w_{0}$ and $w_{\mathrm{a}}$. We characterize the impact of these parameters on concentrations for different halo masses and redshifts. In agreement with previous works, and for all cosmologies studied, we find that there exists a tight correlation between the characteristic densities of dark matter haloes within their scale radii, $r_{-2}$, and the critical density of the Universe at a suitably defined formation time. This finding, when combined with excursion set modelling of halo formation histories, allows us to accurately predict the concentrations of dark matter haloes as a function of mass, redshift, and cosmology. We use our simulations to test the reliability of a number of published models for predicting halo concentration and highlight when they succeed or fail to reproduce the cosmological $c(M,z)$ relation.
△ Less
Submitted 27 July, 2022;
originally announced July 2022.
-
The Bacco Simulation Project: Bacco Hybrid Lagrangian Bias Expansion Model in Redshift Space
Authors:
Marcos Pellejero-Ibanez,
Raul E. Angulo,
Matteo Zennaro,
Jens Stuecker,
Sergio Contreras,
Giovanni Arico,
Francisco Maion
Abstract:
We present an emulator that accurately predicts the power spectrum of galaxies in redshift space as a function of cosmological parameters. Our emulator is based on a 2nd-order Lagrangian bias expansion that is displaced to Eulerian space using cosmological $N$-body simulations. Redshift space distortions are then imprinted using the non-linear velocity field of simulated particles and haloes. We b…
▽ More
We present an emulator that accurately predicts the power spectrum of galaxies in redshift space as a function of cosmological parameters. Our emulator is based on a 2nd-order Lagrangian bias expansion that is displaced to Eulerian space using cosmological $N$-body simulations. Redshift space distortions are then imprinted using the non-linear velocity field of simulated particles and haloes. We build the emulator using a forward neural network trained with the simulations of the BACCO project, which covers an 8-dimensional parameter space including massive neutrinos and dynamical dark energy. We show that our emulator provides unbiased cosmological constraints from the monopole, quadrupole, and hexadecapole of a mock galaxy catalogue that mimics the BOSS-CMASS sample down to nonlinear scales ($k\sim0.6$[$h/$Mpc]$^{3}$). This work opens up the possibility of robustly extracting cosmological information from small scales using observations of the large-scale structure of the Universe.
△ Less
Submitted 27 July, 2022; v1 submitted 13 July, 2022;
originally announced July 2022.
-
Constraining the Baryonic Feedback with Cosmic Shear Using the DES Year-3 Small-Scale Measurements
Authors:
A. Chen,
G. Aricò,
D. Huterer,
R. Angulo,
N. Weaverdyck,
O. Friedrich,
L. F. Secco,
C. Hernández-Monteagudo,
A. Alarcon,
O. Alves,
A. Amon,
F. Andrade-Oliveira,
E. Baxter,
K. Bechtol,
M. R. Becker,
G. M. Bernstein,
J. Blazek,
A. Brandao-Souza,
S. L. Bridle,
H. Camacho,
A. Campos,
A. Carnero Rosell,
M. Carrasco Kind,
R. Cawthon,
C. Chang
, et al. (117 additional authors not shown)
Abstract:
We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package \texttt{Baccoemu} to accelerate the evaluation of the baryonic nonlinear matter power spectrum. We design our ana…
▽ More
We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package \texttt{Baccoemu} to accelerate the evaluation of the baryonic nonlinear matter power spectrum. We design our analysis pipeline to focus on the constraints of the baryonic suppression effects, utilizing the implication given by a principal component analysis on the Fisher forecasts. Our constraint on the baryonic effects can then be used to better model and ameliorate the effects of baryons in producing cosmological constraints from the next generation large-scale structure surveys. We detect the baryonic suppression on the cosmic shear measurements with a $\sim 2 σ$ significance. The characteristic halo mass for which half of the gas is ejected by baryonic feedback is constrained to be $M_c > 10^{13.2} h^{-1} M_{\odot}$ (95\% C.L.). The best-fit baryonic suppression is $\sim 5\%$ at $k=1.0 {\rm Mpc}\ h^{-1}$ and $\sim 15\%$ at $k=5.0 {\rm Mpc} \ h^{-1}$. Our findings are robust with respect to the assumptions about the cosmological parameters, specifics of the baryonic model, and intrinsic alignments.
△ Less
Submitted 17 June, 2022;
originally announced June 2022.
-
Modelling galaxy clustering in redshift space with a Lagrangian bias formalism and $N$-body simulations
Authors:
Marcos Pellejero-Ibanez,
Jens Stuecker,
Raul E. Angulo,
Matteo Zennaro,
Sergio Contreras,
Giovanni Arico
Abstract:
Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describ…
▽ More
Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describe redshift-space distortions even down to small scales. This model corresponds to a second-order perturbative Lagrangian bias expansion which is advected to Eulerian space employing a displacement field extracted from $N$-body simulations. Eulerian coordinates are then transformed into redshift space by directly employing simulated velocity fields augmented with nuisance parameters capturing various possible satellite fractions and intra-halo small-scale velocities. We quantify the accuracy of our approach against samples of physically-motivated mock galaxies selected according to either Stellar Mass (SM) or Star Formation Rate (SFR) at multiple abundances and at $z=0$ and $1$. We find our model describes the monopole, quadrupole, and hexadecapole of the galaxy-power spectra down to scales of $k\approx 0.6 [h/$Mpc] within the accuracy of our simulations. This approach could pave the way to significantly increase the amount of cosmological information to be extracted from future galaxy surveys.
△ Less
Submitted 17 September, 2021;
originally announced September 2021.
-
Accelerating Large-Scale-Structure data analyses by emulating Boltzmann solvers and Lagrangian Perturbation Theory
Authors:
Giovanni Aricò,
Raul E. Angulo,
Matteo Zennaro
Abstract:
The linear matter power spectrum is an essential ingredient in all theoretical models for interpreting large-scale-structure observables. Although Boltzmann codes such as CLASS or CAMB are very efficient at computing the linear spectrum, the analysis of data usually requires $10^4$-$10^6$ evaluations, which means this task can be the most computationally expensive aspect of data analysis. Here, we…
▽ More
The linear matter power spectrum is an essential ingredient in all theoretical models for interpreting large-scale-structure observables. Although Boltzmann codes such as CLASS or CAMB are very efficient at computing the linear spectrum, the analysis of data usually requires $10^4$-$10^6$ evaluations, which means this task can be the most computationally expensive aspect of data analysis. Here, we address this problem by building a neural network emulator that provides the linear theory (total and cold) matter power spectrum in about one millisecond with 0.2% (0.5%) accuracy over redshifts $z \le 3$ ($z \le 9$), and scales $10^{-4} \le k \, [h {\rm Mpc^{-1}}] < 50$. We train this emulator with more than 200,000 measurements, spanning a broad cosmological parameter space that includes massive neutrinos and dynamical dark energy. We show that the parameter range and accuracy of our emulator is enough to get unbiased cosmological constraints in the analysis of a Euclid-like weak lensing survey. Complementing this emulator, we train 15 other emulators for the cross-spectra of various linear fields in Eulerian space, as predicted by 2nd-order Lagrangian Perturbation theory, which can be used to accelerate perturbative bias descriptions of galaxy clustering. Our emulators are specially designed to be used in combination with emulators for the nonlinear matter power spectrum and for baryonic effects, all of which are publicly available at http://www.dipc.org/bacco.
△ Less
Submitted 17 June, 2022; v1 submitted 29 April, 2021;
originally announced April 2021.
-
Non-universality of the mass function: dependence on the growth rate and power spectrum shape
Authors:
Lurdes Ondaro-Mallea,
Raul E. Angulo,
Matteo Zennaro,
Sergio Contreras,
Giovanni Aricò
Abstract:
The abundance of dark matter haloes is one of the key probes of the growth of structure and expansion history of the Universe. Theoretical predictions for this quantity usually assume that, when expressed in a certain form, it depends only on the mass variance of the linear density field. However, cosmological simulations have revealed that this assumption breaks, leading to 10-20% systematic effe…
▽ More
The abundance of dark matter haloes is one of the key probes of the growth of structure and expansion history of the Universe. Theoretical predictions for this quantity usually assume that, when expressed in a certain form, it depends only on the mass variance of the linear density field. However, cosmological simulations have revealed that this assumption breaks, leading to 10-20% systematic effects. In this paper we employ a specially-designed suite of simulations to further investigate this problem. Specifically, we carry out cosmological N-body simulations where we systematically vary growth history at a fixed linear density field, or vary the power spectrum shape at a fixed growth history. We show that the halo mass function generically depends on these quantities, thus showing a clear signal of non-universality. Most of this effect can be traced back to the way in which the same linear fluctuation grows differently into the nonlinear regime depending on details of its assembly history. With these results, we propose a parameterization with explicit dependence on the linear growth rate and power spectrum shape. Using an independent suite of simulations, we show that this fitting function accurately captures the mass function of haloes over cosmologies spanning a vast parameter space, including massive neutrinos and dynamical dark energy. Finally, we employ this tool to improve the accuracy of so-called cosmology-rescaling methods and show they can deliver 2% accurate predictions for the halo mass function over the whole range of currently viable cosmologies.
△ Less
Submitted 17 February, 2021;
originally announced February 2021.
-
The BACCO simulation project: biased tracers in real space
Authors:
Matteo Zennaro,
Raul E. Angulo,
Marcos Pellejero-Ibáñez,
Jens Stücker,
Sergio Contreras,
Giovanni Aricò
Abstract:
We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advecte…
▽ More
We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advected to Eulerian space using the displacement field of numerical simulations. The cosmology-dependence is captured thanks to a cosmology-rescaling algorithm. We show that our emulator is capable of describing the power spectrum of galaxy formation simulations for a sample mimicking that of a typical Emission-Line survey at $z \sim 1$ with an accuracy of $1-2\%$ up to nonlinear scales $k \sim 0.7 h \mathrm{Mpc}^{-1}$.
△ Less
Submitted 2 November, 2023; v1 submitted 28 January, 2021;
originally announced January 2021.
-
The BACCO Simulation Project: A baryonification emulator with Neural Networks
Authors:
Giovanni Aricò,
Raul E. Angulo,
Sergio Contreras,
Lurdes Ondaro-Mallea,
Marcos Pellejero-Ibañez,
Matteo Zennaro
Abstract:
We present a neural-network emulator for baryonic effects in the non-linear matter power spectrum. We calibrate this emulator using more than 50,000 measurements in a 15-dimensional parameters space, varying cosmology and baryonic physics. Baryonic physics is described through a baryonification algorithm, that has been shown to accurately capture the relevant effects on the power spectrum and bisp…
▽ More
We present a neural-network emulator for baryonic effects in the non-linear matter power spectrum. We calibrate this emulator using more than 50,000 measurements in a 15-dimensional parameters space, varying cosmology and baryonic physics. Baryonic physics is described through a baryonification algorithm, that has been shown to accurately capture the relevant effects on the power spectrum and bispectrum in state-of-the-art hydrodynamical simulations. Cosmological parameters are sampled using a cosmology-rescaling approach including massive neutrinos and dynamical dark energy. The specific quantity we emulate is the ratio between matter power spectrum with baryons and gravity-only, and we estimate the overall precision of the emulator to be 1-2%, at all scales 0.01 < k < 5 h/Mpc, and redshifts 0 < z < 1.5. We also obtain an accuracy of 1-2%, when testing the emulator against a collection of 74 different cosmological hydrodynamical simulations and their respective gravity-only counterparts. We show also that only one baryonic parameter, namely Mc, which set the gas fraction retained per halo mass, is enough to have accurate and realistic predictions of the baryonic feedback at a given epoch. Our emulator will become publicly available in http://www.dipc.org/bacco.
△ Less
Submitted 30 November, 2020;
originally announced November 2020.
-
Simultaneous modelling of matter power spectrum and bispectrum in the presence of baryons
Authors:
Giovanni Aricò,
Raul E. Angulo,
Carlos Hernández-Monteagudo,
Sergio Contreras,
Matteo Zennaro
Abstract:
We demonstrate that baryonification algorithms, which displace particles in gravity-only simulations according to physically-motivated prescriptions, can simultaneously capture the impact of baryonic physics on the 2 and 3-point statistics of matter. Specifically, we show that our implementation of a baryonification algorithm jointly fits the changes induced by baryons on the power spectrum and eq…
▽ More
We demonstrate that baryonification algorithms, which displace particles in gravity-only simulations according to physically-motivated prescriptions, can simultaneously capture the impact of baryonic physics on the 2 and 3-point statistics of matter. Specifically, we show that our implementation of a baryonification algorithm jointly fits the changes induced by baryons on the power spectrum and equilateral bispectrum on scales up to k < 5 h/Mpc and redshifts z<2, as measured in six different cosmological hydrodynamical simulations. The accuracy of our fits are typically 1% for the power spectrum, and for the equilateral and squeezed bispectra, which somewhat degrades to 3% for simulations with extreme feedback prescriptions. Our results support the physical assumptions underlying baryonification approaches, and encourage their use in interpreting weak gravitational lensing and other cosmological observables.
△ Less
Submitted 29 September, 2020;
originally announced September 2020.
-
Tomographic Constraints on Gravity from Angular Redshift Fluctuations in the Late Universe
Authors:
Carlos Hernández-Monteagudo,
Jonás Chaves-Montero,
Raúl E. Angulo,
Giovanni Aricò
Abstract:
Fluctuations in sky maps of the galaxy redshifts, dubbed as angular redshift fluctuations (ARF), contain precise information about the growth rate of structures and the nature of gravity in the Universe. Specifically, ARF constrain the combination of cosmological parameters $H/H_0\,fσ_8(z)$, while being an intrinsically tomographic probe and largely insensitive to many observational systematic err…
▽ More
Fluctuations in sky maps of the galaxy redshifts, dubbed as angular redshift fluctuations (ARF), contain precise information about the growth rate of structures and the nature of gravity in the Universe. Specifically, ARF constrain the combination of cosmological parameters $H/H_0\,fσ_8(z)$, while being an intrinsically tomographic probe and largely insensitive to many observational systematic errors, all this without requiring the assumption of any redshift-to-distance relation under a given fiducial cosmology. We present the first cosmological constraints derived from ARF by using BOSS LOWZ+CMASS DR12 galaxy samples, obtaining 7\%-accurate constraints on $H/H_0 fσ_8(z)$ at more than 20 redshifts over the range $z \in [0.26,0.72]$. Our best-fitting value is $10\%$ larger, but compatible at the $1.4σ$ level, than the $Λ$CDM expectation set by {\it Planck} observations of the Cosmic Microwave Background (CMB) radiation. Our tomographic measurements, combined with these CMB data, provides one of the strongest constraints on the gravity index $γ$, $γ=0.44^{+0.09}_{-0.07}$, which lies within $2σ$ from the prediction of General Relativity ($γ_{\rm GR}\simeq 0.55$).
△ Less
Submitted 23 March, 2021; v1 submitted 13 May, 2020;
originally announced May 2020.
-
The BACCO Simulation Project: Exploiting the full power of large-scale structure for cosmology
Authors:
Raul E. Angulo,
Matteo Zennaro,
Sergio Contreras,
Giovanni Aricò,
Marcos Pellejero-Ibañez,
Jens Stücker
Abstract:
We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $\sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear m…
▽ More
We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $\sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear mass power spectrum over the redshift range $0<z<1.5$ and over scales $10^{-2} < k/(h Mpc^{-1} ) < 5$ for 800 points in an 8-dimensional cosmological parameter space. For an efficient interpolation of the results, we build an emulator and compare its predictions against several widely-used methods. Over the whole range of scales considered, we expect our predictions to be accurate at the 2\% level for parameters in the minimal $Λ$ CDM model and to 3\% when extended to dynamical dark energy and massive neutrinos. We make our emulator publicly available under http://www.dipc.org/bacco
△ Less
Submitted 6 July, 2021; v1 submitted 13 April, 2020;
originally announced April 2020.
-
3\%-accurate predictions for the clustering of dark matter, haloes and subhaloes, over a wide range of cosmologies and scales
Authors:
S. Contreras,
R. E. Angulo,
M. Zennaro,
G. Aricó,
M. Pellejero-Ibañez
Abstract:
Predicting the spatial distribution of objects as a function of cosmology is an essential ingredient for the exploitation of future galaxy surveys. In this paper we show that a specially-designed suite of gravity-only simulations together with cosmology-rescaling algorithms can provide the clustering of dark matter, haloes, and subhaloes with high precision. Specifically, with only 3 $N$-body simu…
▽ More
Predicting the spatial distribution of objects as a function of cosmology is an essential ingredient for the exploitation of future galaxy surveys. In this paper we show that a specially-designed suite of gravity-only simulations together with cosmology-rescaling algorithms can provide the clustering of dark matter, haloes, and subhaloes with high precision. Specifically, with only 3 $N$-body simulations we obtain the power spectrum of dark matter at $z=0$ and $z=1$ to better than 3\% precision for essentially all currently viable values of 8 cosmological parameters, including massive neutrinos and dynamical dark energy, and over the whole range of scales explored, 0.03 < $k/h^{-1}Mpc$ < 5. This precision holds at the same level for mass-selected haloes and for subhaloes selected according to their peak maximum circular velocity. As an initial application of these predictions, we successfully constrain $Ω_{\rm m}$, $σ_8$, and the scatter in subhalo-abundance-matching employing the projected correlation function of mock SDSS galaxies.
△ Less
Submitted 22 September, 2020; v1 submitted 9 January, 2020;
originally announced January 2020.
-
Cosmological parameter estimation via iterative emulation of likelihoods
Authors:
Marcos Pellejero-Ibañez,
Raul E. Angulo,
Giovanni Aricó,
Matteo Zennaro,
Sergio Contreras,
Jens Stücker
Abstract:
The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterativ…
▽ More
The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterative Gaussian emulation of the target likelihood function. This requires a minimal number of likelihood evaluations and naturally accommodates for stochasticity in theoretical models. We apply the algorithm to estimate 9 parameters from the monopole and quadrupole of a mock power spectrum in redshift space. We obtain accurate posterior distribution functions with approximately 100 times fewer likelihood evaluations than an affine invariant MCMC, roughly independently from the dimensionality of the problem. We anticipate that our parameter estimation algorithm will accelerate the adoption of more accurate theoretical models in data analysis, enabling more comprehensive exploitation of cosmological observables.
△ Less
Submitted 23 January, 2020; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Modelling the large scale structure of the Universe as a function of cosmology and baryonic physics
Authors:
Giovanni Aricò,
Raul E. Angulo,
Carlos Hernández-Monteagudo,
Sergio Contreras,
Matteo Zennaro,
Marcos Pellejero-Ibañez,
Yetli Rosas-Guevara
Abstract:
We present and test a framework that models the three-dimensional distribution of mass in the Universe as a function of cosmological and astrophysical parameters. Our approach combines two different techniques: a rescaling algorithm that modifies the cosmology of gravity-only N-body simulations, and a baryonification algorithm which mimics the effects of astrophysical processes induced by baryons,…
▽ More
We present and test a framework that models the three-dimensional distribution of mass in the Universe as a function of cosmological and astrophysical parameters. Our approach combines two different techniques: a rescaling algorithm that modifies the cosmology of gravity-only N-body simulations, and a baryonification algorithm which mimics the effects of astrophysical processes induced by baryons, such as star formation and AGN feedback. We show how this approach can accurately reproduce the effects of baryons on the matter power spectrum of various state-of-the-art hydro-dynamical simulations (EAGLE, Illustris, Illustris-TNG, Horizon-AGN, and OWLS,Cosmo-OWLS and BAHAMAS), to percent level from very large down to small, highly nonlinear scales, k= 5 h/Mpc, and from z=0 up to z=2. We highlight that, thanks to the heavy optimisation of the algorithms, we can obtain these predictions for arbitrary baryonic models and cosmology (including massive neutrinos and dynamical dark energy models) with an almost negligible CPU cost. Therefore, this approach is efficient enough for cosmological data analyses. With these tools in hand we explore the degeneracies between cosmological and astrophysical parameters in the nonlinear mass power spectrum. Our findings suggest that after marginalising over baryonic physics, cosmological constraints inferred from weak gravitational lensing should be moderately degraded.
△ Less
Submitted 17 January, 2020; v1 submitted 19 November, 2019;
originally announced November 2019.
-
How to add massive neutrinos to your $Λ$CDM simulation -- extending cosmology rescaling algorithms
Authors:
Matteo Zennaro,
Raúl E. Angulo,
Giovanni Aricò,
Sergio Contreras,
Marcos Pellejero-Ibáñez
Abstract:
Providing accurate predictions for the spatial distribution of matter and luminous tracers in the presence of massive neutrinos is an important task, given the imminent arrival of highly accurate large-scale structure observations. In this work, we address this challenge by extending cosmology-rescaling algorithms to massive neutrino cosmologies. In this way, a $Λ$CDM simulation can be modified to…
▽ More
Providing accurate predictions for the spatial distribution of matter and luminous tracers in the presence of massive neutrinos is an important task, given the imminent arrival of highly accurate large-scale structure observations. In this work, we address this challenge by extending cosmology-rescaling algorithms to massive neutrino cosmologies. In this way, a $Λ$CDM simulation can be modified to provide nonlinear structure formation predictions in the presence a hot component of arbitrary mass, and, if desired, to include non-gravitational modifications to the clustering of matter on large scales. We test the accuracy of the method by comparing its predictions to a suite of simulations carried out explicitly including a neutrino component in its evolution equations. We find that, for neutrino masses in the range $M_ν\in [0.06, 0.3] ~ \mathrm{eV}$ the matter power spectrum is recovered to better than $1\%$ on all scales $k<2~h~\mathrm{Mpc}^{-1}$. Similarly, the halo mass function is predicted at a few percent level over the range $M_{\rm halo} \in [10^{12}, 10^{15}] ~ h^{-1} ~ \mathrm{M}_{\odot}$, and so do also the multipoles of the galaxy 2-point correlation function in redshift space over $r \in [0.1, 200] ~ h^{-1} ~ \mathrm{Mpc}$. We provide parametric forms for the necessary transformations, as a function of $Ω_{\rm m}$ and $Ω_ν$ for various target redshifts.
△ Less
Submitted 3 December, 2019; v1 submitted 21 May, 2019;
originally announced May 2019.