No Arabic abstract
In the far future of an accelerating LCDM cosmology, the cosmic web of large-scale structure consists of a set of increasingly isolated halos in dynamical equilibrium. We examine the approach of collisionless dark matter to hydrostatic equilibrium using a large N-body simulation evolved to scale factor a = 100, well beyond the vacuum--matter equality epoch, a_eq ~ 0.75, and 53/h Gyr into the future for a concordance model universe (Omega_m ~ 0.3, Omega_Lambda ~ 0.7). The radial phase-space structure of halos -- characterized at a < a_eq by a pair of zero-velocity surfaces that bracket a dynamically active accretion region -- simplifies at a > 10 a_eq when these surfaces merge to create a single zero-velocity surface, clearly defining the halo outer boundary, rhalo, and its enclosed mass, mhalo. This boundary approaches a fixed physical size encompassing a mean interior density ~ 5 times the critical density, similar to the turnaround value in a classical Einstein-deSitter model. We relate mhalo to other scales currently used to define halo mass (m200, mvir, m180b) and find that m200 is approximately half of the total asymptotic cluster mass, while m180b follows the evolution of the inner zero velocity surface for a < 2 but becomes much larger than the total bound mass for a > 3. The radial density profile of all bound halo material is well fit by a truncated Hernquist profile. An NFW profile provides a somewhat better fit interior to r200 but is much too shallow in the range r200 < r < rhalo.
Employing hydrodynamic simulations of structure formation in a LCDM cosmology, we study the history of cosmic star formation from the dark ages at redshift z~20 to the present. In addition to gravity and ordinary hydrodynamics, our model includes radiative heating and cooling of gas, star formation, supernova feedback, and galactic winds. By making use of a comprehensive set of simulations on interlocking scales and epochs, we demonstrate numerical convergence of our results on all relevant halo mass scales, ranging from 10^8 to 10^15 Msun/h. The predicted density of cosmic star formation is broadly consistent with measurements, given observational uncertainty. From the present epoch, it gradually rises by about a factor of ten to a peak at z~5-6, which is beyond the redshift range where it has been estimated observationally. 50% of the stars are predicted to have formed by redshift z~2.1, and are thus older than 10.4 Gyr, while only 25% form at redshifts lower than z~1. The mean age of all stars at the present is about 9 Gyr. Our model predicts a total stellar density at z=0 of Omega_*=0.004, corresponding to about 10% of all baryons being locked up in long-lived stars, in agreement with recent determinations of the luminosity density of the Universe. We determine the multiplicity function of cosmic star formation as a function of redshift; i.e. the distribution of star formation with respect to halo mass. We also briefly examine possible implications of our predicted star formation history for reionisation of hydrogen in the Universe. We find that the star formation rate predicted by the simulations is sufficient to account for hydrogen reionisation by z~6, but only if a high escape fraction close to unity is assumed. (abridged)
We study the concentration of dark matter halos and its evolution in N-body simulations of the standard LCDM cosmology. The results presented in this paper are based on 4 large N-body simulations with about 10 billion particles each: the Millennium-I and II, Bolshoi, and MultiDark simulations. The MultiDark (or BigBolshoi) simulation is introduced in this paper. This suite of simulations with high mass resolution over a large volume allows us to compute with unprecedented accuracy the concentration over a large range of scales (about six orders of magnitude in mass), which constitutes the state-of-the-art of our current knowledge on this basic property of dark matter halos in the LCDM cosmology. We find that there is consistency among the different simulation data sets. We confirm a novel feature for halo concentrations at high redshifts: a flattening and upturn with increasing mass. The concentration c(M,z) as a function of mass and the redshift and for different cosmological parameters shows a remarkably complex pattern. However, when expressed in terms of the linear rms fluctuation of the density field sigma(M,z), the halo concentration c(sigma) shows a nearly-universal simple U-shaped behaviour with a minimum at a well defined scale at sigma=0.71. Yet, some small dependences with redshift and cosmology still remain. At the high-mass end (sigma < 1) the median halo kinematic profiles show large signatures of infall and highly radial orbits. This c-sigma(M,z) relation can be accurately parametrized and provides an analytical model for the dependence of concentration on halo mass. When applied to galaxy clusters, our estimates of concentrations are substantially larger -- by a factor up to 1.5 -- than previous results from smaller simulations, and are in much better agreement with results of observations. (abridged)
We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies.
We present a comparison of the statistical properties of dark matter halo merger trees extracted from the Millennium Simulation with Extended Press-Schechter (EPS) formalism and the related GALFORM Monte-Carlo method for generating ensembles of merger trees. The volume, mass resolution and output frequency make the Millennium Simulation a unique resource for the study of the hierarchical growth of structure. We construct the merger trees of present day friends-of-friends groups and calculate a variety of statistics that quantify the masses of their progenitors as a function of redshift; accretion rates; and the redshift distribution of their most recent major merger. We also look in the forward direction and quantify the present day mass distribution of halos into which high redshift progenitors of a specific mass become incorporated. We find that EPS formalism and its Monte-Carlo extension capture the qualitative behaviour of all these statistics but, as redshift increases they systematically underestimate the masses of the most massive progenitors. This shortcoming is worst for the Monte-Carlo algorithm. We present a fitting function to a scaled version of the progenitor mass distribution and show how it can be used to make more accurate predictions of both progenitor and final halo mass distributions.
We use a combination of a cosmological N-body simulation of the concordance Lambda cold dark matter (LCDM) paradigm and a semi-analytic model of galaxy formation to investigate the spin development of central supermassive black holes (BHs) and its relation to the BH host galaxy properties. In order to compute BH spins, we use the alpha-model of Shakura & Sunyaev and consider the King et al. warped disc alignment criterion. The orientation of the accretion disc is inferred from the angular momentum of the source of accreted material, which bears a close relationship to the large-scale structure in the simulation. We find that the final BH spin depends almost exclusively on the accretion history and only weakly on the warped disc alignment. The main mechanisms of BH spin-up are found to be gas cooling processes and disc instabilities, a result that is only partially compatible with Monte-Carlo models where the main spin-up mechanisms are major mergers and disc instabilities; the latter results are reproduced when implementing randomly oriented accretion discs in our model. Regarding the BH population, we find that more massive BHs, which are hosted by massive ellipticals, have higher spin values than less-massive BHs, hosted by spiral galaxies. We analyse whether gas accretion rates and BH spins can be used as tracers of the radio loudness of active galactic nuclei (AGN). We find that the current observational indications of an increasing trend of radio-loud AGN fractions with stellar and BH mass can be easily obtained when placing lower limits on the BH spin, with a minimum influence from limits on the accretion rates; a model with random accretion disc orientations is unable to reproduce this trend. (ABRIDGED)