Ellipsoidal Collapse Model Research Articles

Dark matter halo mass functions and density profiles from mass and energy cascade

Halo abundance and structure play a central role for modeling structure formation and evolution. Without relying on a spherical or ellipsoidal collapse model, we analytically derive the halo mass function and cuspy halo density (inner slope of −4/3) based on the mass and energy cascade theory in dark matter flow. The hierarchical halo structure formation leads to halo or particle random walk with a position-dependent waiting time tau _g. First, the inverse mass cascade from small to large scales leads to the halo random walk in mass space with tau _gpropto m_h^{-lambda }, where m_h is the halo mass and lambda is a halo geometry parameter with predicted value of 2/3. The corresponding Fokker-Planck solution for halo random walk in mass space gives rise to the halo mass function with a power-law behavior on small scale and exponential decay on large scale. This can be further improved by considering two different lambda for haloes below and above a critical mass scale m_h^*, i.e. a double-lambda halo mass function. Second, a double-gamma density profile can be derived based on the particle random walk in 3D space with a position-dependent waiting time tau _g propto Phi (r)^{-1} propto r^{-gamma }, where Phi is the gravitational potential and r is the particle distance to halo center. Theory predicts gamma =2/3 that leads to a cuspy density profile with an inner slope of −4/3, consistent with the predicted scaling laws from energy cascade. The Press-Schechter mass function and Einasto density profile are just special cases of proposed models. The small scale permanence can be identified due to the scale-independent rate of mass and energy cascade, where density profiles of different halo masses and redshifts converge to the -4/3 scaling law (rho _h propto r^{-4/3}) on small scales. Theory predicts the halo number density scales with halo mass as propto m_h^{-1.9}, while the halo mass density scales as propto m_h^{4/9}. Results were compared against the Illustris simulations. This new perspective provides a theory for nearly universal halo mass functions and density profiles.

Voids and haloes-in-voids statistics as a probe of the expansion history of the Universe

ABSTRACT Structures in the Universe are arranged into the cosmic web. Distributions, statistics, and evolutions of the structures can be used as probes for cosmological models. We investigate the number density of voids and dark matter haloes-in-voids in the excursion set theory (EST). We study the Markov and non-Markov frameworks of EST in both spherical and ellipsoidal collapse models. Afterward, we compare the number density of voids and haloes-in-voids in the standard lambda cold dark matter and the reconstructed model. The reconstructed model is a model-independent reconstruction based on background observations. This work explores the effects of the collapse model barrier in the different EST frameworks on the statistics of voids and the statistics of haloes-in-voids. Finally, we find the hint that cosmological models can be distinguished by the number density of haloes-in-voids in the 1.0−2.5 redshift range. The maximum difference is observed in z ∼ 1.9.

Primordial black hole merger rate in ellipsoidal-collapse dark matter halo models

We study the merger rate of primordial black holes (PBHs) in the ellipsoidal-collapse model of halo to explain the dark matter abundance by the PBH merger estimated from the gravitational waves detections via the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO) detectors. We indicate that the PBH merger rate within each halo for the ellipsoidal models is more significant than for the spherical models. We specify that the PBH merger rate per unit time and per unit volume for the ellipsoidal-collapse halo models is about one order of magnitude higher than the corresponding spherical models. Moreover, we calculate the evolution of the PBH total merger rate as a function of redshift. The results indicate that the evolution for the ellipsoidal halo models is more sensitive than spherical halo models, as expected from the models. Finally, we present a constraint on the PBH abundance within the context of ellipsoidal and spherical models. By comparing the results with the aLIGO mergers during the third observing run, we show that the merger rate in the ellipsoidal-collapse halo models falls within the aLIGO window, while the same result is not valid for the spherical-collapse ones. Furthermore, we compare the total merger rate of PBHs in terms of their fraction in the ellipsoidal-collapse halo models for several masses of PBHs. The results suggest that the total merger rate of PBHs changes inversely with their masses. We also estimate the relation between the fraction of PBHs and their masses in the ellipsoidal-collapse halo model and show it for a narrow mass distribution of PBHs. The outcome shows that the constraint inferred from the PBH merger rate for the ellipsoidal-collapse halo models can be potentially stronger than the corresponding result obtained for the spherical-collapse ones.

Evaluating backreaction with the ellipsoidal collapse model

We evaluate the effect of structure formation on the average expansion rate with a statistical treatment where density peaks and troughs are modelled as homogeneous ellipsoids. This extends earlier work that used spherical regions. We find that the shear and the presence of filamentary and planar structures have only a small impact on the results. The expansion rate times the age of the universe Ht increases from 2/3 to 0.83 at late times, in order of magnitude agreement with observations, although the change is slower and takes longer than in the real universe. We discuss shortcomings that have to be addressed for this and similar statistical models in the literature to develop into realistic quantitative treatment of backreaction.

Evolution and statistics of non-sphericity of dark matter halos from cosmological N-body simulation

Abstract We revisit the non-sphericity of cluster-mass-scale halos from cosmological N-body simulation on the basis of triaxial modeling. In order to understand the difference between the simulation results and the conventional ellipsoidal collapse model (EC), we first consider the evolution of individual simulated halos. The major difference between EC and the simulation becomes appreciable after the turnaround epoch. Moreover, it is sensitive to the individual evolution history of each halo. Despite such strong dependence on individual halos, the resulting non-sphericity of halos exhibits weak but robust mass dependence in a statistical fashion; massive halos are more spherical up to the turnaround, but gradually become less spherical by z = 0. This is clearly inconsistent with the EC prediction: massive halos are usually more spherical. In addition, at z = 0, inner regions of the simulated halos are less spherical than outer regions; that is, the density distribution inside the halos is highly inhomogeneous and therefore not self-similar (concentric ellipsoids with the same axis ratio and orientation). This is also inconsistent with the homogeneous density distribution that is commonly assumed in EC. Since most of previous fitting formulae for the probability distribution function (PDF) of the axis ratio of triaxial ellipsoids have been constructed under the self-similarity assumption, they are not accurate. Indeed, we compute the PDF of the projected axis ratio a1/a2 directly from the simulation data without the self-similarity assumption, and find that it is very sensitive to the assumption. The latter needs to be carefully taken into account in direct comparison with observations, and therefore we provide an empirical fitting formula for the PDF of a1/a2. Our preliminary analysis suggests that the derived PDF of a1/a2 roughly agrees with the current weak-lensing observations. More importantly, the present results will be useful for future exploration of the non-sphericity of clusters in X-ray and optical observations.

Velocity and mass bias in the distribution of dark matter haloes

The non-linear, scale-dependent bias in the mass distribution of galaxies and the underlying dark matter is a key systematic affecting the extraction of cosmological parameters from galaxy clustering. Using 95 million halos from the Millennium-XXL N-body simulation, we find that the mass bias is scale independent only for $k<0.1 h{\rm Mpc}^{-1}$ today ($z=0$) and for $k<0.2 h{\rm Mpc}^{-1}$ at $z=0.7$. We test analytic halo bias models against our simulation measurements and find that the model of Tinker et al. 2005 is accurate to better then 5% at $z=0$. However, the simulation results are better fit by an ellipsoidal collapse model at $z=0.7$. We highlight, for the first time, another potentially serious systematic due to a sampling bias in the halo velocity divergence power spectra which will affect the comparison between observations and any redshift space distortion model which assumes dark matter velocity statistics with no velocity bias. By measuring the velocity divergence power spectra for different sized halo samples, we find that there is a significant bias which increases with decreasing number density. This bias is approximately 20% at $k=0.1h$Mpc$^{-1}$ for a halo sample of number density $\bar{n} = 10^{-3} (h/$Mpc$)^3$ at both $z=0$ and $z=0.7$ for the velocity divergence auto power spectrum. Given the importance of redshift space distortions as a probe of dark energy and the on-going major effort to advance models for the clustering signal in redshift space, our results show this velocity bias introduces another systematic, alongside scale-dependent halo mass bias, which cannot be neglected.

The formation of CDM haloes – I. Collapse thresholds and the ellipsoidal collapse model

In the excursion set approach to structure formation initially spherical regions of the linear density field collapse to form haloes of mass $M$ at redshift $z_{\rm id}$ if their linearly extrapolated density contrast, averaged on that scale, exceeds some critical threshold, $\delta_{\rm c}(z_{\rm id})$. The value of $\delta_{\rm c}(z_{\rm id})$ is often calculated from the spherical or ellipsoidal collapse model, which provide well-defined predictions given auxiliary properties of the tidal field at a given location. We use two cosmological simulations of structure growth in a $\Lambda$ cold dark matter scenario to quantify $\delta_{\rm c}(z_{\rm id})$, its dependence on the surrounding tidal field, as well as on the shapes of the Lagrangian regions that collapse to form haloes at $z_{\rm id}$. Our results indicate that the ellipsoidal collapse model provides an accurate description of the mean dependence of $\delta_{\rm c}(z_{\rm id})$ on both the strength of the tidal field and on halo mass. However, for a given $z_{\rm id}$, $\delta_{\rm c}(z_{\rm id})$ depends strongly on the halo's characteristic formation redshift: the earlier a halo forms, the higher its initial density contrast. Surprisingly, the majority of haloes forming $today$ fall below the ellipsoidal collapse barrier, contradicting the model predictions. We trace the origin of this effect to the non-spherical shapes of Lagrangian haloes, which arise naturally due to the asymmetry of the linear tidal field. We show that a modified collapse model, that accounts for the triaxial shape of protohaloes, provides a more accurate description of the measured minimum overdensities of recently collapsed objects.

The formation of cold dark matter haloes – II. Collapse time and tides

We use two cosmological simulations of structure formation in the LambdaCDM scenario to study the evolutionary histories of dark-matter haloes and to characterize the Lagrangian regions from which they form. We focus on haloes identified at redshift z_id=0 and show that the classic ellipsoidal collapse model systematically overestimates their collapse times. If one imposes that halo collapse takes place at z_id, this model requires starting from a significantly lower linear density contrast than what is measured in the simulations at the locations of halo formation. We attempt to explain this discrepancy by testing two key assumptions of the model. First, we show that the tides felt by collapsing haloes due to the surrounding large-scale structure evolve non-linearly. Although this effect becomes increasingly important for low-mass haloes, accounting for it in the ellipsoidal collapse model only marginally improves the agreement with N-body simulations. Second, we track the time evolution of the physical volume occupied by forming haloes and show that, after turnaround, it generally stabilizes at a well-defined redshift, z_c>z_id, contrary to the basic assumption of extended Press-Schechter theory based on excursion sets. We discuss the implications of this result for understanding the origin of the mass-dependence and scatter in the linear threshold for halo formation. Finally, we show that, when tuned for collapse at z_c, a modified version of the ellipsoidal collapse model that also accounts for the triaxial nature of protohaloes predicts their linear density contrast in an unbiased way.

Ellipticity and prolaticity of the initial gravitational-shear field at the position of density maxima

Dark-matter haloes are supposed to form at the positions of maxima in the initial matter density field. The gravitational-shear field's ellipticity and prolaticity that serve as input for the ellipsoidal-collapse model, however, are derived from a distribution that does not take the additional maximum constraint into account. In this article, I quantify the variations of the most probable and the expected values of the ellipticity and the prolaticity when considering this additional constraint as well as the implications for the ellipsoidal-collapse model. Based on the statistics of Gaussian random fields, it is possible to set up a joint distribution for the eigenvalues of the gravitational-shear tensor and the matter density that incorporates the maximum constraint by invoking a vanishing first derivative and a negative definite second derivative of the density field into the calculation. In the density range relevant for cosmological structure formation, both the most probable and the expected value of the ellipticity calculated from the standard distribution used in the literature are about 3-8 per cent higher compared to the ones calculated under the additional assumption of a density maximum. Additionally, the analogous quantities for the prolaticity do not vanish but acquire slightly positive values in the range of $10^{-3}$-$10^{-2}$. For large overdensities, the predictions from both distributions converge. The values for $\delta_\mathrm{c}$ and $\Delta_\mathrm{v}$ derived from the ellipsoidal-collapse model using the standard distribution for the initial ellipticity and prolaticity are up to 4 and 6 per cent higher, respectively, than those obtained taking the additional maximum constraint into account in the range of $10^{13}$-$10^{15}~h^{-1}~\mathrm{M}_\odot$ in mass and 0-2 in redshift.

Halo shapes, initial shear field, and cosmic web

The ellipsoidal collapse model, combined with the excursion set theory, allows one to estimate the shapes of dark matter halos as seen in high-resolution numerical simulations. The same theoretical framework predicts a quasi-universal behaviour for the conditional axis ratio distributions at later times, set by initial conditions and unaltered by non-linear evolution. The formalism for halo shapes is also useful in making the connection with the initial shear field of the cosmic web, which plays a crucial role in the formation of large-scale structures. The author has briefly discussed the basic aspects of the modelling, as well as the implications of a new formula for the constrained eigenvalues of the initial shear field, given the fact that positions are peaks or dips in the corresponding density field – and not random locations. This formula leads to a new generalized excursion set algorithm for peaks in Gaussian random fields. The results highlighted, here, are relevant for a number of applications, especially for weak lensing studies and for devising algorithms to find and classify structures in the cosmic web.

Path integral approach for merger rates of dark matter haloes

We use path integrals in order to estimate merger rates of dark matter haloes using the Extended Press-Schechter approximation (EPS) for the Spherical Collapse (SC) and the Ellipsoidal Collapse (EC) models. Merger rates have been calculated for masses in the range $10^{10}M_{\odot}\mathrm{h}^{-1}$ to $10^{14}M_{\odot}\mathrm{h}^{-1}$ and for redshifts $z$ in the range 0 to 3. A detailed comparison between these models is presented. Path approach gives a better agreement with the exact solutions for constrained distributions than the approach of \cite{shto02}. Although this improvement seems not to be very large, our results show that the path approach is a step to the right direction. Differences between the two widely used barriers, spherical and ellipsoidal, depend crucially on the mass of the descendant halo. These differences become larger for decreasing mass of the descendant halo. The use of additional terms in the expansion used in the path approach, other improvements as well as detailed comparisons with the predictions of N-body simulations, that could improve our understanding about the important issue of structure formation, are under study.

Modelling the shapes of the largest gravitationally bound objects

We combine the physics of the ellipsoidal collapse model with the excursion set theory to study the shapes of dark matter halos. In particular, we develop an analytic approximation to the nonlinear evolution that is more accurate than the Zeldovich approximation; we introduce a planar representation of halo axis ratios, which allows a concise and intuitive description of the dynamics of collapsing regions and allows one to relate the final shape of a halo to its initial shape; we provide simple physical explanations for some empirical fitting formulae obtained from numerical studies. Comparison with simulations is challenging, as there is no agreement about how to define a non-spherical gravitationally bound object. Nevertheless, we find that our model matches the conditional minor-to-intermediate axis ratio distribution rather well, although it disagrees with the numerical results in reproducing the minor-to-major axis ratio distribution. In particular, the mass dependence of the minor-to-major axis distribution appears to be the opposite to what is found in many previous numerical studies, where low-mass halos are preferentially more spherical than high-mass halos. In our model, the high-mass halos are predicted to be more spherical, consistent with results based on a more recent and elaborate halo finding algorithm, and with observations of the mass dependence of the shapes of early-type galaxies. We suggest that some of the disagreement with some previous numerical studies may be alleviated if we consider only isolated halos.

Toward a Universal Formulation of the Halo Mass Function

We compute the dark matter halo mass function using the excursion set formalism for a diffusive barrier with linearly drifting average which captures the main features of the ellipsoidal collapse model. We evaluate the non-Markovian corrections due to the sharp filtering of the linear density field in real space with a path-integral method. We find an unprecedented agreement with N-body simulation data with deviations ≲5% over the range of masses probed by the simulations. This indicates that the excursion set in combination with a realistic modeling of the collapse threshold can provide a robust estimation of the halo mass function.

Merger rates of dark matter haloes: a comparison between EPS and N-body results

We calculate merger rates of dark matter haloes using the Extended Press-Schechter approximation (EPS) for the Spherical Collapse (SC) and the Ellipsoidal Collapse (EC) models. Merger rates have been calculated for masses in the range $10^{10}M_{\odot}\mathrm{h}^{-1}$ to $10^{14}M_{\odot}\mathrm{h}^{-1}$ and for redshifts $z$ in the range 0 to 3 and they have been compared with merger rates that have been proposed by other authors as fits to the results of N-body simulations. The detailed comparison presented here shows that the agreement between the analytical models and N-body simulations depends crucially on the mass of the descendant halo. For some range of masses and redshifts either SC or EC models approximate satisfactory the results of N-body simulations but for other cases both models are less satisfactory or even bad approximations. We showed, by studying the parameters of the problem that a disagreement --if it appears-- does not depend on the values of the parameters but on the kind of the particular solution used for the distribution of progenitors or on the nature of EPS methods. Further studies could help to improve our understanding about the physical processes during the formation of dark matter haloes.

STATISTICAL PROPERTIES OF SUPERCLUSTER-LIKE FILAMENTS FROM COSMOLOGICAL SIMULATIONS

We study large-scale structures from numerical simulations, paying particular attention to supercluster-like structures. A grid-density-contour based algorithm is adopted to locate connected groups. With the increase of the linking density threshold, the foam- like cosmic web is subsequently broken into individual supercluster-like groups and further halos which are in accordance to groups with the linking density threshold {\rho}/{\rho}= 1 + {\delta} = 80. By analyzing sets of cosmological simulations with varying cosmological parameters, we find that an universal mass function exists not only for halos but also for low-density supercluster-like groups until the linking density threshold decreases to a density where the global percolation of large-scale structures occurs. We further show that the mass functions of different groups can be well described by the Jenkins form with the parameters being dependent on the linking density threshold. On the other hand, these low- density supercluster-like groups cannot be directly associated with the predictions from the excursion set theory with effective barriers obtained from dynamical collapse models, and the peak exclusion effect must be taken into account. Including such an effect, the mass function of groups with the linking density threshold 1 + {\delta} = 16 is in good agreements with that from the excursion set theory with a nearly flat effective barrier. A simplified analysis of the ellipsoidal collapse model indicates that the barrier for collapses along two axes to form filaments is approximately flat in scales, thus we define groups identified with 1 + {\delta} = 16 as filaments. We then further study the halo-filament conditional mass function and the filament-halo conditional mass function, and compare them with the predictions from the two-barrier excursion set theory. The shape statistics for filaments are also presented.

THE LARGE-SCALE BIAS OF DARK MATTER HALOS: NUMERICAL CALIBRATION AND MODEL TESTS

We measure the clustering of dark matter halos in a large set of collisionless cosmological simulations of the flat LCDM cosmology. Halos are identified using the spherical overdensity algorithm, which finds the mass around isolated peaks in the density field such that the mean density is Delta times the background. We calibrate fitting functions for the large scale bias that are adaptable to any value of Delta we examine. We find a ~6% scatter about our best fit bias relation. Our fitting functions couple to the halo mass functions of Tinker et. al. (2008) such that bias of all dark matter is normalized to unity. We demonstrate that the bias of massive, rare halos is higher than that predicted in the modified ellipsoidal collapse model of Sheth, Mo, & Tormen (2001), and approaches the predictions of the spherical collapse model for the rarest halos. Halo bias results based on friends-of-friends halos identified with linking length 0.2 are systematically lower than for halos with the canonical Delta=200 overdensity by ~10%. In contrast to our previous results on the mass function, we find that the universal bias function evolves very weakly with redshift, if at all. We use our numerical results, both for the mass function and the bias relation, to test the peak-background split model for halo bias. We find that the peak-background split achieves a reasonable agreement with the numerical results, but ~20% residuals remain, both at high and low masses.

Triaxial collapse and virialisation of dark-matter haloes

We reconsider the ellipsoidal-collapse model and extend it in two ways: We modify the treatment of the external gravitational shear field, introducing a hybrid model in between linear and non-linear evolution, and we introduce a virialisation criterion derived from the tensor virial theorem to replace the ad-hoc criterion employed so far. We compute the collapse parameters delta_c and Delta_v and find that they increase with ellipticity e and decrease with prolaticity p. We marginalise them over the appropriate distribution of e and p and show the marginalised results as functions of halo mass and virialisation redshift. While the hybrid model for the external shear gives results very similar to those obtained from the non-linear model, ellipsoidal collapse changes the collapse parameters typically by (20...50)%, in a way increasing with decreasing halo mass and decreasing virialisation redshift. We qualitatively confirm the dependence on mass and virialisation redshift of a fitting formula for delta_c, but find noticeable quantitative differences in particular at low mass and high redshift. The derived mass function is in good agreement with mass functions recently proposed in the literature.

THE HALO MASS FUNCTION FROM EXCURSION SET THEORY. II. THE DIFFUSING BARRIER

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/2\sigma^2\}$ to $\exp\{-a\delta_c^2/2\sigma^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]

COLLAPSE BARRIERS AND HALO ABUNDANCE: TESTING THE EXCURSION SET ANSATZ

Our heuristic understanding of the abundance of dark matter halos centers around the concept of a density threshold, or "barrier", for gravitational collapse. If one adopts the ansatz that regions of the linearly evolved density field smoothed on mass scale M with an overdensity that exceeds the barrier will undergo gravitational collapse into halos of mass M, the corresponding abundance of such halos can be estimated simply as a fraction of the mass density satisfying the collapse criterion divided by the mass M. The key ingredient of this ansatz is therefore the functional form of the collapse barrier as a function of mass M or, equivalently, of the variance sigma^2(M). Several such barriers based on the spherical, Zel'dovich, and ellipsoidal collapse models have been extensively discussed. Using large scale cosmological simulations, we show that the relation between the linear overdensity and the mass variance for regions that collapse to form halos by the present epoch resembles expectations from dynamical models of ellipsoidal collapse. However, we also show that using such a collapse barrier with the excursion set ansatz predicts a halo mass function inconsistent with that measured directly in cosmological simulations. This inconsistency demonstrates a failure of the excursion set ansatz as a physical model for halo collapse. We discuss implications of our results for understanding the collapse epoch for halos as a function of mass, and avenues for improving consistency between analytical models for the collapse epoch and the results of cosmological simulations.

Merger history trees of dark matter haloes in moving barrier models

We present an algorithm for generating merger histories of dark matter haloes. The algorithm is based on the excursion set approach with moving barriers whose shape is motivated by the ellipsoidal collapse model of halo formation. In contrast to most other merger-tree algorithms, ours takes discrete steps in mass rather than time. This allows us to quantify effects which arise from the fact that outputs from numerical simulations are usually in discrete time bins. In addition, it suggests a natural set of scaling variables for describing the abundance of halo progenitors; this scaling is not as general as that associated with a spherical collapse. We test our algorithm by comparing its predictions with measurements in numerical simulations. The progenitor mass fractions and mass functions are in good agreement, as is the predicted scaling law. We also test the formation-redshift distribution, the mass distribution at formation, and the redshift distribution of the most recent major merger; all are in reasonable agreement with N-body simulation data, over a broad range of masses and redshifts. Finally, we study the effects of sampling in discrete time snapshots. In all cases, the improvement over algorithms based on the spherical collapse assumption is significant.