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  .     

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.     

The redshift evolution of clustering in the Hubble Deep Field

We present a correlation function analysis for the catalogue of photometric redshifts obtained from the Hubble Deep Field image by Fernandez-Soto, Lanzetta & Yahil. By dividing the catalogue into redshift bins of width Δ z =0.4 we measured the angular correlation function w ( θ ) as a function of redshift up to z ∼4.8. From these measurements we derive the trend of the correlation length r ₀. We find that r ₀( z ) is roughly constant with look-back time up to z ≃2, and then increases to higher values at z ≳2.4. We estimate the values of r ₀, assuming ξ ( r , z )=[ r r ₀( z )]^{− γ} , γ =1.8 and various geometries. For Ω₀=1 we find r ₀( z =3)≃7.00±4.87  h ⁻¹ Mpc, in good agreement with the values obtained from analysis of the Lyman break galaxies.     

Peculiar velocities and the mean density parameter

We study the peculiar velocity field inferred from the Mark III spirals using a new method of analysis. We estimate optimal values of Tully–Fisher scatter and zero-point offset, and we derive the three-dimensional rms peculiar velocity ( σ _v ) of the galaxies in the samples analysed. We check our statistical analysis using mock catalogues derived from numerical simulations of cold dark matter (CDM) models considering measurement uncertainties and sampling variations. Our best determination for the observations is σ _v =(660±50) km s⁻¹. We use the linear theory relation between σ _v , the density parameter Ω, and the galaxy correlation function ξ ( r ) to infer the quantity     , where b is the linear bias parameter of optical galaxies and the uncertainties correspond to bootstrap resampling and an estimated cosmic variance added in quadrature. Our findings are consistent with the results of cluster abundances and redshift-space distortion of the two-point correlation function. These statistical measurements suggest a low value of the density parameter Ω∼0.4 if optical galaxies are not strongly biased tracers of mass.     

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.     

Power spectrum shape from peculiar velocity data

We constrain the velocity power spectrum shape parameter Γ in linear theory using the nine bulk flow and shear moments estimated from four recent peculiar velocity surveys. For each survey, a likelihood function for Γ was found after marginalizing over the power spectrum amplitude  σ₈Ω^0.6_m  using constraints obtained from comparisons between redshift surveys and peculiar velocity data. In order to maximize the accuracy of our analyses, the velocity noise σ_* was estimated directly for each survey. A statistical analysis of the differences between the values of the moments estimated from different surveys showed consistency with theoretical predictions, suggesting that all the surveys investigated reflect the same large-scale flows. The peculiar velocity surveys were combined into a composite survey yielding the constraint  Γ= 0.13^+0.09_−0.05  . This value is lower than, but consistent with, values obtained using redshift surveys and cosmic microwave background data.     

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 model of slingshot prominences in rapidly rotating stars

We analyse the spatial clustering properties of the ROSAT All-Sky Survey (RASS) 1 Bright Sample, an X-ray flux-limited catalogue of galaxy clusters selected from the southern part of the survey. The two-point correlation function ( r ) of the whole sample is well fitted (in an Einsteinde Sitter model) by the power law =( r r ₀) , with and (95.4 per cent confidence level with one fitting parameter). We use the RASS 1 Bright Sample as a first application of a theoretical model that aims to predict the clustering properties of X-ray clusters in flux-limited surveys for different cosmological scenarios. The model uses the theoretical and empirical relations between mass, temperature and X-ray cluster luminosity, and fully accounts for the redshift evolution of the underlying dark matter clustering and cluster bias factor. The comparison between observational results and theoretical predictions shows that the Einsteinde Sitter models display too low a correlation length, while models with a matter density parameter _0m=0.3 (with or without a cosmological constant) are successful in reproducing the observed clustering. The dependence of the correlation length r ₀ on the X-ray limiting flux and luminosity of the sample is generally consistent with the predictions of all our models. Quantitative agreement is however only reached for _0m=0.3 models. The model presented here can be reliably applied to future deeper X-ray cluster surveys: the study of the clustering properties will provide a useful complementary tool to the traditional cluster abundance analyses used to constrain the cosmological parameters.     

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.     

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.     

Measuring Ω0 using cluster evolution

The evolution of the abundance of galaxy clusters depends sensitively on the value of the cosmological density parameter, Ω₀. Recent ASCA data are used to quantify this evolution as measured by the cluster X-ray temperature function. A χ² minimization fit to the cumulative temperature function, as well as a maximum-likelihood estimate (which requires additional assumptions about cluster luminosities), leads to the estimate Ω₀ ≈ 0.45 ± 0.25 (1σ statistical error). Various systematic uncertainties are considered, none of which significantly enhances the probability that Ω₀ = 1. These conclusions hold for models with or without a cosmological constant, i.e., with Λ₀ = 0 or Λ₀ = 1 − Ω₀. The statistical uncertainties are at least as large as any of the individual systematic errors that have been considered here, suggesting that additional temperature measurements of distant clusters will allow an improvement in this estimate. An alternative method that uses the highest redshift clusters to place an upper limit on Ω₀ is also presented and tentatively applied, with the result that Ω₀  1 can be ruled out at the 98 per cent confidence level. Whilst this method does not require a well-defined statistical sample of distant clusters, there are still modelling uncertainties that preclude a firmer conclusion at this time.     

The Durham/UKST Galaxy Redshift Survey – IV. Redshift-space distortions in the two-point correlation function

We have investigated the redshift-space distortions in the optically selected Durham/UKST Galaxy Redshift Survey using the two-point galaxy correlation function perpendicular and parallel to the observer's line of sight, ξ(σ, π). On small, non-linear scales we observe an elongation of the constant ξ(σ, π) contours in the line-of-sight direction. This is a result of the galaxy velocity dispersion and is the common 'Finger of God' effect seen in redshift surveys. Our result for the one-dimensional pairwise rms velocity dispersion is 〈 w ²〉^1/2=416±36 km s⁻¹, which is consistent with those from recent redshift surveys and canonical values, but inconsistent with SCDM or LCDM models. On larger, linear scales we observe a compression of the ξ(σ, π) contours in the line-of-sight direction. This is caused by the infall of galaxies into overdense regions, and the Durham/UKST data favours a value of (Ω^0.6/ b )∼0.5, where Ω is the mean mass density of the Universe and b is the linear bias factor that relates the galaxy and mass distributions. Comparison with other optical estimates yields consistent results, with the conclusion that the data do not favour an unbiased critical-density universe.     

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

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.     

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.     

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.     

Constraints on σ8 from galaxy clustering in N-body simulations and semi-analytic models

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 σ₈ 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  σ₈= 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 σ₈ 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.     

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.     

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  .