Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userModelling the stochasticity of high-redshift halo bias
60
0
0.0
(
0
)
Added by
Ainulnabilah Nasirudin
Publication date
2019
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
A very large dynamic range with simultaneous capture of both large- and small-scales in the simulations of cosmic structures is required for correct modelling of many cosmological phenomena, particularly at high redshift. This is not always available, or when it is, it makes such simulations very expensive. We present a novel sub-grid method for modelling low-mass ($10^5,M_odotleq M_{rm halo}leq 10^9,M_odot$) haloes, which are otherwise unresolved in large-volume cosmological simulations limited in numerical resolution. In addition to the deterministic halo bias that captures the average property, we model its stochasticity that is correlated in time. We find that the instantaneous binned distribution of the number of haloes is well approximated by a log-normal distribution, with overall amplitude modulated by this temporal correlation bias. The robustness of our new scheme is tested against various statistical measures, and we find that temporally correlated stochasticity generates mock halo data that is significantly more reliable than that from temporally uncorrelated stochasticity. Our method can be applied for simulating processes that depend on both the small- and large-scale structures, especially for those that are sensitive to the evolution history of structure formation such as the process of cosmic reionization. As a sample application, we generate a mock distribution of medium-mass ($ 10^{8} leq M/M_{odot} leq 10^{9}$) haloes inside a 500 Mpc$,h^{-1}$, $300^3$ grid simulation box. This mock halo catalogue bears a reasonable statistical agreement with a halo catalogue from numerically-resolved haloes in a smaller box, and therefore will allow a very self-consistent sets of cosmic reionization simulations in a box large enough to generate statistically reliable data.
rate research
Read More
We study the relationship between dark-matter haloes and matter in the MIP $N$-body simulation ensemble, which allows precision measurements of this relationship, even deeply into voids. What enables this is a lack of discreteness, stochasticity, and exclusion, achieved by averaging over hundreds of possible sets of initial small-scale modes, while holding fixed large-scale modes that give the cosmic web. We find (i) that dark-matter-halo formation is greatly suppressed in voids; there is an exponential downturn at low densities in the otherwise power-law matter-to-halo density bias function. Thus, the rarity of haloes in voids is akin to the rarity of the largest clusters, and their abundance is quite sensitive to cosmological parameters. The exponential downturn appears both in an excursion-set model, and in a model in which fluctuations evolve in voids as in an open universe with an effective $Omega_m$ proportional to a large-scale density. We also find that (ii) haloes typically populate the average halo-density field in a super-Poisson way, i.e. with a variance exceeding the mean; and (iii) the rank-order-Gaussianized halo and dark-matter fields are impressively similar in Fourier space. We compare both their power spectra and cross-correlation, supporting the conclusion that one is roughly a strictly-increasing mapping of the other. The MIP ensemble especially reveals how halo abundance varies with `environmental quantities beyond the local matter density; (iv) we find a visual suggestion that at fixed matter density, filaments are more populated by haloes than clusters.
We study the effect of large-scale tidal fields on internal halo properties using a set of N-body simulations. We measure significant cross-correlations between large-scale tidal fields and several non-scalar halo properties: shapes, velocity dispersion, and angular momentum. Selection effects that couple to these non-scalar halo properties can produce anisotropic clustering even in real-space. We investigate the size of this effect and show that it can produce a non-zero quadrupole similar in size to the one generated by linear redshift-space distortions (RSD). Finally, we investigate the clustering properties of halos identified in redshift-space and find enormous deviations from the standard linear RSD model, again caused by anisotropic assembly bias. These effects could contaminate the values of cosmological parameters inferred from the observed redshift-space clustering of galaxies, groups, or 21cm emission from atomic hydrogen, if their selection depends on properties affected by halo assembly bias. We briefly discuss ways in which this effect can be measured in existing and future large-scale structure surveys.
The simplest stochastic halo formation models assume that the traceless part of the shear field acts to increase the initial overdensity (or decrease the underdensity) that a protohalo (or protovoid) must have if it is to form by the present time. Equivalently, it is the difference between the overdensity and (the square root of the) shear that must be larger than a threshold value. To estimate the effect this has on halo abundances using the excursion set approach, we must solve for the first crossing distribution of a barrier of constant height by the random walks associated with the difference, which is now (even for Gaussian initial conditions) a non-Gaussian variate. The correlation properties of such non-Gaussian walks are inherited from those of the density and the shear, and, since they are independent processes, the solution is in fact remarkably simple. We show that this provides an easy way to understand why earlier heuristic arguments about the nature of the solution worked so well. In addition to modelling halos and voids, this potentially simplifies models of the abundance and spatial distribution of filaments and sheets in the cosmic web.
It has been recently shown that any halo velocity bias present in the initial conditions does not decay to unity, in agreement with predictions from peak theory. However, this is at odds with the standard formalism based on the coupled fluids approximation for the coevolution of dark matter and halos. Starting from conservation laws in phase space, we discuss why the fluid momentum conservation equation for the biased tracers needs to be modified in accordance with the change advocated in Baldauf, Desjacques & Seljak (2014). Our findings indicate that a correct description of the halo properties should properly take into account peak constraints when starting from the Vlasov-Boltzmann equation.
Small-scale density fluctuations can significantly affect reionization but are typically modelled quite crudely. Unresolved fluctuations in numerical simulations and analytical calculations are included using a gas clumping factor, typically assumed to be independent of the local environment. In Paper I, we presented an improved, local density-dependent model for the sub-grid gas clumping. Here we extend this using an empirical stochastic model based on the results from high-resolution numerical simulations which fully resolve all relevant fluctuations. Our model reproduces well both the mean density-clumping relation and its scatter. We applied our stochastic model, along with the mean clumping one and the Paper I deterministic model, to create a large-volume realisation of the clumping field, and used these in radiative transfer simulations of cosmic reionization. Our results show that the simplistic mean clumping model delays reionization compared to local density-dependent models, despite producing fewer recombinations overall. This is due to the very different spatial distribution of clumping, resulting in much higher photoionization rates in the latter cases. The mean clumping model produces smaller HII regions throughout most of the reionization, but those percolate faster at late times. It also causes a significant delay in the 21-cm fluctuations peak and yields lower non-Gaussianity and many fewer bright pixels in the PDF distribution. The stochastic density-dependent model shows relatively minor differences from the deterministic one, mostly concentrated around overlap, where it significantly suppresses the 21-cm fluctuations, and at the bright tail of the 21-cm PDFs, where it produces noticeably more bright pixels.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
