Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userExcursion Set Theory for Correlated Random Walks
475
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a new method to compute the first crossing distribution in excursion set theory for the case of correlated random walks. We use a combination of the path integral formalism of Maggiore & Riotto, and the integral equation solution of Zhang & Hui, and Benson et al. to find a numerically robust and convenient algorithm to derive the first crossing distribution in terms of a perturbative expansion around the limit of an uncorrelated random walk. We apply this methodology to the specific case of a Gaussian random density field filtered with a Gaussian smoothing function. By comparing our solutions to results from Monte Carlo calculations of the first crossing distribution we demonstrate that our method accurate for power spectra $P(k)propto k^n$ for $n=1$, becoming less accurate for smaller values of $n$. It is therefore complementary to the method of Musso & Sheth, which will therefore be more useful for standard $Lambda$CDM power spectra. Our approach is quite general, and can be adapted to other smoothing functions, and also to non-Gaussian density fields.
rate research
Read More
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.
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.
In excursion set theory the computation of the halo mass function is mapped into a first-passage time process in the presence of a barrier, which in the spherical collapse model is a constant and in the ellipsoidal collapse model is a fixed function of the variance of the smoothed density field. However, N-body simulations show that dark matter halos grow through a mixture of smooth accretion, violent encounters and fragmentations, and modeling halo collapse as spherical, or even as ellipsoidal, is a significant oversimplification. We propose that some of the physical complications inherent to a realistic description of halo formation can be included in the excursion set theory framework, at least at an effective level, by taking into account that the critical value for collapse is not a fixed constant $delta_c$, as in the spherical collapse model, nor a fixed function of the variance $sigma$ of the smoothed density field, as in the ellipsoidal collapse model, but rather is itself a stochastic variable, whose scatter reflects a number of complicated aspects of the underlying dynamics. Solving the first-passage time problem in the presence of a diffusing barrier we find that the exponential factor in the Press-Schechter mass function changes from $exp{-delta_c^2/2sigma^2}$ to $exp{-adelta_c^2/2sigma^2}$, where $a=1/(1+D_B)$ and $D_B$ is the diffusion coefficient of the barrier. The numerical value of $D_B$, and therefore the corresponding value of $a$, depends among other things on the algorithm used for identifying halos. We discuss the physical origin of the stochasticity of the barrier and we compare with the mass function found in N-body simulations, for the same halo definition.[Abridged]
We compute the effect of primordial non-Gaussianity on the halo mass function, using excursion set theory. In the presence of non-Gaussianity the stochastic evolution of the smoothed density field, as a function of the smoothing scale, is non-markovian and beside local terms that generalize Press-Schechter (PS) theory, there are also memory terms, whose effect on the mass function can be computed using the formalism developed in the first paper of this series. We find that, when computing the effect of the three-point correlator on the mass function, a PS-like approach which consists in neglecting the cloud-in-cloud problem and in multiplying the final result by a fudge factor close to 2, is in principle not justified. When computed correctly in the framework of excursion set theory, in fact, the local contribution vanishes (for all odd-point correlators the contribution of the image gaussian cancels the Press-Schechter contribution rather than adding up), and the result comes entirely from non-trivial memory terms which are absent in PS theory. However it turns out that, in the limit of large halo masses, where the effect of non-Gaussianity is more relevant, these memory terms give a contribution which is the the same as that computed naively with PS theory, plus subleading terms depending on derivatives of the three-point correlator. We finally combine these results with the diffusive barrier model developed in the second paper of this series, and we find that the resulting mass function reproduces recent N-body simulations with non-Gaussian initial conditions, without the introduction of any ad hoc parameter.
A sizeable level of non-Gaussianity in the primordial cosmological perturbations may be induced by a large trispectrum, i.e. by a large connected four-point correlation function. We compute the effect of a primordial non-Gaussian trispectrum on the halo mass function, within excursion set theory. We use the formalism that we have developed in a previous series of papers and which allows us to take into account the fact that, in the presence of non-Gaussianity, the stochastic evolution of the smoothed density field, as a function of the smoothing scale, is non-markovian. In the large mass limit, the leading-order term that we find agrees with the leading-order term of the results found in the literature using a more heuristic Press-Schecther (PS)-type approach. Our approach however also allows us to evaluate consistently the subleading terms, which depend not only on the four-point cumulant but also on derivatives of the four-point correlator, and which cannot be obtained within non-Gaussian extensions of PS theory. We perform explicitly the computation up to next-to-leading order.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
