Universal profile of dark matter haloes and the spherical infall model

I propose a modification of the spherical infall model for the evolution of density fluctuations with initially Gaussian probability distribution and scale-free power spectra in the Einsteinde Sitter universe as developed by Hoffman & Shaham. I introduce a generalized form of the initial density distribution around an overdense region and cut it off at half the interpeak separation, accounting in this way for the presence of the neighbouring fluctuations. Contrary to the original predictions of Hoffman & Shaham, the resulting density profiles within virial radii no longer have a power-law shape, but their steepness increases with distance. The profiles of haloes of galactic mass are well fitted by the universal profile formula of changing slope obtained as a result of N -body simulations by Navarro, Frenk & White. The trend of steeper profiles for smaller masses and higher spectral indices is also reproduced. The agreement between the model and simulations is better for smaller masses and lower spectral indices, which suggests that galaxies form mainly by accretion, while formation of clusters involves merging.     

The real- and redshift-space density distribution functions for large-scale structure in the spherical collapse approximation

We use the spherical collapse (SC) approximation to derive expressions for the smoothed redshift-space probability distribution function (PDF), as well as the p -order hierarchical amplitudes S _p , in both real and redshift space. We compare our results with numerical simulations, focusing on the     standard CDM model, where redshift distortions are strongest. We find good agreement between the SC predictions and the numerical PDF in real space even for     , where σ _L is the linearly evolved rms fluctuation on the smoothing scale. In redshift space, reasonable agreement is possible only for     . Numerical simulations also yield a simple empirical relation between the real-space PDF and the redshift-space PDF: we find that for     , the redshift-space PDF, [ P δ _{( z )}], is, to a good approximation, a simple rescaling of the real-space PDF, P [ δ ], i.e.,     where σ and σ _{( z )} are the real-space and redshift-space rms fluctuations, respectively. This result applies well beyond the validity of linear perturbation theory, and it is a good fit for both the standard CDM model and the ΛCDM model. It breaks down for SCDM at     , but provides a good fit to the ΛCDM models for σ _L as large as 0.8.     

Old high-redshift galaxies and primordial density fluctuation spectra

We have discovered a population of extremely red galaxies at z  ≃ 1.5 which have apparent stellar ages of ≳ 3 Gyr, based on detailed spectroscopy in the rest-frame ultraviolet. In order for galaxies to have existed at the high collapse redshifts indicated by these ages, there must be a minimum level of power in the density fluctuation spectrum on galaxy scales. This paper compares the required power with that inferred from other high-redshift populations: damped Lyα absorbers and Lyman-limit galaxies at z  ≃ 3.2. If the collapse redshifts for the old red galaxies are in the range z _c ≃ 6–8, there is general agreement between the various tracers on the required inhomogeneity on 1-Mpc scales. This level of small-scale power requires the Lyman-limit galaxies to be approximately ν ≃ 3.0 fluctuations, implying a very large bias parameter b  ≃ 6. If the collapse redshifts of the red galaxies are indeed in the range z _c = 6–8 required for power spectrum consistency, their implied ages at z  ≃ 1.5 are between 3 and 3.8 Gyr for essentially any model universe of current age 14 Gyr. The age of these objects as deduced from gravitational collapse thus provides independent support for the ages estimated from their stellar populations. Such early-forming galaxies are rare, and their contribution to the cosmological stellar density is consistent with an extrapolation to higher redshifts of the star formation rate measured at z  < 5; there is no evidence for a general era of spheroid formation at extreme redshifts.     

The power spectrum of rich clusters of galaxies on large spatial scales

We present an analysis of the redshift-space power spectrum, P ( k ), of rich clusters of galaxies based on an automated cluster catalogue selected from the APM Galaxy Survey. We find that P ( k ) can be approximated by a power law, P ( k )∝ kⁿ , with n ≈−1.6 over the wavenumber range 0.04< k <0.1 h Mpc⁻¹. Over this range of wavenumbers, the APM cluster power spectrum has the same shape as the power spectra measured for optical and IRAS galaxies. This is consistent with a simple linear bias model in which different tracers have the same power spectrum as that of the mass distribution, but shifted in amplitude by a constant biasing factor. On larger scales, the power spectrum of APM clusters flattens and appears to turn over on a scale k ∼0.03 h Mpc⁻¹. We compare the power spectra estimated from simulated APM cluster catalogues with those estimated directly from cubical N -body simulation volumes, and find that the APM cluster survey should give reliable estimates of the true power spectrum at wavenumbers k ≳0.02 h Mpc⁻¹. These results suggest that the observed turnover in the power spectrum may be a real feature of the cluster distribution, and that we have detected the transition to a near-scale-invariant power spectrum implied by observations of anisotropies in the cosmic microwave background radiation. The scale of the turnover in the cluster power spectrum is in good agreement with the scale of the turnover observed in the power spectrum of APM galaxies.     

Power-spectrum normalization from the local abundance of rich clusters of galaxies

The number density of rich galaxy clusters still provides the most robust way of normalizing the power spectrum of dark matter perturbations on scales relevant to large-scale structure. We revisit this constraint in the light of several recent developments: (1) the availability of well-defined samples of local clusters with relatively accurate X-ray temperatures; (2) new theoretical mass functions for dark matter haloes, which provide a good fit to large numerical simulations; (3) more accurate mass–temperature relations from larger catalogues of hydrodynamical simulations; (4) the requirement to consider closed as well as open and flat cosmologies to obtain full multiparameter likelihood constraints for CMB and SNe studies. We present a new sample of clusters drawn from the literature and use this sample to obtain improved results on σ ₈, the normalization of the matter power spectrum on scales of 8  h ⁻¹ Mpc, as a function of the matter density and cosmological constant in a universe with general curvature. We discuss our differences with previous work, and the remaining major sources of uncertainty. Final results on the normalization, approximately independent of power spectrum shape, can be expressed as constraints on σ at an appropriate cluster normalization scale R _Cl. We provide fitting formulas for R _Cl and σ ( R _Cl) for general cosmologies, as well as for σ ₈ as a function of cosmology and shape parameter Γ. For flat models we find approximately σ ₈≃(0.495_−0.037^+0.034)Ω_M^−0.60 for Γ=0.23, where the error bar is dominated by uncertainty in the mass–temperature relation.     

The apparent and intrinsic shape of the APM galaxy clusters

We estimate the distribution of intrinsic shapes of APM galaxy clusters from the distribution of their apparent shapes. We measure the projected cluster ellipticities using two alternative methods. The first method is based on moments of the discrete galaxy distribution while the second is based on moments of the smoothed galaxy distribution. We study the performance of both methods using Monte Carlo cluster simulations covering the range of APM cluster distances and including a random distribution of background galaxies. We find that the first method suffers from severe systematic biases, whereas the second is more reliable. After excluding clusters dominated by substructure and quantifying the systematic biases in our estimated shape parameters, we recover a corrected distribution of projected ellipticities. We use the non-parametric kernel method to estimate the smooth apparent ellipticity distribution, and numerically invert a set of integral equations to recover the corresponding distribution of intrinsic ellipticities under the assumption that the clusters are either oblate or prolate spheroids. The prolate spheroidal model fits the APM cluster data best.     

The spherical infall model in a cosmological background density field

We discuss the influence of the cosmological background density field on the spherical infall model. The spherical infall model has been used in the PressSchechter formalism to evaluate the number abundance of clusters of galaxies, as well as to determine the density parameter of the Universe from the infalling flow. Therefore, the understanding of collapse dynamics plays a key role for extracting cosmological information. Here, we consider a modified version of the spherical infall model. We derive the mean field equations from the Newtonian fluid equations, in which the influence of cosmological background inhomogeneity is incorporated into the averaged quantities as the backreaction . By calculating the averaged quantities explicitly, we obtain simple expressions and find that, in the case of a scale-free power spectrum, density fluctuations with a negative spectral index make the infalling velocities slow. This suggests that we underestimate the density parameter when using the simple spherical infall model. In cases with the index n >0, the effect of background inhomogeneity could be negligible and the spherical infall model becomes a good approximation for infalling flows. We also present a realistic example with a cold dark matter power spectrum. In this case, the mean infall tends to be slow owing to the anisotropic random velocity.     

The alignment of clusters using large-scale simulations

The alignment of clusters of galaxies with their nearest neighbours and between clusters within a supercluster is investigated using simulations of 512³ dark matter particles for ΛCDM and τ CDM cosmological models. Strongly significant alignments are found for separations of up to 15  h ⁻¹ Mpc in both cosmologies, but for the ΛCDM model the alignments extend up to separations of 30  h ⁻¹ Mpc. The effect is strongest for nearest neighbours, but is not significant enough to be useful as an observational discriminant between cosmologies. As a check of whether this difference in alignments is present in other cosmologies, smaller simulations with 256³ particles are investigated for four different cosmological models. Because of poor number statistics, only the standard CDM model shows indications of having different alignments from the other models.     

Star formation in galaxies falling into clusters along supercluster-scale filaments

With the help of a statistical parameter derived from optical spectra, we show that the current star formation rate of a galaxy, falling into a cluster along a supercluster filament, is likely to undergo a sudden enhancement before the galaxy reaches the virial radius of the cluster. From a sample of 52 supercluster-scale filaments of galaxies joining a pair of rich clusters of galaxies within the two-degree Field Redshift Survey region, we find a significant enhancement of star formation, within a narrow range between ∼2 and  3  h ⁻¹₇₀ Mpc  of the centre of the cluster into which the galaxy is falling. This burst of star formation is almost exclusively seen in the fainter dwarf galaxies  ( M _B ≥−20)  . The relative position of the peak does not depend on whether the galaxy is a member of a group or not, but non-group galaxies have on average a higher rate of star formation immediately before falling into a cluster. From the various trends, we conclude that the predominant process responsible for this rapid burst is the close interaction with other galaxies falling into the cluster along the same filament, if the interaction occurs before the gas reservoir of the galaxy gets stripped off due to the interaction with the intracluster medium.     

Constraining the cosmological baryon density with X-ray clusters

We study the properties of X-ray galaxy clusters in four cold dark matter models with different baryon fractions Ω_BM, ranging from 5 to 20 per cent. By using an original three-dimensional hydrodynamic code based on the piecewise parabolic method, we run simulations on a box with a size of 64  h ⁻¹ Mpc and we identify the clusters by selecting the peaks in the X-ray luminosity field. We analyse these mock catalogues by computing the mass function, the luminosity function, the temperature distribution and the luminosity–temperature relation. By comparing the predictions of the different models to a series of recent observational results, we find that only the models with low baryonic content agree with the data, while models with larger baryon fraction are well outside the 1σ error bars. In particular, the analysis of the luminosity functions, both bolometric and in the energy band [0.5–2] keV, requires Ω_BM ≲ 0.05 when we fix the values h  = 0.5 and n  = 0.8 for the Hubble parameter and the primordial spectral index, respectively. Moreover we find that, independently of the cosmological scenario, all the considered quantities have very little redshift evolution, particularly between z  = 0.5 and 0.     

The impact of dust on the scaling properties of galaxy clusters

We investigate the effect of dust on the scaling properties of galaxy clusters based on hydrodynamic N -body simulations of structure formation. We have simulated five dust models plus radiative cooling and adiabatic models using the same initial conditions for all runs. The numerical implementation of dust was based on the analytical computations of Montier & Giard. We set up dust simulations to cover different combinations of dust parameters that make evident the effects of size and abundance of dust grains. Comparing our radiative plus dust cooling runs with a purely radiative cooling simulation, we find that dust has an impact on cluster scaling relations. It mainly affects the normalization of the scalings (and their evolution), whereas it introduces no significant differences in their slopes. The strength of the effect critically depends on the dust abundance and grain size parameters as well as on the cluster scaling. Indeed, cooling due to dust is effective in the cluster regime and has a stronger effect on the 'baryon driven' statistical properties of clusters such as   L _X– M , Y – M , S – M   scaling relations. Major differences, relative to the radiative cooling model, are as high as 25 per cent for the   L _X– M   normalization, and about 10 per cent for the Y – M and S – M normalizations at redshift zero. On the other hand, we find that dust has almost no impact on the 'dark matter driven'   T _mw– M   scaling relation. The effects are found to be dependent in equal parts on both dust abundances and grain size distributions for the scalings investigated in this paper. Higher dust abundances and smaller grain sizes cause larger departures from the radiative cooling (i.e. with no dust) model.     

Cosmological perturbation theory and the spherical collapse model — III. The velocity divergence field and the Ω dependence

Cosmological perturbation theory (PT) is a useful tool to study the cumulants of the density and velocity fields in the large-scale structure of the Universe. In Papers I and II of this series we saw that the spherical collapse (SC) model provides the exact solution to PT at tree-level and gives a good approximation to the loop corrections (next-to-leading orders), indicating negligible tidal effects. Here, we derive predictions for the (smoothed) cumulants of the velocity divergence field θ ≡ ▽ ⊙  v for an irrotational fluid in the SC model. By comparing these with the exact analytic results of Scoccimarro &38; Frieman, it is shown that, at least for the unsmoothed case, the loop corrections to the cumulants of θ are dominated by tidal effects. However, most of the tidal contribution seems to cancel out when computing the hierarchical ratios, T _J  = 〈θ ^J 〉 / 〈θ²〉  ^{J −1}. We also extend the work presented in Papers I and II to give predictions for the cumulants of the density and velocity divergence fields in non-flat spaces. In particular, we show the equivalence between the spherically symmetric solution to the equations of motion in the SC model (given in terms of the density) and that of the Lagrangian PT approach (given in terms of the displacement field). It is shown that the Ω dependence is very weak for both cosmic fields even at one loop (a 10 per cent effect at most), except for the overall factor f (Ω) that couples to the velocity divergence.     

Cosmological perturbation theory and the spherical collapse model — II. Non-Gaussian initial conditions

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.     

Evolution of peaks in weakly non-linear density field and dark halo profiles

Using the two-point Edgeworth series up to second order in the linear rms density fluctuation we construct the weakly non-linear conditional probability distribution function for the density field around an overdense region. This requires calculating the two-point analogues of the skewness parameter S ₃. We test the dependence of the two-point skewness on distance from the peak for scale-free power spectra and Gaussian smoothing. The statistical features of such a conditional distribution are given as the values obtained within linear theory corrected by the terms that arise as a result of weakly non-linear evolution. The expected density around the peak is found to be always below the linear prediction while its dispersion is always larger than in the linear case. For large enough overdensities the weakly non-linear corrections can be more significant than the peak constraint introduced by Bardeen et al. We apply these results to the spherical model of collapse as developed by Hoffman & Shaham and find that in general the effect of weakly non-linear interactions is to decrease the scale from which a peak gathers mass and therefore also the mass itself. In the case of an open universe this results in steepening of the final profile of the virialized proto-object.