Sunyaev–Zel'dovich cluster survey simulations for Planck

We examine the ability of the future Planck mission to provide a catalogue of galaxy clusters observed via their Sunyaev–Zel'dovich (SZ) distortion in the cosmic microwave background (CMB). For this purpose we produce full-sky SZ maps based on N -body simulations and scaling relations between cluster properties for several cosmological models. We extrapolate the N -body simulations by a mass function to high redshifts in order to obtain a realistic SZ background. The simulated Planck observations include, besides the thermal and kinematic SZ effects, contributions from the primordial CMB, extragalactic point sources as well as Galactic dust, free–free and synchrotron emission. A harmonic-space maximum-entropy method is used to separate the SZ signal from contaminating components in combination with a cluster detection algorithm based on thresholding and flux integration to identify clusters and to obtain their fluxes. We estimate a survey sensitivity limit (depending on the quality of the recovered cluster flux) and provide cluster survey completeness and purity estimates. We find that, given our modelling and detection algorithm, Planck will reliably detect at least several thousands of clusters over the full sky. The exact number depends on the particular cosmological model (up to 10 000 cluster detections in a concordance ΛCDM model with  σ₈= 0.9  ). We show that the Galaxy does not significantly affect the cluster detection. Furthermore, the dependence of the thermal SZ power spectrum on the matter variance on scales of  8 h ⁻¹  Mpc and the quality of its reconstruction by the employed method are investigated. Our simulations suggest that the Planck cluster sample will not only be useful as a basis for follow-up observations, but also will have the ability to provide constraints on cosmological parameters.     

The Sunyaev–Zel'dovich effect contribution to WMAP: a cross-correlation between WMAP and ROSAT

We cross-correlate WMAP and ROSAT diffuse X-ray background maps and look for common features in both data sets. We use the power spectrum of the product maps and the cross-power spectrum to highlight a possible correlation. The power spectrum of the product maps does not detect any correlation and the cross-power spectrum does not show any significant deviation from zero. We explore different explanations for this lack of correlation. A universe with a low value of  σ₈  could naturally explain the lack of correlation. We also discuss the systematic effects that can affect this result, in particular the subtraction of some cluster signal from the ROSAT diffuse maps, which could significantly suppress the correlation signal. These systematic effects considerably reduce the significance of our constraints on the cosmological model. When we include the systematic effects, we find a weaker constraint on  σ₈  , allowing models with values as large as  σ₈= 1  (for  Ω_m= 0.3  ) to be consistent with the lack of correlation. To illustrate the capabilities of the method with future high-quality data, we show how from the correlation signal it should be possible to predict the level of contamination of the Sunyaev–Zel'dovich effect on the power spectrum of the cosmic microwave background. Within the systematic errors, we find evidence that this contribution is negligible for WMAP and is expected to be small in experiments like ACBAR or CBI, but can be important for future high-resolution experiments.     

Constraints on structure formation models from the Sunyaev–Zel'dovich effect

In the context of cold dark matter (CDM) cosmological models, we have simulated images of the brightness temperature fluctuations in the cosmic microwave background (CMB) sky owing to the Sunyaev–Zel'dovich (S–Z) effect in a cosmological distribution of clusters. We compare the image statistics with recent ATCA limits on arcmin-scale CMB anisotropy. The S–Z effect produces a generically non-Gaussian field and we compute the variance in the simulated temperature-anisotropy images, after convolution with the ATCA beam pattern, for different cosmological models. All the models are normalized to the 4-yr COBE data. We find an increase in the simulated-sky temperature variance with increase in the cosmological density parameter Ω₀. A comparison with the upper limits on the sky variance set by the ATCA appears to rule out our closed-universe model: low-Ω₀ open-universe models are preferred. The result is independent of any present day observations of σ ₈.     

Sunyaev–Zel'dovich distortion from early galactic winds

We consider the distortion in the cosmic microwave background (CMB) resulting from galactic winds at high redshift. Winds outflowing from galaxies have been hypothesized to be possible sources of metals in the intergalactic medium, which is known to have been enriched to 10^−2.5 Z_⊙ at z ∼3. We model these winds as functions of mass of the parent galaxy and redshift, assuming that they activate at a common initial redshift, z _in, and calculate the mean y -distortion and the angular power spectrum of the distortion in the CMB. We find that the thermal Sunyaev–Zel'dovich (SZ) effect resulting from the winds is consistent with previous estimates. The distortion arising from the kinetic SZ (kSZ) effect is, however, found to be more important than the thermal SZ (tSZ) effect. We find that the distortion resulting from galactic winds is an important contribution to the power spectrum of distortion at very small angular scales ( l ∼10⁴). We also find that the power spectrum resulting from clustering dominates the Poisson power spectrum for l ≤(4–5)×10⁵. We show explicitly how the combined power spectrum from wind dominates over that of clusters at 217 GHz, relevant for PLANCK . We also show how these constraints change when the efficiency of the winds is varied.     

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.     

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.     

Cosmological parameters from velocities, cosmic microwave background and supernovae

We compare and combine likelihood functions of the cosmological parameters Ω_m, h and σ ₈, from peculiar velocities, cosmic microwave background (CMB) and type Ia supernovae. These three data sets directly probe the mass in the Universe, without the need to relate the galaxy distribution to the underlying mass via a 'biasing' relation. We include the recent results from the CMB experiments BOOMERANG and MAXIMA-1. Our analysis assumes a flat Λ cold dark matter (ΛCDM) cosmology with a scale-invariant adiabatic initial power spectrum and baryonic fraction as inferred from big-bang nucleosynthesis. We find that all three data sets agree well, overlapping significantly at the 2 σ level. This therefore justifies a joint analysis, in which we find a joint best-fitting point and 95 per cent confidence limits of     (0.17,0.39),     (0.64,0.86) and     (0.98,1.37). In terms of the natural parameter combinations for these data     (0.40,0.73),     (0.16,0.27). Also for the best-fitting point,     and the age of the Universe is 13.2 Gyr.     

The location of cosmic microwave background peaks in a universe with dark energy

The locations of the peaks of the cosmic microwave background (CMB) spectrum are sensitive indicators of cosmological parameters, yet there is no known analytic formula which accurately describes their dependence on them. We parametrize the location of the peaks as   l _m = l _A( m - φ _m )  , where l _A is the analytically calculable acoustic scale and m labels the peak number. Fitting formulae for the phase shifts φ _m for the first three peaks and the first trough are given. It is shown that in a wide range of parameter space, the acoustic scale l _A can be retrieved from actual CMB measurements of the first three peaks within 1 per cent accuracy. This can be used to speed up likelihood analysis. We describe how the peak shifts can be used to distinguish between different models of dark energy.     

Implications of bias evolution on measurements of the integrated Sachs–Wolfe effect: errors and biases in parameter estimation

The subject of this paper is a quantification of the impact of uncertainties in bias and bias evolution on the interpretation of measurements of the integrated Sachs–Wolfe (ISW) effect, in particular on the estimation of cosmological parameters. We carry out a Fisher matrix analysis for quantifying the degeneracies between the parameters of a dark energy cosmology and bias evolution, for the combination of the PLANCK microwave sky survey with the EUCLID main galaxy sample, where bias evolution   b ( a ) = b ₀+ (1 − a ) b_a   is modelled with two parameters b ₀ and   b_a   . Using a realistic bias model introduces a characteristic suppression of the ISW spectrum on large angular scales, due to the altered distance-weighting functions. The errors in estimating cosmological parameters if the data with evolving bias is interpreted in the framework of cosmologies with constant bias are quantified in an extended Fisher formalism. We find that the best-fitting values of all parameters are shifted by an amount comparable to the statistical accuracy: the estimation bias in units of the statistical accuracy amounts to 1.19 for Ω_m, 0.27 for σ₈ and 0.72 for w for bias evolution with   b_a = 1  . Leaving   b_a   open as a free parameter deteriorates the statistical accuracy, in particular on Ω_m and w .     

The Australia Telescope search for cosmic microwave background anisotropy

In an attempt to detect cosmic microwave background (CMB) anisotropy on arcmin scales, we have made an 8.7-GHz image of a sky region with a resolution of 2 arcmin and high surface brightness sensitivity using the Australia Telescope Compact Array (ATCA) in an ultracompact configuration. The foreground discrete-source confusion was estimated from observations with higher resolution at the same frequency and in a scaled array at a lower frequency. Following the subtraction of the foreground confusion, the field shows no features in excess of the instrument noise. This limits the CMB anisotropy flat-band power to Q _flat < 23.6 μ K with 95 per cent confidence; the ATCA filter function (which is available at the website www.atnf.csiro.au/Research/cmbr/cmbr_atca.html) F _l in multipole l -space peaks at l _eff = 4700 and has half-maximum values at l  = 3350 and 6050.     

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.     

Observational tests of x-matter models

We study gravitational lensing statistics, matter power spectra and the angular power spectra of the cosmic microwave background (CMB) radiation in x-matter models. We adopt an equation of state of x-matter which can express a wide range of matter from pressureless dust to the cosmological constant. A new ingredient in this model is the sound speed of the x-component, in addition to the equation of state w ₀ =  p _x0/ρ_x0. Except for the cosmological constant case, the perturbations of x-matter itself are considered. Our primary interest is in the effect of non-zero sound speed on the structure formation and the CMB spectra. It is found that there exist parameter ranges where x-matter models are consistent with all current observations. The x-matter generally leaves imprints in the CMB anisotropy and the matter power spectrum, which should be detectable in future observations.     

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.     

Lensed CMB temperature and polarization maps from the Millennium Simulation

We have constructed the first all-sky cosmic microwave background (CMB) temperature and polarization lensed maps based on a high-resolution cosmological N -body simulation, the Millennium Simulation (MS). We have exploited the lensing potential map obtained using a previously developed map-making procedure which integrates along the line-of-sight the MS dark matter distribution by stacking and randomizing the simulation boxes up to   z = 127  , and which semi-analytically supplies the large-scale power in the angular lensing potential that is not correctly sampled by the N -body simulation. The lensed sky has been obtained by properly modifying the latest version of the LensPix code to account for the MS structures. We have also produced all-sky lensed maps of the so-called  ψ _E   and  ψ _B   potentials, which are directly related to the electric and magnetic types of polarization. The angular power spectra of the simulated lensed temperature and polarization maps agree well with semi-analytic estimates up to   l ≤ 2500  , while on smaller scales we find a slight excess of power which we interpret as being due to non-linear clustering in the MS. We also observe how non-linear lensing power in the polarized CMB is transferred to large angular scales by suitably misaligned modes in the CMB and the lensing potential. This work is relevant in view of the future CMB probes, as a way to analyse the lensed sky and disentangle the contribution from primordial gravitational waves.     

Neural networks and the separation of cosmic microwave background and astrophysical signals in sky maps

We implement an independent component analysis (ICA) algorithm to separate signals of different origin in sky maps at several frequencies. Owing to its self-organizing capability, it works without prior assumptions on either the frequency dependence or the angular power spectrum of the various signals; rather, it learns directly from the input data how to identify the statistically independent components, on the assumption that all but, at most, one of the components have non-Gaussian distributions.
We have applied the ICA algorithm to simulated patches of the sky at the four frequencies (30, 44, 70 and 100 GHz) used by the Low Frequency Instrument of the European Space Agency's Planck satellite. Simulations include the cosmic microwave background (CMB), the synchrotron and thermal dust emissions, and extragalactic radio sources. The effects of the angular response functions of the detectors and of instrumental noise have been ignored in this first exploratory study. The ICA algorithm reconstructs the spatial distribution of each component with rms errors of about 1 per cent for the CMB, and 10 per cent for the much weaker Galactic components. Radio sources are almost completely recovered down to a flux limit corresponding to ≃0.7 σ _CMB, where σ _CMB is the rms level of the CMB fluctuations. The signal recovered has equal quality on all scales larger than the pixel size. In addition, we show that for the strongest components (CMB and radio sources) the frequency scaling is recovered with per cent precision. Thus, algorithms of the type presented here appear to be very promising tools for component separation. On the other hand, we have been dealing here with a highly idealized situation. Work to include instrumental noise, the effect of different resolving powers at different frequencies and a more complete and realistic characterization of astrophysical foregrounds is in progress.     

The orbits of open clusters in the Galaxy

We present and analyse the kinematics and orbits for a sample of 488 open clusters (OCs) in the Galaxy. The velocity ellipsoid for our present sample is derived as  (σ _U , σ _V , σ _W ) = (28.7, 15.8, 11.0) km s⁻¹  which represents a young thin-disc population. We also confirm that the velocity dispersions increase with the age of a cluster subsample. The orbits of OCs are calculated with three Galactic gravitational potential models. The errors of orbital parameters are also calculated considering the intrinsic variation of the orbital parameters and the effects of observational uncertainties. The observational uncertainties dominate the errors of derived orbital parameters. The vertical motions of clusters calculated using different Galactic disc models are rather different. The observed radial metallicity gradient of clusters is derived with a slope of   b =−0.070 ± 0.011   dex kpc⁻¹. The radial metallicity gradient of clusters based on their apogalactic distances is also derived with a slope of   b =−0.082 ± 0.014   dex kpc⁻¹. The distribution of derived orbital eccentricities for OCs is very similar to that derived for the field population of dwarfs and giants in the thin disc.     

Cross-correlation between WMAP and 2MASS: non-Gaussianity induced by the SZ effect

We study the non-Gaussianity induced by the Sunyaev–Zel'dovich (SZ) effect in cosmic microwave background (CMB) fluctuation maps. If a CMB map is contaminated by the SZ effect of galaxies or galaxy clusters, the CMB maps should have similar non-Gaussian features to the galaxy and cluster fields. Using the WMAP data and 2MASS galaxy catalogue, we show that the non-Gaussianity of the 2MASS galaxies is imprinted on WMAP maps. The signature of non-Gaussianity can be seen with the fourth-order cross-correlation between the wavelet variables of the WMAP maps and 2MASS clusters. The intensity of the fourth-order non-Gaussian features is found to be consistent with the contamination of the SZ effect of 2MASS galaxies. We also show that this non-Gaussianity can not be seen by the high-order autocorrelation of the WMAP . This is because the SZ signals in the autocorrelations of the WMAP data generally are weaker than the WMAP –2MASS cross-correlations by a factor f ², which is the ratio between the powers of the SZ-effect map and the CMB fluctuations on the scale considered. Therefore, the ratio of high-order autocorrelations of CMB maps to cross-correlations of the CMB maps and galaxy field would be effective to constrain the powers of the SZ effect on various scales.     

Decay of the vacuum energy into cosmic microwave background photons

We examine the possibility of the decay of the vacuum energy into a homogeneous distribution of a thermalized cosmic microwave background (CMB), which is characteristic of an adiabatic vacuum energy decay into photons. It is shown that observations of the primordial density fluctuation spectrum, obtained from CMB and galaxy distribution data, restrict the possible decay rate. When photon creation due to an adiabatic vacuum energy decay takes place, the standard linear temperature dependence   T ( z ) = T ₀(1 + z )  is modified, where T ₀ is the present CMB temperature, and can be parametrized by a modified CMB temperature dependence     . From the observed CMB and galaxy distribution data, a strong limit on the maximum value of the decay rate is obtained by placing a maximum value  β_max≃ 3.4 × 10⁻³  on the β parameter.     

CMB temperature polarization correlation and primordial gravitational waves

We examine the use of the TE cross-correlation power spectrum of the cosmic microwave background (CMB) as a complementary test to detect primordial gravitational waves (PGWs). The first method used is based on the determination of the lowest multipole, ℓ₀, where the TE power spectrum,   C ^TE_ℓ  , first changes sign. The second method uses Wiener filtering on the CMB TE data to remove the density perturbations contribution to the TE power spectrum. In principle this leaves only the contribution of PGWs. We examine two toy experiments (one ideal and another more realistic) to see their ability to constrain PGWs using the TE power spectrum alone. We found that an ideal experiment, one limited only by cosmic variance, can detect PGWs with a ratio of tensor to scalar metric perturbation power spectra   r = 0.3  at 99.9 per cent confidence level using only the TE correlation. This value is comparable with current constraints obtained by the Wilkinson Microwave Anisotropy Probe based on the 2σ upper limits to the B-mode amplitude. We demonstrate that to measure PGWs by their contribution to the TE cross-correlation power spectrum in a realistic ground-based experiment when real instrumental noise is taken into account, the tensor-to-scalar ratio, r , should be approximately three times larger.     

Kinematics of O-B5 giants

To study the kinematics of O-B5 giant stars (luminosity class III), 290 non-Gould belt stars with proper motions taken from the Hipparcos catalogue are used, of which 107 have radial velocities taken from other sources. Semidefinite programming solves for the kinematical parameters and the coefficients of the velocity ellipsoid. The condition that both solutions must yield the same solar velocity is enforced. The results obtained are reasonable: solar velocity of 13.83 ± 0.17 km s⁻¹; Oort's constants, in units of km s⁻¹ kpc⁻¹, A = 16.08 ± 0.72 and   B =−10.74 ± 0.65,  implying a rotational velocity of 228.0 ± 21.4 km s⁻¹ if we take the distance to the Galactic Centre as 8.5 ± 1.1 kpc; velocity dispersions, in units of km s⁻¹, of  σ _x = 32.44 ± 5.04, σ _y = 26.16 ± 2.75, σ _z = 18.71 ± 2.39  with a vertex deviation of