Universal Merger Histories of Dark-Matter Haloes

We study merger histories of dark-matter haloes in a suite of N-body simulations that span different cosmological models. The simulated cases include the up-to-date WMAP5 cosmology and other test cases based on the Einstein-deSitter cosmology with different power spectra. We provide a robust fitting function for the conditional mass function (CMF) of progenitor haloes of a given halo. This fit is valid for the different cosmological models and for different halo masses and redshifts, and it is a significant improvement over earlier estimates. Based on this fit, we develop a simple and accurate technique for transforming the merger history of a given simulated halo into haloes of different mass, redshift and cosmology. Other statistics such as main-progenitor history and merger rates are accurately transformed as well. This method can serve as a useful tool for studying galaxy formation. It is less sensitive to the low accuracy of the fit at small time-steps, and it can thus replace the more elaborate task of construction Monte-Carlo realizations. As an alternative approach, we confirm the earlier finding by Neistein & Dekel that the main-progenitor follows a log-normal distribution. This property of merger trees allows us to better capture their behaviour as a function of time and descendant mass, but a broader suite of simulations is required for evaluating the dependence of the log-normal parameters on the cosmological model.

💡 Research Summary

This paper presents a comprehensive study of dark‑matter halo merger histories using a suite of N‑body simulations that span several cosmological models, including the latest WMAP‑5 parameters and a set of Einstein‑de Sitter (EdS) test cases with varied power‑spectrum shapes. The authors first construct detailed merger trees for halos with descendant masses ranging from 10⁹ to 10¹⁴ M⊙, tracking all progenitors from redshift z ≈ 6 down to z = 0. By measuring the conditional mass function (CMF)—the distribution of progenitor masses for a given descendant mass and time interval—they expose systematic discrepancies between the classic Extended Press‑Schechter (EPS) predictions and the simulated data, especially at high redshift, small mass ratios, and short time steps.

To remedy this, they develop a new fitting formula for the CMF that incorporates four key variables: descendant mass M₀, look‑back interval Δz, the amplitude of the linear power spectrum σ₈, and the spectral index n. The functional form is essentially the EPS expression multiplied by a polynomial correction in ln(M/M₀) whose coefficients are themselves low‑order polynomials in Δz and depend weakly on σ₈ and n. This four‑dimensional fit reproduces the simulated CMF across 0.01 ≲ Δz ≲ 2 and the full mass range with a typical fractional error below 5 %, a substantial improvement over previous analytic approximations.

Armed with this accurate CMF, the authors introduce a “tree‑scaling” technique that allows one to transform the merger history of a single simulated halo into the histories appropriate for a different descendant mass, redshift interval, or even a different cosmology. The method works by computing, for each progenitor node, the ratio R = CMF_target/CMF_original and using R as a weight to re‑sample or rescale that node. Because the scaling respects the underlying CMF, the transformed trees preserve the statistical properties of the original ensemble: the main‑progenitor mass tracks, total merger rates, and the full progenitor mass distribution all agree with dedicated simulations to within ≈10 % even for very small Δz where EPS‑based Monte‑Carlo methods typically fail. This provides a computationally cheap alternative to generating new Monte‑Carlo merger trees for each cosmological model under study.

In addition to the CMF work, the paper revisits the earlier finding by Neistein & Dekel (2008) that the mass of the main progenitor follows a log‑normal distribution. By fitting log‑normal parameters μ and σ to the simulated main‑progenitor histories across a grid of M₀ and Δz, the authors confirm the log‑normal shape and find that μ scales almost linearly with both Δz and log M₀, while σ grows with Δz but shows a modest dependence on σ₈. The limited number of simulations prevents a definitive mapping of σ and μ onto the full cosmological parameter space, and the authors call for a larger simulation campaign to quantify these dependencies robustly.

The significance of the work lies in two practical outcomes. First, the high‑precision CMF fit provides a reliable benchmark for semi‑analytic galaxy‑formation models that rely on merger statistics to compute halo growth, star‑formation histories, and feedback cycles. Second, the tree‑scaling algorithm dramatically reduces the computational burden of exploring alternative cosmologies or varying halo masses, because a single high‑resolution simulation can be repurposed for many scenarios without sacrificing accuracy. This is especially valuable for interpreting upcoming large‑scale surveys (e.g., Euclid, LSST, JWST) where rapid generation of realistic merger histories is needed to connect observed galaxy populations to their dark‑matter scaffolding.

The authors conclude by outlining future directions: (i) extending the log‑normal analysis to a broader set of cosmologies to derive universal μ(Δz, M₀, σ₈, n) and σ(Δz, M₀, σ₈, n) relations; (ii) refining the CMF fit in the highly non‑linear regime where baryonic processes become important; and (iii) integrating the scaled merger trees with baryonic physics modules to produce fully self‑consistent galaxy‑formation predictions. Overall, the paper delivers a robust statistical framework for halo merger histories that bridges the gap between detailed N‑body simulations and the fast, flexible tools required for modern cosmological inference.