Fast generation of mock galaxy catalogues in modified gravity models with COLA

Fiorini, Bartolomeo; Koyama, Kazuya; Izard, Albert; Winther, Hans A.; Wright, Bill S.; Li, Baojiu

doi:10.48550/arxiv.2106.05197

Cited by 4 publications

(3 citation statements)

References 51 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…In spite of such a massive improvement in speed over those latter codes, it is worthwhile to note that mg-glam is not an approximate code: it solves the full Poisson and MG equations, and its accuracy is only limited by the resolution of the PM grid used, which can be specified by users based on their particular scientific objectives. This makes it different from fast approximate simulation codes such as those [105][106][107][108] based on the COmoving Lagrangian Acceleration method (cola) [109].…”

Section: Introductionmentioning

confidence: 99%

Fast full N-body simulations of generic modified gravity: derivative coupling models

Hernández‐Aguayo

Ruan

et al. 2022

J. Cosmol. Astropart. Phys.

View full text Add to dashboard Cite

We present mg-glam, a code developed for the very fast production of full N-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for two broad classes of MG models with derivative coupling terms — the Vainshtein- and Kmouflage-type models — which respectively features the Vainshtein and Kmouflage screening mechanism. Derived from the parallel particle-mesh code glam, mg-glam incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For Kmouflage, we have proposed a new algorithm for the relaxation solver, and run the first simulations of the model to understand its cosmological behaviour. In a companion paper, we describe versions of this code developed for conformally-coupled MG models, including several variants of f(R) gravity, the symmetron model and coupled quintessence. Altogether, mg-glam has so far implemented the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous (over two orders of magnitude) speedup when comparing its running time with earlier N-body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance. mg-glam is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and accurate emulators for ongoing and future galaxy surveys.

show abstract

Section: Introductionmentioning

confidence: 99%

Fast full N-body simulations of generic modified gravity: derivative coupling models

Hernández‐Aguayo

Ruan

et al. 2022

J. Cosmol. Astropart. Phys.

View full text Add to dashboard Cite

show abstract

“…In spite of such a massive improvement in speed over those latter codes, it is worthwhile to note that -is not an approximate code: it solves the full Poisson and MG equations, and its accuracy is only limited by the resolution of the PM grid used, which can be specified by users based on their particular scientific objectives. This makes it different from fast approximate simulation codes such as those [103][104][105][106] based on the COmoving Lagrangian Acceleration method ( ) [107]. This paper is organised as follows.…”

Section: Introductionmentioning

confidence: 99%

Fast full N-body simulations of generic modified gravity: conformal coupling models

Ruan,

Hernández-Aguayo,

et al. 2021

Preprint

Self Cite

View full text Add to dashboard Cite

We present -, a code developed for the very fast production of full N -body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for three classes of MG models with conformal coupling terms: the f (R) gravity, symmetron and coupled quintessence models. Derived from the parallel particle-mesh code , -incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For f (R) gravity, we have included new variants to diversify the model behaviour, and we have tailored the relaxation algorithms to these to maintain high computational efficiency. In a companion paper, we describe versions of this code developed for derivative coupling MG models, including the Vainshteinand K-mouflage-type models.can model the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous speedup of a factor of more than a hundred compared with earlier N -body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance.is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and the production of accurate emulators for ongoing and future galaxy surveys.

show abstract

“…Specifically, in this work, we will consider the first of these examples, the COLA approach. This method has been well-validated, and is known to give a good agreement on quasi non-linear scales in ΛCDM and beyond-ΛCDM cosmologies when comparing its prediction for the matter power spectrum to the ones from full N-body simulations [12][13][14][15][16]. Additionally, a new avenue using the COLA method was presented in [17].…”

Section: Introductionmentioning

confidence: 99%

Enabling matter power spectrum emulation in beyond-$Λ$CDM cosmologies with COLA

Brando¹,

Fiorini²,

Koyama³

et al. 2022

Preprint

Self Cite

View full text Add to dashboard Cite

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 ≤ z ≤ 3, we find that the COLA method is in excellent agreement with the two emulators for scales up to k ∼ 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 ≤ 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.

show abstract

Fast generation of mock galaxy catalogues in modified gravity models with COLA

Cited by 4 publications

References 51 publications

Fast full N-body simulations of generic modified gravity: derivative coupling models

Fast full N-body simulations of generic modified gravity: derivative coupling models

Fast full N-body simulations of generic modified gravity: conformal coupling models

Enabling matter power spectrum emulation in beyond-$Λ$CDM cosmologies with COLA

Contact Info

Product

Resources

About