首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 64 毫秒
1.
We study the distribution function (DF) of dark matter particles in haloes of mass range  1014–1015 M⊙  . In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Λ cold dark matter (ΛCDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L , a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form     . This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (β0 and  β∞  ) and the scale of transition between them ( L 0). The energy part   f E ( E )  is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density   Q =ρ( r )/σ3( r )  .  相似文献   

2.
  相似文献   

3.
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  (Ω0, λ 0)=(1.0,0.0)  , the open model with  (Ω0, λ 0)=(0.3,0.0)  and the lambda model, which is a flat model with a cosmological constant  (Ω0, λ 0)=(0.3,0.7)  .
The effects of varying the Hubble parameter, H 0, the power spectrum shape parameter, Γ, and the cluster mass normalization, σ 8, 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.  相似文献   

4.
This is the third paper in a series which combines N -body simulations and semi-analytic modelling to provide a fully spatially resolved simulation of the galaxy formation and clustering processes. Here we extract mock redshift surveys from our simulations: a cold dark matter model with either Ω0=1 ( τ CDM) or Ω0=0.3 and Λ=0.7 (ΛCDM). We compare the mock catalogues with the northern region (CfA2N) of the Center for Astrophysics (CfA) Redshift Surveys. We study the properties of galaxy groups and clusters identified using standard observational techniques, and also the relation of these groups to real virialized systems. Most features of CfA2N groups are reproduced quite well by both models with no obvious dependence on Ω0. Redshift‐space correlations and pairwise velocities are also similar in the two cosmologies. The luminosity functions predicted by our galaxy formation models depend sensitively on the treatment of star formation and feedback. For the particular choices of Paper I they agree poorly with the CfA survey. To isolate the effect of this discrepancy on our mock redshift surveys, we modify galaxy luminosities in our simulations to reproduce the CfA luminosity function exactly. This adjustment improves agreement with the observed abundance of groups, which depends primarily on the galaxy luminosity density, but other statistics, connected more closely with the underlying mass distribution, remain unaffected. Regardless of the luminosity function adopted, modest differences with observation remain. These can be attributed to the presence of the 'Great Wall' in the CfA2N. It is unclear whether the greater coherence of the real structure is a result of cosmic variance, given the relatively small region studied, or reflects a physical deficiency of the models.  相似文献   

5.
We generate mock galaxy catalogues for a grid of different cosmologies, using rescaled N -body simulations in tandem with a semi-analytic model run using consistent parameters. Because we predict the galaxy bias, rather than fitting it as a nuisance parameter, we obtain an almost pure constraint on σ8 by comparing the projected two-point correlation function we obtain to that from the Sloan Digital Sky Survey (SDSS). A systematic error arises because different semi-analytic modelling assumptions allow us to fit the r -band luminosity function equally well. Combining our estimate of the error from this source with the statistical error, we find  σ8= 0.97 ± 0.06  . We obtain consistent results if we use galaxy samples with a different magnitude threshold, or if we select galaxies by b J-band rather than r -band luminosity and compare to data from the 2dF Galaxy Redshift Survey (2dFGRS). Our estimate for σ8 is higher than that obtained for other analyses of galaxy data alone, and we attempt to find the source of this difference. We note that in any case, galaxy clustering data provide a very stringent constraint on galaxy formation models.  相似文献   

6.
This paper presents a fast, economical particle-multiple-mesh N -body code optimized for large- N modelling of collisionless dynamical processes, such as black hole wandering or bar–halo interactions, occurring within isolated galaxies. The code has been specially designed to conserve linear momentum. Despite this, it also has variable softening and an efficient block-time-step scheme: the force between any pair of particles is calculated using the finest mesh that encloses them both (respecting Newton's third law) and is updated only on the longest time-step of the two (which conserves momentum). For realistic galaxy models with   N ≳ 106  , it is faster than the fastest comparable tree code by factors ranging from ∼2 (using single time-steps) to ∼10 (multiple time-steps in a concentrated galaxy).  相似文献   

7.
We study the effect of primordial non-Gaussianity on the development of large-scale cosmic structure using high-resolution N -body simulations. In particular, we focus on the topological properties of the 'cosmic web', quantitatively characterized by the Minkowski functionals (MFs), for models with quadratic non-linearities with different values of the usual non-Gaussianity parameter f NL. In the weakly non-linear regime (the amplitude of mass density fluctuations σ0 < 0.1), we find that analytic formulae derived from perturbation theory agree with the numerical results within a few per cent of the amplitude of each MF when | f NL| < 1000. In the non-linear regime, the detailed behaviour of the MFs as functions of threshold density deviates more strongly from the analytical curves, while the overall amplitude of the primordial non-Gaussian effect remains comparable to the perturbative prediction. When smaller-scale information is included, the influence of primordial non-Gaussianity becomes increasingly significant statistically due to decreasing sample variance. We find that the effect of the primordial non-Gaussianity with  | f NL| = 50  is comparable to the sample variance of mass density fields with a volume of 0.125( h −1 Gpc)3 when they are smoothed by Gaussian filter at a scale of 5  h −1 Mpc. The detectability of this effect in actual galaxy surveys will strongly depend on residual uncertainties in cosmological parameters and galaxy biasing.  相似文献   

8.
We present a simple and efficient method to set up spherical structure models for N -body simulations with a multimass technique. This technique reduces by a substantial factor the computer run time needed in order to resolve a given scale as compared to single-mass models. It therefore allows to resolve smaller scales in N -body simulations for a given computer run time. Here, we present several models with an effective resolution of up to  1.68 × 109  particles within their virial radius which are stable over cosmologically relevant time-scales. As an application, we confirm the theoretical prediction by Dehnen that in mergers of collisionless structures like dark matter haloes always the cusp of the steepest progenitor is preserved. We model each merger progenitor with an effective number of particles of approximately 108 particles. We also find that in a core–core merger the central density approximately doubles whereas in the cusp–cusp case the central density only increases by approximately 50 per cent. This may suggest that the central regions of flat structures are better protected and get less energy input through the merger process.  相似文献   

9.
Using N -body simulations with a large set of massless test particles, we compare the predictions of two theories of violent relaxation, the well-known Lynden-Bell theory and the more recent theory by Nakamura. We derive 'weakened' versions of both the theories in which we use the whole equilibrium coarse-grained distribution function     as a constraint instead of the total energy constraint. We use these weakened theories to construct expressions for the conditional probability   Ki (τ)  that a test particle initially at the phase-space coordinate τ would end-up in the i th macro-cell at equilibrium. We show that the logarithm of the ratio   Rij (τ) ≡ Ki (τ)/ Kj (τ)  is directly proportional to the initial phase-space density   f 0(τ)  for the Lynden-Bell theory and inversely proportional to   f 0(τ)  for the Nakamura theory. We then measure   Rij (τ)  using a set of N -body simulations of a system undergoing a gravitational collapse to check the validity of the two theories of violent relaxation. We find that both the theories are at odds with the numerical results, both qualitatively and quantitatively.  相似文献   

10.
A revision of Stodółkiewicz's Monte Carlo code is used to simulate evolution of large star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. A survey of the evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods (Fokker–Planck, Monte Carlo and N -body). The initial rapid mass loss, resulting from stellar evolution of the most massive stars, causes expansion of the whole cluster and eventually leads to the disruption of less bound systems ( W 0=3). Models with larger W 0 survive this phase of evolution and then undergo core collapse and subsequent post-collapse expansion, like isolated models. The expansion phase is eventually reversed when tidal limitation becomes important. The results presented are the first major step in the direction of simulating evolution of real globular clusters by means of the Monte Carlo method.  相似文献   

11.
In the presence of a strong   m = 2  component in a rotating galaxy, the phase-space structure near corotation is shaped to a large extent by the invariant manifolds of the short-period family of unstable periodic orbits terminating at L 1 or L 2. The main effect of these manifolds is to create robust phase correlations among a number of chaotic orbits large enough to support a spiral density wave outside corotation. The phenomenon is described theoretically by soliton-like solutions of a Sine–Gordon equation. Numerical examples are given in an N -body simulation of a barred spiral galaxy. In these examples, we demonstrate how the projection of unstable manifolds in configuration space reproduces essentially the entire observed bar–spiral pattern.  相似文献   

12.
We analyse a   z < 0.1  galaxy sample from the Sloan Digital Sky Survey focusing on the variation in the galaxy colour bimodality with stellar mass     and projected neighbour density Σ, and on measurements of the galaxy stellar mass functions. The characteristic mass increases with environmental density from about  1010.6  to     (Kroupa initial mass function,   H 0= 70  ) for Σ in the range  0.1–10 Mpc−2  . The galaxy population naturally divides into a red and blue sequence with the locus of the sequences in colour–mass and colour–concentration indices not varying strongly with environment. The fraction of galaxies on the red sequence is determined in bins of 0.2 in  log Σ  and     bins). The red fraction   f r   generally increases continuously in both Σ and     such that there is a unified relation:     . Two simple functions are proposed which provide good fits to the data. These data are compared with analogous quantities in semi-analytical models based on the Millennium N -body simulation: the Bower et al. and Croton et al. models that incorporate active galactic nucleus feedback. Both models predict a strong dependence of the red fraction on stellar mass and environment that is qualitatively similar to the observations. However, a quantitative comparison shows that the Bower et al. model is a significantly better match; this appears to be due to the different treatment of feedback in central galaxies.  相似文献   

13.
We use three-integral models to infer the distribution function (DF) of the boxy E3–E4 galaxy NGC 1600 from surface brightness and line-profile data on the minor and major axes. We assume axisymmetry and that the mass-to-light ratio is constant in the central ∼1 R e. Stars in the resulting gravitational potential move mainly on regular orbits. We use an approximate third integral K from perturbation theory and write the DF as a sum of basis functions in the three integrals E , L z and K . We then fit the projected moments of these basis functions to the kinematic observables and deprojected density, using a non-parametric algorithm. The deduced dynamical structure is radially anisotropic, with σ θ σ r ≈ σ φ σ r ≈0.7 on the major axis. Both on the minor axis and near the centre the velocity distribution is more isotropic; thus the model is flattened by equatorial radial orbits. The kinematic data are fitted without the need for a central black hole; the central mass determined previously from ground-based data therefore overestimates the actual black-hole mass. The mass-to-light ratio of the stars is M L V =6  h 50. The anisotropy structure of NGC 1600 with a radially anisotropic main body and more nearly isotropic centre is similar to that found recently in NGC 1399, 2434, 3379 and 6703, suggesting that this pattern may be common amongst massive elliptical galaxies. We discuss a possible merger origin of NGC 1600 in the light of these results.  相似文献   

14.
We estimate the two- and three-dimensional power spectra, P 2( K ) and P 3( k ), of the galaxy distribution by applying a maximum likelihood estimator to pixel maps of the APM Galaxy Survey. The analysis provides optimal estimates of the power spectra and of their covariance matrices if the fluctuations are assumed to be Gaussian. Our estimates of P 2( K ) and P 3( k ) are in good agreement with previous work, but we find that the errors at low wavenumbers have been underestimated in some earlier studies. If the galaxy power spectrum is assumed to have the same shape as the mass power spectrum, then the APM maximum likelihood P 3( k ) estimates at k ≤0.19  h  Mpc−1 constrain the amplitude and shape parameter of a scale-invariant CDM model to lie within the 2 σ ranges 0.74≤( σ 8)g≤1.28 and 0.06≤Γ≤0.46 . Using the Galactic extinction estimates of Schlegel, Finkbeiner & Davis, we show that Galactic obscuration has a negligible effect on galaxy clustering over most of the area of the APM Galaxy Survey.  相似文献   

15.
The stochasticity in the distribution of dark haloes in the cosmic density field is reflected in the distribution function   P V ( N h| δ m)  , which gives the probability of finding N h haloes in a volume V with mass density contrast δ m. We study the properties of this function using high-resolution N -body simulations, and find that   P V ( N h| δ m)  is significantly non-Poisson. The ratio between the variance and the mean goes from ∼1 (Poisson) at  1+ δ m≪1  to <1 (sub-Poisson) at  1+ δ m∼1  to >1 (super-Poisson) at  1+ δ m≫1  . The mean bias relation is found to be well described by halo bias models based on the Press–Schechter formalism. The sub-Poisson variance can be explained as a result of halo exclusion, while the super-Poisson variance at high δ m may be explained as a result of halo clustering. A simple phenomenological model is proposed to describe the behaviour of the variance as a function of δ m. Galaxy distribution in the cosmic density field predicted by semi-analytic models of galaxy formation shows similar stochastic behaviour. We discuss the implications of the stochasticity in halo bias to the modelling of higher order moments of dark haloes and of galaxies.  相似文献   

16.
We present measurements of the clustering properties of galaxies in the field of redshift range 0.5 ≲ z ≲ 1.5 Ultra Steep Spectrum radio sources selected from the Sydney University Molonglo Sky Survey and the National Radio Astronomy Observatories Very Large Array Sky Survey. Galaxies in these USS fields were identified in deep near-infrared observations, complete down to   K s= 20  , using the IRIS2 instrument at the Anglo-Australian Telescope. We used the redshift distribution of   K s < 20  galaxies taken from Cimatti et al. (2002) to constrain the correlation length r 0. We find a strong correlation signal of galaxies with   K s < 20  around our USS sample. A comoving correlation length   r 0= 14.0 ± 2.8  h −1 Mpc  and γ= 1.98 ± 0.15 are derived in a flat cosmological model universe.
We compare our findings with those obtained in a cosmological N -body simulation populated with galform semi-analytic galaxies. We find that clusters of galaxies with masses in the range   M = 1013.4–14.2  h −1 M⊙  have a cluster–galaxy cross-correlation amplitude comparable to those found between the USS hosts and galaxies. These results suggest that distant radio galaxies are excellent tracers of galaxy overdensities and pinpoint the progenitors of present day rich clusters of galaxies.  相似文献   

17.
We test an analytic model for the two-point correlations of galaxy clusters in redshift space using the Hubble volume N -body simulations. The correlation function of clusters shows no enhancement along the line of sight, owing to the lack of any virialized structures in the cluster distribution. However, the distortion of the clustering pattern arising from coherent bulk motions is clearly visible. The distribution of cluster peculiar motions is well described by a Gaussian, except in the extreme high-velocity tails. The simulations produce a small but significant number of clusters with large peculiar motions. The form of the redshift-space power spectrum is strongly influenced by errors in measured cluster redshifts in extant surveys. When these errors are taken into account, the model reproduces the power spectrum recovered from the simulation to an accuracy of 15 per cent or better over a decade in wavenumber. We compare our analytic predictions with the power spectrum measured from the APM cluster redshift survey. The cluster power spectrum constrains the amplitude of density fluctuations, as measured by the linear rms variance in spheres of radius 8  h −1 Mpc, denoted by σ 8. When combined with the constraints on σ 8 and the density parameter Ω derived from the local abundance of clusters, we find a best-fitting cold dark matter model with     and     , for a power spectrum shape that matches that measured for galaxies. However, for the best-fitting value of Ω and given the value of Hubble's constant from recent measurements, the assumed shape of the power spectrum is incompatible with the most readily motivated predictions from the cold dark matter paradigm.  相似文献   

18.
Possible orbital histories of the Sgr dwarf galaxy are explored. A special-purpose N -body code is used to construct the first models of the Milky Way–Sgr dwarf system in which both the Milky Way and the Sgr dwarf are represented by full N -body systems and followed for a Hubble time. These models are used to calibrate a semi-analytic model of the Sgr dwarf's orbit that enables us to explore a wider parameter space than is accessible to the N -body models. We conclude that the extant data on the Sgr dwarf are compatible with a wide range of orbital histories. At one extreme the Sgr dwarf initially possesses ∼1011 M⊙ and starts from a Galactocentric distance R D(0)≳200 kpc. At the other extreme the Sgr dwarf starts with ∼109 M⊙ and R D(0)∼60 kpc, similar to its present apocentric distance. In all cases the Sgr dwarf is initially dark matter dominated and the current velocity dispersion of the Sgr dwarf's dark matter is tightly constrained to be 21±2 km s−1. This number is probably compatible with the smaller measured dispersion of the Sgr dwarf's stars because of (i) the dynamical difference between dark and luminous matter, and (ii) velocity anisotropy.  相似文献   

19.
correlator of the galaxy density field Q 21 is examined from the point of view of biasing. It is shown that, to leading order, it depends on two biasing parameters b b 2, and on q 21, the underlying cumulant correlator of the mass. As the skewness Q 3 has analogous properties, the slope of the correlation function −γ, Q 3 and Q 21 uniquely determine the bias parameter on a particular scale to be b  = γ/6( Q 21 −  Q 3), when working in the context of gravitational instability with Gaussian initial conditions. Thus on large scales, easily accessible with the future Sloan Digital Sky Survey and the 2 Degree Field Survey, it will be possible to extract b b 2 from simple counts-in-cells measurements. Moreover, the higher order cumulants, Q N , successively determine the higher order biasing parameters. From these it is possible to predict higher order cumulant correlators as well. Comparison of the predictions with the measurements will provide internal consistency checks on the validity of the assumptions in the theory, most notably perturbation theory of the growth of fluctuations by gravity and Gaussian initial conditions. Since the method is insensitive Ω, it can be successfully combined with results from velocity fields, which determine Ω0.6/b, to measure the total density parameter in the Universe.  相似文献   

20.
In Paper I of this series, we introduced the spherical collapse (SC) approximation in Lagrangian space as a way of estimating the cumulants ξ J of density fluctuations in cosmological perturbation theory (PT). Within this approximation, the dynamics is decoupled from the statistics of the initial conditions, so we are able to present here the cumulants for generic non-Gaussian initial conditions, which can be estimated to arbitrary order including the smoothing effects. The SC model turns out to recover the exact leading-order non-linear contributions up to terms involving non-local integrals of the J -point functions. We argue that for the hierarchical ratios S J , these non-local terms are subdominant and tend to compensate each other. The resulting predictions show a non-trivial time evolution that can be used to discriminate between models of structure formation. We compare these analytic results with non-Gaussian N -body simulations, which turn out to be in very good agreement up to scales where σ ≲ 1.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号