No Arabic abstract
Halo merger trees describe the hierarchical mass assembly of dark matter haloes, and are the backbone for modeling galaxy formation and evolution. Merger trees constructed using Monte Carlo algorithms based on the extended Press-Schechter (EPS) formalism are complementary to those extracted from N-body simulations, and have the advantage that they are not trammeled by limited numerical resolution and uncertainties in identifying (sub)haloes and linking them between snapshots. This paper compares multiple EPS-based merger tree algorithms to simulation results using four diagnostics: progenitor mass function (PMF), mass assembly history (MAH), merger rate per descendant halo, and the unevolved subhalo mass function (USMF). In general, algorithms based on spherical collapse yield major-merger rates that are too high by a factor of two, resulting in MAHs that are systematically offset. Assuming ellipsoidal collapse solves most of these issues, but the particular algorithm investigated here that incorporates ellipsoidal collapse dramatically overpredicts the minor-merger rate for massive haloes. The only algorithm in our comparison that yields MAHs, merger rates, and USMFs in good agreement with simulations, is that by Parkinson et al. (2008). However, this is not a true EPS-based algorithm as it draws its progenitor masses from a PMF calibrated against simulations, rather than `predicted by EPS. Finally we emphasize that the benchmarks used to test the EPS algorithms are obtained from simulations and are hampered by significant uncertainties themselves. In particular, MAHs and halo merger rates obtained from simulations by different authors reveal discrepancies that easily exceed 50 percent, even when based on the same simulation. Given this status quo, merger trees constructed using the Parkinson et al. algorithm are as accurate as those extracted from N-body simulations.
We present a new Monte-Carlo algorithm to generate merger trees describing the formation history of dark matter halos. The algorithm is a modification of the algorithm of Cole et al (2000) used in the GALFORM semi-analytic galaxy formation model. As such, it is based on the Extended Press-Schechter theory and so should be applicable to hierarchical models with a wide range of power spectra and cosmological models. It is tuned to be in accurate agreement with the conditional mass functions found in the analysis of merger trees extracted from the LCDM Millennium N-body simulation. We present a comparison of its predictions not only with these conditional mass functions, but also with additional statistics of the Millennium Simulation halo merger histories. In all cases we find it to be in good agreement with the Millennium Simulation and thus it should prove to be a very useful tool for semi-analytic models of galaxy formation and for modelling hierarchical structure formation in general. We have made our merger tree generation code and code to navigate the trees available at http://star-www.dur.ac.uk/~cole/merger_trees .
We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos.
The development of methods and algorithms to solve the $N$-body problem for classical, collisionless, non-relativistic particles has made it possible to follow the growth and evolution of cosmic dark matter structures over most of the Universes history. In the best studied case $-$ the cold dark matter or CDM model $-$ the dark matter is assumed to consist of elementary particles that had negligible thermal velocities at early times. Progress over the past three decades has led to a nearly complete description of the assembly, structure and spatial distribution of dark matter haloes, and their substructure in this model, over almost the entire mass range of astronomical objects. On scales of galaxies and above, predictions from this standard CDM model have been shown to provide a remarkably good match to a wide variety of astronomical data over a large range of epochs, from the temperature structure of the cosmic background radiation to the large-scale distribution of galaxies. The frontier in this field has shifted to the relatively unexplored subgalactic scales, the domain of the central regions of massive haloes, and that of low-mass haloes and subhaloes, where potentially fundamental questions remain. Answering them may require: (i) the effect of known but uncertain baryonic processes (involving gas and stars), and/or (ii) alternative models with new dark matter physics. Here we present a review of the field, focusing on our current understanding of dark matter structure from $N$-body simulations and on the challenges ahead.
We study the empirical relation between an astronomical objects angular momentum $J$ and mass $M$, $J=beta M^alpha$, the $J-M$ relation, using N-body simulations. In particular, we investigate the time evolution of the $J-M$ relation to study how the initial power spectrum and cosmological model affect this relation, and to test two popular models of its origin - mechanical equilibrium and tidal torque theory. We find that in the $Lambda$CDM model, $alpha$ starts with a value of $sim 1.5$ at high redshift $z$, increases monotonically, and finally reaches $5/3$ near $z=0$, whereas $beta$ evolves linearly with time in the beginning, reaches a maximum and decreases, and stabilizes finally. A three-regime scheme is proposed to understand this newly observed picture. We show that the tidal torque theory accounts for this time evolution behaviour in the linear regime, whereas $alpha=5/3$ comes from the virial equilibrium of haloes. The $J-M$ relation in the linear regime contains the information of the power spectrum and cosmological model. The $J-M$ relations for haloes in different environments and with different merging histories are also investigated to study the effects of a halos non-linear evolution. An updated and more complete understanding of the $J-M$ relation is thus obtained.
Recent high-resolution N-body CDM simulations indicate that nonsingular three-parameter models such as the Einasto profile perform better than the singular two-parameter models, e.g. the Navarro, Frenk and White, in fitting a wide range of dark matter haloes. While many of the basic properties of the Einasto profile have been discussed in previous studies, a number of analytical properties are still not investigated. In particular, a general analytical formula for the surface density, an important quantity that defines the lensing properties of a dark matter halo, is still lacking to date. To this aim, we used a Mellin integral transform formalism to derive a closed expression for the Einasto surface density and related properties in terms of the Fox H and Meijer G functions, which can be written as series expansions. This enables arbitrary-precision calculations of the surface density and the lensing properties of realistic dark matter halo models. Furthermore, we compared the Sersic and Einasto surface mass densities and found differences between them, which implies that the lensing properties for both profiles differ.