On the effects of resolution in dissipationless cosmological simulations     

Alexander Knebe  rey V. Kravtsov  †  Stefan Gottlöber  Anatoly A. Klypin 《Monthly notices of the Royal Astronomical Society》2000,317(3):630-648

We present a study of numerical effects in dissipationless cosmological simulations. The numerical effects are evaluated and studied by comparing the results of a series of 64³-particle simulations of varying force resolution and number of time-steps, performed using three of the N -body techniques currently used for cosmological simulations: the Particle–Mesh (PM), the Adaptive Particle–Particle–Particle–Mesh (AP³M) and the newer Adaptive Refinement Tree (ART) codes. This study can therefore be interesting both as an analysis of numerical effects and as a systematic comparison of different codes.
We find that the AP³M and the ART codes produce similar results given that convergence is reached within the code type. We also find that numerical effects may affect the high-resolution simulations in ways that have not been discussed before. In particular, our study revealed the presence of two-body scattering, the effects of which can be greatly amplified by inaccuracies in time integration. This process appears to affect the correlation function of matter, the mass function, the inner density of dark matter haloes and other statistics at scales much larger than the force resolution, although different statistics may be affected in a different fashion. We discuss the conditions for which strong two-body scattering is possible and discuss the choice of the force resolution and integration time-step. Furthermore, we discuss recent claims that simulations with force softening smaller than the mean interparticle separation are not trustworthy and argue that this claim is incorrect in general, and applies only to the phase-sensitive statistics. Our conclusion is that, depending on the choice of mass and force resolution and the integration time-step, a force resolution as small as 0.01 of the mean interparticle separation can be justified.  相似文献   

The distribution of ejected subhaloes and its implication for halo assembly bias     

Huiyuan Wang  H. J. Mo  Y. P. Jing 《Monthly notices of the Royal Astronomical Society》2009,396(4):2249-2256

Using a high-resolution cosmological N -body simulation, we identify the ejected population of subhaloes, which are haloes at redshift   z = 0  but were once contained in more massive 'host' haloes at high redshifts. The fraction of the ejected subhaloes in the total halo population of the same mass ranges from 9 to 4 per cent for halo masses from  ∼10¹¹  to  ∼10¹²  h ⁻¹ M_⊙  . Most of the ejected subhaloes are distributed within four times the virial radius of their hosts. These ejected subhaloes have distinct velocity distribution around their hosts in comparison to normal haloes. The number of subhaloes ejected from a host of given mass increases with the assembly redshift of the host. Ejected subhaloes in general reside in high-density regions, and have a much higher bias parameter than normal haloes of the same mass. They also have earlier assembly times, so that they contribute to the assembly bias of dark matter haloes seen in cosmological simulations. However, the assembly bias is not dominated by the ejected population, indicating that large-scale environmental effects on normal haloes are the main source for the assembly bias.  相似文献   

On the origin of cold dark matter halo density profiles     

Yu Lu  H. J. Mo  Neal Katz  Martin D. Weinberg 《Monthly notices of the Royal Astronomical Society》2006,368(4):1931-1940

N -body simulations predict that cold dark matter (CDM) halo-assembly occurs in two phases: (i) a fast-accretion phase with a rapidly deepening potential well; and (ii) a slow-accretion phase characterized by a gentle addition of mass to the outer halo with little change in the inner potential well. We demonstrate, using one-dimensional simulations, that this two-phase accretion leads to CDM haloes of the Navarro, Frenk & White (NFW) form and provides physical insight into the properties of the mass-accretion history that influence the final profile. Assuming that the velocities of CDM particles are effectively isotropized by fluctuations in the gravitational potential during the fast-accretion phase, we show that gravitational collapse in this phase leads to an inner profile  ρ( r ) ∝ r ⁻¹  . Slow accretion on to an established potential well leads to an outer profile with  ρ( r ) ∝ r ⁻³  . The concentration of a halo is determined by the fraction of mass that is accreted during the fast-accretion phase. Using an ensemble of realistic mass-accretion histories, we show that the model predictions of the dependence of halo concentration on halo formation time and, hence, the dependence of halo concentration on halo mass, and the distribution of halo concentrations all match those found in cosmological N -body simulations. Using a simple analytic model that captures much of the important physics, we show that the inner   r ⁻¹  profile of CDM haloes is a natural result of hierarchical mass assembly with an initial phase of rapid accretion.  相似文献   

The Ω dependence in equations of motion     

Adi Nusser  & Jörg M. Colberg 《Monthly notices of the Royal Astronomical Society》1998,294(3):457-464

The equations of motion governing the evolution of a collisionless gravitating system of particles in an expanding universe can be cast in a form which is almost independent of the cosmological density parameter, Ω, and the cosmological constant, Λ. The new equations are expressed in terms of a time variable τ=ln D , where D is the linear rate of growth of density fluctuations. The dependence on the density parameter is proportional to ε=Ω^−0.2−1 times the difference between the peculiar velocity (with respect to τ) of particles and the gravity field (minus the gradient of the potential); or, before shell-crossing, times the sum of the density contrast and the velocity divergence. In a one-dimensional collapse or expansion, the equations are fully independent of Ω and Λ before shell crossing. In the general case, the effect of this weak Ω dependence is to enhance the rate of evolution of density perturbations in dense regions. In a flat universe with Λ7ne;0, this enhancement is less pronounced than in an open universe with Λ=0 and the same Ω. Using the spherical collapse model, we find that the increase of the rms density fluctuations in a low-Ω universe relative to that in a flat universe with the same linear normalization is ∼0.01ε(Ω)〈δ³〉, where δ is the density field in the flat universe. The equations predict that the smooth average velocity field scales like Ω^0.6, while the local velocity dispersion (rms value) scales, approximately, like Ω^0.5. High-resolution N -body simulations confirm these results and show that density fields, when smoothed on scales slightly larger than clusters, are insensitive to the cosmological model. Haloes in an open model simulation are more concentrated than haloes of the same M /Ω in a flat model simulation.  相似文献   

Microlensing in dark matter haloes     

Janne Holopainen  Chris Flynn  Alexander Knebe  Stuart P. Gill  Brad K. Gibson 《Monthly notices of the Royal Astronomical Society》2006,368(3):1209-1222

Using eight dark matter haloes extracted from fully self-consistent cosmological N -body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of massive astronomical compact halo objects (MACHOs),   f = 1.0  , and a single MACHO mass is   m _M= 1.0 M_⊙  , the simulations yield mean values of  τ= 4.7^+5.0_−2.2× 10⁻⁷  and  Γ= 1.6^+1.3_−0.6× 10⁻⁶  events star⁻¹ yr⁻¹. We find that triaxiality and substructure can have major effects on the measured values so that τ and Γ values of up to three times the mean can be found. If we fit our values of τ and Γ to the MACHO collaboration observations, we find   f = 0.23^+0.15_−0.13  and   m _M= 0.44^+0.24_−0.16  . Five out of the eight haloes under investigation produce f and m _M values mainly concentrated within these bounds.  相似文献   

Resolving cosmic structure formation with the Millennium-II Simulation     

Michael Boylan-Kolchin  Volker Springel  Simon D. M. White  Adrian Jenkins  Gerard Lemson 《Monthly notices of the Royal Astronomical Society》2009,398(3):1150-1164

We present the Millennium-II Simulation (MS-II), a very large N -body simulation of dark matter evolution in the concordance Λ cold dark matter (ΛCDM) cosmology. The MS-II assumes the same cosmological parameters and uses the same particle number and output data structure as the original Millennium Simulation (MS), but was carried out in a periodic cube one-fifth the size  (100  h ⁻¹ Mpc)  with five times better spatial resolution (a Plummer equivalent softening of  1.0  h ⁻¹ kpc  ) and with 125 times better mass resolution (a particle mass of  6.9 × 10⁶  h ⁻¹ M_⊙  ). By comparing results at MS and MS-II resolution, we demonstrate excellent convergence in dark matter statistics such as the halo mass function, the subhalo abundance distribution, the mass dependence of halo formation times, the linear and non-linear autocorrelations and power spectra, and halo assembly bias. Together, the two simulations provide precise results for such statistics over an unprecedented range of scales, from haloes similar to those hosting Local Group dwarf spheroidal galaxies to haloes corresponding to the richest galaxy clusters. The 'Milky Way' haloes of the Aquarius Project were selected from a lower resolution version of the MS-II and were then resimulated at much higher resolution. As a result, they are present in the MS-II along with thousands of other similar mass haloes. A comparison of their assembly histories in the MS-II and in resimulations of 1000 times better resolution shows detailed agreement over a factor of 100 in mass growth. We publicly release halo catalogues and assembly trees for the MS-II in the same format within the same archive as those already released for the MS.  相似文献   

Testing tidal-torque theory – I. Spin amplitude and direction     

Cristiano Porciani  Avishai Dekel  Yehuda Hoffman 《Monthly notices of the Royal Astronomical Society》2002,332(2):325-338

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N -body simulation with thousands of well-resolved haloes. The protohaloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the protohaloes match, on average, the spin amplitudes of the virialized haloes of today, if linear growth is assumed until ∼ t ₀/3, or  55–70  per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ∼50° in the TTT prediction at t ₀, such that the linear spatial correlations of spins on scales ≥1  h ⁻¹ Mpc are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second-order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.  相似文献   

The formation and evolution of clusters of galaxies in different cosmogonies     

A. Huss  B. Jain  M. Steinmetz 《Monthly notices of the Royal Astronomical Society》1999,308(4):1011-1031

The formation of galaxy clusters in hierarchically clustering universes is investigated by means of high-resolution N -body simulations. The simulations are performed using a newly developed multimass scheme which combines a PM code with a high-resolution N -body code. Numerical effects resulting from time-stepping and gravitational softening are investigated, as well as the influence of the simulation box size and of the assumed boundary conditions. Special emphasis is laid on the formation process and the influence of various cosmological parameters. Cosmogonies with massive neutrinos are also considered. Differences between clusters in the same cosmological model seem to dominate over differences caused by differing background cosmogony. The cosmological model can alter the time evolution of cluster collapse, but the merging pattern remains fairly similar, e.g. the number of mergers and the mass ratio of mergers. The gross properties of a halo, such as its size and total angular momentum, also evolve in a similar manner for all cosmogonies, and can be described using analytical models. It is shown that the density distribution of a halo shows a characteristic radial dependence which follows a power law with a slope of =1 at small radii and =3 at large radii, independent of the background cosmogony or the considered redshift. The shape of the density profiles follows the generic form proposed by Navarro et al. for all hierarchically clustering scenarios, and retains very little information about the formation process or the cosmological model. Only the central matter concentration of a halo is correlated with the formation time and therefore the corresponding cosmogony. We emphasize the role of non-radial motions of the halo particles in the evolution of the density profile.  相似文献   

Investigating cosmological weak lensing with the ray-bundle method     

C. J.Fluke  R. L.Webster  D. J.Mortlock 《Monthly notices of the Royal Astronomical Society》2002,331(1):180-196

Using the ray-bundle method for calculating gravitational lens magnifications, we outline a method by which the magnification probability may be determined specifically in the weak lensing limit for cosmological models obtained from N -body simulations.
16 different models are investigated, which are variations on three broad classes of cold dark matter model: the standard model with  (Ω₀, λ ₀)=(1.0,0.0)  , the open model with  (Ω₀, λ ₀)=(0.3,0.0)  and the lambda model, which is a flat model with a cosmological constant  (Ω₀, λ ₀)=(0.3,0.7)  .
The effects of varying the Hubble parameter, H ₀, the power spectrum shape parameter, Γ, and the cluster mass normalization, σ ₈, are studied. It is shown that there is no signature of these parameters in the weak lensing magnification distributions. The magnification probability distributions are also shown to be independent of the numerical parameters such as the lens mass and simulation box size in the N -body simulations.  相似文献   

Modelling the redshift-space distortion of galaxy clustering     

Steve Hatton  & Shaun Cole 《Monthly notices of the Royal Astronomical Society》1998,296(1):10-20

  相似文献   

The Adaptive TreePM: an adaptive resolution code for cosmological N-body simulations     

J. S. Bagla  Nishikanta Khandai 《Monthly notices of the Royal Astronomical Society》2009,396(4):2211-2227

Cosmological N -body simulations are used for a variety of applications. Indeed progress in the study of large-scale structures and galaxy formation would have been very limited without this tool. For nearly 20 yr the limitations imposed by computing power forced simulators to ignore some of the basic requirements for modelling gravitational instability. One of the limitations of most cosmological codes has been the use of a force softening length that is much smaller than the typical interparticle separation. This leads to departures from collisionless evolution that is desired in these simulations. We propose a particle-based method with an adaptive resolution where the force softening length is reduced in high-density regions while ensuring that it remains well above the local interparticle separation. The method, called the Adaptive TreePM (ATreePM), is based on the TreePM code. We present the mathematical model and an implementation of this code, and demonstrate that the results converge over a range of options for parameters introduced in generalizing the code from the TreePM code. We explicitly demonstrate collisionless evolution in collapse of an oblique plane wave. We compare the code with the fixed resolution TreePM code and also an implementation that mimics adaptive mesh refinement methods and comment on the agreement and disagreements in the results. We find that in most respects the ATreePM code performs at least as well as the fixed resolution TreePM in highly overdense regions, from clustering and number density of haloes to internal dynamics of haloes. We also show that the adaptive code is faster than the corresponding high-resolution TreePM code.  相似文献   

Mass selection bias in galaxy cluster peculiar velocities from the kinetic Sunyaev–Zel'dovich effect     

Alan C. Peel 《Monthly notices of the Royal Astronomical Society》2006,365(4):1191-1202

Upcoming surveys for galaxy clusters using the Sunyaev–Zel'dovich effect are potentially sensitive enough to create a peculiar velocity catalogue. The statistics of these peculiar velocities are sensitive to cosmological parameters. We develop a method to explore parameter space using N -body simulations in order to quantify dark matter halo velocity statistics which will be useful for cluster peculiar velocity observations. We show that mass selection bias from a kinetic Sunyaev–Zel'dovich velocity catalogue forecasts rms peculiar velocities with a much more complicated  Ω_m  dependency than suggested by linear perturbation theory. In addition, we show that both two-point functions for velocities disagree with linear theory predictions out to  ∼40  h ⁻¹ Mpc  separations. A pedagogical appendix is included developing linear theory notation with respect to the two-point peculiar velocities functions.  相似文献   

Galactic halo cusp–core: tidal compression in mergers     

Avishai Dekel  Jonathan Devor  Guy Hetzroni 《Monthly notices of the Royal Astronomical Society》2003,341(1):326-342

We explain in simple terms how the build-up of dark haloes by merging compact satellites, as in the cold dark matter (CDM) cosmology, inevitably leads to an inner cusp of density profile  ρ∝ r ^−α  with  α≳ 1  , as seen in cosmological N -body simulations. A flatter halo core with  α < 1  exerts on the satellites tidal compression in all directions, which prevents the deposit of stripped satellite material in the core region. This makes the satellite orbits decay from the radius where  α∼ 1  to the halo centre with no local tidal mass transfer, and thus causes a rapid steepening of the inner profile to  α > 1  . These tidal effects, the resultant steepening of the profile to a cusp, and the stability of this cusp to tandem mergers with compact satellites are demonstrated using N -body simulations. The transition at  α∼ 1  is then addressed using toy models in the limiting cases of impulse and adiabatic approximations and using tidal radii for satellites on radial and circular orbits. In an associated paper, we address the subsequent slow convergence from either side to an asymptotic stable cusp with  α≳ 1  . Our analysis thus implies that an inner cusp is enforced when small haloes are typically more compact than larger haloes, as in the CDM scenario, such that enough satellite material makes it intact into the inner halo and is deposited there. We conclude that a necessary condition for maintaining a flat core, as indicated by observations, is that the inner regions of the CDM satellite haloes be puffed up by about 50 per cent such that when they merge into a larger halo they would be disrupted outside the halo core. This puffing up could be due to baryonic feedback processes in small haloes, which may be stimulated by the tidal compression in the halo cores.  相似文献   

Extended Press–Schechter theory and the density profiles of dark matter haloes     

Nicos Hiotelis † 《Monthly notices of the Royal Astronomical Society》2003,344(1):149-155

An inside–out model for the formation of haloes in a hierarchical clustering scenario is studied. The method combines the picture of the spherical infall model and a modification of the extended Press–Schechter theory. The mass accretion rate of a halo is defined to be the rate of its mass increase due to minor mergers. The accreted mass is deposited at the outer shells without changing the density profile of the halo inside its current virial radius. We applied the method to a flat Λ-cold dark matter universe. The resulting density profiles are compared with analytical models proposed in the literature, and a very good agreement is found. A trend is found of the inner density profile to become steeper for larger halo mass, which also results from recent N -body simulations. Additionally, present-day concentrations as well as their time evolution are derived and it is shown that they reproduce the results of large cosmological N -body simulations.  相似文献   

The universal rotation curve of spiral galaxies – II. The dark matter distribution out to the virial radius     

P. Salucci  A. Lapi  C. Tonini  G. Gentile  I. Yegorova  U. Klein 《Monthly notices of the Royal Astronomical Society》2007,378(1):41-47

In the current ΛCDM cosmological scenario, N -body simulations provide us with a universal mass profile, and consequently a universal equilibrium circular velocity of the virialized objects, as galaxies. In this paper we obtain, by combining kinematical data of their inner regions with global observational properties, the universal rotation curve of disc galaxies and the corresponding mass distribution out to their virial radius. This curve extends the results of Paper I, concerning the inner luminous regions of Sb–Im spirals, out to the edge of the galaxy haloes.  相似文献   

Constraining quasar host halo masses with the strength of nearby Lyα forest absorption     

Young-Rae Kim  Rupert A. C. Croft 《Monthly notices of the Royal Astronomical Society》2008,387(1):377-386

Using cosmological hydrodynamic simulations, we measure the mean transmitted flux in the Lyα forest for quasar sightlines that pass near a foreground quasar. We find that the trend of absorption with pixel quasar separation distance can be fitted using a simple power-law form including the usual correlation function parameters r ₀ and γ, so that     . From the simulations, we find the relation between r ₀ and quasar host mass, and formulate this as a way to estimate quasar host dark matter halo masses, quantifying uncertainties due to cosmological and IGM parameters, and redshift errors. With this method, we examine data for ∼9000 quasars from the Sloan Digital Sky Survey (SDSS) Data Release 5, assuming that the effect of ionizing radiation from quasars (the so-called transverse proximity effect) is unimportant (no evidence for it is seen in the data). We find that the best-fitting host halo mass for SDSS quasars with mean redshift z = 3 and absolute G -band magnitude −27.5 is  log  M /M_⊙= 12.68^+0.81_−0.67  . We also use the Lyman-Break Galaxy (LBG) and Lyα forest data of Adelberger et al. in a similar fashion to constrain the halo mass of LBGs to be  log₁₀  M /M_⊙= 11.41^+0.54_−0.59  , a factor of ∼20 lower than the bright quasars. In addition, we study the redshift distortions of the Lyα forest around quasars, using the simulations. We use the quadrupole to monopole ratio of the quasar Lyα forest correlation function as a measure of the squashing effect. We find its dependence on halo mass difficult to measure, but find that it may be useful for constraining cosmic geometry.  相似文献   

Dark matter halo profiles in scale-free cosmologies     

Steffen R. Knollmann  Chris Power  Alexander Knebe 《Monthly notices of the Royal Astronomical Society》2008,385(2):545-552

We explore the dependence of the central logarithmic slope of dark matter halo density profiles α on the spectral index n of the linear matter power spectrum P ( k ) using cosmological N -body simulations of scale-free models [i.e. P ( k ) ∝ k ⁿ ]. These simulations are based on a set of clear, reproducible and physically motivated criteria that fix the appropriate starting and stopping times for runs, and allow one to compare haloes across models with different spectral indices and mass resolutions. For each of our simulations we identify samples of well-resolved haloes in dynamical equilibrium and we analyse their mass profiles. By parametrizing the mass profile using a 'generalized' Navarro, Frenk & White profile in which the central logarithmic slope α is allowed to vary while preserving the r ⁻³ asymptotic form at large radii, we obtain preferred central slopes for haloes in each of our models. There is a strong correlation between α and n , such that α becomes shallower as n becomes steeper. However, if we normalize our mass profiles by r ₋₂, the radius at which the logarithmic slope of the density profile is −2, we find that these differences are no longer present. This is apparent if we plot the maximum slope     as a function of r / r ₋₂– we find that the profiles are similar for haloes forming in different n models. This reflects the importance of concentration, and reveals that the concentrations of haloes forming in steep- n cosmologies tend to be smaller than those of haloes forming in shallow- n cosmologies. We conclude that there is no evidence for convergence to a unique central asymptotic slope, at least on the scales that we can resolve.  相似文献   

Evolution of the mass function of dark matter haloes     

Darren Reed  Jeffrey Gardner  Thomas Quinn  Joachim Stadel  Mark Fardal  George Lake  Fabio Governato 《Monthly notices of the Royal Astronomical Society》2003,346(2):565-572

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 15, in a  50 h ⁻¹ 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 & Tormen mass function provides an excellent match to all of our data except for redshifts of 10 and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 per cent for the  10¹⁰–10¹¹ 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 & 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 per cent. We test the sensitivity of our results to force accuracy, starting redshift and halo-finding algorithm.  相似文献   

Mass-to-light ratio gradients in early-type galaxy haloes   总被引：1，自引：0，他引：1  

N. R. Napolitano  M. Capaccioli  A. J. Romanowsky  N. G. Douglas  M. R. Merrifield  K. Kuijken  M. Arnaboldi  O. Gerhard  K. C. Freeman 《Monthly notices of the Royal Astronomical Society》2005,357(2):691-706

Owing to the fact that the near future should see a rapidly expanding set of probes of the halo masses of individual early-type galaxies, we introduce a convenient parameter for characterizing the halo masses from both observational and theoretical results:  ∇_ℓϒ  , the logarithmic radial gradient of the mass-to-light ratio. Using halo density profiles from Λ-cold dark matter (CDM) simulations, we derive predictions for this gradient for various galaxy luminosities and star formation efficiencies  ε_SF  . As a pilot study, we assemble the available  ∇_ℓϒ  data from kinematics in early-type galaxies – representing the first unbiased study of halo masses in a wide range of early-type galaxy luminosities – and find a correlation between luminosity and  ∇_ℓϒ  , such that the brightest galaxies appear the most dark-matter dominated. We find that the gradients in most of the brightest galaxies may fit in well with the ΛCDM predictions, but that there is also a population of fainter galaxies whose gradients are so low as to imply an unreasonably high star formation efficiency  ε_SF > 1  . This difficulty is eased if dark haloes are not assumed to have the standard ΛCDM profiles, but lower central concentrations.  相似文献   

Large-scale non-Gaussian mass function and halo bias: tests on N-body simulations     

M. Grossi  L. Verde  C. Carbone  K. Dolag  E. Branchini  F. Iannuzzi  S. Matarrese  L. Moscardini 《Monthly notices of the Royal Astronomical Society》2009,398(1):321-332

  相似文献