scholarly journals Anisotropic separate universe simulations

2020 ◽  
Vol 496 (1) ◽  
pp. 483-496 ◽  
Author(s):  
Shogo Masaki ◽  
Takahiro Nishimichi ◽  
Masahiro Takada

ABSTRACT The long-wavelength coherent overdensity and tidal force, which are not direct observables for a finite-volume survey, affect time evolution of cosmic structure formation and therefore clustering observables through the mode coupling. In this paper, we develop an ‘anisotropic’ separate universe (SU) simulation technique to simulate large-scale structure formation taking into account the effect of large-scale tidal force into the anisotropic expansion of local background. We modify the treepmN-body simulation code to implement the anisotropic SU simulations, and then study the ‘response’ function of matter power spectrum that describes how the matter power spectrum responds to the large-scale tidal effect as a function of wavenumber and redshift for a given global cosmology. We test and validate the SU simulation results from the comparison with the perturbation theory predictions and the results from high-resolution particle-mesh simulations. We find that the response function displays characteristic scale dependencies over the range of scales down to non-linear scales, up to k ≃ 6 h Mpc−1.

Author(s):  
Jens Stücker ◽  
Andreas S Schmidt ◽  
Simon D M White ◽  
Fabian Schmidt ◽  
Oliver Hahn

Abstract We present anisotropic ‘separate universe’ simulations which modify the N-body code gadget4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber k = 7 h Mpc−1. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value at z = 0, RK ≈ 0.5 for k ≥ 3 h Mpc−1 up to the scale where we consider our simulations reliable, k ≤ 7 h Mpc−1. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.


2020 ◽  
Vol 500 (1) ◽  
pp. 1018-1028
Author(s):  
Shogo Masaki ◽  
Takahiro Nishimichi ◽  
Masahiro Takada

ABSTRACT To generate initial conditions for cosmological N-body simulations, one needs to prepare a uniform distribution of simulation particles, the so-called pre-initial condition (pre-IC). The standard method to construct the pre-IC is to place the particles on the lattice grids evenly spaced in the three-dimensional spatial coordinates. However, even after the initial displacement of each particle according to cosmological perturbations, the particle distribution remains to display an artificial anisotropy. Such an artefact causes systematic effects in simulations at later time until the evolved particle distribution sufficiently erases the initial anisotropy. In this paper, we study the impacts of the pre-IC on the anisotropic separate universe simulation, where the effect of large-scale tidal field on structure formation is taken into account using the anisotropic expansion in a local background (simulation volume). To quantify the impacts, we compare the simulations employing the standard grid pre-IC and the glass one, where the latter is supposed to suppress the initial anisotropy. We show that the artificial features in the grid pre-IC simulations are seen until z ∼ 9, while the glass pre-IC simulations appear to be stable and accurate over the range of scales we study. From these results we find that a coupling of the large-scale tidal field with matter clustering is enhanced compared to the leading-order prediction of perturbation theory in the quasi-non-linear regime in the redshift range 5 ≲ z ≲ 15, indicating the importance of tidal field on structure formation at such high redshifts, e.g. during the epoch of reionization.


2020 ◽  
Vol 498 (2) ◽  
pp. 2887-2911 ◽  
Author(s):  
Simon Foreman ◽  
William Coulton ◽  
Francisco Villaescusa-Navarro ◽  
Alexandre Barreira

ABSTRACT The large-scale clustering of matter is impacted by baryonic physics, particularly active galactic nuclei (AGN) feedback. Modelling or mitigating this impact will be essential for making full use of upcoming measurements of cosmic shear and other large-scale structure probes. We study baryonic effects on the matter bispectrum, using measurements from a selection of state-of-the-art hydrodynamical simulations: IllustrisTNG, Illustris, EAGLE, and BAHAMAS. We identify a low-redshift enhancement of the bispectrum, peaking at $k\sim 3\,h\, {\rm Mpc}^{-1}$, which is present in several simulations, and discuss how it can be associated to the evolving nature of AGN feedback at late times. This enhancement does not appear in the matter power spectrum, and therefore represents a new source of degeneracy breaking between two- and three-point statistics. In addition, we provide physical interpretations for other aspects of these measurements, and make initial comparisons to predictions from perturbation theory, empirical fitting formulas, and the response function formalism. We publicly release our measurements (including estimates of their uncertainty due to sample variance) and bispectrum measurement code as resources for the community.


2020 ◽  
Vol 498 (3) ◽  
pp. 3403-3419
Author(s):  
Sebastian Bohr ◽  
Jesús Zavala ◽  
Francis-Yan Cyr-Racine ◽  
Mark Vogelsberger ◽  
Torsten Bringmann ◽  
...  

ABSTRACT We propose two effective parameters that fully characterize galactic-scale structure formation at high redshifts (z ≳ 5) for a variety of dark matter (DM) models that have a primordial cutoff in the matter power spectrum. Our description is within the recently proposed ETHOS framework and includes standard thermal warm DM (WDM) and models with dark acoustic oscillations (DAOs). To define and explore this parameter space, we use high-redshift zoom-in simulations that cover a wide range of non-linear scales from those where DM should behave as CDM (k ∼ 10 h Mpc−1), down to those characterized by the onset of galaxy formation (k ∼ 500 h Mpc−1). We show that the two physically motivated parameters hpeak and kpeak, the amplitude and scale of the first DAO peak, respectively, are sufficient to parametrize the linear matter power spectrum and classify the DM models as belonging to effective non-linear structure formation regions. These are defined by their relative departure from cold DM (kpeak → ∞) and WDM (hpeak = 0) according to the non-linear matter power spectrum and halo mass function. We identify a region where the DAOs still leave a distinct signature from WDM down to z = 5, while a large part of the DAO parameter space is shown to be degenerate with WDM. Our framework can then be used to seamlessly connect a broad class of particle DM models to their structure formation properties at high redshift without the need of additional N-body simulations.


2020 ◽  
Vol 499 (2) ◽  
pp. 1769-1787
Author(s):  
Anaëlle Halle ◽  
Takahiro Nishimichi ◽  
Atsushi Taruya ◽  
Stéphane Colombi ◽  
Francis Bernardeau

ABSTRACT The power spectrum response function of the large-scale structure of the Universe describes how the evolved power spectrum is modified by a small change in initial power through non-linear mode coupling of gravitational evolution. It was previously found that the response function for the coupling from small to large scales is strongly suppressed in amplitude, especially at late times, compared to predictions from perturbation theory (PT) based on the single-stream approximation. One obvious explanation for this is that PT fails to describe the dynamics beyond shell crossing. We test this idea by comparing measurements in N-body simulations to prescriptions based on PT but augmented with adaptive smoothing to account for the formation of non-linear structures of various sizes in the multistream regime. We first start with one-dimensional (1D) cosmology, where the Zel’dovich approximation provides the exact solution in the single-stream regime. Similarly to the three-dimensional (3D) case, the response function of the large-scale modes exhibits a strong suppression in amplitude at small scales that cannot be explained by the Zel’dovich solution alone. However, by performing adaptive smoothing of initial conditions to identify haloes of different sizes and solving approximately post-collapse dynamics in the three-stream regime, agreement between theory and simulations drastically improves. We extend our analyses to the 3D case using the pinocchio algorithm, in which similar adaptive smoothing is implemented on the Lagrangian PT fields to identify haloes and is combined with a spherical halo prescription to account for post-collapse dynamics. Again, a suppression is found in the coupling between small- and large-scale modes and the agreement with simulations is improved.


2021 ◽  
Vol 2021 (12) ◽  
pp. 044
Author(s):  
G. Parimbelli ◽  
G. Scelfo ◽  
S.K. Giri ◽  
A. Schneider ◽  
M. Archidiacono ◽  
...  

Abstract We investigate and quantify the impact of mixed (cold and warm) dark matter models on large-scale structure observables. In this scenario, dark matter comes in two phases, a cold one (CDM) and a warm one (WDM): the presence of the latter causes a suppression in the matter power spectrum which is allowed by current constraints and may be detected in present-day and upcoming surveys. We run a large set of N-body simulations in order to build an efficient and accurate emulator to predict the aforementioned suppression with percent precision over a wide range of values for the WDM mass, Mwdm, and its fraction with respect to the totality of dark matter, fwdm. The suppression in the matter power spectrum is found to be independent of changes in the cosmological parameters at the 2% level for k≲ 10 h/Mpc and z≤ 3.5. In the same ranges, by applying a baryonification procedure on both ΛCDM and CWDM simulations to account for the effect of feedback, we find a similar level of agreement between the two scenarios. We examine the impact that such suppression has on weak lensing and angular galaxy clustering power spectra. Finally, we discuss the impact of mixed dark matter on the shape of the halo mass function and which analytical prescription yields the best agreement with simulations. We provide the reader with an application to galaxy cluster number counts.


2020 ◽  
Vol 501 (1) ◽  
pp. 833-852
Author(s):  
Toshiki Kurita ◽  
Masahiro Takada ◽  
Takahiro Nishimichi ◽  
Ryuichi Takahashi ◽  
Ken Osato ◽  
...  

ABSTRACT We use a suite of N-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the ΛCDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the E-mode field of halo shapes with the matter/halo distribution, $P_{\delta E}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those of the auto-power spectrum of the E-mode, $P^{(\ell)}_{EE}(k)$, based on the E/B-mode decomposition. The IA power spectra have non-vanishing amplitudes over the linear to non-linear scales, and the large-scale amplitudes at k ≲ 0.1 h−1 Mpc are related to the matter power spectrum via a constant coefficient (AIA), similar to the linear bias parameter of galaxy or halo density field. We find that the cross- and auto-power spectra PδE and PEE at non-linear scales, k ≳ 0.1 h−1 Mpc, show different k-dependences relative to the matter power spectrum, suggesting a violation of the non-linear alignment model commonly used to model contaminations of cosmic shear signals. The IA power spectra exhibit baryon acoustic oscillations, and vary with halo samples of different masses, redshifts, and cosmological parameters (Ωm, S8). The cumulative signal-to-noise ratio for the IA power spectra is about 60 per cent of that for the halo density power spectrum, where the super-sample covariance is found to give a significant contribution to the total covariance. Thus our results demonstrate that the IA power spectra of galaxy shapes, measured from imaging and spectroscopic surveys for an overlapping area of the sky, can be used to probe the underlying matter power spectrum, the primordial curvature perturbations, and cosmological parameters, in addition to the standard galaxy density power spectrum.


Author(s):  
D V Gomez-Navarro ◽  
A J Mead ◽  
A Aviles ◽  
A de la Macorra

Abstract Some beyond ΛCDM cosmological models have dark-sector energy densities that suffer phase transitions. Fluctuations entering the horizon during such a transition can receive enhancements that ultimately show up as a distinctive bump in the power spectrum relative to a model with no phase transition. In this work, we study the non-linear evolution of such signatures in the matter power spectrum and correlation function using N-body simulations, perturbation theory and hmcode- a halo-model based method. We focus on modelling the response, computed as the ratio of statistics between a model containing a bump and one without it, rather than in the statistics themselves. Instead of working with a specific theoretical model, we inject a parametric family of Gaussian bumps into otherwise standard ΛCDM spectra. We find that even when the primordial bump is located at linear scales, non-linearities tend to produce a second bump at smaller scales. This effect is understood within the halo model due to a more efficient halo formation. In redshift space these nonlinear signatures are partially erased because of the damping along the line-of-sight direction produced by non-coherent motions of particles at small scales. In configuration space, the bump modulates the correlation function reflecting as oscillations in the response, as it is clear in linear Eulerian theory; however, they become damped because large scale coherent flows have some tendency to occupy regions more depleted of particles. This mechanism is explained within Lagrangian Perturbation Theory and well captured by our simulations.


2020 ◽  
Vol 101 (6) ◽  
Author(s):  
Julian B. Muñoz ◽  
Cora Dvorkin ◽  
Francis-Yan Cyr-Racine

Sign in / Sign up

Export Citation Format

Share Document