scholarly journals N-body simulations for parametrized modified gravity

2020 ◽  
Vol 497 (2) ◽  
pp. 1885-1894
Author(s):  
Farbod Hassani ◽  
Lucas Lombriser

ABSTRACT We present MG-evolution, an N-body code simulating the cosmological structure formation for parametrized modifications of gravity. It is built from the combination of parametrized linear theory with a parametrization of the deeply non-linear cosmological regime extrapolated from modified spherical collapse computations that cover the range of known screening mechanisms. We test MG-evolution, which runs at the speed of conventional ΛCDM simulations, against a suit of existing exact model-specific codes, encompassing linearized and chameleon f(R) gravity as well as the normal branch of the Dvali–Gabadadz–Porrati braneworld model, hence covering both large-field value and large-derivative screening effects. We compare the non-linear power spectra produced by the parametrized and model-specific approaches over the full range of scales set by the box size and resolution of our simulations, k = (0.05 − 2.5) $h\, \mathrm{Mpc}^{-1}$, and for two redshift slices, z = 0 and z = 1. We find sub-percent to one-percent level recovery of all the power spectra generated with the model-specific codes for the full range of scales. MG-evolution can be used for generalized and accurate tests of gravity and dark energy with the increasing wealth of high-precision cosmological survey data becoming available over the next decade.

Author(s):  
Myles A Mitchell ◽  
Christian Arnold ◽  
César Hernández-Aguayo ◽  
Baojiu Li

Abstract We study the effects of two popular modified gravity theories, which incorporate very different screening mechanisms, on the angular power spectra of the thermal (tSZ) and kinematic (kSZ) components of the Sunyaev-Zel’dovich effect. Using the first cosmological simulations that simultaneously incorporate both screened modified gravity and a complete galaxy formation model, we find that the tSZ and kSZ power spectra are significantly enhanced by the strengthened gravitational forces in Hu-Sawicki f(R) gravity and the normal-branch Dvali-Gabadadze-Porrati model. Employing a combination of non-radiative and full-physics simulations, we find that the extra baryonic physics present in the latter acts to suppress the tSZ power on angular scales l ≳ 3000 and the kSZ power on all tested scales, and this is found to have a substantial effect on the model differences. Our results indicate that the tSZ and kSZ power can be used as powerful probes of gravity on large scales, using data from current and upcoming surveys, provided sufficient work is conducted to understand the sensitivity of the constraints to baryonic processes that are currently not fully understood.


2021 ◽  
Vol 503 (4) ◽  
pp. 5638-5645
Author(s):  
Gábor Rácz ◽  
István Szapudi ◽  
István Csabai ◽  
László Dobos

ABSTRACT The classical gravitational force on a torus is anisotropic and always lower than Newton’s 1/r2 law. We demonstrate the effects of periodicity in dark matter only N-body simulations of spherical collapse and standard Lambda cold dark matter (ΛCDM) initial conditions. Periodic boundary conditions cause an overall negative and anisotropic bias in cosmological simulations of cosmic structure formation. The lower amplitude of power spectra of small periodic simulations is a consequence of the missing large-scale modes and the equally important smaller periodic forces. The effect is most significant when the largest mildly non-linear scales are comparable to the linear size of the simulation box, as often is the case for high-resolution hydrodynamical simulations. Spherical collapse morphs into a shape similar to an octahedron. The anisotropic growth distorts the large-scale ΛCDM dark matter structures. We introduce the direction-dependent power spectrum invariant under the octahedral group of the simulation volume and show that the results break spherical symmetry.


2021 ◽  
Author(s):  
Rui feng Zheng ◽  
Jia ming Shi ◽  
Taotao Qiu

Abstract It is well known that primordial black hole (PBH) can be generated in inflation process of the early universe, especially when the inflaton field has some non-trivial features that could break the slow-roll condition. In this paper, we investigate a toy model of inflation with bumpy potential, which has one or several bumps. We found that potential with multi-bump can give rise to power spectra with multi peaks in small-scale region, which can in turn predict the generation of primordial black holes in various mass ranges. We also consider the two possibilities of PBH formation by spherical collapse and elliptical collapse. And discusses the scalar-induced gravitational waves (SIGWs) generated by the second-order scalar perturbations.


2021 ◽  
Vol 10 (6) ◽  
Author(s):  
Matthias Bartelmann ◽  
Johannes Dombrowski ◽  
Sara Konrad ◽  
Elena Kozlikin ◽  
Robert Lilow ◽  
...  

We use the recently developed Kinetic Field Theory (KFT) for cosmic structure formation to show how non-linear power spectra for cosmic density fluctuations can be calculated in a mean-field approximation to the particle interactions. Our main result is a simple, closed and analytic, approximate expression for this power spectrum. This expression has two parameters characterising non-linear structure growth which can be calibrated within KFT itself. Using this self-calibration, the non-linear power spectrum agrees with results obtained from numerical simulations to within typically \lesssim10\,\%≲10% up to wave numbers k\lesssim10\,h\,\mathrm{Mpc}^{-1}k≲10hMpc−1 at redshift z = 0z=0. Adjusting the two parameters to optimise agreement with numerical simulations, the relative difference to numerical results shrinks to typically \lesssim 5\,\%≲5%. As part of the derivation of our mean-field approximation, we show that the effective interaction potential between dark-matter particles relative to Zel’dovich trajectories is sourced by non-linear cosmic density fluctuations only, and is approximately of Yukawa rather than Newtonian shape.


1996 ◽  
Vol 282 (3) ◽  
pp. 767-778 ◽  
Author(s):  
A. N. Taylor ◽  
A. J. S. Hamilton
Keyword(s):  

2021 ◽  
Author(s):  
Soham Sheth ◽  
Francois McKee ◽  
Kieran Neylon ◽  
Ghazala Fazil

Abstract We present a novel reservoir simulator time-step selection approach which uses machine-learning (ML) techniques to analyze the mathematical and physical state of the system and predict time-step sizes which are large while still being efficient to solve, thus making the simulation faster. An optimal time-step choice avoids wasted non-linear and linear equation set-up work when the time-step is too small and avoids highly non-linear systems that take many iterations to solve. Typical time-step selectors use a limited set of features to heuristically predict the size of the next time-step. While they have been effective for simple simulation models, as model complexity increases, there is an increasing need for robust data-driven time-step selection algorithms. We propose two workflows – static and dynamic – that use a diverse set of physical (e.g., well data) and mathematical (e.g., CFL) features to build a predictive ML model. This can be pre-trained or dynamically trained to generate an inference model. The trained model can also be reinforced as new data becomes available and efficiently used for transfer learning. We present the application of these workflows in a commercial reservoir simulator using distinct types of simulation model including black oil, compositional and thermal steam-assisted gravity drainage (SAGD). We have found that history-match and uncertainty/optimization studies benefit most from the static approach while the dynamic approach produces optimum step-sizes for prediction studies. We use a confidence monitor to manage the ML time-step selector at runtime. If the confidence level falls below a threshold, we switch to traditional heuristic method for that time-step. This avoids any degradation in the performance when the model features are outside the training space. Application to several complex cases, including a large field study, shows a significant speedup for single simulations and even better results for multiple simulations. We demonstrate that any simulation can take advantage of the stored state of the trained model and even augment it when new situations are encountered, so the system becomes more effective as it is exposed to more data.


2018 ◽  
Vol 27 (04) ◽  
pp. 1830003 ◽  
Author(s):  
David F. Mota

We review the effects of modified gravity theories, in particular the symmetron and [Formula: see text] gravity, on the nonlinear regime of structure formation. In particular, we investigate the velocity dispersion of galaxy clusters as a function of the halo masses, how the matter power spectra depend on the coupling, range, and screening scale of the fifth force, and on possible ways of detecting violations of the equivalence principle using the mass inferred via lensing methods versus the mass inferred via dynamical methods. Furthermore, we show how one could use different voids statistics as one of the most promising probes of modified gravity.


2019 ◽  
Vol 488 (2) ◽  
pp. 2121-2142 ◽  
Author(s):  
M Cataneo ◽  
L Lombriser ◽  
C Heymans ◽  
A J Mead ◽  
A Barreira ◽  
...  

ABSTRACT We present a general method to compute the non-linear matter power spectrum for dark energy (DE) and modified gravity scenarios with per cent-level accuracy. By adopting the halo model and non-linear perturbation theory, we predict the reaction of a lambda cold dark matter (ΛCDM) matter power spectrum to the physics of an extended cosmological parameter space. By comparing our predictions to N-body simulations we demonstrate that with no-free parameters we can recover the non-linear matter power spectrum for a wide range of different w0–wa DE models to better than 1 per cent accuracy out to k ≈ 1 $h \,{\rm Mpc}^{-1}$. We obtain a similar performance for both DGP and f(R) gravity, with the non-linear matter power spectrum predicted to better than 3 per cent accuracy over the same range of scales. When including direct measurements of the halo mass function from the simulations, this accuracy improves to 1 per cent. With a single suite of standard ΛCDM N-body simulations, our methodology provides a direct route to constrain a wide range of non-standard extensions to the concordance cosmology in the high signal-to-noise non-linear regime.


Sign in / Sign up

Export Citation Format

Share Document