Search

We present the Baryon Pasted (BP) X-ray and thermal Sunyaev-Zel'dovich (tSZ) maps derived from the half-sky Uchuu Lightcone simulation. These BP-Uchuu maps are constructed using more than $75$ million dark matter halos with masses $M_{500c} \geq 10^{13} M_\odot$ within the redshift range $0 \leq z \leq 2$. A distinctive feature of our BP-Uchuu Lightcone maps is their capability to assess the influence of both extrinsic and intrinsic scatter caused by triaxial gaseous halos and internal gas characteristics, respectively, at the map level. We show that triaxial gas drives substantial scatter in X-ray luminosities of clusters and groups, accounting for nearly half of the total scatter in core-excised measurements. Additionally, scatter in the thermal pressure and gas density profiles of halos enhances the X-ray and SZ power spectra, leading to biases in cosmological parameter estimates. These findings are statistically robust due to the extensive sky coverage and large halo sample in the BP-Uchuu maps., Comment: 16 pages, 12 figures, submitted to ApJ. The BP-Uchuu maps are publicly available

The recent tension in the value of the cosmological parameter $S_8 \equiv \sigma_8(\Omega_M/0.3)^{1/2}$, which represents the amplitude of the matter density fluctuations of the universe, has not been resolved. In this work, we present constraints on $S_8$ with the X-ray angular power spectra of clusters and groups measured with the half-sky map from the eROSITA All Sky Survey data release 1 (eRASS1). Thanks to the extensive sky coverage of eRASS1, its power spectrum achieves unprecedented precision compared to previous measurements. Using a well-calibrated, physical halo gas model that includes astrophysics of feedback and non-thermal pressure support, we obtain $S_8 = 0.80^{+0.02}_{-0.01}$ with 1$\sigma$ uncertainty that is competitive against other cosmological probes. Our derived $S_8$ value is smaller than the primary CMB measurements from Planck, but still consistent to within $1\sigma$. We also obtain constraints on the astrophysics of feedback, non-thermal pressure, equation of state in cluster cores, and gas boundaries in clusters and groups. We discuss how additional X-ray observations, and cosmological surveys in microwave and optical, will further improve the cosmological constraints with the angular power spectrum. Our work demonstrates that the angular power spectrum of clusters and groups is a promising probe of both cosmology and astrophysics., Comment: 12 pages, 4 figures, submitted to ApJ

We present a study on the inference of cosmological and astrophysical parameters using stacked galaxy cluster profiles. Utilizing the CAMELS-zoomGZ simulations, we explore how various cluster properties--such as X-ray surface brightness, gas density, temperature, metallicity, and Compton-y profiles--can be used to predict parameters within the 28-dimensional parameter space of the IllustrisTNG model. Through neural networks, we achieve a high correlation coefficient of 0.97 or above for all cosmological parameters, including $\Omega_{\rm m}$, $H_0$, and $\sigma_8$, and over 0.90 for the remaining astrophysical parameters, showcasing the effectiveness of these profiles for parameter inference. We investigate the impact of different radial cuts, with bins ranging from $0.1R_{200c}$ to $0.7R_{200c}$, to simulate current observational constraints. Additionally, we perform a noise sensitivity analysis, adding up to 40\% Gaussian noise (corresponding to signal-to-noise ratios as low as 2.5), revealing that key parameters such as $\Omega_{\rm m}$, $H_0$, and the IMF slope remain robust even under extreme noise conditions. We also compare the performance of full radial profiles against integrated quantities, finding that profiles generally lead to more accurate parameter inferences. Our results demonstrate that stacked galaxy cluster profiles contain crucial information on both astrophysical processes within groups and clusters and the underlying cosmology of the universe. This underscores their significance for interpreting the complex data expected from next-generation surveys and reveals, for the first time, their potential as a powerful tool for parameter inference., Comment: Submitted to ApJ

A systematic comparison of the models of the circumgalactic medium (CGM) and their observables is crucial to understanding the predictive power of the models and constraining physical processes that affect the thermodynamics of CGM. This paper compares four analytic CGM models: precipitation, isentropic, cooling flow, and baryon pasting models for the hot, volume-filling CGM phase, all assuming hydrostatic or quasi-hydrostatic equilibrium. We show that for fiducial parameters of the CGM of a Milky-Way (MW) like galaxy ($\rm M_{vir} \sim 10^{12}~M_{\odot}$ at $z\sim 0$), the thermodynamic profiles -- entropy, density, temperature, and pressure -- show most significant differences between different models at small ($r\lesssim 30$ kpc) and large scales ($r\gtrsim 100$ kpc) while converging at intermediate scales. The slope of the entropy profile, which is one of the most important differentiators between models, is $\approx 0.8$ for the precipitation and cooling flow models, while it is $\approx0.6$ and 0 for the baryon pasting and isentropic models, respectively. We make predictions for various observational quantities for an MW mass halo for the different models, including the projected Sunyaev-Zeldovich (SZ) effect, soft X-ray emission (0.5--2 keV), dispersion measure, and column densities of oxygen ions (OVI, OVII, and OVIII) observable in absorption., Comment: Comments: 15 pages, 4 figures, 3 tables, accepted for publication in MNRAS

Galaxy formation models within cosmological hydrodynamical simulations contain numerous parameters with non-trivial influences over the resulting properties of simulated cosmic structures and galaxy populations. It is computationally challenging to sample these high dimensional parameter spaces with simulations, particularly for halos in the high-mass end of the mass function. In this work, we develop a novel sampling and reduced variance regression method, CARPoolGP, which leverages built-in correlations between samples in different locations of high dimensional parameter spaces to provide an efficient way to explore parameter space and generate low variance emulations of summary statistics. We use this method to extend the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) to include a set of 768 zoom-in simulations of halos in the mass range of $10^{13} - 10^{14.5} M_\odot\,h^{-1}$ that span a 28-dimensional parameter space in the IllustrisTNG model. With these simulations and the CARPoolGP emulation method, we explore parameter trends in the Compton $Y-M$, black hole mass-halo mass, and metallicity-mass relations, as well as thermodynamic profiles and quenched fractions of satellite galaxies. We use these emulations to provide a physical picture of the complex interplay between supernova and active galactic nuclei feedback. We then use emulations of the $Y-M$ relation of massive halos to perform Fisher forecasts on astrophysical parameters for future Sunyaev-Zeldovich observations and find a significant improvement in forecasted constraints. We publicly release both the simulation suite and CARPoolGP software package., Comment: The manuscript was submitted to arxiv after receiving and responding to comments from the first referee report

We present an investigation into a hitherto unexplored systematic that affects the accuracy of galaxy cluster mass estimates with weak gravitational lensing. Specifically, we study the covariance between the weak lensing signal, $\Delta\Sigma$, and the "true" cluster galaxy number count, $N_{\rm gal}$, as measured within a spherical volume that is void of projection effects. By quantifying the impact of this covariance on mass calibration, this work reveals a significant source of systematic uncertainty. Using the MDPL2 simulation with galaxies traced by the SAGE semi-analytic model, we measure the intrinsic property covariance between these observables within the 3D vicinity of the cluster, spanning a range of dynamical mass and redshift values relevant for optical cluster surveys. Our results reveal a negative covariance at small radial scales ($R \lesssim R_{\rm 200c}$) and a null covariance at large scales ($R \gtrsim R_{\rm 200c}$) across most mass and redshift bins. We also find that this covariance results in a $2-3\%$ bias in the halo mass estimates in most bins. Furthermore, by modeling $N_{\rm gal}$ and $\Delta\Sigma$ as multi-(log)-linear equations of secondary halo properties, we provide a quantitative explanation for the physical origin of the negative covariance at small scales. Specifically, we demonstrate that the $N_{\rm gal}$-$\Delta\Sigma$ covariance can be explained by the secondary properties of halos that probe their formation history. We attribute the difference between our results and the positive bias seen in other works with (mock)-cluster finders to projection effects. These findings highlight the importance of accounting for the covariance between observables in cluster mass estimation, which is crucial for obtaining accurate constraints on cosmological parameters.

High-resolution spectroscopy in soft X-rays will open a new window to map multiphase gas in galaxy clusters and probe physics of the intracluster medium (ICM), including chemical enrichment histories, circulation of matter and energy during large-scale structure evolution, stellar and black hole feedback, halo virialization, and gas mixing processes. An eV-level spectral resolution, large field-of-view, and effective area are essential to separate cluster emissions from the Galactic foreground and efficiently map the cluster outskirts. Several mission concepts that meet these criteria have been proposed recently, e.g., LEM, HUBS, and SuperDIOS. This theoretical study explores what information on ICM physics could be recovered with such missions and the associated challenges. We emphasize the need for a comprehensive comparison between simulations and observations to interpret the high-resolution spectroscopic observations correctly. Using Line Emission Mapper (LEM) characteristics as an example, we demonstrate that it enables the use of soft X-ray emission lines (e.g., O VII/VIII and Fe-L complex) from the cluster outskirts to measure the thermodynamic, chemical, and kinematic properties of the gas up to $r_{200}$ and beyond. By generating mock observations with full backgrounds, analysing their images/spectra with observational approaches, and comparing the recovered characteristics with true ones from simulations, we develop six key science drivers for future missions, including the exploration of multiphase gas in galaxy clusters (e.g., temperature fluctuations, phase-space distributions), metallicity, ICM gas bulk motions and turbulence power spectra, ICM-cosmic filament interactions, and advances for cluster cosmology., Comment: 24 pages, 26 figures, submitted to MNRAS. Comments are welcome

Known as the "Missing Baryon Problem", about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations used a novel stacking analysis and detected an OVII absorption line toward the sightline of a luminous quasar, hinting that the missing baryons may reside in the WHIM. To explore how the properties of the OVII absorption line depend on feedback physics, we compare the observational results with predictions obtained from the Cosmology and Astrophysics with MachinE Learning (CAMEL) Simulation suite. CAMELS consists of cosmological simulations with state-of-the-art supernova (SN) and active galactic nuclei (AGN) feedback models from the IllustrisTNG and SIMBA simulations, with varying strengths. We find that the simulated OVII column densities are higher in the outskirts of galaxies than in the large-scale WHIM, but they are consistently lower than those obtained in the Chandra observations, for all feedback runs. We establish that the OVII distribution is primarily sensitive to changes in the SN feedback prescription, whereas changes in the AGN feedback prescription have minimal impact. We also find significant differences in the OVII column densities between the IllustrisTNG and SIMBA runs. We conclude that the tension between the observed and simulated OVII column densities cannot be explained by the wide range of feedback models implemented in CAMELS., Comment: 11 pages, 9 figures. Matching to the published version on MNRAS

The eROSITA Final Equatorial Depth Survey (eFEDS), with a sky area of 140 square degrees with depth equivalent to the equatorial patch of the final eROSITA all-sky survey, represents the largest continuous non-full-sky X-ray fields to-date, making it the premier data set for measuring the angular power spectrum. In this work, we measure the X-ray angular power spectrum of galaxy clusters and groups in the eFEDS field. We show that the measured power spectrum is consistent with past observations, including the ROSAT All Sky Survey, and the Chandra COSMOS and Bootes fields. The predictions of cluster gas halo model that is calibrated from Chandra observations is also consistent with the eFEDS power spectrum. While the eFEDS does not have large enough sky coverage to provide meaningful cosmological constraints, we predict that the X-ray power spectrum from the cycle 4 of the eROSITA All-Sky Survey (eRASS4) will provide constraints on $\Omega_M$ and $\sigma_8$ at the 10% level., Comment: 10 pages, 4 Figures, accepted for publication in MNRAS. Updated arxiv version matched to the published version. Major change: tensions at small-scale power are resolved after point source mask correction

The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$\alpha$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}., Comment: 18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at https://www.camel-simulations.org

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io., Comment: 17 pages, 1 figure. Third paper of a series of four. Hundreds of thousands of labeled 2D maps and 3D grids from thousands of simulated universes publicly available at https://camels-multifield-dataset.readthedocs.io

Accounting for the triaxial shapes of galaxy clusters will become important in the context of upcoming cosmological surveys. We show that, compared to the gas density distribution, the cluster gravitational potential can be better characterised by a simple 3D model and is more robust against fluctuations. Perturbations in the gas density distribution can have a substantial influence on the derived thermodynamic properties, while cluster potentials are smooth and well-approximated by a spheroidal model. We use a statistical sample of 85 galaxy clusters from a large cosmological hydrodynamical simulation to investigate cluster shapes as a function of radius. In particular, we examine the shape of isodensity and isopotential shells and analyze how it is affected by the choice of component (gas vs. potential), substructure removal (for the gas density) and the definition of the computation domain (interior vs. shells). We find that the orientation and axis ratios of gas isodensity contours are degenerate with the presence of substructures and unstable against fluctuations. We observe that, as the derived cluster shape depends on the method used for removing the substructures, thermodynamic properties extracted from e.g. X-ray emissivity profiles suffer from this additional, often underestimated bias. In contrast, the shape reconstruction of the potential is largely unaffected by these factors and converges towards simple geometric models for both relaxed and dynamically active clusters. The observation that cluster potentials are better represented by simple geometrical models and reconstructed with a low level of systematics for both dynamically active and relaxed clusters suggests that by characterising galaxy clusters by their potential rather than by their mass, dynamically active and relaxed clusters could be combined in cosmological studies, improving statistics and lowering scatter., Comment: Updated with referee's comments, revised version submitted to A&A, 16 pages, 14 figures

Recent advances in simulations and observations of galaxy clusters suggest that there exists a physical outer boundary of massive cluster-size dark matter haloes. In this work, we investigate the locations of the outer boundaries of dark matter and gas around cluster-size dark matter haloes, by analyzing a sample of 65 massive dark matter halos extracted from the Omega500 zoom-in hydrodynamical cosmological simulations. We show that the location of accretion shock is offset from that of the dark matter splashback radius, contrary to the prediction of the self-similar models. The accretion shock radius is larger than all definitions of the splashback radius in the literature by 20-100%. The accretion shock radius defined using the steepest drop in the entropy pressure profiles is approximately 2 times larger than the splashback radius defined by the steepest slope in the dark matter density profile, and it is ~1.2 times larger than the edge of the dark matter phase-space structure. We discuss implications of our results for multi-wavelength studies of galaxy clusters., Comment: 7 pages, 8 figures, accepted in MNRAS

The non-spherical shapes of dark matter and gas distributions introduce systematic uncertainties that affect observable-mass relations and selection functions of galaxy groups and clusters. However, the triaxial gas distributions depend on the non-linear physical processes of halo formation histories and baryonic physics, which are challenging to model accurately. In this study we explore a machine learning approach for modelling the dependence of gas shapes on dark matter and baryonic properties. With data from the IllustrisTNG hydrodynamical cosmological simulations, we develop a machine learning pipeline that applies \pkg{XGBoost}, an implementation of gradient boosted decision trees, to predict radial profiles of gas shapes from halo properties. We show that \pkg{XGBoost} models can accurately predict gas shape profiles in dark matter haloes. We also explore model interpretability with \pkg{SHAP}, a method that identifies the most predictive properties at different halo radii. We find that baryonic properties best predict gas shapes in halo cores, whereas dark matter shapes are the main predictors in the halo outskirts. This work demonstrates the power of interpretable machine learning in modelling observable properties of dark matter haloes in the era of multi-wavelength cosmological surveys., Comment: 14 pages, 8 figures, Accepted to MNRAS

The eROSITA X-ray telescope on board the Spectrum-Roentgen-Gamma (SRG) mission will measure the position and properties of about 100,000 clusters of galaxies and 3 million active galactic nuclei over the full sky. To study the statistical properties of this ongoing survey, it is key to estimate the selection function accurately. We create a set of full sky light-cones using the MultiDark and UNIT dark matter only N-body simulations. We present a novel method to predict the X-ray emission of galaxy clusters. Given a set of dark matter halo properties (mass, redshift, ellipticity, offset parameter), we construct an X-ray emissivity profile and image for each halo in the light-cone. We follow the eROSITA scanning strategy to produce a list of X-ray photons on the full sky. We predict scaling relations for the model clusters, which are in good agreement with the literature. The predicted number density of clusters as a function of flux also agrees with previous measurements. Finally, we obtain a scatter of 0.21 (0.07, 0.25) for the X-ray luminosity -- mass (temperature -- mass, luminosity -- temperature) model scaling relations. We provide catalogues with the model photons emitted by clusters and active galactic nuclei. These catalogues will aid the eROSITA end to end simulation flow analysis and in particular the source detection process and cataloguing methods., Comment: 15 pages, 11 figures, accepted in the open journal of astrophysics

The shape of dark matter haloes plays a critical role in constraining cosmology with upcoming large-scale structure surveys. In this paper, we study the correlations between the triaxial shapes and formation histories in dark matter haloes in the MultiDark Planck 2 N-body cosmological simulation. We find that halo ellipticity is strongly correlated with halo properties that serve as proxies of halo formation history, such as halo concentration and the peak-centroid offset. In particular, the correlation between halo ellipticity and halo concentration is nearly independent of the halo density peak height. We present a simple model for the correlation between halo ellipticity and concentration using conditional abundance matching, and provide fitting formulae for the multi-dimensional distributions of triaxial halo shape as a function of halo peak height. We apply our halo shape model to gauge the effects of halo ellipticity and orientation bias on the excess surface mass density profiles in cluster-size halos. Our model should be useful for exploring the impact of triaxial halo shape on cosmological constraints in upcoming weak lensing surveys of galaxy clusters., Comment: 9 pages, 7 figures. Matched to the published version on MNRAS

Upcoming multi-wavelength astronomical surveys will soon discover all massive galaxy clusters and provide unprecedented constraints on cosmology and cluster astrophysics. In this paper, we investigate the constraining power of the multi-band cluster surveys, through a joint analysis of three observables associated with clusters of galaxies, including thermal Sunyaev-Zel'dovich (tSZ) effect in cosmic microwave background (CMB), X-ray emission of ionized gas, and gravitational weak lensing effect of background galaxies by the cluster's gravitational potential. We develop a theoretical framework to predict and interpret two-point correlation statistics among the three observables using a semi-analytic model of intracluster medium (ICM) and halo-based approach. In this work, we show that the auto- and cross-angular power spectra in tSZ, X-ray and lensing statistics from upcoming missions (eROSITA, CMB-S4, and LSST) can help break the degeneracy between cosmology and ICM physics. These correlation statistics are less sensitive to selection biases, and are able to probe ICM physics in distant, faint and small clusters that are otherwise difficult to be detected individually. We show that the correlation statistics are able to provide cosmological constraints comparable to the conventional cluster abundance measurements, while constraining cluster astrophysics at the same time. Our results indicate that the correlation statistics can significantly enhance the scientific returns of upcoming multi-wavelength cluster surveys., Comment: 21 pages, 12 figures, 4 tables, accepted for publication in MNRAS

We use a statistical sample of galaxy clusters from a large cosmological $N$-body$+$hydrodynamics simulation to examine the relation between morphology, or shape, of the X-ray emitting intracluster medium (ICM) and the mass accretion history of the galaxy clusters. We find that the mass accretion rate (MAR) of a cluster is correlated with the ellipticity of the ICM. The correlation is largely driven by material accreted in the last $\sim 4.5$~Gyr, indicating a characteristic time-scale for relaxation of cluster gas. Furthermore, we find that the ellipticity of the outer regions ($R\sim R_{\rm 500c}$) of the ICM is correlated with the overall MAR of clusters, while ellipticity of the inner regions ($\lesssim 0.5 R_{\rm 500c}$) is sensitive to recent major mergers with mass ratios of $\geq 1:3$. Finally, we examine the impact of variations in cluster mass accretion history on the X-ray observable-mass scaling relations. We show that there is a {\it continuous\/} anti-correlation between the residuals in the $T_x-M$ relation and cluster MARs, within which merging and relaxed clusters occupy extremes of the distribution rather than form two peaks in a bi-modal distribution, as was often assumed previously. Our results indicate the systematic uncertainties in the X-ray observable-mass relations can be mitigated by using the information encoded in the apparent ICM ellipticity., Comment: 10 pages, 7 figures, comments welcome

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System., Comment: 5 pages + references; Submitted to the Astro2020 call for science white papers

Gravitational lensing is a powerful probe of the mass distribution of galaxy clusters and cosmology. However, accurate measurements of the cluster mass profiles are limited by uncertainties in cluster astrophysics. In this work, we present a physically motivated model of baryonic effects on the cluster mass profiles, which self-consistently takes into account the impact of baryons on the concentration as well as mass accretion histories of galaxy clusters. We calibrate this model using the Omega500 hydrodynamical cosmological simulations of galaxy clusters with varying baryonic physics. Our model will enable us to simultaneously constrain cluster mass, concentration, and cosmological parameters using stacked weak lensing measurements from upcoming optical cluster surveys., Comment: 12 pages, 7 figures, 2 tables, accepted for publication in MNRAS

The Hitomi X-ray satellite has provided the first direct measurements of the plasma velocity dispersion in a galaxy cluster. It finds a relatively "quiescent" gas with a line-of-sight velocity dispersion ~ 160 km/s, at 30 kpc to 60 kpc from the cluster center. This is surprising given the presence of jets and X-ray cavities that indicates on-going activity and feedback from the active galactic nucleus (AGN) at the cluster center. Using a set of mock Hitomi observations generated from a suite of state-of-the-art cosmological cluster simulations, and an isolated but higher resolution simulation of gas physics in the cluster core, including the effects of cooling and AGN feedback, we examine the likelihood of Hitomi detecting a cluster with the observed velocities. As long as the Perseus has not experienced a major merger in the last few gigayears, and AGN feedback is operating in a "gentle" mode, we reproduce the level of gas motions observed by Hitomi. The frequent mechanical AGN feedback generates net line-of-sight velocity dispersions ~100-200 km/s, bracketing the values measured in the Perseus core. The large-scale velocity shear observed across the core, on the other hand, is generated mainly by cosmic accretion such as mergers. We discuss the implications of these results for AGN feedback physics and cluster cosmology and progress that needs to be made in both simulations and observations, including a Hitomi re-flight and calorimeter-based instruments with higher spatial resolution., Comment: 10 pages, 3 figures, accepted for publication in ApJ; updated with minor changes to match the accepted version

Recent statistical X-ray measurements of the intracluster medium (ICM) indicate that gas temperature profiles in the outskirts of galaxy clusters deviate from self-similar evolution. Using a mass-limited sample of galaxy clusters from cosmological hydrodynamical simulations, we show that the departure from self-similarity can be explained by non-thermal gas motions driven by mergers and accretion. Contrary to previous claims, gaseous substructures only play a minor role in the temperature evolution in cluster outskirts. A careful choice of halo overdensity definition in self-similar scaling mitigates these departures. Our work highlights the importance of non-thermal gas motions in ICM evolution and the use of galaxy clusters as cosmological probes., Comment: 8 pages, 8 figures

The thermal Sunyaev-Zel'dovich (tSZ) effect signal is widely recognized as a robust mass proxy of galaxy clusters with small intrinsic scatter. However, recent observational calibration of the tSZ scaling relation using weak lensing (WL) mass exhibits considerably larger scatter than the intrinsic scatter predicted from numerical simulations. This raises a question as to whether we can realize the full statistical power of ongoing and upcoming tSZ-WL observations of galaxy clusters. In this work, we investigate the origin of observed scatter in the tSZ-WL scaling relation, using mock maps of galaxy clusters extracted from cosmological hydrodynamic simulations. We show that the inferred intrinsic scatter from mock tSZ-WL analyses is considerably larger than the intrinsic scatter measured in simulations, and comparable to the scatter in the observed tSZ-WL relation. We show that this enhanced scatter originates from the combination of the projection of correlated structures along the line of sight and the uncertainty in the cluster radius associated with WL mass estimates, causing the amplitude of the scatter to depend on the covariance between tSZ and WL signals. We present a statistical model to recover the unbiased cluster scaling relation and cosmological parameter by taking into account the covariance in the tSZ-WL mass relation from multi-wavelength cluster surveys., Comment: 13 pages, 10 figures, 3 tables, accepted for publication in MNRAS

Non-thermal pressure in galaxy clusters leads to underestimation of the mass of galaxy clusters based on hydrostatic equilibrium with thermal gas pressure. This occurs even for dynamically relaxed clusters that are used for calibrating the mass-observable scaling relations. We show that the analytical model for non-thermal pressure developed in Shi & Komatsu 2014 can correct for this so-called 'hydrostatic mass bias', if most of the non-thermal pressure comes from bulk and turbulent motions of gas in the intracluster medium. Our correction works for the sample average irrespective of the mass estimation method, or the dynamical state of the clusters. This makes it possible to correct for the bias in the hydrostatic mass estimates from X-ray surface brightness and the Sunyaev-Zel'dovich observations that will be available for clusters in a wide range of redshifts and dynamical states., Comment: 9 pages, 8 figures, published in MNRAS

Gas motions in galaxy clusters play important roles in determining the properties of the intracluster medium (ICM) and in the constraint of cosmological parameters via X-ray and Sunyaev-Zel'dovich effect observations of galaxy clusters. The Hitomi measurements of gas motions in the core of the Perseus Cluster have provided insights into the physics in galaxy clusters. The XARM mission, equipped with the Resolve X-ray micro-calorimeter, will continue Hitomi's legacy by measuring ICM motions through Doppler shifting and broadening of emission lines in a larger number of galaxy clusters, and at larger radii. In this work, we investigate how well we can measure bulk and turbulent gas motions in the ICM with XARM, by analyzing mock XARM simulations of galaxy clusters extracted from cosmological hydrodynamic simulations. We assess how photon counts, spectral fitting methods, multiphase ICM structure, deprojections, and region selection affect the measurements of gas motions. We first show that XARM is capable of recovering the underlying spherically averaged turbulent and bulk velocity profiles for dynamically relaxed clusters to within $\sim 50\%$ with a reasonable amount of photon counts in the X-ray emission lines. We also find that there are considerable azimuthal variations in the ICM velocities, where the velocities measured in a single azimuthal direction can significantly deviate from the true value even in dynamically relaxed systems. Such variation must be taken into account when interpreting data and developing observing strategies. We will discuss the prospect of using the upcoming XARM mission to measure non-thermal pressure and to correct for the hydrostatic mass bias of galaxy clusters. Our results are broadly applicable for future X-ray missions, such as Athena and Lynx., Comment: 14 pages, 12 figures, PASJ accepted (typos corrected)

Galaxy clusters exhibit remarkable self-similar behavior which allows us to establish simple scaling relationships between observable quantities and cluster masses, making galaxy clusters useful cosmological probes. Recent X-ray observations suggest that self-similarity may be broken in the outskirts of galaxy clusters. In this work, we analyze a mass-limited sample of massive galaxy clusters from the Omega500 cosmological hydrodynamic simulation to investigate the self-similarity of the diffuse X-ray emitting intracluster medium (ICM) in the outskirts of galaxy clusters. We find that the self-similarity of the outer ICM profiles is better preserved if they are normalized with respect to the mean density of the universe, while the inner profiles are more self-similar when normalized using the critical density. However, the outer ICM profiles as well as the location of accretion shock around clusters are sensitive to their mass accretion rate, which causes the apparent breaking of self-similarity in cluster outskirts. We also find that the collisional gas does not follow the distribution of collisionless dark matter perfectly in the infall regions of galaxy clusters, leading to 10% departures in the gas-to-dark matter density ratio from the cosmic mean value. Our results have a number implications for interpreting observations of galaxy clusters in X-ray and through the Sunyaev-Zel'dovich effect and their application to cluster cosmology., Comment: 12 pages, 10 figures. Updated to match published version on ApJ

The analysis of X-ray and Sunyaev-Zel\'dovich measurements of the intracluster medium (ICM) assumes that electrons are in thermal equilibrium with ions in the plasma. However, electron-ion equilibration timescales can be comparable to the Hubble time in the low density galaxy cluster outskirts, leading to differences between the electron and ion temperatures. This temperature difference can lead to systematic biases in cluster mass estimates and mass-observable scaling relations. To quantify the impact of non-equilibrium electrons on the ICM profiles in cluster outskirts, we use a high resolution cosmological simulation with a two-temperature model assuming the Spitzer equilibration timescale for the electrons. First, we show how the radial profile of this temperature bias depends on both the mass and mass accretion rate of the cluster; the bias is most pronounced in the most massive and most rapidly accreting clusters. For the most extreme case in our sample, we find that the bias is of order 10% at half of the cluster virial radius and increases to 40% at the edge of the cluster. We also find that gas in filaments is less susceptible to the non-equilibrium effect, leading to azimuthal variations at large cluster-centric radii. By analyzing mock Chandra observations of simulated clusters, we show that such azimuthal variations can be probed with deep X-ray observations. Finally, the mass-dependent temperature bias introduces biases in hydrostatic mass and cluster temperature, which has implications for cluster-based cosmological inferences. We provide a mass-dependent model for the temperature bias profile which can be useful for modeling the effect of electron-ion equilibration in galaxy clusters., Comment: 7 pages, 4 figures, Updated to match published version on ApJ

Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10% to 30%. The largest bias originates by both violation of hydrostatic equilibrium and an additional temperature bias caused by inhomogeneities in the X-ray emitting intra-cluster medium (ICM). To elucidate on this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle-hydrodynamics (SPH) and adaptive-mesh-refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in no-radiative simulations, while it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with the observed one in the central regions of clusters. Around R_500 the temperature inhomogeneities of the SPH simulations can generate twice the typical hydrostatic-equilibrium mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii., Comment: 13 pages, 12 figures, 4 tables

Cosmological constraints from X-ray and microwave observations of galaxy clusters are subjected to systematic uncertainties. Non-thermal pressure support due to internal gas motions in galaxy clusters is one of the major sources of astrophysical uncertainties. Using a mass-limited sample of galaxy clusters from a high-resolution hydrodynamical cosmological simulation, we characterize the non-thermal pressure fraction profile and study its dependence on redshift, mass, and mass accretion rate. We find that the non-thermal pressure fraction profile is universal across redshift when galaxy cluster radii are defined with respect to the mean matter density of the universe instead of the commonly used critical density. We also find that the non-thermal pressure is predominantly radial, and the gas velocity anisotropy profile exhibits strong universality when galaxy cluster radii are defined with respect to the mean matter density of the universe. However, we find that the non-thermal pressure fraction is strongly dependent on the mass accretion rate of the galaxy cluster. We provide fitting formulae for the universal non-thermal pressure fraction and velocity anisotropy profiles of gas in galaxy clusters, which should be useful in modeling astrophysical uncertainties pertinent to using galaxy clusters as cosmological probes., Comment: 9 pages, 6 figures; Accepted by ApJ

The study of galaxy cluster outskirts has emerged as one of the new frontiers in extragalactic astrophysics and cosmology with the advent of new observations in X-ray and microwave. However, the thermodynamic properties and chemical enrichment of this diffuse and azimuthally asymmetric component of the intracluster medium (ICM) are still not well understood. This work, for the first time, systematically explores potential observational biases in these regions. To assess X-ray measurements of galaxy cluster properties at large radii ($>{R}_{500c}$), we use mock Chandra analyses of cosmological galaxy cluster simulations. The pipeline is identical to that used for Chandra observations, but the biases discussed in this paper are relevant for all X-ray observations outside of ${R}_{500c}$. We find the following from our analysis: (1) filament regions can contribute as much as $50\%$ at $R_{200c}$ to the emission measure; (2) X-ray temperatures and metal abundances from model fitted mock X-ray spectra in a multi-temperature ICM respectively vary to the level of $10\%$ and $50\%$; (3) resulting density profiles vary to within $10\%$ out to $R_{200c}$, and gas mass, total mass, and baryon fractions all vary to within a few percent; (4) the bias from a metal abundance extrapolated a factor of five higher than the true metal abundance results in total mass measurements biased high by $20\%$ and total gas measurements biased low by $10\%$; and (5) differences in projection and dynamical state of a cluster can lead to gas density slope measurements that differ by a factor of $15\%$ and $30\%$, respectively. The presented results can partially account for some of the recent gas profile measurements in cluster outskirts by, e.g., Suzaku. Our findings are pertinent to future X-ray cosmological constraints from cluster outskirts., Comment: 13 pages, 11 figures, Accepted to ApJ

We use ROSAT All Sky Survey (RASS) broadband X-ray images and the optical clusters identified from SDSS DR7 to estimate the X-ray luminosities around $\sim 65,000$ candidate clusters with masses $\ga 10^{13}\msunh$ based on an Optical to X-ray (OTX) code we develop. We obtain a catalogue with X-ray luminosity for each cluster. This catalog contains 817 clusters (473 at redshift $z\le 0.12$) with $S/N> 3$ in X-ray detection. We find about $65\%$ of these X-ray clusters have their most massive member located near the X-ray flux peak; for the rest $35\%$, the most massive galaxy is separated from the X-ray peak, with the separation following a distribution expected from a NFW profile. We investigate a number of correlations between the optical and X-ray properties of these X-ray clusters, and find that: the cluster X-ray luminosity is correlated with the stellar mass (luminosity) of the clusters, as well as with the stellar mass (luminosity) of the central galaxy and the mass of the halo, but the scatter in these correlations is large. Comparing the properties of X-ray clusters of similar halo masses but having different X-ray luminosities, we find that massive halos with masses $\ga 10^{14}\msunh$ contain a larger fraction of red satellite galaxies when they are brighter in X-ray. ... A cluster catalog containing the optical properties of member galaxies and the X-ray luminosity is available at {\it http://gax.shao.ac.cn/data/Group.html}., Comment: 13 pages, 9 figures, accepted for publication in MNRAS