Abstract
This paper investigates the hierarchy of baryon physics assembly bias relations obtained from state-of-the-art hydrodynamic simulations with respect to the underlying cosmic web spanned by ...the dark matter field. Using the Bias Assignment Method we find that nonlocal bias plays a central role. We classify the cosmic web based on the invariants of the curvature tensor defined not only by the gravitational potential, but especially by the overdensity, as small-scale clustering becomes important in this context. First, the gas density bias relation can be directly mapped onto the dark matter density field to high precision exploiting the strong correlation between them. In a second step, the neutral hydrogen is mapped based on the dark matter and the gas density fields. Finally, the temperature is mapped based on the previous quantities. This permits us to statistically reconstruct the baryon properties within the same simulated volume finding percent precision in the two-point statistics and compatible results in the three-point statistics, in general within 1
σ
, with respect to the reference simulation (with 5–6 orders of magnitude less computing time). This paves the path to establish the best setup for the construction of mocks probing the intergalactic medium for the generation of such key ingredients in the statistical analysis of large forthcoming missions such as DESI, Euclid, J-PAS, and WEAVE.
ABSTRACT
This work investigates the connection between the cosmic web and the halo distribution through the gravitational potential at the field level. We combine three fields of research, cosmic web ...classification, perturbation theory expansions of the halo bias, and halo (galaxy) mock catalogue making methods. In particular, we use the invariants of the tidal field and the velocity shear tensor as generating functions to reproduce the halo number counts of a reference catalogue from full gravity calculations, populating the dark matter field on a mesh well into the non-linear regime ($3\, h^{-1}\, {\rm Mpc}$ scales). Our results show an unprecedented agreement with the reference power spectrum within 1 per cent up to $k=0.72\, h\, {\rm Mpc}^{-1}$. By analysing the three-point statistics on large scales (configurations of up to $k=0.2\, h\, {\rm Mpc}^{-1}$), we find evidence for non-local bias at the 4.8σ confidence level, being compatible with the reference catalogue. In particular, we find that a detailed description of tidal anisotropic clustering on large scales is crucial to achieve this accuracy at the field level. These findings can be particularly important for the analysis of the next generation of galaxy surveys in mock galaxy production.
BAM: bias assignment method to generate mock catalogues Balaguera-Antolínez, A; Kitaura, Francisco-Shu; Pellejero-Ibáñez, Marcos ...
Monthly Notices of the Royal Astronomical Society Letters,
02/2019, Letnik:
483, Številka:
1
Journal Article
Recenzirano
Odprti dostop
ABSTRACT We present BAM: a novel Bias Assignment Method envisaged to generate mock catalogues. Combining the statistics of dark matter tracers from a high-resolution cosmological N-body simulation ...and the dark matter density field calculated from down-sampled initial conditions using efficient structure formation solvers, we extract the halo-bias relation on a mesh of a $3\, h^{-1}$ Mpc cell side resolution as a function of properties of the dark matter density field (e.g. local density, cosmic web type), automatically including stochastic, deterministic, local and non-local components. We use this information to sample the halo density field, accounting for ignored dependencies through an iterative process. By construction, our approach reaches ${\sim } 1\, {\rm {per\,cent}}$ accuracy in the majority of the k-range up to the Nyquist frequency without systematic deviations for power spectra (about $k \sim 1\, h$ Mpc−1) using either particle mesh or Lagrangian perturbation theory based solvers. When using phase-space mapping to compensate the low resolution of the approximate gravity solvers, our method reproduces the bispectra of the reference within $10\, {\rm {per\,cent}}$ precision studying configurations tracing the quasi-non-linear regime. BAM has the potential to become a standard technique to produce mock halo and galaxy catalogues for future galaxy surveys and cosmological studies being highly accurate, efficient and parameter free.
ABSTRACT
We present cosmic birth (COSMological Initial Conditions from Bayesian Inference Reconstructions with THeoretical models): an algorithm to reconstruct the primordial and evolved cosmic ...density fields from galaxy surveys on the light-cone. The displacement and peculiar velocity fields are obtained from forward modelling at different redshift snapshots given some initial cosmic density field within a Gibbs-sampling scheme. This allows us to map galaxies, observed in a light-cone, to a single high redshift and hereby provide tracers and the corresponding survey completeness in Lagrangian space including tetrahedral tessellation mapping. These Lagrangian tracers in turn permit us to efficiently obtain the primordial density field, making the cosmic birth code general to any structure formation model. Our tests are restricted for the time being to augmented Lagrangian perturbation theory. We show how to robustly compute the non-linear Lagrangian bias from clustering measurements in a numerical way, enabling us to get unbiased dark matter field reconstructions at initial cosmic times. We also show that we can accurately recover the information of the dark matter field from the galaxy distribution based on a detailed simulation. Novel key ingredients to this approach are a higher order Hamiltonian-sampling technique and a non-diagonal Hamiltonian mass matrix. This technique could be used to study the Eulerian galaxy bias from galaxy surveys and could become an ideal baryon acoustic reconstruction technique. In summary, this method represents a general reconstruction technique, including in a self-consistent way a survey mask, non-linear and non-local bias, and redshift-space distortions, with an efficiency about 10 times superior to previous comparable methods.
Abstract
In this paper we demonstrate that the information encoded in one single (sufficiently large) N-body simulation can be used to reproduce arbitrary numbers of halo catalogues, using ...approximated realisations of dark matter density fields with different initial conditions. To this end we use as a reference one realisation (from an ensemble of 300) of the Minerva N-body simulations and the recently published Bias Assignment Method to extract the local and non-local bias linking the halo to the dark matter distribution. We use an approximate (and fast) gravity solver to generate 300 dark matter density fields from the down-sampled initial conditions of the reference simulation and sample each of these fields using the halo-bias and a kernel, both calibrated from the arbitrarily chosen realisation of the reference simulation. We show that the power spectrum, its variance and the three-point statistics are reproduced within $\sim 2\%$ (up to k ∼ 1.0 h Mpc−1), $\sim 5-10\%$ and $\sim 10\%$, respectively. Using a model for the real space power spectrum (with three free bias parameters), we show that the covariance matrices obtained from our procedure lead to parameter uncertainties that are compatible within $\sim 10\%$ with respect to those derived from the reference covariance matrix, and motivate approaches that can help to reduce these differences to $\sim 1\%$. Our method has the potential to learn from one simulation with moderate volumes and high-mass resolution and extrapolate the information of the bias and the kernel to larger volumes, making it ideal for the construction of mock catalogues for present and forthcoming observational campaigns such as Euclid or DESI.
Context.
Devising fast and accurate methods of predicting the Lyman-
α
forest at the field level, avoiding the computational burden of running large-volume cosmological hydrodynamic simulations, is ...of fundamental importance to quickly generate the massive set of simulations needed by the state-of-the-art galaxy and Ly
α
forest spectroscopic surveys.
Aims.
We present an improved analytical model to predict the Ly
α
forest at the field level in redshift space from the dark matter field, expanding upon the widely used Fluctuating Gunn-Peterson Approximation (FGPA). Instead of assuming a unique universal relation over the whole considered cosmic volume, we introduce a dependence on the cosmic web environment (knots, filaments, sheets, and voids) in the model, thereby effectively accounting for nonlocal bias. Furthermore, we include a detailed treatment of velocity bias in the redshift space distortion modeling, allowing the velocity bias to be cosmic-web-dependent.
Methods.
We first mapped the dark matter field from real to redshift space through a particle-based relation including velocity bias, depending on the cosmic web classification of the dark matter field in real space. We then formalized an appropriate functional form for our model, building upon the traditional FGPA and including a cutoff and a boosting factor mimicking a threshold and inverse-threshold bias effect, respectively, with model parameters depending on the cosmic web classification in redshift space. Eventually, we fit the coefficients of the model via an efficient Markov chain Monte Carlo scheme.
Results.
We find evidence for a significant difference between the same model parameters in different environments, suggesting that for the investigated setup the simple standard FGPA is not able to adequately predict the Ly
α
forest in the different cosmic web regimes. We reproduce the summary statistics of the reference cosmological hydrodynamic simulation that we use for comparison, yielding an accurate mean transmitted flux, probability distribution function, 3D power spectrum, and bispectrum. In particular, we achieve maximum deviation and average deviation accuracy in the Ly
α
forest 3D power spectrum of ∼3% and ∼0.1% up to
k
∼ 0.4
h
Mpc
−1
, and ∼5% and ∼1.8% up to
k
∼ 1.4
h
Mpc
−1
.
Conclusions.
Our new model outperforms previous analytical efforts to predict the Ly
α
forest at the field level in all the probed summary statistics, and has the potential to become instrumental in the generation of fast accurate mocks for covariance matrices estimation in the context of current and forthcoming Ly
α
forest surveys.
We explore the effects of a positive cosmological constant on astrophysical and cosmological configurations described by a polytropic equation of state. We derive the conditions for equilibrium and ...stability of such configurations and consider some astrophysical examples where our analysis may be relevant. We show that in the presence of the cosmological constant the isothermal sphere is not a viable astrophysical model since the density in this model does not go asymptotically to zero. The cosmological constant implies that, for a polytropic index smaller than 5, the central density has to exceed a certain minimal value in terms of the vacuum density in order to guarantee the existence of a finite-size object. We examine such configurations together with effects of Λ in other exotic possibilities, such as neutrino and boson stars, and we compare our results to N-body simulations. The astrophysical properties and configurations found in this paper are specific features resulting from the existence of a dark energy component. Hence, if found in nature they would be an independent probe of a cosmological constant, complementary to other observations.
The cosmic web from perturbation theory Kitaura, F.-S.; Sinigaglia, F.; Balaguera-Antolínez, A. ...
Astronomy & astrophysics,
03/2024, Letnik:
683
Journal Article
Recenzirano
Odprti dostop
Context. Analysing the large-scale structure (LSS) in the Universe with galaxy surveys demands accurate structure formation models. Such models should ideally be fast and have a clear theoretical ...framework in order to rapidly scan a variety of cosmological parameter spaces without requiring large training data sets. Aims. This study aims to extend Lagrangian perturbation theory (LPT), including viscosity and vorticity, to reproduce the cosmic evolution from dark matter N -body calculations at the field level. Methods. We extend LPT to a Eulerian framework, which we dub eALPT. An ultraviolet regularisation through the spherical collapse model provided by Augmented LPT turns out to be crucial at low redshifts. This iterative method enables modelling of the stress tensor and introduces vorticity. The eALPT model has two free parameters apart from the choice of cosmology, redshift snapshots, cosmic volume, and the number of particles. Results. We find that compared to N -body solvers, the cross-correlation of the dark matter distribution increases at k = 1 h Mpc −1 and z = 0 from ∼55% with the Zel’dovich approximation (∼70% with ALPT), to ∼95% with the three-timestep eALPT, and the power spectra show percentage accuracy up to k ≃ 0.3 h Mpc −1 .