-
K-mouflage at high k: extending the reach of $\texttt{Hi-COLA}$
Authors:
Ashim Sen Gupta,
Bartolomeo Fiorini,
Tessa Baker
Abstract:
The $\texttt{Hi-COLA}$ code is an efficient dark matter simulation suite that flexibly handles the Horndeski family of modified gravity models. In this work we extend the scope of $\texttt{Hi-COLA}$ to accommodate Horndeski theories with K-mouflage screening, allowing for the computation of matter power spectra in the non-linear regime in these models. We explore the boost of the dark matter power…
▽ More
The $\texttt{Hi-COLA}$ code is an efficient dark matter simulation suite that flexibly handles the Horndeski family of modified gravity models. In this work we extend the scope of $\texttt{Hi-COLA}$ to accommodate Horndeski theories with K-mouflage screening, allowing for the computation of matter power spectra in the non-linear regime in these models. We explore the boost of the dark matter power spectrum relative to GR-$Λ$CDM in K-mouflage gravity, and also discuss how large-scale structure computations change between the Einstein and Jordan frames. A dissection of the relative contributions of the modified background, linear growth, fifth force, and the conformal factor (a new inclusion to $\texttt{Hi-COLA}$) to the boost factor is presented. The ability of $\texttt{Hi-COLA}$ to run with general Horndeski models and multiple screening mechanisms makes it an ideal tool for testing gravity with upcoming galaxy survey data.
△ Less
Submitted 30 June, 2024;
originally announced July 2024.
-
Matter Power Spectra in Modified Gravity: A Comparative Study of Approximations and $N$-Body Simulations
Authors:
Benjamin Bose,
Ashim Sen Gupta,
Bartolomeo Fiorini,
Guilherme Brando,
Farbod Hassani,
Tessa Baker,
Lucas Lombriser,
Baojiu Li,
Cheng-Zong Ruan,
Cesar Hernandez-Aguayo,
Luis Atayde,
Noemi Frusciante
Abstract:
Testing gravity and the concordance model of cosmology, $Λ$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and…
▽ More
Testing gravity and the concordance model of cosmology, $Λ$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and K-mouflage. In particular, we compare the predictions from full $N$-body simulations, two $N$-body codes with approximate time integration schemes, a parametrised modified $N$-body implementation and the analytic halo model reaction approach. We find the modification to the $Λ$CDM spectrum is in $2\%$ agreement for $z\leq1$ and $k\leq 1~h/{\rm Mpc}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.
△ Less
Submitted 19 June, 2024;
originally announced June 2024.
-
The Stability of the BAO Linear Point under Modified Gravity
Authors:
Jaemyoung Jason Lee,
Bartolomeo Fiorini,
Farnik Nikakhtar,
Ravi K. Sheth
Abstract:
Baryon Acoustic Oscillations (BAOs) are crucial in cosmological analysis, providing a standard ruler, as well as constraints on dark energy. In General Relativity models, the BAO Linear Point - the midpoint between the dip and the peak in the correlation function - has been shown to be rather robust to evolution and redshift space distortions. We show that this remains true even when the gravity m…
▽ More
Baryon Acoustic Oscillations (BAOs) are crucial in cosmological analysis, providing a standard ruler, as well as constraints on dark energy. In General Relativity models, the BAO Linear Point - the midpoint between the dip and the peak in the correlation function - has been shown to be rather robust to evolution and redshift space distortions. We show that this remains true even when the gravity model is not General Relativity, at least for $f(R)$ and DGP gravity models which have the same expansion history as the standard $Λ$CDM. For the Linear Point to be able to distinguish between modified gravity (MG) and $Λ$CDM, survey volumes of order tens of cubic Gpc are required.
△ Less
Submitted 19 June, 2024; v1 submitted 13 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.
-
Euclid preparation. Modelling spectroscopic clustering on mildly nonlinear scales in beyond-$Λ$CDM models
Authors:
Euclid Collaboration,
B. Bose,
P. Carrilho,
M. Marinucci,
C. Moretti,
M. Pietroni,
E. Carella,
L. Piga,
B. S. Wright,
F. Vernizzi,
C. Carbone,
S. Casas,
G. D'Amico,
N. Frusciante,
K. Koyama,
F. Pace,
A. Pourtsidou,
M. Baldi,
L. F. de la Bella,
B. Fiorini,
C. Giocoli,
L. Lombriser,
N. Aghanim,
A. Amara,
S. Andreon
, et al. (207 additional authors not shown)
Abstract:
We investigate the approximations needed to efficiently predict the large-scale clustering of matter and dark matter halos in beyond-$Λ$CDM scenarios. We examine the normal branch of the Dvali-Gabadadze-Porrati model, the Hu-Sawicki $f(R)$ model, a slowly evolving dark energy, an interacting dark energy model and massive neutrinos. For each, we test approximations for the perturbative kernel calcu…
▽ More
We investigate the approximations needed to efficiently predict the large-scale clustering of matter and dark matter halos in beyond-$Λ$CDM scenarios. We examine the normal branch of the Dvali-Gabadadze-Porrati model, the Hu-Sawicki $f(R)$ model, a slowly evolving dark energy, an interacting dark energy model and massive neutrinos. For each, we test approximations for the perturbative kernel calculations, including the omission of screening terms and the use of perturbative kernels based on the Einstein-de Sitter universe; we explore different infrared-resummation schemes, tracer bias models and a linear treatment of massive neutrinos; we employ two models for redshift space distortions, the Taruya-Nishimishi-Saito prescription and the Effective Field Theory of Large-Scale Structure. This work further provides a preliminary validation of the codes being considered by Euclid for the spectroscopic clustering probe in beyond-$Λ$CDM scenarios. We calculate and compare the $χ^2$ statistic to assess the different modelling choices. This is done by fitting the spectroscopic clustering predictions to measurements from numerical simulations and perturbation theory-based mock data. We compare the behaviour of this statistic in the beyond-$Λ$CDM cases, as a function of the maximum scale included in the fit, to the baseline $Λ$CDM case. We find that the Einstein-de Sitter approximation without screening is surprisingly accurate for all cases when comparing to the halo clustering monopole and quadrupole obtained from simulations. Our results suggest that the inclusion of multiple redshift bins, higher-order multipoles, higher-order clustering statistics (such as the bispectrum) and photometric probes such as weak lensing, will be essential to extract information on massive neutrinos, modified gravity and dark energy.
△ Less
Submitted 11 July, 2024; v1 submitted 22 November, 2023;
originally announced November 2023.
-
Fast production of cosmological emulators in modified gravity: the matter power spectrum
Authors:
Bartolomeo Fiorini,
Kazuya Koyama,
Tessa Baker
Abstract:
We test the convergence of fast simulations based on the COmoving Lagrangian Acceleration (COLA) method for predictions of the matter power spectrum, specialising our analysis in the redshift range $1 \le z \le 1.65$, relevant to high-redshift spectroscopic galaxy surveys. We then focus on the enhancement of the matter power spectrum in modified gravity (MG), the boost factor, using the Dvali-Gaba…
▽ More
We test the convergence of fast simulations based on the COmoving Lagrangian Acceleration (COLA) method for predictions of the matter power spectrum, specialising our analysis in the redshift range $1 \le z \le 1.65$, relevant to high-redshift spectroscopic galaxy surveys. We then focus on the enhancement of the matter power spectrum in modified gravity (MG), the boost factor, using the Dvali-Gabadadze-Porrati (DGP) theory as a test case but developing a general approach that can be applied to other MG theories. After identifying the minimal simulation requirements for accurate DGP boost factors, we design and produce a COLA simulation suite that we use to train a neural network emulator for the DGP boost factor. Using MG-AREPO simulations as a reference, we estimate the emulator accuracy to be of $\sim 3\%$ up to $k=5 \, h {\rm Mpc}^{-1}$ at $0 \leq z \leq 2$. We make the emulator publicly available at: https://github.com/BartolomeoF/nDGPemu.
△ Less
Submitted 20 November, 2023; v1 submitted 9 October, 2023;
originally announced October 2023.
-
Testing Gravity on Cosmological Scales: Theoretical Predictions with the COLA Method
Authors:
Bartolomeo Fiorini
Abstract:
One of the main objectives of stage IV galaxy surveys is to constrain gravity on cosmological scales. To this end, it is crucial to make accurate theoretical predictions in the nonlinear regime of structure formation in order to maximise the scientific return. This is possible at a relatively low computational cost thanks to COLA simulations, an approximate and much faster alternative to full $N$-…
▽ More
One of the main objectives of stage IV galaxy surveys is to constrain gravity on cosmological scales. To this end, it is crucial to make accurate theoretical predictions in the nonlinear regime of structure formation in order to maximise the scientific return. This is possible at a relatively low computational cost thanks to COLA simulations, an approximate and much faster alternative to full $N$-body simulations. In this thesis, we focus on two modified gravity theories, $f(R)$ and nDGP, and present an analysis of how COLA simulations can be combined with empirical models for the galaxy-halo connection to produce realistic mock galaxy catalogues in modified gravity. We then use the resulting galaxy catalogues to study the effects of modified gravity and validate COLA simulations for several probes of the large-scale structure like the galaxy power spectrum, bispectrum and the void-galaxy cross-correlation function. Finally, we propose the COLA method for the extension of cosmological emulators to modified gravity theories and showcase its potential by training an emulator for the nDGP boost factor of the matter power spectrum, i.e., the ratio of the power spectrum in nDGP gravity with that in general relativity.
△ Less
Submitted 10 March, 2023;
originally announced March 2023.
-
$\texttt{Hi-COLA}$: Fast, approximate simulations of structure formation in Horndeski gravity
Authors:
Bill S. Wright,
Ashim Sen Gupta,
Tessa Baker,
Georgios Valogiannis,
Bartolomeo Fiorini
Abstract:
We introduce $\texttt{Hi-COLA}$, a code designed to run fast, approximate $\textit{N}$-body simulations of non-linear structure formation in reduced Horndeski gravity. Given an input Lagrangian, $\texttt{Hi-COLA}$ dynamically constructs the appropriate field equations and consistently solves for the cosmological background, linear growth, and screened fifth force of that theory. Hence…
▽ More
We introduce $\texttt{Hi-COLA}$, a code designed to run fast, approximate $\textit{N}$-body simulations of non-linear structure formation in reduced Horndeski gravity. Given an input Lagrangian, $\texttt{Hi-COLA}$ dynamically constructs the appropriate field equations and consistently solves for the cosmological background, linear growth, and screened fifth force of that theory. Hence $\texttt{Hi-COLA}$ is a general, adaptable, and useful tool that allows the mildly non-linear regime of many Horndeski theories to be investigated for the first time, at low computational cost. In this work, we first describe the screening approximations and simulation setup of $\texttt{Hi-COLA}$ for theories with Vainshtein screening. We validate the code against traditional $\textit{N}$-body simulations for cubic Galileon gravity, finding $2.5\%$ agreement up to $k_{\rm max}=1.2~h/{\rm Mpc}$. To demonstrate the flexibility of $\texttt{Hi-COLA}$, we additionally run the first simulations of an extended shift-symmetric gravity theory. We use the consistency and modularity of $\texttt{Hi-COLA}$ to dissect how the modified background, linear growth, and screened fifth force all contribute to departures from $Λ$CDM in the non-linear matter power spectrum. $\texttt{Hi-COLA}$ can be found at https://github.com/Hi-COLACode/Hi-COLA .
△ Less
Submitted 22 March, 2023; v1 submitted 4 September, 2022;
originally announced September 2022.
-
Studying large-scale structure probes of modified gravity with COLA
Authors:
Bartolomeo Fiorini,
Kazuya Koyama,
Albert Izard
Abstract:
We study the effect of two Modified Gravity (MG) theories, $f(R)$ and nDGP, on three probes of large-scale structure, the real space power spectrum estimator $Q_0$, bispectrum and voids, and validate fast approximate COLA simulations against full $N$-body simulations for the prediction of these probes. We find that using the first three even multipoles of the redshift space power spectrum to estim…
▽ More
We study the effect of two Modified Gravity (MG) theories, $f(R)$ and nDGP, on three probes of large-scale structure, the real space power spectrum estimator $Q_0$, bispectrum and voids, and validate fast approximate COLA simulations against full $N$-body simulations for the prediction of these probes. We find that using the first three even multipoles of the redshift space power spectrum to estimate $Q_0$ is enough to reproduce the MG boost factors of the real space power spectrum for both halo and galaxy catalogues. By analysing the bispectrum and reduced bispectrum of Dark Matter (DM), we show that the strong MG signal present in the DM bispectrum is mainly due to the enhanced power spectrum. We warn about adopting screening approximations in simulations as this neglects non-linear contributions that can source a significant component of the MG bispectrum signal at the DM level, but we argue that this is not a problem for the bispectrum of galaxies in redshift space where the signal is dominated by the non-linear galaxy bias. Finally, we perform void-finding on our galaxy mock catalogues by the ZOBOV watershed algorithm. To apply a linear model for Redshift-Space Distortion (RSD) in the void-galaxy cross-correlation function, we first examine the effects of MG on the void profiles entering into the RSD model. We find relevant MG signals in the integrated-density, velocity dispersion and radial velocity profiles in the nDGP theory. Fitting the RSD model for the linear growth rate, we recover the linear theory prediction in an nDGP model, which is larger than the $Λ$CDM prediction at the $3 σ$ level. In $f(R)$ theory we cannot naively compare the results of the fit with the linear theory prediction as this is scale-dependent, but we obtain results that are consistent with the $Λ$CDM prediction.
△ Less
Submitted 21 November, 2022; v1 submitted 2 August, 2022;
originally announced August 2022.
-
Enabling matter power spectrum emulation in beyond-$Λ$CDM cosmologies with COLA
Authors:
Guilherme Brando,
Bartolomeo Fiorini,
Kazuya Koyama,
Hans A. Winther
Abstract:
We compare and validate COLA (COmoving Lagrangian Acceleration) simulations against existing emulators in the literature, namely Bacco and Euclid Emulator 2. Our analysis focuses on the non-linear response function, i.e., the ratio between the non-linear dark matter power spectrum in a given cosmology with respect to a pre-defined reference cosmology, which is chosen to be the Euclid Emulator 2 re…
▽ More
We compare and validate COLA (COmoving Lagrangian Acceleration) simulations against existing emulators in the literature, namely Bacco and Euclid Emulator 2. Our analysis focuses on the non-linear response function, i.e., the ratio between the non-linear dark matter power spectrum in a given cosmology with respect to a pre-defined reference cosmology, which is chosen to be the Euclid Emulator 2 reference cosmology in this paper. We vary three cosmological parameters, the total matter density, the amplitude of the primordial scalar perturbations and the spectral index. By comparing the COLA non-linear response function with those computed from each emulator in the redshift range $0 \leq z \leq 3$, we find that the COLA method is in excellent agreement with the two emulators for scales up to $k \sim 1 \ h$/Mpc as long as the deviations of the matter power spectrum from the reference cosmology are not too large. We validate the implementation of massive neutrinos in our COLA simulations by varying the sum of neutrino masses to three different values, $0.0$ eV, $0.058$ eV and $0.15$ eV. We show that all three non-linear prescriptions used in this work agree at the $1\%$ level at $k \leq 1 \ h$/Mpc. We then introduce the Effective Field Theory of Dark Energy in our COLA simulations using the $N$-body gauge method. We consider two different modified gravity models in which the growth of structure is enhanced or suppressed at small scales, and show that the response function with respect to the change of modified gravity parameters depends weakly on cosmological parameters in these models.
△ Less
Submitted 1 September, 2022; v1 submitted 21 March, 2022;
originally announced March 2022.
-
Fast generation of mock galaxy catalogues in modified gravity models with COLA
Authors:
Bartolomeo Fiorini,
Kazuya Koyama,
Albert Izard,
Hans A. Winther,
Bill S. Wright,
Baojiu Li
Abstract:
We investigate the viability of producing galaxy mock catalogues with COmoving Lagrangian Acceleration (COLA) simulations in Modified Gravity (MG) models employing the Halo Occupation Distribution (HOD) formalism. In this work, we focus on two theories of MG: $f(R)$ gravity with the chameleon mechanism, and a braneworld model (nDGP) that incorporates the Vainshtein mechanism. We use a suite of ful…
▽ More
We investigate the viability of producing galaxy mock catalogues with COmoving Lagrangian Acceleration (COLA) simulations in Modified Gravity (MG) models employing the Halo Occupation Distribution (HOD) formalism. In this work, we focus on two theories of MG: $f(R)$ gravity with the chameleon mechanism, and a braneworld model (nDGP) that incorporates the Vainshtein mechanism. We use a suite of full $N$-body simulations in MG as a benchmark to test the accuracy of COLA simulations. At the level of Dark Matter (DM), we show that COLA accurately reproduces the matter power spectrum up to $k \sim 1 h {\rm Mpc}^{-1}$, while it is less accurate in reproducing the velocity field. To produce halo catalogues, we find that the ROCKSTAR halo-finder does not perform well with COLA simulations. On the other hand, using a simple Friends-of-Friends (FoF) finder and an empirical mass conversion from FoF to spherical over-density masses, we are able to produce halo catalogues in COLA that are in good agreement with those in $N$-body simulations. To consider the effects of the MG fifth force on the halo profile, we derive simple fitting formulae for the concentration-mass and the velocity dispersion-mass relations that we calibrate using ROCKSTAR halo catalogues in $N$-body simulations. We then use these results to extend the HOD formalism to modified gravity simulations in COLA. We use an HOD model with five parameters that we tune to obtain galaxy catalogues in redshift space. We find that despite the great freedom of the HOD model, MG leaves characteristic imprints in the redshift space power spectrum multipoles and these features are well captured by the COLA galaxy catalogues.
△ Less
Submitted 23 August, 2021; v1 submitted 9 June, 2021;
originally announced June 2021.
-
Cosmic variance of $H_0$ in light of forthcoming high-redshift surveys
Authors:
Giuseppe Fanizza,
Bartolomeo Fiorini,
Giovanni Marozzi
Abstract:
Forthcoming surveys will extend the understanding of cosmological large scale structures up to unprecedented redshift. According to this perspective, we present a fully relativistic framework to evaluate the impact of stochastic inhomogeneities on the determination of the Hubble constant. To this aim, we work within linear perturbation theory and relate the fluctuations of the luminosity distance-…
▽ More
Forthcoming surveys will extend the understanding of cosmological large scale structures up to unprecedented redshift. According to this perspective, we present a fully relativistic framework to evaluate the impact of stochastic inhomogeneities on the determination of the Hubble constant. To this aim, we work within linear perturbation theory and relate the fluctuations of the luminosity distance-redshift relation, in the Cosmic Concordance model, to the intrinsic uncertainty associated to the measurement of $H_0$ from high-redshift surveys ($0.15\le z\le3.85$). We first present the detailed derivation of the luminosity distance-redshift relation 2-point correlation function and then provide analytical results for all the involved relativistic effects, such as peculiar velocity, lensing, time delay and (integrated) Sachs-Wolfe, and their angular spectra. Hence, we apply our analytical results to the study of high-redshift Hubble diagram, according to what has been recently claimed in literature. Following the specific of Euclid Deep Survey and LSST, we conclude that the cosmic variance associated with the measurement of the Hubble constant is at most of 0.1 %. Our work extends the analysis already done in literature for closer sources, where only peculiar velocity has been taken into account. We then conclude that deep surveys will provide an estimation of the $H_0$ which will be more precise than the one obtained from local sources, at least in regard of the intrinsic uncertainty related to a stochastic distribution of inhomogeneities.
△ Less
Submitted 13 September, 2021; v1 submitted 24 February, 2021;
originally announced February 2021.