Showing 1–50 of 68 results for author: Lau, E

Masses of Sunyaev-Zel'dovich Galaxy Clusters Detected by The Atacama Cosmology Telescope: Stacked Lensing Measurements with Subaru HSC Year 3 data

Authors: Masato Shirasaki, Cristóbal Sifón, Hironao Miyatake, Erwin Lau, Zhuowen Zhang, Neta Bahcall, Mark Devlin, Jo Dunkley, Arya Farahi, Matt Hilton, Yen-Ting Lin, Daisuke Nagai, Suzanne T. Staggs, Tomomi Sunayama, David Spergel, Edward J. Wollack

Abstract: We present a stacked lensing analysis of 96 galaxy clusters selected by the thermal Sunyaev-Zel'dovich (SZ) effect in maps of the cosmic microwave background (CMB). We select foreground galaxy clusters with a $5σ$-level SZ threshold in CMB observations from the Atacama Cosmology Telescope, while we define background source galaxies for the lensing analysis with secure photometric redshift cuts in… ▽ More We present a stacked lensing analysis of 96 galaxy clusters selected by the thermal Sunyaev-Zel'dovich (SZ) effect in maps of the cosmic microwave background (CMB). We select foreground galaxy clusters with a $5σ$-level SZ threshold in CMB observations from the Atacama Cosmology Telescope, while we define background source galaxies for the lensing analysis with secure photometric redshift cuts in Year 3 data of the Subaru Hyper Suprime Cam survey. We detect the stacked lensing signal in the range of $0.1 < R\, [h^{-1}\mathrm{Mpc}] < 100$ in each of three cluster redshift bins, $0.092<z\le0.445$, $0.445<z\le0.695$, and $0.695<z\le1.180$, with 32 galaxy clusters in each bin. The cumulative signal-to-noise ratios of the lensing signal are $14.6$, $12.0$, and $6.6$, respectively. Using a halo-based forward model, we then constrain statistical relationships between the mass inferred from the SZ observation (i.e. SZ mass) and the total mass derived from our stacked lensing measurements. At the average SZ mass in the cluster sample ($2.1-2.4\times10^{14}\, h^{-1}M_\odot$), our likelihood analysis shows that the average total mass differs from the SZ counterpart by a factor of $1.3 \pm 0.2$, $1.6 \pm 0.2$, and $1.6 \pm 0.3$ ($68\%$) in the aforementioned redshift ranges, respectively. Our limits are consistent with previous lensing measurements, and we find that the cluster modeling choices can introduce a $1σ$-level difference in our parameter inferences. △ Less

Submitted 12 July, 2024; v1 submitted 11 July, 2024; originally announced July 2024.

Comments: 34 pages, 17 figures

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

OLIMPO: a Balloon-Borne SZE Imager to Probe ICM Dynamics and the WHIM

Authors: Jack Sayers, Camille Avestruz, Ritoban Basu Thakur, Elia Stefano Battistelli, Esra Bulbul, Federico Caccioti, Fabio Columbro, Alessandro Coppolecchia, Scott Cray, Giuseppe D'Alessandro, Paolo de Bernardis, Marco de Petris, Shaul Hanany, Luca Lamagna, Erwin Lau, Silvia Masi, Allesandro Paiella, Giorgio Pettinari, Francesco Piacentini, Eitan Rapaport, Larry Rudnick, Irina Zhuravleva, John ZuHuone

Abstract: OLIMPO is a proposed Antarctic balloon-borne Sunyaev-Zel'dovich effect (SZE) imager to study gas dynamics associated with structure formation along with the properties of the warm-hot intergalactic medium (WHIM) residing in the connective filaments. During a 25 day flight OLIMPO will image a total of 10 z~0.05 galaxy clusters and 8 bridges at 145, 250, 365, and 460 GHz at an angular resolution of… ▽ More OLIMPO is a proposed Antarctic balloon-borne Sunyaev-Zel'dovich effect (SZE) imager to study gas dynamics associated with structure formation along with the properties of the warm-hot intergalactic medium (WHIM) residing in the connective filaments. During a 25 day flight OLIMPO will image a total of 10 z~0.05 galaxy clusters and 8 bridges at 145, 250, 365, and 460 GHz at an angular resolution of 1.0'-3.3'. The maps will be significantly deeper than those planned from CMB-S4 and CCAT-P, and will have excellent fidelity to the large angular scales of our low-z targets, which are difficult to probe from the ground. In combination with X-ray data from eROSITA and XRISM we will transform our current static view of galaxy clusters into a full dynamic picture by measuring the internal intra-cluster medium (ICM) velocity structure with the kinematic SZE, X-ray spectroscopy, and the power spectrum of ICM fluctuations. Radio observations from ASKAP and MeerKAT will be used to better understand the connection between ICM turbulence and shocks with the relativistic plasma. Beyond the cluster boundary, we will combine thermal SZE maps from OLIMPO with X-ray imaging from eROSITA to measure the thermodynamics of the WHIM residing in filaments, providing a better understanding of its properties and its contribution to the total baryon budget. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: From the proceedings of the mm Universe 2023

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

Circumgalactic Medium on the Largest Scales: Detecting X-ray Absorption Lines with Large-Area Microcalorimeters

Authors: Akos Bogdan, Ildar Khabibullin, Orsolya Kovacs, Gerrit Schellenberger, John ZuHone, Joseph Burchett, Klaus Dolag, Eugene Churazov, William Forman, Christine Jones, Caroline Kilbourne, Ralph Kraft, Erwin Lau, Maxim Markevitch, Dan McCammon, Daisuke Nagai, Dylan Nelson, Anna Ogorzalek, Benjamin Oppenheimer, Arnab Sarkar, Yuanyuan Su, Nhut Truong, Sylvain Veilleux, Stephan Vladutescu-Zopp, Irina Zhuravleva

Abstract: The circumgalactic medium (CGM) plays a crucial role in galaxy evolution as it fuels star formation, retains metals ejected from the galaxies, and hosts gas flows in and out of galaxies. For Milky Way-type and more massive galaxies, the bulk of the CGM is in hot phases best accessible at X-ray wavelengths. However, our understanding of the CGM remains largely unconstrained due to its tenuous natur… ▽ More The circumgalactic medium (CGM) plays a crucial role in galaxy evolution as it fuels star formation, retains metals ejected from the galaxies, and hosts gas flows in and out of galaxies. For Milky Way-type and more massive galaxies, the bulk of the CGM is in hot phases best accessible at X-ray wavelengths. However, our understanding of the CGM remains largely unconstrained due to its tenuous nature. A promising way to probe the CGM is via X-ray absorption studies. Traditional absorption studies utilize bright background quasars, but this method probes the CGM in a pencil beam, and, due to the rarity of bright quasars, the galaxy population available for study is limited. Large-area, high spectral resolution X-ray microcalorimeters offer a new approach to exploring the CGM in emission and absorption. Here, we demonstrate that the cumulative X-ray emission from cosmic X-ray background sources can probe the CGM in absorption. We construct column density maps of major X-ray ions from the Magneticum simulation and build realistic mock images of nine galaxies to explore the detectability of X-ray absorption lines arising from the large-scale CGM. We conclude that the OVII absorption line is detectable around individual massive galaxies at the $3σ-6σ$ confidence level. For Milky Way-type galaxies, the OVII and OVIII absorption lines are detectable at the $\sim\,6σ$ and $\sim\,3σ$ levels even beyond the virial radius when co-adding data from multiple galaxies. This approach complements emission studies, does not require additional exposures, and will allow probing of the baryon budget and the CGM at the largest scales. △ Less

Submitted 8 June, 2023; originally announced June 2023.

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

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

Line Emission Mapper (LEM): Probing the physics of cosmic ecosystems

Authors: Ralph Kraft, Maxim Markevitch, Caroline Kilbourne, Joseph S. Adams, Hiroki Akamatsu, Mohammadreza Ayromlou, Simon R. Bandler, Marco Barbera, Douglas A. Bennett, Anil Bhardwaj, Veronica Biffi, Dennis Bodewits, Akos Bogdan, Massimiliano Bonamente, Stefano Borgani, Graziella Branduardi-Raymont, Joel N. Bregman, Joseph N. Burchett, Jenna Cann, Jenny Carter, Priyanka Chakraborty, Eugene Churazov, Robert A. Crain, Renata Cumbee, Romeel Dave , et al. (85 additional authors not shown)

Abstract: The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will answer the outstanding questions of the Universe's structure formation. It will also provide transformative new observing capabilities for every area of astrophysics, and to heliophysics and planetary physics as well. LEM's main goal is a comprehensive look at the physics of galaxy formation, including stellar and black-hole… ▽ More The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will answer the outstanding questions of the Universe's structure formation. It will also provide transformative new observing capabilities for every area of astrophysics, and to heliophysics and planetary physics as well. LEM's main goal is a comprehensive look at the physics of galaxy formation, including stellar and black-hole feedback and flows of baryonic matter into and out of galaxies. These processes are best studied in X-rays, and emission-line mapping is the pressing need in this area. LEM will use a large microcalorimeter array/IFU, covering a 30x30' field with 10" angular resolution, to map the soft X-ray line emission from objects that constitute galactic ecosystems. These include supernova remnants, star-forming regions, superbubbles, galactic outflows (such as the Fermi/eROSITA bubbles in the Milky Way and their analogs in other galaxies), the Circumgalactic Medium in the Milky Way and other galaxies, and the Intergalactic Medium at the outskirts and beyond the confines of galaxies and clusters. LEM's 1-2 eV spectral resolution in the 0.2-2 keV band will make it possible to disentangle the faintest emission lines in those objects from the bright Milky Way foreground, providing groundbreaking measurements of the physics of these plasmas, from temperatures, densities, chemical composition to gas dynamics. While LEM's main focus is on galaxy formation, it will provide transformative capability for all classes of astrophysical objects, from the Earth's magnetosphere, planets and comets to the interstellar medium and X-ray binaries in nearby galaxies, AGN, and cooling gas in galaxy clusters. In addition to pointed observations, LEM will perform a shallow all-sky survey that will dramatically expand the discovery space. △ Less

Submitted 12 April, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: 18 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity. v2: All-sky survey figure expanded, references fixed. v3: Added energy resolution measurements for prototype detector array. v4: Author list and reference fixes

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

Cluster outskirts and their connection to the cosmic web

Authors: Stephen Walker, Erwin Lau

Abstract: We review the latest developments in our X-ray observational and theoretical understanding of the outskirts of galaxy clusters, and their connection to the cosmic web. The faint cluster outskirts are challenging regions to observe in X-rays, requiring highly sensitive telescopes with low and stable background levels. We present our latest understanding of the thermodynamic profiles of clusters in… ▽ More We review the latest developments in our X-ray observational and theoretical understanding of the outskirts of galaxy clusters, and their connection to the cosmic web. The faint cluster outskirts are challenging regions to observe in X-rays, requiring highly sensitive telescopes with low and stable background levels. We present our latest understanding of the thermodynamic profiles of clusters in the outskirts, and the biases that gas clumping and non-thermal pressure support can introduce. Features in the outskirts due to merging activity are discussed, along with the chemical enrichment of the outskirts ICM. We describe future prospects for X-ray observations to explore further out in the cluster outskirts and probe their connections to the cosmic web. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 39 pages. Accepted for publication. This Chapter will appear in the Section "Galaxy Clusters" (Section Editors: E. Pointecouteau, E. Rasia, A. Simionescu) of the "Handbook of X-ray and Gamma-ray Astrophysics" (Editors in chief: C. Bambi and A. Santangelo)

The Circumgalactic Medium from the CAMELS Simulations: Forecasting Constraints on Feedback Processes from Future Sunyaev-Zeldovich Observations

Authors: Emily Moser, Nicholas Battaglia, Daisuke Nagai, Erwin Lau, Luis Fernando Machado Poletti Valle, Francisco Villaescusa-Navarro, Stefania Amodeo, Daniel Angles-Alcazar, Greg L. Bryan, Romeel Dave, Lars Hernquist, Mark Vogelsberger

Abstract: The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Lea… ▽ More The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS), that varies four different feedback parameters of two previously existing hydrodynamical simulations, IllustrisTNG and SIMBA. We capture the dependencies of SZ radial profiles on these feedback parameters with an emulator, calculate their derivatives, and forecast future constraints on these feedback parameters from upcoming experiments. We find that for a DESI-like (Dark Energy Spectroscopic Instrument) galaxy sample observed by the Simons Observatory all four feedback parameters are able to be constrained (some within the $10\%$ level), indicating that future observations will be able to further restrict the parameter space for these sub-grid models. Given the modeled galaxy sample and forecasted errors in this work, we find that the inner SZ profiles contribute more to the constraining power than the outer profiles. Finally, we find that, despite the wide range of AGN feedback parameter variation in the CAMELS simulation suite, we cannot reproduce the tSZ signal of galaxies selected by the Baryon Oscillation Spectroscopic Survey as measured by the Atacama Cosmology Telescope. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Comments: 23 pages, 7 figures, comments welcome

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

Voyage through the Hidden Physics of the Cosmic Web

Authors: A. Simionescu, S. Ettori, N. Werner, D. Nagai, F. Vazza, H. Akamatsu, C. Pinto, J. de Plaa, N. Wijers, D. Nelson, E. Pointecouteau, G. W. Pratt, D. Spiga, G. Vacanti, E. Lau, M. Rossetti, F. Gastaldello, V. Biffi, E. Bulbul, M. J. Collon, J. W. den Herder, D. Eckert, F. Fraternali, B. Mingo, G. Pareschi , et al. (5 additional authors not shown)

Abstract: The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons… ▽ More The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons. This will be the first mission capable to reach the accretion shocks located several times farther than the virial radii of galaxy clusters, and reveal the out-of-equilibrium parts of the intra-cluster medium which are live witnesses to the physics of cosmic accretion. It will also enable a view of the thermodynamics, kinematics, and chemical composition of the circumgalactic medium in galaxies with masses similar to the Milky Way, at the same level of detail that $Athena$ will unravel for the virialized regions of massive galaxy clusters, delivering a transformative understanding of the evolution of those galaxies in which most of the stars and metals in the Universe were formed. Finally, the proposed X-ray satellite will connect the dots of the large-scale structure by mapping, at high spectral resolution, as much as 100% of the diffuse gas hotter than $10^6$ K that fills the filaments of the cosmic web at low redshifts, down to an over-density of 1, both in emission and in absorption against the ubiquitous cosmic X-ray background, surveying at least 1600 square degrees over 5 years in orbit. This requires a large effective area (~10 m$^2$ at 1 keV) over a large field of view ($\sim1$ deg$^2$), a megapixel cryogenic microcalorimeter array providing integral field spectroscopy with a resolving power $E/ΔE$ = 2000 at 0.6 keV and a spatial resolution of 5 arcsec in the soft X-ray band, and a low and stable instrumental background ensuring high sensitivity to faint, extended emission. △ Less

Submitted 26 April, 2021; v1 submitted 5 August, 2019; originally announced August 2019.

Comments: White paper submitted in response to ESA's Voyage 2050 Call. Accepted for publication in Experimental Astronomy

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

Multi-scale analysis of turbulence evolution in the density stratified intracluster medium

Authors: Xun Shi, Daisuke Nagai, Erwin Lau

Abstract: The diffuse hot medium inside clusters of galaxies typically exhibits turbulent motions whose amplitude increases with radius, as revealed by cosmological hydrodynamical simulations. However, its physical origin remains unclear. It could either be due to an excess injection of turbulence at large radii, or faster turbulence dissipation at small radii. We investigate this by studying the time evolu… ▽ More The diffuse hot medium inside clusters of galaxies typically exhibits turbulent motions whose amplitude increases with radius, as revealed by cosmological hydrodynamical simulations. However, its physical origin remains unclear. It could either be due to an excess injection of turbulence at large radii, or faster turbulence dissipation at small radii. We investigate this by studying the time evolution of turbulence in the intracluster medium (ICM) after major mergers, using the Omega500 non-radiative hydrodynamical cosmological simulations. By applying a novel wavelet analysis to study the radial dependence of the ICM turbulence spectrum, we discover that faster turbulence dissipation in the inner high density regions leads to the increasing turbulence amplitude with radius. We also find that the ICM turbulence at all radii decays in two phases after a major merger: an early fast decay phase followed by a slow secular decay phase. The buoyancy effects resulting from the ICM density stratification becomes increasingly important during turbulence decay, as revealed by a decreasing turbulence Froude number $Fr \sim \mathcal{O}(1)$. Our results indicate that the stronger density stratification and smaller eddy turn-over time are the likely causes of the faster turbulence dissipation rate in the inner regions of the cluster. △ Less

Submitted 4 September, 2018; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: 8 pages, 7 figures, accepted to MNRAS

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

Cold Fronts and Shocks Formed by Gas Streams in Galaxy Clusters

Authors: Elad Zinger, Avishai Dekel, Yuval Birnboim, Daisuke Nagai, Erwin Lau, Andrey. V. Kravtsov

Abstract: Cold Fronts and shocks are hallmarks of the complex intra-cluster medium (ICM) in galaxy clusters. They are thought to occur due to gas motions within the ICM and are often attributed to galaxy mergers within the cluster. Using hydro-cosmological simulations of clusters of galaxies, we show that collisions of inflowing gas streams, seen to penetrate to the very centre of about half the clusters, o… ▽ More Cold Fronts and shocks are hallmarks of the complex intra-cluster medium (ICM) in galaxy clusters. They are thought to occur due to gas motions within the ICM and are often attributed to galaxy mergers within the cluster. Using hydro-cosmological simulations of clusters of galaxies, we show that collisions of inflowing gas streams, seen to penetrate to the very centre of about half the clusters, offer an additional mechanism for the formation of shocks and cold fronts in cluster cores. Unlike episodic merger events, a gas stream inflow persists over a period of several Gyrs and it could generate a particular pattern of multiple cold fronts and shocks. △ Less

Submitted 9 February, 2018; v1 submitted 17 September, 2016; originally announced September 2016.

Comments: Accepted to MNRAS, January 2018

The Galaxy Cluster Merger Catalog: An Online Repository of Mock Observations from Simulated Galaxy Cluster Mergers

Authors: J. A. ZuHone, K. Kowalik, E. Ohman, E. Lau, D. Nagai

Abstract: We present the "Galaxy Cluster Merger Catalog." This catalog provides an extensive suite of mock observations and related data for N-body and hydrodynamical simulations of galaxy cluster mergers and clusters from cosmological simulations. These mock observations consist of projections of a number of important observable quantities in several different wavebands as well as along different lines of… ▽ More We present the "Galaxy Cluster Merger Catalog." This catalog provides an extensive suite of mock observations and related data for N-body and hydrodynamical simulations of galaxy cluster mergers and clusters from cosmological simulations. These mock observations consist of projections of a number of important observable quantities in several different wavebands as well as along different lines of sight through each simulation domain. The web interface to the catalog consists of easily browseable images over epoch and projection direction, as well as download links for the raw data and a JS9 interface for interactive data exploration. The data is presented within a consistent format so that comparison between simulations is straightforward. All of the data products are provided in the standard FITS file format. Data is being stored on the yt Hub (http://hub.yt), which allows for remote access and analysis using a Jupyter notebook server. Future versions of the catalog will include simulations from a number of research groups and a variety of research topics related to the study of interactions of galaxy clusters with each other and with their member galaxies. The catalog is located at http://gcmc.hub.yt. △ Less

Submitted 22 December, 2017; v1 submitted 13 September, 2016; originally announced September 2016.

Comments: 12 pages, 6 figures, version accepted to ApJ and with additional authors; data at http://gcmc.hub.yt

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)

nIFTy galaxy cluster simulations I: dark matter & non-radiative models

Authors: Federico Sembolini, Gustavo Yepes, Frazer R. Pearce, Alexander Knebe, Scott T. Kay, Chris Power, Weiguang Cui, Alexander M. Beck, Stefano Borgani, Claudio Dalla Vecchia, Romeel Davé, Pascal Jahan Elahi, Sean February, Shuiyao Huang, Alex Hobbs, Neal Katz, Erwin Lau, Ian G. McCarthy, Giuseppe Murante, Daisuke Nagai, Kaylea Nelson, Richard D. A. Newton, Ewald Puchwein, Justin I. Read, Alexandro Saro , et al. (2 additional authors not shown)

Abstract: We have simulated the formation of a galaxy cluster in a $Λ$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanc… ▽ More We have simulated the formation of a galaxy cluster in a $Λ$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods. △ Less

Submitted 20 March, 2015; originally announced March 2015.

Comments: 21 pages, 13 figures, 4 tables - submitted to MNRAS

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