Complex astrophysical systems often exhibit low-scatter relations between observable properties (e.g., luminosity, velocity dispersion, oscillation period). These scaling relations illuminate the underlying physics, and can provide observational tools for estimating masses and distances. Machine learning can provide a fast and systematic way to search for new scaling relations (or for simple extensions to existing relations) in abstract high-dimensional parameter spaces. We use a machine learning tool called symbolic regression (SR), which models patterns in a dataset in the form of analytic equations. We focus on the Sunyaev-Zeldovich flux−cluster mass relation ( Y SZ − M ), the scatter in which affects inference of cosmological parameters from cluster abundance data. Using SR on the data from the IllustrisTNG hydrodynamical simulation, we find a new proxy for cluster mass which combines Y SZ and concentration of ionized gas ( c gas ): M ∝ Y conc 3/5 ≡ Y SZ 3/5 (1 − A c gas ). Y conc reduces the scatter in the predicted M by ∼20 − 30% for large clusters ( M ≳ 10 14 h −1 M ⊙ ), as compared to using just Y SZ . We show that the dependence on c gas is linked to cores of clusters exhibiting larger scatter than their outskirts. Finally, we test Y conc on clusters from CAMELS simulations and show that Y conc is robust against variations in cosmology, subgrid physics, and cosmic variance. Our results and methodology can be useful for accurate multiwavelength cluster mass estimation from upcoming CMB and X-ray surveys like ACT, SO, eROSITA and CMB-S4.
more »
« less
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.
more »
« less
- PAR ID:
- 10409892
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 522
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2628-2643
- Size(s):
- p. 2628-2643
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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$$.more » « less
-
Abstract Upcoming cosmic shear analyses will precisely measure the cosmic matter distribution at low redshifts. At these redshifts, the matter distribution is affected by galaxy formation physics, primarily baryonic feedback from star formation and active galactic nuclei. Employing measurements from theMagneticumandIllustrisTNGsimulations and a dark matter + baryon (DMB) halo model, this paper demonstrates that Sunyaev-Zel'dovich (SZ) effect observations of galaxy clusters, whose masses have been calibrated using weak gravitational lensing, can constrain the baryonic impact on cosmic shear with statistical and systematic errors subdominant to the measurement errors of DES-Y3 and LSST-Y1, with systematic errors on S8and Ωmreaching 10% and 50% of the statistical errors, respectively. For LSST-Y6 and Roman surveys, these systematic errors increase to 150% and 100% of the statistical errors, indicating the necessity for further model developments for future surveys. We further dissect the contributions from different scales and halos with different masses to cosmic shear, highlighting the dominant role of SZ clusters at scales critical for cosmic shear analyses. These findings suggest a promising avenue for future joint analyses of Cosmic Microwave Background (CMB) and lensing surveys.more » « less
-
Abstract We show the improvement to cosmological constraints from galaxy cluster surveys with the addition of cosmic microwave background (CMB)-cluster lensing data. We explore the cosmological implications of adding mass information from the 3.1σdetection of gravitational lensing of the CMB by galaxy clusters to the Sunyaev–Zel’dovich (SZ) selected galaxy cluster sample from the 2500 deg2SPT-SZ survey and targeted optical and X-ray follow-up data. In the ΛCDM model, the combination of the cluster sample with the Planck power spectrum measurements prefers . Adding the cluster data reduces the uncertainty on this quantity by a factor of 1.4, which is unchanged whether the 3.1σCMB-cluster lensing measurement is included or not. We then forecast the impact of CMB-cluster lensing measurements with future cluster catalogs. Adding CMB-cluster lensing measurements to the SZ cluster catalog of the ongoing SPT-3G survey is expected to improve the expected constraint on the dark energy equation of statewby a factor of 1.3 toσ(w) = 0.19. We find the largest improvements from CMB-cluster lensing measurements to be forσ8, where adding CMB-cluster lensing data to the cluster number counts reduces the expected uncertainty onσ8by respective factors of 2.4 and 3.6 for SPT-3G and CMB-S4.more » « less
-
ABSTRACT JKCS041 ($z=1.8$) is one of the most distant galaxy cluster systems known, seen when the Universe was less than 4 billion years old. Recent Sunyaev–Zeldovich (SZ) observations show a temperature decrement that is less than expected based on mass estimates of the system from X-ray, weak gravitational lensing, and galaxy richness measurements. In this paper, we seek to explain the observables – in particular the low SZ decrement and single SZ peak, the projected offset between the X-ray and SZ peaks of $$\approx$$220 kpc, the gas mass measurements and the lensing mass estimate. We use the gamer-2 hydrodynamic code to carry out idealized numerical simulations of cluster mergers and compare resulting synthetic maps with the observational data. Generically, a merger process is necessary to reproduce the observed offset between the SZ and X-ray peaks. From our exploration of parameter space, seen a few tenths of a Gyr after first core passage, two components with total mass of $$\approx 2\times 10^{14} \,\text{M}_\odot$$, mass ratio of $$\approx$$2:3, gas fraction of $0.05-0.1$, and Navarro, Frenk and White mass density profile concentrations c$$\approx$$ 5 are scenarios that are consistent with the observational data. For consistency with the SZ and X-ray measurements, our simulations exclude total mass in excess of $$\approx 3\times 10^{14} {\rm M}_{\odot }$$, primarily based on the SZ signal. The mass ratio is constrained by the SZ–X-ray offset and magnitude of the SZ signal, ruling out systems with equal and vastly different masses.more » « less
An official website of the United States government
