Arxiv:1810.00890V2 [Astro-Ph.CO] 4 Jan 2019 Bologna, Italy 2 Stephen Walker Et Al

Noname manuscript No. (will be inserted by the editor)

The Physics of Galaxy Cluster Outskirts

Stephen Walker · Aurora Simionescu · Daisuke Nagai · Nobuhiro Okabe · Dominique Eckert · Tony Mroczkowski · Hiroki Akamatsu · Stefano Ettori · Vittorio Ghirardini

Received: date / Accepted: date

Stephen Walker Astrophysics Science Division, X-ray Astrophysics Laboratory, Code 662, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA E-mail: [email protected] Aurora Simionescu SRON, Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands Institute of Space and Astronautical Science (ISAS), JAXA, 3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa, 252-5210, Japan Daisuke Nagai Department of Physics, Yale University, PO Box 208101, New Haven, CT, USA Yale Center for Astronomy and Astrophysics, PO Box 208101, New Haven, CT, USA Nobuhiro Okabe Department of Physical Science, Hiroshima University, 1-3-1 Kagamiyama, Higashi- Hiroshima, Hiroshima 739-8526, Japan Dominique Eckert Max-Planck-Institut fürextraterrestrische Physik, Giessenbachstrasse 1, 85748 Garching, Germany Tony Mroczkowski European Southern Observatory (ESO), Karl-Schwarzschild-Str. 2, D-85748 Garching, Ger- many Hiroki Akamatsu SRON, Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands Stefano Ettori INAF, Osservatorio di Astrofisica e Scienza dello Spazio, via Pietro Gobetti 93/3, 40129 Bologna, Italy INFN, Sezione di Bologna, viale Berti Pichat 6/2, I-40127 Bologna, Italy Vittorio Ghirardini Dipartimento di Fisica e Astronomia Universit di Bologna, Via Piero Gobetti, 93/2, 40129 Bologna, Italy INAF, Osservatorio di Astrofisica e Scienza dello Spazio, via Pietro Gobetti 93/3, 40129 arXiv:1810.00890v2 [astro-ph.CO] 4 Jan 2019 Bologna, Italy 2 Stephen Walker et al.

Abstract As the largest virialized structures in the universe, galaxy clusters continue to grow and accrete matter from the cosmic web. Due to the low gas density in the outskirts of clusters, measurements are very challenging, requiring extremely sensitive telescopes across the entire electromagnetic spectrum. Observations using X-rays, the Sunyaev-Zeldovich eﬀect, and weak lensing and galaxy distributions from the optical band, have over the last decade helped to unravel this exciting new frontier of cluster astrophysics, where the infall and virialization of matter takes place. Here, we review the current state of the art in our observational and theoretical understanding of cluster outskirts, and discuss future prospects for exploration using newly planned and proposed observatories. Keywords Galaxy clusters intracluster matter ·

1 What are the physical processes unique to cluster outskirts?

In recent years, the exploration of galaxy cluster outskirts has emerged as one of the new frontiers for studying cluster astrophysics and cosmology. It is a new territory for understanding the physics of intracluster medium (ICM) and intergalactic medium (IGM), and it is particularly important for the interpretation of recent X-ray and SZ observations and cluster mass estimates based on these observations. In the hierarchical structure formation model, galaxy clusters grow and evolve through a series of mergers and accretion from the surrounding large-scale structure in their outer envelope. Hydrodynamical cosmological simulations predict that accretion and mergers are ubiquitous and important for cluster formation. The accretion physics leaves distinctive marks especially in the outer envelopes of galaxy clusters, giving rise to caustics in the dark matter (DM) density proﬁle, an inhomogeneous gas density distribution, internal bulk and turbulent gas motions, and non-equilibrium electrons in the ICM. We describe all these eﬀects in detail in the subsections that follow.

1.1 Physically signiﬁcant scales and radii in the outskirts

Recent works have shown that accreting DM forms the so-called “splashback radius” (Diemer and Kravtsov 2014), which defines a natural boundary for a DM halo. The splashback radius (Rsp) arises from the recently accreted DM particles piling up at the apocenter of the halo as they re-emerge from the halo center after first infall (Adhikari et al. 2014; Shi 2016b). This process creates a caustic (see the dashed curves on the left panels of Fig.1) that can be identified as a sharp drop in the DM density profiles and is reflected also in the phase space structure of the DM (Mansfield et al. 2017; Diemer et al. 2017). Since the splashback process is governed by the infall velocity of DM particles and the gravitational potential of the halo, the splashback radius also depends on the mass accretion rate (MAR) of the DM halo (More et al. 2015), as shown in the right panel of Fig.1. The Physics of Galaxy Cluster Outskirts 3

3.0 Mansﬁeld+17 model ﬁt

Rsp 2.5 Rsh

m 2.0 200 R/R 1.5

1.0

0.5 0 2 4 6 Γ (0.0 z 0.5) 200m ≤ ≤

Fig. 1 Left panels: Gas entropy (top left), DM density (top right), gas temperature (bottom left) and gas pressure (bottom right) maps of the simulated cluster extracted from the non-radiative Omega500 hydrodynamical cosmological simulation. The inner dashed lines indicate the splashback shell computed using the method from Mansfield et al.(2017), whereas the outer white lines indicate the shock shell found by the discontinuous jump in entropy as well as pressure. Note that the accretion shocks are located outside of the splashback radius. Right panels: The splashback radius Rsp (red points) and the accretion shock radius Rsh (blue points), normalized by the halo radius R200m, plotted as a function of the mass accretion rate (MAR) of the cluster-size halos, Γ200m. The solid line is the Rsp Γ relation from DM-only cosmological simulations (Mansfield et al. 2017). The dashed− line indicates the fitted relation of Rsh/Rsp measured in the Omega500 simulation. Figures taken from Aung et al. in preparation and reprinted with permission.

The dynamics of accreting gas, on the other hand, is distinct from that of collisionless DM, because the collisional gas is shock heated during its first infall (Lau et al. 2015; Shi 2016a). The accretion shock (with radius Rsh) thus defines the spatial extent of the hot gas in DM halos (see the solid curves on the left panels in Fig.1) and plays an important role in a variety of cluster science, such as modeling of the ICM thermodynamics for X-ray and microwave surveys. This accretion shock is referred to as an “external shock” in the previous studies of shocks in galaxy clusters (e.g. Ryu et al. 2003; Molnar et al. 2009). The external accretion shock is due to the infall of low density gas from the void regions onto the cluster potential, and is different from “internal shocks” due to mergers and filamentary accretion. It is the external accretion shock that forms the outer boundary of the gas halo. Because the accretion shock directly relates to the infall process, it also serves as an important probe for mass assembly of the cluster gas. In practice, the mass of the cluster is usually quoted within an overdensity 3 radius r∆ as M∆ (4π/3)∆ρcrit(z)r , where ρcrit(z) is the critical density of ≡ ∆ the universe at the redshift of the cluster, and r∆ is the radius within which the mean enclosed density is ∆ ρcrit(z). ∆ = 500 roughly corresponds to the radius out to which the current× Chandra and XMM-Newton cluster observa- 4 Stephen Walker et al.

tions can reliably measure gas density and temperature profiles of the ICM. As such, any regions between R500c and the shock radius Rsh are considered the outskirts of galaxy clusters, where the ratios of various radius definitions are approximately given by R500c : R200c : R200m : Rsp : Rsh = 1 : 1.4 : 3 : 4 : 6. Perhaps the most common reference radius in cluster outskirts studies is R200c, often denoted more simply as r200 and referred to as the virial radius. The use of r200 for the virial radius stems from early calculations which considered the collapse of a spherical top hat perturbation. After initially expanding along with the rest of the universe, the gravitational force of such a perturbation is sufficient to slow its expansion, decoupling it from the universe’s expansion, and allowing it to eventually stop expanding (at its turn around radius) and collapse. Eventually the collapsing perturbation reaches virial equilibrium within its virial radius, which is equal to half of the turn around radius. For an Einstein-de Sitter universe (Ωm = 1) this virial radius is r 2 r178, which for simplicity has historically been rounded to r200. For 18π ≈ a ΛCDM universe (Ωm = 0.3, ΩΛ = 0.7), the virial radius found from this spherical top-hat calculation is actually r100=1.36r200. However, as the use of r200 to refer to the virial radius remains prevalent in the literature, it is this value that we shall use in this review. When scaling by the mean mass density of the universe rather than the critical density of the universe, the subscript letter changes from c to m. In this review, when it is not explicitly specified, the radii and masses are scaled by the critical density of the universe.

1.2 Radial Proﬁles of the Thermodynamical Properties in the ICM

15 −1 Figure2 shows the radial ICM profiles of a M200c = 1.1 10 h M cluster at z = 0 from the nIFTy radiative simulation cluster comparison× project (Sem- bolini et al. 2016). Starting from the same initial condition, these simulations were performed using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET) and modeling hydrodynamics with full subgrid physics of galaxy formation, including radiative cooling, star formation and thermal ac- tive galactic nucleus (AGN) feedback. This demonstrates that, while there are still significant variations in the simulated ICM profiles in the inner regions (r < R2500c), the prediction of the ICM profiles in the outskirts of clusters (r > R500c) is remarkably robust to the numerical schemes as well as the still uncertain cluster core physics. Figure3 shows the evolution of the ICM profiles in the Omega500 adiabatic cosmological simulation performed by using the ART code (Lau et al. 2015). This figure illustrates that the inner ICM profiles are self-similar when the profiles are normalized with respect to the critical density of the universe, while the ICM profiles in the outskirts are self-similar when they are normalized with respect to the mean mass density of the universe. This trend arises from the fact that the inner regions have formed early when the universe was still matter dominated (ΩM 1), while the cluster outskirts have accreted ≈ The Physics of Galaxy Cluster Outskirts 5

3 2 2 1 1 0

residuals 0 residuals -1 -1 -2 1016

1015 1000 ] ) ] 2 3 cm -2 Mpc -3 14 /(h 10 [keV h SUN Ramses-AGN Ramses-AGN 2/3 e M ArepoIL

-1 100 ArepoIL G3-X G3-X [ h G2-X G2-X G3-PESPH

GAS 13 G3-PESPH S = T/n

ρ 10 G3-Magnet.

G3-Magnet. G3-OWLS G3-OWLS ArepoSH ArepoSH G3-MUSIC G3-MUSIC G2-MUSICpi G2-MUSICpi 10 1012

0.1 1.0 0.1 1.0 R [h-1 Mpc] R [h-1 Mpc]

1.0 3 0.5 2 0.0 1 -0.5

residuals -1.0 residuals 0 -1.5 -1 Ramses-AGN Ramses-AGN 15 ArepoIL 1.5 ArepoIL G3-X G3-X G2-X G2-X G3-PESPH 15 G3-PESPH G3-Magnet. G3-Magnet. G3-OWLS )] ArepoSH m G3-OWLS Ω ArepoSH G3-MUSIC / b G3-MUSIC

G2-MUSICpi Ω 10 1.0 G2-MUSICpi K]

10 7 (< R)/( gas T [10 kT [ keV ]

5 (< R) [f 0.5 gas ϒ

0 0 0.0 0.1 1.0 0.1 1.0 R [h-1 Mpc] R [h-1 Mpc]

Fig. 2 Radial ICM proﬁles at z = 0 from the nIFTy radiative simulation cluster comparison project, showing gas density (top-left), entropy (top-right), temperature (bottom-left), and gas fraction (bottom-right). The top panel shows the diﬀerence between each simulation and the reference G3-MUSIC simulation. The dashed line indicates to R2500c and the dotted line to R500c for the reference G3-MUSIC values. Figures taken from Sembolini et al.(2016), reprinted with permission.

materials recently and their mass accretion rate is aﬀected by the repulsive force of dark energy.

1.3 Gas clumping

Simulations of clusters (Nagai and Lau 2011; Roncarelli et al. 2013) predict that in the outskirts the gas density distribution becomes increasing inhomoge- nous (i.e. ‘clumpy’). One might expect the survivability of these gas clumps to thermal conduction and hydrodynamic instabilities during the infall to depend on details of the ICM microphysics – therefore, in addition to being a source of bias that must be addressed and quantiﬁed, the level of gas clumping 6 Stephen Walker et al.

z = 0.0 10 10 z = 0.5 z = 1.0 c m 200

z = 1.5 200 /K /K ) )

r 1 1 r ( ( K K

0.1 0.1 0.4 0.4 0.2 0.2 0.0 0.0 0.2 0.2 −0.4 −0.4 −0.6 −0.6 − 0.1 1 − 0.1 1 r/R200c r/R200m

Fig. 3 Evolution of the average gas entropy profiles at z = 0, 0.5, 1.0, and 1.5 from the Omega500 adiabatic cosmological simulation. The lower subplot in each panel shows the deviations of the profiles from different z outputs relative to that of the z = 0.5 clusters. The left panel shows the profiles of cluster halos defined using the critical density, while the right panel shows the profiles of the same cluster halos defined using the mean density. The shaded regions indicate the 1σ scatter around the mean gas profiles for the z = 0 clusters. Figures taken from Lau et al.(2015), reprinted with permission.

can ultimately hold important clues about the physics of large-scale structure growth. In X-ray analyses of galaxy clusters, it is commonly assumed that the ICM is a single phase medium characterized by one temperature and gas density within each radial bin. If this approximation holds, then the thermodynamical properties of the ICM inferred from X-ray studies should follow closely the radial trends predicted by the numerical simulations shown in the previous section. However, if the ICM is clumpy, the gas density inferred from the X- p ray surface brightness is overestimated by C(r), with

2 2 C ρ / ρgas . (1) ≡ h gasi h i

Although the bias in the inferred ICM mass is moderate (< 10%) in the high- pressured regions in the interior of galaxy clusters (Mathiesen et al. 1999), gas clumping can become significant in the outskirts (r > r200) and could serve as a major source of systematic bias in X-ray measurements of ICM profiles (Nagai and Lau 2011; Roncarelli et al. 2013). Fig.4 shows that the clumping factor of the X-ray emitting gas (T > 106 K) expected from numerical simulations is C 1.3 at r = r200, and it increases with radius, reaching C 5 at r = 2r200. ∼ 2/3 ∼ From the definition of entropy S Te/ne , the overestimate of gas density due to clumping causes an underestimate≡ of the observed entropy profile by C(r)1/3 if the clumps have the same temperature as the ambient medium. The Physics of Galaxy Cluster Outskirts 7

Fig. 4 Left panel: the mass dependence of the median clumping factor proﬁles of gas with 6 14 1 T > 10 K. Results are shown for the present-day (z = 0) clusters with M200 > 10 h− M 14 1 (solid lines) and M200 < 10 h− M (dashed lines), respectively, in the runs including cooling and star formation (CSF, thick) and non-radiative (NR, thin) runs. Right panel: the redshift dependence of the median clumping factor proﬁles in the CSF (thick) and NR (thin) runs at z = 0 (solid lines) and z = 0.6 (dashed lines). The errorbars on the solid lines indicate the interquartile range. Figures taken from Nagai and Lau(2011), reprinted with permission.

1.4 Bulk motions, turbulence and deviations from hydrostatic equilibrium

Merging and accretion events generate a significant level of bulk and turbulent gas motions within the ICM and affect its thermodynamic structure (see the left panel of Figure5). Hydrodynamical simulations predict that the non- thermal pressure provided by these gas motions increases from 10% at ≈ r = R500c to 30% at r = 2R500c (Lau et al. 2009; Vazza et al. 2009; Battaglia et al. 2012) on≈ average. The non-thermal pressure fraction also depends on the mass accretion rate (Nelson et al. 2014; Shi et al. 2015) which gives rise to a scatter for cluster-size halos with any given mass (see the right panel of Fig. 5). Understanding the non-thermal pressure in galaxy clusters is especially important because X-ray and SZ observations typically measure only the thermal pressure of the gas. The non-thermal pressure, if neglected, introduces biases in our physical understanding of the ICM profiles as well as the hydrostatic mass estimation. The increasing kinetic energy fraction at larger radii in particular can lead to suppressed temperatures and flatter entropy profiles in cluster outskirts. Gas motions also contribute to the support against gravity which leads to biases in cluster mass estimates based on the assumption of hydrostatic equilibrium (HSE) (Rasia et al. 2006; Nagai et al. 2007; Nelson et al. 2012, 2014a; Shi et al. 2016). The HSE mass bias remains one of the primary sources of systematic uncertainties in the calibration of cluster observable-mass rela- 8 Stephen Walker et al.

1e7 1.5 Relaxed 8

1.0 7

0.5 6

5 0.0

Mpc/h 4 Temperature (K) Temperature −0.5 3

−1.0 2

Velocity 1 −1.5 −1.5 −1.0 −0.5 0.0 0.5 1.0 1.5 Mpc/h 500 km/s

Fig. 5 Left panel: Gas velocity field (shown by the white vectors) in a slice through a galaxy cluster simulation after subtracting the global mean bulk velocity of the cluster. The underlying colour map represents the gas temperature of this relaxed cluster. The black circles show r2500, r500 and r200. Figure taken from Nelson et al.(2014a), reprinted with permission. Right panel: The profile of non-thermal pressure fraction, f P /P nth ≡ nth total (where the total pressure, Ptotal, is the sum of the thermal pressure, Pth, and the non- thermal pressure, Pnth), as a function of the cluster-centric radius scaled with r200m at z = 0. Modelled and simulated profiles of an early growth sub-sample (Γ200m < 1.8, red lines) and a late growth sub-sample (Γ200m > 2.7,lines) are shown. This plot compares the predicted profiles from the analytic model of Shi and Komatsu(2014) (modeled) with the profiles for a cosmological hydrodynamics simulation from Nelson et al.(2014) (simulated), showing good agreement between the two methods. The lines and the shaded regions are the mean profile of the sample and the 16/84 percentiles, respectively. Figure taken from Shi et al.(2015), reprinted with permission.

tions and cosmological parameters derived from clusters (Allen et al. 2008; Vikhlinin et al. 2009; Planck Collaboration et al. 2013a).

1.5 Electron-Ion Equilibration in Cluster Outskirts

X-ray and SZ measurements of the hot gas in galaxy clusters are often interpreted under a number of simplifying assumptions on the poorly understood physics of the ICM. One such assumption is that the electrons are in thermal equilibrium with the surrounding ions; i.e., Te = Ti. However, this assumption is untenable in the low-density outskirts of galaxy clusters (see Fig.6) due to the extended timescale for electrons to reach equilibrium via Coulomb collisions (Avestruz et al. 2015, and references therein). When an electron-ion plasma passes through a shock, most of the kinetic energy goes into heating the heavier ions, causing Ti Te. After the shock, electrons and ions slowly equilibrate via Coulomb interactions, each converging to the mean gas temperature, Tgas = (neTe + niTi)/(ne + ni), over a typical electron-ion equilibration timescale, tei. The evolution of the electron temperature is given by,

dTe Ti Te = − (γ 1)Te ( v) , (2) dt tei − − ∇ · The Physics of Galaxy Cluster Outskirts 9

1.0 1.00

0.9 0.95 0.8 0.90 0.7 s a 0.6 g 0.85 T / 0.5 e T 0.80 0.4 0.75 14 1 0.7

Fig. 6 Left panels: Sliced maps of Te/Tgas for the most relaxed and least massive cluster from the Omega500 non-radiative simulation. The colorbar indicates Te/Tgas. The dimension for each panel is 15.6 h 1Mpc 15.6 h 1Mpc, with depth of 7.6 h 1kpc. The circle in dashed − × − − line shows R200m of the cluster. Between 1.0 < R/R200m < 1.5, Te/Tgas is close to unity in the filamentary gas with high momentum flux entering the cluster well, but is significantly less than unity in the diffuse gas at the same cluster-centric radii. Right panels: Black lines show the average profiles of the mass-weighted temperature bias binned by mass for a sample of 65 simulated clusters at z = 0 and their mainline progenitor data at z = 0.5, 1, and 1.5 prior to binning the profiles by mass. The number of clusters in each mass bin is denoted in parentheses in the legend label. Red lines show the fitting function (equation4) with input masses that correspond to the average mass in each bin. The shaded regions show the scatter of the cluster in each mass bin around the average profiles. Figures taken from Avestruz et al.(2015), reprinted with permission. where the second term accounts for adiabatic compression heating and cooling. For the fully ionized ICM, including contributions from both protons and He++, the timescale for equilibration is given by

7 3/2 8 Te/10 K tei 6.3 10 yr −5 −3 (3) ≈ × (ni/10 cm ) (ln Λ/40)

(Spitzer 1962). In the outskirts of galaxy clusters, the collision rate of electrons and protons becomes longer than the age of the universe. The resulting deviation between electron and ion temperatures is larger in more massive and less relaxed systems, ranging from 5% in relaxed clusters to 30% for clusters undergoing major mergers (Fig.6, right panel). Avestruz et al.(2015) have found a simple analytic ﬁtting function for the proﬁle of the temperature bias (shown as the red lines in the right hand panel of Fig.6):

−a Te(R) (x/xt) be(R) = , (4) −b a/b ≡ Tgas(R) (1 + (x/xt) ) 10 Stephen Walker et al.

where x = R/R200m and the parameters which account for the mass dependence are,

14 −1 [−0.1798,0.3448,−0.006] [xt, a, b] = [0.629, 0.086, 2.851] M200m/10 h M × These results are for non-radiative runs of the Omega500 simulations. Avestruz et al.(2016) have found that performing the same simulations with cooling, star formation and feedback included has little effect on the shape of the ICM temperature profiles. The presence of non-equilibrium electrons leads to a significant suppression of the SZ effect signal at large cluster-centric radius. The suppression of the electron pressure also leads to an underestimate of the hydrostatic mass. Merger-driven, internal shocks may also generate significant populations of non-equilibrium electrons in the cluster core, leading to a 5% bias on the integrated SZ mass proxy during cluster mergers (Rudd and Nagai 2009). Note, however, that the Spitzer timescale adopted in these simulations is an upper limit of the true equilibration time, as there are other physical mech- anisms such as plasma instabilities that can couple the temperatures of the different ion species. For example, since the gyroradii are considerably smaller than the mean free path, it is possible that the electrons and ions equilibrate near-instantaneously via interactions mediated by the magnetic field. Future studies of the difference between the electron and ion temperatures will shed new insights into the microphysics of how the gas is heated in galaxy clusters.

2 What tools do we have for investigating cluster outskirts physics?

2.1 Weak lensing and galaxy distribution

The weak gravitational lensing (WL) effect leads to small distortions of the apparent shapes of background galaxies located behind a massive object (e.g. Bartelmann and Schneider 2001). WL analysis, which utilizes these coherent distortions in a statistical way, is a powerful technique to reconstruct the mass distribution from the central cluster region all the way to the cluster outskirts, without invoking assumptions about the dynamical state. Notably, the advent of wide-field cameras installed on large ground-based telescopes (e.g., Subaru/Suprime-cam, Subaru/HSC, and Canada-France-Hawaii Tele- scope) enables us to measure the cluster mass distribution beyond R500c. Since cluster outskirts are significantly affected by mass accretion flows, WL analysis provides us with an important direct probe of the mass distribution of galaxy clusters and their surrounding environments. Recently, Sereno et al.(2018) showed that the lensing signal can be detected out to a remarkable distance of 30 Mpc around a massive SZ-selected cluster halo, with a measured shear that∼ is very large and marginally consistent with ΛCDM predictions. Along with a tremendous improvement of the observing instruments, wide field surveys, such as Hyper Suprime-Cam Subaru Strategic Program (HSC- SSP; Aihara et al. 2018), Canada-France-Hawaii Telescope Legacy Survey The Physics of Galaxy Cluster Outskirts 11

(CFHTLS; Shan et al. 2012), the Dark Energy Survey (DES; Dark Energy Sur- vey Collaboration et al. 2016), and Sloan Digital Sky Survey (SDSS; Eisenstein et al. 2011) can and will compile a considerable number of galaxy clusters and measure an average mass profile stacked over the sample, with little statistical uncertainty. Due to their relatively collisionless nature, member galaxies are also expected to trace dark matter (e.g. Okabe and Umetsu 2008), providing an excellent proxy for the mass distribution even when WL measurements are difficult. The complementary techniques of galaxy and WL mass distributions provide us with a unique opportunity to map cluster outskirts, especially in a < < radial range of R500c r R200m. Since it is challenging∼ ∼ to detect diffuse gas beyond r200 through X-ray and SZ measurements, WL and galaxy distributions are a compelling route for understanding the physics of cluster outskirts around R200m, roughly corresponding to the predictions of the splashback radius (e.g. Oguri and Hamana 2011; Diemer and Kravtsov 2014; Adhikari et al. 2014; More et al. 2015).

2.2 X-ray Observations

The hot, tenuous ICM is generally well described as a collisionally ionized, optically thin plasma. Particularly in the low-density regions located at the cluster outskirts, optical depth effects are expected to be negligible, enabling a straightforward interpretation of X-ray cluster spectra. Three fundamental processes contribute to this X-ray emission. At the high temperatures (>2 keV) typical for the plasma found in massive galaxy clusters, the dominant of these is thermal bremsstrahlung (free-free emission); additional contributions to the spectrum arise from recombination radiation caused by the capture of an electron by an ion following ionization (free-bound emission), and from deexcitation radiation of an electron changing the quantum level in an ion (bound-bound emission). Both the shape of the bremsstrahlung continuum and ratios between the strengths of various emission lines can be used as a diagnostic for the plasma temperature. It is important to note that all the radiative processes listed above depend on the encounter between an electron and an ion, implying that the X-ray emissivity is always proportional to the square of the plasma density. The strength of the X-ray signal therefore de- creases rapidly as the density drops from 10−1 cm−3 in the cluster cores to 10−5 cm−3 near the virial radii. ∼ ∼ Hence, a detailed understanding of the instrumental and cosmic backgrounds, and their associated systematic uncertainties, is of utmost impor- tance for enabling robust measurements of the properties of the X-ray faint cluster outskirts. X-ray satellites located in a low-Earth orbit (LEO) are at an advantage for these studies, because the Earth’s magnetic field shields the detectors onboard, resulting in a reduced and more stable particle background. Among the X-ray telescopes in LEO, the work horses of cluster outskirts studies to date have been the Röntgensatellit (ROSAT), launched in 1990 and led 12 Stephen Walker et al.

by the German Aerospace Center, and Suzaku, launched in 2005 by the Japan Aerospace Exploration Agency. While the limited energy band and spectral resolution of ROSAT only allow robust measurements of the emissivity and gas density (e.g. Neumann and Arnaud 1999; Ettori and Fabian 1999), Suzaku has provided the first X-ray spectroscopic constraints on the gas temperature and metallicity near the virial radii of galaxy clusters (Reiprich et al. 2009; George et al. 2009; Bautz et al. 2009; Fujita et al. 2008; Simionescu et al. 2011). In terms of the cosmic foregrounds and backgrounds, X-ray spectroscopy with Suzaku has largely benefited from complementary data sets obtained from the ROSAT All-Sky Survey (RASS), which played a crucial role in constraining the soft X-ray foreground from the Galactic halo, and from Chandra, whose exquisite spatial resolution allowed for the identification and removal of distant point sources contaminating the cluster outskirts signal. The exciting physics near galaxy clusters’ virial radii initially hinted at by early Suzaku observations sparked a new effort to improve the background modeling of other telescopes not located in LEO; despite their highly elliptical orbit and more variable particle background, Chandra and XMM-Newton also succeeded in placing spectroscopic constraints on the properties of the ICM in the outskirts of some systems (Urban et al. 2011a; Bonamente et al. 2013a). In addition to spectroscopy, deep, high-resolution imaging with Chandra has provided complementary information on the structure of the material being accreted onto clusters of galaxies. Likewise, the magnitude of surface brightness fluctuations in ROSAT and XMM-Newton images have provided constraints on the surviving clumps of gas falling into the ICM from the surrounding large-scale structure.

2.3 Cluster Outskirts via the Thermal Sunyaev–Zeldovich Eﬀect

The thermal Sunyaev–Zeldovich effect (Sunyaev and Zeldovich 1972, SZ effect) has long been regarded as a promising tool for probing the cluster outskirts. This is due to its linear dependence on thermal electron pressure, as compared to the X-ray which diminishes as density squared; the surface brightness of the thermal SZ effect is proportional to the Compton y param- 2 R eter, where y = σT /mec Pe(`)d` (see e.g. Carlstrom et al.(2002)). Stacked and ensemble-averaged results using e.g. the Atacama Cosmology Telescope (ACT), Bolocam, the Planck satellite probe of the cosmic microwave background (CMB), and the South Pole Telescope (SPT) have begun to deliver on this promise by placing competitive constraints on the outer slopes of cluster pressure profiles (e.g. Planck Collaboration et al. 2013b; Sayers et al. 2013, 2016; Romero et al. 2017, 2018), though much work remains for extending this to individual systems, lower masses, and higher redshifts. Assuming the Arnaud et al.(2010) pressure profile, one can predict the levels of Compton y 10−7– 10−6 one would need to directly image the outskirts of a typical massive∼ galaxy cluster (see Figure7). As a comparison, the Planck SZ measurement of the massive cluster Abell 2319 (Ghirardini et al. The Physics of Galaxy Cluster Outskirts 13

2018a) is shown in the right hand panel of Figure7. Also for comparison, the stacked Planck SZ measurement of the filamentary material between galaxies was found to be at a level of y 10−8 (Tanimura et al. 2017; de Graaff et al. 2017), consistent with the Arnaud∼ et al.(2010) pressure profile prediction at 6 R500 in Figure7. × Direct imaging of the cluster outskirts through the SZ effect will be more challenging than stacked or radially-averaged measurements, and not simply because longer integration times are required. The 1-σ confusion limit due to the cosmic infrared background (CIB), composed of high-z dusty galaxies that are bright at submillimeter and millimeter wavelengths, is expected to dominate the RMS in single-band SZ measurements at levels of y 10−6. For typical beam sizes of a few arcseconds to ten arcminutes, the confusion∼ is nearly independent of beam size (see e.g. Bétherminet al.(2017) and the discussion of confusion in the SZ review in these proceedings). However, due to the spatial correlation of sources, smaller beams are strongly favored as the higher resolution can more easily separate the contributions of spatially- correlated sources. Higher resolution instruments also provide the opportunity to average over many beams to reach lower RMS noise levels than a confusion- limited measurement with a broad beam, while multi-chroic/spectroscopic instruments can directly fit for the CMB and CIB contributions. Thus the trend toward higher spatial and spectral resolutions holds great potential for direct SZ imaging of cluster outskirts. The works cited above as well as those discussed in Section 2.4.1 have shown that stacked and radially-averaged SZ measurements in the cluster outskirts can be used in combination with lensing or X-ray constraints to infer the level of non-thermal pressure support (Siegel et al. 2018), the level to which thermal HSE is valid (Ettori et al. 2018), and may tentatively indicate a detection of the long-sought accretion shock (Hurier et al. 2017). As SZ instrumentation progresses (see Section 4.1.4), SZ measurements will improve constraints on ICM turbulence and non-thermal pressure, allow tests of the nature of clumping in the outskirts, and provide better constraints on the accretion shock jump conditions. For a review of the many astrophysical processes SZ measurements can probe, see Mroczkowski et al.(2018) in these proceedings.

2.4 Multi-wavelength synergies

2.4.1 X-ray+SZ

The advent of the new generation of SZ instruments has opened the exciting prospect of combining the information obtained from X-ray and SZ techniques. Indeed, when the density nX can be reconstructed from the soft X-ray surface brightness and the thermal pressure PSZ can be directly measured from the SZ eﬀect, the two quantities can be directly combined through the ideal gas equation to recover the gas temperature kT = PSZ /nX and the entropy K = −5/3 PSZ nX . The X-ray/SZ method has the advantage of bypassing the use 14 Stephen Walker et al.

10-5

10-6

10-7

10-8

101 102 103 Radius (kpc)

Fig. 7 Left panel: Model Compton y profile of a 3 1014M cluster at z = 0.5 assuming an × Arnaud et al.(2010) pressure profile. The black dashed line marks R500, and the red dashed line marks 6 R500. Right panel: Planck measured Compton y profile for Abell 2319, from Ghirardini et× al. 2018b, reprinted with permission. The dotted and dashed vertical lines show the location of R500 and 2 R500 respectively. Note that this cluster is more massive than the toy model shown in the× left panel, which explains the higher Compton y value at R500.

of spectroscopic X-ray information, which requires high signal-to-noise data and is prone to potential systematics in the background-dominated regime of cluster outskirts. Conversely, the joint X-ray/SZ technique only requires a detection of the X-ray signal in a soft energy band where the sensitivity of current X-ray instruments peaks, thus allowing our measurements to extend to larger radii. The combination of X-ray and SZ data thus represents a promising avenue towards accurate measurements of the cluster properties out to the virial radius. The technique of inferring temperature from X-ray+SZ surface brightness was pioneered two decades ago (as discussed in Mroczkowski et al. (2018)), and was applied in Eckert et al.(2013a,b) for a sample of 18 clusters with available ROSAT/PSPC and Planck SZ data, extending the methodology into the cluster outskirts.

The XMM-Newton Cluster Outskirts Project (X-COP, Eckert et al. 2017) applies a similar technique on a sample of nearby clusters with high signal- to-noise ratio in the Planck survey. The project provides deep XMM-Newton mapping for 13 massive clusters covering almost uniformly the region extending out to 2R500c R100c. After the pilot studies targeting A2142 (Tchernin et al. 2016) and A2319≈ (Ghirardini et al. 2018b) that demonstrated the potential of XMM-Newton and Planck to reconstruct jointly the thermodynamic properties with high precision all the way out to R100c, Ghirardini et al.(2018a) presented the final results on the profiles of the most relevant thermodynamic quantities (gas density, temperature, pressure and entropy) recovered over two decades in radius, assessing the level of their “universality” (once rescaled only for the mass within a given overdensity) and the amount of their intrinsic scatter. The combination of these X-ray/SZ based profiles allowed also to constrain the hydrostatic mass with a relative (statistical) median error of 3 % at R500c and 6% at r200 (Ettori et al. 2018). The Physics of Galaxy Cluster Outskirts 15

2.4.2 X-ray+WL

As discussed in Section 1.4, the fraction of non-thermal pressure due to gas motions compared to the thermal pressure is expected to increase significantly towards the cluster outskirts, introducing uncertainties in the mass determi- nation when using X-ray data alone and relying on the assumption of HSE. In addition, the possible presence of gas clumping can affect the X-ray gas density and consequently also lead to biases in the mass and overdensity radius. On the other hand, gravitational lensing studies are complementary to X-ray measurements, because lensing observables do not require any assumptions on the cluster dynamical state or ICM microphysics. WL mass estimates are, however, sensitive to assumptions about the 3D shapes and halo orientations (e.g. Oguri et al. 2005) and substructures (e.g. Okabe et al. 2014a) in the cluster gravitational potential, as well as any other large-scale structure between the lensed sources and the observer (e.g. Hoekstra 2003). A comparison of X-ray and WL masses for a large number of clusters is a powerful route to test the validity of the HSE assumption in the cluster outskirts. The two independent measurements indirectly constrain non-thermal pressure components involved in turbulence and/or bulk motions and their radial dependence. As discussed in Section 1.2, the ICM temperature and number density profiles are expected to follow ‘universal’ profiles when normalized by the mass parameters, as follows:

M E(z) a f (M , r/r ) = n E(z)2 ∆ g (r/r ) = n g (r/r ) (5) n ∆ ∆ 0 14 −1 n ∆ ∗ n ∆ 10 h70 M M E(z) b f (M , r/r ) = T ∆ g (r/r ) = T g (r/r ). T ∆ ∆ 0 14 −1 T ∆ ∗ T ∆ 10 h70 M

Here, M∆ is the enclosed mass within r∆ and a and b are slopes of the mass dependence (a = 0 and b = 2/3 for a self-similar solution); gn and gT are the radial functions for the temperature and density scaled by r∆, respectively. Given the WL masses and X-ray observables, one can simultaneously fit the gas density and temperature profiles with WL masses, self-consistently taking into account their uncertainties. The joint likelihood for the number density and the temperature profiles is given by

X T −1 2 ln = ln(det(Cij)) + v C vij, (6) − L ij ij i,j ln(n(˜r)) ln(fn(M∆, r˜)) v = − , ln(T (˜r)) ln(fT (M∆, r˜)) − C = Cstat + Cint

Here, i and j denote the i-th cluster and j-th radial bins, respectively. Cstat is the error covariance matrix for the data vector v, Cint is the intrinsic covariance andr ˜ = r/r∆. When one considers a single observable, the oﬀ-diagonal 16 Stephen Walker et al.

101 Mi-5 logh<-19.43 Mi-5 logh<-20.43 Mi-5 logh<-21.43

] 1

2 10 − 101

Mpc 0 2 10 h )[

R 0 ( 10 g

Σ 0 1 1 1 10 Vpeak > 135kms− Vpeak > 175kms− Vpeak > 280kms− 1 Observed Observed 10− Observed

1 0 1 1 0 1 1 0 1 10− 10 10 10− 10 10 10− 10 10

0.2 − 0.4 − 0.6 − R 0.8 − log 1.0

/d − g 1.2 − 1.4 log Σ − d 1.6 − 200m 200m 200m

1.8 R R R − 2.0 1 0 1 1 0 1 1 0 1 − 10− 10 10 10− 10 10 10− 10 10 1 1 1 R [h− Mpc] R [h− Mpc] R [h− Mpc]

Fig. 8 Top panels: the surface number density profiles of SDSS photometric galaxies with different magnitude thresholds around the entire redMaPPer cluster sample with z [0.1, 0.33] and richness λ [20, 100] , are shown using symbols with error bars. The dashed∈ lines correspond to (sub)halo∈ surface density profiles in the MultiDark-Planck II simulation around clusters. Vpeak is the maximum circular velocity of the dark matter subhalos for the galaxies corresponding to the different magnitude thresholds used. Bottom panels: the logarithmic slope of the surface density profiles are shown using solid and dashed lines for the observed galaxy and the subhalo surface density profiles, respectively. The observed slope of the surface density profile has a shape which is similar to that expected from simulations. Figure taken from More et al.(2016), reprinted with permission. elements in C are set to be zero. The pressure and entropy profiles can be computed by a combination of these functions. A tremendous advantage of this method is that it does not need to assume or fix the radial slopes of X- ray observables based on theoretical predictions. The best-fit slopes directly provide us with the radial changes of the investigated X-ray quantities (Okabe et al. 2014b).

3 What are the latest constraints on the physical processes that shape the cluster outskirts?

3.1 Splashback radius

Recent studies using galaxy and WL mass distributions have provided a deeper understanding of the splashback radius that is a physical boundary of galaxy clusters (as explained in Section 1.1). More et al.(2016) discovered significant evidence for the splashback radius using the projected number density profiles The Physics of Galaxy Cluster Outskirts 17 of photometric galaxies (Fig.8). They computed the stacked surface number densities around each cluster from the redMaPPer cluster catalog (Rykoff et al. 2014; Rozo and Rykoff 2014) obtained from the Sloan Digital Sky Survey, as a function of comoving projected distance from the cluster center. The surface densities show a clear, sharp steepening around scales of 1h−1 Mpc ∼ 0.9r200m (Figure8), in good agreement with predictions of the splashback radius(e.g.∼ Diemer and Kravtsov 2014; Adhikari et al. 2014; More et al. 2015) in three-dimensional space. More et al.(2016) also subdivided their galaxy cluster sample into two subsamples, namely high and low galaxy concentrations, and found a 6.6σ difference in the clustering amplitudes of the two surface density profiles. Zu et al. (2017) and Busch and White(2017) have pointed out that such an apparent assembly bias would be caused by the classification based on projected posi- tions of cluster members, and that it does not appear in the three-dimensional distribution of matter around clusters. Additional studies using photometric galaxies (galaxies which have only photometric redshifts) have been carried out (e.g. Baxter et al. 2017; Nishizawa et al. 2018). Baxter et al.(2017) found using SDSS data that the logarithmic slope of halo profile of the cluster (which they separate from the surrounding infalling material) reaches about 5 over a narrow range of radius, and galaxy colors become abruptly redder within− the location of the steepening of the density profile, suggesting at the quenching of star formation. Nishizawa et al. (2018) investigated whether the splashback radius could be seen using the HSC-SSP data, however the relatively small survey area of 200 deg2 makes it difficult to detect because of statistical uncertainty. They∼ found that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirts of clusters. Umetsu and Diemer(2017) have carried out a stacked weak-lensing analysis for 16 clusters from the Cluster Lensing And Supernova survey with Hubble (CLASH) program. Since the data area is limited within the field-of-view of the Subaru/Suprime-Cam ( 0.25 deg2), the study provided a lower limit on ∼ 3D the splashback radius of the clusters at the median redshift z 0.35, of Rsp > 3D −1 ' 0.89r200m, corresponding to Rsp > 1.83 h Mpc. Chang et al.(2017) measured the galaxy number density and weak-lensing mass profiles around redMaPPer galaxy clusters in the first year Dark En- ergy Survey (DES) data. They found strong evidence of a splashback-like −1 steepening of the galaxy density profile, Rsp = 1.16 0.08 h Mpc, and the first detection of a splashback-like steepening in± weak lensing signals, −1 Rsp = 1.28 0.18 h Mpc. Chang et al.(2017) show that a change of slope in the weak-lensing± signal around the splashback-like radius is weaker than that in the projected galaxy distribution. This is because lensing signals correspond to Σ¯(< R) Σ(R), where Σ(R) is the projected mass density at the projected distance R−and Σ¯(< R) is the average of the projected mass density within R. They found that the two independent measurements of Rsp are in good agreement. 18 Stephen Walker et al.

A2199 PKS0745-191 A2204 A1914 A2744 A2667 A3376 RXCJ0605 A3667 Antlia Perseus Centaurus A399&401 HydraA A1644 A1246 A478 A2255 A1413 A1689 A85 A2319 ESO3060170 CIZA 2242 A2029 A3158 A3266 A773 Coma A2142 A1068 Virgo A644 RXJ1159 Subaru RXC1825 ZW1215 A383 A1835 Zw2341 A1795 UGC03957 Chandra

A133 XMM+Planck Suzaku A1750

Fig. 9 Venn diagram showing the coverage of cluster outskirts by Suzaku, Chandra, Sub- aru weak lensing, and joint XMM/Planck observations. For each cluster measurements of thermodynamic proﬁles (density, temperature, entropy) are available out to at least r200. The individual works behind each cluster’s analysis are tabulated in table1. ≈

3.2 Thermal properties of the gas

The recent advances in our ability to measure the thermodynamic properties of the ICM in the outskirts of galaxy clusters have come through two different avenues. With Suzaku (2005-2015), X-ray astronomy for the first time had a low and stable background, high effective area, CCD spectral resolution X- ray satellite capable of accurate spectroscopy at and beyond r200. This allowed direct temperature and thus entropy measurements of the ICM near the cluster virial radii to be achieved through X-ray spectral fitting. With Planck (2009- 2013), highly sensitive SZ observations of a large sample of clusters has been compiled, allowing the ICM pressure to be measured out to r200. By combining these SZ pressure profiles with X-ray surface brightness data from ROSAT and XMM, independent measurements of the ICM properties in the cluster outskirts have been achieved. While early progress concentrated on exploring the outskirts of individual systems, over the last 5 years these studies have expanded to samples of clusters and lower mass groups, allowing general trends in the cluster population to be studied. The current samples of clusters observed are summarized in the Venn diagram in Fig.9, which shows the coverage available from different observatories. For each cluster in this diagram, thermodynamic profiles (temperature, density, entropy and pressure) can be extracted out to at least The Physics of Galaxy Cluster Outskirts 19

Table 1 List of references for studies of the outskirts of individual galaxy clusters shown in Fig.9.

Object name References Object name References Perseus Simionescu et al.(2011) Abell 1835 Okabe et al.(2014b) Simionescu et al.(2012) Ichikawa et al.(2013), Urban et al.(2014) Bonamente et al.(2013b) Coma Simionescu et al.(2013) Morandi et al.(2012) Akamatsu et al.(2013a) Abell 2029 Walker et al.(2012c) Virgo Urban et al.(2011b) Abell 2142 Akamatsu et al.(2011) Simionescu et al.(2017) Tchernin et al.(2016) Hydra A Okabe et al.(2014b) Abell 2204 Reiprich et al.(2009) Sato et al.(2012) Abell 2255 Ghirardini et al.(2018b) Centaurus Walker et al.(2013) Abell 2319 Ghirardini et al.(2017) PKS0745 Walker et al.(2012a) Abell 2744 Ibaraki et al.(2014) Sculptor Sato et al.(2010) Abell 3158 Ghirardini et al.(2018b) Antlia Wong et al.(2016) Abell 3266 Ghirardini et al.(2018b) Abell 85 Ghirardini et al.(2018b) Abell 3376 Akamatsu et al.(2012b) Ichinohe et al.(2015) Abell 3667 Akamatsu et al.(2012a) Abell 133 Morandi and Cui(2014) CIZA2242.8 Akamatsu et al.(2013b) A399/401 Akamatsu et al.(2017) Jee et al.(2015) Abell 478 Okabe et al.(2014b) Okabe et al.(2015) Abell 644 Ghirardini et al.(2018b) ZwCl2341.1 Akamatsu et al.(2013b) Abell 1246 Sato et al.(2014) ZW1215 Ghirardini et al.(2018b) Abell 1413 Hoshino et al.(2010) RXJ1159+5531 Su et al.(2015) Abell 1644 Ghirardini et al.(2018b) Humphrey et al.(2012) Abell 1689 Okabe et al.(2014b) RXC1825 Ghirardini et al.(2018b) Kawaharada et al.(2010) ESO 3060170 Su et al.(2013) Abell 1750 Bulbul et al.(2016) UGC 03957 Th¨olken et al.(2016) Abell 1795 Bautz et al.(2009)

r200. Table1 tabulates all of the works available on the outskirts of these objects.≈ Of particular interest is the radial profile of the gas entropy K(r), which is related to the gas temperature (kT ) and the electron density (ne) by K = 2/3 kT/ne . Numerical simulations of the formation of galaxy clusters (Voit et al. 2005) which include only gravitational physics have found that as clusters grow their entropy increases radially as a powerlaw. When clusters are scaled by the self similar entropy K500, and the scale radius R500c, these simulations predict that the entropy profile should follow the equation K(r)/K500 = 1.1 1.47(r/R500c) . Since these simulations include only gravitational physics, any observed deviations from this baseline entropy profile are interpreted as evidence for non-gravitational physics. Studies of the central regions of clusters have found the entropy to exceed this baseline prediction (e.g. Pratt et al. 2010), most likely due to energy dissipated into the ICM through AGN feedback and mergers. By comparing the outskirts entropy with the baseline profile, we can begin to search for the presence of physics beyond just simple gravitational collapse. Samples of the X-ray brightest clusters observed by Suzaku in the cluster outskirts have been compiled to investigate general trends. In early attempts 20 Stephen Walker et al.

(Walker et al. 2012b), where the entropy profile shapes were simply compared with one another, it was found that most Suzaku clusters had entropy profiles following a roughly similar shape, flattening at around R500c before unphys- ically turning over. In Walker et al.(2013), the outskirts for a sample of 13 X-ray bright Suzaku clusters were compiled, and their temperature, density and entropy profiles were compared to the predicted profiles from simulations which follow the baseline entropy profile from purely gravitational collapse (Voit et al. 2005) and the universal pressure profile (Arnaud et al. 2010). These profiles were appropriately scaled by mass and redshift to factor out the self similar behaviour. Here we expand upon the analysis of Walker et al.(2013) by including all the clusters studied by Suzaku in the outskirts published to date, which increases the sample size to 24 clusters (Fig. 10 and Fig. 11). The Suzaku entropy profiles (Fig. 10, top left) rise from the core, before flattening and intersecting the baseline entropy profile around R500c. Between R500c and r200 the entropy profiles remain flat. For 11 objects, entropy measurements beyond r200 are possible. When the sample is divided into mass bins (Fig. 10, top 14 right), we see that clusters with M500 > 2 10 M tend to lie well below the baseline prediction for the entropy profile× in these regions (by a factor of 14 around 0.5). In contrast, low mass clusters and groups with M500 < 2 10 M × tend to have entropies consistent with the baseline profile beyond r200, though the number of such low mass systems in the sample is small. The entropy excess over the baseline prediction at small radii is much more pronounced in these low-mass systems, presumably because, in shallower gravitational potential wells, the impact of AGN heating is relatively more important than in high- mass clusters. The observed pressure profiles for non-merging clusters tend to be in rea- sonable agreement with the universal pressure profile derived in Arnaud et al. (2010) from XMM-Newton observations of the central regions of clusters combined with simulations (the solid black line in the bottom panels of Fig. 10). They also agree well with the range found by stacking Planck clusters in Planck Collaboration et al.(2013b), shown by the shaded cyan region. For clusters where merging is taking place (blue points), the pressure profiles tend to lie above the universal pressure profile, though the scatter is large. When the self similarly scaled temperature profiles are compared to the predicted profile (Fig. 11, top panels) the temperatures out to r200 are generally in good agreement with expectations (the solid black line). The self similarly scaled density profiles on the other hand (Fig. 11, bottom panels) tend to be higher than the expected profile. This causes the entropy profiles to lie preferentially below the baseline entropy profile between R500c and r200. This tendency for the gas density to be higher than predicted is qualitatively consistent with the cluster outskirts beyond R500c being increasingly clumpy. If we have cooler, denser, unresolved gas clumps in pressure equilibrium with the hotter, more diffuse intracluster medium, this could account for the high densities and low temperatures (and hence low entropies) that have been ob- The Physics of Galaxy Cluster Outskirts 21

Fig. 10 Figures adapted and expanded from Walker et al.(2013) to include all the Suzaku outskirts results to date. Top panels: Entropy profiles scaled by self-similar entropy at R500c, and compared to the baseline entropy profile (black line). Bottom panels: The self-similarly scaled Suzaku pressure profiles are compared to the universal pressure profile (black line) and the range of the Planck cluster pressures profiles (cyan region). In the left hand column all the clusters are plotted with individual plot symbols. In the right hand column they are 14 14 grouped into clusters with M500 > 2 10 M ,M500 < 2 10 M , and merging clusters. × × served. This is discussed further in section 3.3.1 on observational constraints on gas clumping. 14 For the massive clusters (M500 > 2 10 M ) where a temperature mea- × surement is possible outside r200, the scaled temperatures tend to lie below the expectations from simulations (Fig. 11, top right), driving the entropy down even further (Fig. 10, top right). By contrast, the scaled temperature mea- 14 surements for low mass systems (M500 < 2 10 M ) and merging systems × outside r200 tend to be higher, and in better agreement with theoretical expec- 22 Stephen Walker et al.

Fig. 11 Figures adapted and expanded from Walker et al.(2013) to include all the Suzaku outskirts results to date. Top panels: The self-similarly scaled Suzaku temperature profiles in the outskirts (outside 0.4r200). Bottom panels: Self-similarly scaled Suzaku density profiles in the outskirts. The solid black curves on the left panels show the predicted temperature and density profile shapes obtained by assuming both the baseline entropy profile and the universal pressure profile to hold true. In the left hand column all the clusters are plotted with individual plot symbols. In the right hand column they are grouped into clusters with 14 14 M500 > 2 10 M (red), M500 < 2 10 M (black), and merging clusters (blue). In the bottom right× panel we compare the Suzaku× densities to those obtained by stacking ROSAT (grey) and Chandra (shaded cyan) clusters.

tations. Larger samples of groups will need to be assembled in the future to determine whether this trend is signiﬁcant. Suzaku proved to be very adept at measuring the outskirts of groups, and even managed to measure temperatures in very poor groups well outside the cluster core. For instance Nugent et al. (2017) obtained temperatures at r330 and r680 for NGC 3402 and NGC 5129 The Physics of Galaxy Cluster Outskirts 23

30:00.0 N bkg. cluster

43:00:00.0 r200

30:00.0 E

42:00:00.0 M86 W M90 E 30:00.0 M60 IC 310 M58

41:00:00.0

NGC4325

30:00.0

40:00:00.0 30 arcmin

655 kpc 39:30:00.0 30:00.0 25:00.0 3:20:00.0 15:00.0 10:00.0

0 1 4 9 20 41 83 166 335 669 1333

103 ) 2 m c

V e k (

y p o r t n E 102

cold fronts

10-1 100 Radius (Mpc)

Fig. 12 Results from the Large Suzaku Projects on the outskirts of Coma (Simionescu et al. 2013, left), Perseus (Urban et al. 2014, centre) and Virgo (Simionescu et al. 2017, right), reprinted with permission. The Suzaku mosaics are shown in the top row. The bottom row shows the entropy profiles (the Coma figures include the density and pressure profiles as well). For Perseus and Virgo, the shape of the entropy profile agrees with the average best fit shape of the sample of clusters studied in Walker et al.(2012b) (shown by the black curve).

13 13 respectively, which have M200 values of just 1.164 10 M and 2.38 10 × × M respectively. The Suzaku density proﬁles, meanwhile, are generally in good agreement with the proﬁles obtained by stacking ROSAT (Eckert et al. 2012) and Chan- dra (Morandi et al. 2015) data, as shown in the bottom right of Fig. 11. Merg- ing systems, marked in blue, tend to have higher densities than non merging systems in the outskirts. Large, deep Suzaku programs covering the most nearby, X-ray bright clusters (Perseus – Simionescu et al. 2011; Urban et al. 2014, Centaurus – Walker et al. 2013, Coma – Simionescu et al. 2013, and Virgo – Simionescu et al. 2017), have provided the most detailed view yet of individual systems (Fig. 12). The azimuthal variations seen in Perseus in particular have shown how dramatically a ‘relaxed’ cluster can deviate from spherical symmetry, with large scale gas sloshing (Simionescu et al. 2012; Walker et al. 2018) reaching out to roughly half the virial radius, and possibly much further out. In Perseus this large scale sloshing appears to be a continuation of the sloshing behaviour seen in the cluster core (Walker et al. 2017). These result show that the gas sloshing caused by minor mergers commonly seen in cluster cores (Markevitch and Vikhlinin 2007, ZuHone et al. 2011), which results in a spiral pattern of 24 Stephen Walker et al.

alternating cold fronts, can extend out to extremely large radii. Such large scale gas sloshing out to half the virial radius has also been identified in Abell 2142 (Rossetti et al. 2013), RXJ2014.8-2430 (Walker et al. 2014) and Abell 1763 (Douglass et al. 2018). The highest signal to noise in the outskirts with Suzaku is achievable for the Perseus cluster, the brightest cluster in the X-ray sky, whose large angular extent allows detailed thermodynamic profiles to be studied even with Suzaku’s large PSF (a half power diameter of 2 arcmin). From an initial two strips from the core to the outskirts (Simionescu∼ et al. 2011), the mosaic has grown to have 8 strips equally spaced around the cluster (Urban et al. 2014) (Fig. 12, centre column). The azimuthally averaged entropy profile for Perseus follows the same general trend for the sample of clusters studied in Walker et al.(2013), namely flattening at R500c, and remaining below the baseline entropy profile outside r200 (Fig. 12, centre bottom panel). When the density profiles are considered, these again tend to lie systematically above the predicted level. For the Virgo and Coma clusters, on the other hand, the entropy at r200 is similar to the predictions from the baseline entropy profile, while it fails to match the theoretical expectations at lower radii, presumably at least in part because of the unrelaxed dynamical state of these systems. Interestingly, the shape of the entropy profile in Virgo agrees with the empirical shape found in Walker et al.(2012b) by fitting the profiles of a sample of X-ray bright clusters, but with a higher overall normalization; instead of following the baseline entropy power-law shape at small radii and flattening out below it at large radii, in Virgo, the entropy is above the baseline at almost all radii, and flattens out to meet the expected profile at r200 (Fig. 12, lower right panel). Since one cannot rule out the possibility that the assumed baselines of radial profiles of X-ray observables may mislead or bias our understanding of thermodynamics in the cluster outskirts, Okabe et al.(2014b) have carried out a joint X-ray and WL analysis for four clusters which were observed by both the Suzaku satellite and the Subaru telescope. Because WL mass measurements are free from cluster dynamical states, this is a direct and powerful way to compare two independent observables. By simultaneously fitting the radial profiles of X-ray observables with WL masses, they found that the temperature sharply drops outside 0.5r500 and the density slope is slightly shallower (Figure 13), resulting in∼ an entropy profile that is flat in the cluster outskirts. The normalized pressure profile (normalized using equations5) is in good agreement with the Planck pressure profile. As the Suzaku data sometimes only cover part of the cluster azimuth at the virial radius, it is possible that some of the measurements are partially affected by accretion flows from the surrounding large-scale structure, which produce azimuthal variations of the outskirts temperature. Kawaharada et al.(2010) and Ichikawa et al.(2013) have investigated the correlation between outskirts temperatures obtained with Suzaku and the galaxy number density from the SDSS DR7 catalog for the very massive clusters A1835 and A1689. The number densities of galaxies were computed by slices of photometric redshifts in order to include red and blue galaxies. They found spatial temperature variations in The Physics of Galaxy Cluster Outskirts 25

10-1 100 10-1 100 102 D 101 n/n T/T 100 100 Hydra A A478 -1 10 A1689 A1835 10-2 10-1 − − 102 101 P/P K/K 100 100

10-1 rr r 10-2 10-1 10-1 100 10-1 100 − − 10-1 100 10-1 100 r/r r/r

Fig. 13 Figures taken from Okabe et al.(2014b). Normalized radial profiles of gas number density (top-left), temperature (top-right), pressure (bottom-left) and entropy (bottom- right), obtained from joint analysis of the observed gas density and temperature profiles for four clusters which also have weak lensing data. The normalizations n and T are given ∗ ∗ by equations5, with P and K derived from these using P = n T and K = T /n2/3. The lower subpanels show∗ the deviations∗ from the respective∗ best-fit∗ ∗ profiles.∗ The∗ overall∗ amplitude of deviations is different from cluster to cluster, reflecting the intrinsic scatter between the X-ray observables and weak-lensing mass. The black vertical bars denote r2500, r1000, and r500, from left to right.

< < the cluster outskirts (R500c r rvir) and a correlation between the outskirts temperatures and galaxy number∼ ∼ density surrounding the clusters (14). This suggests that cluster outskirts are directly affected by mass accretion flows from the surrounding large-scale environment. In order to fully understand the environmental effects and cluster thermal evolution, more data to cover large samples of clusters with full azimuthal coverage are vitally important. Using ROSAT data alone, Eckert et al.(2012) explored the density profiles of a sample of 31 local clusters, having full azimuthal coverage out to 26 Stephen Walker et al.

L − L −

′

− − - - - - - L ,

Fig. 14 Comparison of ratio of the temperature deviations (∆T ) from the azimuthally < < averaged temperature (TAll) in the cluster outskirts (r500 r rvir) and number density contrast δ = ∆n/ n of large-scale structure for the region∼ of (1–2)∼ rvir (left) and (2–4)rvir (right). Red circlesh andi blue diamonds represent the results for Abell 1835 and Abell 1689, respectively. Errors are the 68% CL uncertainty. Spearman’s rank correlation coeﬃcient, rs, gives that theses correlations do not occur by accident more than 90% (left) and 67% (right) conﬁdent. Figure is taken from Ichikawa et al.(2013), reprinted with permission.

r200, consisting of 14 cool core and 17 non-cool cores (shown in the bottom right hand panel of Fig. 11 as the solid pink region). They found a significant difference between the densities in the outskirts of cool core and non-cool cores: cool cores tended to have systemically lower densities in the outskirts than non-cool cores. When the observed density profiles were compared to both non-radiative numerical simulations and those including cooling and star formation, the densities in the outskirts were found to be higher than the simulations. The azimuthal scatter in the X-ray surface brightness was found to increase out into the outskirts, reaching around 60-70 percent at r200, with no significant difference between the scatter between cool core and non-cool core clusters in the 0.7-1.2r200 radial range. The ROSAT analysis was further developed in Eckert et al.(2015), in which it was found that the azimuthal median of the surface brightness provides a better way to recover the true underlying density profiles of clusters in the presence of inhomogeneities (such as gas clumping). Further recent progress has been made through the XMM Cluster Outskirts project (X-COP, Eckert et al. 2017), which uses deep XMM observations of 13 clusters with the highest signal to noise in the Planck SZ survey, combining X-ray surface brightness profiles with Planck SZ derived pressure profiles to obtain measurements of the temperature, entropy and mass profiles. Tchernin et al.(2016) presented results for the first cluster from the X-COP project, Abell 2142 (top panels of Fig. 15). Combining the Planck pressure profile with the non-clumping corrected density profile from XMM results in an entropy profile that flattens beyond R500c, similar to the Suzaku results. The Physics of Galaxy Cluster Outskirts 27

1.51e-04

27.7 X-ray/SZ, photometric median NE 7.58e-05 NW X-ray/SZ, photometric mean 27.6 X-ray spectroscopic data 3.80e-05 Purely gravitational collapse model 27.5 1

27.4 1.89e-05 500 K/K 27.3 9.49e-06

27.2 4.74e-06

Declination 27.1

2.34e-06 27.0 1 R/R 1.2 500

26.9 1.15e-06 th 1 SE /K 0.8 26.8 obs 0.6 SW 5.48e-07 K 0.4 26.7 0.2 2.50e-07 1 240.2 240.0 239.9 239.7 239.6 239.4 239.3 239.1 R/R Right ascension 500 1 00e 07

Fig. 15 The top two panels show figures from Tchernin et al.(2016), reprinted with permission. Top left: The XMM mosaic for A2142. Top right: Entropy profile obtained by combining the X-ray surface brightness profile from XMM-Newton with the Planck pressure profile for A2142. When clumping is not corrected for, the entropy profile flattens (cyan profile). When the median X-ray surface brightness profile is used to clumping correct the X-ray data (green profile), the entropy profile is in better agreement with the baseline entropy profile shown as the solid black line. The bottom two panels show figures from Ghirardini et al.(2018b), reprinted with permission. Bottom left:The XMM-Newton mosaics for the XCOP clusters. Bottom right: Clumping corrected entropy profiles for the X-COP sample.

When clumping is taken into account, using the azimuthal median density profile following Eckert et al.(2015), the resulting decrease in the outskirts density causes the entropy profile to continue to rise, in better agreement with the Voit et al.(2005) baseline entropy profile. Ghirardini et al.(2017) then presented results from Abell 2319 (the strongest Planck detection in X-COP with a signal to noise of 51). Again, when the non- clumping corrected density profile is combined with the Planck pressure profile, an entropy profile which flattens beyond r500 is found. However, in this case, even when the densities are clumping corrected, the entropy profile remains flat beyond r500. The pressure profile is also flatter and higher than the average Planck profile, and the universal pressure profile. These results suggest that, in this particular system, a non-thermal pressure support of about 40 percent of the total pressure is present at r200. Ghirardini et al.(2018b) have studied the thermodynamical profiles for the full X-COP sample consisting of 13 clusters, shown in the bottom panels of Fig. 28 Stephen Walker et al.

15. For 10 of the clusters these profiles extend out to r200, while the remaining 3 extend out to around 1.2r200. When clumping is taken into account using the azimuthal median technique (Eckert et al. 2015), the entropy profiles are mostly consistent with the Voit baseline entropy profile except for Abell 2319, which remains flat at r200 because of the high level of non-thermal pressure discussed in Ghirardini et al.(2017). The entropy profile for Abell 644 also appears flat, but the errors are large enough to be consistent with the baseline entropy profile.

3.3 Virialization of the accreted gas

3.3.1 Clumping

As described in section 3.2, accurate measurements of the thermodynamic profiles in the outskirts of clusters require us to take into account the observational bias resulting from gas clumping. Whilst clumping presents a complicating factor in measuring the properties of the bulk of the intracluster medium, being able to accurately measure the level of gas clumping also provides a powerful test of numerical simulations of cluster physics in the outskirts. As a result, much work, both theoretical and observational, has been directed at trying to understand the level of clumping in the cluster outskirts. In the first Suzaku measurements of the outskirts of the Perseus cluster, Simionescu et al.(2011), found the gas mass fraction to rise well above the mean cosmic baryon fraction in the outskirts. The authors estimated the mass profile near the virial radius by extrapolating a parametrized NFW profile constrained primarily by the data at smaller radii, in order to avoid significant bias due to the possible deviations from HSE discussed in Section 1.4. It was concluded that gas clumping must be biasing high the measurements of the gas density at the cluster’s edges. Because the X-ray emission is proportional to the square of the gas density, unresolved dense clumps will increase the observed X-ray luminosity and bias the inferred gas density high by a factor of √C, where C is the ‘clumping factor’ defined in Equation1. Since √C is what is most directly found from observations, this is the quantity we will discuss. By requiring that the true gas mass fraction equal the mean cosmic baryon fraction, an early estimate of this ‘clumping factor’ was obtained, which rose steadily from from unity at 0.5r200, to √C 4 at r200 (Simionescu et al. 2011). This early result was, however, based on' only one arm of the Perseus mosaic, which probed only a small fraction (around 5 percent) of the cluster azimuth. Walker et al.(2013) assembled a sample of the X-ray brightest clusters studied in the outskirts with Suzaku and explored their self-similar scaled thermodynamic profiles. As already shown in Fig. 10, the entropy profiles of the clusters lie systematically below the baseline entropy profile for purely gravitational hierarchical structure formation from around 0.5r200 outwards. Deviations from this baseline entropy profile must be due to non-gravitational The Physics of Galaxy Cluster Outskirts 29

R200 2 1.6 Data 1.8 V13 1.5 R13 NR

R13 CSF+AGN 1.4 1.6 R13 Residual

1.3 C 1.4 C 1.2

1.2 1.1

1 1

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.9 103 R/R500 R [kpc]

Fig. 16 Top left: Modified figure from Walker et al.(2013) showing the clumping overestimate √C needed to bring the observed Suzaku density profiles into agreement with the density profile expected from theory if the pressure and entropy agree with their baseline profiles (merging clusters have been excluded from this plot). Top right: clumping level along the 8 different arms in the Perseus Cluster Suzaku mosaic estimated by Urban et al. (2014) using the same method. Bottom left: Clumping factors obtained in Eckert et al. (2015) by comparing the median and mean X-ray surface brightness profiles of ROSAT clusters (black points). The surves show the predictions of several sets of numerical simulations: non-radiative AMR (red, Vazza et al. 2013), non-radiative SPH (blue), SPH including cooling, star formation and AGN feedback (green), and “residual” clumping after masking substructures (dashed blue; all from Roncarelli et al. 2013). Bottom right: Clumping in A2142 determined by Tchernin et al.(2016) comparing the median and mean X-ray surface brightness in XMM-Newton X-COP data. The dashed and dot-dashed vertical lines show r500 and r200 respectively, while the inner triple dot-dashed line shows the radius out to which large scale gas sloshing reaches. Figures reprinted with permission.

processes. If it is assumed that gas clumping is the only observational bias affecting our measurements of the outskirts, it is possible to estimate the level of gas clumping by calculating the overestimate of the gas density needed to bring the entropy profiles back into agreement with the baseline entropy profile. Walker et al.(2013) calculated these estimated clumping factor profiles for the Suzaku clusters. Again, √C tends to increase from unity at 0.5r200 ∼ up to around a factor of 1.5-2 at r200. 30 Stephen Walker et al.

Outside r200, the estimate of √C increases to much higher values, ranging from 2-10. However, when Walker et al.(2013) compare the temperature and density profiles against their theoretical predications (Fig. 11), they find that outside r200 the entropy decrement below the baseline entropy profile is driven not only by gas clumping but also by the gas temperatures being lower than expected. This will therefore lead to an overestimate of the gas clumping factor if all of the entropy decrement below the baseline entropy profile is assumed to be caused by an overestimate of the gas density. When the gas densities alone are compared to theoretical expectations (Fig. 11, bottom panels), we see that the densities typically lie at a constant level above the theoretical expectation, with an overestimate (i.e. a √C estimate) in the range √C 1-2. This clumping factor estimate is plotted in the top left hand panel of Fig.≈ 16. A similar method of estimating the level of gas clumping by solving for the density overestimate that would bring both the observed entropy and pressure in agreement with their expected baseline profiles was applied to the full 8-arms Suzaku Key Project mosaic of the Perseus Cluster by Urban et al. (2014). They found that the northwestern arm for which clumping was measured in Simionescu et al.(2011) showed the second highest clumping factor of the eight directions considered (see top right panel of Figure 16), and that the average clumping level over the 8-arms azimuth had a lower value of √C 1.6 ≈ at r200, in good agreement with other clusters in the Walker et al.(2013) sample. For the Coma and Virgo clusters also covered by Suzaku Key and Large Projects, the situation is complicated by the disturbed dynamical nature of these systems. The Planck and Suzaku pressure profiles for the Coma Cluster appear to be in agreement with each other out to r200, and no significant over- estimation of the X-ray gas density through clumping is required by the data between R500c and r200 (Simionescu et al. 2013). In the Virgo Cluster, the average pressure profile obtained from Suzaku does lie higher than the Planck measurement for the same system, concomitant with a marked flattening of the entropy profile; qualitatively this suggests clumping must be important in the Virgo Cluster but, due to the unrelaxed nature of the system, a quanti- tative estimation of the clumping level and disentangling the contributions of clumping versus non-thermal pressure is difficult (Simionescu et al. 2017). With the higher spatial resolution and full azimuthal coverage afforded by ROSAT, Eckert et al.(2015) found that is possible to estimate the clumping factor by comparing the median X-ray surface brightness profiles of clusters with the mean profile. The azimuthal median can be optimally estimated by adaptively binning the X-ray image with a target binning of 20 counts per bin. The median of the surface brightness distribution in each annulus was found to be insensitive to the presence of outliers. The resulting clumping factors for the ROSAT clusters are shown in the bottom left hand panel of Fig. 16. From 0.5r200 to r200 these rise gradually from unity to a range of √C 1.1-1.5, agreeing with the range of √C 1-2 found in the Suzaku estimates≈ obtained ≈ using entropy profile arguments. Outside r200 the ROSAT clumping factors are in the range √C 1-1.6. This range is reasonably consistent with the range of clumping factor≈ that is needed to bring the observed Suzaku densities into The Physics of Galaxy Cluster Outskirts 31

agreement with the theoretical expectations (Fig. 16, top left hand panel). Eckert et al.(2015) also compared the measured clumping factor profiles with the predictions of several sets of numerical simulations (Roncarelli et al. 2013; Vazza et al. 2013). Models including baryonic physics (cooling, star formation and feedback processes) provide an excellent match to the observed clumping factor (see also Planelles et al. 2017). Conversely, non-radiative simulations tend to overpredict the level of clumping in the ICM. The combined effect of cooling and AGN feedback is responsible for the difference between the various baryonic physics models. Cooling removes the most structured gas phase from the X-ray emitting phase in the central regions of infalling halos, whereas AGN feedback pushes the gas away from and smoothes the density profiles, thus lowering the level of inhomogeneities in the ICM. Using this method of comparing the median and mean X-ray surface brightness profiles, the X-COP project has produced constraints for individual clusters, using XMM-Newton observations of the outer regions of clusters which have the highest signal to noise ratio in the Planck cluster catalog. The profile for √C obtained for Abell 2142 is shown in the bottom right hand panel of Fig. 16. This profile again increases from unity around R500c, reaching √C 1.4 ' just outside r200. With deep high spatial resolution observations of the cluster outskirts using Chandra, it is possible to put direct constraints on the level of gas clumping by exploring 2-dimensional fluctuations in the X-ray surface brightness, and searching for excess variance over the Poisson noise. Morandi et al.(2013) employed this method to study gas clumping in Chandra observations of Abell 1835, and later in Morandi and Cui(2014) used the same method on deep observations of the outskirts of Abell 133. Morandi et al.(2017) also used this method to study clumping in the outskirts of the galaxy group NGC 2563. In all cases a clumping factor of around 2-3 (implying √C=1.4-1.7) at r200 is found. Due to the need to bin the data to achieve sufficient signal to noise, gas clumping on very small scales (less than 10kpc in the case of Abell 1835) cannot be constrained using this method, so these estimates are lower limits to the true level of gas clumping. In principle it should be possible to directly image gas clumps with suffi- ciently deep observations having the requisite sensitivity. For example, in Ichi- nohe et al.(2015) three possible gas clumps have been identified in the merging cluster Abell 85 using Chandra observations. These gas clumps, whose X-ray luminosities are in the range 1-2 1040 erg s−1, show features of ram-pressure stripping. × It is, however, difficult to discover such X-ray faint diffuse emission from gas clumps with existing instruments because it requires deep mosaicked observations. One of the most efficient observation strategies is to use optical and/or weak-lensing information. Subaru weak-lensing mass maps covering most of the Coma cluster have revealed the presence of many dark matter subhalos 12 with masses above a few 10 M (Okabe et al. 2014a). Tangential shear profiles stacked over the dark matter subhalos show sharp drops outside specific radii, suggesting that the subhalos are tidally truncated and located inside the 32 Stephen Walker et al. clusters. Thus, they are remnants of group-size objects. Dark matter subhalos are associated with known optical groups. X-ray follow-up studies (Andrade- Santos et al. 2013; Sasaki et al. 2015, 2016) using Chandra, XMM-Newton and Suzaku discovered evidence that diffuse X-ray emission is associated with some of the subhalos. Dark matter subhalos are also associated with the known X-ray gas clump of NGC 4839 which shows ram-pressure striping. Despite early disagreements and uncertainties, it is remarkable that the levels of clumping inferred from very different methods (Suzaku thermodynamic profiles, the mean versus median X-ray brightness in ROSAT and XMM- Newton, and direct Chandra imaging) seem to give consistent ranges for the gas clumping factor in recent works. These measurements agree well with the level of clumping seen in numerical simulations. In order to fully understand the interplay between gas clumps and the ICM, multi-wavelength approaches combining X-ray, optical and weak-lensing are essential to future studies.

3.3.2 Non-equilibrium ionization

As described in section 1.5, the ICM in the cluster outskirts can be in a non-equilibrium ionization state, with the ion temperatures being higher than the electron temperatures. Since our X-ray spectroscopic measurements using CCDs only allow us to observe the electron temperature at present, this could cause our temperature measurements in the outskirts to be biased low compared to the true gas temperature. Simulations indicate that this effect could lead to up to a 10-15 percent underestimate of the gas temperature at R200c for massive clusters (see Fig.6, where R200c is roughly 0.5R200m). It is possible that this may contribute to the lower than expected temperatures found by Suzaku beyond R200c in massive clusters (Fig. 11, top panels). In early works, Hoshino et al.(2010) and Akamatsu et al.(2012a) at- tempted to examine the possibility of such an ion-electron two temperature structure in the clusters A1413 and A3667 respectively using Suzaku data. In these works the authors compared the observed entropy profile shapes with the baseline r1.1 powerlaw shape. If it is assumed that any flattening from this profile shape is due to the gas temperature being underestimated due to non-equilibrium ionization, then this can provide a basic upper limit on the magnitude of this effect. This approach does however omit to factor in the effects of gas clumping on flattening the entropy profile shape. At present, no firm observation of the presence of non-equilibrium ionization conditions has been possible in the outskirts of clusters.

3.3.3 Non-thermal pressure support

The presence of non-thermal pressure in cluster outskirts induces an additional pressure contribution to balance gravity and biases total masses estimated assuming hydrostatic equilibrium towards low values1 (Lapi et al. 2010; Kawa-

1 A detailed discussion of hydrostatic mass reconstruction techniques and the corresponding systematics and biases is described in another review of this series (Pratt et al.). The Physics of Galaxy Cluster Outskirts 33

harada et al. 2010; Ichikawa et al. 2013; Ettori et al. 2013a; Fusco-Femiano and Lapi 2014). The level of this non-thermal pressure support can be evaluated by comparing measurements of the total mass obtained through the hydrostatic equilibrium equation and via other independent techniques (like gravitational lensing, virial theorem) or assuming that hydrodynamical cosmological simulations provide reliable estimates of the gas mass fraction in clusters. Subaru WL observations can cover the entire regions of galaxy clusters and can thus constrain M∆ where ∆ 2500 100. From a comparison of Suzaku ∼ − and Subaru masses at r200 for a small sample of four clusters, Okabe et al. (2014b) found that the mass bias (the mean mass ratio between hydrostatic equilibrium masses and WL masses) is on average 1 b 0.5 at ∆ = 200. − ∼ Due to its dependence on the gas and hydrostatic mass, the gas mass fraction reconstructed using the assumption of hydrostatic equilibrium is highly sensitive to the level of non-thermal pressure. Since the most massive clusters are expected to be fair archives of the universal baryon budget (e.g. White et al. 1993; Evrard 1997; Nagai et al. 2007; Sembolini et al. 2016), the effective level of non-thermal pressure support can be estimated by comparing hydrostatic measurements of the gas fraction with the cosmic value (Ωb/Ωm = 0.156 0.003, Planck Collaboration 2015) properly corrected for any baryon depletion± and contribution from cold (i.e. not observable in X-rays, like stars in galaxies) baryons. Generally speaking, the effects of gas clumping and non-thermal pressure on the gas fraction both induce a positive bias in the measurements of fgas in cluster outskirts. The interpretation of excess gas fractions alone is thus not unique. Recently, Eckert et al.(2018) applied the azimuthal median technique (see Sect. 3.3.1) to clean the estimated gas density profiles from the effects of clumping, which allowed the authors to disentangle clumping and non-thermal pressure support. The clumping-corrected profiles were used to reconstruct hydrostatic mass profiles with a precision of < 5% at R500c for all X-COP clusters (Ettori et al. 2018). In the left-hand panel of Fig. 17 we show the measured hydrostatic gas fractions for the 12 X-COP clusters (Eckert et al. 2018). With one exception, all the measurements lie very close to the expected hot gas fraction, which implies that the average level of non-thermal pressure in the sample is low. The corresponding levels of non-thermal pressure are shown in the right-hand panel. The median non-thermal pressure fractions in the sample are 6% and 10% at R500c and r200, respectively. A detailed comparison of the mass estimates (Ettori et al. 2018) also shows that X-COP hydrostatic masses agree with WL-based measurements for the same clusters within 10%, which sup- ports this conclusion. In the CLASH sample, Siegel et al.(2018) reached a similar conclusion using a combination of WL, X-ray and SZ data and set an upper limit of 11% on the level of non-thermal pressure at r500. 34 Stephen Walker et al.

0.24 Nelson+14 Universal fgas, 500 The300 A1644 fgas, 500 0.5 A1795 f 0.21 gas, 200 A2029 A2142 A2319 0.4 A2255 0.18 A2319 A3158 t E

o A3266

S / b m A2142 t A2255 H A644 A85 P

, 0.3 0.15 A3158 / A85 s A2029 A1795 T a RXC1825 RXC1825 A3266 N

g A644 A1644 f

P ZW1215 0.12 0.2

ZW1215

0.09 0.1

0.06 0.0 1015 0.6 0.8 1.0 1.2 1.4 1.6 M500, HSE [M ] R/R500

Fig. 17 Constraints on non-thermal pressure support in cluster outskirts in the X-COP sample (reproduced from Eckert et al. 2018), reprinted with permission. The left-hand panel shows the hydrostatic gas fraction corrected for gas clumping estimated at R500c (blue points) and r200 (red). The green shaded area shows the cosmic baryon fraction from Planck (Planck Collaboration 2015) whereas the magenta line and shaded area represents the expected hot gas fraction and uncertainties after accounting for baryon depletion and stellar mass. In the right-hand panel we show the inferred level of non-thermal pressure acting against gravity for the 12 X-COP clusters. The curves and shaded areas show the median non-thermal pressure and dispersion from two diﬀerent sets of numerical simulations (Nelson et al.(2014b), blue; Rasia et al. private communication, green).

4 What progress can future missions bring?

The last decade has brought about a transformation in our understanding of the outskirts of galaxy clusters. Owing to the extremely low density of the gas in the outskirts, these measurements remain highly challenging. Here we conclude by reviewing the capabilities of new planned and proposed missions across the electromagnetic spectrum, which will continue to unlock the physics in the outskirts of galaxy clusters.

4.1 Ongoing and planned missions

4.1.1 eROSITA

The extended Roentgen Survey with an Imaging Telescope Array (eROSITA) is the main instrument on the upcoming Russian Spektrum-Roentgen-Gamma (SRG) mission slated for launch in 2019 in an orbit around L2 (Predehl et al. 2010). With its large grasp, eROSITA will perform the first all sky X-ray imaging survey since ROSAT in the early 1990s, achieving a sensitivity 30 times higher than that of ROSAT in the soft X-ray band (0.5-2.0keV). It∼ will perform the first ever all sky X-ray imaging survey at high X-ray energies (2.0-10keV). It is expected that eROSITA will detect all galaxy clusters with 14 masses greater than 3 10 M in the observable universe (Pillepich et al. 2012). The dramatic∼ increase× in sample size afforded by eROSITA promises a transformational advance in our ability to use galaxy clusters as probes The Physics of Galaxy Cluster Outskirts 35 of cosmological parameters. Additionally, the eROSITA telescopes feature a high effective area to focal length ratio, which effectively reduces the particle background rate per unit area (which dominates the background level above around 2 keV) by about an order of magnitude compared to other instruments located on a high orbit (XMM-Newton, Chandra). The combination of wide field of view, moderately high angular resolution and low particle background will make eROSITA very sensitive to low surface-brightness regions such as cluster outskirts.

4.1.2 XRISM/XARM

The X-Ray Imaging and Spectroscopy Mission (XRISM), formerly XARM (the X-ray Astronomy Recovery Mission, Tashiro et al. 2018) is an upcoming NASA/JAXA collaboration (with contributions from ESA) to replace the ill- fated Hitomi satellite that was lost after only 37 days of operation in 2016. XRISM centers on its high spectral resolution X-ray microcalorimeter (Re- solve), which is identical to that flown on Hitomi. Running in parallel is a wide field X-ray imager (Xtend), with a field of view around 4 times that of Suzaku. XRISM is currently scheduled to launch in 2021. Importantly for cluster outskirts research, XRISM will launch into the same low Earth orbit as Suzaku, providing it with a low and stable background level. While observations of the cluster outer edges with the Resolve instrument would likely be prohibitively expensive in terms of the required exposure time, extensive studies of the thermodynamical properties in the outskirts of clusters can be carried out using Xtend, particularly in the later, cryogen-free phase of the mission. Xtends large field of view (40 40 arcmin) will provide a much more complete coverage of cluster outskirts,× minimizing the need for multiple mosaicked observations that represented the standard observing strategy for Suzaku. In addition, XRISM will also have a better point spread function than Suzaku (1.2 arcmin HPD compared to 2 arcmin). This means that more spatially resolved proles will be attainable, allowing us to produce moderate spatial resolution maps of cluster temperature, entropy and metal abundance to the virial radius, with high azimuthal coverage.

4.1.3 Athena

The Advanced Telescope for High Energy Astrophysics (ATHENA) is a major X-ray observatory being developed by ESA for launch in the early 2030s (the launch date is currently 2031). ATHENA will use new mirror technologies (silicon pore optics) to achieve a huge effective area of 1.4m2 at 1 keV with an angular resolution of 5 arcseconds for a low mass. The combination of a large effective area and moderate angular resolution allows ATHENA to avoid photon starvation, low sensitivity to point sources and source confusion. The unique cluster science enabled by ATHENA’s capabilities is discussed in detail in Ettori et al.(2013b), and briefly summarized below. 36 Stephen Walker et al.

The ATHENA mirror serves a high spectral and spatial resolution microcalorimeter (the X-ray Integral Field Unit, X-IFU), which will achieve an energy resolution of 2.5 eV with a spatial resolution of 5 arcseconds (over 10 times better than XRISM) over a 5 arcmin field of view (Barret et al. 2018). The effective area will be 150 times better than XRISM at soft X-ray energies (0.6 keV), and over 12 times better at 6 keV, near the Fe-K complex. This will allow high spatial resolution mapping of gas motions and turbulence in clusters throughout their volumes for the first time, and will reveal the contributions of non-thermal pressure support in the outskirts of clusters. By measuring the ion temperature directly from the line widths, it will be possible to characterize the election-ion two temperature structure in the outskirts. The high level of sensitivity of the X-IFU will enable measurements of X-ray line absorption and emission from both the cluster outskirts and the warm-hot intergalactic medium (WHIM) permeating large-scale structure filaments. ATHENA also features a wide field X-ray imager (WFI), with a field of view of 40 40 arcmin (Meidinger et al. 2017). The combination of a large field of view and× unprecedented effective area will allow the outskirts of clusters to be routinely unlocked for X-ray spectroscopy, allowing temperature and entropy measurements out to the virial radius and beyond for large samples of clusters. The WFI instrument will feature a low and stable instrumental background, such that it will be extremely well suited for the study of cluster outskirts. The capabilities of ATHENA/WFI are highlighted in Fig. 18, where we show a simulated temperature profile for a local galaxy cluster with a mean temperature of 8 keV at z = 0.1, i.e. similar to Abell 2142 and PKS 0745. This fiducial cluster is assumed to follow Universal profiles of temperature and emission measure in the radial range [0 1.2]r200. We used the Reiprich et al.(2013) parametrization of the temperature− profile, which was calculated from a linear fit to a compilation of Suzaku temperature profiles in the radial range [0.3 1.15]r200. For the Universal emission measure profile, we used the mean self-similar− scaled emission measure determined by Eckert et al.(2012) from a sample of 31 nearby clusters observed with ROSAT/PSPC (see Table C.1 of that paper). The metal abundance of the cluster was assumed to be flat with a mean value of 0.3Z . Our simulations show that 100ks observations of local clusters with the WFI will be able to measure the gas temperature with a precision of just 3% at r200. Given the expected stable instrumental background level (reproducibility of 2%), systematics in the reconstruction are expected to be minimal.

4.1.4 Upcoming SZ instrumentation

The current capabilities and science goals for many of the current and future instruments for SZ observations are discussed in detail in the review ‘Astro- physics with the Spatially and Spectrally Resolved Sunyaev–Zeldovich Eﬀects’ in these proceedings (Mroczkowski et al. 2018). We brieﬂy summarize some of the upcoming survey experiments and instrumentation here, with most of the major facilities listed in Table2. There are two broad categories for SZ TheThe Hot and Physics Energetic Universe: of The Galaxy Astrophysics ofCluster galaxy groups and Outskirtsclusters with Athena+ 37

9 Simulated WFI data Abell 2142 8 spectrum at Input profile

R200 with 7 Background +/-2% = 22’ Athena+ WFI 200 6 R

5 kT [keV] 4 3 2 1 Abell 2142 0 0.2 0.4 0.6 0.8 1 1.2 R/R200

Fig.Figure 2: The 18 combined X-IFU and WFI exposures will provide a temperature profile that will improve by an order of magnitude Left: the present constraintscapabilities on the thermodynamical properties of ATHENA/WFI of the ICM out to R500 and beyond. This in plot cluster presents a XMM outskirts.mosaic XMM-Newton image of of A2142 at z=0.09 (right) and a 100 k-sec WFI exposure of a region of 10 arcmin2 across the virial radius (left): the gas emissivity, A2142temperature and( metallicityTchernin are constrained et withal. relative 2016 uncertainties). of The2, 3 and 18green per cent, respectively, circle at 90%shows confidence (see the position of r200 220. The inset also Appendix). 2 ∼ shows a simulated 100ks spectrum for a 10 arcmin region around r200; the red curve is the source spectrum (reproduced from Ettori et al. 2013b), reprinted with permission. Right: High-resolution spectroscopy from the Athena+ X-IFU will resolve the emission complex around 6.7 keV arising from SimulatedHe-like Fe lines and their 100ks satellites, allowing ATHENA/WFI us to use these lines as strong temperature and independent temperature profile diagnostic of for a fiducial 8 keV galaxy cluster at the ICM in the range 107-108 K (Porquet et al. 2010). This measurement is particularly important to verify the zstandard= 0assumption.1. The that the inputICM is a tenuous temperature gas in collisional ionisation profile equilibrium. is This shown may not in fact in be blackthe and the simulated WFI data as case in the cluster outskirts, due to the low density and short timescales. In particular, the estimate of the intensity ratio thebetweenred Fe XXV data and Fe XXVI points. Kα line emissions The trace upperthe ionization andequilibrium lower state that can blue be altered dotted by ongoing lines show the approximate level of mergers (e.g. Akahori & Yoshikawa 2010). In the shocked regions of the ICM, departures are also expected from the systematicsassumed Maxwellian electron expected distribution. In case in of supra-thermal the measurements emission, this effect can be resolved assuming via enhanced a 2% reproducibility of the particle background.equivalent widths of the Fe XXV j-satellite (Kaastra et al. 2009). Through spatially resolved high-resolution spectroscopy, Athena+ will measure intra-cluster gas motions and turbulence, showing how the baryonic gas evolves in the dark matter potential wells, and probe the true nature and physical state of the X-ray emitting plasma. instrumentation:4. HOW AND WHEN survey WAS THE instruments ENERGY IN built THE forICM studies of the CMB, and user GENERATED AND DISTRIBUTED? instrumentsWithin the cluster atmosphere, to it has whichbeen recognised scientiststhat non-gravitational processes, can particularly propose galaxy feedback observations. There are currently from outflows created by supernovae (SN) and super-massive black holes, must play a fundamental role in the history 2of all major massive galaxies SZas well as survey the evolution of groups instruments and clusters as a whole. F ineedback operation, from AGN emission and the Advanced ACTpol camera stellar winds is likely (i) to produce the galaxy red sequence, (ii) to provide the extra energy required to keep large quantities of gas in cluster cores from cooling all the way down to molecular clouds, (iii) to account for the energy (i.e. onentropy) the excess observed 6-meter in the central Atacama regions of the group and Cosmology cluster ICM (e.g. Sun et al. Telescope2009; Pratt et al. 2010). (Thornton et al. 2016; Crowley Understanding the process of energy injection, and whether the energy was introduced early in the formation of the etfirst haloes al. (with 2018 further ),consequences and for SPT3G galaxy formation history), on or the gradually 10-meter over time by AGN feedback, South SN Pole Telescope (Benson et al. driven galactic winds, or an as-yet unknown physical process, is crucial to our understanding of structure formation and 2014evolution. The; Bender specific case of AGN et feedback al. in 2016 local (z < 0.5); Andersonsystems is discussed further et in Croston, al. 2018Sanders et a.). Both achieve arcminute res- (2013, Athena+ supporting paper) The various feedback processes, coupled with cooling, affect the overall ∼ olutionPage 4 in their SZ decrement bands, roughly centered at 90 and 150 GHz. Future upgrades include the Simons Observatory and CMB-S4, both of which can improve results with stacking of the large cluster samples they will at- tain (> 105 candidates spanning all redshifts at which bona fide clusters have formed; see The Simons Observatory Collaboration et al.(2018); Abazajian et al.(2016)). Targeted observations at subarcminute resolution with, for example, upgrades to ALMA (Kitayama et al. 2016; Basu et al. 2016), MUSTANG-2 (Dicker et al. 2014), NIKA2 (Catalano et al. 2018) and new instruments such as TolTEC (Austermann et al. 2018; Bryan et al. 2018) will enable more detailed mapping of the thermal and kinematic SZ effect. Spectrally resolved measurements with the intensity mapping experiments TIME (Crites et al. 2014) and CONCERTO (Lagache 2018) will allow better constraints on the product of velocity and electron opacity of the cluster gas through the kinematic SZ effect and on the mass-weighted gas temperature through relativistic corrections to the thermal SZ effect (Itoh et al. 1998; Chluba et al. 2012; Erler et al. 2018). Thus the next generation instruments are set to strengthen the role of the SZ effects in probing the ICM properties.

4.1.5 Upcoming optical/WL

On-going and future imaging surveys with limiting magnitudes less than about 25 are summarized in Table3. Deep and wide imaging surveys enable us to 38 Stephen Walker et al.

Name Frequencies Location Type (GHz) (S/T) Adv.ACTpola 30–270 6-m in Atacama Desert, Chile S SPT3Gb 95/150/220 South Pole Station, Antarctica S ALMAc 35–950 66 elements in the Atacama Desert, Chile T MUSTANG-2d 90 100-m Green Bank Telescope, USA T NIKA2e 150/260 IRAM 30-m, Spain T

TolTECf 150/220/270 50-m Large Millimeter Telescope, Mexico S/T CCAT-primeg 200-1000 6-m in Atacama Desert, Chile S ∼ Simons Obs.h 27-280 6-m in Atacama Desert, Chile S CONCERTOi 200-360 12-m APEX, Atacama Desert, Chile S/T TIMEj 183–336 12-meter ARO on Kitt Peak S/T

AtLASTk 35–950 Concept 50-m in Atacama Desert, Chile S/T ∼ CMB-S4l 30-300 Concept array of 1-m & 6-m telescopes S CSSTm 90–400 Concept 30-m in Atacama Desert, Chile S ∼ Table 2 A list highlighting a few current and near-future instruments capable of measuring the SZ effect. Currently fielded instruments are separated from upcoming and proposed instruments by horizontal lines. The primary observing mode – survey or targeted observation – is designated by an ‘S’ or ‘T.’ a) The 3rd generation of the Atacama Cosmology Telescope, known as the “advanced polarimetry” generation of ACT (Adv.ACTpol, Hilton et al. 2018). b) 3rd generation of the South Pole Telescope (SPT), a CMB survey instrument Benson et al. (2014); Bender et al.(2016). c) Atacama Large Millimeter/Submillimeter Array (ALMA), a 66 element interferometer on the Chajnantor plateau, 5000 meters above sea level (A.S.L). The first SZ results are reported in Kitayama et al.(2016); Basu et al.(2016). d) The 2nd generation Multiplexed SQUID/TES Array at Ninety GHz (MUSTANG-2, Dicker et al. 2014). e) The 2nd generation New IRAM KIDs Array (NIKA2, Catalano et al. 2018; Adam et al. 2018). f) TolTEC Austermann et al.(2018); Bryan et al.(2018). g) CCAT-prime 6-meter mm/submm survey telescope currently being constructed on Cerro Chajnantor, at 5600 meters A.S.L. Parshley et al.(2018); Stacey et al.(2018). h) The Simons Observatory (S.O.) is constructing a new, large aperature 6-meter CMB survey telescope to be added to its facilities in Cerro Toco, at 5200 meters A.S.L. (The Simons Observatory Collaboration et al. 2018). i) The CarbON CII line in post-rEionization and ReionizaTiOn (CONCERTO) epoch project, planned for the 12-meter Atacama Pathfinder Experiment (APEX) telescope at 5100 meters A.S.L. Lagache 2018). j) The Tomographic intensity mapping experiment (TIME, Crites et al. 2014), designed to measure [CII] emission at high redshifts. k) The At- acama Large Aperture Submillimeter/millimeter Telescope (AtLAST), a concept 50-meter class telescope currently under study by the community (Bertoldi 2018; Mroczkowski and Noroozian 2018). l) The Chajnantor Submm Survey Telescope (CSST), a concept 30-meter class survey telescope that would be led by Caltech (Padin 2014). m) CMB-S4 is a concept DOE ‘Stage IV’ dark energy / cosmology project soon to be proposed, and will likely include several 1–6 meter telescopes located at the South Pole and in the Atacama Desert (Abaza- jian et al. 2016). A 10-meter component to improve CMB-S4’s cluster finding capabilities is also under consideration. systematically carry out the science of cluster galaxies and weak-lensing. The KiDS, DES and HSC-SSP surveys are underway and have published initial science outputs. The LSST will have first light in 2020, and start the operations phase in 2022. We also mention other imaging surveys; CIFS and Euclid. The CIFS survey using the CFHT will cover 5000 deg2 with the r-band and 10000 deg2 with the u-band by Jan. 31st∼ 2020. The depth will be 24.1 mag∼ and 23.6 The Physics of Galaxy Cluster Outskirts 39

Survey Telescope Depth Area Seeing Status r (5 σ) [deg2] [arcsec] KiDS VST 25.2 1,500 0.7 Underway DES DES 25.0 5,000 0.9 Underway HSC-SSP Subaru 26.0 1,400 0.7 Underway LSST LSST 26.9 20,000 0.7 First light 2020

Table 3 Major optical imaging surveys. The depth is defined by the 5σ aperture magnitude within 200 diameter. mag, which is defined by 10σ aperture magnitudes for point sources within 200 diameter, for the r-band and the u-band, respectively. The Euclid satellite will be launched in 2021 and cover 15,000 deg2 in the near-infrared bands with a limiting magnitude of 24 mag. Its deep layer covers 40 deg2 with two magnitudes deeper than∼ that of the wide one. Euclid will be able to detect 14 all clusters with M200 > 2 10 M at 3-σ significance out to redshift 2. It is expected to find around× 60,000 clusters with 0.21. These deep and wide optical surveys will give us unique opportunities to discover a considerably large number of galaxy clusters. Furthermore, combined with cluster catalogs which will be obtained by future X-ray, SZ, optical or WL projects, the surveys will be able to constrain well the splashback radii and their redshifts and mass dependences using optical galaxy distributions and/or tangential shear profiles.

4.2 Mission concepts

4.2.1 Advanced X-ray Imaging Satellite: AXIS

The Advanced X-ray Imaging Satellite (AXIS) is a NASA Probe Mission concept (mission cost less than 1 billion USD) being considered for the 2020 decadal survey (Mushotzky 2018). The concept consists of a CCD imager served by a large effective area telescope (at least 10 times the effective area of Chandra) over a much larger bandpass of 0.1-16 keV, with high spatial resolution over a large field of view of around 25 25 arcmins (0.4 arcsec HPD on axis, increasing only slightly to 1 arcsec at× the edge of the field of view). It would fly in a low earth orbit, providing it with a detector background around 4-5 times lower than that of Chandra. For background limited observations, the figure of merit for the detectability of low surface brightness X-ray emission (effective area/(focal length2 background)) is a factor of 65 better than Chandra, and would be the highest× of any X-ray observatory ever flown (see Fig. 19). AXIS’s high spatial resolution is designed to allow it to complement ATHENA (in much the same way that Chandra’s high spatial resolution complements the lower spatial resolution but higher effective area of XMM-Newton). In the right hand panel of Fig. 19 we show the improvement in sensitivity of AXIS over Chandra in simulated 100ks observations of the 40 Stephen Walker et al.

0.5-5keV Chandra 0.5-5keV AXIS

r200 r200

Fig. 19 Left:Figures of merit for the detectability of low surface brightness X-ray emission, given by: eﬀective area/(focal length2 background), reprinted with permission. Right: Comparing 100ks simulated exposures of× cluster outskirts with Chandra (left) and AXIS (right). AXIS’s high sensitivity would allow it to routinely probe the 1-2r200 region and to view the connections between cosmic web ﬁlaments and clusters. Figures taken from Mushotzky(2018).

outskirts of a massive (8 keV) cluster at redshift 0.1. With AXIS the 1–2r200 region of clusters would be routinely opened up to X-ray spectroscopy and imaging.

4.2.2 Lynx

Lynx is a NASA large class X-ray mission concept being considered for the 2020 decadal survey (Gaskin et al. 2017). Lynx would feature the largest X- ray telescope ever built with a collecting area of 2m2 at 1 keV, over a 10 arcmin radius ﬁeld of view. Its PSF would be similar to Chandra’s on-axis (0.5 arcsec HPD), but will be sustained at less than 1 arcsec HPD over the whole ﬁeld of view. In addition to the imaging detector, an X-ray grating spectrometer and an X-ray microcalorimeter are planned for high spectral resolution work. Similar to Athena, Lynx would launch to L2. Together with Athena, Lynx’s high resolution spectroscopy capabilities can bring important further progress in understanding the dynamics and electron-ion equilibration balance in cluster outskirts.

4.2.3 HUBS

The Hot Universe Baryon Surveyor (HUBS2) is a Chinese led concept focusing on looking for the ‘missing baryons’ in intergalactic and circumgalactic space. The mission concept aims for a large effective area ( 1000cm2), 1 arcmin angular resolution and in particular a large field of view∼ ( 1 square degree) with high spectral resolution (2 eV at 0.6 keV) in the soft∼ X-ray band (0.1 to 2 keV). A superconducting Transition Edge Sensor (TES) will be employed ∼ 2 http://heat.tsinghua.edu.cn/~hubs/ The Physics of Galaxy Cluster Outskirts 41 as the main spectrometer, which will have a similar number of pixels to the Athena X-IFU but with a larger absorber. With a low Earth and low inclina- tion orbit giving it a low instrumental background, the mission aims to study the OVII (0.57 keV) and OVIII (0.65 keV) lines in both emission and absorption. HUBS plans to perform three different observational programs: deep pointings, a medium survey, and an all-sky survey. As a secondary objective, HUBS will also investigate wider science cases including the hot interstellar medium, the diffuse X-ray background, supernova remnants, as well as charge exchange processes in the solar system.

4.2.4 Super-DIOS

Super-DIOS (Diffuse Intergalactic Oxygen Surveyor) is a Japanese led concept (Ohashi et al. 2018) for a high spectral resolution X-ray spectroscopy mission in the early 2030s, which aims to detect the WHIM through its redshifted emission lines from OVII, OVIII and other ions. The mission concept features a high energy resolution of 2 eV at 1 keV over a large field of view of 30 arcmin, with a spatial resolution of 10 arcsec over a large effective area ( 2000cm2 at 1.0 keV). The number of pixels for the TES detector would be∼ 10 times that of Athena. The science with Super-DIOS has its base in the DIOS mission concept (Ohashi et al. 2016) and will extend it further. The main purposes of the mission concept are as follows: 1. Detection and survey of the inter/circum-galactic medium in the redshift range z < 0.3 over 5 5 degrees. 2. Characterization of the× ICM around the virial and splash-back radii. 3. Probing the metal synthesis in the early Universe using X-ray afterglows of gamma-ray bursts as background light. The large effective area, low background, high spectral resolution, and large field of view of both HUBS and Super-DIOS are ideal for detailed studies of the thermal, dynamical, and chemical structure in the outskirts of nearby clusters of galaxies, and the surrounding large-scale structure filaments that connect to them.

4.2.5 AtLAST

The Atacama Large Aperture Submillimeter/millimeter Telescope (AtLAST) is a 50-meter class wide ﬁeld of view (> 1 deg2) single dish telescope currently under study by the community.3 This would provide coverage of SZ/CMB bands, complementing ALMA both by mapping large regions > 105 faster than ALMA can, and by recovering spatial scales from 5 arcseconds× to a few degrees. ALMA is currently limited to recover < 10 scales,∼ while expansion to lower frequencies will only improve this by a factor∼ of 2.5 . AtLAST will be able to measure cluster masses through CMB lensing, while× also simultaneously providing (sub)millimeter constraints. Such constraints are important

3 http://atlast-telescope.org 42 Stephen Walker et al.

for reliable background subtraction in sensitive kinetic and thermal SZ studies, as discussed in the SZ review chapter in these proceedings.

5 Summary

5.1 What have we learned about the physics of galaxy cluster outskirts in the last 5 years? ∼ The last 5 years have seen great progress in our understanding of cluster outskirts,∼ both in theory and in observation. Simulations have helped to re- fine the way we define the outer boundaries of clusters by showing there to be a ‘splashback radius’ in the collisionless dark matter distribution, which lies within the accretion shock of the collisional gas. Significant observational evidence has now been found to support the prediction of a splashback radius. A growing diversity of methods for exploring the intracluster medium out to r200 has flourished and matured, combining X-rays, the Sunyaev Zeldovich effect and weak lensing. Suzaku has probed the outskirts of 35 galaxy clusters out to at least r200, from low mass groups to massive clusters, allowing samples to be compiled and trends to emerge for the cluster population as a whole. The density profiles tend to be higher than theoretical expectations, consistent with the expected gas clumping bias. For 11 of these clusters measurements are possible beyond r200. The Suzaku results suggest that galaxy groups and galaxy clusters may have different behaviours outside r200, with clusters tending to have steeper normalized temperature profiles and lower normalized entropies than groups, however a larger sample is needed to confirm this trend. A growing sample of clusters (13) has also been studied in the outskirts by combining Planck SZ observations of gas pressure with XMM-Newton observations of gas density. These studies have also found gas clumping to play an important role in cluster outskirts, which must be taken into account to find the true thermodynamic profiles of the intracluster medium. Extremely deep observations of cluster outskirts with Chandra have also revealed the presence of excess surface brightness fluctuations due to gas clumping. These different methods of probing the cluster outskirts (Suzaku, XMM/Planck and Chan- dra) have now converged on similar values for this level of gas clumping, and agree with theoretical expectations. Evidence for non-thermal pressure support in the outskirts of some clusters has been found, reaching a level of 10 percent at r200, and has been successfully disentangled from the gas clumping bias. This helps to deepen our understanding of possible biases in the measurement of cluster masses out to large radius. Detailed studies of nearby clusters such as Perseus have shown that ‘relaxed’ clusters can deviate substantially from spherical symmetry, with cold fronts resulting from gas sloshing extending out to at least half the cluster virial radius, far further out than previously realized. The Physics of Galaxy Cluster Outskirts 43

5.2 What are the most promising areas for future research in the physics of cluster outskirts?

The advent of the high throughput, high spectral resolution X-ray microcalorimeter aboard ATHENA will herald a new era in cluster outskirts research. For the first time we will be able to characterize the ICM velocity field, and be able to assess the level of turbulence and coherent bulk motion of the ICM, allowing us to accurately probe and map the level of non-thermal pressure support in cluster outskirts. The high effective areas of new X-ray observatories will dramatically improve our ability to measure metal abundances in cluster outskirts, possibly allowing the resolved contributions from single elements, to study the enrichment processes that are occurring in the accreting regions. We will be able to explore the preferred routes of mass accretion, and explore the connections between clusters and the surrounding cosmic web. Studies of cluster outskirts will be able to extend out to very large samples, reaching down to very small groups and out to high redshift. We will then be able to explore how the cluster outskirts has evolved with cosmic time, and for the cluster population as a whole. Increasingly powerful missions studying the SZ effect will be able to attack the problem of observing the accretion shocks around clusters.

Acknowledgements S.W. was supported by an appointment to the NASA Postdoctoral Program at the Goddard Space Flight Center, administered by the Universities Space Re- search Association through a contract with NASA. A.S. gratefully acknowledges support by the Women In Science Excel (WISE) programme of the Netherlands Organisation for Scientific Research (NWO). D.N. acknowledges Yale University for granting a triennial leave and the Max-Planck-Institut fürAstrophysik for hospitality when this work was completed. N.O. was supported by the Funds for the Development of HumanResources in Science and Technology under MEXT, Japan and Core Research for Energetic Universe in Hiroshima University (the MEXT program for promoting the enhancement of research universities, Japan). T.M. is supported for scientific activities by ESO’s Directorate for Science. H.A. acknowledges the support of NWO via a Veni grant. S.E. acknowledges financial contribution from the contracts NARO15 ASI-INAF I/037/12/0, ASI 2015-046-R.0 and ASI-INAF n.2017-14-H.0. SRON is supported financially by NWO, the Netherlands Organization for Scientific Research.

References

K.N. Abazajian, P. Adshead, Z. Ahmed, S.W. Allen, D. Alonso, K.S. Arnold, C. Bacci- galupi, J.G. Bartlett, N. Battaglia, B.A. Benson, C.A. Bischoﬀ, J. Borrill, V. Buza, E. Calabrese, R. Caldwell, J.E. Carlstrom, C.L. Chang, T.M. Crawford, F.-Y. Cyr-Racine, F. De Bernardis, T. de Haan, S. di Serego Alighieri, J. Dunkley, C. Dvorkin, J. Errard, G. Fabbian, S. Feeney, S. Ferraro, J.P. Filippini, R. Flauger, G.M. Fuller, V. Glusce- vic, D. Green, D. Grin, E. Grohs, J.W. Henning, J.C. Hill, R. Hlozek, G. Holder, W. Holzapfel, W. Hu, K.M. Huﬀenberger, R. Keskitalo, L. Knox, A. Kosowsky, J. Kovac, E.D. Kovetz, C.-L. Kuo, A. Kusaka, M. Le Jeune, A.T. Lee, M. Lilley, M. Loverde, M.S. Madhavacheril, A. Mantz, D.J.E. Marsh, J. McMahon, P.D. Meerburg, J. Meyers, A.D. Miller, J.B. Munoz, H.N. Nguyen, M.D. Niemack, M. Peloso, J. Peloton, L. Pogosian, C. Pryke, M. Raveri, C.L. Reichardt, G. Rocha, A. Rotti, E. Schaan, M.M. Schmittfull, D. Scott, N. Sehgal, S. Shandera, B.D. Sherwin, T.L. Smith, L. Sorbo, G.D. Starkman, 44 Stephen Walker et al.

K.T. Story, A. van Engelen, J.D. Vieira, S. Watson, N. Whitehorn, W.L. Kimmy Wu, CMB-S4 Science Book, First Edition. ArXiv e-prints (2016) R. Adam, A. Adane, P.A.R. Ade, P. André,A. Andrianasolo, H. Aussel, A. Beelen, A. Benoˆıt, A. Bideaud, N. Billot, O. Bourrion, A. Bracco, M. Calvo, A. Catalano, G. Coiffard, B. Comis, M. De Petris, F.-X. Désert,S. Doyle, E.F.C. Driessen, R. Evans, J. Goupy, C. Kramer, G. Lagache, S. Leclercq, J.-P. Leggeri, J.-F. Lestrade, J.F. Mac´ıas-Pérez, P. Mauskopf, F. Mayet, A. Maury, A. Monfardini, S. Navarro, E. Pascale, L. Perotto, G. Pisano, N. Ponthieu, V. Revéret,A. Rigby, A. Ritacco, C. Romero, H. Roussel, F. Ruppin, K. Schuster, A. Sievers, S. Triqueneaux, C. Tucker, R. Zylka, The NIKA2 large- field-of-view millimetre continuum camera for the 30 m IRAM telescope. A&A 609, 115 (2018). doi:10.1051/0004-6361/201731503 S. Adhikari, N. Dalal, R.T. Chamberlain, Splashback in accreting dark matter halos. ArXiv e-prints (2014) H. Aihara, N. Arimoto, R. Armstrong, S. Arnouts, N.A. Bahcall, S. Bickerton, J. Bosch, K. Bundy, P.L. Capak, J.H.H. Chan, M. Chiba, J. Coupon, E. Egami, M. Enoki, F. Finet, H. Fujimori, S. Fujimoto, H. Furusawa, J. Furusawa, T. Goto, A. Goulding, J.P. Greco, J.E. Greene, J.E. Gunn, T. Hamana, Y. Harikane, Y. Hashimoto, T. Hattori, M. Hayashi, Y. Hayashi, K.G. He lminiak,R. Higuchi, C. Hikage, P.T.P. Ho, B.-C. Hsieh, K. Huang, S. Huang, H. Ikeda, M. Imanishi, A.K. Inoue, K. Iwasawa, I. Iwata, A.T. Jaelani, H.-Y. Jian, Y. Kamata, H. Karoji, N. Kashikawa, N. Katayama, S. Kawanomoto, I. Kayo, J. Koda, M. Koike, T. Kojima, Y. Komiyama, A. Konno, S. Koshida, Y. Koyama, H. Kusakabe, A. Leauthaud, C.-H. Lee, L. Lin, Y.-T. Lin, R.H. Lupton, R. Mandelbaum, Y. Matsuoka, E. Medezinski, S. Mineo, S. Miyama, H. Miyatake, S. Miyazaki, R. Momose, A. More, S. More, Y. Moritani, T.J. Moriya, T. Morokuma, S. Mukae, R. Murata, H. Murayama, T. Nagao, F. Nakata, M. Niida, H. Niikura, A.J. Nishizawa, Y. Obuchi, M. Oguri, Y. Oishi, N. Okabe, S. Okamoto, Y. Okura, Y. Ono, M. Onodera, M. Onoue, K. Osato, M. Ouchi, P.A. Price, T.-S. Pyo, M. Sako, M. Sawicki, T. Shibuya, K. Shimasaku, A. Shimono, M. Shirasaki, J.D. Silverman, M. Simet, J. Speagle, D.N. Spergel, M.A. Strauss, Y. Sugahara, N. Sugiyama, Y. Suto, S.H. Suyu, N. Suzuki, P.J. Tait, M. Takada, T. Takata, N. Tamura, M.M. Tanaka, M. Tanaka, M. Tanaka, Y. Tanaka, T. Terai, Y. Terashima, Y. Toba, N. Tominaga, J. Toshikawa, E.L. Turner, T. Uchida, H. Uchiyama, K. Umetsu, F. Uraguchi, Y. Urata, T. Usuda, Y. Utsumi, S.-Y. Wang, W.-H. Wang, K.C. Wong, K. Yabe, Y. Yamada, H. Yamanoi, N. Yasuda, S. Yeh, A. Yonehara, S. Yuma, The Hyper Suprime-Cam SSP Survey: Overview and survey design. PASJ 70, 4 (2018). doi:10.1093/pasj/psx066 H. Akamatsu, A. Hoshino, Y. Ishisaki, T. Ohashi, K. Sato, Y. Takei, N. Ota, X-Ray Study of the Outer Region of Abell 2142 with Suzaku. PASJ 63, 1019–1033 (2011). doi:10.1093/pasj/63.sp3.S1019 H. Akamatsu, J. de Plaa, J. Kaastra, Y. Ishisaki, T. Ohashi, M. Kawaharada, K. Nakazawa, Properties of the Intracluster Medium of Abell 3667 Observed with Suzaku XIS. PASJ 64, 49 (2012a). doi:10.1093/pasj/64.3.49 H. Akamatsu, M. Takizawa, K. Nakazawa, Y. Fukazawa, Y. Ishisaki, T. Ohashi, X-Ray View of the Shock Front in the Merging Cluster Abell 3376 with Suzaku. PASJ 64, 67 (2012b). doi:10.1093/pasj/64.4.67 H. Akamatsu, S. Inoue, T. Sato, K. Matsusita, Y. Ishisaki, C.L. Sarazin, Suzaku X-Ray Observations of the Accreting NGC 4839 Group of Galaxies and a Radio Relic in the Coma Cluster. PASJ 65, 89 (2013a). doi:10.1093/pasj/65.4.89 H. Akamatsu, S. Inoue, T. Sato, K. Matsusita, Y. Ishisaki, C.L. Sarazin, Suzaku X-Ray Observations of the Accreting NGC 4839 Group of Galaxies and a Radio Relic in the Coma Cluster. PASJ 65, 89 (2013b). doi:10.1093/pasj/65.4.89 H. Akamatsu, M. Mizuno, N. Ota, Y.-Y. Zhang, R.J. van Weeren, H. Kawahara, Y. Fukazawa, J.S. Kaastra, M. Kawaharada, K. Nakazawa, T. Ohashi, H.J.A. Röttgering, M. Takizawa, J. Vink, F. Zandanel, Suzaku observations of the merging galaxy cluster Abell 2255: The northeast radio relic. A&A 600, 100 (2017). doi:10.1051/0004- 6361/201628400 S.W. Allen, D.A. Rapetti, R.W. Schmidt, H. Ebeling, R.G. Morris, A.C. Fabian, Improved constraints on dark energy from Chandra X-ray observations of the largest relaxed galaxy clusters. MNRAS 383, 879–896 (2008). doi:10.1111/j.1365-2966.2007.12610.x The Physics of Galaxy Cluster Outskirts 45

A.J. Anderson, P.A.R. Ade, Z. Ahmed, J.E. Austermann, J.S. Avva, P.S. Barry, R.B. Thakur, A.N. Bender, B.A. Benson, L.E. Bleem, K. Byrum, J.E. Carlstrom, F.W. Carter, T. Cecil, C.L. Chang, H.M. Cho, J.F. Cliche, T.M. Crawford, A. Cukierman, E.V. Denison, T. de Haan, J. Ding, M.A. Dobbs, D. Dutcher, W. Everett, A. Foster, R.N. Gannon, A. Gilbert, J.C. Groh, N.W. Halverson, A.H. Harke-Hosemann, N.L. Har- rington, J.W. Henning, G.C. Hilton, G.P. Holder, W.L. Holzapfel, N. Huang, K.D. Irwin, O.B. Jeong, M. Jonas, T. Khaire, L. Knox, A.M. Kofman, M. Korman, D. Kubik, S. Kuhlmann, N. Kuklev, C.L. Kuo, A.T. Lee, E.M. Leitch, A.E. Lowitz, S.S. Meyer, D. Michalik, J. Montgomery, A. Nadolski, T. Natoli, H. Nguyen, G.I. Noble, V. Novosad, S. Padin, Z. Pan, J. Pearson, C.M. Posada, A. Rahlin, C.L. Reichardt, J.E. Ruhl, L.J. Saunders, J.T. Sayre, I. Shirley, E. Shirokoff, G. Smecher, J.A. Sobrin, A.A. Stark, K.T. Story, A. Suzuki, Q.Y. Tang, K.L. Thompson, C. Tucker, L.R. Vale, K. Vanderlinde, J.D. Vieira, G. Wang, N. Whitehorn, V. Yefremenko, K.W. Yoon, M.R. Young, SPT- 3G: A Multichroic Receiver for the South Pole Telescope. Journal of Low Temperature Physics (2018). doi:10.1007/s10909-018-2007-z F. AndradeSantos, P.E.J. Nulsen, R.P. Kraft, W.R. Forman, C. Jones, E. Churazov, A. Vikhlinin, Dark Matter Subhalos and the X-Ray Morphology of the Coma Cluster. ApJ 766, 107 (2013). doi:10.1088/0004-637X/766/2/107 M. Arnaud, G.W. Pratt, R. Piffaretti, H. Böhringer, J.H. Croston, E. Pointecouteau, The universal galaxy cluster pressure profile from a representative sample of nearby systems (REXCESS) and the YSZ -M500 relation. A&A 517, 92 (2010). doi:10.1051/0004- 6361/200913416 J.E. Austermann, J.A. Beall, S.A. Bryan, B. Dober, J. Gao, G. Hilton, J. Hubmayr, P. Mauskopf, C.M. McKenney, S.M. Simon, J.N. Ullom, M.R. Vissers, G.W. Wilson, Millimeter-Wave Polarimeters Using Kinetic Inductance Detectors for TolTEC and Be- yond. Journal of Low Temperature Physics (2018). doi:10.1007/s10909-018-1949-5 C. Avestruz, D. Nagai, E.T. Lau, Stirred, not Clumped: Evolution of Temperature Profiles in the Outskirts of Galaxy Clusters. ArXiv e-prints (2016) C. Avestruz, D. Nagai, E.T. Lau, K. Nelson, Non-equilibrium Electrons in the Outskirts of Galaxy Clusters. ApJ 808, 176 (2015). doi:10.1088/0004-637X/808/2/176 D. Barret, T. Lam Trong, J.-W. den Herder, L. Piro, M. Cappi, J. Houvelin, R. Kelley, J.M. Mas-Hesse, K. Mitsuda, S. Paltani, G. Rauw, A. Rozanska, J. Wilms, S. Bandler, M. Bar- bera, X. Barcons, E. Bozzo, M.T. Ceballos, I. Charles, E. Costantini, A. Decourchelle, R. den Hartog, L. Duband, J.-M. Duval, F. Fiore, F. Gatti, A. Goldwurm, B. Jackson, P. Jonker, C. Kilbourne, C. Macculi, M. Mendez, S. Molendi, P. Orleanski, F. Pajot, E. Pointecouteau, F. Porter, G.W. Pratt, D. Prêle,L. Ravera, K. Sato, J. Schaye, K. Shinozaki, T. Thibert, L. Valenziano, V. Valette, J. Vink, N. Webb, M. Wise, N. Ya- masaki, F. Douchin, J.-M. Mesnager, B. Pontet, A. Pradines, G. Branduardi-Raymont, E. Bulbul, M. Dadina, S. Ettori, A. Finoguenov, Y. Fukazawa, A. Janiuk, J. Kaastra, P. Mazzotta, J. Miller, G. Miniutti, Y. Naze, F. Nicastro, S. Scioritino, A. Simonescu, J.M. Torrejon, B. Frezouls, H. Geoffray, P. Peille, C. Aicardi, J. André,C. Daniel, A. Clénet, C. Etcheverry, E. Gloaguen, G. Hervet, A. Jolly, A. Ledot, I. Paillet, R. Schmisser, B. Vella, J.-C. Damery, K. Boyce, M. Dipirro, S. Lotti, D. Schwander, S. Smith, B.-J. Van Leeuwen, H. van Weers, N. Clerc, B. Cobo, T. Dauser, C. Kirsch, E. Cucchetti, M. Eckart, P. Ferrando, L. Natalucci, The ATHENA X-ray Integral Field Unit (X-IFU), in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 10699, 2018, p. 106991. doi:10.1117/12.2312409 M. Bartelmann, P. Schneider, Weak gravitational lensing. Phys. Rep. 340, 291–472 (2001). doi:10.1016/S0370-1573(00)00082-X K. Basu, M. Sommer, J. Erler, D. Eckert, F. Vazza, B. Magnelli, F. Bertoldi, P. Tozzi, ALMA-SZ Detection of a Galaxy Cluster Merger Shock at Half the Age of the Universe. ApJ 829, 23 (2016). doi:10.3847/2041-8205/829/2/L23 N. Battaglia, J.R. Bond, C. Pfrommer, J.L. Sievers, On the Cluster Physics of Sunyaev- Zel’dovich and X-ray Surveys III: Measurement Biases and Cosmological Evolution of Gas and Stellar Mass Fractions. ArXiv e-prints (2012) M.W. Bautz, E.D. Miller, J.S. Sanders, K.A. Arnaud, R.F. Mushotzky, F.S. Porter, K. Hayashida, J.P. Henry, J.P. Hughes, M. Kawaharada, K. Makashima, M. Sato, T. 46 Stephen Walker et al.

Tamura, Suzaku Observations of Abell 1795: Cluster Emission to r200. PASJ 61, 1117 (2009) E. Baxter, C. Chang, B. Jain, S. Adhikari, N. Dalal, A. Kravtsov, S. More, E. Rozo, E. Rykoff, R.K. Sheth, The Halo Boundary of Galaxy Clusters in the SDSS. ApJ 841, 18 (2017). doi:10.3847/1538-4357/aa6ff0 A.N. Bender, P.A.R. Ade, A.J. Anderson, J. Avva, Z. Ahmed, K. Arnold, J.E. Austermann, R. Basu Thakur, B.A. Benson, L.E. Bleem, K. Byrum, J.E. Carlstrom, F.W. Carter, C.L. Chang, H.M. Cho, J.F. Cliche, T.M. Crawford, A. Cukierman, D.A. Czaplewski, J. Ding, R. Divan, T. de Haan, M.A. Dobbs, D. Dutcher, W. Everett, A. Gilbert, J.C. Groh, R. Guyser, N.W. Halverson, A. Harke-Hosemann, N.L. Harrington, K. Hattori, J.W. Henning, G.C. Hilton, W.L. Holzapfel, N. Huang, K.D. Irwin, O. Jeong, T. Khaire, M. Korman, D. Kubik, C.L. Kuo, A.T. Lee, E.M. Leitch, S. Lendinez, S.S. Meyer, C.S. Miller, J. Montgomery, A. Nadolski, T. Natoli, H. Nguyen, V. Novosad, S. Padin, Z. Pan, J. Pearson, C.M. Posada, A. Rahlin, C.L. Reichardt, J.E. Ruhl, B.R. Saliwanchik, J.T. Sayre, J.A. Shariff, I. Shirley, E. Shirokoff, G. Smecher, J. Sobrin, L. Stan, A.A. Stark, K. Story, A. Suzuki, Q.Y. Tang, K.L. Thompson, C. Tucker, K. Vanderlinde, J.D. Vieira, G. Wang, N. Whitehorn, V. Yefremenko, K.W. Yoon, Integrated performance of a frequency domain multiplexing readout in the SPT-3G receiver, in Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VIII. Proc. SPIE, vol. 9914, 2016, p. 99141. doi:10.1117/12.2232146 B.A. Benson, P.A.R. Ade, Z. Ahmed, S.W. Allen, K. Arnold, J.E. Austermann, A.N. Ben- der, L.E. Bleem, J.E. Carlstrom, C.L. Chang, H.M. Cho, J.F. Cliche, T.M. Crawford, A. Cukierman, T. de Haan, M.A. Dobbs, D. Dutcher, W. Everett, A. Gilbert, N.W. Halverson, D. Hanson, N.L. Harrington, K. Hattori, J.W. Henning, G.C. Hilton, G.P. Holder, W.L. Holzapfel, K.D. Irwin, R. Keisler, L. Knox, D. Kubik, C.L. Kuo, A.T. Lee, E.M. Leitch, D. Li, M. McDonald, S.S. Meyer, J. Montgomery, M. Myers, T. Natoli, H. Nguyen, V. Novosad, S. Padin, Z. Pan, J. Pearson, C. Reichardt, J.E. Ruhl, B.R. Saliwanchik, G. Simard, G. Smecher, J.T. Sayre, E. Shirokoff, A.A. Stark, K. Story, A. Suzuki, K.L. Thompson, C. Tucker, K. Vanderlinde, J.D. Vieira, A. Vikhlinin, G. Wang, V. Yefremenko, K.W. Yoon, SPT-3G: a next-generation cosmic microwave background polarization experiment on the South Pole telescope, in Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VII. Proc. SPIE, vol. 9153, 2014, p. 91531. doi:10.1117/12.2057305 F. Bertoldi, The Atacama Large Aperture Submm/mm Telescope (AtLAST) Project, in Atacama Large-Aperture Submm/mm Telescope (AtLAST), 2018, p. 3. doi:10.5281/zenodo.1158842 M. Béthermin,H.-Y. Wu, G. Lagache, I. Davidzon, N. Ponthieu, M. Cousin, L. Wang, O. Doré,E. Daddi, A. Lapi, The impact of clustering and angular resolution on far- infrared and millimeter continuum observations. A&A 607, 89 (2017). doi:10.1051/0004- 6361/201730866 M. Bonamente, D. Landry, B. Maughan, P. Giles, M. Joy, J. Nevalainen, Chandra X- ray observations of Abell 1835 to the virial radius. MNRAS 428, 2812–2823 (2013a). doi:10.1093/mnras/sts202 M. Bonamente, D. Landry, B. Maughan, P. Giles, M. Joy, J. Nevalainen, Chandra X- ray observations of Abell 1835 to the virial radius. MNRAS 428, 2812–2823 (2013b). doi:10.1093/mnras/sts202 S. Bryan, J. Austermann, D. Ferrusca, P. Mauskopf, J. McMahon, A. Montana, S. Simon, G. Novak, D. Sanchez-Arguelles, G. Wilson, for the TolTEC Team, Optical Design of the TolTEC Millimeter-wave Camera. ArXiv e-prints (2018) E. Bulbul, M. Markevitch, A. Foster, E. Miller, M. Bautz, M. Loewenstein, S.W. Randall, R.K. Smith, Searching for the 3.5 keV Line in the Stacked Suzaku Observations of Galaxy Clusters. ApJ 831, 55 (2016). doi:10.3847/0004-637X/831/1/55 P. Busch, S.D.M. White, Assembly bias and splashback in galaxy clusters. MNRAS 470, 4767–4781 (2017). doi:10.1093/mnras/stx1584 J.E. Carlstrom, G.P. Holder, E.D. Reese, Cosmology with the Sunyaev-Zel’dovich Effect. ARA&A 40, 643–680 (2002) A. Catalano, R. Adam, P.A.R. Ade, P. André,H. Aussel, A. Beelen, A. Benoˆıt,A. Bideaud, N. Billot, O. Bourrion, M. Calvo, B. Comis, M. De Petris, F.-X. Désert,S. Doyle, The Physics of Galaxy Cluster Outskirts 47

E.F.C. Driessen, J. Goupy, C. Kramer, G. Lagache, S. Leclercq, J.-F. Lestrade, J.F. Mac´ıas-Pérez,P. Mauskopf, F. Mayet, A. Monfardini, E. Pascale, L. Perotto, G. Pisano, N. Ponthieu, V. Revéret,A. Ritacco, C. Romero, H. Roussel, F. Ruppin, K. Schuster, A. Sievers, S. Triqueneaux, C. Tucker, R. Zylka, E. Barria, G. Bres, P. Camus, P. Chanthib, G. Donnier-Valentin, O. Exshaw, G. Garde, A. Gerardin, J.-P. Leggeri, F. Levy-Bertrand, C. Guttin, C. Hoarau, M. Grollier, J.-L. Mocellin, G. Pont, H. Rodenas, O. Tissot, G. Galvez, D. John, H. Ungerechts, S. Sanchez, P. Mellado, M. Munoz, F. Pierfederici, J. Penalver, S. Navarro, G. Bosson, J.-L. Bouly, J. Bouvier, C. Geraci, C. Li, J. Menu, N. Ponchant, S. Roni, S. Roudier, J.P. Scordillis, D. Tourres, C. Vescovi, A. Barbier, D. Billon-Pierron, A. Adane, A. Andrianasolo, A. Bracco, G. Coiffard, R. Evans, A. Maury, A. Rigby, The NIKA2 Instrument at 30-m IRAM Telescope: Performance and Results. Journal of Low Temperature Physics (2018). doi:10.1007/s10909-018-1884-5 C. Chang, E. Baxter, B. Jain, C. Sánchez, S. Adhikari, T.N. Varga, Y. Fang, E. Rozo, E.S. Rykoff, A. Kravtsov, D. Gruen, E.M. Huff, M. Jarvis, A.G. Kim, J. Prat, N. MacCrann, T. McClintock, A. Palmese, D. Rapetti, R.P. Rollins, S. Samuroff, E. Sheldon, M.A. Troxel, R.H. Wechsler, Y. Zhang, J. Zuntz, T.M.C. Abbott, F.B. Abdalla, S. Allam, J. Annis, K. Bechtol, A. Benoit-Lévy,G.M. Bernstein, D. Brooks, E. Buckley-Geer, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, C.B. D’Andrea, L.N. da, C. Davis, S. Desai, H.T. Diehl, J.P. Dietrich, A. Drlica-Wagner, T.F. Eifler, B. Flaugher, P. Fosalba, J. Frieman, J. Garc´ıa-Bellido,E. Gaztanaga, D.W. Gerdes, R.A. Gruendl, J. Gschwend, G. Gutierrez, K. Honscheid, D.J. James, T. Jeltema, E. Krause, K. Kuehn, O. Lahav, M. Lima, M. March, J.L. Marshall, P. Martini, P. Melchior, F. Menanteau, R. Miquel, J.J. Mohr, B. Nord, R.L.C. Ogando, A.A. Plazas, E. Sanchez, V. Scarpine, R. Schindler, M. Schubnell, I. Sevilla-Noarbe, R.C. Smith, M. Smith, M. Soares-Santos, F. Sobreira, E. Suchyta, M.E.C. Swanson, G. Tarle, J. Weller, The Splashback Feature around DES Galaxy Clusters: Galaxy Density and Weak Lensing Profiles. ArXiv e-prints (2017) J. Chluba, D. Nagai, S. Sazonov, K. Nelson, A fast and accurate method for computing the Sunyaev-Zel’dovich signal of hot galaxy clusters. MNRAS 426, 510–530 (2012). doi:10.1111/j.1365-2966.2012.21741.x A.T. Crites, J.J. Bock, C.M. Bradford, T.C. Chang, A.R. Cooray, L. Duband, Y. Gong, S. Hailey-Dunsheath, J. Hunacek, P.M. Koch, C.T. Li, R.C. O’Brient, T. Prouve, E. Shirokoff, M.B. Silva, Z. Staniszewski, B. Uzgil, M. Zemcov, The TIME-Pilot intensity mapping experiment, in Millimeter, Submillimeter, and Far-Infrared Detec- tors and Instrumentation for Astronomy VII. Proc. SPIE, vol. 9153, 2014, p. 91531. doi:10.1117/12.2057207 K.T. Crowley, J.E. Austermann, S.K. Choi, S.M. Duff, P.A. Gallardo, S.-P.P. Ho, J. Hub- mayr, B.J. Koopman, F. Nati, M.D. Niemack, M. Salatino, S.M. Simon, S.T. Staggs, J.R. Stevens, J.N. Ullom, E.M. Vavagiakis, E.J. Wollack, Advanced ACTPol TES Device Parameters and Noise Performance in Fielded Arrays. ArXiv e-prints (2018) Dark Energy Survey Collaboration, T. Abbott, F.B. Abdalla, J. Aleksić,S. Allam, A. Amara, D. Bacon, E. Balbinot, M. Banerji, K. Bechtol, A. Benoit-Lévy,G.M. Bernstein, E. Bertin, J. Blazek, C. Bonnett, S. Bridle, D. Brooks, R.J. Brunner, E. Buckley-Geer, D.L. Burke, G.B. Caminha, D. Capozzi, J. Carlsen, A. Carnero-Rosell, M. Carollo, M. Carrasco-Kind, J. Carretero, F.J. Castander, L. Clerkin, T. Collett, C. Conselice, M. Crocce, C.E. Cunha, C.B. D’Andrea, L.N. da Costa, T.M. Davis, S. Desai, H.T. Diehl, J.P. Dietrich, S. Dodelson, P. Doel, A. Drlica-Wagner, J. Estrada, J. Etherington, A.E. Evrard, J. Fabbri, D.A. Finley, B. Flaugher, R.J. Foley, P. Fosalba, J. Frieman, J. Garc´ıa-Bellido,E. Gaztanaga, D.W. Gerdes, T. Giannantonio, D.A. Goldstein, D. Gruen, R.A. Gruendl, P. Guarnieri, G. Gutierrez, W. Hartley, K. Honscheid, B. Jain, D.J. James, T. Jeltema, S. Jouvel, R. Kessler, A. King, D. Kirk, R. Kron, K. Kuehn, N. Kuropatkin, O. Lahav, T.S. Li, M. Lima, H. Lin, M.A.G. Maia, M. Makler, M. Manera, C. Maraston, J.L. Marshall, P. Martini, R.G. McMahon, P. Melchior, A. Merson, C.J. Miller, R. Miquel, J.J. Mohr, X. Morice-Atkinson, K. Naidoo, E. Neilsen, R.C. Nichol, B. Nord, R. Ogando, F. Ostrovski, A. Palmese, A. Papadopoulos, H.V. Peiris, J. Peoples, W.J. Percival, A.A. Plazas, S.L. Reed, A. Refregier, A.K. Romer, A. Roodman, A. Ross, E. Rozo, E.S. Rykoff, I. Sadeh, M. Sako, C. Sánchez, E. Sanchez, B. Santiago, V. Scarpine, M. Schubnell, I. Sevilla-Noarbe, E. Sheldon, M. Smith, R.C. Smith, M. Soares- Santos, F. Sobreira, M. Soumagnac, E. Suchyta, M. Sullivan, M. Swanson, G. Tarle, J. 48 Stephen Walker et al.

Thaler, D. Thomas, R.C. Thomas, D. Tucker, J.D. Vieira, V. Vikram, A.R. Walker, R.H. Wechsler, J. Weller, W. Wester, L. Whiteway, H. Wilcox, B. Yanny, Y. Zhang, J. Zuntz, The Dark Energy Survey: more than dark energy - an overview. MNRAS 460, 1270–1299 (2016). doi:10.1093/mnras/stw641 A. de Graaff, Y.-C. Cai, C. Heymans, J.A. Peacock, Missing baryons in the cosmic web revealed by the Sunyaev-Zel’dovich effect. ArXiv e-prints (2017) S.R. Dicker, P.A.R. Ade, J. Aguirre, J.A. Brevik, H.M. Cho, R. Datta, M.J. Devlin, B. Dober, D. Egan, J. Ford, P. Ford, G. Hilton, K.D. Irwin, B.S. Mason, P. Marganian, M. Mello, J.J. McMahon, T. Mroczkowski, M. Rosenman, C. Tucker, L. Vale, S. White, M. Whitehead, A.H. Young, MUSTANG 2: A Large Focal Plane Array for the 100 m Green Bank Telescope. Journal of Low Temperature Physics 176, 808–814 (2014). doi:10.1007/s10909-013-1070-8 B. Diemer, A.V. Kravtsov, Dependence of the Outer Density Profiles of Halos on Their Mass Accretion Rate. ApJ 789, 1 (2014). doi:10.1088/0004-637X/789/1/1 B. Diemer, P. Mansfield, A.V. Kravtsov, S. More, The Splashback Radius of Halos from Particle Dynamics. II. Dependence on Mass, Accretion Rate, Redshift, and Cosmology. ApJ 843, 140 (2017). doi:10.3847/1538-4357/aa79ab E.M. Douglass, E.L. Blanton, S.W. Randall, T.E. Clarke, L.O.V. Edwards, Z. Sabry, J.A. ZuHone, The Megaparsec-scale Gas-sloshing Spiral in the Remnant Cool Core Cluster Abell 1763. ApJ 868, 121 (2018). doi:10.3847/1538-4357/aae9e7 D. Eckert, F. Vazza, S. Ettori, S. Molendi, D. Nagai, E.T. Lau, M. Roncarelli, M. Rossetti, S.L. Snowden, F. Gastaldello, The gas distribution in the outer regions of galaxy clusters. A&A 541, 57 (2012). doi:10.1051/0004-6361/201118281 D. Eckert, S. Molendi, F. Vazza, S. Ettori, S. Paltani, The X-ray/SZ view of the virial region. I. Thermodynamic properties. A&A 551, 22 (2013a). doi:10.1051/0004-6361/201220402 D. Eckert, S. Ettori, S. Molendi, F. Vazza, S. Paltani, The X-ray/SZ view of the virial region. II. Gas mass fraction. A&A 551, 23 (2013b). doi:10.1051/0004-6361/201220403 D. Eckert, M. Roncarelli, S. Ettori, S. Molendi, F. Vazza, F. Gastaldello, M. Rossetti, Gas clumping in galaxy clusters. MNRAS 447, 2198–2208 (2015). doi:10.1093/mnras/stu2590 D. Eckert, M. Gaspari, M.S. Owers, E. Roediger, S. Molendi, F. Gastaldello, S. Paltani, S. Ettori, T. Venturi, M. Rossetti, L. Rudnick, Deep Chandra observations of the stripped galaxy group falling into Abell 2142. A&A 605, 25 (2017). doi:10.1051/0004- 6361/201730555 D. Eckert, V. Ghirardini, S. Ettori, E. Rasia, V. Biffi, E. Pointecouteau, M. Rossetti, S. Molendi, F. Vazza, F. Gastaldello, M. Gaspari, S. De Grandi, S. Ghizzardi, H. Bourdin, C. Tchernin, M. Roncarelli, Non-thermal pressure support in X-COP galaxy clusters. ArXiv e-prints (2018) D.J. Eisenstein, D.H. Weinberg, E. Agol, H. Aihara, C. Allende Prieto, S.F. Anderson, J.A. Arns, E.´ Aubourg, S. Bailey, E. Balbinot, et al., SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way, and Extra-Solar Planetary Systems. AJ 142, 72 (2011). doi:10.1088/0004-6256/142/3/72 J. Erler, K. Basu, J. Chluba, F. Bertoldi, Planck’s view on the spectrum of the Sunyaev- Zeldovich effect. MNRAS 476, 3360–3381 (2018). doi:10.1093/mnras/sty327 S. Ettori, A.C. Fabian, ROSAT PSPC observations of 36 high-luminosity clusters of galaxies: constraints on the gas fraction. MNRAS 305, 834–848 (1999). doi:10.1046/j.1365- 8711.1999.02460.x S. Ettori, A. Donnarumma, E. Pointecouteau, T.H. Reiprich, S. Giodini, L. Lovisari, R.W. Schmidt, Mass Profiles of Galaxy Clusters from X-ray Analysis. Space Sci. Rev. 177, 119–154 (2013a). doi:10.1007/s11214-013-9976-7 S. Ettori, G.W. Pratt, J. de Plaa, D. Eckert, J. Nevalainen, E.S. Battistelli, S. Borgani, J.H. Croston, A. Finoguenov, J. Kaastra, M. Gaspari, F. Gastaldello, M. Gitti, S. Molendi, E. Pointecouteau, T.J. Ponman, T.H. Reiprich, M. Roncarelli, M. Rossetti, J.S. Sanders, M. Sun, G. Trinchieri, F. Vazza, M. Arnaud, H. Böringher,F. Brighenti, H. Dahle, S. De Grandi, J.J. Mohr, A. Moretti, S. Schindler, The Hot and Energetic Universe: The astrophysics of galaxy groups and clusters. ArXiv e-prints (2013b) S. Ettori, V. Ghirardini, D. Eckert, E. Pointecouteau, F. Gastaldello, M. Sereno, M. Gaspari, S. Ghizzardi, M. Roncarelli, M. Rossetti, Hydrostatic mass profiles in X-COP galaxy The Physics of Galaxy Cluster Outskirts 49

clusters. ArXiv e-prints (2018) A.E. Evrard, The intracluster gas fraction in X-ray clusters - Constraints on the clustered mass density. MNRAS 292, 289 (1997). doi:10.1093/mnras/292.2.289 Y. Fujita, N. Tawa, K. Hayashida, M. Takizawa, H. Matsumoto, N. Okabe, T.H. Reiprich, High Metallicity of the X-Ray Gas Up to the Virial Radius of a Binary Cluster of Galaxies: Evidence of Galactic Superwinds at High-Redshift. PASJ 60, 343 (2008) R. Fusco-Femiano, A. Lapi, Entropy Flattening, Gas Clumping, and Turbulence in Galaxy Clusters. ApJ 783, 76 (2014). doi:10.1088/0004-637X/783/2/76 J.A. Gaskin, R. Allured, S.R. Bandler, S. Basso, M.W. Bautz, M.F. Baysinger, M.P. Biskach, T.M. Boswell, P.D. Capizzo, K.-W. Chan, M.M. Civitani, L.M. Cohen, V. Cotroneo, J.M. Davis, C.T. DeRoo, M.J. DiPirro, A. Dominguez, L.L. Fabisinski, A.D. Falcone, E. Figueroa-Feliciano, J.C. Garcia, K.E. Gelmis, R.K. Heilmann, R.C. Hop- kins, T. Jackson, K. Kilaru, R.P. Kraft, T. Liu, R.S. McClelland, R.L. McEntaffer, K.S. McCarley, J.A. Mulqueen, F. Ozel,¨ G. Pareschi, P.B. Reid, R.E. Riveros, M.A. Rodriguez, J.W. Rowe, T.T. Saha, M.L. Schattenburg, A.R. Schnell, D.A. Schwartz, P.M. Solly, R.M. Suggs, S.G. Sutherlin, D.A. Swartz, S. Trolier-McKinstry, J.H. Tutt, A. Vikhlinin, J. Walker, W. Yoon, W.W. Zhang, Lynx Mission concept status, in Soci- ety of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 10397, 2017, p. 103970. doi:10.1117/12.2273911 M.R. George, A.C. Fabian, J.S. Sanders, A.J. Young, H.R. Russell, X-ray observations of the galaxy cluster PKS0745-191: to the virial radius, and beyond. MNRAS 395, 657–666 (2009). doi:10.1111/j.1365-2966.2009.14547.x V. Ghirardini, S. Ettori, S. Amodeo, R. Capasso, M. Sereno, On the evolution of the entropy and pressure profiles in X-ray luminous galaxy clusters at z greater than 0.4. A&A 604, 100 (2017). doi:10.1051/0004-6361/201630209 V. Ghirardini, D. Eckert, S. Ettori, E. Pointecouteau, S. Molendi, M. Gaspari, M. Rossetti, S. De Grandi, M. Roncarelli, H. Bourdin, P. Mazzotta, E. Rasia, F. Vazza, The universal thermodynamic properties of the intracluster medium over two decades in radius in the X-COP sample. ArXiv e-prints (2018a) V. Ghirardini, S. Ettori, D. Eckert, S. Molendi, F. Gastaldello, E. Pointecouteau, G. Hurier, H. Bourdin, The XMM Cluster Outskirts Project (X-COP): Thermodynamic properties of the intracluster medium out to R200 in Abell 2319. A&A 614, 7 (2018b). doi:10.1051/0004-6361/201731748 M. Hilton, M. Hasselfield, C. Sifón,N. Battaglia, S. Aiola, V. Bharadwaj, J.R. Bond, S.K. Choi, D. Crichton, R. Datta, M.J. Devlin, J. Dunkley, R. Dünner,P.A. Gallardo, M. Gralla, A.D. Hincks, S.-P.P. Ho, J. Hubmayr, K.M. Huffenberger, J.P. Hughes, B.J. Koopman, A. Kosowsky, T. Louis, M.S. Madhavacheril, T.A. Marriage, L. Maurin, J. McMahon, H. Miyatake, K. Moodley, S. Næss, F. Nati, L. Newburgh, M.D. Niemack, M. Oguri, L.A. Page, B. Partridge, B.L. Schmitt, J. Sievers, D.N. Spergel, S.T. Staggs, H. Trac, A. van Engelen, E.M. Vavagiakis, E.J. Wollack, The Atacama Cosmology Tele- scope: The Two-season ACTPol Sunyaev-Zeldovich Effect Selected Cluster Catalog. ApJS 235, 20 (2018). doi:10.3847/1538-4365/aaa6cb H. Hoekstra, How well can we determine cluster mass profiles from weak lensing? MNRAS 339, 1155–1162 (2003). doi:10.1046/j.1365-8711.2003.06264.x A. Hoshino, J.P. Henry, K. Sato, H. Akamatsu, W. Yokota, S. Sasaki, Y. Ishisaki, T. Ohashi, M. Bautz, Y. Fukazawa, N. Kawano, A. Furuzawa, K. Hayashida, N. Tawa, J.P. Hughes, M. Kokubun, T. Tamura, X-Ray Temperature and Mass Measurements to the Virial Radius of Abell 1413 with Suzaku. PASJ 62, 371–389 (2010). doi:10.1093/pasj/62.2.371 P.J. Humphrey, D.A. Buote, F. Brighenti, H.M.L.G. Flohic, F. Gastaldello, W.G. Mathews, Tracing the Gas to the Virial Radius (R 100) in a Fossil Group. ApJ 748, 11 (2012). doi:10.1088/0004-637X/748/1/11 G. Hurier, R. Adam, U. Keshet, First detection of a virial shock with SZ data: implication for the mass accretion rate of Abell 2319. ArXiv e-prints (2017) Y. Ibaraki, N. Ota, H. Akamatsu, Y.-Y. Zhang, A. Finoguenov, Suzaku study of gas properties along filaments of A2744. A&A 562, 11 (2014). doi:10.1051/0004-6361/201322806 K. Ichikawa, K. Matsushita, N. Okabe, K. Sato, Y.-Y. Zhang, A. Finoguenov, Y. Fujita, Y. Fukazawa, M. Kawaharada, K. Nakazawa, T. Ohashi, N. Ota, M. Takizawa, T. Tamura, 50 Stephen Walker et al.

K. Umetsu, Suzaku Observations of the Outskirts of A1835: Deviation from Hydrostatic Equilibrium. ApJ 766, 90 (2013). doi:10.1088/0004-637X/766/2/90 Y. Ichinohe, N. Werner, A. Simionescu, S.W. Allen, R.E.A. Canning, S. Ehlert, F. Mernier, T. Takahashi, The growth of the galaxy cluster Abell 85: mergers, shocks, stripping and seeding of clumping. MNRAS 448, 2971–2986 (2015). doi:10.1093/mnras/stv217 N. Itoh, Y. Kohyama, S. Nozawa, Relativistic Corrections to the Sunyaev-Zeldovich Effect for Clusters of Galaxies. ApJ 502, 7–15 (1998). doi:10.1086/305876 M.J. Jee, A. Stroe, W. Dawson, D. Wittman, H. Hoekstra, M. Brüggen, H. Röttgering, D. Sobral, R.J. van Weeren, MC 2: Constraining the Dark Matter Distribution of the Violent Merging Galaxy Cluster CIZA J2242.8+5301 by Piercing through the Milky Way. ApJ 802, 46 (2015). doi:10.1088/0004-637X/802/1/46 M. Kawaharada, N. Okabe, K. Umetsu, M. Takizawa, K. Matsushita, Y. Fukazawa, T. Hamana, S. Miyazaki, K. Nakazawa, T. Ohashi, Suzaku Observation of A1689: Anisotropic Temperature and Entropy Distributions Associated with the Large-scale Structure. ApJ 714, 423–441 (2010). doi:10.1088/0004-637X/714/1/423 T. Kitayama, S. Ueda, S. Takakuwa, T. Tsutsumi, E. Komatsu, T. Akahori, D. Iono, T. Izumi, R. Kawabe, K. Kohno, H. Matsuo, N. Ota, Y. Suto, M. Takizawa, K. Yoshikawa, The Sunyaev-Zel’dovich effect at 5”: RX J1347.5-1145 imaged by ALMA. PASJ 68, 88 (2016). doi:10.1093/pasj/psw082 G. Lagache, Continuum and Line Emission Simulation of Star-Forming Galaxies and De- velopment of a New Sub-mm Inte, in Atacama Large-Aperture Submm/mm Telescope (AtLAST), 2018, p. 22. doi:10.5281/zenodo.1159045 A. Lapi, R. Fusco-Femiano, A. Cavaliere, Probing the astrophysics of cluster outskirts. A&A 516, 34 (2010). doi:10.1051/0004-6361/201014218 E.T. Lau, A.V. Kravtsov, D. Nagai, Residual Gas Motions in the Intracluster Medium and Bias in Hydrostatic Measurements of Mass Profiles of Clusters. ApJ 705, 1129–1138 (2009). doi:10.1088/0004-637X/705/2/1129 E.T. Lau, D. Nagai, C. Avestruz, K. Nelson, A. Vikhlinin, Mass Accretion and its Effects on the Self-similarity of Gas Profiles in the Outskirts of Galaxy Clusters. ApJ 806, 68 (2015). doi:10.1088/0004-637X/806/1/68 P. Mansfield, A.V. Kravtsov, B. Diemer, Splashback Shells of Cold Dark Matter Halos. ApJ 841, 34 (2017). doi:10.3847/1538-4357/aa7047 M. Markevitch, A. Vikhlinin, Shocks and cold fronts in galaxy clusters. Phys. Rep. 443, 1–53 (2007). doi:10.1016/j.physrep.2007.01.001 B. Mathiesen, A.E. Evrard, J.J. Mohr, The Effects of Clumping and Substructure on Intr- acluster Medium Mass Measurements. ApJ 520, 21–24 (1999) N. Meidinger, M. Barbera, V. Emberger, M. Fürmetz,M. Manhart, J. Müller-Seidlitz,K. Nandra, M. Plattner, A. Rau, W. Treberspurg, The Wide Field Imager instrument for Athena, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 10397, 2017, p. 103970. doi:10.1117/12.2271844 S.M. Molnar, N. Hearn, Z. Haiman, G. Bryan, A.E. Evrard, G. Lake, Accretion Shocks in Clusters of Galaxies and Their SZ Signature from Cosmological Simulations. ApJ 696, 1640–1656 (2009). doi:10.1088/0004-637X/696/2/1640 A. Morandi, W. Cui, Measuring the gas clumping in Abell 133. MNRAS 437, 1909–1917 (2014). doi:10.1093/mnras/stt2021 A. Morandi, D. Nagai, W. Cui, Non-parametric method for measuring gas inhomogeneities from X-ray observations of galaxy clusters. MNRAS 436, 1123–1129 (2013). doi:10.1093/mnras/stt1636 A. Morandi, M. Limousin, J. Sayers, S.R. Golwala, N.G. Czakon, E. Pierpaoli, E. Jullo, J. Richard, S. Ameglio, X-ray, lensing and Sunyaev-Zel’dovich triaxial analysis of Abell 1835 out to R200. MNRAS 425, 2069–2082 (2012). doi:10.1111/j.1365- 2966.2012.21196.x A. Morandi, M. Sun, W. Forman, C. Jones, The galaxy cluster outskirts probed by Chandra. MNRAS 450, 2261–2278 (2015). doi:10.1093/mnras/stv660 A. Morandi, M. Sun, J. Mulchaey, D. Nagai, M. Bonamente, Gas distribution and clumpiness in the galaxy group NGC 2563. MNRAS 469, 2423–2433 (2017). doi:10.1093/mnras/stx1031 The Physics of Galaxy Cluster Outskirts 51

S. More, B. Diemer, A.V. Kravtsov, The Splashback Radius as a Physical Halo Boundary and the Growth of Halo Mass. ApJ 810, 36 (2015). doi:10.1088/0004-637X/810/1/36 S. More, H. Miyatake, M. Takada, B. Diemer, A.V. Kravtsov, N.K. Dalal, A. More, R. Murata, R. Mandelbaum, E. Rozo, E.S. Rykoff, M. Oguri, D.N. Spergel, Detection of the Splashback Radius and Halo Assembly Bias of Massive Galaxy Clusters. ApJ 825, 39 (2016). doi:10.3847/0004-637X/825/1/39 T. Mroczkowski, O. Noroozian, AtLAST Instrumentation Considerations and Overview, in Atacama Large-Aperture Submm/mm Telescope (AtLAST), 2018, p. 26. doi:10.5281/zenodo.1159053 T. Mroczkowski, D. Nagai, K. Basu, J. Chluba, J. Sayers, R. Adam, E. Churazov, A. Crites, L. Di Mascolo, D. Eckert, J. Macias-Perez, F. Mayet, L. Perotto, E. Pointecouteau, C. Romero, F. Ruppin, E. Scannapieco, J. ZuHone, Astrophysics with the Spatially and Spectrally Resolved Sunyaev-Zeldovich Effects: A Millimetre/Submillimetre Probe of the Warm and Hot Universe. ArXiv e-prints (2018) R. Mushotzky, AXIS: a probe class next generation high angular resolution x-ray imaging satellite, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 10699, 2018, p. 1069929. doi:10.1117/12.2310003 D. Nagai, E.T. Lau, Gas Clumping in the Outskirts of ΛCDM Clusters. ApJ 731, 10 (2011). doi:10.1088/2041-8205/731/1/L10 D. Nagai, A.V. Kravtsov, A. Vikhlinin, Effects of Galaxy Formation on Thermodynamics of the Intracluster Medium. ApJ 668, 1–14 (2007). doi:10.1086/521328 K. Nelson, E.T. Lau, D. Nagai, Hydrodynamic Simulation of Non-thermal Pressure Profiles of Galaxy Clusters. ApJ 792, 25 (2014). doi:10.1088/0004-637X/792/1/25 K. Nelson, D.H. Rudd, L. Shaw, D. Nagai, Evolution of the Merger-induced Hydrostatic Mass Bias in Galaxy Clusters. ApJ 751, 121 (2012). doi:10.1088/0004-637X/751/2/121 K. Nelson, E.T. Lau, D. Nagai, D.H. Rudd, L. Yu, Weighing Galaxy Clusters with Gas. II. On the Origin of Hydrostatic Mass Bias in ΛCDM Galaxy Clusters. ApJ 782, 107 (2014a). doi:10.1088/0004-637X/782/2/107 K. Nelson, E.T. Lau, D. Nagai, D.H. Rudd, L. Yu, Weighing Galaxy Clusters with Gas. II. On the Origin of Hydrostatic Mass Bias in ΛCDM Galaxy Clusters. ApJ 782, 107 (2014b). doi:10.1088/0004-637X/782/2/107 D.M. Neumann, M. Arnaud, Regularity in the X-ray surface brightness profiles of galaxy clusters and the M-T relation. A&A 348, 711–727 (1999) A.J. Nishizawa, M. Oguri, T. Oogi, S. More, T. Nishimichi, M. Nagashima, Y.-T. Lin, R. Mandelbaum, M. Takada, N. Bahcall, J. Coupon, S. Huang, H.-Y. Jian, Y. Komiyama, A. Leauthaud, L. Lin, H. Miyatake, S. Miyazaki, M. Tanaka, First results on the cluster galaxy population from the Subaru Hyper Suprime-Cam survey. II. Faint end color- magnitude diagrams and radial profiles of red and blue galaxies at 0.1 < z < 1.1. PASJ 70, 24 (2018). doi:10.1093/pasj/psx106 J.M. Nugent, X. Dai, M. Sun, $Suzaku$ Measurements of Hot Halo Emission at Outskirts for Two Poor Galaxy Groups: NGC 3402 and NGC 5129. arXiv e-prints, 1704–08266 (2017) M. Oguri, T. Hamana, Detailed cluster lensing profiles at large radii and the impact on cluster weak lensing studies. MNRAS 414, 1851–1861 (2011). doi:10.1111/j.1365- 2966.2011.18481.x M. Oguri, M. Takada, K. Umetsu, T. Broadhurst, Can the Steep Mass Profile of A1689 Be Explained by a Triaxial Dark Halo? ApJ 632, 841–846 (2005). doi:10.1086/452629 T. Ohashi, Y. Ishisaki, Y. Ezoe, S. Yamada, G. Kuromaru, S. Suzuki, Y. Tawara, I. Mitsuishi, Y. Babazaki, K. Mitsuda, N.Y. Yamasaki, Y. Takei, R. Yamamoto, T. Hayashi, N. Ota, R.L. Kelley, K. Sakai, DIOS: the dark baryon exploring mission, in Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray. Proc. SPIE, vol. 9905, 2016, p. 99051. doi:10.1117/12.2232274 T. Ohashi, Y. Ishisaki, Y. Ezoe, S. Yamada, R. Hayakawa, K. Nunomura, K. Sato, Y. Tawara, I. Mitsuishi, K. Ohtsuka, K. Mitsuda, N.Y. Yamasaki, T. Kikuchi, T. Hayashi, H. Mu- ramatsu, Y. Nakashima, N. Ota, K. Osato, Y. Ichinohe, M.E. Eckart, S.R. Bandler, R.L. Kelley, C.A. Kilbourne, Super DIOS: future x-ray spectroscopic mission to search for dark baryons, in Society of Photo-Optical Instrumentation Engineers (SPIE) Con- 52 Stephen Walker et al.

ference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 10699, 2018, p. 1069928. doi:10.1117/12.2313122 N. Okabe, K. Umetsu, Subaru Weak Lensing Study of Seven Merging Clusters: Distributions of Mass and Baryons. PASJ 60, 345–375 (2008) N. Okabe, T. Futamase, M. Kajisawa, R. Kuroshima, Subaru Weak-lensing Survey of Dark Matter Subhalos in the Coma Cluster: Subhalo Mass Function and Statistical Properties. ApJ 784, 90 (2014a). doi:10.1088/0004-637X/784/2/90 N. Okabe, K. Umetsu, T. Tamura, Y. Fujita, M. Takizawa, Y.-Y. Zhang, K. Matsushita, T. Hamana, Y. Fukazawa, T. Futamase, M. Kawaharada, S. Miyazaki, Y. Mochizuki, K. Nakazawa, T. Ohashi, N. Ota, T. Sasaki, K. Sato, S. Tam, Universal profiles of the intracluster medium from Suzaku X-ray and Subaru weak-lensing observations*. PASJ 66, 99 (2014b). doi:10.1093/pasj/psu075 N. Okabe, H. Akamatsu, J. Kakuwa, Y. Fujita, Y. Zhang, M. Tanaka, K. Umetsu, Radio relics tracing the projected mass distribution in CIZA J2242.8+5301*. PASJ (2015). doi:10.1093/pasj/psv085 S. Padin, Inexpensive mount for a large millimeter-wavelength telescope. Appl. Opt. 53, 4431–4439 (2014). doi:10.1364/AO.53.004431 S.C. Parshley, J. Kronshage, J. Blair, T. Herter, M. Nolta, G.J. Stacey, A. Bazarko, F. Bertoldi, R. Bustos, D.B. Campbell, S. Chapman, N. Cothard, M. Devlin, J. Erler, M. Fich, P.A. Gallardo, R. Giovanelli, U. Graf, S. Gramke, M.P. Haynes, R. Hills, M. Limon, J.G. Mangum, J. McMahon, M.D. Niemack, T. Nikola, M. Omlor, D.A. Riechers, K. Steeger, J. Stutzki, E.M. Vavagiakis, CCAT-prime: a novel telescope for sub-millimeter astronomy, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 10700, 2018, p. 107005. doi:10.1117/12.2314046 A. Pillepich, C. Porciani, T.H. Reiprich, The X-ray cluster survey with eRosita: forecasts for cosmology, cluster physics and primordial non-Gaussianity. MNRAS 422, 44–69 (2012). doi:10.1111/j.1365-2966.2012.20443.x Planck Collaboration, Planck 2015 results. XIII. Cosmological parameters. ArXiv e-prints (2015) Planck Collaboration, P.A.R. Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, M. Ash- down, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A.J. Banday, et al., Planck 2013 results. XX. Cosmology from Sunyaev-Zeldovich cluster counts. ArXiv e-prints (2013a) Planck Collaboration, P.A.R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, F. Atrio- Barandela, J. Aumont, C. Baccigalupi, A. Balbi, A.J. Banday, et al., Planck intermediate results. V. Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect. A&A 550, 131 (2013b). doi:10.1051/0004-6361/201220040 S. Planelles, D. Fabjan, S. Borgani, G. Murante, E. Rasia, V. Biffi, N. Truong, C. Ragone- Figueroa, G.L. Granato, K. Dolag, E. Pierpaoli, A.M. Beck, L.K. Steinborn, M. Gaspari, Pressure of the hot gas in simulations of galaxy clusters. MNRAS 467, 3827–3847 (2017). doi:10.1093/mnras/stx318 G.W. Pratt, M. Arnaud, R. Piffaretti, H. Böhringer, T.J. Ponman, J.H. Croston, G.M. Voit, S. Borgani, R.G. Bower, Gas entropy in a representative sample of nearby X-ray galaxy clusters (REXCESS): relationship to gas mass fraction. A&A 511, 85 (2010). doi:10.1051/0004-6361/200913309 P. Predehl, R. Andritschke, H. Böhringer,W. Bornemann, H. Bräuninger,H. Brunner, M. Brusa, W. Burkert, V. Burwitz, N. Cappelluti, E. Churazov, K. Dennerl, J. Eder, J. Elbs, M. Freyberg, et al., eROSITA on SRG, in Space Telescopes and Instrumentation 2010: Ultraviolet to Gamma Ray. Proc. SPIE, vol. 7732, 2010, p. 77320. doi:10.1117/12.856577 E. Rasia, S. Ettori, L. Moscardini, P. Mazzotta, S. Borgani, K. Dolag, G. Tormen, L.M. Cheng, A. Diaferio, Systematics in the X-ray cluster mass estimators. MNRAS 369, 2013–2024 (2006) T.H. Reiprich, D.S. Hudson, Y. Zhang, K. Sato, Y. Ishisaki, A. Hoshino, T. Ohashi, N. Ota, Y. Fujita, Suzaku measurement of Abell 2204’s intracluster gas temperature profile out to 1800 kpc. A&A 501, 899–905 (2009). doi:10.1051/0004-6361/200810404 T.H. Reiprich, K. Basu, S. Ettori, H. Israel, L. Lovisari, S. Molendi, E. Pointecouteau, M. Roncarelli, Outskirts of Galaxy Clusters. Space Sci. Rev. 177, 195–245 (2013). doi:10.1007/s11214-013-9983-8 The Physics of Galaxy Cluster Outskirts 53

C.E. Romero, B.S. Mason, J. Sayers, T. Mroczkowski, C. Sarazin, M. Donahue, A. Baldi, T.E. Clarke, A.H. Young, J. Sievers, S.R. Dicker, E.D. Reese, N. Czakon, M. Devlin, P.M. Korngut, S. Golwala, Galaxy Cluster Pressure Profiles as Determined by Sunyaev Zel’dovich Effect Observations with MUSTANG and Bolocam. II. Joint Analysis of 14 Clusters. ApJ 838, 86 (2017). doi:10.3847/1538-4357/aa643f C. Romero, M. McWilliam, J.-F. Mac´ıas-Pérez,R. Adam, P. Ade, P. André,H. Aussel, A. Beelen, A. Benoˆıt,A. Bideaud, N. Billot, O. Bourrion, M. Calvo, A. Catalano, G. Coif- fard, B. Comis, M. de Petris, F.-X. Désert,S. Doyle, J. Goupy, C. Kramer, G. Lagache, S. Leclercq, J.-F. Lestrade, P. Mauskopf, F. Mayet, A. Monfardini, E. Pascale, L. Perotto, G. Pisano, N. Ponthieu, V. Revéret,A. Ritacco, H. Roussel, F. Ruppin, K. Schuster, A. Sievers, S. Triqueneaux, C. Tucker, R. Zylka, A multi-instrument non-parametric reconstruction of the electron pressure profile in the galaxy cluster CLJ1226.9+3332. A&A 612, 39 (2018). doi:10.1051/0004-6361/201731599 M. Roncarelli, S. Ettori, S. Borgani, K. Dolag, D. Fabjan, L. Moscardini, Large-scale inhomogeneities of the intracluster medium: improving mass estimates using the observed azimuthal scatter. MNRAS 432, 3030–3046 (2013). doi:10.1093/mnras/stt654 M. Rossetti, D. Eckert, S. De Grandi, F. Gastaldello, S. Ghizzardi, E. Roediger, S. Molendi, Abell 2142 at large scales: An extreme case for sloshing? A&A 556, 44 (2013). doi:10.1051/0004-6361/201321319 E. Rozo, E.S. Rykoff, redMaPPer II: X-Ray and SZ Performance Benchmarks for the SDSS Catalog. ApJ 783, 80 (2014). doi:10.1088/0004-637X/783/2/80 D.H. Rudd, D. Nagai, Non-equilibrium Electrons and the Sunyaev-Zel’dovich Effect of Galaxy Clusters. accepted to the ApJL (2009) E.S. Rykoff, E. Rozo, M.T. Busha, C.E. Cunha, A. Finoguenov, A. Evrard, J. Hao, B.P. Koester, A. Leauthaud, B. Nord, M. Pierre, R. Reddick, T. Sadibekova, E.S. Sheldon, R.H. Wechsler, redMaPPer. I. Algorithm and SDSS DR8 Catalog. ApJ 785, 104 (2014). doi:10.1088/0004-637X/785/2/104 D. Ryu, H. Kang, E. Hallman, T.W. Jones, Cosmological Shock Waves and Their Role in the Large-Scale Structure of the Universe. ApJ 593, 599–610 (2003). doi:10.1086/376723 T. Sasaki, K. Matsushita, K. Sato, N. Okabe, Suzaku Observations of Subhalos in the Coma Cluster. ApJ 806, 123 (2015). doi:10.1088/0004-637X/806/1/123 T. Sasaki, K. Matsushita, K. Sato, N. Okabe, X-ray observations of a subhalo associated with the NGC 4839 group infalling toward the Coma cluster. PASJ 68, 85 (2016). doi:10.1093/pasj/psw078 K. Sato, R.L. Kelley, Y. Takei, T. Tamura, N.Y. Yamasaki, T. Ohashi, A. Gupta, M. Galeazzi, Study of the Intracluster and Intergalactic Medium in the Sculptor Super- cluster with Suzaku. PASJ 62, 1423–1433 (2010). doi:10.1093/pasj/62.6.1423 K. Sato, K. Matsushita, N.Y. Yamasaki, S. Sasaki, T. Ohashi, Temperature and entropy profiles to the virial radius of the Abell 1246 cluster observed with Suzaku. PASJ 66, 85 (2014). doi:10.1093/pasj/psu061 T. Sato, T. Sasaki, K. Matsushita, E. Sakuma, K. Sato, Y. Fujita, N. Okabe, Y. Fukazawa, K. Ichikawa, M. Kawaharada, K. Nakazawa, T. Ohashi, N. Ota, M. Takizawa, T. Tamura, Suzaku Observations of the Hydra A Cluster out to the Virial Radius. PASJ 64, 95 (2012) J. Sayers, N.G. Czakon, A. Mantz, S.R. Golwala, S. Ameglio, T.P. Downes, P.M. Koch, K.-Y. Lin, B.J. Maughan, S.M. Molnar, L. Moustakas, T. Mroczkowski, E. Pierpaoli, J.A. Shitanishi, S. Siegel, K. Umetsu, N. Van der Pyl, Sunyaev-Zel’dovich-measured Pressure Profiles from the Bolocam X-Ray/SZ Galaxy Cluster Sample. ApJ 768, 177 (2013). doi:10.1088/0004-637X/768/2/177 J. Sayers, S.R. Golwala, A.B. Mantz, J. Merten, S.M. Molnar, M. Naka, G. Pailet, E. Pier- paoli, S.R. Siegel, B. Wolman, A Comparison and Joint Analysis of Sunyaev-Zel’dovich Effect Measurements from Planck and Bolocam for a Set of 47 Massive Galaxy Clusters. ApJ 832, 26 (2016). doi:10.3847/0004-637X/832/1/26 F. Sembolini, P.J. Elahi, F.R. Pearce, C. Power, A. Knebe, S.T. Kay, W. Cui, G. Yepes, A.M. Beck, S. Borgani, D. Cunnama, R. Davé,S. February, S. Huang, N. Katz, I.G. McCarthy, G. Murante, R.D.A. Newton, V. Perret, E. Puchwein, A. Saro, J. Schaye, R. Teyssier, nIFTy galaxy cluster simulations - II. Radiative models. MNRAS 459, 2973–2991 (2016). doi:10.1093/mnras/stw800 54 Stephen Walker et al.

M. Sereno, C. Giocoli, L. Izzo, F. Marulli, A. Veropalumbo, S. Ettori, L. Moscardini, G. Covone, A. Ferragamo, R. Barrena, A. Streblyanska, Gravitational lensing detection of an extremely dense environment around a galaxy cluster. Nature Astronomy (2018). doi:10.1038/s41550-018-0508-y H. Shan, J.-P. Kneib, C. Tao, Z. Fan, M. Jauzac, M. Limousin, R. Massey, J. Rhodes, K. Thanjavur, H.J. McCracken, Weak Lensing Measurement of Galaxy Clusters in the CFHTLS-Wide Survey. ApJ 748, 56 (2012). doi:10.1088/0004-637X/748/1/56 X. Shi, Locations of accretion shocks around galaxy clusters and the ICM properties: insights from self-similar spherical collapse with arbitrary mass accretion rates. MNRAS 461, 1804–1815 (2016a). doi:10.1093/mnras/stw1418 X. Shi, The outer proﬁle of dark matter haloes: an analytical approach. MNRAS 459, 3711– 3720 (2016b). doi:10.1093/mnras/stw925 X. Shi, E. Komatsu, Analytical model for non-thermal pressure in galaxy clusters. MNRAS 442, 521–532 (2014). doi:10.1093/mnras/stu858 X. Shi, E. Komatsu, K. Nelson, D. Nagai, Analytical model for non-thermal pressure in galaxy clusters - II. Comparison with cosmological hydrodynamics simulation. MNRAS 448, 1020–1029 (2015). doi:10.1093/mnras/stv036 X. Shi, E. Komatsu, D. Nagai, E.T. Lau, Analytical model for non-thermal pressure in galaxy clusters - III. Removing the hydrostatic mass bias. MNRAS 455, 2936–2944 (2016). doi:10.1093/mnras/stv2504 S.R. Siegel, J. Sayers, A. Mahdavi, M. Donahue, J. Merten, A. Zitrin, M. Meneghetti, K. Umetsu, N.G. Czakon, S.R. Golwala, M. Postman, P.M. Koch, A.M. Koekemoer, K.- Y. Lin, P. Melchior, S.M. Molnar, L. Moustakas, T.K. Mroczkowski, E. Pierpaoli, J. Shitanishi, Constraints on the Mass, Concentration, and Nonthermal Pressure Support of Six CLASH Clusters from a Joint Analysis of X-Ray, SZ, and Lensing Data. ApJ 861, 71 (2018). doi:10.3847/1538-4357/aac5f8 A. Simionescu, S.W. Allen, A. Mantz, N. Werner, Y. Takei, R.G. Morris, A.C. Fabian, J.S. Sanders, P.E.J. Nulsen, M.R. George, G.B. Taylor, Baryons at the Edge of the X-ray- Brightest Galaxy Cluster. Science 331, 1576 (2011). doi:10.1126/science.1200331 A. Simionescu, N. Werner, O. Urban, S.W. Allen, A.C. Fabian, J.S. Sanders, A. Mantz, P.E.J. Nulsen, Y. Takei, Large-scale Motions in the Perseus Galaxy Cluster. ApJ 757, 182 (2012). doi:10.1088/0004-637X/757/2/182 A. Simionescu, N. Werner, O. Urban, S.W. Allen, A.C. Fabian, A. Mantz, K. Matsushita, P.E.J. Nulsen, J.S. Sanders, T. Sasaki, T. Sato, Y. Takei, S.A. Walker, Thermodynamics of the Coma Cluster Outskirts. ApJ 775, 4 (2013). doi:10.1088/0004-637X/775/1/4 A. Simionescu, N. Werner, A. Mantz, S.W. Allen, O. Urban, Witnessing the growth of the nearest galaxy cluster: thermodynamics of the Virgo Cluster outskirts. MNRAS 469, 1476–1495 (2017). doi:10.1093/mnras/stx919 L. Spitzer, Physics of Fully Ionized Gases (Physics of Fully Ionized Gases, New York: In- terscience (2nd edition), 1962, ???, 1962) G.J. Stacey, M. Aravena, K. Basu, N. Battaglia, B. Beringue, F. Bertoldi, J.R. Bond, P. Breysse, R. Bustos, S. Chapman, D.T. Chung, N. Cothard, J. Erler, M. Fich, S. Foreman, P. Gallardo, R. Giovanelli, U.U. Graf, M.P. Haynes, R. Herrera-Camus, T.L. Herter, R. Hloˇzek,D. Johnstone, L. Keating, B. Magnelli, D. Meerburg, J. Meyers, N. Murray, M. Niemack, T. Nikola, M. Nolta, S.C. Parshley, D.A. Riechers, P. Schilke, D. Scott, G. Stein, J. Stevens, J. Stutzki, E.M. Vavagiakis, M.P. Viero, CCAT-Prime: science with an ultra-wideﬁeld submillimeter observatory on Cerro Chajnantor, in Soci- ety of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 10700, 2018, p. 107001. doi:10.1117/12.2314031 Y. Su, R.E. White III, E.D. Miller, Suzaku Observations of the X-Ray Brightest Fossil Group ESO 3060170. ApJ 775, 89 (2013). doi:10.1088/0004-637X/775/2/89 Y. Su, D. Buote, F. Gastaldello, F. Brighenti, The Entire Virial Radius of the Fossil Cluster RX J1159+5531: I. Gas Properties. ApJ 805, 104 (2015). doi:10.1088/0004- 637X/805/2/104 R.A. Sunyaev, Y.B. Zeldovich, The observation of relic radiation as a test of the nature of x-ray radiation from the clusters of galaxies. Comments Astrophys. Space Phys. 4, 173 (1972) The Physics of Galaxy Cluster Outskirts 55

H. Tanimura, G. Hinshaw, I.G. McCarthy, L. Van Waerbeke, Y.-Z. Ma, A. Mead, A. Hojjati, T. Tröster,A Search for Warm/Hot Gas Filaments Between Pairs of SDSS Luminous Red Galaxies. ArXiv e-prints (2017) M. Tashiro, H. Maejima, K. Toda, R. Kelley, L. Reichenthal, J. Lobell, R. Petre, M. Guainazzi, E. Costantini, M. Edison, R. Fujimoto, M. Grim, K. Hayashida, J.-W. den Herder, Y. Ishisaki, S. Paltani, K. Matsushita, K. Mori, G. Sneiderman, Y. Takei, Y. Ter- ada, H. Tomida, H. Akamatsu, L. Angelini, Y. Arai, H. Awaki, I. Babyk, A. Bamba, P. Barfknecht, K. Barnstable, T. Bialas, B. Blagojevic, J. Bonafede, C. Brambora, L. Bren- neman, G. Brown, K. Brown, L. Burns, E. Canavan, T. Carnahan, M. Chiao, B. Comber, L. Corrales, C. de Vries, J. Dercksen, M. Diaz-Trigo, T. Dillard, M. DiPirro, C. Done, T. Dotani, K. Ebisawa, M. Eckart, T. Enoto, Y. Ezoe, C. Ferrigno, Y. Fukazawa, Y. Fujita, A. Furuzawa, L. Gallo, S. Graham, L. Gu, K. Hagino, K. Hamaguchi, I. Hatsukade, D. Hawes, T. Hayashi, C. Hegarty, N. Hell, J. Hiraga, E. Hodges-Kluck, M. Holland, A. Hornschemeier, A. Hoshino, Y. Ichinohe, R. Iizuka, K. Ishibashi, M. Ishida, K. Ishikawa, K. Ishimura, B. James, T. Kallman, E. Kara, S. Katsuda, S. Kenyon, C. Kilbourne, M. Kimball, T. Kitaguti, S. Kitamoto, S. Kobayashi, T. Kohmura, S. Koyama, A. Kubota, M. Leutenegger, T. Lockard, M. Loewenstein, Y. Maeda, L. Marbley, M. Markevitch, H. Matsumoto, K. Matsuzaki, D. McCammon, B. McNamara, J. Miko, E. Miller, J. Miller, K. Minesugi, I. Mitsuishi, T. Mizuno, H. Mori, K. Mukai, H. Murakami, R. Mushotzky, H. Nakajima, H. Nakamura, S. Nakashima, K. Nakazawa, C. Natsukari, K. Nigo, Y. Nishioka, K. Nobukawa, M. Nobukawa, H. Noda, H. Odaka, M. Ogawa, T. Ohashi, M. Ohno, M. Ohta, T. Okajima, A. Okamoto, M. Onizuka, N. Ota, M. Ozaki, P. Plucin- sky, F.S. Porter, K. Pottschmidt, K. Sato, R. Sato, M. Sawada, H. Seta, K. Shelton, Y. Shibano, M. Shida, M. Shidatsu, P. Shirron, A. Simionescu, R. Smith, K. Someya, Y. Soong, Y. Suagawara, A. Szymkowiak, H. Takahashi, T. Tamagawa, T. Tamura, T. Tanaka, Y. Terashima, Y. Tsuboi, M. Tsujimoto, H. Tsunemi, T. Tsuru, H. Uchida, H. Uchiyama, Y. Ueda, S. Uno, T. Walsh, S. Watanabe, B. Williams, R. Wolfs, M. Wright, S. Yamada, H. Yamaguchi, K. Yamaoka, N. Yamasaki, S. Yamauchi, M. Yamauchi, K. Yanagase, T. Yaqoob, S. Yasuda, N. Yoshioka, J. Zabala, Z. Irina, Concept of the X-ray Astronomy Recovery Mission, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 10699, 2018, p. 1069922. doi:10.1117/12.2309455 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, The XMM Cluster Outskirts Project (X-COP): Physical conditions of Abell 2142 up to the virial radius. A&A 595, 42 (2016). doi:10.1051/0004-6361/201628183 The Simons Observatory Collaboration, P. Ade, J. Aguirre, Z. Ahmed, S. Aiola, A. Ali, D. Alonso, M.A. Alvarez, K. Arnold, P. Ashton, et al., The Simons Observatory: Science goals and forecasts. ArXiv e-prints (2018) S. Thölken, L. Lovisari, T.H. Reiprich, J. Hasenbusch, X-ray analysis of the galaxy group UGC 03957 beyond R200 with Suzaku. A&A 592, 37 (2016). doi:10.1051/0004- 6361/201527608 R.J. Thornton, P.A.R. Ade, S. Aiola, F.E. Angilè,M. Amiri, J.A. Beall, D.T. Becker, H.- M. Cho, S.K. Choi, P. Corlies, K.P. Coughlin, R. Datta, M.J. Devlin, S.R. Dicker, R. Dünner,J.W. Fowler, A.E. Fox, P.A. Gallardo, J. Gao, E. Grace, M. Halpern, M. Hasselfield, S.W. Henderson, G.C. Hilton, A.D. Hincks, S.P. Ho, J. Hubmayr, K.D. Irwin, J. Klein, B. Koopman, D. Li, T. Louis, M. Lungu, L. Maurin, J. McMahon, C.D. Munson, S. Naess, F. Nati, L. Newburgh, J. Nibarger, M.D. Niemack, P. Niraula, M.R. Nolta, L.A. Page, C.G. Pappas, A. Schillaci, B.L. Schmitt, N. Sehgal, J.L. Sievers, S.M. Simon, S.T. Staggs, C. Tucker, M. Uehara, J. van Lanen, J.T. Ward, E.J. Wollack, The Atacama Cosmology Telescope: The Polarization-sensitive ACTPol Instrument. ApJS 227, 21 (2016). doi:10.3847/1538-4365/227/2/21 K. Umetsu, B. Diemer, Lensing Constraints on the Mass Profile Shape and the Splashback Radius of Galaxy Clusters. ApJ 836, 231 (2017). doi:10.3847/1538-4357/aa5c90 O. Urban, N. Werner, A. Simionescu, S.W. Allen, H. Böhringer, X-ray spectroscopy of the Virgo Cluster out to the virial radius. MNRAS 414, 2101–2111 (2011a). doi:10.1111/j.1365-2966.2011.18526.x O. Urban, N. Werner, A. Simionescu, S.W. Allen, H. Böhringer, X-ray spectroscopy 56 Stephen Walker et al.

of the Virgo Cluster out to the virial radius. MNRAS 414, 2101–2111 (2011b). doi:10.1111/j.1365-2966.2011.18526.x O. Urban, A. Simionescu, N. Werner, S.W. Allen, S. Ehlert, I. Zhuravleva, R.G. Morris, A.C. Fabian, A. Mantz, P.E.J. Nulsen, J.S. Sanders, Y. Takei, Azimuthally resolved X-ray spectroscopy to the edge of the Perseus Cluster. MNRAS 437, 3939–3961 (2014). doi:10.1093/mnras/stt2209 F. Vazza, G. Brunetti, A. Kritsuk, R. Wagner, C. Gheller, M. Norman, Turbulent motions and shocks waves in galaxy clusters simulated with adaptive mesh refinement. A&A 504, 33–43 (2009). doi:10.1051/0004-6361/200912535 F. Vazza, D. Eckert, A. Simionescu, M. Brüggen,S. Ettori, Properties of gas clumps and gas clumping factor in the intra-cluster medium. MNRAS 429, 799–814 (2013). doi:10.1093/mnras/sts375 A. Vikhlinin, A.V. Kravtsov, R.A. Burenin, H. Ebeling, W.R. Forman, A. Hornstrup, C. Jones, S.S. Murray, D. Nagai, H. Quintana, A. Voevodkin, Chandra Cluster Cos- mology Project III: Cosmological Parameter Constraints. ApJ 692, 1060–1074 (2009). doi:10.1088/0004-637X/692/2/1060 G.M. Voit, S.T. Kay, G.L. Bryan, The baseline intracluster entropy profile from gravitational structure formation. MNRAS 364, 909–916 (2005) S.A. Walker, A.C. Fabian, J.S. Sanders, Large-scale gas sloshing out to half the virial radius in the strongest cool core REXCESS galaxy cluster, RXJ2014.8-2430. MNRAS 441, 31–35 (2014). doi:10.1093/mnrasl/slu040 S.A. Walker, A.C. Fabian, J.S. Sanders, M.R. George, Further X-ray observations of the galaxy cluster PKS 0745-191 to the virial radius and beyond. MNRAS 424, 1826–1840 (2012a). doi:10.1111/j.1365-2966.2012.21282.x S.A. Walker, A.C. Fabian, J.S. Sanders, M.R. George, Galaxy cluster outskirts: a universal entropy profile for relaxed clusters? MNRAS 427, 45–49 (2012b). doi:10.1111/j.1745- 3933.2012.01342.x S.A. Walker, A.C. Fabian, J.S. Sanders, M.R. George, Y. Tawara, X-ray observations of the galaxy cluster Abell 2029 to the virial radius. MNRAS 422, 3503–3515 (2012c). doi:10.1111/j.1365-2966.2012.20860.x S.A. Walker, A.C. Fabian, J.S. Sanders, A. Simionescu, Y. Tawara, X-ray exploration of the outskirts of the nearby Centaurus cluster using Suzaku and Chandra. MNRAS 432, 554–569 (2013). doi:10.1093/mnras/stt497 S.A. Walker, J. Hlavacek-Larrondo, M. Gendron-Marsolais, A.C. Fabian, H. Intema, J.S. Sanders, J.T. Bamford, R. van Weeren, Is there a giant Kelvin-Helmholtz instabil- ity in the sloshing cold front of the Perseus cluster? MNRAS 468, 2506–2516 (2017). doi:10.1093/mnras/stx640 S.A. Walker, J. ZuHone, A. Fabian, J. Sanders, The split in the ancient cold front in the Perseus cluster. Nature Astronomy 2, 292–296 (2018). doi:10.1038/s41550-018-0401-8 S.D.M. White, U.G. Briel, J.P. Henry, X-ray archaeology in the Coma cluster. MNRAS 261, 8–12 (1993) K.-W. Wong, J.A. Irwin, D.R. Wik, M. Sun, C.L. Sarazin, Y. Fujita, T.H. Reiprich, Suzaku X-ray Observations of the Nearest Non-Cool Core Cluster, Antlia: Dynamically Young but with Remarkably Relaxed Outskirts. ApJ 829, 49 (2016). doi:10.3847/0004- 637X/829/1/49 Y. Zu, R. Mandelbaum, M. Simet, E. Rozo, E.S. Rykoff, On the level of cluster assembly bias in SDSS. MNRAS 470, 551–560 (2017). doi:10.1093/mnras/stx1264 J.A. ZuHone, M. Markevitch, D. Lee, Sloshing of the Magnetized Cool Gas in the Cores of Galaxy Clusters. ApJ 743, 16 (2011). doi:10.1088/0004-637X/743/1/16