-
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 19 January, 2024; v1 submitted 22 November, 2023;
originally announced November 2023.
-
Constraints on modified gravity from the BOSS galaxy survey
Authors:
Lorenzo Piga,
Marco Marinucci,
Guido D'Amico,
Massimo Pietroni,
Filippo Vernizzi,
Bill S. Wright
Abstract:
We develop a pipeline to set new constraints on scale-independent modified gravity, from the galaxy power spectrum in redshift space of BOSS DR12. The latter is modelled using the effective field theory of large-scale structure up to 1-loop order in perturbation theory. We test our pipeline on synthetic and simulated data, to assess systematic biases on the inferred cosmological parameters due to…
▽ More
We develop a pipeline to set new constraints on scale-independent modified gravity, from the galaxy power spectrum in redshift space of BOSS DR12. The latter is modelled using the effective field theory of large-scale structure up to 1-loop order in perturbation theory. We test our pipeline on synthetic and simulated data, to assess systematic biases on the inferred cosmological parameters due to marginalization and theoretical errors, and we apply it to the normal branch of the DGP model with a $Λ$CDM background. When applied to synthetic data and cosmological simulations, we observe biased posteriors due to the strong degeneracy between the nDGP parameter $Ω_{\rm rc}$ and the primordial amplitude of fluctuations $A_s$. Fixing the latter to the Planck central value, we obtain a posterior distribution with $Ω_{\rm rc}\lesssim 0.65$ at 95$\%$ C.L., under the assumption of a flat prior on $\log_{10} Ω_{\rm rc}$. This upper bound, however, depends strongly on the prior on $Ω_{\rm rc}$. To alleviate this effect, we provide an upper bound based on the Bayes factor between the nDGP model and $Λ$CDM model, which gives $Ω_{\rm rc}\lesssim 0.2$ at 95$\%$ C.L..
△ Less
Submitted 19 March, 2023; v1 submitted 22 November, 2022;
originally announced November 2022.
-
Euclid: Modelling massive neutrinos in cosmology -- a code comparison
Authors:
J. Adamek,
R. E. Angulo,
C. Arnold,
M. Baldi,
M. Biagetti,
B. Bose,
C. Carbone,
T. Castro,
J. Dakin,
K. Dolag,
W. Elbers,
C. Fidler,
C. Giocoli,
S. Hannestad,
F. Hassani,
C. Hernández-Aguayo,
K. Koyama,
B. Li,
R. Mauland,
P. Monaco,
C. Moretti,
D. F. Mota,
C. Partmann,
G. Parimbelli,
D. Potter
, et al. (111 additional authors not shown)
Abstract:
The measurement of the absolute neutrino mass scale from cosmological large-scale clustering data is one of the key science goals of the Euclid mission. Such a measurement relies on precise modelling of the impact of neutrinos on structure formation, which can be studied with $N$-body simulations. Here we present the results from a major code comparison effort to establish the maturity and reliabi…
▽ More
The measurement of the absolute neutrino mass scale from cosmological large-scale clustering data is one of the key science goals of the Euclid mission. Such a measurement relies on precise modelling of the impact of neutrinos on structure formation, which can be studied with $N$-body simulations. Here we present the results from a major code comparison effort to establish the maturity and reliability of numerical methods for treating massive neutrinos. The comparison includes eleven full $N$-body implementations (not all of them independent), two $N$-body schemes with approximate time integration, and four additional codes that directly predict or emulate the matter power spectrum. Using a common set of initial data we quantify the relative agreement on the nonlinear power spectrum of cold dark matter and baryons and, for the $N$-body codes, also the relative agreement on the bispectrum, halo mass function, and halo bias. We find that the different numerical implementations produce fully consistent results. We can therefore be confident that we can model the impact of massive neutrinos at the sub-percent level in the most common summary statistics. We also provide a code validation pipeline for future reference.
△ Less
Submitted 8 August, 2023; v1 submitted 22 November, 2022;
originally announced November 2022.
-
$\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.
-
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.
-
On the road to percent accuracy V: the non-linear power spectrum beyond $Λ$CDM with massive neutrinos and baryonic feedback
Authors:
Benjamin Bose,
Bill S. Wright,
Matteo Cataneo,
Alkistis Pourtsidou,
Carlo Giocoli,
Lucas Lombriser,
Ian G. McCarthy,
Marco Baldi,
Simon Pfeifer,
Qianli Xia
Abstract:
In the context of forthcoming galaxy surveys, to ensure unbiased constraints on cosmology and gravity when using non-linear structure information, percent-level accuracy is required when modelling the power spectrum. This calls for frameworks that can accurately capture the relevant physical effects, while allowing for deviations from $Λ$CDM. Massive neutrino and baryonic physics are two of the mo…
▽ More
In the context of forthcoming galaxy surveys, to ensure unbiased constraints on cosmology and gravity when using non-linear structure information, percent-level accuracy is required when modelling the power spectrum. This calls for frameworks that can accurately capture the relevant physical effects, while allowing for deviations from $Λ$CDM. Massive neutrino and baryonic physics are two of the most relevant such effects. We present an integration of the halo model reaction frameworks for massive neutrinos and beyond-$Λ$CDM cosmologies. The integrated halo model reaction, combined with a pseudo power spectrum modelled by HMCode2020 is then compared against $N$-body simulations that include both massive neutrinos and an $f(R)$ modification to gravity. We find that the framework is 4% accurate down to at least $k\approx 3 \, h/{\rm Mpc}$ for a modification to gravity of $|f_{\rm R0}|\leq 10^{-5}$ and for the total neutrino mass $M_ν\equiv \sum m_ν\leq 0.15$ eV. We also find that the framework is 4% consistent with EuclidEmulator2 as well as the Bacco emulator for \B{most of the considered} $νw$CDM cosmologies down to at least $k \approx 3 \, h$/Mpc. Finally, we compare against hydrodynamical simulations employing HMCode2020's baryonic feedback modelling on top of the halo model reaction. For $νΛ$CDM cosmologies we find 2% accuracy for $M_ν\leq 0.48$eV down to at least $k\approx 5h$/Mpc. Similar accuracy is found when comparing to $νw$CDM hydrodynamical simulations with $M_ν= 0.06$eV. This offers the first non-linear, theoretically general means of accurately including massive neutrinos for beyond-$Λ$CDM cosmologies, and further suggests that baryonic, massive neutrino and dark energy physics can be reliably modelled independently.
△ Less
Submitted 25 September, 2021; v1 submitted 25 May, 2021;
originally announced May 2021.
-
Investigating the degeneracy between modified gravity and massive neutrinos with redshift-space distortions
Authors:
Bill S. Wright,
Kazuya Koyama,
Hans A. Winther,
Gong-Bo Zhao
Abstract:
There is a well known degeneracy between the enhancement of the growth of large-scale structure produced by modified gravity models and the suppression due to the free-streaming of massive neutrinos at late times. This makes the matter power-spectrum alone a poor probe to distinguish between modified gravity and the concordance $Λ$CDM model when neutrino masses are not strongly constrained. In thi…
▽ More
There is a well known degeneracy between the enhancement of the growth of large-scale structure produced by modified gravity models and the suppression due to the free-streaming of massive neutrinos at late times. This makes the matter power-spectrum alone a poor probe to distinguish between modified gravity and the concordance $Λ$CDM model when neutrino masses are not strongly constrained. In this work, we investigate the potential of using redshift-space distortions (RSD) to break this degeneracy when the modification to gravity is scale-dependent in the form of Hu-Sawicki $f(R)$. We find that if the linear growth rate can be recovered from the RSD signal, the degeneracy can be broken at the level of the dark matter field. However, this requires accurate modelling of the non-linearities in the RSD signal, and we here present an extension of the standard perturbation theory-based model for non-linear RSD that includes both Hu-Sawicki $f(R)$ modified gravity and massive neutrinos.
△ Less
Submitted 3 August, 2020; v1 submitted 27 February, 2019;
originally announced February 2019.
-
Constraining the time variation of Newton's constant $G$ with gravitational-wave standard sirens and supernovae
Authors:
Wen Zhao,
Bill S. Wright,
Baojiu Li
Abstract:
The intrinsic peak luminosity of Type Ia supernovae (SNIa) depends on the value of Newton's gravitational constant $G$, through the Chandrasekhar mass $M_{\rm Ch}\propto G^{-3/2}$. If the luminosity distance can be independently determined, the SNIa can be treated as a tracker to constrain the possible time variation of $G$ in different redshift ranges. The gravitational-wave (GW) standard sirens,…
▽ More
The intrinsic peak luminosity of Type Ia supernovae (SNIa) depends on the value of Newton's gravitational constant $G$, through the Chandrasekhar mass $M_{\rm Ch}\propto G^{-3/2}$. If the luminosity distance can be independently determined, the SNIa can be treated as a tracker to constrain the possible time variation of $G$ in different redshift ranges. The gravitational-wave (GW) standard sirens, caused by the coalescence of binary neutron stars, provide a model-independent way to measure the distance of GW events, which can be used to determine the luminosity distances of SNIa by interpolation, provided the GW and SNIa samples have similar redshift ranges. We demonstrate that combining the GW observations of third-generation detectors with SNIa data provides a powerful and model-independent way to measure $G$ in a wide redshift range, which can constrain the ratio $G/G_0$, where $G$ and $G_0$ are respectively the values in the redshift ranges $z>0.1$ and $z<0.1$, at the level of $1.5\%$.
△ Less
Submitted 18 October, 2018; v1 submitted 9 April, 2018;
originally announced April 2018.
-
Type Ia supernovae, standardisable candles, and gravity
Authors:
Bill S. Wright,
Baojiu Li
Abstract:
Type Ia supernovae (SNIe) are generally accepted to act as standardisable candles, and their use in cosmology led to the first confirmation of the as yet unexplained accelerated cosmic expansion. Many of the theoretical models to explain the cosmic acceleration assume modifications to Einsteinian General Relativity which accelerate the expansion, but the question of whether such modifications also…
▽ More
Type Ia supernovae (SNIe) are generally accepted to act as standardisable candles, and their use in cosmology led to the first confirmation of the as yet unexplained accelerated cosmic expansion. Many of the theoretical models to explain the cosmic acceleration assume modifications to Einsteinian General Relativity which accelerate the expansion, but the question of whether such modifications also affect the ability of SNIe to be standardisable candles has rarely been addressed. This paper is an attempt to answer this question. For this we adopt a semi-analytical model to calculate SNIe light curves in non-standard gravity. We use this model to show that the average rescaled intrinsic peak luminosity -- a quantity that is assumed to be constant with redshift in standard analyses of Type Ia supernova (SNIa) cosmology data -- depends on the strength of gravity in the supernova's local environment because the latter determines the Chandrasekhar mass -- the mass of the SNIa's white dwarf progenitor right before the explosion. This means that SNIe are no longer standardisable candles in scenarios where the strength of gravity evolves over time, and therefore the cosmology implied by the existing SNIa data will be different when analysed in the context of such models. As an example, we show that the observational SNIa cosmology data can be fitted with both a model where $(Ω_{\rm M}, Ω_Λ)=(0.62, 0.38)$ and Newton's constant $G$ varies as $G(z)=G_0(1+z)^{-1/4}$ and the standard model where $(Ω_{\rm M}, Ω_Λ)=(0.3, 0.7)$ and $G$ is constant, when the Universe is assumed to be flat.
△ Less
Submitted 29 October, 2018; v1 submitted 19 October, 2017;
originally announced October 2017.
-
COLA with massive neutrinos
Authors:
Bill S. Wright,
Hans A. Winther,
Kazuya Koyama
Abstract:
The effect of massive neutrinos on the growth of cold dark matter perturbations acts as a scale-dependent Newton's constant and leads to scale-dependent growth factors just as we often find in models of gravity beyond General Relativity. We show how to compute growth factors for $Λ$CDM and general modified gravity cosmologies combined with massive neutrinos in Lagrangian perturbation theory for us…
▽ More
The effect of massive neutrinos on the growth of cold dark matter perturbations acts as a scale-dependent Newton's constant and leads to scale-dependent growth factors just as we often find in models of gravity beyond General Relativity. We show how to compute growth factors for $Λ$CDM and general modified gravity cosmologies combined with massive neutrinos in Lagrangian perturbation theory for use in COLA and extensions thereof. We implement this together with the grid-based massive neutrino method of Brandbyge and Hannestad in $\texttt{MG-PICOLA}$ and compare COLA simulations to full $\it N$-body simulations of $Λ$CDM and $f(R)$ gravity with massive neutrinos. Our implementation is computationally cheap if the underlying cosmology already has scale-dependent growth factors and it is shown to be able to produce results that match $\it N$-body to percent level accuracy for both the total and CDM matter power-spectra up to $k\lesssim 1 h/$Mpc.
△ Less
Submitted 19 October, 2017; v1 submitted 23 May, 2017;
originally announced May 2017.
-
COLA with scale-dependent growth: applications to screened modified gravity models
Authors:
Hans A. Winther,
Kazuya Koyama,
Marc Manera,
Bill S. Wright,
Gong-Bo Zhao
Abstract:
We present a general parallelized and easy-to-use code to perform numerical simulations of structure formation using the COLA (COmoving Lagrangian Acceleration) method for cosmological models that exhibit scale-dependent growth at the level of first and second order Lagrangian perturbation theory. For modified gravity theories we also include screening using a fast approximate method that covers a…
▽ More
We present a general parallelized and easy-to-use code to perform numerical simulations of structure formation using the COLA (COmoving Lagrangian Acceleration) method for cosmological models that exhibit scale-dependent growth at the level of first and second order Lagrangian perturbation theory. For modified gravity theories we also include screening using a fast approximate method that covers all the main examples of screening mechanisms in the literature. We test the code by comparing it to full simulations of two popular modified gravity models, namely $f(R)$ gravity and nDGP, and find good agreement in the modified gravity boost-factors relative to $Λ$CDM even when using a fairly small number of COLA time steps.
△ Less
Submitted 19 July, 2017; v1 submitted 2 March, 2017;
originally announced March 2017.