Dark Matter Halo Merger Research Articles

Super-Eddington accretion in high-redshift black holes and the emergence of jetted AGN

Abstract In this work we study the co-evolution of central black holes (BHs) and host galaxies by utilizing an advanced iteration of the DELPHI semi-analytical model of galaxy formation and evolution. Based on dark matter halo merger trees spanning the redshift range from z = 20 to z = 4, it now incorporates essential components such as gas heating and cooling, cold and hot BH accretion, jet and radiative AGN feedback. We show how different BH growth models impact quasar and galaxy observables at z ≥ 5, providing predictions that will help discriminate between super-Eddington and Eddington-limited accretion models: despite being both consistent with observed properties of SMBHs and their host galaxies at z ∼ 5 − 7, they become very clearly distinguishable at higher redshift and in the intermediate mass regime. We find that the super-Eddington model, unlike the Eddington-limited scenario, predicts a gap in the BH mass function corresponding to the intermediate-mass range 104 M⊙ < Mbh < 106 M⊙. Additionally, it predicts black holes up to two orders of magnitude more massive for the same stellar mass at z = 9. The resulting velocity dispersion – BH mass relation at z ≥ 5 is consistent with local measurements, suggesting that its slope and normalisation are independent of redshift. Depending on the Eddington ratio, we also model the emergence of AGN jets, predicting their duty cycle across as a function of BH mass and their potential impact on the observed number density distribution of high-redshift AGN in the hard X-ray band.

The effect of pressure-anisotropy-driven kinetic instabilities on magnetic field amplification in galaxy clusters

Context. The intracluster medium (ICM) is the low-density diffuse gas that fills the space between galaxies within galaxy clusters. It is primarily composed of magnetized plasma, which reaches virial temperatures of up to 108 K, probably due to mergers of subhalos. Under these conditions, the plasma is weakly collisional and therefore has an anisotropic pressure tensor with respect to the local direction of the magnetic field. This triggers very fast, Larmor-scale, pressure-anisotropy-driven kinetic instabilities that alter magnetic field amplification. Aims. We aim to study magnetic field amplification through a turbulent, small-scale dynamo, including the effects of the kinetic instabilities, during the evolution of a typical massive galaxy cluster. A specific aim of this work is to establish a redshift limit from which a dynamo has to start to amplify the magnetic field up to equipartition with the turbulent velocity field at redshift z = 0. Methods. We implemented one-dimensional radial profiles for various plasma quantities for merger trees generated with the modified GALFORM algorithm. We assumed that turbulence is driven by successive mergers of dark matter halos and constructed effective models for the Reynolds number Reeff dependence on the magnetic field in three different magnetization regimes (unmagnetized, magnetized “kinetic”, and magnetized “fluid”), including the effects of kinetic instabilities. The magnetic field growth rate is calculated for the different Reeff models. Results. The model results in a higher magnetic field growth rate at higher redshift. For all scenarios considered in this study, to reach equipartition at z = 0, it is sufficient for the amplification of the magnetic field to start at redshift zstart ≈ 1.5 and above. The time to reach equipartition can be significantly shorter in cases with systematically smaller turbulent forcing scales and for the highest Reeff models. Conclusions. The origin of magnetic fields in the weakly collisional ICM can be explained by the small-scale turbulent dynamo, provided that the dynamo process starts beyond a given redshift. Merger trees are useful tools for studying the evolution of magnetic fields in weakly collisional plasmas, and could also be used to constrain the different stages of the dynamo that could potentially be observed by future radio telescopes.

A Global Semianalytic Model of the First Stars and Galaxies Including Dark Matter Halo Merger Histories

We present a new self-consistent semianalytic model of the first stars and galaxies to explore the high-redshift (z ≥ 15) Population III (PopIII) and metal-enriched star formation histories. Our model includes the detailed merger history of dark matter halos generated with Monte Carlo merger trees. We calibrate the minimum halo mass for PopIII star formation from recent hydrodynamical cosmological simulations that simultaneously include the baryon–dark matter streaming velocity, Lyman–Werner (LW) feedback, and molecular hydrogen self-shielding. We find an overall increase in the resulting star formation rate density (SFRD) compared to calibrations based on previous simulations (e.g., the PopIII SFRD is over an order of magnitude higher at z = 35−15). We evaluate the effect of the halo-to-halo scatter in this critical mass and find that it increases the PopIII stellar mass density by a factor ∼1.5 at z ≥ 15. Additionally, we assess the impact of various semianalytic/analytic prescriptions for halo assembly and star formation previously adopted in the literature. For example, we find that models assuming smooth halo growth computed via abundance matching predict SFRDs similar to the merger tree model for our fiducial model parameters, but that they may underestimate the PopIII SFRD in cases of strong LW feedback. Finally, we simulate subvolumes of the Universe with our model both to quantify the reduction in total star formation in numerical simulations due to a lack of density fluctuations on spatial scales larger than the simulation box, and to determine spatial fluctuations in SFRD due to the diversity in halo abundances and merger histories.

ACACIA: a new method to produce on-the-fly merger trees in the ramses code

ABSTRACT The implementation of ACACIA, a new algorithm to generate dark matter halo merger trees with the Adaptive Mesh Refinement code RAMSES, is presented. The algorithm is fully parallel and based on the Message Passing Interface. As opposed to most available merger tree tools, it works on the fly during the course of the N-body simulation. It can track dark matter substructures individually using the index of the most bound particle in the clump. Once a halo (or a sub-halo) merges into another one, the algorithm still tracks it through the last identified most bound particle in the clump, allowing to check at later snapshots whether the merging event was definitive, or whether it was only temporary, with the clump only traversing another one. The same technique can be used to track orphan galaxies that are not assigned to a parent clump anymore because the clump dissolved due to numerical overmerging. We study in detail the impact of various parameters on the resulting halo catalogues and corresponding merger histories. We then compare the performance of our method using standard validation diagnostics, demonstrating that we reach a quality similar to the best available and commonly used merger tree tools. As a proof of concept, we use our merger tree algorithm together with a parametrized stellar-mass-to-halo-mass relation and generate a mock galaxy catalogue that shows good agreement with observational data.

SatGen: a semi-analytical satellite galaxy generator – I. The model and its application to Local-Group satellite statistics

ABSTRACT We present a semi-analytical model of satellite galaxies, SatGen, which can generate large statistical samples of satellite populations for a host halo of desired mass, redshift, and assembly history. The model combines dark matter (DM) halo merger trees, empirical relations for the galaxy–halo connection, and analytical prescriptions for tidal effects, dynamical friction, and ram-pressure stripping. SatGen emulates cosmological zoom-in hydrosimulations in certain aspects. Satellites can reside in cored or cuspy DM subhaloes, depending on the halo response to baryonic physics that can be formulated from hydrosimulations and physical modelling. The subhalo profile and the stellar mass and size of a satellite evolve depending on its tidal mass-loss and initial structure. The host galaxy can include a baryonic disc and a stellar bulge, each described by a density profile that allows analytic satellite orbit integration. SatGen complements simulations by propagating the effect of halo response found in simulated field galaxies to satellites (not properly resolved in simulations) and outperforms simulations by sampling the halo-to-halo variance of satellite statistics and overcoming artificial disruption due to insufficient resolution. As a first application, we use the model to study satellites of Milky Way (MW)- and M31-sized hosts, making it emulate simulations of bursty star formation and of smooth star formation, respectively, and to experiment with a disc potential in the host halo. We find that our model reproduces the observed satellite statistics reasonably well. Different physical recipes make a difference in satellite abundance and spatial distribution at the 25 per cent level, not large enough to be distinguished by current observations given the halo-to-halo variance. The MW/M31 disc depletes satellites by ${\sim } 20{{\ \rm per\ cent}}$ and has a subtle effect of diversifying the internal structure of satellites, which is important for alleviating certain small-scale problems. We discuss the conditions for a massive satellite to survive in MW/M31.

Self-consistent Semianalytic Modeling of Feedback during Primordial Star Formation and Reionization

Abstract We present a new semianalytic model of the formation of the first stars. Our method takes dark matter halo merger trees (including three-dimensional spatial information) from cosmological N-body simulations as input and applies analytic prescriptions to compute both the Population III and metal-enriched star formation histories. We have developed a novel method to accurately compute the major feedback processes affecting Population III star formation: H2 photodissociation from Lyman–Werner (LW) radiation, suppression of star formation due to inhomogeneous reionization, and metal enrichment via supernova winds. Our method utilizes a grid-based approach relying on fast Fourier transforms to rapidly track the LW intensity, ionization fraction, and metallicity in three dimensions throughout the simulation box. We present simulations for a wide range of astrophysical model parameters from z ≈ 30 to 6. Initially long-range LW feedback and local metal enrichment and reionization feedback dominate. However, for z ≲ 15 we find that the star formation rate density (SFRD) of Population III stars is impacted by the combination of external metal enrichment (metals from one halo polluting other pristine halos) and inhomogeneous reionization. We find that the interplay of these processes is particularly important for the Population III SFRD at z ≲ 10. Reionization feedback delays star formation long enough for metal bubbles to reach halos that would otherwise form Population III stars. Including these effects can lead to more than an order-of-magnitude decrease in the Population III SFRD at z = 6 compared to LW feedback alone.

Living with Neighbors. II. Statistical Analysis of Flybys and Mergers of Dark Matter Halos in Cosmological Simulations

Abstract We present a statistical analysis of the flybys of dark matter halos compared to mergers, using cosmological N-body simulations. We mainly focus on gravitationally interacting target halos with mass of 1010.8–1013.0 h −1 M ⊙, and their neighbors are counted only when the mass ratio is 1:3–3:1 and the distance is less than the sum of the virial radii of target and neighbor. The neighbors are divided into the flyby or merger samples if the pair’s total energy is greater or smaller, respectively, than the capture criterion with consideration of dynamical friction. The main results are as follows: (a) the flyby fraction increases by up to a factor of 50 with decreasing halo mass and by up to a factor of 400 with increasing large-scale density, while the merger fraction does not show any significant dependencies on these two parameters; (b) the redshift evolution of the flyby fraction is twofold, increasing with redshift at 0 < z < 1 and remaining constant at z > 1, while the merger fraction increases monotonically with redshift at z = 0 ∼ 4; (c) Multiple interactions with two or more neighbors are on average flyby-dominated, and their fraction has a mass and environment dependence similar to that for the flyby fraction; and (d) Given that flybys substantially outnumber mergers toward z = 0 (by a factor of five) and the multiple interactions are flyby-dominated, the flyby’s contribution to galactic evolution is stronger than ever at the present epoch, especially for less massive halos and in the higher density environment. We propose a scenario that connects the evolution of the flyby and merger fractions to the hierarchical structure formation process.

Formation of globular cluster systems – II. Impact of the cut-off of the cluster initial mass function

Observations of young star clusters reveal that the high-mass end of the cluster initial mass function (CIMF) deviates from a pure power-law and instead truncates exponentially. We investigate the effects of this truncation on the formation of globular cluster (GC) systems by updating our analytic model for cluster formation and evolution, which is based on dark matter halo merger trees coupled to empirical galactic scaling relations, and has been shown in previous work to match a wide array of observational data. The cutoff masses of $M_c=10^{6.5} M_{\odot}$ or $10^{7}M_{\odot}$ match many scaling relations: between the GC system mass and host halo mass, between the average metallicity of the GC system and host halo mass, and the distribution of cluster masses. This range of $M_c$ agrees with indirect measurements from extragalactic GC systems. Models with $M_c<10^{6.5}M_{\odot}$ cannot reproduce the observed GC metallicity and mass distributions in massive galaxies. The slope of the mass-metallicity relation for metal-poor clusters (blue tilt) for all $M_c$ models is consistent with observations within their errors, when measured using the same method. We introduce an alternative, more robust fitting method, which reveals a trend of increasing tilt slope for lower $M_c$. In our model the blue tilt arises because the metal-poor clusters form in relatively low-mass galaxies which lack sufficient cold gas to sample the CIMF at highest masses. Massive blue clusters form in progressively more massive galaxies and inherit their higher metallicity. The metal-rich clusters do not exhibit such a tilt because they form in significantly more massive galaxies, which have enough cold gas to fully sample the CIMF.

Descendants of the first stars: the distinct chemical signature of second-generation stars

Extremely metal-poor (EMP) stars in the Milky Way (MW) allow us to infer the properties of their progenitors by comparing their chemical composition to the metal yields of the first supernovae. This method is most powerful when applied to mono-enriched stars, i.e. stars that formed from gas that was enriched by only one previous supernova. We present a novel diagnostic to identify this subclass of EMP stars. We model the first generations of star formation semi-analytically, based on dark matter halo merger trees that yield MW-like halos at the present day. Radiative and chemical feedback are included self-consistently and we trace all elements up to zinc. Mono-enriched stars account for only $\sim 1\%$ of second generation stars in our fiducial model and we provide an analytical formula for this probability. We also present a novel analytical diagnostic to identify mono-enriched stars, based on the metal yields of the first supernovae. This new diagnostic allows us to derive our main results independently from the specific assumptions made regarding Pop III star formation, and we apply it to a set of observed EMP stars to demonstrate its strengths and limitations. Our results may provide selection criteria for current and future surveys and therefore contribute to a deeper understanding of EMP stars and their progenitors.

Exploring the cosmic evolution of habitability with galaxy merger trees

We combine inferred galaxy properties from a semi-analytic galaxy evolution model incorporating dark matter halo merger trees with new estimates of supernova and gamma ray burst rates as a function of metallicity from stellar population synthesis models incorporating binary interactions. We use these to explore the stellar mass fraction of galaxies irradiated by energetic astrophysical transients and its evolution over cosmic time, and thus the fraction which is potentially habitable by life like our own. We find that 18 per cent of the stellar mass in the Universe is likely to have been irradiated within the last 260 Myr, with GRBs dominating that fraction. We do not see a strong dependence of irradiated stellar mass fraction on stellar mass or richness of the galaxy environment. We consider a representative merger tree as a Local Group analogue, and find that there are galaxies at all masses which have retained a high habitable fraction (>40 per cent) over the last 6 Gyr, but also that there are galaxies at all masses where the merger history and associated star formation have rendered galaxies effectively uninhabitable. This illustrates the need to consider detailed merger trees when evaluating the cosmic evolution of habitability.

The globular cluster-dark matter halo connection.

I present a simple phenomenological model for the observed linear scaling of the stellar mass in old globular clusters (GCs) with z = 0 halo mass in which the stellar mass in GCs scales linearly with progenitor halo mass at z = 6 above a minimum halo mass for GC formation. This model reproduces the observed M GCs-M halo relation at z = 0 and results in a prediction for the minimum halo mass at z = 6 required for hosting one GC: M min(z = 6) = 1.07 × 109 M⊙. Translated to z = 0, the mean threshold mass is M halo(z = 0) ≈ 2 × 1010 M⊙. I explore the observability of GCs in the reionization era and their contribution to cosmic reionization, both of which depend sensitively on the (unknown) ratio of GC birth mass to present-day stellar mass, ξ. Based on current detections of z ≳ 6 objects with M 1500<-17, values of ξ > 10 are strongly disfavoured; this, in turn, has potentially important implications for GC formation scenarios. Even for low values of ξ, some observed high-z galaxies may actually be GCs, complicating estimates of reionization-era galaxy ultraviolet luminosity functions and constraints on dark matter models. GCs are likely important reionization sources if 5 ≲ ξ ≲ 10. I also explore predictions for the fraction of accreted versus in situ GCs in the local Universe and for descendants of systems at the halo mass threshold of GC formation (dwarf galaxies). An appealing feature of the model presented here is the ability to make predictions for GC properties based solely on dark matter halo merger trees.

The segregation of baryons and dark matter during halo assembly

The standard galaxy formation theory assumes that baryons and dark matter are initially well mixed before becoming segregated due to radiative cooling. We use non-radiative hydrodynamical simulations to explicitly examine this assumption and find that baryons and dark matter can also be segregated due to different characteristics of gas and dark matter during the buildup of the halo. As a result, baryons in many haloes do not originate from the same Lagrangian region as the dark matter. When using the fraction of corresponding dark matter and gas particles in the initial conditions (the ‘paired fraction’) as a proxy of the dark matter and gas segregation strength of a halo, on average about 25 per cent of the baryonic and dark matter of the final halo are segregated in the initial conditions. This is at odds with the assumption of the standard galaxy formation model. A consequence of this effect is that the baryons and dark matter of the same halo initially experience different tidal torques and thus their angular momentum vectors are often misaligned. The degree of the misalignment is largely preserved during later halo assembly and can be understood with the tidal torque theory. The result challenges the precision of some semi-analytical approaches that utilize dark matter halo merger trees to infer properties of gas associated with dark matter haloes.

LATE-STAGE GALAXY MERGERS IN COSMOS TO z ∼ 1

The role of major mergers in galaxy and black hole formation is not well-constrained. To help address this, we develop an automated method to identify late-stage galaxy mergers before coalescence of the galactic cores. The resulting sample of mergers is distinct from those obtained using pair-finding and morphological indicators. Our method relies on median-filtering of high-resolution images to distinguish two concentrated galaxy nuclei at small separations. This method does not rely on low surface brightness features to identify mergers, and is therefore reliable to high redshift. Using mock images, we derive statistical contamination and incompleteness corrections for the fraction of late-stage mergers. The mock images show that our method returns an uncontaminated () sample of mergers with projected separations between 2.2 and 8 out to . We apply our new method to a magnitude-limited () sample of 44,164 galaxies from the COSMOS HST/ACS catalog. Using a mass-complete sample with and , we find of systems are late-stage mergers. Correcting for incompleteness and contamination, the fractional merger rate increases strongly with redshift as , in agreement both with earlier studies and with dark matter halo merger rates. Separating the sample into star-forming and quiescent galaxies shows that the merger rate for star-forming galaxies increases strongly with redshift, , while the merger rate for quiescent galaxies is consistent with no evolution, . The merger rate also becomes steeper with decreasing stellar mass. Limiting our sample to galaxies with spectroscopic redshifts from zCOSMOS, we find that the star formation rates and X-ray selected active galactic nucleus (AGN) activity in likely late-stage mergers are higher by factors of relative to those of a control sample. Combining our sample with more widely separated pairs, we find that of star formation and of AGN activity are triggered by close encounters () or mergers, providing additional evidence that major mergers are not the only channels for star formation and black hole growth.

EFFECTS OF LARGE-SCALE ENVIRONMENT ON THE ASSEMBLY HISTORY OF CENTRAL GALAXIES

We examine whether large-scale environment affects the mass assembly history of their central galaxies. To facilitate this, we constructed dark matter halo merger trees from a cosmological N-body simulation and calculated the formation and evolution of galaxies using a semi-analytic method. We confirm earlier results that smaller halos show a notable difference in formation time with a mild dependence on large-scale environment. However, using a semi-analytic model, we found that on average the growth rate of the stellar mass of central galaxies is largely insensitive to large-scale environment. Although our results show that the star formation rate (SFR) and the stellar mass of central galaxies in smaller halos are slightly affected by the assembly bias of halos, those galaxies are faint, and the difference in the SFR is minute, and therefore it is challenging to detect it in real galaxies given the current observational accuracy. Future galaxy surveys, such as the BigBOSS experiment and the Large Synoptic Survey Telescope, which are expected to provide observational data for fainter objects, will provide a chance to test our model predictions.

Particle ejection during mergers of dark matter halos

Dark matter halos are built from accretion and merging. During merging some of the dark matter particles may be ejected with velocities higher than the escape velocity. We use both N-body simulations and single-particle smooth-field simulations to demonstrate that rapid changes to the mean field potential are responsible for such ejection, and in particular that dynamical friction plays no significant role in it. Studying a range of minor mergers, we find that typically between 5–15% of the particles from the smaller of the two merging structures are ejected. We also find that the ejected particles originate essentially from the small halo, and more specifically are particles in the small halo which pass later through the region in which the merging occurs.

Dark-matter halo mergers as a fertile environment for low-mass Population III star formation

While Population III stars are typically thought to be massive, pathways towards lower-mass Pop III stars may exist when the cooling of the gas is particularly enhanced. A possible route is enhanced HD cooling during the merging of dark-matter halos. The mergers can lead to a high ionization degree catalysing the formation of HD molecules and may cool the gas down to the cosmic microwave background (CMB) temperature. In this paper, we investigate the merging of mini-halos with masses of a few 10$^5$ M$_\odot$ and explore the feasibility of this scenario. We have performed three-dimensional cosmological hydrodynamics calculations with the ENZO code, solving the thermal and chemical evolution of the gas by employing the astrochemistry package KROME. Our results show that the HD abundance is increased by two orders of magnitude compared to the no-merging case and the halo cools down to $\sim$60 K triggering fragmentation. Based on Jeans estimates the expected stellar masses are about 10 M$_\odot$. Our findings show that the merging scenario is a potential pathway for the formation of low-mass stars.

The reionization of He II and the temperature evolution of the intergalactic medium

A number of observations suggest that He II in the intergalactic medium (IGM) was fully ionized at z ∼ 3, probably by quasi-stellar objects (QSOs). Here we construct a simple model of a QSO to study the reionization of He II and the corresponding thermal evolution of the IGM. We assume that QSOs are triggered by major mergers of dark matter halos, and the luminosity evolution of individual QSOs is described by an initial accretion stage with a constant Eddington ratio and then a power-law decay driven by long term disk evolution or fueling. Once a QSO is triggered, it immediately ionizes its surrounding area as an ionized bubble. The resulting changes in size and volume of the bubble are determined by the luminosity evolution of the central QSO. With the emergence of more and more bubbles, they eventually overlap each other and finally permeate the whole universe. During the He II reionization, the IGM temperature increases due to the photoheating by the ionization processes. Applying the bubble model and considering various heating and cooling mechanisms, we trace the thermal evolution of the IGM and obtain the average IGM temperature as a function of redshift, which is very consistent with observations. The increase in IGM temperature due to the reionization of He II may be determined more accurately in the future, which may put robust constraints on the QSO model and the physics of He II reionization.

The clustering of QSOs and the dark matter halos that host them

The spatial clustering of QSOs is an important measurable quantity which can be used to infer the properties of dark matter halos that host them. We construct a simple QSO model to explain the linear bias of QSOs measured by recent observations and explore the properties of dark matter halos that host a QSO. We assume that major mergers of dark matter halos can lead to the triggering of QSO phenomena, and the evolution of luminosity for a QSO generally shows two accretion phases, i.e., initially having a constant Eddington ratio due to the self-regulation of the accretion process when supply is sufficient, and then declining in rate with time as a power law due to either diminished supply or long term disk evolution. Using a Markov Chain Monte Carlo method, the model parameters are constrained by fitting the observationally determined QSO luminosity functions (LFs) in the hard X-ray and in the optical band simultaneously. Adopting the model parameters that best fit the QSO LFs, the linear bias of QSOs can be predicted and then compared with the observational measurements by accounting for various selection effects in different QSO surveys. We find that the latest measurements of the linear bias of QSOs from both the SDSS and BOSS QSO surveys can be well reproduced. The typical mass of SDSS QSOs at redshift 1.5 < z < 4.5 is ∼ (3 – 6) × 1012 h−1 M⊙ and the typical mass of BOSS QSOs at z ∼ 2.4 is ∼ 2 × 1012 h−1 M⊙. For relatively faint QSOs, the mass distribution of their host dark matter halos is wider than that of bright QSOs because faint QSOs can be hosted in both big halos and smaller halos, but bright QSOs are only hosted in big halos, which is part of the reason for the predicted weak dependence of the linear biases on the QSO luminosity.

On the scatter in the relation between stellar mass and halo mass: random or halo formation time dependent?

The empirical HOD model of Wang et al. 2006 fits, by construction, both the stellar mass function and correlation function of galaxies in the local Universe. In contrast, the semi-analytical models of De Lucia & Blazoit 2007 (DLB07) and Guo et al. 2011 (Guo11), built on the same dark matter halo merger trees than the empirical model, still have difficulties in reproducing these observational data simultaneously. We compare the relations between the stellar mass of galaxies and their host halo mass in the three models, and find that they are different. When the relations are rescaled to have the same median values and the same scatter as in Wang et al., the rescaled DLB07 model can fit both the measured galaxy stellar mass function and the correlation function measured in different galaxy stellar mass bins. In contrast, the rescaled Guo11 model still over-predicts the clustering of low-mass galaxies. This indicates that the detail of how galaxies populate the scatter in the stellar mass -- halo mass relation does play an important role in determining the correlation functions of galaxies. While the stellar mass of galaxies in the Wang et al. model depends only on halo mass and is randomly distributed within the scatter, galaxy stellar mass depends also on the halo formation time in semi-analytical models. At fixed value of infall mass, galaxies that lie above the median stellar mass -- halo mass relation reside in haloes that formed earlier, while galaxies that lie below the median relation reside in haloes that formed later. This effect is much stronger in Guo11 than in DLB07, which explains the over-clustering of low mass galaxies in Guo11. Our results illustrate that the assumption of random scatter in the relation between stellar and halo mass as employed by current HOD and abundance matching models may be problematic in case a significant assembly bias exists in the real Universe.

Dark matter halo merger histories beyond cold dark matter – I. Methods and application to warm dark matter

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.