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
Characterizing hydrostatic mass bias with mock-X
ABSTRACT Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we compare simulations to observations in the correct way. Here we introduce the mock-X analysis framework, a multiwavelength tool that generates synthetic images from cosmological simulations and derives halo properties via observational methods. We detail our methods for generating optical, Compton-y and X-ray images. Outlining our synthetic X-ray image analysis method, we demonstrate the capabilities of the framework by exploring hydrostatic mass bias for the IllustrisTNG, BAHAMAS, and MACSIS simulations. Using simulation derived profiles we find an approximately constant bias b ≈ 0.13 with cluster mass, independent of hydrodynamical method, or subgrid physics. However, the hydrostatic bias derived from synthetic observations is mass-dependent, increasing to b = 0.3 for the most massive clusters. This result is driven by a single temperature fit to a spectrum produced by gas with a wide temperature distribution in quasi-pressure equilibrium. The spectroscopic temperature and mass estimate are biased low by cooler gas dominating the emission, due to its quadratic density dependence. The bias and the scatter in estimated mass remain independent of the numerical method and subgrid physics. Our results are consistent with current observations and future surveys will contain sufficient samples of massive clusters to confirm the mass dependence of the hydrostatic bias.
more »
« less
- NSF-PAR ID:
- 10283622
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 506
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 2533 to 2550
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT Dynamically relaxed galaxy clusters have long played an important role in galaxy cluster studies because it is thought their properties can be reconstructed more precisely and with less systematics. As relaxed clusters are desirable, there exist a plethora of criteria for classifying a galaxy cluster as relaxed. In this work, we examine 9 commonly used observational and theoretical morphological metrics extracted from $54\, 000$mock-X synthetic X-ray images of galaxy clusters taken from the IllustrisTNG, BAHAMAS, and MACSIS simulation suites. We find that the simulated criteria distributions are in reasonable agreement with the observed distributions. Many criteria distributions evolve as a function of redshift, cluster mass, numerical resolution, and subgrid physics, limiting the effectiveness of a single relaxation threshold value. All criteria are positively correlated with each other, however, the strength of the correlation is sensitive to redshift, mass, and numerical choices. Driven by the intrinsic scatter inherent to all morphological metrics and the arbitrary nature of relaxation threshold values, we find the consistency of relaxed subsets defined by the different metrics to be relatively poor. Therefore, the use of relaxed cluster subsets introduces significant selection effects that are non-trivial to resolve.more » « less
-
Abstract The hot gas that constitutes the intracluster medium (ICM) has been studied at X-ray and millimeter/submillimeter wavelengths (Sunyaev–Zel’dovich effect) for decades. Fast radio bursts (FRBs) offer an additional method of directly measuring the ICM and gas surrounding clusters via observables such as dispersion measure (DM) and Faraday rotation measure. We report the discovery of two FRB sources detected with the Deep Synoptic Array whose host galaxies belong to massive galaxy clusters. In both cases, the FRBs exhibit excess extragalactic DM, some of which likely originate in the ICM of their respective clusters. FRB 20220914A resides in the galaxy cluster A2310 at z = 0.1125 with a projected offset from the cluster center of 520 ± 50 kpc. The host of a second source, FRB 20220509G, is an elliptical galaxy at z = 0.0894 that belongs to the galaxy cluster A2311 at the projected offset of 870 ± 50 kpc. These sources represent the first time an FRB has been localized to a galaxy cluster. We combine our FRB data with archival X-ray, Sunyaev–Zel'dovich (SZ), and optical observations of these clusters in order to infer properties of the ICM, including a measurement of gas temperature from DM and y SZ of 0.8–3.9 keV. We then compare our results to massive cluster halos from the IllustrisTNG simulation. Finally, we describe how large samples of localized FRBs from future surveys will constrain the ICM, particularly beyond the virial radius of clusters.more » « less
-
more » « less
ABSTRACT We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq (1\!-\!3)\times 10^{15}{\, \rm M_\odot }$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\!-\!0.2R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(\sigma /m)= 0.5\!-\!1\, {\rm cm}^2\, {\rm g}^{-1}$ than CDM and SIDM with $(\sigma /m)=0.1\, {\rm cm}^2\, {\rm g}^{-1}$. Our statistical analysis indicates that the latter two models are disfavoured at the $68{{\ \rm per\ cent}}$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.
-
Abstract We provide the first combined cosmological analysis of the South Pole Telescope (SPT) and Planck cluster catalogs. The aim is to provide an independent calibration for Planck scaling relations, exploiting the cosmological constraining power of the SPT-SZ cluster catalog and its dedicated weak lensing (WL) and X-ray follow-up observations. We build a new version of the Planck cluster likelihood. In the ν Λ CDM scenario, focusing on the mass slope and mass bias of Planck scaling relations, we find α SZ = 1.49 − 0.10 + 0.07 and 1 − b SZ = 0.69 − 0.14 + 0.07 , respectively. The results for the mass slope show a ∼4 σ departure from the self-similar evolution, α SZ ∼ 1.8. This shift is mainly driven by the matter density value preferred by SPT data, Ω m = 0.30 ± 0.03, lower than the one obtained by Planck data alone, Ω m = 0.37 − 0.06 + 0.02 . The mass bias constraints are consistent both with outcomes of hydrodynamical simulations and external WL calibrations, (1 − b ) ∼ 0.8, and with results required by the Planck cosmic microwave background cosmology, (1 − b ) ∼ 0.6. From this analysis, we obtain a new catalog of Planck cluster masses M 500 . We estimate the ratio between the published Planck M SZ masses and our derived masses M 500 , as a “measured mass bias,” 1 − b M . We analyze the mass, redshift, and detection noise dependence of 1 − b M , finding an increasing trend toward high redshift and low mass. These results mimic the effect of departure from self-similarity in cluster evolution, showing different dependencies for the low-mass, high-mass, low- z , and high- z regimes.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/stab1276
