Cosmological constraints from the observed angular cross-power spectrum between Sunyaev-Zel’dovich and X-ray surveys
Key Words.:large-scale structures, galaxy clusters, cosmic microwave background, intracluster medium
We present the first detection of the cross-correlation angular power spectrum between the thermal Sunyaev-Zel’dovich (tSZ) effect and the X-ray emission over the full sky. The tSZ effect and X-rays are produced by the same hot gas within groups and clusters of galaxies, which creates a naturally strong correlation between them that can be used to boost the joint signal and derive cosmological parameters. We computed the correlation between the ROSAT All Sky Survey in the 0.5-2 keV energy band and the tSZ effect reconstructed from six Planck all-sky frequency maps between 70 and 545 GHz. We detect a significant correlation over a wide range of angular scales. In the range , the cross-correlation of X-rays to tSZ is detected at an overall significance of 28 . As part of our systematic study, we performed a multi-frequency modelling of the AGN contamination and the correlation between cosmic infra-red background and X-rays. Taking advantage of the strong dependence of the cross-correlation signal on the amplitude of the power spectrum, we constrained , where modelling uncertainties dominate statistical and systematic uncertainties. We also derived constraints on the mass indices of scaling relations between the halo mass and X-ray luminosity, L-M, and SZ signal, Y-M, , and on the indices of the extra-redshift evolution, .
Galaxy clusters are the largest virialized structures in the Universe. They are excellent tracers of the matter distribution, and their
abundance can be used to constrain the cosmological model in an independent way.
Galaxy clusters can be identified in the optical bands as concentrations of galaxies (see e.g. Abell et al., 1989; Gladders & Yee, 2005; Koester et al., 2007; Rykoff et al., 2013), they can be observed in X-rays by the bremsstrahlung emission
produced by the ionized intra-cluster medium (ICM) (see e.g. Bohringer et al., 2000; Ebeling et al., 2000, 2001; Böhringer et al., 2001). The same hot ICM also creates a distortion in the black-body spectrum of the cosmic
microwave background (CMB) through the thermal Sunyaev-Zel’dovich (tSZ) effect (Sunyaev & Zeldovich, 1969, 1972), an inverse-Compton scattering between the CMB photons and the ionized electrons in the ICM.
Recent catalogues based on measurements of the tSZ have been produced from Planck (Planck Collaboration early results. VIII, 2011; Planck Collaboration results. XXIX, 2014), ACT (Marriage et al., 2011), and SPT (Bleem et al., 2014) data.
The expected number of galaxy clusters is extremely sensitive to cosmological parameters, especially to the normalization of the matter power spectrum, or its fluctuations on 8 Mpc scales, . Galaxy-cluster surveys have been used to constrain the cosmological parameters to a high level of precision (The Dark Energy Survey Collaboration, 2005; Vanderlinde et al., 2010; Sehgal et al., 2011; Böhringer et al., 2014; Planck Collaboration results. XX, 2014).
Recent tSZ surveys, spanning wide areas in the sky, enable measuring the angular power spectrum of the tSZ signal over a wide range of multipoles and using it to constrain the cosmological model (Planck Collaboration results. XXI, 2014).
The present constraints on cosmological parameters from tSZ power spectrum are limited by our knowledge of the relation between the total mass and the tSZ flux in galaxy
clusters (Planck Collaboration early results. XI, 2011; Hasselfield et al., 2013; Benson et al., 2013).
An important factor that limits the accuracy in the determination of cosmological parameters is the contamination of the tSZ by other astrophysical emissions, mainly radio, infra-red point sources,
and cosmic infra-red background (Dunkley et al., 2011; Shirokoff et al., 2011; Reichardt et al., 2012; Sievers et al., 2013; Planck Collaboration results. XXI, 2014).
Performing a power spectrum analysis in
the X-ray domain is difficult because the X-ray power spectrum is dominated by the signal from active galactic nuclei (AGN).
Previous works have only been able to set upper limits on the cosmological model from the X-ray power spectrum alone (Diego et al., 2003).
Since the X-ray emission and tSZ effect are produced by the same hot gas in galaxy clusters, we expect a strong correlation signal between tSZ and X-ray surveys. The correlation between tSZ and X-rays may reduce the contamination and systematic effects that affect each sample separately and set tighter limits on the cosmological model.
Recent studies (see e.g. Hurier et al., 2014) have shown that the higher sensitivity and resolution of the Planck mission will probably allow us to detect a significant cross-correlation of X-rays and tSZ.
The correlation between the tSZ signal and the X-ray emission has already been used in several studies. Statistical measurements of the tSZ effect have been obtained from WMAP
data by stacking the temperature anisotropies toward known X-ray clusters (Afshordi et al., 2005; Lieu et al., 2006; Afshordi et al., 2007; Atrio-Barandela et al., 2008; Diego & Partridge, 2010; Komatsu et al., 2011; Melin et al., 2011; Planck Collaboration early results. X, 2011) or by computing angular cross-power
spectra between cluster catalogues and temperature anisotropy maps from WMAP (Fosalba et al., 2003; Bennett et al., 2003; Myers et al., 2004; Hernández-Monteagudo et al., 2004, 2006; Hinshaw et al., 2007) and more recently from Planck (Hajian et al., 2013).
Attempts have also been made to directly cross-correlate WMAP temperature maps with ROSAT All Sky Survey (RASS) (Diego et al., 2003), but this was limited by the sensitivity and resolution of the WMAP experiment.
The paper is organized as follows: In Sect. 2 we detail our modelling procedure for the cross-correlation of X-rays and tSZ. In Sect. 3 we present the Planck and ROSAT data we used for this analysis. In Sect. 4 we present our methodology for extracting the tSZ-X angular cross-correlation power spectrum and show the first significant detection of the tSZ-X cross-correlation power spectrum from tSZ and X-ray full-sky data. In Sect. 5.1 we estimate the uncertainty levels produced by data and modelling limitations. In Sect. 6 we show the constraints for cosmological and scaling law parameters. We discuss our results in Sect.7.
2 Modelling the tSZ-Xray cross-correlation
2.1 Thermal Sunyaev-Zel’dovich efffect
The tSZ effect produces a small spectral distortion in the black-body spectrum of the CMB (Sunyaev & Zeldovich, 1969, 1972). Its intensity is related to the integral of the pressure along the line of sight, or more specifically, the Compton parameter, which is defined as
where is the distance element along the line of sight
and and are the electron number density and the temperature.
In units of CMB temperature, the contribution of the tSZ effect to the sub-millimeter sky intensity for a given observation frequency is given by
Neglecting relativistic corrections, we have , with . This function is equal to 0 around GHz, it is negative at lower frequencies and positive for higher frequencies. The characteristic signal of the spectral distortion from the ICM can be effectively used to directly isolate the projected pressure distribution in galaxy clusters from multi-frequency maps.
2.2 X-ray emission from galaxy clusters
The ionized gas in the intra-cluster medium produces an X-ray emission via bremsstrahlung. This radiation is proportional to the square of the electron density. The energy spectrum of the X-ray emission from a galaxy cluster mainly depends on the temperature, , of the ICM and to a lesser extent on the metallicity, , of the gas. From an observational point of view, the X-ray spectrum also depends on the redshift, . The observed count rate at low energy depends on the column density of neutral hydrogen, , on the line of sight. In this work, we model the galaxy clusters emission using a metal model (Mewe et al., 1985). We refer to Hurier et al. (2014) for a more detailed description of the X-ray emission modelling.
2.3 Cross-correlation between tSZ effect and X-ray emission
The angular cross-power spectrum of tSZ effect and X-ray count map reads
with and the coefficients from the spherical harmonics decomposition of tSZ map and X-ray count map. This equation is applied without the detector transfer function and only for the full sky. Nevertheless, when computed for partial regions of the sky, we take the mask and beam effects into account, see Sect. 4. To model this cross-correlation and the auto-correlation power spectra, we assume the following general expression
where is the Poissonian contribution and is the two-halo term.
The Poissonian term can be computed by assuming the square of the Fourier transform of normalized tSZ and X-ray projected profiles, weighted by the mass function and the respective tSZ effect flux and X-ray count-rate of galaxy clusters (see e.g. Cole & Kaiser, 1988; Komatsu & Seljak, 2002, for a derivation of the tSZ auto-correlation angular power spectrum),
where and are the average X-ray count-rate and tSZ flux that depends on and . They are given by scaling relations (see Hurier et al., 2014). is the mass function of the
dark matter halo (we considered here the fitting formula of Tinker et al. (2008)), and is
the co-moving volume element.
The factor accounts for the bias produced by the scatter of scaling relations (see Hurier et al., 2014).
The Fourier transform of a 3D profile projected across the line of sight on the sphere reads , where is either tSZ or X-ray 3D profile, , , is the scale radius of the profile.
The two-halo term corresponds to large-scale fluctuations of the dark matter field that induce correlations in the cluster distribution across the sky. It can be computed as (see e.g. Komatsu & Kitayama, 1999; Diego & Majumdar, 2004; Taburet et al., 2011)
with the time dependent linear bias that relates the matter power spectrum, , to the power spectrum of the cluster distribution across the sky. Following Mo & White (1996) and Komatsu & Kitayama (1999), we adopt
with , is the linear growth factor and is the over-density threshold for spherical collapse.
We stress that the two-halo term is negligible for tSZ-X-ray cross-correlation purposes.
The amplitude of the spectrum follows the dependencies (Hurier et al., 2014)
where , , are the normalizations of the - and - relations, , and are the power-law indices ( and ), is the bias between the mass estimated from X-ray measurements and the true matter mass of galaxy clusters. In the following we consider (Planck Collaboration results. XX, 2014, and references therein). and is the normalization of the mass function, for which we consider an uncertainty of 10%.
We list the value and uncertainty on parameters of the - and - scaling relations in Table. 1. We define
as a single parameter that accounts for the overall amplitude
of the tSZ-X power spectrum.
We used the Planck nominal mission dataset (Planck Collaboration results. I, 2014) available at Planck Legacy Archive (PLA111http://www.sciops.esa.int).
We considered frequencies from 70 to 857 GHz. The two lowest frequency channels, at 30 and 44 GHz, were not considered because
their angular resolution is too poor.
We assumed that the Planck beams can be well approximated by circular Gaussian beams. We considered FWHM values from Planck Collaboration results. VII (2014).
For the tSZ transmission in Planck spectral bandpasses, we used the values provided in Planck Collaboration results. IX (2014).
We used the ROSAT all-sky survey (RASS) public data222ftp://ftp.xray.mpe.mpg.de/rosat/archive/, which cover 99.8% of the sky, including 97% that have an exposure time longer than 100s (Voges et al., 1999).
X-ray photons with an energy below 0.5 keV were not considered to reduce the impact of absorption (considering the
photoelectric cross-section from Morrison & McCammon (1983) and the typical value from Kalberla et al. (2005)).
Then, we constructed a full-sky map of the photon count rate in the energy range 0.5-2.0 keV from each ROSAT photon event file and exposure map.
We projected each event over the sky using HEALPix (Górski et al., 2005) pixelisation scheme at a resolution of .
The ROSAT exposure maps were reprojected using a nearest-neighbour interpolation on an HEALPix grid with .
Thus, the reprojected full-sky RASS has a resolution of 1.7 arcmin (size of the HEALPix pixels for ). Below, we account for the loss of power produced by the convolution of the RASS data with the pixel window function.
We have checked that reprojecting ROSAT number count images provides similar results for our analysis as would be achieved with ROSAT photon event files.
We also constructed a mask to avoid contamination by galactic or point source emissions in Planck data.
Considering that the Planck 857 GHz channel is a good tracer of the thermal dust emission, we chose to mask all regions that present an emission above 3 333 is defined as the unit in which a black body spectrum at 2.725 K is flat with respect to the frequency. at this frequency (see Planck Collaboration results. IX, 2014, for unit convention and conversion).
We masked all sources from the Planck compact source catalogue (Planck Collaboration results. XXVIII, 2014) detected in at least one frequency with a signal-to-noise
ratio (S/N) above 5.
We also masked all regions with an exposure of 0 seconds in the RASS survey. After applying these cuts, we kept about 60% of the sky for the analysis. We present the resulting mask in Fig. 1. We also considered masks of 20% and 40% for robustness checks with different thresholds for the thermal dust emission.
We chose to exclude the Planck official CMB and foreground maps from the cleaning of frequency maps. This choice was motivated by the contamination of the Planck official CMB maps by tSZ residuals (see e.g. Bobin et al., 2014), which would bias our analysis. We also show in Sect. 4.1 that galactic foregrounds can be kept under control using thermal dust cleaning and a galactic mask. As a consequence, we did not perform CO (Planck Collaboration results. XIII, 2013) or low-frequency galactic emission (Planck Collaboration results. XII, 2014) cleaning, because it would have increased the noise level without significantly improving the results.
4 Dectecting the tSZ-Xray cross-correlation
The tSZ-X angular cross-power spectrum can be estimated from frequency maps or from a tSZ Compton parameter map.
In this section, we explore the two approaches and discuss advantages and drawbacks of each of them.
In each case, we binned the power spectra and corrected the cross-spectra for beam and mask effects. The beam was corrected for by dividing the power-spectra by the beam transfer function in space. We deconvolved the power spectra by the mask-induced mixing matrix, which also accounts for the covered sky fraction. For a detailed description of the mask correction see Tristram et al. (2005).
4.1 Analysis for each frequency channel
We computed the angular cross-power spectrum, , between Planck frequency maps from 70 to 857 GHz and the RASS full-sky map in the 0.5-2.0 keV energy band.
At a given frequency, the sky signal is dominated by CMB below 217 GHz and by thermal dust emission above 217 GHz.
We first corrected for CMB contamination by applying the transformation
where is the tSZ transmission in the Planck map at the frequency, .
The term is very small because the tSZ transmission in the Planck 217 GHz channel is almost null.
This cleaning process prevents contamination at calibration uncertainty level (Planck Collaboration results. VIII, 2014) by processes that follow the same spectral energy distribution (SED) as the CMB black-body radiation (e.g. the kinetic SZ effect).
Then we corrected for thermal dust contamination by applying the transformation
The value of was estimated from a linear fit of the lowest multipoles of the power spectra of ()
because the thermal dust emission dominates the cross-power spectra at low multipoles.
The factor is negligible because of the ratio between tSZ and thermal dust SEDs at high frequencies.
spectra are referred to as the cleaned cross-power spectra in the following.
In Fig. 2, we present the cleaned cross-power spectra between the RASS map and the Planck maps at 70, 100, 143, 353, and 545 GHz.
We computed the of each spectrum compared to a null correlation. We obtain , , , , and for 31 degrees of freedom from to at 70, 100, 143, 353, and 545 GHz, respectively. We also display the cross-correlation at 217 GHz without CMB and dust subtraction, and with dust correction alone. After dust correction, the cross-spectrum at 217 GHz is compatible with zero, with for 31 degree of freedom.
We compared these cross-spectra with the theoretical prediction assuming the Planck-CMB best-fitting CDM cosmology (, , and km/s/Mpc) (Planck Collaboration results. XVI, 2014). We observe a significant correlation that follows a tSZ SED. However, we note that the correlation signal appears to be significantly weaker than the prediction. We also note a significant contamination in the cross-power spectra at 100 and 143 GHz. This contamination can be seen at high multi-pole values ( > 1000) and is interpreted as a correlation between radio and X-ray AGNs that compensates the negative correlation between the tSZ and the X-ray emission from galaxy clusters at those frequencies.
The angular cross-power spectrum at 545 GHz significantly exceeds the theoretical prediction. This excess can be interpreted as a contamination by the correlation between X-ray emission and the cosmic infra-red background.
Figure 2 illustrates the main advantages of using a multi-frequency measurement of the tSZ-X cross-power spectrum. Indeed, with multi-frequency measurements, we can easily distinguish contamination sources from the tSZ-X signal by using the differences in SEDs.
4.2 Compton parameter map based analysis
We used a second approach to extract the tSZ-X cross-power spectrum.
In this case, we built a tSZ Compton parameter map (-map) using the MILCA method (Hurier et al., 2013) on the Planck maps from at 100 to 857 GHz.
We verified that including frequencies from 30 to 70 GHz does not change the results.
The tSZ-X angular cross-power spectrum is obtained by directly cross-correlating the reprojected RASS full-sky map and the -map, using 69% of the sky.
In Fig. 3, we present the derived cross-power spectrum and compare it with the theoretical prediction assuming the Planck-CMB best-fitting cosmology.
We also computed the of the spectrum with respect to a null correlation. We derive for 33 degrees of freedom from to .
Similarly to the multi-frequency analysis, we observe that the measured power spectrum is significantly weaker than the prediction.
However, the shape of the measured power spectrum agrees with the shape of the prediction.
As the -map reconstruction is performed both in pixel and frequency domains, it allows us to extract the tSZ signal at a higher S/N than a linear combination of power spectra.
However, the identification and estimate of contamination sources is more complicated.
Consequently, the two approaches are complementary and were used together to check the robustness of the results.
4.3 Robustness of the detection
In this section we verify that our signal is produced by galaxy cluster and not by AGNs.
We estimated the contamination from known AGN using X-ray catalogues such as the ROSAT Bright Survey (RBS, Fischer et al., 1998).
First, we projected the RBS sources on a full-sky map.
We computed the cross-spectrum between the -map and the projected RBS map considering only sources flagged as AGNs. Then, we performed the same analysis for sources flagged as clusters.
Uncertainties were computed following the approach presented in Sect. 5.1.
We present the obtained cross-spectra in Fig. 4. We do not observe any significant contribution from this AGNs sample. These spectra illustrate that our signal is dominated by X-ray emission from galaxy clusters and not from X-ray AGNs.
We also estimated the contribution of known tSZ and X-ray galaxy clusters to the tSZ-X cross-power spectrum. To do so, we built several tSZ maps for different subsamples of galaxy clusters.
We considered the Planck SZ catalogue (Planck Collaboration early results. VIII, 2011; Planck Collaboration results. XXIX, 2014), the Planck cosmo sample (Planck Collaboration results. XX, 2014), and the MCXC
(Piffaretti et al., 2011).
For the Planck SZ catalogue we used confirmed galaxy clusters (861 clusters) with fluxes and radius estimated from 2D likelihoods provided in Planck Collaboration results. XXIX (2014) and a universal pressure profile (Arnaud et al., 2010). For X-rays clusters, we predicted the tSZ flux assuming a scaling relation between and from Planck Collaboration early results. XI (2011) and X-ray deduced values for .
We projected each cluster on an oversampled grid with a pixel size of 0.1 (e.g. drizzling) to avoid flux loss during the projection. Then we convolved the oversampled map with a beam FWHM of 10 arc minutes. We reprojected the oversampled map on a HEALpix full-sky map with 1.7 arc minute pixels using a nearest-neighbour interpolation.
Then, we computed the angular cross-power spectrum between the tSZ template maps and the ROSAT full-sky count-rate map. Uncertainties were computed following the approach presented in in Sect. 5.1.
In Fig. 5, we present the derived angular power spectra for each sample of clusters. The known X-ray clusters from the MCXC contribute about 75% of the total power at low- (). Confirmed Planck tSZ clusters contribute about 50% of the power of the tSZ-X spectrum at low- (). The Planck cosmo sample contributes about 1% at low-.
Known clusters contribute essentially at low-. At high-, the tSZ-X spectrum presents contribution from low-mass and high- galaxy clusters. Planck -confirmed clusters and MCXC clusters only contribute 20% and 40% of the total power at , as expected considering the contribution from undetected galaxy clusters.
These spectra demonstrate that the full tSZ-X cross-power spectrum contains additional informations on cosmology relevant to tSZ number count studies.
5 Estimating bias and uncertainties
We discuss in detail the main sources of uncertainties and biases in the measurement of the tSZ-X cross-power spectrum.
5.1 Statistical uncertainties
The cleaned angular cross-power spectra are constructed using a linear combination of three cross-power spectra, , at different frequencies, see Sect. 4.1. Consequently, we need to estimate variances and covariances of these spectra. The variance, , of can be expressed as
and the covariance, , between and reads
Then, the variance of can be computed as
It can be useful to combine constraints from the different cleaned power spectra. To do so, we computed the covariance between and as
We propagated the whole covariance matrix through binning and mask deconvolution processes.
In Fig. 6, we present the uncertainties on and the correlation factor between and for the cleaned cross-power spectra at 70, 100, 143, 353, and 545 GHz.
In this figure and in the diagonal plots, we observe that at low frequencies the uncertainty level increases with the multipole. This is produced by the deconvolution from the instrumental resolution. We observe a flattening in the uncertainty level at high- for the highest frequencies. This flattening indicates that the signal is dominated by the instrumental noise.
At low- the correlation factors are close to 1 or -1, indicating an almost total correlation between uncertainties below at 70, 100, 143, and 353 GHz. This indicates that the uncertainties at these multipole values are dominated by cosmic variance of the tSZ-X cross-spectrum.
The instrumental noise domination at high- also appears in the correlation factor between frequencies. When the noise dominates the uncertainty budget, these coefficients become zero.
For a tSZ-X power spectrum estimated from a -map, the variance, , can be directly estimated as
The MILCA method is tailored to minimize instrumental noise and the variance of other astrophysical emissions, therefore we derive lower uncertainties (Fig. 3) than using the previous approach (Fig. 2).
We also probate these uncertainty through the binning and mask deconvolution processes.
5.2 Bias sources
There are several possible sources of bias in the tSZ-X cross-power spectrum.
The thermal dust emission from the Milky Way can produce bias
because it is correlated with the absorption of the X-rays.
By masking the galactic plane and considering X-rays in the range [0.5,2.0] keV, we ensure that this bias is under control (see Morrison & McCammon, 1983, for absorption cross-sections).
We have tested several galactic cuts (20, 40, and 60%) and found no significant variations in our measure (below 0.2 on the tSZ-X cross-power spectra). This bias can thus be safely neglected.
There are also sources of biases that are correlated with the cluster spatial distribution over the sky.
Indeed, all extra-galactic emissions are correlated with the matter distribution.
On the X-ray side, we have mainly two components, the X-ray emission from clusters and from AGNs (noted AGN hereafter).
On the microwave side, we have the tSZ effect, the radio emission from AGNs (noted AGN hereafter), and the cosmic infra-red background (CIB).
Consequently, we have five cross-spectra from extra-galactic objects that can bias the measurement.
The extrapolated radio emission from AGNs is weak at Planck frequencies (e.g. Planck Collaboration results. XXIX, 2014), thus it does not produce a significant bias. In addition, the small number of un-masked high-flux X-ray AGNs prevents a significant contamination in terms of tSZ-AGN correlation. As a consequence, we considered that tSZ-AGN, AGN-X, and CIB-AGN cross-correlations are weaker than the AGN-AGN correlations. This leaves the tSZ-X, AGN-AGN, and CIB-X correlations as potentially significant contributions for the microwave-X correlation.
For the multi-frequency estimation of the tSZ-X correlation (see Sect. 4.1), this sources of bias can be separated from the tSZ-X signal using their SED.
However, this is not possible in the -map approach (see Sect. 4.2).
To estimate the AGN-AGN level of contamination in the tSZ-X angular cross-power obtained from the reprojected RASS full sky map and the -map,
we assumed that radio catalogues are accurate tracers of the AGN population that may bias our measurement.
The NRAO VLA Sky Survey (NVSS) (Condon et al., 1998) is a 1.4 GHz continuum
survey covering the entire sky north of Dec. The
associated catalogue of
discrete sources contains over 1.8 million radio sources. South of
Dec and at galactic latitudes ,
the Sydney University Molonglo Sky Survey (SUMSS) (Mauch et al., 2003, 2008)
is a 843 MHz continuum survey also providing a radio source catalogue.
SUMSS and NVSS have similar sensitivities and angular resolutions, and
combined, they cover the whole sky.
We extrapolated the SUMSS sources at 1.4 GHz assuming a spectral index of -1 in intensity units. Then, we built a map of NVSS sources weighted by their flux and filled the region Dec with extrapolated SUMSS sources. We computed the cross-power spectra between the radio-sources map and the -map, , and the RASS full-sky map, .
Thus, we computed the estimate of the AGN-AGN contamination in the spectrum as
In Fig. 7 we present the derived spectrum. We fitted with a constant considering that AGN are point-like sources. Neglecting their clustering444AGN clustering can be safely neglected, the AGN power over the sky is dominated by a small number of objects. Similarly to what we observed for galaxy clusters, the clustering term is smaller
than the Poissonian term., we derive cts.s.arcmin.sr (red line in Fig. 7).
Finally, we corrected the -RASS cross-spectrum for this bias and propagated the related uncertainties.
|X-rays 70-545 GHz||Eq. 18||0.804||0.805||0.006||0.800 - 0.809||0.001||1.22|
|X-rays 70 GHz||Eq. 17||0.782||0.784||0.015||0.770 - 0.798||0.020||1.51|
|X-rays 100 GHz||Eq. 17||0.798||0.799||0.007||0.793 - 0.805||0.010||0.94|
|X-rays 143 GHz||Eq. 17||0.799||0.799||0.006||0.794 - 0.804||0.010||1.24|
|X-rays 353 GHz||Eq. 17||0.815||0.816||0.008||0.809 - 0.823||0.015||0.61|
|X-rays 70 GHz||Eq. 18||0.791||0.790||0.020||0.771 - 0.809||0.001||1.51|
|X-rays 100 GHz||Eq. 18||0.815||0.816||0.010||0.807 - 0.824||0.001||0.85|
|X-rays 143 GHz||Eq. 18||0.806||0.806||0.007||0.800 - 0.812||0.001||1.20|
|X-rays 353 GHz||Eq. 18||0.808||0.809||0.008||0.801 - 0.816||0.001||0.61|
|X-rays 545 GHz||Eq. 18||0.803||0.808||0.018||0.790 - 0.822||0.001||1.35|
|X-rays -map||Eq. 19||0.804||0.804||0.003||0.801 - 0.806||0.002||1.15|
We explored three approaches to set cosmological constraints using the measurement of the tSZ-X cross-power spectrum. First, we considered individual frequencies, neglecting all sources of bias. Then, we considered the constraint from each frequency considering a multi-frequency adjustment for the biases. Finally, we considered the cross-spectrum with the tSZ -map. In the following, we fit the data considering from 50 to 2000. We computed the likelihood functions assuming Gaussian uncertainties over a grid for , , and , considering fixed values for other parameters. Then, we marginalized the likelihood to express it as a function of the degeneracy relation between parameters, . For cosmological constraints, we used fixed values of and and propagated the uncertainties on these values to . In the following, uncertainties are given for 68% confidence level. We verified that we derived compatible results within the error bars using different galactic masks (see Sect. 3). This demonstrates that our results are not significantly affected by galactic foreground contamination.
6.1 Cosmological constraints per frequency channel
Considering a single-frequency approach, we estimated for each frequency individually, excluding 545 GHz, which is contaminated. We assumed the model
The derived likelihoods are presented in Fig. 8 and the best-fitting values are summarized in Table. 2.
Our best-fitting values for increase with frequency
with this simple modelling.
This behaviour is produced by the contamination of AGN and CIB. The radio-loud AGNs produce an excess of correlation at low frequency that compensates for the anti-correlation between tSZ and X-ray emission.
At high frequency, the CIB-X correlation produces an excess of correlation that biases the correlation between tSZ and X-ray emission.
The values deduced from spectra at 100 GHz agree well with those
at 143 GHz, which present the lowest bias level.
6.2 Cosmological constraints from the multi-frequency approach
To prevent contamination produced by AGNs and CIB in the fit, we modelled the measured spectra, , as follows:
where (Gispert et al., 2000) and (Planck Collaboration Int. results. VII, 2013) are fixed SED for CIB and radio-loud AGN contaminations.
Given the large uncertainties at 545 GHz for the tSZ-X cross-correlation measurement, the signal is not sensitive to the shape of the CIB-X correlation. Therefore and for simplicity, we assumed that the CIB-X angular cross-power spectrum has the same shape, with respect to , as the tSZ-X cross-correlation.
We fitted for , and the amplitude of the radio and CIB contamination and the tSZ-X angular cross-power spectrum degenerate amplitude defined in Eq. 8. Parameters and were fitted using the five frequencies 70, 100, 143, 353, and 545 GHz; was fitted individually for each frequency.
The derived likelihoods are presented in Fig. 8. Best-fitting values are summarized in Table. 2. The values now agree well, including the value deduced from the 545 GHz spectra. There is no significant frequency dependence for , indicating that our modelling accounts properly for contamination by other astrophysical components.
Then, we adjusted using the model presented in Eq. 18 and all frequencies simultaneously.
Figure 9 presents the likelihood function for and . We obtain a best-fit value of , with statistical and systematic uncertainties, respectively, where systematic uncertainties account for residual contamination.
We estimated the expected value for considering the cosmological constraints from analyses based on CMB (Planck Collaboration results. XVI, 2014) and cluster counts (Planck Collaboration results. XX, 2014). Cluster number counts and CMB power spectra present different degeneracies for cosmological parameters from those of tSZ-X cross-spectrum. This leads to large error bars when is estimated from cluster counts or CMB. Uncertainty levels derived assuming cluster number count and CMB cosmology are described in Hurier et al. (2014).
We stress that these uncertainties are not representative of the uncertainties on . We need to propagate the uncertainty on the modelling parameters to . Consequently, the constraints on cosmological parameters read with statistical, systematic, and modeling uncertainties, respectively. The constraint on gives We observe that the error budget is dominated by modelling uncertainties.
6.3 Cosmological constraints from a -map
We estimated the level of CIB-X contamination in the -RASS power spectrum by propagating the amplitude through the MILCA linear weights. We deduce that this contamination is negligible. Indeed, the CIB only contaminates the tSZ-X correlation at high frequency, where the weights are smaller than at low frequencies.
To reproduce the -RASS spectrum, we assumed the following modedling:
where is estimated as presented in Sect. 5.2.
From this measurement, we obtain and .
In Fig. 9, we present the likelihood for and . The constraints on derived from the tSZ-X power spectrum are consistent with predictions from cluster counts and CMB. Moreover, the constraints derived from the -map are consistent with constraints derived from the multi-frequency analysis presented in Sect. 6.2.
6.4 Constraints on scaling-law parameters
In Fig. 10, we present the constraints on the scaling relation power-law indices from combined tSZ-X and tSZ-tSZ power spectra.
We assumed a fixed pressure profile (best-fit on a GNFW profile from Arnaud et al., 2010) and polytropic index of 1.5 for galaxy clusters.
From the tSZ-Xray cross-spectrum, we derive , which is to be compared with when considering values listed in Table. 1. Figure 10 shows that constraints from the tSZ-X spectrum alone are consistent with previous constraints (see e.g. Pratt et al., 2009; Planck Collaboration early results. XI, 2011) on and .
To explore the possible evolution of the scaling relations with redshift, we assumed a power-law dependence of the form . Then, we fitted for the tSZ and X-ray evolution indices. We derive from the tSZ-X angular cross-power spectrum.
We stress that the constraints on scaling-law indices are degenerate with the assumed density and temperature profiles for galaxy clusters (as discussed in Hurier et al., 2014).
We have performed the first detection of the tSZ-X angular cross-power spectrum directly estimated from X-ray count rate and tSZ signal.
We detected a clearly correlated emission between tSZ-map and X-ray count-rate map in the range at an overall significance of 28 .
We accounted for statistical systematics from residuals (mainly AGN and CIB) and for modelling uncertainties.
Using the amplitude of the tSZ-X angular cross-power spectrum, we derived constraints on cosmological parameters with statistical, systematic, and modelling uncertainties.
These constraints read if we consider from Planck-CMB best-fit cosmology (Planck Collaboration results. XVI, 2014).
We tested the robustness of this result using an analysis based on multiple frequencies and on a tSZ-map.
With a mass bias of 0.2, this result is compatible with constraints based on cluster count and CMB on from Planck Collaboration results. XX (2014) and Planck Collaboration results. XVI (2014). This constraint is also consistent with the cross-correlation measured by Hajian et al. (2013) between the tSZ emission and an X-ray sample of galaxy clusters and with the X-ray cluster number count (Vikhlinin et al., 2009; Böhringer et al., 2014).
Using the shape of the tSZ-X cross-spectrum, we were able to constrain the scaling-law mass and evolution indices of and .
We derived and .
Consistently with previous results from Reichert et al. (2011) for the X-ray scaling relation, we found that the redshift evolution of the scaling law is consistent with self-similar evolution.
The authors thanks A.beelen, M.Arnaud and G.Pratt for useful discussions. We acknowledge the support of the French Agence Nationale de la Recherche under grant ANR-11-BD56-015. This research has made use of the ROSAT all-sky survey data which have been processed at MPE. The development of Planck has been supported by: ESA; CNES and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MICINN and JA (Spain); Tekes, AoF and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); and The development of Planck has been supported by: ESA; CNES and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MICINN and JA (Spain); Tekes, AoF and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); and PRACE (EU).
- Abell et al. (1989) Abell, G. O., Corwin, Jr., H. G., & Olowin, R. P. 1989, ApJS, 70, 1
- Afshordi et al. (2007) Afshordi, N., Lin, Y.-T., Nagai, D., & Sanderson, A. J. R. 2007, MNRAS, 378, 293
- Afshordi et al. (2005) Afshordi, N., Lin, Y.-T., & Sanderson, A. J. R. 2005, ApJ, 629, 1
- Arnaud et al. (2010) Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92
- Atrio-Barandela et al. (2008) Atrio-Barandela, F., Kashlinsky, A., Kocevski, D., & Ebeling, H. 2008, ApJ, 675, L57
- Bennett et al. (2003) Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2003, ApJS, 148, 97
- Benson et al. (2013) Benson, B. A., de Haan, T., Dudley, J. P., et al. 2013, ApJ, 763, 147
- Bleem et al. (2014) Bleem, L. E., Stalder, B., de Haan, T., et al. 2014, ArXiv e-prints
- Bobin et al. (2014) Bobin, J., Sureau, F., Starck, J.-L., Rassat, A., & Paykari, P. 2014, A&A, 563, A105
- Böhringer et al. (2014) Böhringer, H., Chon, G., & Collins, C. A. 2014, ArXiv e-prints
- Böhringer et al. (2001) Böhringer, H., Schuecker, P., Guzzo, L., et al. 2001, A&A, 369, 826
- Bohringer et al. (2000) Bohringer, H., Voges, W., Huchra, J. P., et al. 2000, VizieR Online Data Catalog, 212, 90435
- Cole & Kaiser (1988) Cole, S. & Kaiser, N. 1988, MNRAS, 233, 637
- Condon et al. (1998) Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693
- Diego & Majumdar (2004) Diego, J. M. & Majumdar, S. 2004, MNRAS, 352, 993
- Diego & Partridge (2010) Diego, J. M. & Partridge, B. 2010, MNRAS, 402, 1179
- Diego et al. (2003) Diego, J. M., Silk, J., & Sliwa, W. 2003, MNRAS, 346, 940
- Dunkley et al. (2011) Dunkley, J., Hlozek, R., Sievers, J., et al. 2011, ApJ, 739, 52
- Ebeling et al. (2000) Ebeling, H., Edge, A. C., Allen, S. W., et al. 2000, VizieR Online Data Catalog, 731, 80333
- Ebeling et al. (2001) Ebeling, H., Edge, A. C., & Henry, J. P. 2001, ApJ, 553, 668
- Fischer et al. (1998) Fischer, J.-U., Hasinger, G., Schwope, A. D., et al. 1998, Astronomische Nachrichten, 319, 347
- Fosalba et al. (2003) Fosalba, P., Gaztañaga, E., & Castander, F. J. 2003, ApJ, 597, L89
- Gispert et al. (2000) Gispert, R., Lagache, G., & Puget, J. L. 2000, A&A, 360, 1
- Gladders & Yee (2005) Gladders, M. D. & Yee, H. K. C. 2005, ApJS, 157, 1
- Górski et al. (2005) Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759
- Hajian et al. (2013) Hajian, A., Battaglia, N., Spergel, D. N., et al. 2013, ArXiv e-prints
- Hasselfield et al. (2013) Hasselfield, M., Hilton, M., Marriage, T. A., et al. 2013, J. Cosmology Astropart. Phys., 7, 8
- Hernández-Monteagudo et al. (2004) Hernández-Monteagudo, C., Genova-Santos, R., & Atrio-Barandela, F. 2004, ApJ, 613, L89
- Hernández-Monteagudo et al. (2006) Hernández-Monteagudo, C., Macías-Pérez, J. F., Tristram, M., & Désert, F.-X. 2006, A&A, 449, 41
- Hinshaw et al. (2007) Hinshaw, G., Nolta, M. R., Bennett, C. L., et al. 2007, ApJS, 170, 288
- Hurier et al. (2014) Hurier, G., Aghanim, N., & Douspis, M. 2014, A&A, 568, A57
- Hurier et al. (2013) Hurier, G., Macías-Pérez, J. F., & Hildebrandt, S. 2013, A&A, 558, A118
- Kalberla et al. (2005) Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775
- Koester et al. (2007) Koester, B. P., McKay, T. A., Annis, J., et al. 2007, ApJ, 660, 239
- Komatsu & Kitayama (1999) Komatsu, E. & Kitayama, T. 1999, ApJ, 526, L1
- Komatsu & Seljak (2002) Komatsu, E. & Seljak, U. 2002, MNRAS, 336, 1256
- Komatsu et al. (2011) Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18
- Lieu et al. (2006) Lieu, R., Mittaz, J. P. D., & Zhang, S.-N. 2006, ApJ, 648, 176
- Marriage et al. (2011) Marriage, T. A., Acquaviva, V., Ade, P. A. R., et al. 2011, ApJ, 737, 61
- Mauch et al. (2003) Mauch, T., Murphy, T., Buttery, H. J., et al. 2003, MNRAS, 342, 1117
- Mauch et al. (2008) Mauch, T., Murphy, T., Buttery, H. J., et al. 2008, VizieR Online Data Catalog, 8081, 0
- Melin et al. (2011) Melin, J.-B., Bartlett, J. G., Delabrouille, J., et al. 2011, A&A, 525, A139
- Mewe et al. (1985) Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 1985, A&AS, 62, 197
- Mo & White (1996) Mo, H. J. & White, S. D. M. 1996, MNRAS, 282, 347
- Morrison & McCammon (1983) Morrison, R. & McCammon, D. 1983, ApJ, 270, 119
- Myers et al. (2004) Myers, A. D., Shanks, T., Outram, P. J., Frith, W. J., & Wolfendale, A. W. 2004, MNRAS, 347, L67
- Piffaretti et al. (2011) Piffaretti, R., Arnaud, M., Pratt, G. W., Pointecouteau, E., & Melin, J.-B. 2011, A&A, 534, A109
- Planck Collaboration early results. VIII (2011) Planck Collaboration early results. VIII. 2011, A&A, 536, A8
- Planck Collaboration early results. X (2011) Planck Collaboration early results. X. 2011, A&A, 536, A10
- Planck Collaboration early results. XI (2011) Planck Collaboration early results. XI. 2011, A&A, 536, A11
- Planck Collaboration Int. results. VII (2013) Planck Collaboration Int. results. VII. 2013, A&A, 550, A133
- Planck Collaboration results. I (2014) Planck Collaboration results. I. 2014, A&A, 571, A1
- Planck Collaboration results. IX (2014) Planck Collaboration results. IX. 2014, A&A, 571, A9
- Planck Collaboration results. VII (2014) Planck Collaboration results. VII. 2014, A&A, 571, A7
- Planck Collaboration results. VIII (2014) Planck Collaboration results. VIII. 2014, A&A, 571, A8
- Planck Collaboration results. XII (2014) Planck Collaboration results. XII. 2014, A&A, 571, A12
- Planck Collaboration results. XIII (2013) Planck Collaboration results. XIII. 2013, e-prints ArXiv: 1303.5073
- Planck Collaboration results. XVI (2014) Planck Collaboration results. XVI. 2014, A&A, 571, A16
- Planck Collaboration results. XX (2014) Planck Collaboration results. XX. 2014, A&A, 571, A20
- Planck Collaboration results. XXI (2014) Planck Collaboration results. XXI. 2014, A&A, 571, A21
- Planck Collaboration results. XXIX (2014) Planck Collaboration results. XXIX. 2014, A&A, 571, A29
- Planck Collaboration results. XXVIII (2014) Planck Collaboration results. XXVIII. 2014, A&A, 571, A28
- Pratt et al. (2009) Pratt, G. W., Croston, J. H., Arnaud, M., & Böhringer, H. 2009, A&A, 498, 361
- Reichardt et al. (2012) Reichardt, C. L., Shaw, L., Zahn, O., et al. 2012, ApJ, 755, 70
- Reichert et al. (2011) Reichert, A., Böhringer, H., Fassbender, R., & Mühlegger, M. 2011, A&A, 535, A4
- Rykoff et al. (2013) Rykoff, E. S., Rozo, E., Busha, M. T., et al. 2013, ArXiv e-prints
- Sehgal et al. (2011) Sehgal, N., Trac, H., Acquaviva, V., et al. 2011, ApJ, 732, 44
- Shirokoff et al. (2011) Shirokoff, E., Reichardt, C. L., Shaw, L., et al. 2011, ApJ, 736, 61
- Sievers et al. (2013) Sievers, J. L., Hlozek, R. A., Nolta, M. R., et al. 2013, J. Cosmology Astropart. Phys., 10, 60
- Sunyaev & Zeldovich (1969) Sunyaev, R. A. & Zeldovich, Y. B. 1969, Nature, 223, 721
- Sunyaev & Zeldovich (1972) Sunyaev, R. A. & Zeldovich, Y. B. 1972, Comments on Astrophysics and Space Physics, 4, 173
- Taburet et al. (2011) Taburet, N., Hernández-Monteagudo, C., Aghanim, N., Douspis, M., & Sunyaev, R. A. 2011, MNRAS, 418, 2207
- The Dark Energy Survey Collaboration (2005) The Dark Energy Survey Collaboration. 2005, ArXiv Astrophysics e-prints
- Tinker et al. (2008) Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 688, 709
- Tristram et al. (2005) Tristram, M., Macías-Pérez, J. F., Renault, C., & Santos, D. 2005, MNRAS, 358, 833
- Vanderlinde et al. (2010) Vanderlinde, K., Crawford, T. M., de Haan, T., et al. 2010, ApJ, 722, 1180
- Vikhlinin et al. (2009) Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009, ApJ, 692, 1060
- Voges et al. (1999) Voges, W., Aschenbach, B., Boller, T., et al. 1999, A&A, 349, 389