Cosmological Perturbation Theory and the Spherical Collapse model — I. Gaussian initial conditions

We present a simple and intuitive approximation for solving the perturbation theory (PT) of small cosmic fluctuations. We consider only the spherically symmetric or monopole contribution to the PT integrals, which yields the exact result for tree-graphs (i.e. at leading order). We find that the non-linear evolution in Lagrangian space is then given by a simple local transformation over the initial conditions, although it is not local in Euler space. This transformation is found to be described by the spherical collapse (SC) dynamics, as it is the exact solution in the shearless (and therefore local) approximation in Lagrangian space. Taking advantage of this property, it is straightforward to derive the one-point cumulants, ξ_J, for both the unsmoothed and smoothed density fields to arbitrary order in the perturbative regime. To leading-order this reproduces, and provides us with a simple explanation for, the exact results obtained by Bernardeau. We then show that the SC model leads to accurate estimates for the next corrective terms when compared with the results derived in the exact perturbation theory making use of the loop calculations. The agreement is within a few per cent for the hierarchical ratios S _J  = ξ _J /ξ ^{J −1}₂. We compare our analytic results with N -body simulations, which turn out to be in very good agreement up to scales where σ ≈ 1. A similar treatment is presented to estimate higher order corrections in the Zel'dovich approximation. These results represent a powerful and readily usable tool to produce analytical predictions that describe the gravitational clustering of large-scale structure in the weakly non-linear regime.     

The three-point function in large-scale structure – I. The weakly non-linear regime in N-body simulations

This paper presents a comparison of the predictions for the two- and three-point correlation functions of density fluctuations, ξ and ζ , in gravitational perturbation theory (PT) against large cold dark matter (CDM) simulations. This comparison is made possible for the first time on large weakly non-linear scales (>10  h ⁻¹ Mpc) thanks to the development of a new algorithm for estimating correlation functions for millions of points in only a few minutes. Previous studies in the literature comparing the PT predictions of the three-point statistics with simulations have focused mostly on Fourier space, angular space or smoothed fields. Results in configuration space, such as those presented here, were limited to small scales where leading-order PT gives a poor approximation. Here we also propose and apply a method for separating the first-order and subsequent contributions to PT by combining different output times from the evolved simulations. We find that in all cases there is a regime where simulations do reproduce the leading-order (tree-level) predictions of PT for the reduced three-point function   Q ₃∼ ζ / ξ ²  . For steeply decreasing correlations (such as the standard CDM model) deviations from the tree-level results are important even at relatively large scales, ≃20 Mpc  h ⁻¹. On larger scales ξ goes to zero and the results are dominated by sampling errors. In more realistic models (such as the ΛCDM cosmology) deviations from the leading-order PT become important at smaller scales   r ≃10 Mpc  h ^-1  , although this depends on the particular three-point configuration. We characterize the range of validity of this agreement and show the behaviour of the next-order (one-loop) corrections.     

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.     

Generalized cumulant correlators and hierarchical clustering

The cumulant correlators, C _pq , are statistical quantities that generalize the better-known S _p parameters; the former are obtained from the two-point probability distribution function of the density fluctuations while the latter describe only the one-point distribution. If galaxy clustering develops from Gaussian initial fluctuations and a small-angle approximation is adopted, standard perturbative methods suggest a particular hierarchical relationship of the C _pq for projected clustering data, such as that obtained from the Automatic Plate Measuring (APM) survey. We establish the usefulness of the two-point cumulants for describing hierarchical clustering by comparing such calculations against available measurements from projected catalogues, finding very good agreement. We extend the idea of cumulant correlators to multipoint generalized cumulant correlators (related to the higher-order correlation functions). We extend previous studies in the highly non-linear regime to express the generalized cumulant correlators in terms of the underlying 'tree amplitudes' of hierarchical scaling models. Such considerations lead to a technique for determining these hierarchical amplitudes, to arbitrary order, from galaxy catalogues and numerical simulations. Knowledge of these amplitudes yields important clues about the phenomenology of gravitational clustering. For instance, we show that a three-point cumulant correlator can be used to separate the tree amplitudes up to sixth order. We also combine the particular hierarchical Ansatz of Bernardeau & Schaeffer with extended and hyper-extended perturbation theory to infer values of the tree amplitudes in the highly non-linear regime.     

The effects of anticorrelation on gravitational clustering

We use non-linear scaling relations (NSRs) to investigate the effects arising from the existence of negative correlations on the evolution of gravitational clustering in an expanding universe. It turns out that such anticorrelated regions have important dynamical effects on all scales. In particular, the mere existence of negative values for the linear two-point correlation function ξ¯ _L over some range of scales starting from l = L ₀ implies that the non-linear correlation function is bounded from above at all scales x < L ₀ . This also results in the relation ξ¯   ∝  x ⁻³ , at these scales, at late times, independent of the original form of the correlation function. Current observations do not rule out the existence of negative ξ¯ for 200  h ⁻¹ Mpc≲ ξ¯ ≲1000  h ⁻¹ Mpc; the present work may thus have relevance for the real Universe. The only assumption made in the analysis is the existence of NSR; the results are independent of the form of the NSR as well as of the stable clustering hypothesis.     

Clustering of luminous red galaxies – II. Small-scale redshift-space distortions

This is the second paper of a series where we study the clustering of luminous red galaxies (LRG) in the recent spectroscopic Sloan Digital Sky Survey (SDSS) data release, DR6, which has 75 000 LRG covering over  1 Gpc³  h ⁻³  for  0.15 < z < 0.47  . Here, we focus on modelling redshift-space distortions in  ξ(σ, π)  , the two-point correlation in separate line-of-sight and perpendicular directions, at small scales and in the line-of-sight. We show that a simple Kaiser model for the anisotropic two-point correlation function in redshift space, convolved with a distribution of random peculiar velocities with an exponential form, can describe well the correlation of LRG on all scales. We show that to describe with accuracy the so-called 'fingers-of-God' (FOG) elongations in the radial direction, it is necessary to model the scale dependence of both bias b and the pairwise rms peculiar velocity σ₁₂ with the distance. We show how both quantities can be inferred from the  ξ(σ, π)  data. From   r ≃ 10 Mpc  h ⁻¹  to   r ≃ 1 Mpc  h ⁻¹  , both the bias and σ₁₂ are shown to increase by a factor of 2: from   b = 2  to 4 and from  σ₁₂= 400  to  800 km s⁻¹  . The latter is in good agreement, within a 5 per cent accuracy in the recovered velocities, with direct velocity measurements in dark matter simulations with  Ω_m= 0.25  and  σ₈= 0.85  .     

On the origin of cold dark matter halo density profiles

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.     

Clustering of galaxies at 3.6 μm in the Spitzer Wide-area Infrared Extragalactic legacy survey

We investigate the clustering of galaxies selected in the 3.6 μm band of the Spitzer Wide-area Infrared Extragalactic (SWIRE) legacy survey. The angular two-point correlation function is calculated for 11 samples with flux limits of S _3.6≥ 4–400 μJy, over an 8 deg² field. The angular clustering strength is measured at >5σ significance at all flux limits, with amplitudes of A = (0.49–29) × 10⁻³ at 1°, for a power-law model, A θ^−0.8. We estimate the redshift distributions of the samples using phenomological models, simulations and photometric redshifts, and so derive the spatial correlation lengths. We compare our results with the Galaxies In Cosmological Simulations (GalICS) models of galaxy evolution and with parametrized models of clustering evolution. The GalICS simulations are consistent with our angular correlation functions, but fail to match the spatial clustering inferred from the phenomological models or the photometric redshifts. We find that the uncertainties in the redshift distributions of our samples dominate the statistical errors in our estimates of the spatial clustering. At low redshifts (median z ≤ 0.5), the comoving correlation length is approximately constant,   r ₀= 6.1 ± 0.5  h ⁻¹  Mpc, and then decreases with increasing redshift to a value of 2.9 ± 0.3  h ⁻¹ Mpc for the faintest sample, for which the median redshift is z ∼ 1. We suggest that this trend can be attributed to a decrease in the average galaxy and halo mass in the fainter flux-limited samples, corresponding to changes in the relative numbers of early- and late-type galaxies. However, we cannot rule out strong evolution of the correlation length over  0.5 < z < 1  .     

Power spectrum for the small-scale Universe

The first objects to arise in a cold dark matter (CDM) universe present a daunting challenge for models of structure formation. In the ultra small-scale limit, CDM structures form nearly simultaneously across a wide range of scales. Hierarchical clustering no longer provides a guiding principle for theoretical analyses and the computation time required to carry out credible simulations becomes prohibitively high. To gain insight into this problem, we perform high-resolution  ( N = 720³–1584³)  simulations of an Einstein–de Sitter cosmology where the initial power spectrum is   P ( k ) ∝ k ⁿ ,  with  −2.5 ≤ n ≤− 1  . Self-similar scaling is established for   n =−1  and −2 more convincingly than in previous, lower resolution simulations and for the first time, self-similar scaling is established for an   n =−2.25  simulation. However, finite box-size effects induce departures from self-similar scaling in our   n =−2.5  simulation. We compare our results with the predictions for the power spectrum from (one-loop) perturbation theory and demonstrate that the renormalization group approach suggested by McDonald improves perturbation theory's ability to predict the power spectrum in the quasi-linear regime. In the non-linear regime, our power spectra differ significantly from the widely used fitting formulae of Peacock & Dodds and Smith et al. and a new fitting formula is presented. Implications of our results for the stable clustering hypothesis versus halo model debate are discussed. Our power spectra are inconsistent with predictions of the stable clustering hypothesis in the high- k limit and lend credence to the halo model. Nevertheless, the fitting formula advocated in this paper is purely empirical and not derived from a specific formulation of the halo model.     

The Durham/UKST Galaxy Redshift Survey – III. Large-scale structure via the two-point correlation function

We have investigated the statistical clustering properties of galaxies by calculating the two-point galaxy correlation function from the optically selected Durham/UKST Galaxy Redshift Survey. This survey is magnitude-limited to b _J∼17, contains ∼2500 galaxies sampled at a rate of one-in-three and surveys a ∼4×10⁶ ( h ⁻¹ Mpc)³ volume of space. We have empirically determined the optimal method of estimating the two-point correlation function from just such a magnitude-limited survey. Applying our methods to this survey, we find that our redshift-space results agree well with those from previous optical surveys. In particular, we confirm the previously claimed detections of large-scale power out to ∼40 h ⁻¹ Mpc scales. We compare with two common models of cosmological structure formation and find that our two-point correlation function has power significantly in excess of the standard cold dark matter model in the 10–30 h ⁻¹ Mpc region. We therefore support the observational results of the APM galaxy survey. Given that only the redshift-space clustering can be measured directly, we use standard modelling methods and indirectly estimate the real-space two-point correlation function from the projected two-point correlation function. We then invert this projected correlation function to obtain an estimate of the spatial two-point correlation function in real space. This correlation function in real space has a lower amplitude than that in redshift space, but a steeper slope.     

Variance and skewness in the FIRST survey

We investigate the large-scale clustering of radio sources in the FIRST 1.4-GHz survey by analysing the distribution function ( counts in cells ). We select a reliable sample from the the FIRST catalogue, paying particular attention to the problem of how to define single radio sources from the multiple components listed. We also consider the incompleteness of the catalogue. We estimate the angular two-point correlation function w (θ), the variance Ψ₂ and skewness Ψ₃ of the distribution for the various subsamples chosen on different criteria. Both w (θ) and Ψ₂ show power-law behaviour with an amplitude corresponding to a spatial correlation length of r ₀ ∼ 10  h ⁻¹Mpc. We detect significant skewness in the distribution, the first such detection in radio surveys. This skewness is found to be related to the variance through Ψ₃ =  S ₃(Ψ₂)^α, with α = 1.9 ± 0.1, consistent with the non-linear gravitational growth of perturbations from primordial Gaussian initial conditions. We show that the amplitude of variance and the skewness are consistent with realistic models of galaxy clustering.     

The redshift-space two-point correlation function of ELAIS-S1 galaxies

We investigate the clustering properties of galaxies in the recently completed ELAIS-S1 redshift survey through their spatial two-point autocorrelation function. We used a subsample of the ELAIS-S1 catalogue covering approximately 4 deg² and consisting of 148 objects selected at 15 μm with a flux >0.5 mJy and a redshift   z < 0.5  . We detected a positive signal in the correlation function that in the range of separations  1–10  h ⁻¹ Mpc  is well approximated by a power law with a slope  γ= 1.4 ± 0.25  and a correlation length   s ₀= 5.4 ± 1.2  h ⁻¹ Mpc  , at the 90 per cent significance level. This result is in good agreement with the redshift-space correlation function measured in more local samples of mid-infrared-selected galaxies such as the IRAS Point Source Catalog (PSC z ) redshift survey. This suggests a lack of significant clustering evolution of infrared-selected objects out to   z = 0.5  that is further confirmed by the consistency found between the correlation functions measured in a local  ( z < 0.2)  and a distant  (0.2 < z < 0.5)  subsample of ELAIS-S1 galaxies. We also confirm that optically selected galaxies in the local redshift surveys, especially those of the SDSS sample, are significantly more clustered than infrared objects.     

A measurement of the faint source correlation function in the GOODS and UDF surveys

We present a stable procedure for defining and measuring the two point angular autocorrelation function,   w (θ) =[θ/θ₀( V )]^−Γ  , of faint  (25 < V < 29)  , barely resolved and unresolved sources in the Hubble Space Telescope Great Observatories Origins Deep Survey and Ultra Deep Field data sets. We construct catalogues that include close pairs and faint detections. We show, for the first time, that, on subarcsec scales, the correlation function exceeds unity. This correlation function is well fit by a power law with index  Γ≈ 2.5  and a  θ₀= 10^{−0.1( V −25.8)} arcsec  . This is very different from the values of  Γ≈ 0.7  and  θ₀( r ) = 10^{−0.4( r −21.5)} arcsec  associated with the gravitational clustering of brighter galaxies. This observed clustering probably reflects the presence of giant star-forming regions within galactic-scale potential wells. Its measurement enables a new approach to measuring the redshift distribution of the faintest sources in the sky.     

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.     

ARCS, the Arcminute Radio Cluster-lens Search – I. Selection criteria and initial results

We present the results of an unbiased radio search for gravitational lensing events with image separations between 15 and 60 arcsec, which would be associated with clusters of galaxies with masses >10^13–14 M_⊙. A parent population of 1023 extended radio sources stronger than 35 mJy with stellar optical identifications was selected using the FIRST radio catalogue at 1.4 GHz and the APM optical catalogue. The FIRST catalogue was then searched for companions to the parent sources stronger than 7 mJy and with separation in the range 15 to 60 arcsec. Higher-resolution observations of the resulting 38 lens candidates were made with the VLA at 1.4 and 5 GHz, and with MERLIN at 5 GHz in order to test the lens hypothesis in each case. None of our targets was found to be a gravitational lens system. These results provide the best current constraint on the lensing rate for this angular scale, but improved calculations of lensing rates from realistic simulations of the clustering of matter on the relevant scales are required before cosmologically significant constraints can be derived from this null result. We now have an efficient, tested observational strategy with which it will be possible to make an order-of-magnitude larger unbiased search in the near future.     

Irradiation pressure effects in close binary systems

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 ⁻¹ Mpc, denoted by σ ₈. When combined with the constraints on σ ₈ 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.     

CDM models with a BSI step-like primordial spectrum and a cosmological constant

A class of spatially flat models with cold dark matter (CDM), a cosmological constant and a broken-scale-invariant (BSI) step-like primordial (initial) spectrum of adiabatic perturbations, generated in an exactly solvable inflationary model where the inflaton potential has a rapid change of its first derivative at some point, is confronted with existing observational data on angular fluctuations of the CMB temperature, galaxy clustering and peculiar velocities of galaxies. If we locate the step in the initial spectrum at k  ≃ 0.05  h Mpc⁻¹, where a feature in the spectrum of Abell clusters of galaxies was found that could reflect a property of the initial spectrum, and if the large-scale flat plateau of the spectrum is normalized according to the COBE data, the only remaining parameter of the spectrum is p — the ratio of amplitudes of the metric perturbations between the small-scale and large-scale flat plateaux. Allowed regions in the plane of parameters (Ω = 1 − Ω_Λ,  H ₀) satisfying all data have been found for p lying in the region (0.8–1.7). Especially good agreement of the form of the present power spectrum in this model with the form of the cluster power spectrum is obtained for the inverted step ( p  < 1,  p  = 0.7–0.8), when the initial spectrum has slightly more power on small scales.     

Chandra measurements of the distribution of mass in the luminous lensing cluster Abell 2390

We present spatially resolved X-ray spectroscopy of the luminous lensing cluster Abell 2390, using observations made with the Chandra observatory. The temperature of the X-ray gas rises with increasing radius within the central ∼ 200 kpc of the cluster, and then remains approximately isothermal, with kT =11.5_−1.6^+1.5 keV , out to the limits of the observations at r ∼1.0 Mpc . The total mass profile determined from the Chandra data has a form in good agreement with the predictions from numerical simulations. Using the parametrization of Navarro, Frenk and White, we measure a scale radius r _s∼0.8 Mpc and a concentration parameter c ∼3 . The best-fitting X-ray mass model is in good agreement with independent gravitational lensing results and optical measurements of the galaxy velocity dispersion in the cluster. The X-ray gas to total mass ratio rises with increasing radius with f _gas∼21 per cent at r =0.9 Mpc . The azimuthally averaged 0.3–7.0 keV surface brightness profile exhibits a small core radius and a clear 'break' at r ∼500 kpc , where the slope changes from S _X ∼^∝  r ^−1.5 to S _X ∼^∝  r ^−3.6 . The data for the central region of the cluster indicate the presence of a cooling flow with a mass deposition rate of 200–300 M_⊙ yr⁻¹ and an effective age of 2–3 Gyr .