Showing 1–47 of 47 results for author: Lau, E T

Comparison of Models for the Warm-Hot Circumgalactic Medium around Milky Way-like Galaxies

Authors: Priyanka Singh, Erwin T. Lau, Yakov Faerman, Jonathan Stern, Daisuke Nagai

Abstract: 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 hy… ▽ More 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. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Comments: 15 pages, 4 figures, 3 tables, accepted for publication in MNRAS

Zooming by in the CARPoolGP lane: new CAMELS-TNG simulations of zoomed-in massive halos

Authors: Max E. Lee, Shy Genel, Benjamin D. Wandelt, Benjamin Zhang, Ana Maria Delgado, Shivam Pandey, Erwin T. Lau, Christopher Carr, Harrison Cook, Daisuke Nagai, Daniel Angles-Alcazar, Francisco Villaescusa-Navarro, Greg L. Bryan

Abstract: 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 dev… ▽ More 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. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: The manuscript was submitted to arxiv after receiving and responding to comments from the first referee report

Impact of Property Covariance on Cluster Weak lensing Scaling Relations

Authors: Zhuowen Zhang, Arya Farahi, Daisuke Nagai, Erwin T. Lau, Joshua Frieman, Marina Ricci, Anja von der Linden, Hao-yi Wu

Abstract: 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, $ΔΣ$, 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 thi… ▽ More 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, $ΔΣ$, 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 $ΔΣ$ 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}$-$ΔΣ$ 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. △ Less

Submitted 6 June, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Mapping the Intracluster Medium in the Era of High-resolution X-ray Spectroscopy

Authors: Congyao Zhang, Irina Zhuravleva, Maxim Markevitch, John ZuHone, François Mernier, Veronica Biffi, Ákos Bogdán, Priyanka Chakraborty, Eugene Churazov, Klaus Dolag, Stefano Ettori, William R. Forman, Christine Jones, Ildar Khabibullin, Caroline Kilbourne, Ralph Kraft, Erwin T. Lau, Sheng-Chieh Lin, Daisuke Nagai, Dylan Nelson, Anna Ogorzałek, Elena Rasia, Arnab Sarkar, Aurora Simionescu, Yuanyuan Su , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 24 pages, 26 figures, submitted to MNRAS. Comments are welcome

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

Authors: Amanda Butler Contreras, Erwin T. Lau, Benjamin D. Oppenheimer, Ákos Bogdán, Megan Tillman, Daisuke Nagai, Orsolya E. Kovács, Blakesley Burkhart

Abstract: 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 qua… ▽ More 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. △ Less

Submitted 12 January, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 11 pages, 9 figures. Matching to the published version on MNRAS

Journal ref: MNRAS, Volume 519, Issue 2, February 2023, Pages 2251-2261

The X-ray Angular Power Spectrum of Extended Sources in the eROSITA Final Equatorial Depth Survey

Authors: Erwin T. Lau, Akos Bogdan, Urmila Chadayammuri, Daisuke Nagai, Ralph Kraft, Nico Cappelluti

Abstract: 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 group… ▽ More 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 $Ω_M$ and $σ_8$ at the 10% level. △ Less

Submitted 14 November, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: 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 CAMELS project: public data release

Authors: Francisco Villaescusa-Navarro, Shy Genel, Daniel Anglés-Alcázar, Lucia A. Perez, Pablo Villanueva-Domingo, Digvijay Wadekar, Helen Shao, Faizan G. Mohammad, Sultan Hassan, Emily Moser, Erwin T. Lau, Luis Fernando Machado Poletti Valle, Andrina Nicola, Leander Thiele, Yongseok Jo, Oliver H. E. Philcox, Benjamin D. Oppenheimer, Megan Tillman, ChangHoon Hahn, Neerav Kaushal, Alice Pisani, Matthew Gebhardt, Ana Maria Delgado, Joyce Caliendo, Christina Kreisch , et al. (22 additional authors not shown)

Abstract: 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… ▽ More 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-$α$ 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}. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: 18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at https://www.camel-simulations.org

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

Authors: Francisco Villaescusa-Navarro, Shy Genel, Daniel Angles-Alcazar, Leander Thiele, Romeel Dave, Desika Narayanan, Andrina Nicola, Yin Li, Pablo Villanueva-Domingo, Benjamin Wandelt, David N. Spergel, Rachel S. Somerville, Jose Manuel Zorrilla Matilla, Faizan G. Mohammad, Sultan Hassan, Helen Shao, Digvijay Wadekar, Michael Eickenberg, Kaze W. K. Wong, Gabriella Contardo, Yongseok Jo, Emily Moser, Erwin T. Lau, Luis Fernando Machado Poletti Valle, Lucia A. Perez , et al. (3 additional authors not shown)

Abstract: 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 year… ▽ More 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. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 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

Triaxiality in galaxy clusters: Mass versus Potential reconstructions

Authors: Sebastian Stapelberg, Céline Tchernin, Damaris Hug, Erwin T. Lau, Matthias Bartelmann

Abstract: 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… ▽ More 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. △ Less

Submitted 2 February, 2022; v1 submitted 24 December, 2020; originally announced December 2020.

Comments: Updated with referee's comments, revised version submitted to A&A, 16 pages, 14 figures

Journal ref: A&A 663, A17 (2022)

Shock and Splash: Gas and Dark Matter Halo Boundaries around LambdaCDM Galaxy Clusters

Authors: Han Aung, Daisuke Nagai, Erwin T. Lau

Abstract: 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 co… ▽ More 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. △ Less

Submitted 22 October, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: 7 pages, 8 figures, accepted in MNRAS

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

Authors: Luis Fernando Machado Poletti Valle, Camille Avestruz, David J. Barnes, Arya Farahi, Erwin T. Lau, Daisuke Nagai

Abstract: 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 machin… ▽ More 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. △ Less

Submitted 2 August, 2021; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 14 pages, 8 figures, Accepted to MNRAS

Full-sky photon simulation of clusters and active galactic nuclei in the soft X-rays for eROSITA

Authors: Johan Comparat, Dominique Eckert, Alexis Finoguenov, Robert Schmidt, Jeremy Sanders, Daisuke Nagai, Erwin T. Lau, Florian Kaefer, Florian Pacaud, Nicolas Clerc, Thomas H. Reiprich, Esra Bulbul, Jacob Ider Chitham, Chia-Hsun Chuang, Vittorio Ghirardini, Violeta Gonzalez-Perez, Ghassem Gozaliazl, Charles C. Kirkpatrick, Anatoly Klypin, Andrea Merloni, Kirpal Nandra, Teng Liu, Francisco Prada, Miriam E. Ramos-Ceja, Mara Salvato , et al. (3 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 7 December, 2020; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 15 pages, 11 figures, accepted in the open journal of astrophysics

Journal ref: Volume 3, id 13 2020 in the Open Journal of Astrophysics

Characterizing galaxy clusters by their gravitational potential: systematics of cluster potential reconstruction

Authors: C. Tchernin, E. T. Lau, S. Stapelberg, D. Hug, M. Bartelmann

Abstract: Context. Biases in mass measurements of galaxy clusters are one of the major limiting systematics in constraining cosmology with clusters. Aims. We aim to demonstrate that the systematics associated with cluster gravitational potentials are smaller than the hydrostatic mass bias and that cluster potentials could therefore be a good alternative to cluster masses in cosmological studies. Methods. Us… ▽ More Context. Biases in mass measurements of galaxy clusters are one of the major limiting systematics in constraining cosmology with clusters. Aims. We aim to demonstrate that the systematics associated with cluster gravitational potentials are smaller than the hydrostatic mass bias and that cluster potentials could therefore be a good alternative to cluster masses in cosmological studies. Methods. Using cosmological simulations of galaxy clusters, we compute the biases in the hydrostatic mass (HE mass) and those in the gravitational potential, reconstructed from measurements at X-ray and millimeter wavelengths. In particular, we investigate the effects of the presence of substructures and of non-thermal pressure support on both the HE mass and the reconstructed potential. Results. We find that the bias in the reconstructed potential (6%) is less than that of the HE mass (13%), and that the scatter in the reconstructed potential decreases by about 35% with respect to that in the HE mass. Conclusions. This study shows that characterizing galaxy clusters by their gravitational potential is a promising alternative to using cluster masses in cluster cosmology. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: submitted to the journal A&A, 16 pages, 26 figures

Journal ref: A&A 644, A126 (2020)

Correlations between Triaxial Shapes and Formation History of Dark Matter Haloes

Authors: Erwin T. Lau, Andrew P. Hearin, Daisuke Nagai, Nico Cappelluti

Abstract: 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 hal… ▽ More 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. △ Less

Submitted 27 October, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 9 pages, 7 figures. Matched to the published version on MNRAS

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

Authors: Masato Shirasaki, Erwin T. Lau, Daisuke Nagai

Abstract: 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 cosm… ▽ More 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. △ Less

Submitted 24 October, 2019; v1 submitted 4 September, 2019; originally announced September 2019.

Comments: 21 pages, 12 figures, 4 tables, accepted for publication in MNRAS

Imprints of Mass Accretion History on the Shape of the Intracluster Medium and the $T_X-M$ Relation

Authors: Huanqing Chen, Camille Avestruz, Andrey V. Kravtsov, Erwin T. Lau, Daisuke Nagai

Abstract: 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 driv… ▽ More 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. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: 10 pages, 7 figures, comments welcome

Unveiling the Galaxy Cluster - Cosmic Web Connection with X-ray observations in the Next Decade

Authors: Stephen A. Walker, Daisuke Nagai, A. Simionescu, M. Markevitch, H. Akamatsu, M. Arnaud, C. Avestruz, M. Bautz, V. Biffi, S. Borgani, E. Bulbul, E. Churazov, K. Dolag, D. Eckert, S. Ettori, Y. Fujita, M. Gaspari, V. Ghirardini, R. Kraft, E. T. Lau, A. Mantz, K. Matsushita, M. McDonald, E. Miller, T. Mroczkowski , et al. (13 additional authors not shown)

Abstract: In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect… ▽ More In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 10 pages, 3 figures, Science white paper submitted to the Astro2020 Decadal Survey

Science from an Ultra-Deep, High-Resolution Millimeter-Wave Survey

Authors: Neelima Sehgal, Ho Nam Nguyen, Joel Meyers, Moritz Munchmeyer, Tony Mroczkowski, Luca Di Mascolo, Eric Baxter, Francis-Yan Cyr-Racine, Mathew Madhavacheril, Benjamin Beringue, Gil Holder, Daisuke Nagai, Simon Dicker, Cora Dvorkin, Simone Ferraro, George M. Fuller, Vera Gluscevic, Dongwon Han, Bhuvnesh Jain, Bradley Johnson, Pamela Klaassen, Daan Meerburg, Pavel Motloch, David N. Spergel, Alexander van Engelen , et al. (44 additional authors not shown)

Abstract: 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 woul… ▽ More 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. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: 5 pages + references; Submitted to the Astro2020 call for science white papers

Modelling Baryonic Effects on Galaxy Cluster Mass Profiles

Authors: Masato Shirasaki, Erwin T. Lau, Daisuke Nagai

Abstract: 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 concentra… ▽ More 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. △ Less

Submitted 28 March, 2018; v1 submitted 16 November, 2017; originally announced November 2017.

Comments: 12 pages, 7 figures, 2 tables, accepted for publication in MNRAS

Physical Origins of Gas Motions in Galaxy Cluster Cores: Interpreting Hitomi Observations of the Perseus Cluster

Authors: Erwin T. Lau, Massimo Gaspari, Daisuke Nagai, Paolo Coppi

Abstract: 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 nuc… ▽ More 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. △ Less

Submitted 18 September, 2017; v1 submitted 17 May, 2017; originally announced May 2017.

Comments: 10 pages, 3 figures, accepted for publication in ApJ; updated with minor changes to match the accepted version

The XMM Cluster Outskirts Project (X-COP): Physical conditions to the virial radius of Abell 2142

Authors: C. Tchernin, D. Eckert, S. Ettori, E. Pointecouteau, S. Paltani, S. Molendi, G. Hurier, F. Gastaldello, E. T. Lau, D. Nagai, M. Roncarelli, M. Rossetti

Abstract: Context. Galaxy clusters are continuously growing through the accretion of matter in their outskirts. This process induces inhomogeneities in the gas density distribution (clumping) which need to be taken into account to recover the physical properties of the intracluster medium (ICM) at large radii. Aims. We studied the thermodynamic properties in the outskirts (R > R500) of the massive galaxy cl… ▽ More Context. Galaxy clusters are continuously growing through the accretion of matter in their outskirts. This process induces inhomogeneities in the gas density distribution (clumping) which need to be taken into account to recover the physical properties of the intracluster medium (ICM) at large radii. Aims. We studied the thermodynamic properties in the outskirts (R > R500) of the massive galaxy cluster Abell 2142 by combining the Sunyaev Zel'dovich (SZ) effect with the X-ray signal. Methods. We combined the SZ pressure profile measured by Planck with the XMM-Newton gas density profile to recover radial profiles of temperature, entropy and hydrostatic mass out to 2R500. We used a method that is insensitive to clumping to recover the gas density, and we compared the results with traditional X-ray measurement techniques. Results. When taking clumping into account, our joint SZ/X-ray entropy profile is consistent with the predictions from pure gravitational collapse, whereas a significant entropy flattening is found when the effect of clumping is neglected. The hydrostatic mass profile recovered using joint X-ray/SZ data agrees with that obtained from spectroscopic X-ray measurements and with mass reconstructions obtained through weak lensing and galaxy kinematics. Conclusions. We found that clumping can explain the entropy flattening observed by Suzaku in the outskirts of several clusters. When using a method insensitive to clumping for the reconstruction of the gas density, the thermodynamic properties of Abell 2142 are compatible with the assumption that the thermal gas pressure sustains gravity and that the entropy is injected at accretion shocks, with no need to evoke more exotic physics. Our results highlight the need for X-ray observations with sufficient spatial resolution, and large collecting area, to understand the processes at work in cluster outer regions. △ Less

Submitted 17 June, 2016; originally announced June 2016.

Comments: 22 pages, 32 figures, accepted in the journal A&A

Journal ref: A&A 595, A42 (2016)

Stirred, not Clumped: Evolution of Temperature Profiles in the Outskirts of Galaxy Clusters

Authors: Camille Avestruz, Daisuke Nagai, Erwin T. Lau

Abstract: 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.… ▽ More 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. △ Less

Submitted 27 August, 2019; v1 submitted 5 May, 2016; originally announced May 2016.

Comments: 8 pages, 8 figures

Journal ref: The Astrophysical Journal, Volume 833, Issue 2, article id. 227, 9 pp. (2016)

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

Authors: Masato Shirasaki, Daisuke Nagai, Erwin T. Lau

Abstract: 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… ▽ More 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. △ Less

Submitted 23 May, 2016; v1 submitted 28 March, 2016; originally announced March 2016.

Comments: 13 pages, 10 figures, 3 tables, accepted for publication in MNRAS

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

Authors: Xun Shi, Eiichiro Komatsu, Daisuke Nagai, Erwin T. Lau

Abstract: 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 'hy… ▽ More 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. △ Less

Submitted 1 December, 2015; v1 submitted 15 July, 2015; originally announced July 2015.

Comments: 9 pages, 8 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2016 455 (3): 2936-2944

Constraining hydrostatic mass bias of galaxy clusters with high-resolution X-ray spectroscopy

Authors: Naomi Ota, Daisuke Nagai, Erwin T. Lau

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 27 June, 2018; v1 submitted 9 July, 2015; originally announced July 2015.

Comments: 14 pages, 12 figures, PASJ accepted (typos corrected)

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

Authors: Erwin T. Lau, Daisuke Nagai, Camille Avestruz, Kaylea Nelson, Alexey Vikhlinin

Abstract: 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… ▽ More 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. △ Less

Submitted 13 June, 2015; v1 submitted 19 November, 2014; originally announced November 2014.

Comments: 12 pages, 10 figures. Updated to match published version on ApJ

Journal ref: The Astrophysical Journal 2015, 806(1), 68

Non-Equilibrium Electrons in the Outskirts of Galaxy Clusters

Authors: Camille Avestruz, Daisuke Nagai, Erwin T. Lau, Kaylea Nelson

Abstract: 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… ▽ More 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. △ Less

Submitted 19 August, 2015; v1 submitted 29 October, 2014; originally announced October 2014.

Comments: 7 pages, 4 figures, Updated to match published version on ApJ

Journal ref: 2015, ApJ, 808, 176

Temperature Structure of the Intra-Cluster Medium from SPH and AMR simulations

Authors: Elena Rasia, Erwin T. Lau, Stefano Borgani, Daisuke Nagai, Klaus Dolag, Camille Avestruz, Gian Luigi Granato, Pasquale Mazzotta, Giuseppe Murante, Kaylea Nelson, Cinthia Ragone-Figueroa

Abstract: 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… ▽ More 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. △ Less

Submitted 7 August, 2014; v1 submitted 17 June, 2014; originally announced June 2014.

Comments: 13 pages, 12 figures, 4 tables

Journal ref: 2014, ApJ, 791, 96

The relation between gas density and velocity power spectra in galaxy clusters: qualitative treatment and cosmological simulations

Authors: I. Zhuravleva, E. Churazov, A. A. Schekochihin, E. T. Lau, D. Nagai, M. Gaspari, S. W. Allen, K. Nelson, I. J. Parrish

Abstract: We address the problem of evaluating the power spectrum of the velocity field of the ICM using only information on the plasma density fluctuations, which can be measured today by Chandra and XMM-Newton observatories. We argue that for relaxed clusters there is a linear relation between the rms density and velocity fluctuations across a range of scales, from the largest ones, where motions are domi… ▽ More We address the problem of evaluating the power spectrum of the velocity field of the ICM using only information on the plasma density fluctuations, which can be measured today by Chandra and XMM-Newton observatories. We argue that for relaxed clusters there is a linear relation between the rms density and velocity fluctuations across a range of scales, from the largest ones, where motions are dominated by buoyancy, down to small, turbulent scales: $(δρ_k/ρ)^2 = η_1^2 (V_{1,k}/c_s)^2$, where $δρ_k/ρ$ is the spectral amplitude of the density perturbations at wave number $k$, $V_{1,k}^2=V_k^2/3$ is the mean square component of the velocity field, $c_s$ is the sound speed, and $η_1$ is a dimensionless constant of order unity. Using cosmological simulations of relaxed galaxy clusters, we calibrate this relation and find $η_1\approx 1 \pm 0.3$. We argue that this value is set at large scales by buoyancy physics, while at small scales the density and velocity power spectra are proportional because the former are a passive scalar advected by the latter. This opens an interesting possibility to use gas density power spectra as a proxy for the velocity power spectra in relaxed clusters, across a wide range of scales. △ Less

Submitted 21 April, 2014; originally announced April 2014.

Comments: 6 pages, 3 figures, submitted to ApJ Letters

The relation between gas density and velocity power spectra in galaxy clusters: high-resolution hydrodynamic simulations and the role of conduction

Authors: M. Gaspari, E. Churazov, D. Nagai, E. T. Lau, I. Zhuravleva

Abstract: Exploring the ICM power spectrum can help us to probe the physics of galaxy clusters. Using high-resolution 3D plasma simulations, we study the statistics of the velocity field and its relation with the thermodynamic perturbations. The normalization of the ICM spectrum (density, entropy, or pressure) is linearly tied to the level of large-scale motions, which excite both gravity and sound waves du… ▽ More Exploring the ICM power spectrum can help us to probe the physics of galaxy clusters. Using high-resolution 3D plasma simulations, we study the statistics of the velocity field and its relation with the thermodynamic perturbations. The normalization of the ICM spectrum (density, entropy, or pressure) is linearly tied to the level of large-scale motions, which excite both gravity and sound waves due to stratification. For low 3D Mach number M~0.25, gravity waves mainly drive entropy perturbations, traced by preferentially tangential turbulence. For M>0.5, sound waves start to significantly contribute, passing the leading role to compressive pressure fluctuations, associated with isotropic (or slightly radial) turbulence. Density and temperature fluctuations are then characterized by the dominant process: isobaric (low M), adiabatic (high M), or isothermal (strong conduction). Most clusters reside in the intermediate regime, showing a mixture of gravity and sound waves, hence drifting towards isotropic velocities. Remarkably, regardless of the regime, the variance of density perturbations is comparable to the 1D Mach number. This linear relation allows to easily convert between gas motions and ICM perturbations, which can be exploited by Chandra, XMM data and by the forthcoming Astro-H. At intermediate and small scales (10-100 kpc), the turbulent velocities develop a Kolmogorov cascade. The thermodynamic perturbations act as effective tracers of the velocity field, broadly consistent with the Kolmogorov-Obukhov-Corrsin advection theory. Thermal conduction acts to damp the gas fluctuations, washing out the filamentary structures and steepening the spectrum, while leaving unaltered the velocity cascade. The ratio of the velocity and density spectrum thus inverts the downtrend shown by the non-diffusive models, allowing to probe the presence of significant conductivity in the ICM. △ Less

Submitted 7 July, 2014; v1 submitted 21 April, 2014; originally announced April 2014.

Comments: Accepted by A&A; 15 pages, 10 figures; added insights and references - thank you for the positive feedback

Hydrodynamic Simulation of Non-thermal Pressure Profiles of Galaxy Clusters

Authors: Kaylea Nelson, Erwin T. Lau, Daisuke Nagai

Abstract: 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-therm… ▽ More 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. △ Less

Submitted 16 August, 2014; v1 submitted 17 April, 2014; originally announced April 2014.

Comments: 9 pages, 6 figures; Accepted by ApJ

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

Authors: Camille Avestruz, Erwin T. Lau, Daisuke Nagai, Alexey Vikhlinin

Abstract: 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, systemat… ▽ More 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. △ Less

Submitted 23 October, 2014; v1 submitted 17 April, 2014; originally announced April 2014.

Comments: 13 pages, 11 figures, Accepted to ApJ

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

Authors: Lei Wang, Xiaohu Yang, Shiyin Shen, H. J. Mo, Frank C. van den Bosch, Wentao Luo, Yu Wang, Erwin T. Lau, Q. D. Wang, Xi Kang, Ran Li

Abstract: 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… ▽ More 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}. △ Less

Submitted 28 December, 2013; originally announced December 2013.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

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

Authors: David Harvey, Eric Tittley, Richard Massey, Thomas D. Kitching, Andy Taylor, Simon R. Pike, Scott T. Kay, Erwin T. Lau, Daisuke Nagai

Abstract: 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… ▽ More 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. △ Less

Submitted 20 November, 2013; v1 submitted 7 October, 2013; originally announced October 2013.

Comments: 14 pages, 11 figures

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

Authors: Kaylea Nelson, Erwin T. Lau, Daisuke Nagai, Douglas H. Rudd, Liang Yu

Abstract: 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 hy… ▽ More 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. △ Less

Submitted 10 February, 2014; v1 submitted 29 August, 2013; originally announced August 2013.

Comments: 9 pages, 6 figures, 1 table; Updated to match version accepted by ApJ

Journal ref: ApJ, 182 (2014), 107

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

Authors: Daisuke Nagai, Erwin T. Lau, Camille Avestruz, Kaylea Nelson, Douglas H. Rudd

Abstract: 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 b… ▽ More 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. △ Less

Submitted 1 October, 2013; v1 submitted 8 July, 2013; originally announced July 2013.

Comments: 8 pages, 7 figures, accepted for publication in ApJ

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

Authors: Erwin T. Lau, Daisuke Nagai, Kaylea Nelson

Abstract: 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 larges… ▽ More 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. △ Less

Submitted 1 October, 2013; v1 submitted 17 June, 2013; originally announced June 2013.

Comments: 6 pages, 2 figures, accepted for publication in ApJ

Bias from gas inhomogeneities in the pressure profiles as measured from X-ray and SZ observations

Authors: S. Khedekar, E. Churazov, A. Kravtsov, I. Zhuravleva, E. T. Lau, D. Nagai, R. Sunyaev

Abstract: X-ray observations of galaxy clusters provide emission measure weighted spectra, arising from a range of density and temperature fluctuations in the intra-cluster medium (ICM). This is fitted to a single temperature plasma emission model to provide an estimate of the gas density and temperature, which are sensitive to the gas inhomogeneities. Therefore, X-ray observations yield a potentially biase… ▽ More X-ray observations of galaxy clusters provide emission measure weighted spectra, arising from a range of density and temperature fluctuations in the intra-cluster medium (ICM). This is fitted to a single temperature plasma emission model to provide an estimate of the gas density and temperature, which are sensitive to the gas inhomogeneities. Therefore, X-ray observations yield a potentially biased estimate of the thermal gas pressure, P_X. At the same time Sunyaev-Zeldovich (SZ) observations directly measure the integrated gas pressure, P_SZ. If the X-ray pressure profiles are strongly biased with respect to the SZ, then one has the possibility to probe the gas inhomogeneities, even at scales unresolved by the current generation of telescopes. At the same time, a weak bias has implications for the interchangeable use of mass proxies like Y_SZ and Y_X as cosmological probes. In this paper we investigate the dependence of the bias, defined as b_P(r)=P_X(r)/P_SZ(r)-1, on the characteristics of fluctuations in the ICM taking into account the correlation between temperature and density fluctuations. We made a simple prediction of the irreducible bias in idealised X-ray vs SZ observations using multi-temperature plasma emission model. We also provide a simple fitting form to estimate the bias given the distribution of fluctuations. Analysing a sample of 16 clusters extracted from hydrodynamical simulations, we find that the median value of bias is within +/-3% within R_500, it decreases to -5% at R_500<r<1.5R_500 and then rises back to ~0% at r >~ 2R_500. The scatter of b_P(r) between individual relaxed clusters is small -- at the level of <0.03 within R_500, but turns significantly larger (0.25) and highly skewed at r>~ 1.5R_500. For any relaxed cluster we find |b_P(r)| < 15% within R_500, across different implementations of input physics in the simulations. [abridged] △ Less

Submitted 5 February, 2013; v1 submitted 14 November, 2012; originally announced November 2012.

Comments: 14 pages, 9 figures, 2 tables, accepted in MNRAS

Journal ref: MNRAS 2013

Quantifying properties of ICM inhomogeneities

Authors: I. Zhuravleva, E. Churazov, A. Kravtsov, E. T. Lau, D. Nagai, R. Sunyaev

Abstract: We present a new method to identify and characterize the structure of the intracluster medium (ICM) in simulated galaxy clusters. The method uses the median of gas properties, such as density and pressure, which we show to be very robust to the presence of gas inhomogeneities. In particular, we show that the radial profiles of median gas properties are smooth and do not exhibit fluctuations at loc… ▽ More We present a new method to identify and characterize the structure of the intracluster medium (ICM) in simulated galaxy clusters. The method uses the median of gas properties, such as density and pressure, which we show to be very robust to the presence of gas inhomogeneities. In particular, we show that the radial profiles of median gas properties are smooth and do not exhibit fluctuations at locations of massive clumps in contrast to mean and mode properties. It is shown that distribution of gas properties in a given radial shell can be well described by a log-normal PDF and a tail. The former corresponds to a nearly hydrostatic bulk component, accounting for ~99% of the volume, while the tail corresponds to high density inhomogeneities. We show that this results in a simple and robust separation of the diffuse and clumpy components of the ICM. The FWHM of the density distribution grows with radius and varies from ~0.15 dex in cluster centre to ~0.5 dex at 2r_500 in relaxed clusters. The small scatter in the width between relaxed clusters suggests that the degree of inhomogeneity is a robust characteristic of the ICM. It broadly agrees with the amplitude of density perturbations in the Coma cluster. We discuss the origin of ICM density variations in spherical shells and show that less than 20% of the width can be attributed to the triaxiality of the cluster gravitational potential. As a link to X-ray observations of real clusters we evaluated the ICM clumping factor with and without high density inhomogeneities. We argue that these two cases represent upper and lower limits on the departure of the observed X-ray emissivity from the median value. We find that the typical value of the clumping factor in the bulk component of relaxed clusters varies from ~1.1-1.2 at r_500 up to ~1.3-1.4 at r_200, in broad agreement with recent observations. △ Less

Submitted 24 October, 2012; originally announced October 2012.

Comments: 16 pages, 12 figure, accepted to MNRAS

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

Authors: Erwin T. Lau, Daisuke Nagai, Andrey V. Kravtsov, Alexey Vikhlinin, Andrew R. Zentner

Abstract: 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 stell… ▽ More 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. △ Less

Submitted 5 September, 2012; v1 submitted 10 January, 2012; originally announced January 2012.

Comments: 10 pages, 6 figures, matching the published version in ApJ. Corrected missing reference in the arxiv version

Journal ref: ApJ 2012, Vol 755, 116

The gas distribution in the outer regions of galaxy clusters

Authors: D. Eckert, F. Vazza, S. Ettori, S. Molendi, D. Nagai, E. T. Lau, M. Roncarelli, M. Rossetti, S. L. Snowden, F. Gastaldello

Abstract: We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC… ▽ More We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r200 and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond \sim r500. Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or treating gas clumping are in better agreement with the observed gas distribution. We report for the first time the high-confidence detection of a systematic difference between cool-core and non-cool core clusters beyond \sim 0.3r200, which we explain by a different distribution of the gas in the two classes. Beyond \sim r500, galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. [Abridged] △ Less

Submitted 7 October, 2015; v1 submitted 31 October, 2011; originally announced November 2011.

Comments: The data for the average profiles and individual clusters can be downloaded at: http://www.isdc.unige.ch/~deckert/newsite/The_Planck_ROSAT_project.html

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

Authors: Lei Wang, Xiaohu Yang, Wentao Luo, Erwin T. Lau, Yu Wang, H. J. Mo, Frank C. van den Bosch, Q. D. Wang

Abstract: 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 ar… ▽ More 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. △ Less

Submitted 9 January, 2012; v1 submitted 10 October, 2011; originally announced October 2011.

Comments: 22 pages, 13figures, 1 table, submitted to ApJ

Characterizing Galaxy Clusters with Gravitational Potential

Authors: Erwin T. Lau

Abstract: 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 intri… ▽ More 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. △ Less

Submitted 6 June, 2011; v1 submitted 10 September, 2010; originally announced September 2010.

Comments: 10 pages, 7 figures, 4 tables. Matching the version accepted for publication in the Astrophysical Journal

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

Authors: Laurie D. Shaw, Daisuke Nagai, Suman Bhattacharya, Erwin T. Lau

Abstract: 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, t… ▽ More 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. △ Less

Submitted 8 November, 2010; v1 submitted 10 June, 2010; originally announced June 2010.

Comments: 15 Pages, 9 Figures, updated to match version accepted by ApJ

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

Authors: Erwin T. Lau, Daisuke Nagai, Andrey V. Kravtsov, Andrew R. Zentner

Abstract: 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 s… ▽ More 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. △ Less

Submitted 20 April, 2011; v1 submitted 11 March, 2010; originally announced March 2010.

Comments: 11 pages, 8 figures, 3 tables, updated to match the version accepted for publication in the Astrophysical Journal

Effects of Baryon Dissipation on the Dark Matter Virial Scaling Relation

Authors: Erwin T. Lau, Daisuke Nagai, Andrey V. Kravtsov

Abstract: 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 chan… ▽ More 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. △ Less

Submitted 16 November, 2009; v1 submitted 14 August, 2009; originally announced August 2009.

Comments: 7 pages, 5 figures, accepted by ApJ for publication

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

Authors: Erwin T. Lau, Andrey V. Kravtsov, Daisuke Nagai

Abstract: 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 concentrati… ▽ More 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. △ Less

Submitted 17 September, 2009; v1 submitted 27 March, 2009; originally announced March 2009.

Comments: 11 pages, 10 figures, accepted by ApJ

Journal ref: Astrophys.J.705:1129-1138,2009