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.     

On the errors on Ω0: Monte Carlo simulations of the EMSS cluster sample

We perform Monte Carlo simulations of synthetic EMSS cluster samples, to quantify the systematic errors and the statistical uncertainties on the estimate of Ω₀ derived from fits to the cluster number density evolution and to the X-ray temperature distribution up to z =0.83 . We identify the scatter around the relation between cluster X-ray luminosity and temperature to be a source of systematic error, of the order of Δ_systΩ₀=0.09 , if not properly taken into account in the modelling. After correcting for this bias, our best Ω₀ is 0.66. The uncertainties on the shape and normalization of the power spectrum of matter fluctuations imply relatively large uncertainties on this estimate of Ω₀, of the order of Δ_statΩ₀=0.1 at the 1 σ level. On the other hand, the statistical uncertainties due to the finite size of the high-redshift sample are twice as small. Therefore, what is needed in order to improve the accuracy of Ω₀ estimates based on cluster number density evolution is a more reliable measure of the local temperature function and a better understanding of the cluster observed properties both in the local Universe and at high redshift, that is the relation between cluster mass, temperature and luminosity. This requires detailed observations of X-ray selected cluster samples, in comparison with hydrodynamic simulations including refined physics.     

Local Group velocity versus gravity: the coherence function

In maximum-likelihood analyses of the Local Group (LG) acceleration, the object describing non-linear effects is the coherence function (CF), i.e. the cross-correlation coefficient of the Fourier modes of the velocity and gravity fields. We study the CF both analytically, using perturbation theory, and numerically, using a hydrodynamic code. The dependence of the function on Ω_m and the shape of the power spectrum is very weak. The only cosmological parameter that the CF is strongly sensitive to is the normalization σ ₈ of the underlying density field. A perturbative approximation for the function turns out to be accurate as long as σ ₈ is smaller than about 0.3. For higher normalizations we provide an analytical fit for the CF as a function of σ ₈ and the wavevector. The characteristic decoherence scale which our formula predicts is an order of magnitude smaller than that determined by Strauss et al. This implies that present likelihood constraints on cosmological parameters from analyses of the LG acceleration are significantly tighter than hitherto reported.     

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.     

The XMM Cluster Survey: forecasting cosmological and cluster scaling-relation parameter constraints

We forecast the constraints on the values of  σ₈, Ω_m  and cluster scaling-relation parameters which we expect to obtain from the XMM Cluster Survey (XCS). We assume a flat Λ cold dark matter Universe and perform a Monte Carlo Markov Chain analysis of the evolution of the number density of galaxy clusters that takes into account a detailed simulated selection function. Comparing our current observed number of clusters shows good agreement with predictions. We determine the expected degradation of the constraints as a result of self-calibrating the luminosity–temperature relation (with scatter), including temperature measurement errors, and relying on photometric methods for the estimation of galaxy cluster redshifts. We examine the effects of systematic errors in scaling relation and measurement error assumptions. Using only  ( T , z )  self-calibration, we expect to measure Ω_m to ±0.03 (and  Ω_Λ  to the same accuracy assuming flatness), and σ₈ to ±0.05, also constraining the normalization and slope of the luminosity–temperature relation to ±6 and ±13 per cent (at 1σ), respectively, in the process. Self-calibration fails to jointly constrain the scatter and redshift evolution of the luminosity–temperature relation significantly. Additional archival and/or follow-up data will improve on this. We do not expect measurement errors or imperfect knowledge of their distribution to degrade constraints significantly. Scaling-relation systematics can easily lead to cosmological constraints 2σ or more away from the fiducial model. Our treatment is the first exact treatment to this level of detail, and introduces a new 'smoothed ML' (Maximum Likelihood) estimate of expected constraints.     

Dynamics of the boxy elliptical galaxy NGC 1600

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

Dark energy and the evolution of spherical overdensities

We use the non-linear spherical model in cold dark matter (CDM) cosmologies with dark energy to investigate the effects of dark energy on the growth of structure and the formation of virialized structures. We consider dark energy models with a constant equation-of-state parameter w . For  −1 < w < −1/3  , clusters form earlier and are more concentrated in quintessence than in ΛCDM models, but they form later and are less concentrated than in the corresponding open model with the same matter density and no dark energy. We point out some confusion in the literature around the expression of the collapse factor (ratio of the radius of the sphere at virialization to that at turnaround) derived from the virial theorem. We use the Sheth & Tormen extension of the Press–Schechter framework to calculate the evolution of the cluster abundance in different models and show the sensitivity of the cluster abundance to both the amplitude of the mass fluctuations, σ₈, and the σ₈– w normalization, selected to match either the cosmic microwave background observations or the abundance of X-ray clusters.     

Investigating cosmological weak lensing with the ray-bundle method

Using the ray-bundle method for calculating gravitational lens magnifications, we outline a method by which the magnification probability may be determined specifically in the weak lensing limit for cosmological models obtained from N -body simulations.
16 different models are investigated, which are variations on three broad classes of cold dark matter model: the standard model with  (Ω₀, λ ₀)=(1.0,0.0)  , the open model with  (Ω₀, λ ₀)=(0.3,0.0)  and the lambda model, which is a flat model with a cosmological constant  (Ω₀, λ ₀)=(0.3,0.7)  .
The effects of varying the Hubble parameter, H ₀, the power spectrum shape parameter, Γ, and the cluster mass normalization, σ ₈, are studied. It is shown that there is no signature of these parameters in the weak lensing magnification distributions. The magnification probability distributions are also shown to be independent of the numerical parameters such as the lens mass and simulation box size in the N -body simulations.     

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.     

LoCuSS: the connection between brightest cluster galaxy activity, gas cooling and dynamical disturbance of X-ray cluster cores

We study the distribution of projected offsets between the cluster X-ray centroid and the brightest cluster galaxy (BCG) for 65 X-ray-selected clusters from the Local Cluster Substructure Survey, with a median redshift of   z = 0.23  . We find a clear correlation between X-ray/BCG projected offset and the logarithmic slope of the cluster gas density profile at  0.04 r ₅₀₀(α  ), implying that more dynamically disturbed clusters have weaker cool cores. Furthermore, there is a close correspondence between the activity of the BCG, in terms of detected Hα and radio emission, and the X-ray/BCG offset, with the line-emitting galaxies all residing in clusters with X-ray/BCG offsets of ≤15 kpc. Of the BCGs with  α < −0.85  and an offset <0.02 r ₅₀₀, 96 per cent (23/24) have optical emission and 88 per cent (21/24) are radio active, while none has optical emission outside these criteria. We also study the cluster gas fraction ( f _gas) within r ₅₀₀ and find a significant correlation with X-ray/BCG projected offset. The mean f _gas of the 'small offset' clusters (<0.02 r ₅₀₀) is  0.106 ± 0.005 (σ= 0.03  ) compared to  0.145 ± 0.009 (σ= 0.04  ) for those with an offset >0.02 r ₅₀₀, indicating that the total mass may be systematically underestimated in clusters with larger X-ray/BCG offsets. Our results imply a link between cool core strength and cluster dynamical state consistent with the view that cluster mergers can significantly perturb cool cores, and set new constraints on models of the evolution of the intracluster medium.     

Measuring Ω0 from the entropy evolution of clusters

In this paper we have extended the entropy-driven model of cluster evolution developed by Bower in order to be able to predict the evolution of galaxy clusters for a range of cosmological scenarios. We have applied this model to recent measurements of the evolution of the L _x− T normalization and X-ray luminosity function in order to place constraints on cosmological parameters. We find that these measurements alone do not select a particular cosmological framework. An additional constraint is required on the effective slope of the power spectrum to break the degeneracy that exists between this and the background cosmology. We therefore include a theoretical calculation of the Ω₀ dependence on the power spectrum, based on the cold dark matter paradigm, which infers Ω₀<0.55 (0.1<Ω₀<0.7 for Ω₀+Λ₀=1), at the 95 per cent confidence level. Alternatively, an independent measurement of the slope of the power spectrum from galaxy clustering requires Ω₀<0.6 (Ω₀<0.65 for Ω₀+Λ₀=1), again to 95 per cent confidence. The rate of entropy evolution is insensitive to the values of Ω₀ considered, although it is sensitive to changes in the distribution of the intracluster medium.     

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.     

Stellar velocity dispersion in narrow-line Seyfert 1 galaxies

Several authors have recently explored, for narrow-line Seyfert 1 galaxies (NLS1s), the relationship between black hole mass ( M _BH) and stellar velocity dispersion (σ_*). Their results are more or less in agreement and seem to indicate that NLS1s fill the region below the fit obtained by Tremaine et al., showing a range of σ_* similar to that of Seyfert 1 galaxies, and a lower M _BH. Until now, the [O  iii ] width has been used in place of the stellar velocity dispersion, but some indications have begun to arise against the effectiveness of the gaseous kinematics in representing the bulge potential, at least in NLS1s. Bian & Zhao have stressed the urgency of producing true σ_* measurements. Here, we present new stellar velocity dispersions obtained through direct measurements of the Ca  ii absorption triplet (∼8550 Å) in the nuclei of eight NLS1 galaxies. The resulting σ_* values and a comparison with σ_[O III] confirm our suspicion that [O  iii ] typically overestimates the stellar velocity dispersion. We demonstrate that NLS1s follow the   M _BH–σ_*  relation as Seyfert 1, quasars and non-active galaxies.     

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  .     

Can simulations reproduce the observed temperature–mass relation for clusters of galaxies?

It has become increasingly apparent that traditional hydrodynamical simulations of galaxy clusters are unable to reproduce the observed properties of galaxy clusters, in particular overpredicting the mass corresponding to a given cluster temperature. Such overestimation may lead to systematic errors in results using galaxy clusters as cosmological probes, such as constraints on the density perturbation normalization σ ₈. In this paper we demonstrate that inclusion of additional gas physics, namely radiative cooling and a possible pre-heating of gas prior to cluster formation, is able to bring the temperature–mass relation in the innermost parts of clusters into good agreement with recent determinations by Allen, Schmidt & Fabian using Chandra data.     

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.     

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.     

Redshift and velocity dispersion of the cluster of galaxies around NGC 326

Redshifts of several galaxies thought to be associated with NGC 326 are determined. The results confirm the presence of a cluster and find a mean redshift of     and a line-of-sight velocity dispersion σ _z =599 (+230,−110) km s⁻¹. The velocity dispersion and previously measured X-ray gas temperature of kT ≃1.9 keV are consistent with the cluster σ _z kT relation, and NGC 326 is seen to be a slowly moving member of the cluster.     

Revisiting the baryon fractions of galaxy clusters: a comparison with WMAP 3-yr results

The universal baryonic mass fraction  (Ω_b/Ω_m)  can be sensitively constrained using X-ray observations of galaxy clusters. In this paper, we compare the baryonic mass fraction inferred from measurements of the cosmic microwave background with the gas mass fractions ( f _gas) of a large sample of clusters taken from the recent literature. In systems cooler than 4 keV, f _gas declines as the system temperature decreases. However, in higher temperature systems, f _gas( r ₅₀₀) converges to  ≈(0.12 ± 0.02)( h /0.72)^−1.5  , where the uncertainty reflects the systematic variations between clusters at r ₅₀₀. This is significantly lower than the maximum-likelihood value of the baryon fraction from the recently released Wilkinson Microwave Anisotropy Probe ( WMAP ) 3-yr results. We investigate possible reasons for this discrepancy, including the effects of radiative cooling and non-gravitational heating, and conclude that the most likely solution is that Ω_m is higher than the best-fitting WMAP value (we find  Ω_m= 0.36^+0.11_−0.08  ), but consistent at the 2σ level. Degeneracies within the WMAP data require that σ₈ must also be greater than the maximum likelihood value for consistency between the data sets.     

The STAGES view of red spirals and dusty red galaxies: mass-dependent quenching of star formation in cluster infall

We investigate the properties of optically passive spirals and dusty red galaxies in the A901/2 cluster complex at redshift ∼0.17 using rest-frame near-ultraviolet–optical spectral energy distributions, 24-μm infrared data and Hubble Space Telescope morphologies from the STAGES data set. The cluster sample is based on COMBO-17 redshifts with an rms precision of  σ _cz ≈ 2000 km s⁻¹  . We find that 'dusty red galaxies' and 'optically passive spirals' in A901/2 are largely the same phenomenon, and that they form stars at a substantial rate, which is only four times lower than that in blue spirals at fixed mass. This star formation is more obscured than in blue galaxies and its optical signatures are weak. They appear predominantly in the stellar mass range of  log  M _*/M_⊙=[10, 11]  where they constitute over half of the star-forming galaxies in the cluster; they are thus a vital ingredient for understanding the overall picture of star formation quenching in clusters. We find that the mean specific star formation rate (SFR) of star-forming galaxies in the cluster is clearly lower than in the field, in contrast to the specific SFR properties of blue galaxies alone, which appear similar in cluster and field. Such a rich red spiral population is best explained if quenching is a slow process and morphological transformation is delayed even more. At  log  M _*/M_⊙ < 10  , such galaxies are rare, suggesting that their quenching is fast and accompanied by morphological change. We note that edge-on spirals play a minor role; despite being dust reddened they form only a small fraction of spirals independent of environment.