Lambda Cold Dark Matter ( Delta *LCDM) is now the standard theory of structure formation in the universe. We present the first results from the new Bolshoi dissipationless cosmological Delta *LCDM ...simulation that uses cosmological parameters favored by current observations. The Bolshoi simulation was run in a volume 250 h --1 Mpc on a side using ~8 billion particles with mass and force resolution adequate to follow subhalos down to the completeness limit of V circ = 50 km s--1 maximum circular velocity. Using merger trees derived from analysis of 180 stored time steps we find the circular velocities of satellites before they fall into their host halos. Using excellent statistics of halos and subhalos (~10 million at every moment and ~50 million over the whole history) we present accurate approximations for statistics such as the halo mass function, the concentrations for distinct halos and subhalos, the abundance of halos as a function of their circular velocity, and the abundance and the spatial distribution of subhalos. We find that at high redshifts the concentration falls to a minimum value of about 4.0 and then rises for higher values of halo mass--a new result. We present approximations for the velocity and mass functions of distinct halos as a function of redshift. We find that while the Sheth-Tormen (ST) approximation for the mass function of halos found by spherical overdensity is quite accurate at low redshifts, the ST formula overpredicts the abundance of halos by nearly an order of magnitude by z = 10. We find that the number of subhalos scales with the circular velocity of the host halo as V 1/2 host, and that subhalos have nearly the same radial distribution as dark matter particles at radii 0.3-2 times the host halo virial radius. The subhalo velocity function N(> V sub) scales as V --3 circ. Combining the results of Bolshoi and Via Lactea-II simulations, we find that inside the virial radius of halos with the number of satellites is N(> V sub) = (V sub/58 km s--1)--3 for satellite circular velocities in the range 4 km s--1 < V sub < 150 km s--1.
ABSTRACT
We introduce the Uchuu suite of large high-resolution cosmological N-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion (12 8003) dark matter particles in a box ...of side-length 2.0 $\, h^{-1} \, \rm Gpc$, with particle mass of 3.27 × 108$\, h^{-1}\, \rm M_{\odot }$. The highest resolution simulation, Shin-Uchuu, consists of 262 billion (64003) particles in a box of side-length 140 $\, h^{-1} \, \rm Mpc$, with particle mass of 8.97 × 105$\, h^{-1}\, \rm M_{\odot }$. Combining these simulations, we can follow the evolution of dark matter haloes and subhaloes spanning those hosting dwarf galaxies to massive galaxy clusters across an unprecedented volume. In this first paper, we present basic statistics, dark matter power spectra, and the halo and subhalo mass functions, which demonstrate the wide dynamic range and superb statistics of the Uchuu suite. From an analysis of the evolution of the power spectra, we conclude that our simulations remain accurate from the baryon acoustic oscillation scale down to the very small. We also provide parameters of a mass–concentration model, which describes the evolution of halo concentration and reproduces our simulation data to within 5 per cent for haloes with masses spanning nearly eight orders of magnitude at redshift 0 ≤ z ≤ 14. There is an upturn in the mass–concentration relation for the population of all haloes and of relaxed haloes at z ≳ 0.5, whereas no upturn is detected at z < 0.5. We make publicly available various N-body products as part of Uchuu Data Release 1 on the Skies & Universes site.1 Future releases will include gravitational lensing maps and mock galaxy, X-ray cluster, and active galactic nucleus catalogues.
The cooling of gas in the centers of dark matter halos is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually ...calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. In contrast, halos in the hierarchical structure formation scenarios grow via multiple violent mergers and accretion along filaments, and particle orbits in the halos are highly eccentric. We study the effects of the cooling of gas in the inner regions of halos using high-resolution cosmological simulations that include gas dynamics, radiative cooling, and star formation. We find that the dissipation of gas indeed increases the density of dark matter and steepens its radial profile in the inner regions of halos compared to the case without cooling. For the first time, we test the adiabatic contraction model in cosmological simulations and find that the standard model systematically overpredicts the increase of dark matter density in the inner 5% of the virial radius. We show that the model can be improved by a simple modification of the assumed invariant from M(r)r to M(r) r, where r and r are the current and orbit-averaged particle positions. This modification approximately accounts for orbital eccentricities of particles and reproduces simulation profiles to within 10%-20%. We present analytical fitting functions that accurately describe the transformation of the dark matter profile in the modified model and can be used for interpretation of observations.
Context.
Galaxy clusters are luminous tracers of the most massive dark matter haloes in the Universe. To use them as a cosmological probe, a detailed description of the properties of dark matter ...haloes is required.
Aims.
We characterize how the dynamical state of haloes impacts the dark matter halo mass function at the high-mass end (i.e., for haloes hosting clusters of galaxies).
Methods.
We used the dark matter-only MultiDark suite of simulations and the high-mass objects
M
> 2.7 × 10
13
M
⊙
h
−1
therein. We measured the mean relations of concentration, offset, and spin as a function of dark matter halo mass and redshift. We investigated the distributions around the mean relations. We measured the dark matter halo mass function as a function of offset, spin, and redshift. We formulated a generalized mass function framework that accounts for the dynamical state of the dark matter haloes.
Results.
We confirm the recent discovery of the concentration upturn at high masses and provide a model that predicts the concentration for different values of mass and redshift with one single equation. We model the distributions around the mean values of concentration, offset, and spin with modified Schechter functions. We find that the concentration of low-mass haloes shows a faster redshift evolution compared to high-mass haloes, especially in the high-concentration regime. We find that the offset parameter is systematically smaller at low redshift, in agreement with the relaxation of structures at recent times. The peak of its distribution shifts by a factor of ∼1.5 from
z
= 1.4 to
z
= 0. The individual models are combined into a comprehensive mass function model, which predicts the mass function as a function of spin and offset. Our model recovers the fiducial mass function with ∼3% accuracy at redshift 0 and accounts for redshift evolution up to
z
∼ 1.5.
Results.
This new approach accounts for the dynamical state of the halo when measuring the halo mass function. It offers a connection with dynamical selection effects in galaxy cluster observations. This is key toward precision cosmology using cluster counts as a probe.
We present a study of the spatial distribution of dwarf satellites (or subhalos) in galactic dark matter halos using dissipationless cosmological simulations of the concordance flat cold dark matter ...(CDM) model with vacuum energy. We find that subhalos are distributed anisotropically and are preferentially located along the major axes of the triaxial mass distributions of their hosts. The Kolmogorov-Smirnov probability for drawing our simulated subhalo sample from an isotropic distribution is P sub(KS) 1.5 x 10 super(-4). An isotropic distribution of subhalos is thus not the correct null hypothesis for testing the CDM paradigm. The nearly planar distribution of observed Milky Way (MW) satellites is marginally consistent (probability 0.02) with being drawn randomly from the subhalo distribution in our simulations. Furthermore, if we select the subhalos likely to be luminous, we find a distribution that is consistent with the observed MW satellites. In fact, we show that subsamples of the subhalo population with a centrally concentrated radial distribution that is similar to that of the MW dwarfs typically exhibit a comparable degree of planarity. We explore the origin of the observed subhalo anisotropy and conclude that it is likely due to (1) the preferential accretion of satellites along filaments, often closely aligned with the major axis of the host halo, and (2) evolution of satellite orbits within the prolate, triaxial potentials typical of CDM halos. Agreement between predictions and observations requires the major axis of the outer dark matter halo of the Milky Way to be nearly perpendicular to the disk. We discuss possible observational tests of such disk-halo alignment with current large galaxy surveys.
Hierarchical cold dark matter (CDM) models predict that Milky Way-sized halos contain several hundred dense low-mass dark matter satellites (the substructure), an order of magnitude more than the ...number of observed satellites in the Local Group. If the CDM paradigm is correct, this prediction implies that the Milky Way and Andromeda are filled with numerous dark halos. To understand why these halos failed to form stars and become galaxies, we need to understand their history. We analyze the dynamical evolution of the substructure halos in a high-resolution cosmological simulation of Milky Way-sized halos in the lambda CDM cosmology. We find that about 10% of the substructure halos with the present masses ~10 super(8)-10 super(9) M sub(o) (circular velocities V sub(m) ~ 30 km s super(-1)) had considerably larger masses and circular velocities when they formed at redshifts z > ~ 2. After the initial period of mass accretion in isolation, these objects experience dramatic mass loss because of tidal stripping. Our analysis shows that strong tidal interaction is often caused by actively merging massive neighboring halos, even before the satellites are accreted by their host halo. These results can explain how the smallest dwarf spheroidal galaxies of the Local Group were able to build up a sizable stellar mass in their seemingly shallow potential wells. We propose a new model in which all the luminous dwarf spheroidals in the Local Group are descendants of the relatively massive (> ~10 super(9) M sub(o)) high-redshift systems, in which the gas could cool efficiently by atomic line emission, and which were not significantly affected by the extragalactic ultraviolet radiation. We present a simple galaxy formation model based on the trajectories extracted from the simulation, which accounts for the bursts of star formation after strong tidal shocks and the inefficiency of gas cooling in halos with virial temperatures T sub(vir) ~ 10 super(4) K. Our model reproduces the abundance, spatial distribution, and morphological segregation of the observed Galactic satellites. The results are insensitive to the redshift of reionization.
We analyze the halo occupation distribution (HOD) and two-point correlation function of galaxy-size dark matter halos using high-resolution dissipationless simulations of the concordance flat lambda ...CDM model. The halo samples include both the host halos and the "subhalos," distinct gravitationally bound halos within the virialized regions of larger host systems. We find that the HOD, the probability distribution for a halo of mass M to host a number of subhalos N, is similar to that found in semianalytic and N-body+gasdynamics studies. Its first moment, ang10Nang0 sub(M), has a complicated shape consisting of a step, a shoulder, and a power-law high-mass tail. The HOD can be described by Poisson statistics at high halo masses but becomes sub-Poisson for ang10Nang0 sub(M) ~ 4. We show that the HOD can be understood as a combination of the probability for a halo of mass M to host a central galaxy and the probability to host a given number N sub(s) of satellite galaxies. The former can be approximated by a steplike function, while the latter can be well approximated by a Poisson distribution, fully specified by its first moment. The first moment of the satellite HOD can be well described by a simple power law ang10N sub(s)ang0 proportional M super( beta ) with beta ~ 1 for a wide range of number densities, redshifts, and different power spectrum normalizations. This formulation provides a simple but accurate model for the halo occupation distribution found in simulations. At z = 0, the two-point correlation function (CF) of galactic halos can be well fitted by a power law down to approx100 h super(-1) kpc with an amplitude and slope similar to those of observed galaxies. The dependence of correlation amplitude on the number density of objects is in general agreement with results from the Sloan Digital Sky Survey. At redshifts z > ~ 1, we find significant departures from the power-law shape of the CF at small scales, where the CF steepens because of a more pronounced one-halo component. The departures from the power law may thus be easier to detect in high-redshift galaxy surveys than at the present-day epoch. They can be used to put useful constraints on the environments and formation of galaxies. If the deviations are as strong as indicated by our results, the assumption of the single power law often used in observational analyses of high-redshift clustering is dangerous and is likely to bias the estimates of the correlation length and slope of the correlation function.
The velocity shear tensor: tracer of halo alignment Libeskind, Noam I; Hoffman, Yehuda; Forero-Romero, Jaime ...
Monthly notices of the Royal Astronomical Society,
01/2013, Letnik:
428, Številka:
3
Journal Article
Recenzirano
Odprti dostop
The alignment of dark matter (DM) haloes and the surrounding large-scale structure (LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital angular momentum of subhaloes ...are investigated relative to the LSS using the eigenvectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc h
−1, deep within the non-linear regime. Knots, filaments, sheets and voids are associated with regions that are collapsing along 3, 2, 1 or 0 principal directions simultaneously. Each halo is tagged with a web classification (i.e. knot halo, filament halo, etc.) according to the nature of the collapse at the halo position. The full distribution of shear eigenvalues is found to be substantially different from that tagged to haloes, indicating that the observed velocity shear is significantly biased. We find that larger mass haloes live in regions where the shear is more isotropic, namely the expansion or collapse is more spherical. A correlation is found between the halo shape and the eigenvectors of the shear tensor, with the longest (shortest) axis of the halo shape being aligned with the slowest (fastest) collapsing eigenvector. This correlation is web independent, suggesting that the velocity shear is a fundamental tracer of the halo alignment. A similar result is found for the alignment of halo spin with the cosmic web. It has been shown that high-mass haloes exhibit a spin flip with respect to the LSS: we find that the mass at which this spin flip occurs is web dependent and not universal as suggested previously. Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, providing a possible insight into the highly correlated corotation of the Milky Way's satellite system. The present study suggests that the velocity shear tensor constitutes the natural framework for studying the directional properties of the non-linear LSS and those of haloes and galaxies.
Resonances in barred galaxies Ceverino, D.; Klypin, A.
Monthly notices of the Royal Astronomical Society,
08/2007, Letnik:
379, Številka:
3
Journal Article
Recenzirano
Odprti dostop
The inner parts of many spiral galaxies are dominated by bars. These are strong non-axisymmetric features which significantly affect orbits of stars and dark matter particles. One of the main effects ...is the dynamical resonances between galactic material and the bar. We detect and characterize these resonances in N-body models of barred galaxies by measuring angular and radial frequencies of individual orbits. We found narrow peaks in the distribution of orbital frequencies with each peak corresponding to a specific resonance. We found five different resonances in the stellar disc and two in the dark matter. The corotation resonance (CR) and the inner and outer Lindblad resonances are the most populated. The spatial distributions of particles near resonances are wide. For example, the inner Lindblad resonance is not localized at a given radius. Particles near this resonance are mainly distributed along the bar and span a wide range of radii. On the other hand, particles near the CR are distributed in two broad areas around the two stable Lagrange points. The distribution resembles a wide ring at the corotation radius. Resonances capture disc and halo material in near-resonant orbits. Our analysis of orbits in both N-body simulations and simple analytical models indicates that resonances tend to prevent the dynamical evolution of this trapped material. Only if the bar evolves as a whole, resonances drift through the phase space. In this case particles anchored near resonant orbits track the resonance shift and evolve. The criteria to ensure a correct resonant behaviour discussed by Weinberg and Katz can be achieved with few millions particles because the regions of trapped orbits near resonances are large and evolving.