Academia.edu no longer supports Internet Explorer.
To browse Academia.edu and the wider internet faster and more securely, please take a few seconds to upgrade your browser.
2015, Physical Review D
We perform a study of the nonlinear clustering of matter in the late-forming dark matter (LFDM) scenario in which dark matter results from the transition of a nonminimally coupled scalar field from radiation to collisionless matter. A distinct feature of this model is the presence of a damped oscillatory cutoff in the linear matter power spectrum at small scales. We use a suite of high-resolution N-body simulations to study the imprints of LFDM on the nonlinear matter power spectrum, the halo mass and velocity functions and the halo density profiles. The model largely satisfies highredshift matter power spectrum constraints from Lyman-α forest measurements, while it predicts suppressed abundance of low-mass halos (∼ 10 9 -10 10 h -1 M ) at all redshifts compared to a vanilla ΛCDM model. The analysis of the LFDM halo velocity function shows a better agreement than the ΛCDM prediction with the observed abundance of low-velocity galaxies in the local volume. Halos with mass M 10 11 h -1 M show minor departures of the density profiles from ΛCDM expectations, while smaller-mass halos are less dense, consistent with the fact that they form later than their ΛCDM counterparts.
Monthly Notices of the Royal Astronomical Society, 1999
We develop a model for the growth of dark matter halos and use it to study their evolved density profiles. In this model, halos are spherical and form by quiescent accretion of matter in clumps, called satellites. The halo mass as a function of redshift is given by the mass of the most massive progenitor, and is determined from Monte-Carlo realizations of the merger-history tree. Inside the halo, satellites move under the action of the gravitational force of the halo and a dynamical friction drag force. The associated equation of motion is solved numerically. The energy lost to dynamical friction is transferred to the halo in the form of kinetic energy. As they sink into the halo, satellites continually lose matter as a result of tidal stripping. The stripped matter moves inside the halo free of dynamical friction. The evolved density profiles are steeper than those obtained by assuming that, once they have been accreted onto the parent halo, satellites remain at a fixed distance from the halo center. We find that the final density profile depends mainly on the rate of infall of matter onto the halo. This, in turn, depends on the initial fluctuation field as well as on cosmology. For mass scales where the effective spectral index of the initial density field is less than −1, the model predicts a profile which can only approximately be matched by the one parameter family of curves suggested by Navarro, Frenk and White (1997). For scale-free power-spectra with initial slope n, the density profile within about 1% of the virial radius is ρ ∝ r −β , with 3(3 + n)/(5 + n) ≤ β ≤ 3(3 + n)/(4 + n).
1998
We develop a model for the growth of dark matter halos and use it to study their evolved density profiles. In this model, halos are spherical and form by quiescent accretion of matter in clumps, called satellites. The halo mass as a function of redshift is given by the mass of the most massive progenitor, and is determined from Monte-Carlo realizations of the merger-history tree. Inside the halo, satellites move under the action of the gravitational force of the halo and a dynamical friction drag force. The associated equation of motion is solved numerically. The energy lost to dynamical friction is transferred to the halo in the form of kinetic energy. As they sink into the halo, satellites continually lose matter as a result of tidal stripping. The stripped matter moves inside the halo free of dynamical friction. The evolved density profiles are steeper than those obtained by assuming that, once they have been accreted onto the parent halo, satellites remain at a fixed distance from the halo center. We find that the final density profile depends mainly on the rate of infall of matter onto the halo. This, in turn, depends on the initial fluctuation field as well as on cosmology. For mass scales where the effective spectral index of the initial density field is less than -1, the model predicts a profile which can only approximately be matched by the one parameter family of curves suggested by Navarro, Frenk and White (1997). For scale-free power-spectra with initial slope $n$, the density profile within about 1% of the virial radius is $\rho\propto r^{-\beta}$, with $3(3+n)/(5+n)\le\beta\le 3(3+n)/(4+n)$.
High resolution N-body simulations have all but converged on a common empirical form for the shape of the density profiles of halos, but the full understanding of the underlying physics of halo formation has eluded them so far. We investigate the formation and structure of dark matter halos using analytical and semi-analytical techniques. Our halos are formed via an extended secondary infall model (ESIM); they contain secondary perturbations and hence random tangential and radial motions which affect the halo's evolution at it undergoes shell-crossing and virialization. Even though the density profiles of NFW and ESIM halos are different their phase-space density distributions are the same: ρ/σ 3 ∝ r −α , with α = 1.875 over ∼ 3 decades in radius. We use two approaches to try to explain this "universal" slope: (1) The Jeans equation analysis yields many insights, however, does not answer why α = 1.875. (2) The secondary infall model of the 1960's and 1970's, augmented by "thermal motions" of particles does predict that halos should have α = 1.875. However, this relies on assumptions of spherical symmetry and slow accretion. While for ESIM halos these assumptions are justified, they most certainly break down for simulated halos which forms hierarchically. We speculate that our argument may apply to an "on-average" formation scenario of halos within merger-driven numerical simulations, and thereby explain why α = 1.875 for NFW halos. Thus, ρ/σ 3 ∝ r −1.875 may be a generic feature of violent relaxation.
Astrophysical Journal, 2008
We report a series of high-resolution cosmological N-body simulations designed to explore the formation and properties of dark matter halos with masses close to the damping scale of the primordial power spectrum of density fluctuations. We further investigate the effect that the addition of a random component, v rms , into the particle velocity field has on the structure of halos. We adopted as a fiducial model the Λ warm dark matter (ΛWDM) cosmology with a non-thermal sterile neutrino mass of 0.5 keV. The filtering mass corresponds then to M f = 2.6 × 10 12 h −1 M ⊙ . Halos of masses close to M f were simulated with several million of particles. The results show that, on one hand, the inner density slope of these halos (at radii < ∼ 0.02 the virial radius R v ) is systematically steeper than the one corresponding to the Navarro-Frenk-White (NFW) fit or to the cold dark matter counterpart. On the other hand, the overall density profile (radii larger than 0.02 R v ) is less curved and less concentrated than the NFW fit, with an outer slope shallower than -3. For simulations with v rms , the inner halo density profiles flatten significantly at radii smaller than 2-3 h −1 kpc ( < ∼ 0.010R v − 0.015R v ). A constant density core is not detected in our simulations, with the exception of one halo for which the flat core radius is ≈ 1h −1 kpc. Nevertheless, if "cored" density profiles are used to fit the halo profiles, the inferred core radii are ≈ (0.1 − 0.8)h −1 kpc, in rough agreement with theoretical predictions based on phase-space constrains, and on dynamical models of warm gravitational collapse. A reduction of v rms by a factor of 3 produces a modest decrease in core radii, less than a factor of 1.5. We discuss the extension of our results into several contexts, for example, to the structure of the cold DM micro-halos at the damping scale of this model.
Monthly Notices of The Royal Astronomical Society, 2005
Using six high resolution dissipationless simulations with a varying box size in a flat LCDM universe, we study the mass and redshift dependence of dark matter halo shapes for M_vir = 9.0e11 - 2.0e14, over the redshift range z=0-3, and for two values of sigma_8=0.75 and 0.9. Remarkably, we find that the redshift, mass, and sigma_8 dependence of the mean smallest-to-largest axis ratio of halos is well described by the simple power-law relation <s> = (0.54 +- 0.02)(M_vir/M_*)^(-0.050 +- 0.003), where s is measured at 0.3 R_vir and the z and sigma_8 dependences are governed by the characteristic nonlinear mass, M_*=M_*(z,sigma_8). We find that the scatter about the mean s is well described by a Gaussian with sigma ~ 0.1, for all masses and redshifts. We compare our results to a variety of previous works on halo shapes and find that reported differences between studies are primarily explained by differences in their methodologies. We address the evolutionary aspects of individual halo shapes by following the shapes of the halos through ~100 snapshots in time. We determine the formation scalefactor a_c as defined by Wechsler et al. (2002) and find that it can be related to the halo shape at z = 0 and its evolution over time.
2007
(Context) In a Universe dominated by dark matter, halos are continuously accreting mass (violently or not) and such mechanism affects their dynamical state. (Aims) The evolution of dark matter halos in phase-space, and using the phase-space density indicator Q=rho/sigma^3 as a tracer, is discussed. (Methods) We have performed cosmological N-body simulations from which we have carried a detailed study of the evolution of ~35 dark halos in the interval 0<z<10. (Results)The follow up of individual halos indicates two distinct evolutionary phases. First, an early and fast decrease of Q associated to virialization after the gravitational collapse takes place. The nice agreement between simulated data and theoretical expectations based on the spherical collapse model support such a conjecture. The late and long period where a slow decrease of the phase-space density occurs is related to accretion and merger episodes. The study of some merger events in the phase-space (radial velocity versus radial distance) reveals the formation of structures quite similar to caustics generated in secondary infall models of halo formation. After mixing in phase-space, halos in quasi-equilibrium have flat-topped velocity distributions (negative kurtosis) with respect to Gaussians. The effect is more noticiable for captured satellites and/or substructures than for the host halo.
EAS Publications Series, 2006
The formation and structure of dark matter halos is studied by constrained simulations. A series of experiments of the formation of a 10 12 h −1 M⊙ halo is designed to study the dependence of the density profile on its merging history. We find that the halo growth consist of several quiescent phases intermitted by violent events, with the density well approximated by the NFW profile during the former phases. We find that (1) the NFW scale radius Rs stays constant during the quiescent phase and grows abruptly during the violent one. In contrast, the virial radius grows linearly during the quiescent and abruptly during the violent phases. (2) The central density stays unchanged during the quiescent phase while dropping abruptly during the violent phase, and it does not reflect the formation time of the halo. (3) The clear separation of the evolution of an individual halo into quiescent and violent phases implies that its entire evolution cannot be fitted by simple scaling relations.
Monthly Notices of the Royal Astronomical Society, 2003
We use a high resolution ΛCDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of fifteen, in a 50 h −1 Mpc volume containing 80 million particles. Our low redshift results allow us to probe low σ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth and Tormen mass function provides an excellent match to all of our data except for redshifts of ten and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 percent for the 10 10 − 10 11 M ⊙ haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth and Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high redshift (z ≃ 6) QSOs, which are thought to be hosted by haloes at 5σ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50%. We test the sensitivity of our results to force accuracy, starting redshift, and halo finding algorithm.
The Astrophysical Journal, 2002
Current data suggest that the central mass densities ρ 0 and phase-space densities Q ≡ ρ 0 /σ 3 V of cosmological halos in the present universe are correlated with their velocity dispersions σ V over a very wide range of σ V from less than 10 to more than 1000 km s −1 . Such correlations are an expected consequence of the statistical correlation of the formation epochs of virialized objects in the CDM model with their masses; the smaller-mass halos typically form first and merge to form larger-mass halos later. We have derived the Q − σ V and ρ 0 − σ V correlations for different CDM cosmologies and compared the predicted correlations with the observed properties of a sample of low-redshift halos ranging in size from dwarf spheroidal galaxies to galaxy clusters. Our predictions are generally consistent with the data, with preference for the currently-favored, flat ΛCDM model. Such a comparison serves to test the basic CDM paradigm while constraining the background cosmology and the power-spectrum of primordial density fluctuations, including larger wavenumbers than have previously been constrained.
Monthly Notices of the Royal Astronomical Society, 1998
We examine the properties of dark matter halos within a rich galaxy cluster using a high resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 halos with circular velocities larger than 80 km s −1 within the cluster's virial radius of 2 Mpc. This enables an unprecedented study of the statistical properties of a large sample of dark matter halos evolving in a dense environment. The cumulative fraction of mass attached to these halos varies from 0% at 200 kpc, to 13% at the virial radius. Even at this resolution the overmerging problem persists; halos that pass within 200 kpc of the cluster center are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor halos. The median ratio of apocentric to pericentric radii is 6:1; the orbital distribution is close to isotropic, circular orbits are rare, radial orbits are common. The orbits of halos are unbiased with respect to both position within the cluster and with the orbits of the smooth dark matter background and no velocity bias is detected. The tidal radii of surviving halos are generally well-fit using the simple analytic prediction applied to their orbital pericenters. Halos within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other halos that cause significant mass loss -"galaxy harassment". Mergers between halos do not occur inside the clusters virial radius.
Monthly Notices of the Royal Astronomical Society, 2004
We use numerical simulations in a ΛCDM cosmology to model density profiles in a set of sixteen dark matter haloes with resolutions of up to seven million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10 11 to 10 14 M ⊙ , down to approximately 0.5% of the virial radius, and from redshift 5 to the present, covering a larger range in parameter space than previous works. We confirm that the cusp of the density profile is set at redshifts of two or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates.
Monthly Notices of The Royal Astronomical Society, 2006
Using six high-resolution dissipationless simulations with a varying box size in a flat Lambda cold dark matter (ΛCDM) universe, we study the mass and redshift dependence of dark matter halo shapes for Mvir= 9.0 × 1011− 2.0 × 1014 h−1 M⊙, over the redshift range z= 0–3, and for two values of σ8= 0.75 and 0.9. Remarkably, we find that the redshift, mass and σ8 dependence of the mean smallest-to-largest axis ratio of haloes is well described by the simple power-law relation 〈s〉= (0.54 ± 0.02)(Mvir/M*)−0.050±0.003, where s is measured at 0.3Rvir, and the z and σ8 dependences are governed by the characteristic non-linear mass, M*=M*(z, σ8). We find that the scatter about the mean s is well described by a Gaussian with σ∼ 0.1, for all masses and redshifts. We compare our results to a variety of previous works on halo shapes and find that reported differences between studies are primarily explained by differences in their methodologies. We address the evolutionary aspects of individual halo shapes by following the shapes of the haloes through ∼100 snapshots in time. We determine the formation scalefactor ac as defined by Wechsler et al. and find that it can be related to the halo shape at z= 0 and its evolution over time.
Astrophysical Journal, 2006
The formation and structure of dark matter (DM) halos is studied by means of constrained realizations of Gaussian fields using N-body simulations. A series of experiments of the formation of a 10^{12} Msun halo is designed to study the dependence of the density profile on its merging history. We confirm that the halo growth consists of violent and quiescent phases, with the density well approximated by the Navarro-Frenk-White (NFW) profile during the latter phases. We find that (1) the NFW scale radius R_s stays constant during the quiescent phase and grows abruptly during the violent one. In contrast, the virial radius grows linearly during the quiescent and abruptly during the violent phases. (2) The central density stays unchanged during the quiescent phase while dropping abruptly during the violent phase. (3) The value of \rs reflects the violent merging history of the halo, and depends on the number of violent events and their fractional magnitudes, independent of the time and order of these events. It does not reflect the formation time of the halo. (4) The fractional change in R_s is a nonlinear function of the fractional absorbed kinetic energy within R_s in a violent event.
1997
We use N -body simulations to investigate the structure and dynamical evolution of dark matter halos in clusters of galaxies. Our sample consists of nine massive halos from an Einstein-De Sitter universe with scale free power spectrum and spectral index n = -1. Halos are resolved by 20000 particles each, on average, and have a dynamical resolution of 20-25 kpc, as shown by extensive tests. Large scale tidal fields are included up to a scale L = 150 Mpc using background particles. We find that the halo formation process can be characterized by the alternation of two dynamical configurations: a merging phase and a relaxation phase, defined by their signature on the evolution of the total mass and root mean square (rms) velocity. Halos spend on average one third of their evolution in the merging phase and two thirds in the relaxation phase. Using this definition, we study the density profiles and show how they change during the halo dynamical history. In particular, we find that the average density profiles of our halos are fitted by the Navarro, Frenk & White (1995) analytical model with an rms residual of 17% between the virial radius R v and 0.01R v . The Hernquist (1990) analytical density profiles fits the same halos with an rms residual of 26%. The trend with mass of the scale radius of these fits is marginally consistent with that found by : compared to their results our halos are more centrally concentrated, and the relation between scale radius and halo mass is slightly steeper. We find a moderately large scatter in this relation, due both to dynamical evolution within halos and to fluctuations in the halo population. We analyze the dynamical equilibrium of our halos using the Jeans' equation, and find that on average they are approximately in equilibrium within their virial radius. Finally, we find that the projected mass profiles of our simulated halos are in very good agreement with the profiles of three rich galaxy clusters derived from strong and weak gravitational lensing observations.
Monthly Notices of the Royal Astronomical Society, 1997
We use N -body simulations to investigate the structure and dynamical evolution of dark matter halos in clusters of galaxies. Our sample consists of nine massive halos from an Einstein-De Sitter universe with scale free power spectrum and spectral index n = −1. Halos are resolved by 20000 particles each, on average, and have a dynamical resolution of 20-25 kpc, as shown by extensive tests. Large scale tidal fields are included up to a scale L = 150 Mpc using background particles. We find that the halo formation process can be characterized by the alternation of two dynamical configurations: a merging phase and a relaxation phase, defined by their signature on the evolution of the total mass and root mean square (rms) velocity. Halos spend on average one third of their evolution in the merging phase and two thirds in the relaxation phase. Using this definition, we study the density profiles and show how they change during the halo dynamical history. In particular, we find that the average density profiles of our halos are fitted by the Navarro, Frenk & White (1995) analytical model with an rms residual of 17% between the virial radius R v and 0.01R v . The Hernquist (1990) analytical density profiles fits the same halos with an rms residual of 26%. The trend with mass of the scale radius of these fits is marginally consistent with that found by : compared to their results our halos are more centrally concentrated, and the relation between scale radius and halo mass is slightly steeper. We find a moderately large scatter in this relation, due both to dynamical evolution within halos and to fluctuations in the halo population. We analyze the dynamical equilibrium of our halos using the Jeans' equation, and find that on average they are approximately in equilibrium within their virial radius. Finally, we find that the projected mass profiles of our simulated halos are in very good agreement with the profiles of three rich galaxy clusters derived from strong and weak gravitational lensing observations.
The Astrophysical Journal, 2007
The evolution of the phase-space density profile in dark matter (DM) halos is investigated by means of constrained simulations, designed to control the merging history of a given DM halo. Halos evolve through a series of quiescent phases of a slow accretion intermitted by violent events of major mergers. In the quiescent phases the density of the halo closely follows the NFW profile and the phase-space density profile, Q(r), is given by the Taylor & Navarro power law, r −β , where β ≈ 1.9 and stays remarkably stable over the Hubble time. Expressing the phase-space density by the NFW parameters, Q(r) = Q s (r/R s ) −β , the evolution of Q is determined by Q s . We have found that the effective mass surface density within R s , Σ s ≡ ρ s R s , remains constant throughout the evolution of a given DM halo along the main branch of its merging tree. This invariance entails that Q s ∝ R
Monthly Notices of the Royal Astronomical Society, 2014
We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity M ⋆ /M halo constrains the inner (γ) and outer (β) slopes of dark matter density, and the sharpness of transition between the slopes (α), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs.
The Astrophysical Journal, 2015
We present zoom-in N-body + Hydrodynamic simulations of dwarf central galaxies formed in Warm Dark Matter (WDM) halos with masses at present-day of 2 − 4 × 10 10 M ⊙. Two different cases are considered, the first one when halo masses are close to the corresponding half-mode filtering scale, M f (m WDM =1.2 keV) and the second when they are 20 to 30 times the corresponding M f (m WDM = 3.0 keV). The WDM simulations are compared with the respective Cold Dark Matter (CDM) simulations. The dwarfs formed in halos of masses (20 − 30)M f have roughly similar properties and evolution than their CDM counterparts; on the contrary, those formed in halos of masses around M f , are systematically different from their CDM counterparts. As compared to the CDM dwarfs, they assemble the dark and stellar masses later, having mass-weighted stellar ages 1.4-4.8 Gyr younger; their circular velocity profiles are shallower, with maximal velocities 20-60% lower; their stellar distributions are much less centrally concentrated and with larger effective radii, by factors 1.3-3. The WDM dwarfs at the filtering scale (m WDM =1.2 keV) have disk-like structures, and end in most cases with higher gas fractions and lower stellar-to-total mass ratios than their CDM counterparts. The late halo assembly, low halo concentrations, and the absence of satellites of the former with respect to the latter, are at the basis of the differences.
The Astrophysical Journal, 2003
We measure the average gravitational shear profile of 6 massive clusters (M vir ∼ 10 15 M ⊙ ) at z = 0.3 out to a radius ∼ 2h −1 Mpc. The measurements are fitted to a generalized NFW-like halo model ρ(r) with an arbitrary r → 0 slope α. The data are well fitted by such a model with a central cusp with α ∼ 0.9 − 1.6 (68% confidence interval). For the standard-NFW case α = 1.0, we find a concentration parameter c vir that is consistent with recent predictions from high-resolution CDM N-body simulations. Our data are also well fitted by an isothermal sphere model with a softened core. For this model, our 1σ upper limit for the core radius corresponds to a limit σ ⋆ ≤ 0.1cm 2 g −1 on the elastic collision cross-section in a self-interacting dark matter model.
We explore the density profiles of dark matter halos formed in high-resolution cosmological N-body simulations with power spectra of density fluctuations damped below a given length scale R f. The density profile shape of halos with masses smaller than the mass M f corresponding to R f does not significantly differ from that of the cold dark matter (CDM) halos; on the average, this shape is well described by the Navarro, Frenk, & White profile. However, the concentration of halos less massive than M f remains constant as the mass decreases, while for the CDM halos the concentration increases. The cosmogony of the former halos is not hierarchical but due to fragmentation of pancakes, their typical formation epoch being slightly later than that of halos of mass close to M f. The lower concentrations of these halos relative to the CDM ones may be explained by their late formation. If the suppression of power is associated with the free-streaming damping characteristic for warm dark matter (WDM) particles, then particles will have a residual thermal velocity dispersion v th. We have included this velocity dispersion in our simulations and have found that its influence on the halo structure is negligible. Using an analytical formalism, we show that halos formed by monolithic collapse will present a significant soft core only when particles have velocity dispersions much larger than the v th of warm particles. Relative to CDM, in a WDM model with R f ≈ 0.15 − 0.20 Mpc (flat cosmology with Ω Λ = h = 0.7) we find that not only the number of satellite galaxies in Milky-Way systems is in better agreement with the observations, but also the concentrations, the Tully-Fisher relation, and the formation epochs of dwarf galaxies.
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.