Context. Galaxy clusters are an important tool for cosmology, and their detection and characterization are key goals for current and future surveys. Using data from the Wide-field Infrared Survey Explorer (WISE), the Massive and Distant Clusters of WISE Survey (MaDCoWS) located 2839 significant galaxy overdensities at redshifts 0.7 ≲ z ≲ 1.5, which included extensive follow-up imaging from the Spitzer Space Telescope to determine cluster richnesses. Concurrently, the Atacama Cosmology Telescope (ACT) has produced large area millimeter-wave maps in three frequency bands along with a large catalog of Sunyaev-Zeldovich (SZ)-selected clusters as part of its Data Release 5 (DR5). Aims. We aim to verify and characterize MaDCoWS clusters using measurements of, or limits on, their thermal SZ effect signatures. We also use these detections to establish the scaling relation between SZ mass and the MaDCoWS-defined richness. Methods. Using the maps and cluster catalog from DR5, we explore the scaling between SZ mass and cluster richness. We do this by comparing cataloged detections and extracting individual and stacked SZ signals from the MaDCoWS cluster locations. We use complementary radio survey data from the Very Large Array, submillimeter data from Herschel , and ACT 224 GHz data to assess the impact of contaminating sources on the SZ signals from both ACT and MaDCoWS clusters. We use a hierarchical Bayesian model to fit the mass-richness scaling relation, allowing for clusters to be drawn from two populations: one, a Gaussian centered on the mass-richness relation, and the other, a Gaussian centered on zero SZ signal. Results. We find that MaDCoWS clusters have submillimeter contamination that is consistent with a gray-body spectrum, while the ACT clusters are consistent with no submillimeter emission on average. Additionally, the intrinsic radio intensities of ACT clusters are lower than those of MaDCoWS clusters, even when the ACT clusters are restricted to the same redshift range as the MaDCoWS clusters. We find the best-fit ACT SZ mass versus MaDCoWS richness scaling relation has a slope of p 1 = 1.84 −0.14 +0.15 , where the slope is defined as M λ ∝ 15 p 1 and λ 15 is the richness. We also find that the ACT SZ signals for a significant fraction (∼57%) of the MaDCoWS sample can statistically be described as being drawn from a noise-like distribution, indicating that the candidates are possibly dominated by low-mass and unvirialized systems that are below the mass limit of the ACT sample. Further, we note that a large portion of the optically confirmed ACT clusters located in the same volume of the sky as MaDCoWS are not selected by MaDCoWS, indicating that the MaDCoWS sample is not complete with respect to SZ selection. Finally, we find that the radio loud fraction of MaDCoWS clusters increases with richness, while we find no evidence that the submillimeter emission of the MaDCoWS clusters evolves with richness. Conclusions. We conclude that the original MaDCoWS selection function is not well defined and, as such, reiterate the MaDCoWS collaboration’s recommendation that the sample is suited for probing cluster and galaxy evolution, but not cosmological analyses. We find a best-fit mass-richness relation slope that agrees with the published MaDCoWS preliminary results. Additionally, we find that while the approximate level of infill of the ACT and MaDCoWS cluster SZ signals (1–2%) is subdominant to other sources of uncertainty for current generation experiments, characterizing and removing this bias will be critical for next-generation experiments hoping to constrain cluster masses at the sub-percent level.
more »
« less
Emulating Sunyaev–Zeldovich images of galaxy clusters using autoencoders
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.
more » « less- NSF-PAR ID:
- 10397508
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 513
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 333-344
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
ABSTRACT Cosmological simulations are an important theoretical pillar for understanding non-linear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large-volume simulations of cosmic structure formation aiming to better understand physical processes on large scales and to help interpret upcoming large-scale galaxy surveys. We here focus on the full physics box MTNG740 that computes a volume of $740\, \mathrm{Mpc}^3$ with a baryonic mass resolution of $3.1\times ~10^7\, \mathrm{M_\odot }$ using arepo with 80.6 billion cells and the IllustrisTNG galaxy formation model. We verify that the galaxy properties produced by MTNG740 are consistent with the TNG simulations, including more recent observations. We focus on galaxy clusters and analyse cluster scaling relations and radial profiles. We show that both are broadly consistent with various observational constraints. We demonstrate that the SZ-signal on a deep light-cone is consistent with Planck limits. Finally, we compare MTNG740 clusters with galaxy clusters found in Planck and the SDSS-8 RedMaPPer richness catalogue in observational space, finding very good agreement as well. However, simultaneously matching cluster masses, richness, and Compton-y requires us to assume that the SZ mass estimates for Planck clusters are underestimated by 0.2 dex on average. Due to its unprecedented volume for a high-resolution hydrodynamical calculation, the MTNG740 simulation offers rich possibilities to study baryons in galaxies, galaxy clusters, and in large-scale structure, and in particular their impact on upcoming large cosmological surveys.
-
more » « less
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.
-
null (Ed.)Context. Inferences about dark matter, dark energy, and the missing baryons all depend on the accuracy of our model of large-scale structure evolution. In particular, with cosmological simulations in our model of the Universe, we trace the growth of structure, and visualize the build-up of bigger structures from smaller ones and of gaseous filaments connecting galaxy clusters. Aims. Here we aim to reveal the complexity of the large-scale structure assembly process in great detail and on scales from tens of kiloparsecs up to more than 10 Mpc with new sensitive large-scale observations from the latest generation of instruments. We also aim to compare our findings with expectations from our cosmological model. Methods. We used dedicated SRG/eROSITA performance verification (PV) X-ray, ASKAP/EMU Early Science radio, and DECam optical observations of a ~15 deg 2 region around the nearby interacting galaxy cluster system A3391/95 to study the warm-hot gas in cluster outskirts and filaments, the surrounding large-scale structure and its formation process, the morphological complexity in the inner parts of the clusters, and the (re-)acceleration of plasma. We also used complementary Sunyaev-Zeldovich (SZ) effect data from the Planck survey and custom-made Galactic total (neutral plus molecular) hydrogen column density maps based on the HI4PI and IRAS surveys. We relate the observations to expectations from cosmological hydrodynamic simulations from the Magneticum suite. Results. We trace the irregular morphology of warm and hot gas of the main clusters from their centers out to well beyond their characteristic radii, r 200 . Between the two main cluster systems, we observe an emission bridge on large scale and with good spatial resolution. This bridge includes a known galaxy group but this can only partially explain the emission. Most gas in the bridge appears hot, but thanks to eROSITA’s unique soft response and large field of view, we discover some tantalizing hints for warm, truly primordial filamentary gas connecting the clusters. Several matter clumps physically surrounding the system are detected. For the “Northern Clump,” we provide evidence that it is falling towards A3391 from the X-ray hot gas morphology and radio lobe structure of its central AGN. Moreover, the shapes of these X-ray and radio structures appear to be formed by gas well beyond the virial radius, r 100 , of A3391, thereby providing an indirect way of probing the gas in this elusive environment. Many of the extended sources in the field detected by eROSITA are also known clusters or new clusters in the background, including a known SZ cluster at redshift z = 1. We find roughly an order of magnitude more cluster candidates than the SPT and ACT surveys together in the same area. We discover an emission filament north of the virial radius of A3391 connecting to the Northern Clump. Furthermore, the absorption-corrected eROSITA surface brightness map shows that this emission filament extends south of A3395 and beyond an extended X-ray-emitting object (the “Little Southern Clump”) towards another galaxy cluster, all at the same redshift. The total projected length of this continuous warm-hot emission filament is 15 Mpc, running almost 4 degrees across the entire eROSITA PV observation field. The Northern and Southern Filament are each detected at >4 σ . The Planck SZ map additionally appears to support the presence of both new filaments. Furthermore, the DECam galaxy density map shows galaxy overdensities in the same regions. Overall, the new datasets provide impressive confirmation of the theoretically expected structure formation processes on the individual system level, including the surrounding warm-hot intergalactic medium distribution; the similarities of features found in a similar system in the Magneticum simulation are striking. Our spatially resolved findings show that baryons indeed reside in large-scale warm-hot gas filaments with a clumpy structure.more » « less
-
more » « less
ABSTRACT We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg2 SPT-SZ survey on the locations of clusters identified in that same data set. The sample contains 516 clusters with mean mass $\langle M_{\rm 200m}\rangle = 10^{14.9} \, {\rm M}_\odot$ and redshift 〈z〉 = 0.55. We analyse in parallel a set of zoom-in hydrodynamical simulations from the three hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = 1.08 ± 0.09, measured at 3.1σ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = 4.58 ± 1.24 at 2.0σ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyse clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius ($R_{\rm sh,\, acc}$) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, $R_{\rm sh,\, acc}/R_{\rm sp}\gt 2.16 \pm 0.59$.
