scholarly journals Fast numerical method to generate halo catalogues in modified gravity (part I): second-order Lagrangian perturbation theory

2020 ◽  
Vol 493 (1) ◽  
pp. 1153-1164 ◽  
Author(s):  
C Moretti ◽  
S Mozzon ◽  
P Monaco ◽  
E Munari ◽  
M Baldi

ABSTRACT We present a new numerical method to determine second-order Lagrangian displacement fields in presence of modified gravity (MG). We start from the extension of Lagrangian perturbation theory (LPT) to a class of MG models, which can be described by a parametrized Poisson equation. We exploit Fast Fourier transforms to compute the full source term of the differential equation for the second-order Lagrangian displacement field. We compare its mean to the source term computed for specific configurations, for which a k-dependent solution can be found numerically. We choose the configuration that best matches the full source term, thus obtaining an approximate factorization of the second-order displacement field as the space term valid for Λ Cold Dark Matter (ΛCDM) times a k-dependent, second-order growth factor. Such approximation is used to compute second-order displacements for particles. The method is tested against N-body simulations run with standard and f(R) gravity: we rely on the results of a friends-of-friends code run on the N-body snapshots to assign particles to haloes, then compute the halo power spectrum. We find very consistent results for the two gravity theories: second-order LPT (2LPT) allows to recover the N-body halo power spectrum within ∼10 per cent precision to k ∼ 0.2–0.4 h Mpc−1, as well as halo positions. We show that the performance of 2LPT with MG is the same (within 1 per cent) as the one obtained for standard ΛCDM case. This formulation of 2LPT can quickly generate dark matter distributions with f(R) gravity, and can easily be extended to other MG theories.

2020 ◽  
Vol 501 (1) ◽  
pp. L71-L75
Author(s):  
Cornelius Rampf ◽  
Oliver Hahn

ABSTRACT Perturbation theory is an indispensable tool for studying the cosmic large-scale structure, and establishing its limits is therefore of utmost importance. One crucial limitation of perturbation theory is shell-crossing, which is the instance when cold-dark-matter trajectories intersect for the first time. We investigate Lagrangian perturbation theory (LPT) at very high orders in the vicinity of the first shell-crossing for random initial data in a realistic three-dimensional Universe. For this, we have numerically implemented the all-order recursion relations for the matter trajectories, from which the convergence of the LPT series at shell-crossing is established. Convergence studies performed at large orders reveal the nature of the convergence-limiting singularities. These singularities are not the well-known density singularities at shell-crossing but occur at later times when LPT already ceased to provide physically meaningful results.


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 ◽  
Vol 650 ◽  
pp. A113
Author(s):  
Margot M. Brouwer ◽  
Kyle A. Oman ◽  
Edwin A. Valentijn ◽  
Maciej Bilicki ◽  
Catherine Heymans ◽  
...  

We present measurements of the radial gravitational acceleration around isolated galaxies, comparing the expected gravitational acceleration given the baryonic matter (gbar) with the observed gravitational acceleration (gobs), using weak lensing measurements from the fourth data release of the Kilo-Degree Survey (KiDS-1000). These measurements extend the radial acceleration relation (RAR), traditionally measured using galaxy rotation curves, by 2 decades in gobs into the low-acceleration regime beyond the outskirts of the observable galaxy. We compare our RAR measurements to the predictions of two modified gravity (MG) theories: modified Newtonian dynamics and Verlinde’s emergent gravity (EG). We find that the measured relation between gobs and gbar agrees well with the MG predictions. In addition, we find a difference of at least 6σ between the RARs of early- and late-type galaxies (split by Sérsic index and u − r colour) with the same stellar mass. Current MG theories involve a gravity modification that is independent of other galaxy properties, which would be unable to explain this behaviour, although the EG theory is still limited to spherically symmetric static mass models. The difference might be explained if only the early-type galaxies have significant (Mgas ≈ M⋆) circumgalactic gaseous haloes. The observed behaviour is also expected in Λ-cold dark matter (ΛCDM) models where the galaxy-to-halo mass relation depends on the galaxy formation history. We find that MICE, a ΛCDM simulation with hybrid halo occupation distribution modelling and abundance matching, reproduces the observed RAR but significantly differs from BAHAMAS, a hydrodynamical cosmological galaxy formation simulation. Our results are sensitive to the amount of circumgalactic gas; current observational constraints indicate that the resulting corrections are likely moderate. Measurements of the lensing RAR with future cosmological surveys (such as Euclid) will be able to further distinguish between MG and ΛCDM models if systematic uncertainties in the baryonic mass distribution around galaxies are reduced.


2019 ◽  
Vol 485 (4) ◽  
pp. 5474-5489 ◽  
Author(s):  
Mark R Lovell ◽  
Jesús Zavala ◽  
Mark Vogelsberger

Abstract A cut-off in the linear matter power spectrum at dwarf galaxy scales has been shown to affect the abundance, formation mechanism and age of dwarf haloes, and their galaxies at high and low redshifts. We use hydrodynamical simulations of galaxy formation within the ETHOS framework in a benchmark model that has such a cut-off and that has been shown to be an alternative to the cold dark matter (CDM) model that alleviates its dwarf-scale challenges. We show how galaxies in this model form differently to CDM, on a halo-by-halo basis, at redshifts z ≥ 6. We show that when CDM haloes with masses around the ETHOS half-mode mass scale are resimulated with the ETHOS matter power spectrum, they form with 50 per cent less mass than their CDM counterparts due to their later formation times, yet they retain more of their gas reservoir due to the different behaviour of gas and dark matter during the monolithic collapse of the first haloes in models with a galactic-scale cut-off. As a result, galaxies in ETHOS haloes near the cut-off scale grow rapidly between z = 10 and 6 and by z = 6 end up having very similar stellar masses, higher gas fractions and higher star formation rates relative to their CDM counterparts. We highlight these differences by making predictions for how the number of galaxies with old stellar populations is suppressed in ETHOS for both z = 6 galaxies and for gas-poor Local Group fossil galaxies. Interestingly, we find an age gradient in ETHOS between galaxies that form in high- and low-density environments.


1996 ◽  
Vol 466 ◽  
pp. 13 ◽  
Author(s):  
Anatoly Klypin ◽  
Joel Primack ◽  
Jon Holtzman

1993 ◽  
Vol 415 ◽  
pp. L67 ◽  
Author(s):  
Tereasa G. Brainerd ◽  
Jens V. Villumsen

1987 ◽  
Vol 117 ◽  
pp. 367-367
Author(s):  
Rosemary F. G. Wyse ◽  
Bernard J. T. Jones

We present a simple model for the formation of elliptical galaxies, based on a binary clustering hierarchy of dark matter, the chemical enrichment of the gas at each level being controlled by supernovae. The initial conditions for the non-linear phases of galaxy formation are set by the post-recombination power spectrum of density fluctuations. We investigate two models for this power spectrum - the first is a straightforward power law, |δk|2 ∝ kn, and the second is Peeble's analytic approximation to the emergent spectrum in a universe dominated by cold dark matter. The normalisation is chosen such that on some scale, say M ∼ 1012M⊙, the objects that condense out have properties - radius and velocity dispersion - resembling ‘typical’ galaxies. There is some ambiguity in this due to the poorly determined mass-to-light ratio of a typical elliptical galaxy — we look at two normalisations, σ1D ∼ 350kms−1 and σ1D ∼ 140kms−1. The choice determines which of Compton cooling or hydrogen cooling is more important during the galaxy formation period. The non-linear behaviour of the perturbations is treated by the homogeneous sphere approximation.


Sign in / Sign up

Export Citation Format

Share Document