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.     

The Sunyaev–Zel'dovich effect and Faraday rotation contributions of galaxy groups to the CMB angular power spectrum

The Sunyaev–Zel'dovich (SZ) effect and the Faraday rotation from haloes are examined over a wide mass range, including gas condensation and magnetic field evolution. Contributions to the cosmic microwave background (CMB) angular power spectrum are evaluated for galaxy clusters, galaxy groups and galaxies. Smaller mass haloes are found to play a more important role than massive haloes for the B -mode polarization associated with the SZ CMB anisotropies. The B modes from the Faraday rotation dominate the secondary B modes caused by gravitational lensing at  ℓ > 3000  . Measurement of B -mode polarization in combination with the SZ power spectrum can potentially provide important constraints on intracluster magnetic field and gas evolution at early epochs.     

An SZ temperature decrement–X-ray luminosity relation for galaxy clusters

We present the observed relation between Δ T _SZ, the cosmic microwave background (CMB) temperature decrement due to the Sunyaev–Zeldovich (SZ) effect, and L , the X-ray luminosity of galaxy clusters. We discuss this relation in terms of the cluster properties, and show that the slope of the observed Δ T _SZ– L relation is in agreement with both the L – T _e relation based on numerical simulations and X-ray emission observations, and the M _gas– L relation based on observation. The slope of the Δ T _SZ– L relation is also consistent with the M _tot– L relation, where M _tot is the cluster total mass based on gravitational lensing observations. This agreement may be taken to imply a constant gas mass fraction within galaxy clusters, however, there are large uncertainties, dominated by observational errors, associated with these relations. Using the Δ T _SZ– L relation and the cluster X-ray luminosity function, we evaluate the local cluster contribution to arcmin-scale cosmic microwave background anisotropies. The Compton distortion y -parameter produced by galaxy clusters through the SZ effect is roughly two orders of magnitude lower than the current upper limit based on FIRAS observations.     

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.     

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.     

Statistics of the Sunyaev–Zel'dovich effect power spectrum

Using large numbers of simulations of the microwave sky, incorporating the cosmic microwave background (CMB) and the Sunyaev–Zel'dovich (SZ) effect due to clusters, we investigate the statistics of the power spectrum at microwave frequencies between spherical multipoles of 1000 and 10 000. From these virtual sky maps, we find that the spectrum of the SZ effect has a larger standard deviation by a factor of 3 than would be expected from purely Gaussian realizations, and has a distribution that is significantly skewed towards higher values, especially when small map sizes are used. The standard deviation is also increased by around 10 per cent compared to the trispectrum calculation due to the clustering of galaxy clusters. We also consider the effects of including residual point sources and uncertainties in the gas physics. This has implications for the excess power measured in the CMB power spectrum by the Cosmic Background Imager (CBI) and Berkeley–Illinois–Maryland Association (BIMA) experiments. Our results indicate that the observed excess could be explained using a lower value of σ₈ than previously suggested, however the effect is not enough to match  σ₈= 0.825  . The uncertainties in the gas physics could also play a substantial role. We have made our maps of the SZ effect available online.     

Do Type Ia supernovae prove Λ>0?

The evidence for positive cosmological constant Λ from Type Ia supernovae is re-examined.
Both high redshift supernova teams are found to underestimate the effects of host galaxy extinction. The evidence for an absolute magnitude–decay time relation is much weakened if supernovae not observed before maximum light are excluded. Inclusion of such objects artificially suppresses the scatter about the mean relation.
With a consistent treatment of host galaxy extinction and elimination of supernovae not observed before maximum, the evidence for a positive lambda is not very significant  (3–4 σ )  . A factor which may contribute to apparent faintness of high- z supernovae is evolution of the host galaxy extinction with z .
The Hubble diagram using all high- z distance estimates, including SZ clusters and gravitational lens time-delay estimates, does not appear inconsistent with an  Ω_o=1  model.
Although a positive Λ can provide an (albeit physically unmotivated) resolution of the low curvature implied by cosmic microwave background (CMB) experiments and evidence that  Ω_o<1  from large-scale structure, the direct evidence from Type Ia supernovae seems at present to be inconclusive.     

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.     

Bayesian modelling of clusters of galaxies from multifrequency-pointed Sunyaev–Zel'dovich observations

We present a Bayesian approach to modelling galaxy clusters using multi-frequency pointed observations from telescopes that exploit the Sunyaev–Zel'dovich effect. We use the recently developed multinest technique to explore the high-dimensional parameter spaces and also to calculate the Bayesian evidence. This permits robust parameter estimation as well as model comparison. Tests on simulated Arcminute Microkelvin Imager observations of a cluster, in the presence of primary CMB signal, radio point sources (detected as well as an unresolved background) and receiver noise, show that our algorithm is able to analyse jointly the data from six frequency channels, sample the posterior space of the model and calculate the Bayesian evidence very efficiently on a single processor. We also illustrate the robustness of our detection process by applying it to a field with radio sources and primordial CMB but no cluster, and show that indeed no cluster is identified. The extension of our methodology to the detection and modelling of multiple clusters in multi-frequency SZ survey data will be described in a future work.     

Constraints on cosmological parameters from recent measurements of cosmic microwave background anisotropy

A key prediction of cosmological theories for the origin and evolution of structure in the Universe is the existence of a 'Doppler peak' in the angular power spectrum of cosmic microwave background (CMB) fluctuations. We present new results from a study of recent CMB observations which provide the first strong evidence for the existence of a 'Doppler peak' localized in both angular scale and amplitude. This first estimate of the angular position of the peak is used to place a new direct limit on the curvature of the Universe, corresponding to a density of Ω = 0.7^+0.8_−0.5, consistent with a flat universe. Very low-density 'open' universe models are inconsistent with this limit unless there is a significant contribution from a cosmological constant. For a flat standard cold dark matter dominated universe we use our results in conjunction with big bang nucleosynthesis constraints to determine the value of the Hubble constant as H ₀ = 30 − 70 km s⁻¹ Mpc⁻¹ for baryon fractions Ω_b = 0.05 to 0.2. For H ₀ = 50 km s⁻¹ Mpc⁻¹ we find the primordial spectral index of the fluctuations to be n  = 1.1 ± 0.1, in close agreement with the inflationary prediction of n  ≃ 1.0.     

Limits on isocurvature perturbations from non-Gaussianity in WMAP temperature anisotropy

We study the effect of primordial isocurvature perturbations on non-Gaussian properties of cosmic microwave background (CMB) temperature anisotropies. We consider generic forms of the non-linearity of isocurvature perturbations which can be applied to a wide range of theoretical models. We derive analytical expressions for the bispectrum and the Minkowski Functionals for CMB temperature fluctuations to describe the non-Gaussianity from isocurvature perturbations. We find that the isocurvature non-Gaussianity in the quadratic isocurvature model, where the isocurvature perturbation S is written as a quadratic function of the Gaussian variable  σ,  S =σ²−〈σ²〉  , can give the same signal-to-noise ratio as   f _NL= 30  even if we impose the current observational limit on the fraction of isocurvature perturbations contained in the primordial power spectrum α. We give constraints on isocurvature non-Gaussianity from Minkowski Functionals using the Wilkinson Microwave Anisotropy Probe ( WMAP ) 5-year data. We do not find a significant signal of isocurvature non-Gaussianity. For the quadratic isocurvature model, we obtain a stringent upper limit on the isocurvature fraction  α < 0.070  (95 per cent CL) for a scale-invariant spectrum which is comparable to the limit obtained from the power spectrum.     

Sunyaev–Zel'dovich predictions for the Planck Surveyor satellite using the Hubble Volume simulations

We use the billion-particle Hubble Volume simulations to make statistical predictions for the distribution of galaxy clusters that will be observed by the Planck Surveyor satellite through their effect on the cosmic microwave background – the Sunyaev–Zel'dovich (SZ) effect. We utilize the lightcone data sets for both critical density ( τ CDM) and flat low-density (ΛCDM) cosmologies: a 'full-sky' survey out to z ∼0.5 , two 'octant' data sets out to beyond z =1 , and a 100 square degree data set extending to z ∼4 . Making simple, but robust, assumptions regarding both the thermodynamic state of the gas and the detection of objects against an unresolved background, we present the expected number of SZ sources as a function of redshift and angular size, and also as a function of flux (for both the thermal and kinetic effects) for three of the relevant High Frequency Instrument frequency channels. We confirm the expectation that the Planck Surveyor will detect around 5×10⁴ clusters, though the exact number is sensitive to the choice of several parameters including the baryon fraction, and also to the cluster density profile, so that either cosmology may predict more clusters. We also find that the majority of detected sources should be at z <1.5 , and we estimate that around 1 per cent of clusters will be spatially resolved by the Planck Surveyor , though this has a large uncertainty.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

A study of the Sunyaev–Zel'dovich increment using archival SCUBA data

In a search for evidence of the short wavelength increment in the Sunyaev–Zel'dovich (SZ) effect, we have analysed archival galaxy cluster data from the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope, resulting in the most complete pointed survey of clusters at 850 μm to date. SCUBA's 850-μm passband overlaps the peak of the SZ increment. The sample consists of 44 galaxy clusters in the range 0 < z < 1.3. Maps of each of the clusters have been made and sources have been extracted; as an ancillary product, we generate the most thorough galaxy cluster point source list yet from SCUBA. 17 of these clusters are free of obvious active galactic nuclei (AGN) and have data deep enough to provide interesting measurements of the expected SZ signal. Specialized analysis techniques are employed to extract the SZ effect signal from these SCUBA data, including using SCUBA's short wavelength band as an atmospheric monitor and fitting the long wavelength channel to a model of the spatial distribution of each cluster's SZ effect. By explicitly excising the exact cluster centre from our analysis, we demonstrate that emission from galaxies within the cluster does not contaminate our measurement. The SZ amplitudes from our measurements are consistently higher than the amplitudes inferred from low-frequency measurements of the SZ decrement.     

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.     

Measurement of the Sunyaev–Zel'dovich increment in massive galaxy clusters

We have detected the Sunyaev–Zel'dovich (SZ) increment at 850 μm in two galaxy clusters (Cl 0016+16 and MS 1054.4−0321) using the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. Fits to the isothermal β model yield a central Compton y parameter of  (2.2 ± 0.7) × 10⁻⁴  and a central 850-μm flux of  Δ I ₀= 2.2 ± 0.7 mJy beam⁻¹  in Cl 0016. This can be combined with decrement measurements to infer   y = (2.38 ±^0.36_0.34) × 10⁻⁴  and   v _pec= 400±¹⁹⁰⁰₁₄₀₀ km s⁻¹  . In MS 1054 we find a peak 850-μm flux of  Δ I ₀= 2.0 ± 1.0 mJy beam⁻¹  and   y = (2.0 ± 1.0) × 10⁻⁴  . To be successful such measurements require large chop throws and non-standard data analysis techniques. In particular, the 450-μm data are used to remove atmospheric variations in the 850-μm data. An explicit annular model is fit to the SCUBA difference data in order to extract the radial profile, and separately fit to the model differences to minimize the effect of correlations induced by our scanning strategy. We have demonstrated that with sufficient care, SCUBA can be used to measure the SZ increment in massive, compact galaxy clusters.     

The formation and survival of discs in a ΛCDM universe

We study the formation of galaxies in a Λ cold dark matter (ΛCDM) universe using high-resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight isolated haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at   z = 0  , none of our galaxies contains a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The   z = 0  spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at   z = 0  nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at   z ≳ 2  , regardless of their   z = 0  morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic ΛCDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.     

Neutralino annihilations and the gas temperature in the dark ages

Assuming that the dark matter is entirely made up of neutralinos, we re-visit the role of their annihilation on the temperature of diffuse gas in the high-redshift universe  ( z > 10)  , before the formation of luminous structures. We consider neutralinos of particle mass 36 and 100 GeV. The former is able to produce  ∼7  e ⁻ e ⁺  particles per annihilation through the fremionic channel, and the latter ∼53 particles assuming a purely bosonic channel. High-energy   e ⁻ e ⁺  particles up-scatter the cosmic microwave background (CMB) photons into higher energies via the inverse-Compton scattering. The process produces a power-law   e ⁻ e ⁺  energy spectrum of index −1 in the energy range of interest, independent of the initial energy distribution. The corresponding energy spectrum of the up-scattered photons is a power law of index −1/2, if absorption by the gas is not included. The scattered photons photoheat the gas by releasing electrons which deposit a fraction (14 per cent) of their energy as heat into the ambient medium. For uniformly distributed neutralinos, the heating is insignificant. The effect is greatly enhanced by the clumping of neutralinos into dense haloes. We use a time-dependent clumping model which takes into account the damping of density fluctuations on mass-scales smaller than  ∼10⁻⁶ M_⊙  . With this clumping model, the heating mechanism boosts the gas temperature above that of the CMB after a redshift of   z ∼ 30  . By   z ≈ 10  , the gas temperature is nearly 100 times its temperature when no heating is invoked. Similar increase is obtained for the two neutralino masses considered.     

Large-scale cosmic microwave background anisotropies and dark energy

In this paper we investigate the effects of perturbations in a dark energy component with a constant equation of state on large-scale cosmic microwave background (CMB) anisotropies. The inclusion of perturbations increases the large-scale power. We investigate more speculative dark energy models with   w < −1  and find the opposite behaviour. Overall the inclusion of perturbations in the dark energy component increases the degeneracies. We generalize the parametrization of the dark energy fluctuations to allow for an arbitrary constant sound speed, and we show how constraints from CMB experiments change if this is included. Combining CMB with large-scale structure, Hubble parameter and supernovae observations we obtain   w =−1.02 ± 0.16 (1σ)  as a constraint on the equation of state, which is almost independent of the sound speed chosen. With the presented analysis we find no significant constraint on the constant speed of sound of the dark energy component.     

Large-scale cosmic homogeneity from a multifractal analysis of the PSCz catalogue

We investigate the behaviour of galaxy clustering on large scales using the PSC z catalogue. In particular, we ask whether there is any evidence of large-scale fractal behaviour in this catalogue. We find the correlation dimension in this survey varies with scale, consistent with other analyses. For example, our results on small and intermediate scales are consistent those obtained from the QDOT sample, but the larger PSC z sample allows us to extend the analysis out to much larger scales. We find firm evidence that the sample becomes homogeneous at large scales; the correlation dimension of the sample is D ₂=2.992±0.003 for r >30  h ⁻¹ Mpc. This provides strong evidence in favour of a universe that obeys the cosmological principle.