Augmenting astrophysical scaling relations with machine learning : application to reducing the SZ flux-mass scatter
share
Complex systems (stars, supernovae, galaxies, and clusters) often exhibit low scatter relations between observable properties (e.g., luminosity, velocity dispersion, oscillation period, temperature). These scaling relations can illuminate the underlying physics and can provide observational tools for estimating masses and distances. Machine learning can provide a 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 the patterns in a given 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 large clusters (M≳ 10^14 h^-1 M_⊙) at both high and low redshifts, 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 simulations of the CAMELS project and show that Y_conc is robust against variations in cosmology, astrophysics, subgrid physics, and cosmic variance. Our results and methodology can be useful for accurate multiwavelength cluster mass estimation from current and upcoming CMB and X-ray surveys like ACT, SO, SPT, eROSITA and CMB-S4.
READ FULL TEXT