Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Halo Mass Function from Excursion Set Theory. II. The Diffusing Barrier
145
0
0.0
(
0
)
Authors
Michele Maggiore
Ask ChatGPT about the research
No Arabic abstract
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]
rate research
Read More
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 classic method for computing the mass function of dark matter halos is provided by excursion set theory, where density perturbations evolve stochastically with the smoothing scale, and the problem of computing the probability of halo formation is mapped into the so-called first-passage time problem in the presence of a barrier. While the full dynamical complexity of halo formation can only be revealed through N-body simulations, excursion set theory provides a simple analytic framework for understanding various aspects of this complex process. In this series of paper we propose improvements of both technical and conceptual aspects of excursion set theory, and we explore up to which point the method can reproduce quantitatively the data from N-body simulations. In paper I of the series we show how to derive excursion set theory from a path integral formulation. This allows us both to derive rigorously the absorbing barrier boundary condition, that in the usual formulation is just postulated, and to deal analytically with the non-markovian nature of the random walk. Such a non-markovian dynamics inevitably enters when either the density is smoothed with filters such as the top-hat filter in coordinate space (which is the only filter associated to a well defined halo mass) or when one considers non-Gaussian fluctuations. In these cases, beside ``markovian terms, we find ``memory terms that reflect the non-markovianity of the evolution with the smoothing scale. We develop a general formalism for evaluating perturbatively these non-markovian corrections, and in this paper we perform explicitly the computation of the halo mass function for gaussian fluctuations, to first order in the non-markovian corrections due to the use of a tophat filter in coordinate space.
We use the Excursion Set formalism to compute the properties of the halo mass distribution for a stochastic barrier model which encapsulates the main features of the ellipsoidal collapse of dark matter halos. Non-markovian corrections due to the sharp filtering of the linear density field in real space are computed with the path-integral technique introduced by Maggiore & Riotto (2010). Here, we provide a detailed derivation of the results presented in Corasaniti & Achitouv (2011) and extend the mass function analysis to higher redshift. We also derive an analytical expression for the linear halo bias. We find the analytically derived mass function to be in remarkable agreement with N-body simulation data from Tinker et al. (2008) with differences smaller than ~5% over the range of mass probed by the simulations. The excursion set solution from Monte Carlo generated random walks shows the same level of agreement, thus confirming the validity of the path-integral approach for the barrier model considered here. Similarly the analysis of the linear halo bias shows deviations no greater than 20%. Overall these results indicate that the Excursion Set formalism in combination with a realistic modeling of the conditions of halo collapse can provide an accurate description of the halo mass distribution.
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.
We calculate the mass distribution of Primordial Black Holes (PBHs) produced during metric preheating. After inflation, the oscillations of the inflaton at the bottom of its potential source a parametric resonant instability for small-scale scalar perturbations, that may collapse into black holes. After reviewing in a pedagogical way different techniques that have been developed in the literature to compute mass distributions of PBHs, we focus on the excursion-set approach. We derive a Volterra integral equation that is free of a singularity sometimes encountered, and apply it to the case of metric preheating. We find that if the energy density at which the instability stops, $rho_Gamma$, is sufficiently smaller than the one at which inflation ends, $rho_mathrm{end}$, namely if $rho_Gamma^{1/4}/rho_mathrm{end}^{1/4}< 10^{-5}(rho_mathrm{end}^{1/4}/10^{16}mathrm{GeV})^{3/2}$, then PBHs dominate the universe content at the end of the oscillatory phase. This confirms the previous analysis of arXiv:1907.04236 . By properly accounting for the cloud-in-cloud mechanism, we find that the mass distribution is more suppressed at low masses than previously thought, and peaks several orders of magnitude above the Hubble mass at the end of inflation. The peak mass ranges from $10$ g to stellar masses, giving rise to different possible cosmological effects that we discuss.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
