Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSelf-consistency of the Excursion Set Approach
157
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The excursion set approach provides a framework for predicting how the abundance of dark matter halos depends on the initial conditions. A key ingredient of this formalism comes from the physics of halo formation: the specification of a critical overdensity threshold (barrier) which protohalos must exceed if they are to form bound virialized halos at a later time. Another ingredient is statistical, as it requires the specification of the appropriate statistical ensemble over which to average when making predictions. The excursion set approach explicitly averages over all initial positions, thus implicitly assuming that the appropriate ensemble is that associated with randomly chosen positions in space, rather than special positions such as peaks of the initial density field. Since halos are known to collapse around special positions, it is not clear that the physical and statistical assumptions which underlie the excursion set approach are self-consistent. We argue that they are at least for low mass halos, and illustrate by comparing our excursion set predictions with numerical data from the DEUS simulations.
rate research
Read More
Recently, we provided a simple but accurate formula which closely approximates the first crossing distribution associated with random walks having correlated steps. The approximation is accurate for the wide range of barrier shapes of current interest and is based on the requirement that, in addition to having the right height, the walk must cross the barrier going upwards. Therefore, it only requires knowledge of the bivariate distribution of the walk height and slope, and is particularly useful for excursion set models of the massive end of the halo mass function. However, it diverges at lower masses. We show how to cure this divergence by using a formulation which requires knowledge of just one other variable. While our analysis is general, we use examples based on Gaussian initial conditions to illustrate our results. Our formulation, which is simple and fast, yields excellent agreement with the considerably more computationally expensive Monte-Carlo solution of the first crossing distribution, for a wide variety of moving barriers, even at very low masses.
Insight into a number of interesting questions in cosmology can be obtained from the first crossing distributions of physically motivated barriers by random walks with correlated steps. We write the first crossing distribution as a formal series, ordered by the number of times a walk upcrosses the barrier. Since the fraction of walks with many upcrossings is negligible if the walk has not taken many steps, the leading order term in this series is the most relevant for understanding the massive objects of most interest in cosmology. This first term only requires knowledge of the bivariate distribution of the walk height and slope, and provides an excellent approximation to the first crossing distribution for all barriers and smoothing filters of current interest. We show that this simplicity survives when extending the approach to the case of non-Gaussian random fields. For non-Gaussian fields which are obtained by deterministic transformations of a Gaussian, the first crossing distribution is simply related to that for Gaussian walks crossing a suitably rescaled barrier. Our analysis shows that this is a useful way to think of the generic case as well. Although our study is motivated by the possibility that the primordial fluctuation field was non-Gaussian, our results are general. In particular, they do not assume the non-Gaussianity is small, so they may be viewed as the solution to an excursion set analysis of the late-time, nonlinear fluctuation field rather than the initial one. They are also useful for models in which the barrier height is determined by quantities other than the initial density, since most other physically motivated variables (such as the shear) are usually stochastic and non-Gaussian. We use the Lognormal transformation to illustrate some of our arguments.
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.
We provide a simple formula that accurately approximates the first crossing distribution of barriers having a wide variety of shapes, by random walks with a wide range of correlations between steps. Special cases of it are useful for estimating halo abundances, evolution, and bias, as well as the nonlinear counts in cells distribution. We discuss how it can be extended to allow for the dependence of the barrier on quantities other than overdensity, to construct an excursion set model for peaks, and to show why assembly and scale dependent bias are generic even at the linear level.
The excursion set approach uses the statistics of the density field smoothed on a wide range of scales, to gain insight into a number of interesting processes in nonlinear structure formation, such as cluster assembly, merging and clustering. The approach treats the curve defined by the height of the overdensity fluctuation field when changing the smoothing scale as a random walk. The steps of the walks are often assumed to be uncorrelated, so that the walk heights are a Markov process, even though this assumption is known to be inaccurate for physically relevant filters. We develop a model in which the walk steps, rather than heights, are a Markov process, and correlations between steps arise because of nearest neighbour interactions. This model is a particular case of a general class, which we call Markov Velocity models. We show how these can approximate the walks generated by arbitrary power spectra and filters, and, unlike walks with Markov heights, provide a very good approximation to physically relevant models. We define a Markov Velocity Monte Carlo algorithm to generate walks whose first crossing distribution is very similar to that of TopHat-smoothed LCDM walks. Finally, we demonstrate that Markov Velocity walks generically exhibit a simple but realistic form of assembly bias, so we expect them to be useful in the construction of more realistic merger history trees.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
