Baryon oscillations in galaxy and matter power-spectrum covariance matrices

We investigate large-amplitude baryon acoustic oscillations (BAOs) in off-diagonal entries of cosmological power-spectrum covariance matrices. These covariance-matrix BAOs describe the increased attenuation of power-spectrum BAOs caused by upward fluctuations in large-scale power. We derive an analytic approximation to covariance-matrix entries in the BAO regime, and check the analytical predictions using N -body simulations. These BAOs look much stronger than the BAOs in the power spectrum, but seem detectable only at about a 1σ level in gigaparsec-scale galaxy surveys. In estimating cosmological parameters using matter or galaxy power spectra, including the covariance-matrix BAOs can have a several per cent effect on error bar widths for some parameters directly related to the BAOs, such as the baryon fraction. Also, we find that including the numerous galaxies in small haloes in a survey can reduce error bars in these cosmological parameters more than the simple reduction in shot noise might suggest.     

Clustering of galaxy clusters in cold dark matter universes

We use very large cosmological N -body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ( τ CDM) and a flat low density model with a cosmological constant (ΛCDM). Our simulations each use 10⁹ particles to follow the mass distribution within cubes of side 2  h ⁻¹ Gpc ( τ CDM) and 3  h ⁻¹ Gpc (ΛCDM) with a force resolution better than 10⁻⁴ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from     for samples with mean separation     to     for samples with     The lower value here corresponds to τ CDM and the upper to ΛCDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10 per cent using the simple model of Sheth, Mo & Tormen. This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the automated plate measurement (APM) cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for ΛCDM. The observed power spectrum (for a cluster sample of mean separation     ) lies significantly above the predictions of both models.     

Testing deprojection algorithms on mock angular catalogues: evidence for a break in the power spectrum

We produce mock angular catalogues from simulations with different initial power spectra to test methods that recover measures of clustering in three dimensions, such as the power spectrum, variance and higher order cumulants. We find that the statistical properties derived from the angular mock catalogues are in good agreement with the intrinsic clustering in the simulations. In particular, we concentrate on the detailed predictions for the shape of the power spectrum, P ( k ). We find that there is good evidence for a break in the galaxy P ( k ) at scales in the range 0.02< k <0.06 h Mpc⁻¹, using an inversion technique applied to the angular correlation function measured from the APM Galaxy Survey. For variants on the standard cold dark matter (CDM) model, a fit at the location of the break implies Ω h =0.45±0.10, where Ω is the ratio of the total matter density to the critical density, and Hubble's constant is parametrized as H ₀=100 h km s⁻¹ Mpc⁻¹. On slightly smaller, though still quasi-linear scales, there is a feature in the APM power spectrum where the local slope changes appreciably, with the best match to CDM models obtained for Ω h ≃0.2. Hence the location and narrowness of the break in the APM power spectrum combined with the rapid change in its slope on quasi-linear scales cannot be matched by any variant of CDM, including models that have a non-zero cosmological constant or a tilt to the slope of the primordial P ( k ). These results are independent of the overall normalization of the CDM models or any simple bias that exists betwen the galaxy and mass distributions.     

SDSS galaxy clustering: luminosity and colour dependence and stochasticity

Differences in clustering properties between galaxy subpopulations complicate the cosmological interpretation of the galaxy power spectrum, but can also provide insights about the physics underlying galaxy formation. To study the nature of this relative clustering, we perform a counts-in-cells analysis of galaxies in the Sloan Digital Sky Survey in which we measure the relative bias between pairs of galaxy subsamples of different luminosities and colours. We use a generalized  χ²  test to determine if the relative bias between each pair of subsamples is consistent with the simplest deterministic linear bias model, and we also use a maximum likelihood technique to further understand the nature of the relative bias between each pair. We find that the simple, deterministic model is a good fit for the luminosity-dependent bias on scales above  ∼2  h ⁻¹ Mpc  , which is good news for using magnitude-limited surveys for cosmology. However, the colour-dependent bias shows evidence for stochasticity and/or non-linearity which increases in strength towards smaller scales, in agreement with previous studies of stochastic bias. Also, confirming hints seen in earlier work, the luminosity-dependent bias for red galaxies is significantly different from that of blue galaxies: both luminous and dim red galaxies have higher bias than moderately bright red galaxies, whereas the biasing of blue galaxies is not strongly luminosity dependent. These results can be used to constrain galaxy formation models and also to quantify how the colour and luminosity selection of a galaxy survey can impact measurements of the cosmological matter power spectrum.     

Maximum likelihood estimates of the two- and three-dimensional power spectra of the APM Galaxy Survey

We estimate the two- and three-dimensional power spectra, P ₂( K ) and P ₃( 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 ₂( K ) and P ₃( 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 ₃( k ) estimates at k ≤0.19  h  Mpc⁻¹ constrain the amplitude and shape parameter of a scale-invariant CDM model to lie within the 2 σ ranges 0.74≤( σ ₈)_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.     

Simulations of baryon oscillations

The coupling of photons and baryons by Thomson scattering in the early universe imprints features in both the cosmic microwave background (CMB) and matter power spectra. The former have been used to constrain a host of cosmological parameters, the latter have the potential to strongly constrain the expansion history of the universe and dark energy. Key to this program is the means to localize the primordial features in observations of galaxy spectra which necessarily involve galaxy bias, non-linear evolution and redshift space distortions. We present calculations, based on mock catalogs produced from high-resolution N-body simulations, which show the range of behaviors we might expect of galaxies in the real universe. We investigate physically motivated fitting forms which include the effects of non-linearity, galaxy bias and redshift space distortions and discuss methods for analysis of upcoming data. In agreement with earlier work, we find that a survey of several Gpc³ would constrain the sound horizon at z  1 to about 1%.     

Peaks in the cosmological density field: sensitivity to initial power spectrum, redshift distortions and galaxy halo occupation

We investigate the number density of maxima in the cosmological galaxy density field smoothed with a filter as a probe of clustering. In previous work it has been shown that this statistic is closely related to the slope of the linear power spectrum, even when the directly measured power spectrum is non-linear. In the present paper we investigate the sensitivity of the peak number density to various models with differing power spectra, including rolling index models, cosmologies with massive neutrinos and different baryon densities. We find that rolling index models which have given an improved fit to CMB/LSS (cosmic microwave background/large scale structure) data yield a ∼10 per cent difference in peak density compared to the scale invariant case. Models with 0.3 eV neutrinos have effects of similar magnitude and it should be possible to constrain them with data from current galaxy redshift surveys. Baryon oscillations in the power spectrum also give rise to distinctive features in the peak density. These are preserved without modification when measured from the peak density in fully non-linear N -body simulations. Using the simulations, we also investigate how the peak density is modified in the presence of redshift distortions. Redshift distortions cause a suppression of the number of peaks, largely due to fingers of God overlapping in redshift space. We find that this effect can be modelled by using a modification of the input power spectrum. We also study the results when the simulation density field is traced by galaxies obtained by populating haloes with a halo occupation distribution consistent with observations. The peak number density is consistent with that in the dark matter for filter scales  >4  h ⁻¹ Mpc  , for which we find good agreement with the linear theory predictions. In a companion paper we analyse data from the 2dF Galaxy Redshift Survey.     

Scale-dependent bias induced by local non-Gaussianity: a comparison to N-body simulations

We investigate the effect of primordial non-Gaussianity of the local f _NL type on the auto- and cross-power spectra of dark matter haloes using simulations of the Λ cold dark matter cosmology. We perform a series of large N -body simulations of both positive and negative f _NL, spanning the range between 10 and 100. Theoretical models predict a scale-dependent bias correction  Δ b ( k , f _NL)  that depends on the linear halo bias   b ( M )  . We measure the power spectra for a range of halo mass and redshifts covering the relevant range of existing galaxy and quasar populations. We show that auto- and cross-correlation analyses of bias are consistent with each other. We find that for low wavenumbers with   k < 0.03  h  Mpc⁻¹  the theory and the simulations agree well with each other for biased haloes with   b ( M ) > 1.5  . We show that a scale-independent bias correction improves the comparison between theory and simulations on smaller scales, where the scale-dependent effect rapidly becomes negligible. The current limits on f _NL from Slosar et al. come mostly from very large scales   k < 0.01  h  Mpc⁻¹  and, therefore, remain valid. For the halo samples with   b ( M ) < 1.5 − 2  , we find that the scale-dependent bias from non-Gaussianity actually exceeds the theoretical predictions. Our results are consistent with the bias correction scaling linearly with f _NL.     

Measurement of the baryonic acoustic oscillation scale in 21 cm intensity fluctuations during the reionization era

It has recently been suggested that the power spectrum of redshifted 21 cm fluctuations could be used to measure the scale of baryonic acoustic oscillations (BAOs) during the reionization era. The resulting measurements are potentially as precise as those offered by the next generation of galaxy redshift surveys at lower redshift. However, unlike galaxy redshift surveys, which in the linear regime are subject to a scale-independent galaxy bias, the growth of ionized regions during reionization is thought to introduce a strongly scale-dependent relationship between the 21 cm and mass power spectra. We use a seminumerical model for reionization to assess the impact of ionized regions on the precision and accuracy with which the BAO scale could be measured using redshifted 21 cm observations. For a model in which reionization is completed at   z ∼ 6  , we find that the constraints on the BAO scale are not systematically biased at   z ≳ 6.5  . In this scenario, and assuming the sensitivity attainable with a low-frequency array comprising 10 times the collecting area of the Murchison Widefield Array, the BAO scale could be measured to within 1.5 per cent in the range  6.5 ≲ z ≲ 7.5  .     

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.     

Voids in the Las Campanas Redshift Survey versus cold dark matter models

We analysed the distribution of void sizes in the two-dimensional slices of the Las Campanas Redshift Survey (LCRS). Fourteen volume-limited subsamples were extracted from the six slices to cover a large part of the survey and to test the robustness of the results against cosmic variance. Thirteen samples were randomly culled to produce homogeneously selected samples. We then studied the relationship between the cumulative area covered by voids and the void size as a property of the void hierarchy. We found that the distribution of void sizes scaled with the mean galaxy separation, λ . In particular, we found that the size of voids covering half of the area is given by D _med≈ λ +(12±3)  h ⁻¹ Mpc. Next, by employing an environmental density threshold criterion to identify mock galaxies, we were able to extend this analysis to mock samples from dynamical N -body simulations of cold dark matter (CDM) models. To reproduce the observed void statistics, overdensity thresholds of δ _th≈0,…,1 are necessary. We compared standard (SCDM), open (OCDM), vacuum energy dominated (ΛCDM) and broken scale invariant CDM models (BCDM): we found that both the void size distribution and the two-point correlation function provided important and complementary information on the large-scale matter distribution. The dependence of the void statistics on the threshold criterion for the mock galaxy identification showed that the galaxy biasing was more crucial for the void size distribution than were differences between the cosmological models.     

Power spectrum analysis of the ESO Slice Project galaxy redshift survey

We measure the power spectrum of the galaxy distribution in the ESO Slice Project (ESP) galaxy redshift survey. We develop a technique to describe the survey window function analytically, and then deconvolve it from the measured power spectrum using a variant of the Lucy method. We test the whole deconvolution procedure on ESP mock catalogues drawn from large N -body simulations, and find that it is reliable for recovering the correct amplitude and shape of P ( k ) at k >0.065  h  Mpc⁻¹. In general, the technique is applicable to any survey composed of a collection of circular fields with an arbitrary pattern on the sky, as typical of surveys based on fibre spectrographs. The estimated power spectrum has a well-defined power-law shape k ⁿ with n ≃−2.2 for k ≥0.2  h  Mpc⁻¹, and a smooth bend to a flatter shape ( n ≃−1.6) for smaller k . The smallest wavenumber where a meaningful reconstruction can be performed ( k ∼0.06  h  Mpc⁻¹) does not allow us to explore the range of scales where other power spectra seem to show a flattening and hint at a turnover. We also find, by a direct comparison of the Fourier transforms, that the estimate of the two-point correlation function ξ ( s ) is much less sensitive to the effect of a problematic window function, such as that of the ESP, than the power spectrum. Comparison with other surveys shows an excellent agreement with estimates from blue-selected surveys. In particular, the ESP power spectrum is virtually indistinguishable from that of the Durham–UKST survey over the common range of k , an indirect confirmation of the quality of the deconvolution technique applied.     

Evolutionary Stellar Population Synthesis at 2 Å Spectral Resolution

I present an evolutionary stellar population synthesis model which predicts spectral energy distributions, SEDs, for simple old single-metallicity stellar populations, SSPs, in the wavelength intervalsλλ 3856–4476Å and 4795–5465Å at a resolution of FWHM$equals;1.8 $Aring;, on the basis of an extensive empirical spectral library composed of ≈550 stars. The synthesized model spectra can be used to analyse observed galaxy spectra in a very easy and flexible way, allowing us to adapt the theoretical predictions to the characteristics of the data instead of proceeding in the opposite direction (as we must do, for example, when transforming observational data to a particular system of indices at specific resolution/s, such as Lick, which is heavily instrument-dependent).The SSP spectra, with flux-calibrated response curves, can be smoothed to the same resolution as that of the data or to galaxy internal velocity dispersion, allowing us to analyse the observed spectra in their own system, and with no need to correct the index measurements for velocity dispersion. Thus, we are able to use all the information contained in the data, at their higher spectral resolution. Excellent fits are obtained for various metal-rich globular clusters at relatively high resolution, and well known spectral peculiarities such as the strong CN absorption features are detected. When applied to early-type galaxies the model also shows its potential for studying element ratios such as the Mg/Fe overabundance. The model spectra provided robust age indicators for old stellar populations which do not depend on metallicity and therefore have a great potential for solving the age–metallicitydegeneracy of galaxy spectra.     

Testing cosmological structure formation using redshift-space distortions

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for observing the build-up of cosmological structure. In this paper, we develop and test a new statistic based on anisotropies in the measured galaxy power spectrum, which is independent of galaxy bias and matches the matter power spectrum shape on large scales. The amplitude provides a constraint on the derivative of the linear growth rate through   f σ₈(mass)  . This demonstrates that spectroscopic galaxy surveys offer many of the same advantages as weak lensing surveys, in that they both use galaxies as test particles to probe all matter in the Universe. They are complementary as redshift-space distortions probe non-relativistic velocities and therefore the temporal metric perturbations, while weak lensing tests the sum of the temporal and spatial metric perturbations. The degree to which our estimator can be pushed into the non-linear regime is considered and we show that a simple Gaussian damping model, similar to that previously used to model the behaviour of the power spectrum on very small scales, can also model the quasi-linear behaviour of our estimator. This enhances the information that can be extracted from surveys for Λ cold dark matter (ΛCDM) models.     

A full-sky prediction of the Sunyaev–Zeldovich effect from diffuse hot gas in the local universe and the upper limit from the WMAP data

We use the Point Source Catalogue Redshift Survey galaxy redshift catalogue combined with constrained simulations based on the IRAS 1.2-Jy galaxy density field to estimate the contribution of hot gas in the local universe to the Sunyaev–Zeldovich (SZ) effect on a large scale. We produce a full-sky healpix map predicting the SZ effect from clusters as well as diffuse hot gas within  80  h ⁻¹ Mpc  . Performing cross-correlation tests between this map and the WMAP data in pixel, harmonic and wavelet space we can put an upper limit on the effect. We conclude that the SZ effect from diffuse gas in the local universe cannot be detected in current cosmic microwave background (CMB) data and is not a large-scale contaminating factor  (ℓ < 60)  in studies of CMB angular anisotropies. We derive an upper limit for the mean temperature decrement of  Δ T < 0.33 μK  at the 2σ confidence level for the 61-GHz frequency channel. However, for future high-sensitivity experiments observing at a wider range of frequencies, the predicted large-scale SZ effect could be of importance.     

Cosmological implications of the PSCz PDF and its moments

We compare the probability density function (PDF) and its low-order moments (variance and skewness) of the smoothed IRAS Point Source Catalogue Redshift Survey (PSC z ) galaxy density field and of the corresponding simulated PSC z look-alikes, generated from N -body simulations of six different dark matter models: four structure-normalized with     and     , one COBE -normalized, and the old standard cold dark matter model. The galaxy distributions are smoothed with a Gaussian window at three different smoothing scales,     , 10 and 15  h ⁻¹ Mpc. We find that the simulation PSC z look-alike PDFs are sensitive only to the normalization of the power spectrum, probably owing to the shape similarity of the simulated galaxy power spectrum on the relevant scales. We find that the only models that are consistent, at a high significance level, with the observed PSC z PDF are models with a relatively low power spectrum normalization     . From the phenomenologically derived σ ₈–moments relation, fitted from the simulation data, we find that the PSC z moments suggest     .     

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

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