Sunyaev-Zeldovich Surveys Research Articles

Cosmological shocks around galaxy clusters: a coherent investigation with DES, SPT, and ACT

ABSTRACT We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogues from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev–Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around 105 clusters with mass and redshift ranges $10^{13.7} \lt M_{\rm 200m}/\, {\rm M}_\odot \lt 10^{15.5}$ and 0.1 < z < 2, and the total sky coverage of the maps is $\approx 15\, 000 \deg ^2$. We find a clear pressure deficit at R/R200m ≈ 1.1 in SZ profiles around both ACT and SPT clusters, estimated at 6σ significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2–5 compared to predictions, and we discuss possible interpretations of this behaviour. An oriented stacking of clusters – where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogues – shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius.

The SZ flux-mass (Y–M) relation at low-halo masses: improvements with symbolic regression and strong constraints on baryonic feedback

ABSTRACT Feedback from active galactic nuclei (AGNs) and supernovae can affect measurements of integrated Sunyaev–Zeldovich (SZ) flux of haloes (YSZ) from cosmic microwave background (CMB) surveys, and cause its relation with the halo mass (YSZ–M) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with extensive variations in feedback prescriptions. We use a combination of two machine learning tools (random forest and symbolic regression) to search for analogues of the Y–M relation which are more robust to feedback processes for low masses ($M\lesssim 10^{14}\, \mathrm{ h}^{-1} \, \mathrm{ M}_\odot$); we find that simply replacing Y → Y(1 + M*/Mgas) in the relation makes it remarkably self-similar. This could serve as a robust multiwavelength mass proxy for low-mass clusters and galaxy groups. Our methodology can also be generally useful to improve the domain of validity of other astrophysical scaling relations. We also forecast that measurements of the Y–M relation could provide per cent level constraints on certain combinations of feedback parameters and/or rule out a major part of the parameter space of supernova and AGN feedback models used in current state-of-the-art hydrodynamic simulations. Our results can be useful for using upcoming SZ surveys (e.g. SO, CMB-S4) and galaxy surveys (e.g. DESI and Rubin) to constrain the nature of baryonic feedback. Finally, we find that the alternative relation, Y–M*, provides complementary information on feedback than Y–M.

Emulating Sunyaev–Zeldovich images of galaxy clusters using autoencoders

ABSTRACT We develop a machine-learning (ML) algorithm that generates high-resolution thermal Sunyaev–Zeldovich (SZ) maps of novel galaxy clusters given only halo mass and mass accretion rate (MAR). The algorithm uses a conditional variational autoencoder (CVAE) in the form of a convolutional neural network and is trained with SZ maps generated from the IllustrisTNG simulation. Our method can reproduce many of the details of galaxy clusters that analytical models usually lack, such as internal structure and aspherical distribution of gas created by mergers, while achieving the same computational feasibility, allowing us to generate mock SZ maps for over 105 clusters in 30 s on a laptop. We show that the model is capable of generating novel clusters (i.e. not found in the training set) and that the model accurately reproduces the effects of mass and MAR on the SZ images, such as scatter, asymmetry, and concentration, in addition to modelling merging sub-clusters. This work demonstrates the viability of ML-based methods for producing the number of realistic, high-resolution maps of galaxy clusters necessary to achieve statistical constraints from future SZ surveys.

Removing the giants and learning from the crowd: A new SZ power spectrum method and revised Compton y-map analysis

ABSTRACT The Sunyaev–Zeldovich (SZ) effect provides a powerful cosmological probe, which traditionally is approached independently as cluster number count (CNC) or power spectrum (PS) analysis. Here, we devise a new method for analysing the y-map by introducing the survey completeness function, conventionally only used in the CNC analysis, in the yy-PS modelling. This provides a systematic method, based mainly on SZ observables, for obtaining two complementary y-maps, one incorporating detected/resolved clusters and the other relying only on diffuse/unresolved SZ contributions. We use the catalogue of clusters obtained in the Planck CNC analysis to define the completeness function linking these two y-maps. The split depends on the chosen signal-to-noise detection threshold, which we vary in our discussion. We carefully propagate the effect of completeness cuts on the non-Gaussian error contributions in the yy-PS analysis, highlighting the benefits of masking massive clusters. Our analysis of the Planck yy-PS for the unresolved component yields a mass bias of b = 0.15 ± 0.04, consistent with the standard value (b ≈ 0.2), in comparison to b = 0.4 ± 0.05 for the total yy-PS. We find indications for this drift being driven by the CIB-tSZ cross-correlation, which dominantly originates from clusters in the resolved component of the y-map. Another possible explanation is the presence of a mass-dependent bias, which has been theoretically motivated and can be quantified with our novel method. We furthermore find first hints for the presence of the 2-halo terms in the yy-PS. Finally, the proposed method provides a new framework for combining the complementary information of the CNC and PS analyses in upcoming SZ surveys.

Constraining the masses of high-redshift clusters with weak lensing: Revised shape calibration testing for the impact of stronger shears and increased blending

Weak lensing measurements suffer from well-known shear estimation biases, which can be partially corrected for with the use of image simulations. In this work we present an analysis of simulated images that mimic Hubble Space Telescope/Advance Camera for Surveys observations of high-redshift galaxy clusters, including cluster specific issues such as non-weak shear and increased blending. Our synthetic galaxies have been generated to have similar observed properties as the background-selected source samples studied in the real images. First, we used simulations with galaxies placed on a grid to determine a revised signal-to-noise-dependent (S/NKSB) correction for multiplicative shear measurement bias, and to quantify the sensitivity of our KSB+ bias calibration to mismatches of galaxy or PSF properties between the real data and the simulations. Next, we studied the impact of increased blending and light contamination from cluster and foreground galaxies, finding it to be negligible for high-redshift (z > 0.7) clusters, whereas shear measurements can be affected at the ∼1% level for lower redshift clusters given their brighter member galaxies. Finally, we studied the impact of fainter neighbours and selection bias using a set of simulated images that mimic the positions and magnitudes of galaxies in Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) data, thereby including realistic clustering. While the initial SExtractor object detection causes a multiplicative shear selection bias of −0.028 ± 0.002, this is reduced to −0.016 ± 0.002 by further cuts applied in our pipeline. Given the limited depth of the CANDELS data, we compared our CANDELS-based estimate for the impact of faint neighbours on the multiplicative shear measurement bias to a grid-based analysis, to which we added clustered galaxies to even fainter magnitudes based on Hubble Ultra Deep Field data, yielding a refined estimate of ∼ − 0.013. Our sensitivity analysis suggests that our pipeline is calibrated to an accuracy of ∼0.015 once all corrections are applied, which is fully sufficient for current and near-future weak lensing studies of high-redshift clusters. As an application, we used it for a refined analysis of three highly relaxed clusters from the South Pole Telescope Sunyaev-Zeldovich survey, where we now included measurements down to the cluster core (r > 200 kpc) as enabled by our work. Compared to previously employed scales (r > 500 kpc), this tightens the cluster mass constraints by a factor 1.38 on average.

Discovery of a Powerful >1061 erg AGN Outburst in the Distant Galaxy Cluster SPT-CLJ0528-5300

Abstract We present ∼103 ks of Chandra observations of the galaxy cluster SPT-CLJ0528-5300 (SPT0528, z = 0.768). This cluster harbors the most radio-loud (L 1.4GHz = 1.01 × 1033 erg s−1 Hz−1) central active galactic nucleus (AGN) of any cluster in the South Pole Telescope (SPT) Sunyaev–Zeldovich survey with available X-ray data. We find evidence of AGN-inflated cavities in the X-ray emission, which are consistent with the orientation of the jet direction revealed by Australia Telescope Compact Array radio data. The combined probability that two such depressions—each at ∼1.4–1.8σ significance, oriented ∼180° apart and aligned with the jet axis—would occur by chance is 0.1%. At ≳1061 erg, the outburst in SPT0528 is among the most energetic known in the universe, and certainly the most powerful known at z > 0.25. This work demonstrates that such powerful outbursts can be detected even in shallow X-ray exposures out to relatively high redshifts (z ∼ 0.8), providing an avenue for studying the evolution of extreme AGN feedback. The ratio of the cavity power ( P cav = ( 9.4 ± 5.8 ) × 10 45 erg s−1) to the cooling luminosity (L cool = (1.5 ± 0.5) × 1044 erg s−1) for SPT0528 is among the highest measured to date. If, in the future, additional systems are discovered at similar redshifts with equally high P cav/L cool ratios, it would imply that the feedback/cooling cycle was not as gentle at high redshifts as in the low-redshift universe.

A matched filter approach for blind joint detection of galaxy clusters in X-ray and SZ surveys

The combination of X-ray and Sunyaev–Zeldovich (SZ) observations can potentially improve the cluster detection efficiency, when compared to using only one of these probes, since both probe the same medium, the hot ionized gas of the intra-cluster medium. We present a method based on matched multifrequency filters (MMF) for detecting galaxy clusters from SZ and X-ray surveys. This method builds on a previously proposed joint X-ray–SZ extraction method and allows the blind detection of clusters, that is finding new clusters without knowing their position, size, or redshift, by searching on SZ and X-ray maps simultaneously. The proposed method is tested using data from the ROSAT all-sky survey and from the Planck survey. The evaluation is done by comparison with existing cluster catalogues in the area of the sky covered by the deep SPT survey. Thanks to the addition of the X-ray information, the joint detection method is able to achieve simultaneously better purity, better detection efficiency, and better position accuracy than its predecessor Planck MMF, which is based on SZ maps alone. For a purity of 85%, the X-ray–SZ method detects 141 confirmed clusters in the SPT region; to detect the same number of confirmed clusters with Planck MMF, we would need to decrease its purity to 70%. We provide a catalogue of 225 sources selected by the proposed method in the SPT footprint, with masses ranging between 0.7 and 14.5 ×1014 M⊙ and redshifts between 0.01 and 1.2.

Optical Identifications of High-Redshift Galaxy Clusters from the Planck Sunyaev–Zeldovich Survey

We present the results of optical identifications and spectroscopic redshifts measurements for galaxy clusters from 2-nd Planck catalogue of Sunyaev-Zeldovich sources (PSZ2), located at high redshifts, $z\approx0.7-0.9$. We used the data of optical observations obtained with Russian-Turkish 1.5-m telescope (RTT150), Sayan observatory 1.6-m telescope, Calar Alto 3.5-m telescope and 6-m SAO RAS telescope (Bolshoi Teleskop Alt-azimutalnyi, BTA). Spectroscopic redshift measurements were obtained for seven galaxy clusters, including one cluster, PSZ2 G126.57+51.61, from the cosmological sample of PSZ2 catalogue. In central regions of two clusters, PSZ2 G069.39+68.05 and PSZ2 G087.39-34.58, the strong gravitationally lensed background galaxies are found, one of them at redshift $z=4.262$. The data presented below roughly double the number of known galaxy clusters in the second Planck catalogue of Sunyaev-Zeldovich sources at high redshifts, $z\approx0.8$.

Planck2015 results

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

Measuring the dynamical state ofPlanckSZ-selected clusters: X-ray peak – BCG offset

We want to characterize the dynamical state of galaxy clusters detected with the Sunyaev-Zeldovich (SZ) effect by Planck and compare them with the dynamical state of clusters selected in X-rays survey. We analyzed a representative subsample of the Planck SZ catalogue, containing the 132 clusters with the highest signal to noise ratio and characterize their dynamical state using as indicator the projected offset between the peak of the X-ray emission and the position of the Brightest cluster galaxy. We study the distribution of our indicator in our sample and compare it to its distribution in X-ray selected samples (HIFLUGCS, MACS and REXCESS). The distributions are significantly different and the fraction of relaxed objects is smaller in the Planck sample ($52 \pm 4 \%$) than in X-ray samples ($\simeq 74\%$) We interpret this result as an indication of different selection effects affecting X-rays (e.g. "cool core bias") and SZ surveys of galaxy clusters.

COOL CORE BIAS IN SUNYAEV–ZEL’DOVICH GALAXY CLUSTER SURVEYS

Sunyaev-Zeldovich (SZ) surveys find massive clusters of galaxies by measuring the inverse Compton scattering of cosmic microwave background off of intra-cluster gas. The cluster selection function from such surveys is expected to be nearly independent of redshift and cluster astrophysics. In this work, we estimate the effect on the observed SZ signal of centrally-peaked gas density profiles (cool cores) and radio emission from the brightest cluster galaxy (BCG) by creating mock observations of a sample of clusters that span the observed range of classical cooling rates and radio luminosities. For each cluster, we make simulated SZ observations by the South Pole Telescope and characterize the cluster selection function, but note that our results are broadly applicable to other SZ surveys. We find that the inclusion of a cool core can cause a change in the measured SPT significance of a cluster between 0.01% - 10% at z > 0.3, increasing with cuspiness of the cool core and angular size on the sky of the cluster (i.e., decreasing redshift, increasing mass). We provide quantitative estimates of the bias in the SZ signal as a function of a gas density cuspiness parameter, redshift, mass, and the 1.4 GHz radio luminosity of the central AGN. Based on this work, we estimate that, for the Phoenix cluster (one of the strongest cool cores known), the presence of a cool core is biasing the SZ significance high by ~ 6%. The ubiquity of radio galaxies at the centers of cool core clusters will offset the cool core bias to varying degrees.

Probing the primordial power spectrum with cluster number counts

We investigate how well galaxy cluster number counts can constrain the primordial power spectrum. Measurements of the primary anisotropies in the cosmic microwave background may be limited, by the presence of foregrounds from secondary sources, to probing the primordial power spectrum at wave numbers less than about $0.30h\text{ }\text{ }{\mathrm{Mpc}}^{\ensuremath{-}1}$. We break up the primordial power spectrum into a number of nodes and interpolate linearly between each node. This allows us to show that cluster number counts could then extend the constraints on the form of the primordial power spectrum up to wave numbers of about $0.45h\text{ }\text{ }{\mathrm{Mpc}}^{\ensuremath{-}1}$. We estimate combinations of constraints from PLANCK and SPT primary cosmic microwave background and their respective Sunyaev-Zeldovich surveys. We find that their constraining ability is limited by uncertainties in the mass-scaling relations. We also estimate the constraint from clusters detected from a SNAP-like gravitational lensing survey. As there is an unambiguous and simple relationship between the filtered shear of the lensing survey and the cluster mass, it may be possible to obtain much tighter constraints on the primordial power spectrum in this case.

SOUTHERN COSMOLOGY SURVEY. III. QSOs FROM COMBINED GALEX AND OPTICAL PHOTOMETRY

We present catalogs of QSO candidates selected using photometry from GALEX combined with SDSS in the Stripe 82 region and Blanco Cosmology Survey (BCS) near declination -55 degrees. The SDSS region contains ~700 objects with magnitude i < 20 and ~3600 objects with i < 21.5 in a ~60 square degree sky region, while the BCS region contains ~280 objects with magnitude i < 20 and ~2000 objects with i < 21.5 for a 11 square degree sky region that is being observed by three current microwave Sunyaev-Zeldovich surveys. Our QSO catalog is the first one in the BCS region. Deep GALEX exposures (~2000 seconds in FUV and NUV, except in three fields) provide high signal-to-noise photometry in the GALEX bands (FUV, NUV < 24.5 mag). From this data, we select QSO candidates using only GALEX and optical r-band photometry, using the method given by Atlee and Gould (2008). In the Stripe 82 field, 60% (30%) of the GALEX selected QSO's with optical magnitude i<20 (i<21.5) also appear in the Richards et al. (2008) QSO catalog constructed using 5-band optical SDSS photometry. Comparison with the same catalog by Richards et al. shows that the completeness of the sample is approximately 40%(25%). However, for regions of the sky with very low dust extinction, like the BCS 23hr field and the Stripe 82 between 0 and 10 degrees in RA, our completeness is close to 95%, demonstrating that deep GALEX observations are almost as efficient as multi-wavelength observations at finding QSO's. GALEX observations thus provide a viable alternate route to QSO catalogs in sky regions where u-band optical photometry is not available. The full catalog is available at http://www.ice.csic.es/personal/jimenez/PHOTOZ

High-redshift Fanaroff-Riley type II radio sources: large-scale X-ray environment

We investigate the properties of the environment around 20 powerful radio galaxies and quasars at redshifts between 0.45 and 1. Using XMM-Newtonand Chandra observations we probe the spatial distribution and the temperature of the cl uster gas. We find that more than 60 per cent of powerful radio sources in the redshift range of our sample lie in a cluster of X-ray luminosity greater than 10 44 erg s 1 , and all but one of the narrow-line radio galaxies, for which the emission from the nucleus is obscured by a torus, lie in a cluster environment. For broad-line quasars the X-ray emission from the core dominates and it is more difficult to measure the cluster environment. However, within the statistical uncertainties we find no significant difference in the properties of the environment as a function of the orientation to the line of sight of the radio jet. This is in agreement with un ification schemes. Our results have important implications for cluster surveys, as cluste rs around powerful radio sources tend to be excluded from X-ray and Sunyaev-Zeldovich surveys of galaxy clusters, and thus can introduce an important bias in the cluster luminosity function. Most of the radio sources are found close to pressure balance with the environment in which they lie, but the two lowexcitation radio galaxies of the sample are observed to be under-pressured. This may be the first observational indication for the presence of non-radi ative particles in the lobes of some powerful radio galaxies. We find that the clusters around rad io sources in the redshift range of our sample have a steeper entropy-temperature relation than local clusters, and the slope is in agreement with the predictions of self-similar gravitat ional heating models for cluster gas infall. This suggests that selection by AGN finds systems les s affected by AGN feedback than the local average. We speculate that this is because the AGN in our sample are sufficiently luminous and rare that their AGN activity is too recent to have caused the onset of measurable feedback and increased entropy in the clusters, especially in the cooler ones where locally the effect of feedback are expected to be most evident. If this is confirmed by forthcoming X-ray missions it will improve our understanding of the heating and cooling processes in high-redshift galaxy clusters.

Redshift Accuracy Requirements for Future Supernova and Number Count Surveys

We investigate the required redshift accuracy of type Ia supernova and cluster number-count surveys in order for the redshift uncertainties not to contribute appreciably to the dark energy parameter error budget. For the SNAP supernova experiment, we find that, without the assistance of ground-based measurements, individual supernova redshifts would need to be determined to about 0.002 or better, which is a challenging but feasible requirement for a low-resolution spectrograph. However, we find that accurate redshifts for z < 0.1 supernovae, obtained with ground-based experiments, are sufficient to immunize the results against even relatively large redshift errors at high z. For the future cluster number-count surveys such as the South Pole Telescope, Planck or DUET, we find that the purely statistical error in photometric redshift is less important, and that the irreducible, systematic bias in redshift drives the requirements. The redshift bias will have to be kept below 0.001-0.005 per redshift bin (which is determined by the filter set), depending on the sky coverage and details of the definition of the minimal mass of the survey. Furthermore, we find that X-ray surveys have a more stringent required redshift accuracy than Sunyaev-Zeldovich (SZ) effect surveys since they use a shorter lever armmore » in redshift; conversely, SZ surveys benefit from their high redshift reach only so long as some redshift information is available for distant (zgtrsim1) clusters.« less

The hybrid SZ power spectrum: combining cluster counts and SZ fluctuations to probe gas physics

Sunyaev-Zeldovich (SZ) effect from a cosmological distribution of clusters carry information on the underlying cosmology as well as the cluster gas physics. In order to study either cosmology or clusters one needs to break the degeneracies between the two. We present a toy model showing how complementary informations from SZ power spectrum and the SZ flux counts, both obtained from upcoming SZ cluster surveys, can be used to mitigate the strong cosmological influence (especially that of sigma_8) on the SZ fluctuations. Once the strong dependence of the cluster SZ power spectrum on sigma_8 is diluted, the cluster power spectrum can be used as a tool in studying cluster gas structure and evolution. The method relies on the ability to write the Poisson contribution to the SZ power spectrum in terms the observed SZ flux counts. We test the toy model by applying the idea to simulations of SZ surveys.

Expectations for Sunyaev‐Zeldovich Cluster Counts: Mass Function versus X‐Ray Luminosity Function

We present a comparison of the Sunyaev-Zeldovich (SZ) cluster counts predicted by the Press-Schechter (PS) mass function (MF) and the X-ray luminosity function (XLF) of clusters. The employment of the cluster XLF, together with the observationally determined X-ray luminosity-temperature (LX-T) relation, may allow us to estimate the SZ cluster counts in a more realistic manner, although such an empirical approach depends sensitively on our current knowledge of the dynamical properties of intracluster gas and its cosmic evolution. Using both the nonevolving and evolving XLFs of clusters suggested by X-ray observations, we calculate the expectations for SZ surveys of clusters with X-ray luminosity LX ≥ 3 × 1044 ergs s-1 and LX ≥ 1 × 1043 ergs s-1 in the 0.5-2.0 band, respectively. The nonevolving XLF results in a significant excess of SZ cluster counts at high redshifts as compared with the evolving XLF, while a slightly steeper LX-T relation than the observed one is needed to reproduce the distributions of SZ clusters predicted by the standard PS formalism. It is pointed out that uncertainties in the cosmological application of future SZ cluster surveys via the standard PS formalism should be carefully studied.