Source author record

Erwin T. Lau

Erwin T. Lau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2023arXiv

X-ray Absorption Lines in the Warm-Hot Intergalactic Medium: Probing Chandra observations with the CAMEL simulations

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.

preprint2022arXiv

Triaxiality in galaxy clusters: Mass versus Potential reconstructions

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.

preprint2021arXiv

SHAPing the Gas: Understanding Gas Shapes in Dark Matter Haloes with Interpretable Machine Learning

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.

preprint2021arXiv

The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence

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.

preprint2019arXiv

Probing Cosmology and Cluster Astrophysics with Multi-Wavelength Surveys I. Correlation Statistics

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.

preprint2016arXiv

Covariance in the Thermal SZ-Weak Lensing Mass Scaling Relation of Galaxy Clusters

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.

preprint2015arXiv

Analytical model for non-thermal pressure in galaxy clusters - III. Removing the hydrostatic mass bias

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.

preprint2015arXiv

Mass Accretion and its Effects on the Self-Similarity of Gas Profiles in the Outskirts of Galaxy Clusters

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.

preprint2015arXiv

Non-Equilibrium Electrons in the Outskirts of Galaxy Clusters

The analysis of X-ray and Sunyaev-Zeldovich 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.

preprint2014arXiv

Hydrodynamic Simulation of Non-thermal Pressure Profiles of Galaxy Clusters

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.

preprint2014arXiv

Testing X-ray Measurements of Galaxy Cluster Outskirts with Cosmological Simulations

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.

preprint2014arXiv

Weighing Galaxy Clusters with Gas. II. On the Origin of Hydrostatic Mass Bias in LambdaCDM Galaxy Clusters

The use of galaxy clusters as cosmological probes hinges on our ability to measure their masses accurately and with high precision. Hydrostatic mass is one of the most common methods for estimating the masses of individual galaxy clusters, which suffer from biases due to departures from hydrostatic equilibrium. Using a large, mass-limited sample of massive galaxy clusters from a high-resolution hydrodynamical cosmological simulation, in this work we show that in addition to turbulent and bulk gas velocities, acceleration of gas introduces biases in the hydrostatic mass estimate of galaxy clusters. In unrelaxed clusters, the acceleration bias is comparable to the bias due to non-thermal pressure associated with merger-induced turbulent and bulk gas motions. In relaxed clusters, the mean mass bias due to acceleration is small (<3%), but the scatter in the mass bias can be reduced by accounting for gas acceleration. Additionally, this acceleration bias is greater in the outskirts of higher redshift clusters where mergers are more frequent and clusters are accreting more rapidly. Since gas acceleration cannot be observed directly, it introduces an irreducible bias for hydrostatic mass estimates. This acceleration bias places limits on how well we can recover cluster masses from future X-ray and microwave observations. We discuss implications for cluster mass estimates based on X-ray, Sunyaev-Zeldovich effect, and gravitational lensing observations and their impact on cluster cosmology.

preprint2013arXiv

Measuring the X-ray luminosities of SDSS DR7 clusters from RASS

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}.

preprint2013arXiv

On the cross-section of Dark Matter using substructure infall into galaxy clusters

We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of systems whilst removing any inherent line-of-sight projections. In order to interpret this ratio as a cross-section of Dark Matter we derive an analytical description of sub-halo infall which encompasses; the force of the main cluster potential, the drag on a gas sub-halo, a model for Dark Matter self-interactions and the resulting sub-halo drag, the force on the gas and galaxies due to the Dark Matter sub-halo potential, and finally the buoyancy on the gas and Dark Matter. We create mock observations from cosmological simulations of structure formation and find that collisionless Dark Matter becomes physically separated from X-ray gas by up to 20h^-1 kpc. Adding realistic levels of noise, we are able to predict achievable constraints from observational data. Current archival data should be able to detect a difference in the dynamical behaviour of Dark Matter and standard model particles at 6 sigma, and measure the total interaction cross-section sigma/m with 68% confidence limits of +/- 1cm2g^-1. We note that this method is not restricted by the limited number of major merging events and is easily extended to large samples of clusters from future surveys which could potentially push statistical errors to 0.1cm^2g^-1.

preprint2013arXiv

Predicting Merger-Induced Gas Motions in Lambda-CDM Galaxy Clusters

In the hierarchical structure formation model, clusters of galaxies form through a sequence of mergers and continuous mass accretion, which generate significant random gas motions especially in their outskirts where material is actively accreting. Non-thermal pressure provided by the internal gas motions affects the thermodynamic structure of the X-ray emitting intracluster plasma and introduces biases in the physical interpretation of X-ray and Sunyaev-Zeldovich effect observations. However, we know very little about the nature of gas motions in galaxy clusters. The ASTRO-H X-ray mission, scheduled to launch in 2015, will have a calorimeter capable of measuring gas motions in galaxy clusters at the level of <100 km/s. In this work, we predict the level of merger-induced gas motions expected in the Lambda-CDM model using hydrodynamical simulations of galaxy cluster formation. We show that the gas velocity dispersion is larger in more massive clusters, but exhibits a large scatter. We show that systems with large gas motions are morphologically disturbed, while early forming, relaxed groups show a smaller level of gas motions. By analyzing mock ASTRO-H observations of simulated clusters, we show that such observations can accurately measure the gas velocity dispersion out to the outskirts of nearby relaxed galaxy clusters. ASTRO-H analysis of merging clusters, on the other hand, requires multi-component spectral fitting and enables unique studies of substructures in galaxy clusters by measuring both the peculiar velocities and the velocity dispersion of gas within individual sub-clusters.

preprint2013arXiv

Weighing Galaxy Clusters with Gas. I. On the Methods of Computing Hydrostatic Mass Bias

Mass estimates of galaxy clusters from X-ray and Sunyeav-Zel'dovich observations assume the intracluster gas is in hydrostatic equilibrium with their gravitational potential. However, since galaxy clusters are dynamically active objects whose dynamical states can deviate significantly from the equilibrium configuration, the departure from the hydrostatic equilibrium assumption is one of the largest sources of systematic uncertainties in cluster cosmology. In the literature there has been two methods for computing the hydrostatic mass bias based on the Euler and the modified Jeans equations, respectively, and there has been some confusion about the validity of these two methods. The word "Jeans" was a misnomer, which incorrectly implies that the gas is collisionless. To avoid further confusion, we instead refer these methods as "summation" and "averaging" methods respectively. In this work, we show that these two methods for computing the hydrostatic mass bias are equivalent by demonstrating that the equation used in the second method can be derived from taking spatial averages of the Euler equation. Specifically, we identify the correspondences of individual terms in these two methods mathematically and show that these correspondences are valid to within a few percent level using hydrodynamical simulations of galaxy cluster formation. In addition, we compute the mass bias associated with the acceleration of gas and show that its contribution is small in the virialized regions in the interior of galaxy cluster, but becomes non-negligible in the cluster outskirts of massive clusters where materials are being actively accreted. We discuss implications for understanding the origin of hydrostatic mass biases and for characterizing them using hydrodynamical simulations as well as observations.

preprint2012arXiv

Constraining Cluster Physics with the Shape of X-ray Clusters: Comparison of Local X-ray Clusters versus LCDM Clusters

Simulations of cluster formation have demonstrated that condensation of baryons into central galaxies during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder, even at radii as large as half of the virial radius. However, such simulations generally predict stellar fractions within cluster virial radii that are ~2 to 3 times larger than the stellar masses deduced from observations. In this work we compare ellipticity profiles of clusters simulated with and without baryonic cooling to the cluster ellipticity profiles derived from Chandra and ROSAT observations in an effort to constrain the fraction of gas that cools and condenses into the central galaxies within clusters. We find that the observed ellipticity profiles are fairly constant with radius, with an average ellipticity of 0.18 +/- 0.05. The observed ellipticity profiles are in good agreement with the predictions of non-radiative simulations. On the other hand, the ellipticity profiles of the clusters in simulations that include radiative cooling, star formation, and supernova feedback (but no AGN feedback) deviate significantly from the observed ellipticity profiles at all radii. The simulations with cooling overpredict (underpredict) ellipticity in the inner (outer) regions of galaxy clusters. By comparing the simulations with and without cooling, we show that the cooling of gas via cooling flows in the central regions of simulated clusters causes the gas distribution to be more oblate in the central regions, but makes the outer gas distribution more spherical. We find that late-time gas cooling and star formation are responsible for the significantly oblate gas distributions in cluster cores, but the gas shapes outside of cluster cores are set primarily by baryon dissipation at high redshift z > 2.

preprint2012arXiv

Cross identification between X-ray and Optical Clusters of Galaxies in the SDSS DR7 Field

We use the ROSAT all sky survey X-ray cluster catalogs and the optical SDSS DR7 galaxy and group catalogs to cross-identify X-ray clusters with their optical counterparts, resulting in a sample of 201 X-ray clusters in the sky coverage of SDSS DR7. We investigate various correlations between the optical and X-ray properties of these X-ray clusters, and find that the following optical properties are correlated with the X-ray luminosity: the central galaxy luminosity, the central galaxy mass, the characteristic group luminosity ($\propto \Lx^{0.43}$), the group stellar mass ($\propto \Lx^{0.46}$), with typical 1-$σ$ scatter of $\sim 0.67$ in $\log \Lx$. Using the observed number distribution of X-ray clusters, we obtain an unbiased scaling relation between the X-ray luminosity, the central galaxy stellar mass and the characteristic satellite stellar mass as ${\log L_X} = -0.26 + 2.90 [\log (M_{\ast, c} + 0.26 M_{\rm sat}) -12.0]$ (and in terms of luminosities, as ${\log L_X} = -0.15 + 2.38 [\log (L_{c} + 0.72 L_{\rm sat}) -12.0]$). We find that the systematic difference between different halo mass estimations, e.g., using the ranking of characteristic group stellar mass or using the X-ray luminosity scaling relation can be used to constrain cosmology. Comparing the properties of groups of similar stellar mass (or optical luminosities) and redshift that are X-ray luminous or under-luminous, we find that X-ray luminous groups have more faint satellite galaxies and higher red fraction in their satellites. The cross-identified X-ray clusters together with their optical properties are provided in Appendix B.

preprint2011arXiv

Characterizing Galaxy Clusters with Gravitational Potential

We propose a simple estimator for the gravitational potential of cluster-size halos using the temperature and density profiles of the intracluster gas based on the assumptions of hydrostatic equilibrium and spherical symmetry. Using high resolution cosmological simulations of galaxy clusters, we show that the scaling relation between this estimator and the gravitational potential has a small intrinsic scatter of ~8%-15%, and it is insensitive to baryon physics outside the cluster core. The slope and the normalization of the scaling relation vary weakly with redshift, and they are relatively independent of the choice of radial range used and the dynamical states of the clusters. The results presented here provide a possible way for using the cluster potential function as an alternative to the cluster mass function in constraining cosmology using galaxy clusters.

preprint2011arXiv

Shapes of Gas, Gravitational Potential and Dark Matter in Lambda-CDM Clusters

We present analysis of the three-dimensional shape of intracluster gas in clusters formed in cosmological simulations of the Lambda-CDM cosmology and compare it to the shape of dark matter distribution and the shape of the overall isopotential surfaces. We find that in simulations with radiative cooling, star formation and stellar feedback (CSF), intracluster gas outside the cluster core is more spherical compared to non-radiative (NR) simulations, while in the core the gas in the CSF runs is more triaxial and has a distinctly oblate shape. The latter reflects the ongoing cooling of gas, which settles into a thick oblate ellipsoid as it loses thermal energy. The shape of the gas in the inner regions of clusters can therefore be a useful diagnostic of gas cooling. We find that gas traces the shape of the underlying potential rather well outside the core, as expected in hydrostatic equilibrium. At smaller radii, however, the gas and potential shapes differ significantly. In the CSF runs, the difference reflects the fact that gas is partly rotationally supported. Interestingly, we find that in NR simulations the difference between gas and potential shape at small radii is due to random gas motions, which make the gas distribution more spherical than the equipotential surfaces. Finally, we use mock Chandra X-ray maps to show that the differences in shapes observed in three-dimensional distribution of gas are discernible in the ellipticity of X-ray isophotes. Contrasting the ellipticities measured in simulated clusters against observations can therefore constrain the amount of cooling of the intracluster medium and the presence of random gas motions in cluster cores.

preprint2010arXiv

Impact of Cluster Physics on the Sunyaev-Zel'dovich Power Spectrum

We use an analytic model to investigate the theoretical uncertainty on the thermal Sunyaev-Zel'dovich (SZ) power spectrum due to astrophysical uncertainties in the thermal structure of the intracluster medium. Our model accounts for star formation and energy feedback (from supernovae and active galactic nuclei) as well as radially dependent non-thermal pressure support due to random gas motions, the latter calibrated by recent hydrodynamical simulations. We compare the model against X-ray observations of low redshift clusters, finding excellent agreement with observed pressure profiles. Varying the levels of feedback and non-thermal pressure support can significantly change both the amplitude and shape of the thermal SZ power spectrum. Increasing the feedback suppresses power at small angular scales, shifting the peak of the power spectrum to lower ell. On the other hand, increasing the non-thermal pressure support has the opposite effect, significantly reducing power at large angular scales. In general, including non-thermal pressure at the level measured in simulations has a large effect on the power spectrum, reducing the amplitude by 50% at angular scales of a few arcminutes compared to a model without a non-thermal component. Our results demonstrate that measurements of the shape of the power spectrum can reveal useful information on important physical processes in groups and clusters, especially at high-redshift where there exists little observational data. Comparing with the recent South Pole Telescope measurements of the small-scale cosmic microwave background power spectrum, we find our model reduces the tension between the values of sigma_8 measured from the SZ power spectrum and from cluster abundances.

preprint2009arXiv

Effects of Baryon Dissipation on the Dark Matter Virial Scaling Relation

We investigate effects of baryon dissipation on the dark matter virial scaling relation between total mass and velocity dispersion and the velocity bias of galaxies in groups and clusters using self-consistent cosmological simulations. We show that the baryon dissipation increases the velocity dispersion of dark matter within the virial radius by 5% - 10%. The effect is mainly driven by the change in density and gravitational potential in inner regions of cluster, and it is larger in lower mass systems where gas cooling and star formation are more efficient. We also show that the galaxy velocity bias depends on how galaxies are selected. Galaxies selected based on their stellar mass exhibit no velocity bias, while galaxies selected based on their total mass show positive bias of ~10%, consistent with previous results based on collisionless dark matter- only simulations. We further find that observational estimates of galaxy velocity dispersion are unbiased with respect to the velocity dispersion of dark matter, provided galaxies are selected using their stellar masses and and their velocity dispersions are computed with more than twenty most massive galaxies. Velocity dispersions estimated with fewer galaxies, on the other hand, can lead to significant underestimate of dynamical masses. Results presented in this paper should be useful in interpretating high-redshift groups and clusters as well as cosmological constraints derived from upcoming optical cluster surveys.

preprint2009arXiv

Residual Gas Motions in the Intracluster Medium and Bias in Hydrostatic Measurements of Mass Profiles of Clusters

We present analysis of bulk and random gas motions in the intracluster medium using high-resolution Eulerian cosmological simulations of sixteen simulated clusters, including both very relaxed and unrelaxed systems and spanning a virial mass range of 5*10^13 - 2*10^15 Msun/h. We investigate effects of the residual subsonic gas motions on the hydrostatic estimates of mass profiles and concentrations of galaxy clusters. In agreement with previous studies we find that the gas motions contribute up to ~ 5%-15% of the total pressure support in relaxed clusters with contribution increasing with cluster-centric radius. The fractional pressure support is higher in unrelaxed systems. This contribution would not be accounted for in hydrostatic estimates of the total mass profile and would lead to systematic underestimate of mass. We demonstrate that total mass can be recovered accurately if pressure due to gas motions measured in simulations is explicitly taken into account in the equation of hydrostatic equilibrium. Given that the underestimate of mass is increasing at larger radii, where gas is less relaxed and contribution of gas motions to pressure is larger, the total density profile derived from hydrostatic analysis is more concentrated than the true profile. This may at least partially explain some high values of concentrations of clusters estimated from hydrostatic analysis of X-ray data.