Mass Accretion Histories Research Articles

Self-similar mass accretion history in scale-free simulations

ABSTRACT Using a scale-free N-body simulation generated with the abacusN-body code, we test the robustness of halo mass accretion histories via their convergence to self-similarity. We compare two halo finders, rockstar and compaso. We find superior self-similarity in halo mass accretion histories determined using rockstar, with convergence to 5 per cent or better between $\sim\!\! 10^2$ and $10^5$ particles. For compaso, we find weaker convergence over a similar region, with at least 10 per cent between $\sim\!\! 10^2$ and $10^4$ particles. Furthermore, we find that the convergence to self-similarity improves as the simulation evolves, with the largest and deepest regions of convergence appearing after the scale factor quadrupled from the time at which non-linear structures begin to form. With sufficient time evolution, halo mass accretion histories are converged to self-similarity within 5 per cent with as few as $\sim\!\! 70$ particles for compaso and within 2 per cent for as few as $\sim\!\! 30$ particles for rockstar.

Characterizing the assembly of dark matter haloes with protohalo size histories: I. Redshift evolution, relation to descendant haloes, and halo assembly bias

ABSTRACT We propose a novel method to quantify the assembly histories of dark matter haloes with the redshift evolution of the mass-weighted spatial variance of their progenitor haloes, that is, the protohalo size history. We find that the protohalo size history for each individual halo at z ∼ 0 can be described by a double power-law function. The amplitude of the fitting function strongly correlates to the central-to-total stellar mass ratios of descendant haloes. The variation of the amplitude of the protohalo size history can induce a strong halo assembly bias effect for massive haloes. This effect is detectable in observation using the central-to-total stellar mass ratio as a proxy of the protohalo size. The correlation to the descendant central-to-total stellar mass ratio and the halo assembly bias effect seen in the protohalo size are much stronger than that seen in the commonly adopted half-mass formation time derived from the mass accretion history. This indicates that the information loss caused by the compression of halo merger trees to mass accretion histories can be captured by the protohalo size history. Protohalo size thus provides a useful quantity to connect protoclusters across cosmic time and to link protoclusters with their descendant clusters in observations.

Modelling the mass accretion histories of dark matter haloes using a gamma formalism

ABSTRACT We present a physical model of the mass accretion histories (MAH) of haloes in concordance with the observed cosmic star formation rate density (CSFRD). We model the MAHs of dark matter haloes using a Gamma (Γ) functional form: $M_h(T) = \frac{M_0}{f_{0}} \, \times \frac{\gamma (\alpha _h, ~\beta _h \times (T-Th))}{\Gamma (\alpha _h)}$, where M0 is the halo mass at present time, T is time, αh and βh are parameters we explore, and f0 is the percentage of the mass of the halo at z = 0 with respect to the final mass of the halo achieved at T = ∞. We use the MAHs of haloes obtained from cosmological simulations and analytical models to constrain our model. f0 can be described by a power-law ($f_{0} = 1- c \times M_{0}^{d}$). Haloes with small masses have already on average attained most of their final masses. The average < f0 > of haloes in the Universe is >0.95 pointing to the direction that the cosmic MAH/CSFRD is saturated at our era. The average < βh > parameter (the depletion rate of the available dark matter for halo growth) is related to the dynamical time-scales of haloes. The α parameter is a power-law index of M0 and represents the early growth a halo experiences before the expansion of the Universe starts to slow it down. Finally, Th (the time that marks the co-evolution/growth of galaxies and haloes after the big bang) is found to be 150–300 million years.

A Direct Measurement of Galaxy Major and Minor Merger Rates and Stellar Mass Accretion Histories at Z < 3 Using Galaxy Pairs in the REFINE Survey

We measure the role of major and minor mergers in forming the stellar masses of galaxies over redshifts 0 < z < 3 using a combination of ∼3.25 deg2 of the deepest ground-based near-infrared imaging taken to date (Ultra Deep Survey, Ultra-VISTA, and VIDEO) as part of the collated REFINE survey. We measure the pair fraction and merger fractions for galaxy mergers of different mass ratios, and quantify the merger rate with newly measured timescales derived from the Illustris simulation as a function of redshift and merger mass ratio. For a M * > 1011 M ⊙ selection, we find that over 0 < z < 3 major mergers with mass ratios greater than 1:4 occur times on average, while minor mergers down to ratios of 1:10 occur on average times per galaxy. We also quantify the role of major and minor mergers in galaxy formation, whereby the increase in mass due to major mergers is while minor mergers account for an increase of using a M * > 1011 M ⊙ selection. We thus find that major mergers add more stellar mass to galaxies than minor mergers over this epoch. Overall, mergers will more than double the mass of massive galaxies over this epoch when selecting by stellar mass. We however find a lower increase in stellar mass when selecting by a constant number density. Finally, we compare our results to simulations, finding that minor mergers are overpredicted in Illustris and in semi-analytical models, suggesting a mismatch between observations and theory in this fundamental aspect of galaxy assembly.

Grumpy: a simple framework for realistic forward modelling of dwarf galaxies

ABSTRACT We present a simple regulator-type framework designed specifically for modelling formation of dwarf galaxies. Despite its simplicity, when coupled with realistic mass accretion histories of haloes from simulations and reasonable choices for model parameter values, the framework can reproduce a remarkably broad range of observed properties of dwarf galaxies over seven orders of magnitude in stellar mass. In particular, we show that the model can simultaneously match observational constraints on the stellar mass–halo mass relation, as well as observed relations between stellar mass and gas phase and stellar metallicities, gas mass, size, and star formation rate, as well as general form and diversity of star formation histories of observed dwarf galaxies. The model can thus be used to predict photometric properties of dwarf galaxies hosted by dark matter haloes in N-body simulations, such as colours, surface brightnesses, and mass-to-light ratios and to forward model observations of dwarf galaxies. We present examples of such modelling and show that colours and surface brightness distributions of model galaxies are in good agreement with observed distributions for dwarfs in recent observational surveys. We also show that in contrast with the common assumption, the absolute magnitude–halo mass relation is generally predicted to have a non-power law form in the dwarf regime, and that the fraction of haloes that host detectable ultra-faint galaxies is sensitive to reionization redshift (zrei) and is predicted to be consistent with observations for zrei ≲ 9.

Is there enough star formation in simulated protoclusters?

ABSTRACT As progenitors of the most massive objects, protoclusters are key to tracing the evolution and star formation history of the Universe, and are responsible for ${\gtrsim }\, 20$ per cent of the cosmic star formation at $z\, {\gt }\, 2$. Using a combination of state-of-the-art hydrodynamical simulations and empirical models, we show that current galaxy formation models do not produce enough star formation in protoclusters to match observations. We find that the star formation rates (SFRs) predicted from the models are an order of magnitude lower than what is seen in observations, despite the relatively good agreement found for their mass-accretion histories, specifically that they lie on an evolutionary path to become Coma-like clusters at $z\, {\simeq }\, 0$. Using a well-studied protocluster core at $z\, {=}\, 4.3$ as a test case, we find that star formation efficiency of protocluster galaxies is higher than predicted by the models. We show that a large part of the discrepancy can be attributed to a dependence of SFR on the numerical resolution of the simulations, with a roughly factor of 3 drop in SFR when the spatial resolution decreases by a factor of 4. We also present predictions up to $z\, {\simeq }\, 7$. Compared to lower redshifts, we find that centrals (the most massive member galaxies) are more distinct from the other galaxies, while protocluster galaxies are less distinct from field galaxies. All these results suggest that, as a rare and extreme population at high z, protoclusters can help constrain galaxy formation models tuned to match the average population at $z\, {\simeq }\, 0$.

Rapid filamentary accretion as the origin of extended thin discs

ABSTRACT Galactic outflows driven by stellar feedback are crucial for explaining the inefficiency of star formation in galaxies. Although strong feedback can promote the formation of galactic discs by limiting star formation at early times and removing low angular momentum (AM) gas, it is not understood how the same feedback can result in diverse objects such as elliptical galaxies or razor thin spiral galaxies. We investigate this problem using cosmological zoom-in simulations of two galaxies forming within 1012 M⊙ haloes with almost identical mass accretion histories and halo spin parameters. However, the two resulting galaxies end up with very different bulge-to-disc ratios at z = 0. At z &amp;gt; 1.5, the two galaxies feature a surface density of star formation ΣSFR ≃ 10 M⊙ yr−1 kpc−2, leading to strong outflows. After the last starburst episode, both galaxies feature a dramatic gaseous disc growth from 1 to 5 kpc during 1 Gyr, a decisive event we dub ‘the Grand Twirl’. After this event, the evolutionary tracks diverge strongly, with one galaxy ending up as a bulge-dominated galaxy, whereas the other ends up as a disc-dominated galaxy. The origins of this dichotomy are the AM of the accreted gas, and whether it adds constructively to the initial disc angular momentum. The build-up of this extended disc leads to a rapid lowering of ΣSFR by over two orders of magnitude with ΣSFR ≲ 0.1 M⊙ yr−1 kpc−2, in remarkable agreement with what is derived from Milky Way stellar populations. As a consequence, supernovae explosions are spread out and cannot launch galactic outflows anymore, allowing for the persistence of a thin, gently star-forming, extended disc.

The past history of galaxy clusters told by their present neighbours

ABSTRACT Galaxy clusters can play a key role in modern cosmology, provided their evolution is properly understood. However, observed clusters give us only a single timeframe of their dynamical state. Therefore, finding present observable data of clusters that are well correlated to their assembly history constitutes an inestimable tool for cosmology. Former studies correlating environmental descriptors of clusters to their formation history are dominated by halo mass–environment relations. This paper presents a mass-free correlation between the present neighbour distribution of cluster-size haloes and the latter mass assembly history. From the Big Multidark simulation, we extract two large samples of random haloes with masses ranging from Virgo to Coma cluster sizes. Additionally, to find the main environmental culprit for the formation history of the Virgo cluster, we compare the Virgo-size haloes to 200 Virgo-like haloes extracted from simulations that resemble the local Universe. The number of neighbours at different cluster-centric distances permits discriminating between clusters with different mass accretion histories. Similarly to Virgo-like haloes, clusters with numerous neighbours within a distance of about two times their virial radius experience a transition at z ≈ 1 between an active period of mass accretion, relative to the mean, and a quiet history. In contrary, clusters with few neighbours share an opposite trend: from passive to active assembly histories. Additionally, clusters with massive companions within about four times their virial radius tend to have recent active merging histories. Therefore, the radial distribution of cluster neighbours provides invaluable insights into the past history of these objects.

Scatter in Sunyaev–Zel’dovich effect scaling relations explained by inter-cluster variance in mass accretion histories

ABSTRACT X-ray and microwave cluster scaling relations are immensely valuable for cosmological analysis. However, their power is limited by astrophysical systematics that bias mass estimates and introduce additional scatter. Turbulence injected into the intracluster medium via mass assembly contributes substantially to cluster non-thermal pressure support, a significant source of such uncertainties. We use an analytical model to compute the assembly-driven non-thermal pressure profiles of haloes based on Monte Carlo-generated accretion histories. We introduce a fitting function for the average non-thermal pressure fraction profile, which exhibits minimal dependence on redshift at fixed peak height. Using the model, we predict deviations from self-similarity and the intrinsic scatter in the Sunyaev–Zel’dovich effect observable-mass scaling relation (YSZ − M) due solely to inter-cluster variation in mass accretion histories. We study the dependence of YSZ − M on aperture radius, cosmology, redshift, and mass limit. The model predicts $5-9{{\ \rm per\ cent}}$ scatter in YSZ − M at z = 0, increasing as the aperture used to compute YSZ increases from R500c to 5R500c. The predicted scatter lies slightly below that of studies based on non-radiative hydro-simulations, illustrating that assembly history variance is likely responsible for a substantial fraction of scatter in YSZ − M. This should be regarded as a lower bound, which will likely increase with the use of an updated gas density model that incorporates a more realistic response to halo assembly. As redshift increases, YSZ − M deviates more from self-similarity and scatter increases. We show that the YSZ − M residuals correlate strongly with the recent halo mass accretion rate, potentially providing an opportunity to infer the latter.

Elliptical galaxies/stellar halos connection

AbstractCosmological simulations predict that early-type galaxies (ETGs) are the results of extended mass accretion histories. The latter are characterized by different numbers of mergers, mergers’ mass ratios and gas fractions, and timing. Depending on the sequence and nature of these mergers that follow the first phase of the in-situ star formation, these accretion histories may lead to ETGs that have low or high mass halos, and that rotate fast or slow. Since the stellar halos maintain the fossil records of the events that led to their formation, a discontinuity may be in place between the inner regions of ETGs and their outer halos, because the time required for the halos’ stars to exchange their energies and momenta is very long compared with the age of these systems. Exquisite deep photometry and extended spectroscopy for significant samples of ETGs are then used to quantify the occurrence and significance of such a transition in the galaxies’ structural and kinematical parameters. Once this transition radius is measured, its dependency with the effective radius of the galaxies’ light distribution and total stellar masses can be investigated. Such correlations can then be compared with the predictions of accreted, i.e. ex-situ vs. in-situ components from cosmological simulations to validate such models.

The Stellar-to-halo Mass Ratios of Passive and Star-forming Galaxies at z ∼ 2–3 from the SMUVS Survey

Abstract In this work, we use measurements of galaxy stellar mass and two-point angular correlation functions to constrain the stellar-to-halo mass ratios (SHMRs) of passive and star-forming galaxies at z ∼ 2–3, as identified in the Spitzer Matching Survey of the UltraVISTA ultra-deep Stripes. We adopt a sophisticated halo modeling approach to statistically divide our two populations into central and satellite galaxies. For central galaxies, we find that the normalization of the SHMR is greater for our passive population. Through the modeling of Λ cold dark matter halo mass accretion histories, we show that this can only arise if the conversion of baryons into stars was more efficient at higher redshifts and additionally that passive galaxies can be plausibly explained as residing in halos with the highest formation redshifts (i.e., those with the lowest accretion rates) at a given halo mass. At a fixed stellar mass, satellite galaxies occupy host halos with a greater mass than central galaxies, and we find further that the fraction of passive galaxies that are satellites is higher than for the combined population. This, and our derived satellite quenching timescales, combined with earlier estimates from the literature, support dynamical/environmental mechanisms as the dominant process for satellite quenching at z ≲ 3.

Reconstructing the mass accretion histories of nearby red nuggets with their globular cluster systems

AbstractIt is generally recognized that massive galaxies form through a combination of in-situ collapse and ex-situ accretion. The in-situ component forms early, where gas collapse and compaction leads to the formation of massive compact systems (blue and red “nuggets”) seen at z &gt; 1. The subsequent accretion of satellites brings in ex-situ material, growing these nuggets in size and mass to appear as the massive early-type galaxies (ETGs) we see locally. Due to stochasticity in the accretion process, in a few rare cases a red nugget will evolve to the present day having undergone little ex-situ mass accretion. The resulting massive, compact and ancient objects have been termed “relic galaxies”. Detailed stellar population and kinematic analyses are required to characterise these systems. However, an additional crucial aspect lies in determining the fraction of ex-situ mass they have accreted since their formation. Globular cluster systems can be used to constrain this fraction, since the oldest and most metal-poor globular clusters in massive galaxies are primarily an accreted, ex-situ population. Models for the formation of relic galaxies and their globular cluster systems suggest that, due to their early compaction and limited accretion of dark-matter dominated satellites, relic galaxies should have characteristically low dark-matter mass fractions compared to ETGs of the same stellar mass.

Sensitivity of dark matter haloes to their accretion histories

We apply our recently proposed "quadratic genetic modification" approach to generating and testing the effects of alternative mass accretion histories for a single $\Lambda$CDM halo. The goal of the technique is to construct different formation histories, varying the overall contribution of mergers to the fixed final mass. This enables targeted studies of galaxy and dark matter halo formation's sensitivity to the smoothness of mass accretion. Here, we focus on two dark matter haloes, each with four different mass accretion histories. We find that the concentration of both haloes systematically decreases as their merger history becomes smoother. This causal trend tracks the known correlation between formation time and concentration parameters in the overall halo population. At fixed formation time, we further establish that halo concentrations are sensitive to the order in which mergers happen. This ability to study an individual halo's response to variations in its history is highly complementary to traditional methods based on emergent correlations from an extended halo population.

Modelling baryonic effects on galaxy cluster mass profiles

Gravitational lensing is a powerful probe of the mass distribution of galaxy clusters and cosmology. However, accurate measurements of the cluster mass profiles are limited by uncertainties in cluster astrophysics. In this work, we present a physically motivated model of baryonic effects on the cluster mass profiles, which self-consistently takes into account the impact of baryons on the concentration as well as mass accretion histories of galaxy clusters. We calibrate this model using the Omega500 hydrodynamical cosmological simulations of galaxy clusters with varying baryonic physics. Our model will enable us to simultaneously constrain cluster mass, concentration, and cosmological parameters using stacked weak lensing measurements from upcoming optical cluster surveys.

The Galaxy–Halo Connection for as Revealed by the Spitzer Matching Survey of the UltraVISTA Ultra-deep Stripes

Abstract The Spitzer Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS) provides unparalleled depth at 3.6 and 4.5 μm over ∼0.66 deg2 of the COSMOS field, allowing precise photometric determinations of redshift and stellar mass. From this unique data set we can connect galaxy samples, selected by stellar mass, to their host dark matter halos for , filling in a large hitherto unexplored region of the parameter space. To interpret the observed galaxy clustering, we use a phenomenological halo model, combined with a novel method to account for uncertainties arising from the use of photometric redshifts. We find that the satellite fraction decreases with increasing redshift and that the clustering amplitude (e.g., comoving correlation length/large-scale bias) displays monotonic trends with redshift and stellar mass. Applying ΛCDM halo mass accretion histories and cumulative abundance arguments for the evolution of stellar mass content, we propose pathways for the coevolution of dark matter and stellar mass assembly. Additionally, we are able to estimate that the halo mass at which the ratio of stellar-to-halo mass is maximized is at . This peak halo mass is here inferred for the first time from stellar mass-selected clustering measurements at , and it implies a mild evolution of this quantity for , consistent with constraints from abundance-matching techniques.

Einasto Profiles and the Dark Matter Power Spectrum

Abstract We study the mass accretion histories (MAHs) and density profiles of dark matter halos using N-body simulations of self-similar gravitational clustering from scale-free power spectra, P(k)∝kn. We pay particular attention to the density profile curvature, which we characterize using the shape parameter, α, of an Einasto profile. In agreement with previous findings our results suggest that, despite vast differences in their MAHs, the density profiles of virialized halos are remarkably alike. Nonetheless, clear departures from self-similarity are evident: for a given spectral index, α increases slightly but systematically with “peak height”, ν ≡ δsc/σ(M, z), regardless of mass or redshift. More importantly, however, the “α − ν” relation depends on n: the steeper the initial power spectrum, the more gradual the curvature of both the mean MAHs and mean density profiles. These results are consistent with previous findings connecting the shapes of halo mass profiles and MAHs and imply that dark matter halos are not structurally self-similar but, through the merger history, retain a memory of the linear density field from which they form.

The mass–concentration–redshift relation of cold and warm dark matter haloes

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $\Lambda$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.

THE CRITERION OF SUPERNOVA EXPLOSION REVISITED: THE MASS ACCRETION HISTORY

By performing neutrino-radiation hydrodynamic simulations in spherical symmetry (1D) and axial symmetry (2D) with different progenitor models by Woosley & Heger from 12 to 100 M⊙, we find that all 1D runs fail to produce an explosion and several 2D runs succeed. The difference in the shock evolutions for different progenitors can be interpreted by the difference in their mass accretion histories, which are in turn determined by the density structures of progenitors. The mass accretion history has two phases in the majority of the models: the earlier phase, in which the mass accretion rate is high and rapidly decreasing, and the later phase, with a low and almost constant accretion rate. They are separated by the so-called turning point, the origin of which is a change of the accreting layer. We argue that shock revival will most likely occur around the turning point and hence that its location in the plane will be a good measure for the possibility of shock revival: if the turning point lies above the critical curve and the system stays there for a long time, shock revival will obtain. In addition, we develop a phenomenological model to approximately evaluate the trajectories in the plane, which, after calibrating free parameters by a small number of 1D simulations, reproduces the location of the turning point reasonably well by using the initial density structure of progenitor alone. We suggest the application of the phenomenological model to a large collection of progenitors in order to infer without simulations which ones are more likely to explode.

Major mergers going Notts: challenges for modern halo finders

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.

The accretion history of dark matter haloes – II. The connections with the mass power spectrum and the density profile

We explore the relation between the structure and mass accretion histories of dark matter halos using a suite of cosmological simulations. We confirm that the formation time, defined as the time when the virial mass of the main progenitor equals the mass enclosed within the scale radius, correlates strongly with concentration. We provide a semi-analytic model for halo mass history that combines analytic relations with fits to simulations. This model has the functional form, $M(z) = M_{0}(1+z)^{\alpha}e^{\beta z}$, where the parameters $\alpha$ and $\beta$ are directly correlated with concentration. We then combine this model for the halo mass history with the analytic relations between $\alpha$, $\beta$ and the linear power spectrum derived by Correa et al. (2014) to establish the physical link between halo concentration and the initial density perturbation field. Finally, we provide fitting formulas for the halo mass history as well as numerical routines, we derive the accretion rate as a function of halo mass, and we demonstrate how the halo mass history depends on cosmology and the adopted definition of halo mass.