{'fraction_non_alphanumeric': 0.038472212361965366, 'fraction_numerical': 0.026474251821637634, 'mean_word_length': 4.657368190448372, 'pattern_counts': {'":': 0, '<': 10, '<?xml version=': 0, '>': 15, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Astrophysics with the Spatially and Spectrally Resolved Sunyaev-Zeldovich Effects: A Millimetre/Submillimetre Probe of the Warm and Hot Universe

In recent years, observations of the Sunyaev-Zeldovich (SZ) effect have had significant cosmological implications and have begun to serve as a powerful and independent probe of the warm and hot gas that pervades the Universe. As a few pioneering studies have already shown, SZ observations both complement X-ray observations -- the traditional tool for studying the intra-cluster medium -- and bring unique capabilities for probing astrophysical processes at high redshifts and out to the low-density regions in the outskirts of galaxy clusters. Advances in SZ observations have largely been driven by developments in centimetre-, millimetre-, and submillimetre-wave instrumentation on ground-based facilities, with notable exceptions including results from the Planck satellite. Here we review the utility of the thermal, kinematic, relativistic, non-thermal, and polarised SZ effects for studies of galaxy clusters and other large scale structures, incorporating the many advances over the past two decades that have impacted SZ theory, simulations, and observations. We also discuss observational results, techniques, and challenges, and aim to give an overview and perspective on emerging opportunities, with the goal of highlighting some of the exciting new directions in this field.


Introduction
In the hierarchical picture of structure formation, hot, X-ray emitting plasma with temperatures exceeding 10 5 K trace the largest overdensities to form in large scale structure: galaxy clusters, galaxy groups, galaxies, and intergalactic filaments. While X-ray observations are the main tool for probing the emission from these structures, particularly at temperatures 10 7 K (or k B T e 1 keV), it has long been known that the free electrons in the ionised plasma scatter cosmic microwave background (CMB) photons to, on average, higher energies, producing several unique spectral signatures referred to as the Sunyaev-Zeldovich (SZ) effects, the beginnings of which were first calculated nearly 50 years ago (Zeldovich & Sunyaev 1969;Sunyaev & Zeldovich 1970). The thermal SZ effect, which is the dominant effect due to upscattering of CMB photons by the hot intra-cluster medium (ICM), was proposed in Sunyaev & Zel'dovich (1972), shortly after the discoveries of X-ray emission from clusters (Byram et al. 1966;Bradt et al. 1967; see e.g. Sarazin 1986 and other chapters of these proceedings for reviews), and of the CMB itself (Penzias & Wilson 1965), as a test of whether the CMB was truly cosmic in origin.
Over the past half century, observational studies relying on the SZ effects have matured from early detections and imaging of known, massive clusters to using SZ measurements as a powerful tool for detection in wide-field surveys and detailed astrophysical studies. In parallel, many more facets of the SZ effect have been discovered, and our understanding has been refined, while both simulations and deep X-ray observations have transformed our view of the intra-cluster, intra-group, circumgalactic, and warm-hot filamentary structures on cosmic scales. Yet we have only begun to scratch the surface. We take this opportunity to review the advances of the past two decades, since the highly cited and seminal reviews of Birkinshaw (1999) and Carlstrom et al. (2002), and discuss opportunities for the near future and beyond.

A brief history of SZ observations
The first attempts to measure the thermal SZ focused on the decrement, and were performed using single-dish telescopes operating at radio wavelengths near 15 GHz (Pariiskii 1973;Gull & Northover 1976;Lake & Partridge 1977;Birkinshaw et al. 1978). By the mid-1980's, measurements with the Owens Valley Radio Observatory (OVRO) 40-meter started to produce more reliable detections in several tens of hours of integration time (Birkinshaw et al. 1984). Interferometric observations were performed later with the Very Large Array (VLA; Moffet & Birkinshaw 1989) and the Ryle telescope Grainge et al. 1993). These two interferometers also provided the first SZ imaging and purported detections based solely on an SZ signal, including the so-called "dark clusters" (i.e., clusters with no X-ray counterparts), which upon further investigation appeared to be spurious (Richards et al. 1997;Jones et al. 1997). Subsequent higher frequency (30-90 GHz) interferometric measurements with the OVRO 10-meter array and the Berkeley-Illinois-Maryland Array (BIMA; Carlstrom et al. 1996) and its successors the Sunyaev-Zeldovich Array and the Combined Array for mm-wave Astronomy (SZA & CARMA, respectively; Muchovej et al. 2007)), the Arcminute MicroKelvin Interferometer (AMI; Zwart et al. 2008) and the Array for Microwave Background Anisotropy (AMiBA; Lin et al. 2009) began interferometric imaging of large statistical samples of known clusters. Around the same time, photometric instruments using bolometers such as the Sunyaev-Zeldovich Infrared Experiment (SuZIE;Wilbanks et al. 1994) and Diabolo (Désert et al. 1998), mounted at the focus of large single dish telescopes, allowed mm-wave measurements at ∼150 GHz, where the SZ decrement is strongest. The next generation of photometric imaging arrays, including Bolocam (Glenn et al. 1998) and the Atacama Pathfinder Experiment Sunyaev-Zeldovich instrument (APEX-SZ; Schwan et al. 2003), were able to build large statistical samples of cluster images using single dish measurements.
The first measurement of the SZ increment toward a cluster-A2163-was performed by the PRONAOS balloon-borne experiment at 350 GHz. When combined with Diabolo and SuZIE measurements, it delivered the first spectro-photometric characterisation of the SZ spectrum (Lamarre et al. 1998). Soon after, the first successful imaging of the SZ increment was demonstrated using SCUBA on the JCMT (Komatsu et al. 1999). In the following decade, higher-quality SZ increment imaging was made possible using the LABOCA bolometer camera, once again for the cluster A2163 (Nord et al. 2009), while further SCUBA observations enabled the first large statistical study in this part of the spectrum (Zemcov et al. 2007). Shortly thereafter, using the Herschel-SPIRE photometer, the first detection of the SZ effect in the far increment at > 600 GHz was published (Zemcov et al. 2010). Building upon these spectro-photometric measurements, the first high-resolution SZ spectral measurement was performed towards the cluster RX J1347.5-1145 using the Z-Spec instrument (Zemcov et al. 2012). Subsequent multi-band photometric measurements towards a high-velocity sub-cluster in MACS J0717.5+3745 showed the first highsignificance deviation from the classical tSZ effect spectrum (Sayers et al. 2013a;Adam et al. 2017b), confirming the low-significance, 2σ indication of a non-zero kSZ signal reported in Mroczkowski et al. (2012).
The first high-spatial resolution SZ images were produced by the Diabolo camera on the IRAM 30-meter telescope and the NOBA camera on the NRO 45-meter telescope. Targeting RX J1347.5-1145-which quickly superseded A2163 as the SZ target of choice-Diabolo delivered a 5 × 5 map of with ∼20 resolution (Pointecouteau et al. 1999(Pointecouteau et al. , 2001, and NOBA delivered a 2 × 2 maps with ∼13 resolution (Komatsu et al. 1999(Komatsu et al. , 2001. Both maps exhibit an excess SZ signal to the southeast of the X-ray peak and central active galactic nucleus (AGN). This peak was not seen in the X-ray by ROSAT PSPC (Schindler et al. 1995), but was soon confirmed by Chandra (Allen et al. 2001). At the same time, OVRO/BIMA obtained ≈20 resolution in some of their interferometric images, including hints of an offset in the SZ centroid RX J1347.5-1145 to the southeast, but these features were never analysed in detail, and the cluster was instead used for determinations of the angular diameter distance (e.g. Reese et al. 2002;Bonamente et al. 2006). The next significant improvement in angular resolution came with the MUSTANG camera on the 100-meter Green Bank Telescope (GBT; Dicker et al. 2008), which imaged shock-heated gas and pressure substructures in a handful of clusters at ∼ 9 resolution (e.g. Mason et al. 2010;Korngut et al. 2011). As discussed in more detail in § 4, recent interferometric imaging with the Atacama Large Millimeter/Submillimeter Array (ALMA) has been used to probe one of the shock features in the 'El Gordo' cluster at a resolution of 3.5 through observations of the SZ effect (Basu et al. 2016), while ALMA and Atacama Compact Array (ACA, § 6.1.1) further improved constraints on RX J1347.5-1145 (Kitayama et al. 2016;Ueda et al. 2018).
Nearly a decade ago, the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) obtained sufficiently high mapping speeds to detect previously unknown clusters in wide-field surveys based on their thermal SZ effect signals (Staniszewski et al. 2009;Menanteau et al. 2010). Subsequent to these ground-based surveys, the Planck satellite surveyed the full sky in 9 photometric bands spanning the range 30-850 GHz, delivering a final catalogue of roughly 2000 SZ-selected clusters (Planck Collaboration et al. 2016e). The Planck survey data have also been used to measure the SZ effect spectrum with the broadest frequency coverage to date (e.g., Planck Collaboration et al. 2011;Hurier 2016a;Erler et al. 2018).
It must be noted that the backdrop to the developments of the past two decades has been the transformational X-ray observations of the ICM with the Chandra and XMM-Newton X-ray observatories. Understanding the physical processes inside the ICM relates directly to our ability to find galaxy clusters in X-ray or SZ surveys, and to model the ICM accurately to infer cluster masses for cosmological applications. SZ observations -especially deep, targeted observations -are also helpful in probing a wide range of physical phenomena and the imprints of merger events, such as shocks and cold fronts, energy dissipation through turbulence, and random bulk motions within the ICM. It is thus important to continue along this trajectory and develop the next generation cm/mm/submm tools for both surveys and targeted observations, particularly to deliver higher spatial and spectral resolution.

Theoretical and computational developments
As mentioned previously, the multiple manifestations of the SZ effect did not come as one discovery, but as a long series of discoveries and refinements. Much of this progress was made in parallel with observational and instrumental developments, as well as advances in simulations and our understanding of ICM physics. We discuss the many aspects of the SZ effect in § 2, aiming to provide a comprehensive but concise guide for the reader, while the references within indicate the long history of developments.

Organisation of this review
This review is presented as follows. In § 2, we discuss the fundamental theory behind the multiple manifestations of, and nuances to, the SZ effects. In § 3, we discuss how the SZ effects can be used to study the thermodynamics of ICM. In § 4, we present recent advances in our understanding of ICM structure, such as AGN feedback, turbulence, and discontinuities. In § 5, we discuss the practical considerations of performing observations of the SZ effects. In § 6, we highlight several of what we consider to be the main subarcminute resolution observatories today, and discuss a few near-term instruments under development and some long-term instruments and projects under study. In § 7 we conclude with an outlook on what the future may bring.

Overview of the SZ effects
The Sunyaev-Zeldovich (SZ) effects (Sunyaev & Zel'dovich 1972;Sunyaev & Zeldovich 1980b) are caused by the scattering of CMB photons with the free electrons residing in the potential wells of clusters of galaxies and, more broadly speaking, the diffuse plasma at large scales. 1 This leads to several CMB signals in the direction of clusters that can be used to learn about ICM physics and cosmology.
The physics behind the SZ signals is quite simple. Electrons at rest with respect to the isotropic CMB produce no net effect, as the number of photons scattered in and out of the line of sight is the same. However, moving electrons can transfer some of their kinetic energy to the CMB photon field through the Doppler effect. This can be appreciated by studying the Compton scattering relation for the ratio of the scattered to initial photon frequency (e.g., Jauch & Rohrlich 1976) (1) Here β = /c is the speed of the scattering electron with Lorentz factor γ = 1/ 1 − β 2 in units of the speed of light, c; m e is the electron mass; h is the Planck constant; µ and µ are respectively the direction cosines of the incoming and scattered photon with respect to the incoming electron; and µ sc is the corresponding direction cosine between the incoming and scattered photons.
Since the typical energy of the CMB photons is very small when compared with that of the electrons, hν γm e c 2 , one can neglect the corresponding photon recoil correction, which is usually dominant in the classical Compton effect. Hence, only Doppler and aberration terms are relevant (no Lorentz factors appear in the right hand side of Eq. 1), and Klein-Nishina corrections, O(hν/m e c 2 ), can be omitted 2 . The maximal photon energy after the scattering event is, ν max = ν (1 + β)/(1 − β) > ν, for a photon being back-scattered in a head-on collision with the electron. Similarly, the minimal scattered photon energy is ν min = ν (1 − β)/(1 + β) < ν for the initial photon travelling in the same direction as the incoming electron and then being backscattered. However, these scattering events occupy a small phase-space volume, and for any given µ, up-scattering of the CMB photon occurs when µ > µ, or equivalently, when the scattered photon deflects towards the direction of the incoming electron. Furthermore, the angular distributions of the incoming electrons and photons play a crucial role for the net energy exchange, as we will see below.
As these simple arguments already illustrate, scattering by moving electrons leads to a change of the CMB intensity in the direction of galaxy clusters, with the spectral shape of the signal depending on the velocity distribution of the electrons. Thermal electrons, described by an isotropic (relativistic) Maxwell-Boltzmann distribution, give rise to the so-called thermal SZ (tSZ) effect ( § 2.1), while the cluster's bulk motion (δ-function in velocity space) causes the kinematic SZ (kSZ) effect ( § 2.2). Depending on the characteristic speed of the electrons, relativistic corrections can become important, yielding the relativistic SZ (rSZ) effect ( § 2.3). Finally, non-thermal velocity distributions (e.g., in the cocoons of radio galaxies, or turbulence and magnetic fields) can create the non-thermal SZ (ntSZ) effect ( § 2.4).
Although the physical origin of all these SZ signals is due to simple electron scattering, each of them has slightly different spectral and spatial dependence across the cluster. With future multi-frequency observations, covering both the low-(ν 150 GHz) and high-frequency (ν 220 GHz) parts of the CMB blackbody, we are thus in principle able to distinguish them. This will provide an exciting opportunity for extracting valuable information about the structure of the cluster's atmosphere and its gas physics. We now explain each of the SZ signals in turn, highlighting how to compute them and detailing which physical parameters they can inform.  ) from an arbitrary angle, and on average is up-scattered to higher energy (blue) by an electron (black). The largest energy is imparted on the photon when it is scattered into the direction of the incoming electron, and it is minimal when deflected into the direction opposite to the incoming electron. However, on average scattering constellations with 90 • angles between the particles are most relevant for the tSZ. The total momentum in the interaction is conserved, so the electron is essentially undeflected by the interaction.

The thermal SZ effect
As CMB photons pass through regions of hot thermal gas (see schematic representation in Figure 1), inverse Compton scattering moves them from the low-frequency region of the blackbody spectrum towards higher energies. In single-scattering events with electrons at speed β drawn from an isotropic velocity distribution there is no net effect, as the gains and losses average out to leading order, leaving a second order term. The average energy gained by a CMB photon in each scattering is determined by ∆ν/ν (4/3) β 2 4kT e /m e c 2 (e.g., Rybicki & Lightman 1979;Sazonov & Sunyaev 2000). In the last step, we used β 2 th /3 ≈ kT e /m e c 2 for a thermal (nonrelativistic) velocity distribution. Similarly, a narrow photon line broadens by ∆ν/ν (2/3) β 2 2kT e /m e c 2 in each scattering event. In the non-relativistic limit, both effects can be incorporated using the Kompaneets equation (Kompaneets 1956), which when applied to the case of SZ clusters 3 reduces to a simple diffusion equation in frequency-space. This approach was originally used by Zeldovich & Sunyaev (1969) to compute the distinctive tSZ distortion signal, often referred to as the Compton y-type distortion. The corresponding distortion is given by 4 : in terms of the CMB intensity. Here x = hν/k B T CMB ≈ ν/56.8 GHz, with T CMB denoting the temperature of the CMB, and the classical tSZ spectral function g(x) is defined implicitly in Eq.
(2). Assuming ∆I ν /I ν 1, one can use the derivative with respect to temperature of the Planck function to alternatively express the signal in terms of the effective CMB temperature, yielding: The function f (x), defined implicitly above, is the classical tSZ spectrum in terms of ∆T CMB . The change in the effective CMB temperature is proportional to the Compton-y parameter, which depends on the Thomson scattering optical depth, τ e , and temperature of the hot electron gas, T e , as Here σ T is the Thompson cross section, P e = n e k B T e is the pressure due to the electrons and n e is the number density of the electrons. The integral is performed over the proper distance along the line of sight. Thus the magnitude of the tSZ signal is a direct measure of the integrated line of sight pressure. Typical clusters contain electrons with kT e 5 − 10 keV, or kT e /m e c 2 0.01 − 0.02. The central optical depth can reach τ e = n e σ T dl 10 −2 , such that for massive clusters one can expect y 10 −4 (see e.g. the cluster outskirts review in these proceedings). The spectral shape of the y-distortion in terms of CMB intensity for y = 10 −4 is illustrated in Figure 2. The signal manifests itself as a deficit in the number of photons at frequencies below ν null ≈ 217 GHz (≈ 1.4 mm) and an increase above ν null (n.b. photon number is conserved by scattering). In the Rayleigh-Jeans limit (x 1) the change in the effective temperature ∆T reduces to ∆T/T ≈ −2y, while at high frequencies (x 1) one has ∆T/T ≈ y(x − 4).
An important property of the tSZ is its near redshift independence. The spectral shape remains unchanged and the tSZ does not suffer from redshift-dimming 5 . 4 An alternative derivation uses that the superposition of blackbodies with slightly different temperatures, as indeed caused by the scattering process, is no longer a blackbody (Zeldovich et al. 1972). In this case, the y-parameter is related to the temperature dispersion, y = 1 2 ∆T T 2 , induced by Doppler-shifts. 5 The CMB intensity increases towards higher redshifts and thus the tSZ signal starts off at a higher level at z > 0, remaining constant relative to the CMB intensity. Alternatively, this can be understood qualitatively as the scattering events producing a fractional change in the intensity of the CMB that would be constant for an observer at any given epoch. We assumed an optical depth τ e = 10 −2 and an overall Compton parameter y = 10 −4 . The dotted, dark red curve illustrates the shape of the unscattered CMB spectrum, which for comparison was scaled by a factor of 5 × 10 −4 .
This makes the tSZ a unique probe of the large-scale structure in the Universe (see e.g. Sunyaev & Zeldovich 1980a;Rephaeli 1995a;Birkinshaw 1999;Carlstrom et al. 2002;Kitayama 2014, for a number of excellent reviews). It can furthermore in principle be used to measure the expansion rate of the Universe through the combination with X-ray data, exploiting the differing density dependencies in their surface brightness integrals (see Equations 13 and 20) to infer the angular diameter distance to a cluster (Silk & White 1978;Cavaliere et al. 1979;Birkinshaw 1979;Hughes & Birkinshaw 1998;Battistelli et al. 2003).

The kinematic SZ effect
The kinematic SZ (kSZ) effect is due to scattering of CMB photons interacting with free electrons undergoing bulk motion relative to the CMB rest-frame (Sunyaev & Zeldovich 1980b, see the schematic diagram in Figure 3). In contrast to the tSZ, the velocity distribution of the electrons is anisotropic in the case (i.e., mono-directional), such that, upon averaging over all scattering angles of the photons drawn from the isotropic CMB, a linear order Doppler term ∝ β remains. Similar physics play a CMB photon bulk motion free electron hot plasma scattered photon Fig. 3 Diagram for the kinematic SZ. A CMB photon (red) enters the hot ICM (light blue) from an arbitrary angle, and, for this geometry, is up-scattered to higher energy (blue) by an electron (white dot) in a moving ICM. To first order order in β = /c only the line of sight projection of the cluster's bulk motion matters for the corresponding intensity change.
crucial role in the formation of the CMB temperature and polarisation anisotropies (e.g., Sunyaev & Zeldovich 1970;Peebles & Yu 1970;Hu & Sugiyama 1995). The kSZ effect produces a shift in CMB temperature in the direction n of a moving cluster, which can be written as in terms of the effective shift in CMB temperature, or in terms of the CMB intensity shift (where I 0 is defined in Eq. 3). This signal is indistinguishable from that of the hot and cold spots in the primary CMB, unless scaledependent information or correlations with other astronomical data are exploited. To leading order, only the line of sight component, β p, = n · β p = µ p β p , of the cluster's peculiar motion is relevant. The temperature shift is furthermore negative for a line of sight velocity away from the observer, and positive when the cluster approaches the observer (consistent with the convention that the z vector increases with redshift).
The parameter y kSZ (see Eq. 6), is often defined in the literature as a kSZ analogue to the Compton-y parameter (e.g. Ruan et al. 2013). The typical peculiar motions of clusters in the standard cosmological model are expected to be β p few×10 −3 (speed few×100 km s −1 ), such that y kSZ few×10 −5 at the map peak. This is about one order of magnitude smaller than the typical tSZ y-parameter (see Figure 2). By measuring this signal, one can in principle map out the large-scale motions of baryons in the Universe, constraining the growth of structure and testing isotropy/homogeneity of the Universe (Zhang & Stebbins 2011;Hand et al. 2012).

Relativistic corrections to the SZ effect
The tSZ and kSZ signals discussed above were obtained assuming non-relativistic speeds for the electrons. For typical bulk motions β p few × 10 −3 this is rather well justified. However, electrons in a thermal gas at temperature kT e 5 keV have typical speeds β 3kT e /m e c 2 0.1 − 0.2. In this case, the non-relativistic approximation (i.e., only including terms up to second order in β) for the tSZ derivation no longer suffices and (special-)relativistic corrections become relevant (e.g., Wright 1979;Fabbri 1981;Rephaeli 1995b), leading to changes of the tSZ and kSZ spectral shapes. Thus, relativistic corrections to the SZ signals ('rSZ' for short) can in principle provide additional information about the detailed temperature and velocity structure of the ICM, as we discuss below.
Relativistic corrections to both the tSZ and kSZ have been studied in detail using various methods to evaluate the Compton collision term. Direct numerical integration with various levels of analytical reductions (e.g., Wright 1979;Fabbri 1981;Rephaeli 1995b;Pointecouteau et al. 1998;Molnar & Birkinshaw 1999;Enßlin & Kaiser 2000;Dolgov et al. 2001;Nozawa et al. 2009) are generally time-consuming but least prone to errors or loss of precision. For a specified temperature range, the computation can be accelerated using fits to the numerical results Itoh & Nozawa 2004) or pre-computed basis functions (Chluba et al. 2012b). Insight into the physics of the problem can be gained through Taylor-series approximations to various orders in Θ e = kT e /m e c 2 and β p = p /c (Stebbins 1997;Challinor & Lasenby 1998;Sazonov & Sunyaev 1998;Itoh et al. 1998;Nozawa et al. 1998;Shimon & Rephaeli 2004;Nozawa & et al. 2006). However, at kT e 5 keV the expansions start to converge quite slowly, since the width of the scattering kernel quickly rises with T e . In particular in the Wien tail of the CMB spectrum this becomes problematic (e.g., Stebbins 1997;Itoh et al. 1998;Chluba et al. 2012b), as derivatives of an exponential are not well approximated by a sum of exponentials. Another kinematic correction is due to the motion of the observer, which can be added by performing a Lorentz transformation of the SZ signal from the CMB rest frame to the observer's frame (Chluba et al. 2005;Nozawa et al. 2005). This leads to a dipolar modulation of the cluster number counts across the sky (Chluba et al. 2005) and can also be interpreted as a distortion of the CMB dipole spectrum (Balashev et al. 2015).
After this broad-brush overview, let us discuss the physics of relativistic temperature corrections in more detail. One of the important effects is that the average energy shift and broadening per scattering both increase more quickly with temperature than in the non-relativistic limit. Including terms up to fourth order in Θ e , one finds (Sazonov & Sunyaev 2000;Chluba et al. 2012b) for the first two moments of the scattering kernel, which illustrates the effect. At higher temperature, it is thus no longer possible to assume |∆ν/ν| 1 in the scattering event, one of the key assumptions in the derivation of the Kompaneets equation. This implies that both higher order derivatives of the blackbody function 6 and higher order moments of the scattering kernel become important. Overall relativistic temperature corrections cause a broadening of the tSZ signal, with a systematic shift towards higher frequencies, as is illustrated in Figure 2. Through this additional dependence it is in principle possible to directly measure the temperature of the cluster (Wright 1979;Fabbri 1981;Rephaeli 1995b;Pointecouteau et al. 1998). As discussed further in §3.3.2, rSZ temperature determinations have been attempted for individual clusters (e.g., Hansen et al. 2002a;Prokhorov & Colafrancesco 2012;Chluba et al. 2013) and in stacking analyses (Hurier 2016b;Erler et al. 2018), albeit with large errors in both. In the future this could also become possible at the tSZ power spectrum level (Remazeilles et al. 2019) and for cluster number counts (Fan & Wu 2003). The relativistic tSZ in principle can also be used to measure the CMB temperature-redshift relation, which in non-standard cosmologies could depart from the standard T CMB ∝ (1 + z) scaling (e.g., Fabbri et al. 1978;Rephaeli 1980;Battistelli et al. 2002;Luzzi et al. 2009). Even if it is fairly difficult to create a change in the CMB temperature at late times (Chluba 2014) without violating CMB spectral distortion constraints from COBE/FIRAS (Mather et al. 1994;Fixsen et al. 1996), this is an interesting application of the rSZ.
Turning to the physics of relativistic corrections to the kSZ, just as those for the tSZ, additional higher order terms become relevant when evaluating the kernel moments. The temperature corrections to the kSZ (at leading order ∝ β p Θ e ) are illustrated in Figure 2. As with the tSZ effect, these again lead to a broadening and systematic shift of the kSZ signal towards higher frequencies. Mixed kinematic and temperature corrections have been considered for terms up to second order in β p (Sazonov & Sunyaev 1998;Nozawa et al. 1998;Shimon & Rephaeli 2004;Nozawa & et al. 2006;Chluba et al. 2012b). The bulk motion of the cluster breaks the isotropy of the CMB blackbody field at second order inducing a correction to the monopole ∝ β 2 p and quadrupole ∝ β 2 p (3µ 2 p − 1)/2 inside the cluster's rest frame. Through relativistic kSZ one can therefore in principle measure two projections of the cluster's velocity. The remaining azimuthal degeneracy can in principle be broken by considering pSZ ( § 2.5). For this the primordial quadrupole terms have to be carefully sub-6 Derivatives of the blackbody occupation number, n bb = 1/(e x − 1), can be given in closed form using (Chluba et al. 2012b).
tracted . Instead of Taylor series approximations of the Compton collision term, by considering the tSZ effect in the cluster frame for an anisotropic CMB photon field, one can account for the kinematic corrections through the Lorentz transformation, an approach that also ensures the correct interpretation of the optical depth as cluster-frame optical depth (Chluba et al. 2012b). Higher order corrections in β p can be added by using a multipole-dependent Kompaneets equation or anisotropic scattering kernels (Chluba et al. 2012a;, however, these modifications are expected to be negligible. All the aforementioned relativistic corrections (kinematic due to the cluster's and observer's motions as well as those due to temperature) can be efficiently computed using SZpack 7 (Chluba et al. 2012b(Chluba et al. , 2013 with complementary numerical and analytical methods implemented to ensure a large amount of flexibility. SZpack furthermore allows including the effect of spatial variations of the temperature and velocity fields (both along the line of sight and within the beam), which give rise to frequency-dependent morphological changes of the SZ signals (Chluba et al. 2013). In the future, the associated moments of the velocity and temperature field could become direct observables, and by combining with X-ray data we may be able to extract detailed information about the ICM structure.

The non-thermal SZ effect
In the discussion of the tSZ effect we only considered thermal distributions of electron momenta. However, the momentum distribution can be more complex and highly relativistic, e.g., having long power-law tails at high energies (e.g., Enßlin & Kaiser 2000). In this case, we refer to the associated distortions as non-thermal SZ (ntSZ) effect. Assuming a general 8 isotropic momentum distribution f (p), where p is the dimensionless electron momentum (i.e., p = p phys /m e c and γ = 1 + p 2 ) and the distribution f (p) has a normalisation ∞ 0 f (p)p 2 dp = 1, one can write the scattered CMB signal as 9 ∆I ν ≈ I 0 x 3 τ e ∞ 0 f (p)p 2 dp s m (p) −s m (p) P(s, p) n bb (x e s ) − n bb (x) ds.
and is normalised as s m −s m P(s, p) ds = 1. One furthermore finds the first moment, ∆ν/ν = s m −s m (e s − 1) P(s, p) ds = 4p 2 /3. We emphasise that Eq. (9) is only strictly valid when anisotropies in the radiation and velocity fields can be neglected. This assumption can in principle be violated by the presence of magnetic fields (Koch et al. 2003;Gopal & Roychowdhury 2010), pressure anisotropies (Khabibullin et al. 2018), anisotropies in the scattering medium , and kSZ effects, as mentioned above. Equation (9) allows us to compute the scattered signal for a given f (p). If we assume the electron energies follow a relativistic Maxwell-Boltzmann distribution, , we can reproduce the tSZ effect with relativistic temperature corrections. Here K 2 (x) is the modified Bessel function of the second kind, ensuring the correct normalisation of f (p). For power-law distributions, analytic expressions can be given in terms of incomplete β-functions (Enßlin & Kaiser 2000;Colafrancesco et al. 2003). Similarly, for mono-energetic electrons, f δ (p) = δ(p − p 0 )/p 2 , the outer integral becomes trivial and only one numerical integral over the blackbody distribution has to be carried out. In Figure 4, we show the ntSZ spectrum for this case. CMB photons are strongly up-scattered towards higher frequencies once the momentum exceeds p 1 (Enßlin & Kaiser 2000;Colafrancesco et al. 2003;Malu et al. 2017). This effect could allow us to search for the presence of aged radio plasmas and relativistic outflows from AGN using CMB measurements. The total contribution to the Comptonisation of CMB photons caused by relativistic plasmas could reach a level equivalent to y few × 10 −6 (Enßlin & Kaiser 2000;Shimon & Rephaeli 2002). However, the detailed shape of the ntSZ spectrum can become complex since it depends strongly on the energy distribution of the electrons doing the scattering. A direct comparison with the tSZ is thus not as straightforward. The ntSZ could also allow us to shed light on the nature of dark matter and annihilating particles (Colafrancesco 2004;Colafrancesco et al. 2006).

The polarised SZ effects
It is well-known that Thomson scattering of CMB photons by free electrons inside clusters leads to a small polarisation effect (e.g., Sunyaev & Zeldovich 1980b, 1981Gibilisco 1997;Kamionkowski & Loeb 1997;Audit & Simmons 1999;Sazonov & Sunyaev 1999;Challinor et al. 2000;Itoh et al. 2000;Shimon et al. 2009). The physics behind this signal is very similar to the way the primordial CMB E-mode polarisation patterns are created (Bond & Efstathiou 1984;Seljak & Zaldarriaga 1997;, the most important ingredient being the presence of a quadrupole anisotropy in the local radiation field. However, for hot electrons inside clusters the Comptonisation of CMB photon can lead to a distinct frequencydependence of the polarised SZ (pSZ). For this, the origin of the local quadrupole anisotropy plays a crucial role, as we discuss now.
The largest effect is due to the local primordial CMB quadrupole, leading to a polarisation amplitude 0.1 τ e Q in units of the CMB temperature (Kamionkowski & Loeb 1997;Sazonov & Sunyaev 1999). Here Q is the CMB quadrupole moment at the location of the cluster. The signal thus can reach a level of 10 −8 of the primary CMB temperature for rich clusters (τ e 0.01 and Q 10 −5 ), and could allow us to measure the CMB quadrupole at different locations in the Universe, thereby in principle circumventing the cosmic variance limit (Kamionkowski & Loeb 1997;Portsmouth 2004;Bunn 2006;Yasini & Pierpaoli 2016;Meyers et al. 2018). However, the frequency-dependence of this signal is identical to that of the primordial CMB polarisation anisotropies, such that knowledge about the cluster location and its redshift are required.
The second largest pSZ signal is due to second scattering corrections of the thermal ( τ e 2 Θ e ) and kinematic ( τ e 2 β p ) SZ signals (e.g., Sazonov & Sunyaev 1999). In both cases, the spectral-dependence follows that of the re-scattered tSZ and kSZ respectively. The signal related to the second scattering of the tSZ signal can be of similar order of magnitude as the one caused by the primordial CMB quadrupole, but since a scattering-induced quadrupole anisotropy is required (see § 2.6) it vanishes along the centre of the cluster unless asphericities or anisotropies in the medium are present (e.g., Sazonov & Sunyaev 1999;Puy et al. 2000;Shimon et al. 2009). For the second scattering pSZ induced by the kSZ only the tangential component of the clusters motion is relevant.
At second order in β p , a y-type CMB quadrupole is also induced by the cluster's motion. This is the original pSZ effect that was highlighted by Sunyaev & Zeldovich (1980b). Again only the tangential component of the cluster's velocity, β p,⊥ , is relevant, causing a quadrupole pattern that leads to a polarisation signal 0.1 τ e β 2 p,⊥ . This can in principle be used to measure the tangential velocity of the cluster's motion in the plane of the sky. Similarly, internal gas motions (e.g., Chluba 2001;Cooray & Chen 2002;Chluba & Mannheim 2002;Diego et al. 2003;Lavaux et al. 2004;Shimon et al. 2006;Maturi et al. 2007) can lead to complex polarisation patterns.
Yet another physical mechanism that could lead to polarisation of the SZ signal is associated with anisotropic distribution function of electrons (Khabibullin et al. 2018). The ICM is an astrophysical example of weakly collisional plasma where the Larmor frequencies of charged particles greatly exceed their collision frequencies. In such conditions the magnetic moments of particles are conserved between collisions and the evolving magnetic fields or heat fluxes can generate pressure anisotropies of particles. Therefore, the characteristic thermal velocities of electrons can differ along and perpendicular to the direction of the magnetic field, inducing a polarisation pattern in the CMB. The signal scales linearly with the optical depth of the region containing large-scale correlated anisotropy (e.g., along ubiquitous cold fronts in clusters), and with the degree of anisotropy itself. It has the same spectral dependence as the polarisation induced by cluster motion with respect to the CMB frame (kinematic SZ effect polarisation), but can be distinguished by its spatial pattern. The magnitude of the effect is on par with majority of other SZ polarisation signals considered here (for a useful summary, see Table 1 of Khabibullin et al. 2018). An increase of the effective electron collisionality due to plasma instabilities will reduce the effect. Such polarisation, therefore, may be an independent probe of the electron collisionality in the ICM, which is one of the key properties of a high-β weakly-collisional plasma 10 from the point of view of both astrophysics and plasma theory.
Even in the future, the aforementioned pSZ signals will be challenging to extract due to their intrinsic faintness, the limited polarisation purity of existing instrumentation technologies, and the large number of polarised astrophysical signals that could contaminate them (e.g. Sunyaev 1982). Beam depolarisation effects will furthermore render much of the polarisation signal unobservable for unresolved clusters. Nevertheless, the various pSZ signals may provide another avenue forward for detailed studies of ICM structure, and maybe become accessible through stacking analyses.

Multiple scattering effects
Another subdominant correction to the SZ signals is caused by multiple scattering events inside rich clusters (Sunyaev & Zeldovich 1980b;Sazonov & Sunyaev 1999;Molnar & Birkinshaw 1999;Dolgov et al. 2001;Itoh et al. 2001;Colafrancesco et al. 2003;Shimon & Rephaeli 2004). This correction is usually derived in the isotropic scattering approximation (ISA), which assumes that the radiation field locally remains isotropic. In this limit, to leading order, the contribution is suppressed by a factor of 10 y relative to the single-scattering tSZ signal, rendering it a 0.1% correction (Itoh et al. 2001). In the ISA no correction ∝ τ e arises.
However, when the scattering-induced anisotropy in the radiation field is included, the scenario changes slightly. In this case, corrections may become detectable, which even in a constant density sphere are caused by the variations of the photon's path in different directions . Hence, the spectral dependence of the multiple-scattering signal is modified, and a new contribution of order τ e /20 relative to the tSZ effect arises . The net signal depends explicitly on the considered line of sight and structure of the medium. Future measurements of the spatial and frequency dependence of the multiple-scattering SZ signal could thus help in the reconstruction of ICM density and temperature profiles. However, such measurements will remain very challenging for the foreseeable future.

ICM thermodynamics through the SZ effects
3.1 The "universal" pressure profile One of the key insights from modern hydrodynamical cosmological simulations is that the hot X-ray emitting plasma in galaxy clusters exhibits a remarkable degree of self-similarity (e.g., Nagai et al. 2007;Battaglia et al. 2010;Lau et al. 2015), where the ICM pressure profile is well characterised by a generalised Navarro, Frenk, & White (NFW) profile (Navarro et al. 1996). The use of the generalised NFW (gNFW) profile to describe pressure, the integral of which yields the tSZ signal (Eq. 5), was first proposed by Nagai et al. (2007): Here, the parameters (γ, β) are the central slope (r r s ) and outer slope (r r s ), respectively. The parameter α modulates how smoothly the slope changes from γ to β around r s , where x = r/r s , r s = r 500 /c 500 , and r 500 is defined as the radius within which the average overdensity is 500× greater than the critical density of the Universe at that redshift, ρ crit (z). 11 All three of the slopes are highly correlated with the value of the scale radius r s . This analytic profile has been widely used to measure pressure profiles using X-ray and SZ data. Measurements of ICM pressure profiles provide important information about the thermodynamic structure of the ICM, including the effects of AGN feedback, bulk and turbulent motions, substructures, and asphericity of clusters. Beyond r 500 , an increasing level of non-thermal pressure support at the level of 10 -30% is expected, depending on the dynamical state of clusters (e.g., Lau et al. 2009;Battaglia et al. 2012b;Nelson et al. 2014).
Motivated by Nagai et al. (2007), the gNFW parametrisation (given in Eq. 11) was first applied to SZ observations in Mroczkowski et al. (2009) and has since displaced the previously used β-model (Cavaliere & Fusco-Femiano 1976 in nearly all SZ-related pressure profile studies. A number of SZ and X-ray studies have since sought to measure or refine the azimuthally-averaged, "universal" pressure profile of the ICM, often attempting to constrain the slope parameters of the gNFW profile (Arnaud et al. 2010;Plagge et al. 2010;Bonamente et al. 2012;Planck Collaboration et al. 2013a;Sayers et al. 2013b;Eckert et al. 2013b;Adam et al. 2015;Sayers et al. 2016a;Ghirardini et al. 2017;Romero et al. 2017;Bourdin et al. 2017;Ruppin et al. 2018). The gNFW parameters are highly covariant, and thus the interest is often more on the overall profile shape, or the combinations of the parameters, rather than the individual parameter values themselves.
The majority of the above studies find agreement with the "universal" pressure profile presented in Arnaud et al. (2010) (A10), especially within r 500 . Where disagreement has been found within r 500 , systematic uncertainties and sample variance are likely to mitigate the purported tension. Beyond r 500 , studies have relied on, or heavily supplemented X-ray data with, tSZ data (Planck Collaboration et al. 2013a;Sayers et al. 2016b;Ghirardini et al. 2018a, e.g.). At r > r 500 , Ghirardini et al. (2018a) find higher pressure relative to A10, while Sayers et al. (2016b) find lower pressure. As the samples are disparate (notably in redshift range), one might take this as an early indication of evolution in the pressure profile, consistent with the hydrodynamical simulation results reported by Battaglia et al. (2012a).
The gNFW profiles are often fit to binned non-parametric pressure profiles (e.g., Plagge et al. 2010;Sayers et al. 2013b), where Plagge et al. (2010) stacked their sample to recover the non-parametric pressure profiles. With higher resolution and more sensitive SZ observations, non-parametric pressure profiles have been fit to the SZ data of individual clusters (Basu et al. 2010;Sayers et al. 2013b;Romero et al. 2018;Ruppin et al. 2018). An example of recent pressure profile constraints for the jointly-fit, multi-scale data in Romero et al. (2017Romero et al. ( , 2018) is shown in Figure 5.
The methods of constraining pressure profiles from SZ data are non-trivial and affect the potential biases and systematic uncertainties associated with these constraints. As discussed in § 5.2, ground-based observations suffer from atmospheric effects that typically limit the spatial scales recoverable from an observation to the instrumental field of view (FoV). There are several ways to mitigate these filtering effects, the main tools being: parametric and non-parametric forward modelling, deconvolution, and the Abel transform (see Silk & White 1978, for reference). If the noise and filtering behaviours are well characterised, one can attempt image deconvolution (Basu et al. 2010;Sayers et al. 2013b). From the deconvolved images, pressure profiles may be deprojected via methods used in X-ray studies, which include an 'onion-skin' method (multiple shells either jointly-fit or fit through a 'peeling' method) or Abel transform (e.g. as applied in Basu et al. 2010). However, the decon-  Romero et al. 2017Romero et al. , 2018 volution suffers at large angular scales, enough that Sayers et al. (2013b) opted to forward model pressure profiles calculated as power-law interpolations between several radii. Forward modelling requires, in its coarsest form, a grid search (e.g. Romero et al. 2015), or more often a Monte Carlo Markov Chain approach (e.g., Bonamente et al. 2004Bonamente et al. , 2006Olamaie et al. 2013;Sayers et al. 2013b;Ruppin et al. 2018).
Another implication of the FoV restricting the spatial scales recovered is that many SZ instruments do not access a wide dynamic range of angular scales. As highresolution SZ instruments were developed ( § 1.1), investigators soon wanted to combine SZ datasets to maximise the range of angular scales over which pressure profiles could be constrained. These studies (e.g. Sayers et al. 2016b;Romero et al. 2017) serve both to check for agreement of pressure profile shapes between SZ and X-ray data as well as overall agreement among the SZ datasets.
Given the complexity of extracting pressure profiles from SZ data, especially ground-based data, consistency checks from joint fits or identical analyses are critical. In comparing their results to previous results, Sayers et al. (2016b) note that differences in pressure profiles can derive from (1) sample selection, (2) instrumental biases, or (3) biases introduced in data processing. For single clusters, the SZ data currently appear consistent (e.g., Romero et al. 2018;Ruppin et al. 2018), but deeper observations of an overlapping sample between mm-wave instruments may yield differences in individual clusters. Additional comparisons across samples of clusters will serve this end as well.
Several current SZ instruments ( § 6) now recover sufficiently large angular scales so as to provide moderate resolution pressure profile constraints beyond r 500 and thus provide a clearer picture as to whether the pressure profile evolves with redshift.
These instruments also benefit from the complementary X-ray data. When combined with X-ray measurements of the ICM density, the pressure profile can also be used to infer other thermodynamic quantities such as temperature and entropy (see § 3.3.1 for further discussion). However, the temperature derivation is degenerate with cluster geometry, and second-order effects like pressure clumping and helium sedimentation could bias the resulting constraints.

SZ scaling relations
Since the thermal energy content of the cluster is determined primarily by the gravitational potential well of dark matter, the aperture (or volumetrically) integrated tSZ signal, is proportional to the thermal energy content of the ICM. Per beam, the tSZ provides calorimetry of the ICM. Hydrodynamical cosmological simulations suggest that Y SZ serves as one of the most robust total mass proxies for galaxy clusters, with a scatter of about 10% (e.g., Motl et al. 2005;Nagai 2006;Battaglia et al. 2012b;Kay et al. 2012;Krause et al. 2012;Yu et al. 2015). This has motivated construction of a lowscatter core-excised X-ray mass proxy, Y X ≡ M gas T X (Kravtsov et al. 2006). Despite the robustness of Y SZ as a mass proxy, there can be large deviations from self-similarity, particularly in extreme cases. Mergers between clusters of galaxies are the most energetic events in the present day Universe, and therefore not surprisingly have an effect on SZ observables (see Figure 6). The first and second core passages of a major merger event induce transient boosts in the tSZ signal via compression and shock heating of the cluster plasma, which are more pronounced for the maximum Comptonisation parameter y max than for its integrated value Y SZ (< R) (Motl et al. 2005;Poole et al. 2007;Wik et al. 2008). SZ observable-mass scaling relations involving the latter are therefore far less affected by mergers than the former or X-ray proxies such as T X and L X .
The scatter in the integrated SZ observable-mass relations originates from the non-thermal pressure provided by bulk and turbulent gas motions generated by mergers and mass accretion (Yu et al. 2015) and asphericity and substructures (Battaglia et al. 2012b) in the ICM, while the normalisation of the Y SZ − M relation is sensitive to the input cluster astrophysics, such as radiative cooling, star formation, and energy injection from stars and AGN feedback (e.g., Nagai 2006;Battaglia et al. 2012b;Kay et al. 2012).
Recent results for cosmological determinations using SZ-selected clusters and, specifically, their SZ signal as a mass proxy have found good agreement with cosmology inferred by other means, improving constraints on the dark energy equation of state and number of neutrino species, with the dominant systematic being the overall scaling of the SZ signal with cluster mass (e.g. Bocquet et al. 2015;Planck Collaboration et al. 2016c;de Haan et al. 2016;Hilton et al. 2018). Detailed resolved and stacked studies of the average pressure profiles and how Y SZ scales with mass therefore provide crucial tests of the simulations. The SZ effect in a cluster core which is undergoing a merger, from a simulation oriented so that the largest gas motions are predominantly within the line of sight. Certain X-ray quantities are also shown for comparison. Shown quantities are Compton y (top-left), y kSZ (top-middle), projected mass-weighted temperature, k B T mw (top-right), X-ray surface brightness, S X (bottom-left), electron optical depth τ e (bottommiddle), and projected X-ray temperature, k B T X (bottom-right). Note that while cold fronts are clearly seen in projected temperature, X-ray emission, and optical depth, they are essentially invisible in the tSZ signal due to the continuity of thermal pressure across the front. However, the bulk velocity of the gas across the fronts produces a kSZ signal.
Merging activity can induce offsets between the X-ray and tSZ peaks, since the former strongly tracks the maximum gas density (S X ∝ n e 2 ) and the latter is determined by the maximum integrated line-of-sight pressure (Molnar et al. 2012;Zhang et al. 2014), while the peaks in the X-ray and SZ surface brightnesses should be coincident for relaxed clusters. Comparison of observed offsets to simulations can be used to estimate merger parameters such as the mass ratio and relative velocity.
Cluster-cluster and cluster-group mergers often exhibit velocities of several thousand km s −1 during core passage, and can therefore also produce a strong kSZ signal in some circumstances (e.g. if close to the line of sight). If the velocities are large, the kSZ signal may even dominate over the tSZ signal during core passages (Ruan et al. 2013). Multi-frequency observations of these systems are essential in order to distinguish between the two effects. Note, however, that due to the linear dependence of the kSZ effect on velocity, projection effects can complicate the interpretation if multiple velocity components are contributing to the line-of-sight kSZ signal.
3.3 The complementarity of X-ray and SZ measurements X-ray and tSZ measurements are independent and highly complementary probes of the thermodynamics and kinematics of the ICM. Joint X-ray/tSZ studies, for example, allow one to reconstruct the thermodynamic state of the ICM by exploiting the different line of sight density dependences in the surface brightness integrals for each. The X-ray surface brightness S X , commonly expressed in cts arcmin −2 s −1 , is: where n e ( ) and T e ( ) are respectively the electron density and temperature along sight line , Λ ee (T e , Z) (in cts cm 5 s −1 ) is the X-ray emissivity measured by the instrument within the energy band used for the observation, z is the cluster's redshift, and Z is the metallicity. In units of erg cm 5 s −1 , there is an additional factor of (1+z) −1 , consistent with the standard (1 + z) −4 dependence of cosmological dimming. Thus, while the Comptonisation parameter (y ∝ n e T e d ; see Equation 5) is linear in density, X-ray surface brightness varies as the square (S X ∝ n e 2 d ). While X-ray emission suffers from cosmological dimming, it is useful to note that scaling relations in ΛCDM cosmology predict that, for a given mass defined with respect to overdensity, clusters at higher redshift are hotter, denser and therefore more X-ray luminous than their local counterparts. As a result, the observable X-ray flux (at fixed mass) may not decrease strongly at z 1, resembling a way that is similar to the SZ signal (e.g. Churazov et al. 2015).
3.3.1 ICM temperature from joint X-ray/tSZ studies SZ and X-ray observations of clusters can in principle be used in combination to derive "mass-weighted" temperatures by using the estimate of the gas pressure from the SZ and that of the gas density from the X-ray (e.g., Adam et al. 2017a). This technique offers the possibility to estimate the temperature, entropy and HSE mass profiles to high redshift (z 1), where X-ray spectroscopic measurements are very challenging due to very low photon counts, and hence characterise a redshift evolution of the ICM profiles throughout the epoch of cluster formation.
However, some caution should be exercised on interpreting the effective weights of temperatures derived from joint tSZ/X-ray analyses. First, the tSZ signal is proportional to the gas mass multiplied by the mass-weighted temperature, T mw , while the temperature of the hot ICM inferred by fitting the X-ray spectrum with a thermal emission model is a spectroscopic temperature T spec (Mazzotta et al. 2004;Vikhlinin 2006). Hydrodynamical simulations predict that there are discrepancies between T mw and T spec (Mathiesen & Evrard 2001;Nagai et al. 2007;Piffaretti & Valdarnini 2008;Rasia et al. 2014, see also Figure 6 for the difference between the T mw and T spec maps of one of the idealised cluster merger simulations).
In addition, it is well known that the density profile results can be biased in the presence of significant gas clumping, or projection from a triaxial ICM distribution when spherical models are assumed (see e.g. Mathiesen et al. 1999;De Filippis et al. 2005;Simionescu et al. 2011;Nagai & Lau 2011;Bonamente et al. 2012;Limousin et al. 2013;Vazza et al. 2013;Battaglia et al. 2015a;Eckert et al. 2015;Umetsu et al. 2015;Rossetti et al. 2016). If these density biases are not accounted for, then the joint tSZ/X-ray temperature estimate will not yield the sought-after T mw value. While it is possible in principle to use the combination of tSZ/X-ray data, including spectroscopic temperature estimates from high X-ray photon counts, to measure cluster triaxiality out to high redshifts (Sereno et al. 2012), once again the degeneracy between the triaxiality parameters and gas clumping/substructures, as well as the challenges of X-ray spectroscopy at high-z, might limit this method's applicability.
Observationally, this technique was pioneered in the first decade of this millennium using, for example, the combination of Nobeyama Telescope/SCUBA and Chandra data (Kitayama et al. 2004), in joint SZA + Chandra observations (Mroczkowski et al. 2009), and in APEX-SZ + XMM-Newton observations (Nord et al. 2009;Basu et al. 2010), although the statistical uncertainties were still very large. The combination of wide FoV and high sensitivity to the SZ signal brought about by Planck allowed a leap forward in the use of the joint tSZ/X-ray technique, as demonstrated in the analyses by Eckert et al. (2013a,b), who exploited ROSAT and Planck data to probe out to ≈ 1.5 × r 500 in a sample of nearby massive clusters. More recently, this effort has been extended using either the archival Planck or Bolocam datasets. These analyses include 1) a joint analysis of Bolocam and Chandra data for a large sample of 45 clusters probing the thermodynamics out to r 500 (Shitanishi et al. 2018), 2) a detailed analysis using a sub-sample of 6 clusters including data from Bolocam, Chandra, the Hubble Space Telescope (HST), and Hyper Suprime-Cam (HSC) lensing data (Siegel et al. 2018) to probe non-thermal pressure support out to r 500 , and 3) a large effort by the XMM Cluster Outskirts Project (XCOP) to study the thermodynamics, non-thermal pressure support, and outskirts (> r 500 ) of the 13 of the most significant Planck detections (see Eckert et al. 2017aEckert et al. , 2019Ettori et al. 2019;. Considerable effort is now on-going to extend these tools to high resolution (subarcminute) SZ samples. For instance, a legacy project using NIKA2 observations of a large sample (∼ 50) clusters at a resolution of ∼ 15 shows potential for NIKA2/XMM-Newton analyses. This potential is presented in Ruppin et al. (2018), following the developments made with the pathfinder camera NIKA (Adam et al. 2015(Adam et al. , 2016Ruppin et al. 2017), and the first two-dimensional temperature map reconstruction with this method was recently published by Adam et al. (2017a) in MACS J0717.5+3745 (z = 0.55). See Figure 7 for illustrations. The authors also directly compared the temperatures recovered through tSZ/X-ray imaging with the spectroscopic X-ray temperatures measured with Chandra and XMM-Newton. The SZ-derived temperature measurement is about 10% larger than the spectroscopic one from XMM-Newton, which is within the calibration uncertainties of both instruments, but may also include systematic effects driven by assumptions about the gas line-of-sight geometry and clumpiness.

ICM temperatures from the relativistic SZ effect
A more direct, spectroscopic SZ estimate of the ICM temperature is possible using multi-frequency tSZ measurements to separate the rSZ and classical tSZ contributions ( § 2.3). While theoretical progress has been made in computing the relativistic terms accurately (e.g., using SZpack, Chluba et al. 2012b), actual measurements of the cluster temperatures through the relativistic corrections have been more challenging. Only recently, using all-sky data from the Planck satellite, which covers the tSZ spectrum almost entirely, has it been possible to constrain the rSZ spectral distortion Fig. 7 Temperature measurement of the ICM using tSZ (NIKA) + X-ray (XMM-Newton) imaging and comparison to X-ray spectroscopic measurements. Left: SZ+X-ray deprojected temperature profile towards MACS J1423.8+2404 (red dashed region) and comparison to XMM-Newton (purple, filled dots) and Chandra (blue, open diamonds) spectroscopic measurements. Right: SZ+X-ray temperature map of the hot gas toward the galaxy cluster MACS J0717.5+3745 obtained from combining XMM-Newton density and NIKA pressure imaging. Figures from Adam et al. (2016Adam et al. ( , 2017a. in a large stacked sample of clusters (Hurier 2016a). With future ground-and spacebased instruments it is expected that measurement of the rSZ effect will turn into a robust technique for inferring cluster temperatures, and hence their masses, for both astrophysical and cosmological analyses.
The primary requirement for rSZ effect measurements is multi-frequency coverage of the tSZ spectrum spanning both the decrement and increment. Early attempts were made by combining data from several different experiments, e.g., by Hansen et al. (2002b) and Nord et al. (2009), but limited sensitivity did not break the degeneracy between the rSZ and velocity-induced kSZ contributions. Using data from the Z-spec grating spectrometer and Bolocam, Zemcov et al. (2012) could constrain the temperature of the hot cluster RX J1347.5-1145 to high accuracy, but only by neglecting the kSZ contribution. Similarly, a combination of X-ray and SZ measurements were used to quantify higher order rSZ contributions for the Bullet cluster (Prokhorov & Colafrancesco 2012;Chluba et al. 2013).
The situation improved considerably with data from the Planck satellite, which provided well-calibrated all-sky measurements for galaxy clusters in the frequency range 30 − 860 GHz. It is also possible to ignore the kSZ contribution by averaging (stacking) the data from several hundred clusters, as the kSZ signal can be both positive or negative due to the random nature of peculiar motions. This approach was explored by Hurier (2016a) and Erler et al. (2018) to constrain the average y-weighted temperature from cluster samples.
In Figure 8 the rSZ modelling results from Erler et al. (2018) are shown. The left panel is an illustration of the rSZ spectral distortions for a fixed value of Compton-y (10 −4 in this case) and the relative positions of the nine frequency bands of Planck. As can be seen, the entire tSZ spectrum is covered by Planck, although for better sensitivity and higher resolution only the 70 − 860 GHz data were used in Erler et al. (2018). In the right panel of Figure 8, the result of stacking images of 772 Planck clusters are shown. The impact of the kSZ effect is averaged out and the thermal spectrum is constrained to sufficient accuracy to make a measurement of the mean sample temperature to 4.4 +2.1 −2.0 keV. This temperature is found to be slightly lower than the mean mass-weighted X-ray spectroscopic temperature of the sample, (T X = 6.91 ± 0.07 keV), although the tension was only at the level of 1.3σ. This difference can potentially indicate a low level of gas clumping in galaxy clusters which causes the density-squared weighted X-ray temperature to stay above the rSZ-derived y-weighted temperature. Figure 8 also demonstrates the importance of modelling the cluster-centric far infra-red (FIR) emission simultaneously to obtain an accurate estimate of the rSZ signal.
Future high resolution, multi-frequency data can open the possibility of direct rSZ measurements of cluster temperature profiles or post-shock electron temperatures. This rSZ-derived mean projected temperature is very close to Compton y-weighted (Hansen 2004;Kay et al. 2012;Morandi et al. 2013;Erler et al. 2018). 12 Nonetheless, as with X-ray spectroscopic temperatures, rSZ-derived temperatures are not free of systematics. In the case of rSZ, these can be due to the effects of mixing from multiple temperature components along the line of sight, where first and higher moments of the temperature and density distribution can present a bias (Chluba et al. 2013). 12 Ideally one would measure mass-weighted temperature, T mw ≡ ∝ y/τ e , rather than one weighted by Compton-y

Mass Profiles
After deprojection of the X-ray and tSZ surface brightness profiles, the derived quantities can be combined through the ideal gas equation to recover the gas temperature kT e = P SZ /n X and the entropy K = P SZ n −5/3 X , where P SZ is the electron pressure inferred from the tSZ data and n X is the electron number density inferred from the X-ray data. It is then possible to reconstruct the gravitating mass profiles of the ICM under the assumption of hydrostatic equilibrium (HSE) by using the equation (e.g., Zaroubi et al. 1998;Puchwein & Bartelmann 2006;Ameglio et al. 2009): where ρ X = µ e m p n X is the gas density inferred from the X-ray data, and µ e ≈ 1.18 is the average weight per electron and m p is the proton rest mass. The hydrostatic mass is known to be slightly biased (see review on cluster mass estimation in these proceedings).
Over the years, high-quality radial (1-D) constraints on ICM properties have been obtained (e.g., Mroczkowski et al. 2009), the most recent examples of which combined either XMM-Newton or Chandra imaging (or both) and Planck or Bolocam SZ data for various cluster samples (e.g., Tchernin et al. 2016;Ghirardini et al. 2018a,b;Shitanishi et al. 2018). Given the linear dependence of the SZ signal on the gas density, measurements obtained through this technique extend farther in radius than typical spectroscopic X-ray measurements and allow one to probe the relatively unexplored outskirts of local clusters out to their virial radius (Ghirardini et al. 2018b). The wide radial range accessible in this study also brings high-quality measurements of hydrostatic mass profiles (typical statistical uncertainty ∼ 3% at r 500 ) to probe the internal structure of dark matter halos ). In the majority of cases, the underlying mass distribution was found to be well described by an NFW profile.

tSZ power spectrum constraints on ICM thermodynamics
Measurements of the tSZ angular power spectrum from wide-field surveys, as opposed to high-resolution measurements of the pressure fluctuations within selected, high-mass clusters (discussed in §4.5), present a complementary approach to constrain several important astrophysical parameters. Due to the redshift independence of the tSZ effect, this contribution to the CMB power spectrum comes from cluster (and galaxy group) ICM at all redshifts as well as from the unvirialised gas in cosmic filaments.
The original application for the tSZ power spectrum measurement was intended for cosmology, due to its strong dependence on several cosmological parameters such as Ω M , the matter density, or σ 8 , the amplitude of matter power spectral fluctuations on 8 Mpc h −1 scales (e.g., White et al. 1993;Komatsu & Seljak 2002). Originally predicted to be dominated by massive clusters at intermediate redshift (Cole & Kaiser 1988), the main contribution to the tSZ power at the angular scales relevant for ground-based CMB instruments (i.e., ∼ 2000 − 3000) in fact originates from low-mass groups and clusters, and it was realised that the uncertain astrophysics (e.g., non-thermal pressure support and energy feedback) of these low-mass systems at intermediate and high redshifts will make the goal of using the tSZ power spectrum for cosmology more difficult to attain (Shaw et al. 2010;Trac et al. 2011;Battaglia et al. 2012b). If, on the other hand, the cosmological parameters are well known from other probes and can be restricted a priori, then the tSZ power spectrum measurement can be used to constrain ICM astrophysics. It was also recently shown that relativistic temperature corrections affect the inference of cosmological parameters from the SZ power spectrum analysis. These are less important at high-(i.e., 2000) but are relevant on the scales probed by Planck (Remazeilles et al. 2019).
One of the main challenges to using the tSZ power spectrum to infer cluster astrophysics lies in the fact that, for arcminute resolution instruments, the tSZ power is buried under cosmic infrared background (CIB) fluctuations, especially at the smallest angular scales accessible to such low-resolution probes ( 1000, see e.g., Reichardt et al. 2012). Using the bandpower measurements from the SPT-SZ survey and a model for the CIB contribution, Ramos-Ceja et al. (2015) showed that the currentgeneration tSZ power spectrum measurements are already useful in constraining the outer slope of the ICM pressure profile, or departures from self-similar evolution, for fixed cosmology. However, the drawback is that several astrophysical effects can conspire to produce similar changes in the tSZ power such that it will be difficult to break the associated parameter degeneracies (e.g., between the amount of non-thermal pressure support and its redshift evolution). It is also expected that a better knowledge of ICM astrophysics and the tSZ-CIB cross-correlation will help to alleviate some of the current tensions between the Planck and SPT-SZ measurements (e.g., McCarthy et al. 2014, Dolag et al. 2016. Future high-precision, high-resolution measurements of the tSZ and CIB power down to smaller angular scales are therefore required. This will likely be one of the (many) science drivers for future ground-based CMB measurements with large-aperture telescopes, such as those described in §6.2.8 and §6.2.9.

Shock fronts
Shock fronts, both from mergers and cosmic accretion, should be readily observable in the tSZ at high spatial resolution due to their pressure jumps (Markevitch & Vikhlinin 2007;Molnar et al. 2009;Ruan et al. 2013). For an ideal gas with adiabatic index γ, the relation between the pressure in the pre-and post-shock regions can be derived from the standard Rankine-Hugoniot jump condition, and is given by (e.g., Landau & Lifshitz 1959) where M is the sonic Mach number of the shock (M ≡ /c s , where c s is the sound speed), and the instantaneous electron pressure in the post-shock (downstream) and pre-shock (upstream) regions as P post and P pre , respectively. The value of adiabatic index is typically assumed to be that of an ideal monoatomic or fully ionised gas, γ = 5/3, but this can vary if there is a significant energy density contribution from non-thermal particles or magnetic fields, or if there is strong magnetic field parallel to the shock front (Helfer 1953, Sarazin et al. 2016, Guo et al. 2018. Note that equation (15) implicitly assumes a single component fluid with equal electron and ion temperature (T e = T ion ), while the tSZ signal can only be used to infer the electron pressure. This immediately shows that if the post-shock thermalisation is not instantaneous, but instead for example follows Coulomb interactions, then M derived from tSZ measurements will give a lower estimate than from the density compression value, and a similar bias will also show up from the X-ray spectroscopic temperature ratio. Typically, one would need to analyse both X-ray and tSZ data with the same geometrical considerations to get a complete picture of the processes involved in thermalisation following a shock. It can be seen from Equation 15 that the pressure ratio scales roughly as M 2 , unlike the density compression ratio n post n pre , which saturates at a value of 4 at M 5. In addition, two advantages of the tSZ effect position it as very attractive for studying shocks in the ICM. First, the redshift independence of the SZ effect means shock-like features can be studied in clusters from arbitrarily high redshifts, unlike X-ray temperature measurements which has to contend with a rapidly diminishing number of photons. Second, the linear dependence on gas density for the tSZ signal amplitude means that it can be better suited to study merger or accretion shocks in the low density cluster outskirts. A third and practical advantage of tSZ observation compared to X-ray spectroscopy is that for a relatively strong shock in the hot ICM, the postshock temperature will reach a high value (in the range 20 − 40 keV), well beyond the spectral window of most X-ray satellites. This can make modelling of important physical processes like the thermal equilibration timescale very difficult. This is not a concern for tSZ data, particularly when it is modelled jointly with the X-ray density compression. These considerations also bring SZ-based shock measurements within observational reach for the cluster outskirts region, e.g. for shocks associated with radio relics (for a review of diffuse radio emission, see Feretti et al. 2012) or the tentative detection of an accretion shock near the cluster splashback radius (Hurier et al. 2017).
Realising these potentials requires good sensitivity and sufficiently high angular resolution, ideally well below the ∼1 beam which is the typical resolution for groundbased survey telescopes like SPT and ACT. An early example of tSZ imaging of putative shock-heated substructure in the ICM was provided by Kitayama et al. (2004) using NOBA data. The first sub-arcminute detection of a shock feature was provided by Korngut et al. (2011), which used the pathfinder MUSTANG-1 bolometer array on the 100-meter Green Bank Telescope (GBT) to image MACS J0744.8+3927. The SZ-identified feature was then used to guide the region selection for X-ray analysis using Chandra, which provided constraints on the shock Mach number. Illustrations from this measurement are shown in Figure 9. It should be noted that the limited FoV of the MUSTANG-1 camera meant that only the shocked region could be detected, while the rest of the ICM contribution is filtered by the instrumental response. There is a same problem with interferometric observations, such as those performed  Fig. 9 A comparison of MUSTANG measurements of the shock-heated gas in the galaxy clusters MACSJ0744.8+3927 and MACSJ0717.5+3745 with multi-wavelength data from the Chandra X-ray Observatory, the Very Large Array (VLA), and the Hubble Space Telescope (HST). Top-left: MUSTANG tSZ measurement at 90 GHz with effective resolution of 13.5 (due to smoothing). Top-middle: Chandra X-ray image of MACSJ0744.8+3927. Top-right: Chandra X-ray measurement in blue with the tSZ signal overplotted in magenta. Lower Left: Fitting for the shock Mach number in MACSJ0744 with three different methods. The fitting was performed on the Chandra X-ray data using regions selected by the MUSTANG imaging. See Korngut et al. (2011) for details. Lower-right: Multi-wavelength MUS-TANG+Chandra+VLA+HST images of MACSJ0717.5+3745. The background image is HST imaging as part of the Frontier Fields project. The red overlay on both panels is S-band (2-4 GHz) radio data taken with the VLA, showing the radio relic associated with re-accelerated AGN tails. The blue overlay is Chandra X-ray imaging. The peach-coloured overlay is the MUSTANG SZ decrement at > 3σ, where the dominant feature appears to be a pressure overdensity shaping the southern half of the radio relic. All figures apart from the lower right panel are from Korngut et al. (2011). Lower-right panel is courtesy Reinout van Weeren, using MUSTANG data presented in Mroczkowski et al. (2012) and Chandra + VLA data presented in van Weeren et al. (2017). with ALMA, as discussed below. This makes it impractical to image a large number of cluster targets to search for SZ shocks. Typically, X-ray observation are used to overcome the spatial filtering limitation or choose targets with known shock features. Low-resolution Planck all-sky survey data have also been used to study shocks, in the nearby Coma cluster (Planck Collaboration et al. 2013b), in combination with XMM-Newton. An alternative approach has been to target radio relics in merging clusters, which are elongated synchrotron emitting regions thought to be associated with shock  An ALMA tSZ effect measurement of a shock front coincident with a bright radio relic, in the massive El Gordo galaxy cluster at z = 0.87. Left: Multi-wavelength view of the shock, with background image from a Spitzer IRAC mosaic at 3.6 µm, and orange and green contours marking the Chandra Xray 0.5 − 2 keV X-ray emission (see also Botteon et al. 2016) and the 2.1 GHz radio emission (Lindner et al. 2014), respectively. The observed SZ intensity distribution from an ALMA deconvolved image is shown in the zoom-out inset. Right: The shock front as seen from both X-ray and tSZ measurements. The upper panel shows the temperature, where the very faint pre-shock region in this high-z cluster do not provide enough photons to constrain the X-ray spectroscopic temperature. In the middle panel the surface brightness jump in X-ray clearly reveals a shock (or a cold front). Blue-dotted line is the mean X-ray background. In the bottom panel the ALMA SZ data indicates a modulated flux (red dot-dashed line) which is the result of interferometric imaging of a step function-like shock signal (shown with the green solid line for the best-fit model). Figures adapted from Basu et al. (2016). fronts. The first tSZ shock measurement on a radio relic was provided by Erler et al. (2015), also for the Coma cluster.
Perhaps the most promising instrument for detailed tSZ studies of known shock fronts in clusters is the ALMA interferometric array. Despite its very limited FoV (< 1 at 90 GHz; see §5.1), ALMA provides unprecedented sensitivity and resolution ( 3 at 90 GHz) for accurate imaging of shock fronts. With future low-frequency additions to the ALMA suite of receivers (ALMA Bands 1 & 2 in particular; see § 6.1.1) the problem of limited FoV will be partially mitigated. This field is still in its early years, with the first ALMA observation of a merger shock being published by Basu et al. (2016) for the notorious galaxy cluster 'El Gordo' at z = 0.87 ( Figure 10). This high-z target demonstrates several of the advantages mentioned earlier, namely a strong complementarity with X-ray data (here X-ray photon count is too low to make a measurement of the pre-shock temperature) and a measurement of shocks in the cluster outskirts at high-z. The inset image in the left panel of Figure 10 shows an image-equivalent ("CLEANed image") of the ALMA interferometric data, where a wave-like intensity pattern indicates the presence of an underlying tSZ shock that roughly looks like a step function. The accepted approach for constraining a shock (or any other) model in such cases in not via the reconstructed images, which have biased noise properties (not uncorrelated), but rather by using a direct Bayesian analysis of the visibility data (or uv-data). This can also benefit from the combination with Xray data to constrain the large scales inaccessible to ALMA. The salient results of the shock analysis from Basu et al. (2016) are reproduced in the right-hand panels of Fig.10. The remarkable consistency of the shock location from tSZ and X-ray measurements, and the quality of the ALMA observation that took only a fraction of the time spent with Chandra for the same object, signify a coming-of-age of tSZ measurements for such astrophysical applications. Additional ALMA observations of shocks in other galaxy clusters have been made, including the well-known bow shock of the Bullet cluster, and the number of results will grow in the coming years. Using multi-frequency ALMA observations, or by combining data from several instruments, an SZ-only decomposition of cluster shock fronts could soon become a realistic possibility. For example, the kSZ effect will reveal directly the velocity discontinuity at the shock front, although projection effects can limit its observability, since most known shocks are seen edge-on. One interesting application for future multi-frequency SZ measurements will be to determine post-shock temperatures via the rSZ effect, since the high temperatures involved will guarantee a large departure from the classical tSZ signal, potentially producing a detectable ntSZ component. This will provide a temperature estimate that is completely independent of X-ray data (tSZ modelling of temperature requires the X-ray density information). Future large-aperture single-dish telescopes like the AtLAST project ( §6.2.8) promise revolutionary new developments in this field.

Cold fronts
The tSZ view of cold fronts, which are contact discontinuities frequently seen in Xray observations, is quite different from that of shocks ( §4.1). Cold fronts manifest as sharp jumps in X-ray surface brightness and temperature, but in opposite directions that roughly maintain thermal pressure balance. Simulations indicate that such contact discontinuities form through the subsonic bulk motions of gas. Due to the fact that cold fronts as seen in simulations are pressure-continuous, they should be invisible in the tSZ if the pressure is predominately thermal. There are two exceptions to this: a) if the ram pressure of the surrounding medium is large enough or b) if there is a strong magnetic draping layer either directly above or below the front surface. Regarding the latter possibility, these layers may form due to shear amplification of the weak ambient field during cold front formation and evolution. Such layers should be observable as a thin deficit in the tSZ signal in the vicinity of the front as seen in X-rays, but a high angular-resolution instrument will be required to detect them. In general, the thermal pressure gradient across cold fronts is not continuous in simulations, and therefore should be detectable in the tSZ. Also, a cold front with a velocity field at least partially oriented along our line of sight is in principle detectable as a contribution to the kSZ signal. The kSZ effect can also reveal the tangential gas motion across a cold front, where the thermal pressure itself has no discontinuity (i.e., interferometric tSZ measurements will not see a cold front). At the smallest scales the motions of electrons and ions in the weakly magnetised ICM are inherently anisotropic due to the fact that the Larmor radii of the electrons and ions are many orders of magnitude smaller than their respective mean free paths. This can manifest itself as a pressure anisotropy due to evolving magnetic fields via adiabatic invariance and other processes. Khabibullin et al. (2018) showed that such a pressure anisotropy could produce polarisation in the tSZ signal due to the different thermal velocities of the electrons in the directions perpendicular and parallel to the local field line geometry, which they estimated to be produced at a level of ∼10 nK. Such polarisation might be prominent near cold front surfaces where the magnetic field is increasing rapidly due to shear amplification which would in turn drive a local pressure anisotropy.

Probing pressure substructures with tSZ effect
Thanks to the advances of high angular resolution imaging experiments, it is possible to detect and characterise pressure substructures in the ICM through SZ studies. Substructures detected through the tSZ effect are related to gas compression driven by merger events associated with infalling substructures in the ICM. In combination with X-ray observations, tSZ substructure characterisation has proven useful in furthering our understanding of the underlying ICM physics. For instance, SZ observations from the DIABOLO (Pointecouteau et al. 2001) and NOBA (Komatsu et al. 2001;Kitayama et al. 2004) instruments, respectively installed on the IRAM-30m and Nobeyama-45m telescopes, detected an offset between the dense X-ray core and the pressure peak in the cluster RX J1347.5-1145 (z = 0.45), which at that time had been known to be a typical dynamically relaxed cool-core cluster. These observations indicate that RX J1347.5-1145 is an ongoing merger, and this cluster was later studied by many targeted SZ experiments, such as MUSTANG-1 (Mason et al. 2010), CARMA (Plagge et al. 2013), NIKA (Adam et al. 2014), BOLOCAM (Sayers et al. 2016b) and ALMA (Kitayama et al. 2016). For illustration, the ALMA image, at a resolution of 5 , is presented in Figure 11 and compared to X-ray and strong lensing data. These data were used to characterise the cluster core and the overpressure caused by the merging of a sub-cluster in RX J1347.5-1145 (Ueda et al. 2018). In combination with X-ray measurements, these SZ data provided the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. Additionally, a comparison of the offsets between the peaks in the SZ and strong lensing signals was used to infer the self-interaction cross-section of dark matter to be less than 3.7 h −1 cm 2 g −1 .
In addition to the emblematic case of RX J1347.5-1145, many other clusters were observed to study substructures using the tSZ effect. MUSTANG observations of the massive Frontier Fields cluster MACSJ0717.5+3745 (z = 0.55), which contains at least 4 merger components (Ma et al. 2009;van Weeren et al. 2017), revealed a high pressure region, presumably resulting from the merger of the 2 most massive components and likely driving and reshaping the southerly radio jet (see the lower-right panel of Figure 9). The shocked overpressure region is confirmed in X-ray observations, which indicate the presence of hot (> 18 keV) gas. An additional SZ feature is seen to the west, in a region of high X-ray surface brightness and lower temperature, indicating the location of an intact cool core. Ground-based observations using on MUSTANG, NIKA, and NIKA2 have also been instrumental in probing ICM substructures at higher redshifts, with studies ranging from identification of substructures, the impact of substructure on cluster pressure profiles, mass estimates and global morphology, and the study of shocks (see for instance, Korngut et al. 2011;Mroczkowski et al. 2012;Adam et al. 2015;Young et al. 2015;Adam et al. 2016;Ruppin et al. 2017Ruppin et al. , 2018, as well as the discussion in § 3.1 and § 4.1). Notably, the maturity reached in last decade by high spatial resolution SZ imaging has enabled advanced image filtering techniques typically reserved for X-ray to be applied to SZ observations (Adam et al. 2018a). As SZ imaging instrumentation improves and the higher signal to noise images are obtained, we can expect that many of the imaging analysis techniques developed for X-ray studies will be applied to SZ imaging as well.

Probing ICM substructures with resolved kSZ
High-spatial resolution maps of the kSZ signal can in principle be used to measure the internal gas motions in the ICM, including bulk and random gas motions Nagai et al. 2003) and rotational motions (Chluba 2001;Cooray & Chen 2002;Chluba & Mannheim 2002;Baldi et al. 2017Baldi et al. , 2018. Such maps will provide a complementary measure of gas motions to those based on the X-ray surface brightness or pressure fluctuation analyses discussed in § 4.5, as well as those which will eventually be produced by microcalorimeters on future X-ray missions like the X-ray Imaging and Spectroscopy Mission (XRISM; Tashiro et al. (2018)), Athena (Nandra et al. 2013), or Lynx (Gaskin et al. 2016), which will be capable of measuring the shift and width of velocity-broadened spectral lines (Inogamov & Sunyaev 2003). The different density dependence of the two methods could potentially disentangle the effect of density variations along the line of sight from the measurement of the bulk velocity (see Figure 6). While the first detection of the kSZ signal was performed in a statistical manner, by measuring the mean pairwise momentum of optically identified clusters (Hand et al. 2012), the first kSZ detections in a single cluster were performed by resolving the signal towards MACS J0717.5+3745 (Mroczkowski et al. 2012;Sayers et al. 2013a;Adam et al. 2017b). At a redshift of z = 0.55, MACS J0717.5+3745 is a strik-ing example of a massive merging system and thus is an excellent target within which to search for the kSZ signal. It hosts at least four sub-groups known to present exceptionally large relative line-of-sight velocities from optical measurements, as first noted by Ma et al. (2009). Given the large merger velocity (β ≈ 0.01) and the compact nature of the source (a few arcminutes), it is not surprising this has served as the best known source where non-zero velocities can be inferred from the kSZ signal. Indeed, despite several attempts at detecting the kSZ signature in other sources (e.g., Holzapfel et al. 1997;Benson et al. 2003;Sayers et al. 2016bSayers et al. , 2018, the kSZ effect signal remains difficult to extract in individual clusters since it follows the same spectral dependence as the primordial CMB. It can only be identified using the two signals' differing spatial distributions, requiring both high sensitivity and high angular resolution observations. Additionally, as the kSZ signal is fainter than the tSZ signal (see Figure 2), it is more prone to contamination from dust and radio emission (see §5.3 and the discussion in e.g. Sayers et al. 2018). Figure 12 shows the Bolocam (Sayers et al. 2013a) and NIKA (Adam et al. 2017b) measurements of MACS J0717.5+3745. In both cases, these analyses took advantage of the dual-band capabilities of the instruments to separate the tSZ from the kSZ contributions, along with ancillary data from e.g. Herschel. Thanks to the higher angular resolution of NIKA with respect to Bolocam (12 and 18 versus 31 and 58 ), they were able to obtain a resolved map of the kSZ signal from the dual-band measurement. As noted above, Herschel observations were essential to account for the contribution from unresolved dusty galaxies, while Chandra and XMM-Newton X-ray data were crucial to disentangle the density, temperature, and pressure, in order to extract the gas line-of-sight peculiar velocity, under various modelling assumptions. The main bullet-like core of MACS J0717.5+3745 was found to move at velocity of ∼ 3000 km s −1 with respect to the CMB reference frame, while no overall net velocity was detected for the whole cluster. Such a high velocity remains consistent with standard ΛCDM cosmology.
Although still very futuristic, by also adding polarisation information through the pSZ (see § 2.5) one could in principle reconstruct the full 3D bulk and rotational velocity vectors, thus directly probing the large-scale motions of ICM structures. Such measurements will require improvements in sensitivity by several orders of magnitude, as well as better spectral coverage for separating and removing contaminating signals.

ICM turbulence from pressure fluctuations
Azimuthally-averaged radial profiles of thermodynamic quantities -the gas pressure, density, temperature or entropy are important proxies for global properties of a cluster, such as its mass. On the other hand, a great deal of information on the ICM dynamics is encoded in the velocity field and in small-scale fluctuations of thermodynamic quantities. Extensive measurements of the former will become available with future X-ray missions like XRISM, Athena, and Lynx and/or future kSZ imaging, while the latter can already be measured with the currently operating X-ray observatories and SZ facilities. Cosmological simulations predict that even in relaxed clusters the ICM is not at rest, and the amplitudes of the velocity/density/pressure variations are sensitive to many currently unknown characteristics of the ICM.
The first attempt to measure pressure fluctuations in the Coma cluster was done in Schuecker et al. (2004), where X-ray images and temperature maps obtained with XMM-Newton were combined to get a map of pressure fluctuations. A notable observational complication of this approach is that the accurate temperature measurements require a large number of X-ray photon counts, implying that only low angular resolution maps can be constructed. This reduces the spatial dynamic range can be probed. It was later realised that it is much easier to deal with the X-ray images themselves, rather than with the pressure maps, since the X-ray emissivity below ∼ 3 keV is a direct proxy for the density fluctuations. These density fluctuations could be related to the pressure variations in the turbulent eddies, variations of the gravitational potential, sound waves, bubbles of relativistic plasma, or entropy variations due to the gas displaced from its equilibrium position (Churazov et al. 2012). Numerical experiments with driven solenoidal turbulence and with the cosmological simulations have shown that for the large-scale subsonic motions, the latter mechanism (entropy perturbations) is likely the dominant source of density variations, which can be related to the scale-dependent amplitude of gas velocities (e.g., Gaspari & Churazov 2013;Gaspari et al. 2014;Zhuravleva et al. 2014a). It was also found that for the cores of relaxed clusters, the isobaric nature of the dominant perturbations can be confirmed observationally Churazov et al. 2016;Zhuravleva et al. 2016). This (and similar) approaches have been used to estimate the gas density and velocity power spectra in a number of clusters, with the freedom in choosing the underlying "unperturbed" density distribution models being the one of the main sources of uncertainty (e.g., Zhuravleva et al. 2014b;Walker et al. 2015;Eckert et al. 2017b;Zhuravleva et al. 2018).
Unlike the X-ray data, tSZ observations are sensitive only to the electron pressure variations. Therefore, SZ data provide a highly complementary (to X-rays) view of the ICM perturbations, "cleaned" from the pure isobaric (entropy) perturbations. Given that in cool cores the isobaric perturbations seem to dominate, we can expect such cores to appear much smoother in deep SZ imaging than in X-rays, although the bubbles of relativistic plasma (see §4.6) and sound waves generated by the AGN (e.g., Fabian et al. 2006) could still contribute. Based on cosmological simulations (e.g., Lau et al. 2009;Zhuravleva et al. 2013;Nelson et al. 2014), one can expect on average ≈ 10% pressure fluctuations close to but outside the cores of relaxed systems (≈ 20% in disturbed systems), progressively increasing toward the outskirts region. Therefore, the role of non-isobaric modes, involving pressure variations, might become more important than the entropy perturbations at large distances from a cluster core. This is especially true for merging systems, where shocks cause large variations in the thermal pressure.
The first exploratory study of SZ fluctuations (Khatri & Gaspari 2016, see Figure 13) was based on the Coma cluster using Planck data. The Coma cluster has an angular size of a few degrees, close to the peak anisotropy in the CMB power spectrum. Thus the analysis necessitates Planck multi-frequency coverage, in order to separate SZ fluctuations from the CMB primary anisotropies and foregrounds. At the same time, the effective Planck angular resolution is ∼ 5 − 10 , which limits the available dynamic range. The characteristic amplitude of the pressure variations obtained from the Planck data is ∼ 30% on scales ∼ 15 (∼ 400 kpc). Khatri & Gaspari (2016) assumed that these pressure variations and the density fluctuations on smaller scales derived from X-ray images (Churazov et al. 2012) are part of the same spectrum, and concluded that pressure variations dominate in Coma. It should be noted, that similar to the X-ray analysis, the choice of the "unperturbed" global model can affect the results. This is because one must divide by the fitted bulk model to estimate the relative fluctuations, (see e.g. Bonafede et al. 2018;Di Mascolo et al. 2018). In addition, residuals from primary CMB and foreground fluctuations provide an additional source of uncertainty, which -as was done in Khatri & Gaspari (2016) -must be carefully accounted for even when multi-wavelength data are available.
A handy way to combine SZ and X-ray data would be by cross-correlating the deviations seen in SZ and X-ray images relative to the same underlying model used to describe the global distributions of the X-ray and SZ signals (see §4.5 in Churazov et al. 2016). This can be used to differentiate between fluctuations having a different nature, e.g., adiabatic versus isobaric fluctuations. Such an approach would work best if the SZ and X-ray images have comparable angular resolutions and comparable S/N per resolution element. For nearby clusters like Coma, the ability to recover large (∼ 1 degree) scales is also required, which is challenging for ground-based SZ instruments with few arcminute FoVs (see § 5.1 and § 5.2) and limited frequency coverage, which impacts component separation and the removal of astrophysical contamination (see 5.3). On the X-ray side, such widefield data could be provided by forthcoming SRG/eRosita, exploiting its large FoV (Predehl et al. 2016;Predehl 2017). Both Xrays and SZ data suffer from projection effects, which limit the dynamic range of the spatial scales that can be probed. Indeed, when the size of the perturbation is much smaller than the characteristic size of the cluster, the amplitude of the perturbation is attenuated by averaging along the line of sight. The pressure distribution (SZ) is obviously more extended than the distribution of the density squared (X-rays), implying that SZ data suffer from projection even more than X-rays. As we elaborate in our conclusions ( §7), future facilities will have the opportunity to address many of these limitations. Despite the current limitations, it is clear that studying pressure fluctuations using SZ data alone or in combination with X-ray data is an extremely powerful way of studying the dynamical state of the ICM and its physical properties.

The physics of AGN feedback
AGN outbursts are the dominant feedback mechanism preventing runaway cooling in galaxy clusters (Gaspari et al. 2011(Gaspari et al. , 2012Li et al. 2015;Yang et al. 2018), and should produce different tSZ and ntSZ effect signatures, both due to the shocks driven by the outbursts and the different possible compositions of the plasma filling the bubbles. Within cool-core galaxy clusters, AGN outbursts with powers up to 10 46 erg s −1 and enthalpies of 10 58 − 10 62 erg are observed to produce giant cavities (see e.g. Blanton et al. 2010;McNamara & Nulsen 2012;Fabian 2012, for reviews). These outburst driven cavities are seen as low-surface brightness patches in X-ray observations, which are filled with radio emission at low frequencies ( 10 GHz).
In Figure 14, we show the results of magneto-hydrodynamical simulations, reported in Ehlert et al. (2018a), of AGN outbursts including cosmic ray support. These cavities drive shocks, and are thought to be supported by magnetic fields, ultra-hot thermal gas, and cosmic rays, although the relative proportions of these sources of pressure are still highly unknown. By exploiting the unique spectral signatures of various SZ effects discussed in § 2, frequency-resolved SZ observations have the unique potential to constrain the nature of the pressure support in cavities (see e.g., Pfrommer et al. 2005;. For example, as discussed in Pfrommer et al. (2005): 1. If magnetic fields provide significant support to the cavity, it should be close to equipartition with a relativistic pressure component, and the SZ decrement will have < 1/8 the strength of standard nonrelativistic thermal gas. 2. If very hot ( 50 keV) gas provides most of the pressure, the strength of the tSZ effect decrement will be suppressed by more than 25% due to relativistic corrections at fixed y-parameter. 3. If non-thermal relativistic cosmic ray protons or electrons dominate the pressure support, the tSZ effect signal from a bubble will be negligible (see Figure 15).
In most such systems, the simplest assumption of equipartition between relativistic and magnetic support is unlikely, as the internal pressure falls far short of the work required to carve out the cavities by more than an order of magnitude. Recently, Abdulla et al. (2018) reported on a 30 GHz observation of the bubbles in MS0735.6+7421, the highest energy pair of bubbles currently known (McNamara & Nulsen 2007;Vantyghem et al. 2014), using the 23-element CARMA interferometer (see Figure 15). Using a model based on the X-ray results in Vantyghem et al. (2014), they find the joint tSZ signature from the bubble pair is consistent with zero at ∼ 5σ. The right panel in Figure 15 shows this result quantitatively. The lack of tSZ signal in the cavities is modelled in terms of a "cavity suppression factor", which at a value of unity will imply there is minimal tSZ contribution and the cavities are pressure supported by either non-thermal relativistic electrons, thermal electrons from a very hot gas (T e 1000 keV), or magnetic fields. This cavity suppression factor is obtained directly from an MCMC analysis of the CARMA visibility data, assuming a specific cluster geometry, resulting in the two histograms shown in the figure for two different line-of-sight orientations. The level of cavity suppression that can be achieved by a thermal plasma in pressure equilibrium with the surroundings and having a specific temperature, or by a non-thermal electron distribution with certain minimum momentum (here p = βγ, where β = /c and γ = 1/ 1 − β 2 = E e /m e c 2 ), are marked by the different solid and dashed lines, respectively.
Tentatively, the Abdulla et al. (2018) result can be considered as an argument for non-thermal pressure support within AGN cavities. This is because it is difficult to make the case for a 1000 keV thermal plasma filling the cavities since it is well beyond the pair-production energy of electrons. In other words, the CARMA observations have provided indirect evidence for cosmic ray transport by AGN bubbles. The streaming of cosmic ray protons from a central AGN into the dense ICM cores has been theorised to be a gentle source of heating that can solve the cooling flow problem (e.g., Guo et al. 2008, Fujita & Ohira 2012, Pfrommer 2013, Wiener et al. 2013. The transport of cosmic rays by AGN bubbles can aid this streaming process as the bubbles rise through the ambient medium at the sound speed, overcoming the limitation posed by Alfvèn waves (Ruszkowski et al. 2017, Ehlert et al. 2018b).
Clearly, data from single-frequency observations cannot be used to fully model the non-thermal SZ spectrum (see §2.4), as more frequencies are necessary to separate the tSZ, kSZ, and rSZ effect contributions. Future observations combining multi- Fig. 16 Main: Simulated Compton y vs. angular radius for a 10 12 M elliptical galaxy at z = 1. The peaked inner curve corresponds to a model where the halo gas is gravitationally-heated ,and has an NFW profile consistent with gas at the virial temperature. The gas is in hydrostatic balance in a host halo with a total mass 5 × 10 13 M (see Marconi & Hunt 2003). The outer component corresponds to a model in which the halo gas is heated by energy injection at the centre, and has a Sedov profile (i.e. a Sedov-Taylor blast wave; see e.g. Sedov 1946Sedov , 1959. The Sedov profile has an energy input of 3 × 10 61 ergs, which is consistent with typical theoretical feedback estimates for such galaxies (e.g. Scannapieco & Oh 2004). The model is observed at a time of 1 Gyr after the initial energy input. The errors are the forecasted noise when combining 2000 galaxies in the TolTEC (see §6.2.1) wide field survey. The dashed vertical lines correspond to the HWHM beam widths of TolTEC and SPT at 2.1 mm. Inset: Co-added 2.1 mm image of 3394 massive, quiescent 0.5 < z < 1 galaxies from SPT at 150 GHz. Although the angular resolution of SPT is ∼1 and the depth at 2.1 mm is approximately 17 µK-arcmin, co-adding nevertheless allows for a 2σ measurement of the SZ effect from the CGM around galaxies outside of clusters. The image is 8 × 8 , and the black circle represents a 1 radial aperture. Main figure from Scannapieco & Mauskopf in prep. Inset image from Spacek et al. 2016. frequency data with few arcsecond angular resolution can be expected to provide this much sought-after proof of cosmic ray heating in the ICM.

Large-scale impact of AGN feedback with stacked SZ measurements
The SZ effect has the potential to constrain the large-scale impact of feedback from AGN. As introduced in §4.6, feedback processes are thought to play a key role in suppressing the growth of new star-forming disks around isolated massive galaxies (e.g., Scannapieco & Oh 2004;Granato et al. 2004;Croton et al. 2006;Thacker et al. 2006;Di Matteo et al. 2008) and preventing the formation of large amounts of ≤ 1 keV gas in cool-core galaxy clusters (e.g., McNamara & Nulsen 2012;Fabian 2012).
In many theoretical models, energetic AGN outflows due to broad absorptionline winds or radio jets heat the circumgalactic medium and suppress the cooling needed to form further generations of stars and AGN. This suppression is redshift dependent, as the higher-redshift circumgalactic medium is more dense and rapidly radiating. Therefore, at high redshift, a more energetic outflow driven by a larger AGN is required for quenching to be effective. This provides an explanation for the large observed drop in the characteristic mass of star-forming galaxies since z=2 (e.g., Arnouts et al. 2005;Treu et al. 2005;Schaye et al. 2010). However, the details of such feedback remain extremely uncertain, and direct measurements of AGN kinetic energy input around isolated galaxies are extremely difficult (e.g., de Kool et al. 2001;Chartas et al. 2007;Borguet et al. 2013).
Measurements of the tSZ effect offer a method to overcome these difficulties. From Equation (5), and as noted previously, we see that the tSZ signal integrated over the area of sky around any object of interest is a calorimetric measure of electron thermal energy. This means that one can potentially use the tSZ effect to measure the total thermal energy output from AGN into the ICM. Searching for this signal was first suggested by Natarajan & Sigurdsson (1999), who roughly estimated the SZ effect from AGN and argued that it would be detectable even outside of galaxy clusters. More detailed calculations were carried out in e.g. Yamada et al. (1999), Yamada & Fujita (2001), and Chatterjee & Kosowsky (2007), including estimates of the redshift dependence of the signal and its global impact on the CMB. Such global calculations were refined in Scannapieco et al. (2008) and Chatterjee et al. (2008), who used large-scale cosmological simulations to compute the impact of AGN on the CMB power spectrum, and showed that SZ measurements would be able to constrain AGN feedback by appropriately stacking measurements around quasars and earlytype galaxies. In Figure 16, we show the results of a simple estimate of the difference in the radial SZ profile around elliptical galaxies. The predictions show that the signal is expected to be within reach of current and near future instruments (specifically, the figure compares the predictions for TolTEC, discussed in §6.2.1).
Constraints on the widespread impact of AGN through galaxy stacking have recently been obtained by several groups. Chatterjee et al. (2010) found a tentative detection of quasar feedback using the Sloan Digital Sky Survey (SDSS) to select galaxies for use in stacking Wilkinson Microwave Anisotropy Probe (WMAP) data, though the significance of AGN feedback in their measurements is disputed (Ruan et al. 2015). Hand et al. (2011) stacked ACT measurements around 2300 SDSSselected luminous red galaxies and found a ≈ 3σ detection. Planck Collaboration et al. (2013d) investigated the relationship between tSZ signal and stellar mass with significant results above stellar masses ≈ 10 11 M (see also Greco et al. 2015). Gralla et al. (2014) stacked ACT data to make a 5σ detection of the thermal SZ signal from radio galaxies. Ruan et al. (2015) also stacked Planck tSZ maps centred on the locations of 26,686 SDSS quasars and estimated the mean thermal energies in gas surrounding them to be ≥ 10 62 ergs (for a counterexample, see Cen & Safarzadeh 2015). Similarly, Crichton et al. (2016) stacked > 17,000 radio-quiet quasars from the SDSS in ACT data and found 3σ evidence for the presence of associated thermalised gas.
The inset panel of Figure 16 illustrates the results of stacking 150 GHz data around 3394 large passive galaxies with an average stellar mass of 1.5 × 10 11 M located outside of large galaxy clusters in the data publicly released by the SPT collaboration (Spacek et al. 2016). Note that at 150 GHz the SZ signal appears as a decrement, and this distinguishes it from all emission processes that can contaminate the signal. The signal was detected with a noise level of 0.7µK at 150 GHz. These measurements, along with similar stacked measurements from ACT (Spacek

Probing the missing baryons with stacked measurements of the SZ effect
A cosmic web of filamentary structures connecting clusters, groups, and galaxies is an essential prediction of nearly all large-volume cosmological simulations, and is expected to be found at temperatures of 10 5 − 10 7 K (Cen & Ostriker 1999;Shull et al. 2012). The emission is thus expected to peak at energies generally too soft to be observed in the X-ray, and the associated densities indicate these regions should have a very low X-ray surface brightness (Bregman et al. 2009). At the same time, this gas is too warm to be easily detected through absorption along the line of sight to AGN (Nicastro et al. 2008), although most searches so far have focused on this approach. Another method is to cross-correlate measurements of the tSZ and kSZ effect with other tracers of large-scale structure such as lensing maps (Van Waerbeke et al. 2014;Battaglia et al. 2015b;Atrio-Barandela & Mücket 2017), galaxy surveys (Hernández-Monteagudo et al. 2015), or the dispersion measures of fast radio bursts (Muñoz & Loeb 2018). Stacked measurements not only improve the signal to noise by increasing the effective integration time, but can also reduce certain forms of astrophysical contamination, as discussed in §5.3.
Recent stacking analyses of Planck data have provided the first putative direct detection of this material using the the tSZ effect (de Graaff et al. 2017;Tanimura et al. 2019). In both works, large catalogues (∼ 10 6 ) of galaxies were used to identify pairs of galaxies separated by ∼ 10 Mpc on the sky. These pairs were then scaled and stacked symmetrically, and models for the galactic SZ signatures were fit to the signal perpendicular to the filament direction and subtracted as shown in Figure 17. The results of both analyses showed a ∼ 5σ residual between the pairs. If a temperature T e ∼ 10 6 K is assumed, consistent with simulations, a significant fraction of the 'missing' baryons are recovered when accounting for the filling fraction of this cosmic web. We caution that the temperatures in simulations range by 2 orders of magnitude, however, and that the result for the baryon fraction recovered is linearly sensitive to the assumed temperature. A similar stacking analysis was also used in Tanimura et al. (2018) to reveal filamentary structures in supercluster environments.
The Planck all-sky Compton y map has a 10 angular resolution, which dilutes the signal over a larger area than it likely subtends. Future higher resolution SZ/CMB observations will suffer less from beam dilution, improving these constraints. Furthermore, this probe is expected to yield more significant results from stacking and cross-correlation studies, as well as through the analysis of the tSZ power spectrum at low multipoles ( < 100). By stacking in different mass-bins, one may furthermore be able to calibrate the temperature-mass relation using the rSZ (e.g., Erler et al. 2018), thereby providing another handle on the hydrostatic mass bias.

Observational considerations
The SZ effect is an attractive probe of the hot ICM for three main reasons: it has redshift independent surface brightness, observations can be performed from the ground, and the volumetrically-integrated tSZ signature is proportional to thermal energy, making it a useful probe for ICM calorimetry or detection of a cluster with an observable that tracks mass well. With the sensitive new instruments now online or planned for the near future (see §6), the integration times often rival those of X-ray observations, particularly for high-z systems and in the case where X-ray spectroscopy is necessary to perform calorimetric determinations of the cluster gas. However, there are several practical matters to take into account when studying cluster astrophysics using SZ observations. Here we briefly discuss interferometric and single-dish photometric observations, and contamination and confusion from mm/submm astrophysical sources.

Interferometric measurements
As noted in the introduction, interferometric arrays have a long history of use in detecting and imaging the tSZ effect. In contrast to single dish telescopes, interferometers have the appeal that they provide high resolution measurements without the requirement of building large aperture elements (see §5.2), and that the uncorrelated atmospheric signal is naturally filtered from the observation. This comes with the expense that the interferometric array behaves as a frequency-dependent spatial filter, and rarely attains the mapping speed of a single dish.
The point source sensitivity of an interferometer with identical array elements scales as (see e.g. Thompson et al. 2017;Condon & Ransom 2016): where N bl = M(M − 1) is the number of baselines in the observation, A e is the effective area of the array elements (i.e., optical efficiency η times collecting area A c ), T sys is the total system noise temperature including atmospheric and instrumental contributions, ∆ν is the bandwidth, and t is the integration time. Thus large gains for continuum signals like the SZ effect can be made by increasing the number of interferometric array elements, reducing the instrumental noise, and increasing the instrument bandwidth. The most fundamental property of an interferometric observation of the SZ effect (or any extended source) is that it samples the Fourier transform of a portion of the sky multiplied by the response pattern of the antennas (see e.g. Thompson et al. 2017;Condon & Ransom 2016). For an array with homogeneous elements, this response pattern is simply a Gaussian with size equal to the diffraction limit of each antenna. Thus, for an interferometer, the primary beam determines the usable FoV of an observation, rather than the angular resolution of the measurement. This, in turn, is determined by the sampling in 'uv-space' (so-called for the Fourier space coordinates u and v). The average response pattern is called the synthesised beam (often approximated as an elliptical Gaussian), which depends both on the uv-coverage and the weighting scheme used for creating an image from the visibilities. It is therefore cleanest to model the data in uv-space, as discussed for the interferometric data presented in, e.g., § 3.1 and § 4.1, In Figure 18, we show idealised, noiseless mock ACA and ALMA Band 3 and 6 observations (centred at 90 and 220 GHz respectively). These show the ability of ALMA/ACA to probe SZ structure at high resolution, but at the expense of having no sensitivity to scales larger than the full width half maximum width of the primary beam. For ALMA, this is 56 /(ν/90 GHz), much smaller than the typical extents of a massive galaxy cluster at any redshift. We discuss the ALMA/ACA interferometers in detail in §6.1.1.

Photometric measurements using direct-detection instruments
While many pioneering SZ studies were performed using interferometers, there has been a nearly universal shift towards the use of photometric imaging arrays for SZ observations, particularly in the case of wide-field SZ surveys where large mapping speeds are required. Currently, state of the art instruments for targeted observations contain a few thousand individual detectors covering 1-3 separate broadband (∆ν/ν ∼ 20%) spectral channels placed within the atmospheric windows centred at roughly 90, 150, and 260 GHz. Coupled to large aperture telescopes (see e.g. § 6.1.2, § 6.1.3, and § 6.2.1), current and near-future instruments obtain ∼ 10 angular resolution over a ∼ 5 FoV. As angular resolution θ at wavelength λ scales as the inverse of the primary mirror diameter D (see §A.1 for details), achieving much higher resolution than this  Fig. 18 Mock ACA and ALMA observations at 90 and 220 GHz of the simulated cluster of Figure 6. The maps were generated applying a natural weighting scheme to the raw visibilities and without correcting for the incomplete sampling of the uv-space through deconvolution. The resulting synthesised beams for each of the observations are also reported. The top row shows the 90 GHz input map, a mock ACA observation, and a mock ALMA observation from left to right. The bottom row is the same, but for 220 GHz. The simulated observations are noise-free in order to highlight solely the filtering effects of the interferometer. requires either an interferometer ( §5.1) or a much larger (> 100 m) single-dish largeaperture telescope.
To reasonable approximation, most modern photometric arrays are limited by photon statistics (i.e., the intrinsic detector and instrument noise is negligible), and their point source sensitivity scales as: where η is the optical efficiency, T load is the effective temperature of the total photon loading on the detector, ν 0 is the observing band-centre, A c is the collecting area of the single-dish telescope, and N det is the number of detectors. The first term in the numerator represents photon shot noise, while the second represents photon Bose noise, and the two terms tend to be approximately equal for ground-based receivers at relevant SZ frequencies. Comparing Equations (16) and (17), it is clear the noise improves with detector count as √ N det for photometric arrays, while the noise improves with baseline number as √ N bl for interferometers. At the observing frequencies where photometric arrays operate, spatial and temporal variations of the column depth of water vapour in the atmosphere produce a signal with a peak-to-peak amplitude of 1 K over long integrations, many orders of magnitude brighter than the SZ effect. To reasonable approximation, these variations occur within a single layer in the atmosphere with a fixed spatial pattern that moves across the telescope's FoV at an angular rate of ∼ 30 /sec (e.g., Lay & Halverson 2000;Sayers et al. 2010). The power spectrum of these variations rises sharply at large angular scales, and tends to dominate the detector data time-streams at frequencies below ∼ 1 Hz for telescope scan rates slower than or comparable to the atmospheric angular drift speed. Therefore, faster telescope scanning, at least up to the point when it is much faster than the atmospheric drift, mitigates the atmospheric contamination by modulating the astronomical signal to higher frequencies in the data time-streams above the dominant atmospheric signal at 1 Hz. To achieve this fast scanning while keeping the cluster within the instrument's FoV, the telescope is often moved in a lissajous 'daisy' pattern (e.g. Korngut et al. 2011) centred on the cluster, although linear raster scans are still used in some cases (e.g. Adam et al. 2018b).
For a single-dish large-aperture telescope, all of the detectors in the focal plane look through a nearly identical column of atmosphere, and so the atmospheric signal at any given instant in time is approximately common in all of the detectors. As a result, this unwanted signal is generally subtracted via a common-mode template(s), usually in combination with a high-pass filter applied to the data time-streams. The largest astronomical scales recoverable after this common-mode removal are thus comparable to the instrument's FoV. An example of this is given in the discussion of MUSTANG-2 ( §6.1.2).
Looking forward, advancements in detector technologies, particularly kinetic inductance detectors (KIDs), hold the promise of far larger detector counts when coupled to high-throughput telescopes. Such advancements would allow for broader spectral coverage (e.g., 6 bands between 90-400 GHz) and significantly larger FoV ( 1 • ). While there are prospects for wide FoV telescopes with improved mapping speeds and ability to recovered larger angular scales (see §6.2.8), there is no immediate prospect for the larger telescope apertures (> 100 meters) required for substantial gains in angular resolution. For that, an interferometer is required ( §5.1).

Astrophysical contamination
As the sensitivity of SZ effect measurements has improved, it has become increasingly important to mitigate contamination from unwanted astrophysical sources coincident on the sky with the cluster.
Radio sources: Early SZ effect observations, mainly at ≤ 30 GHz, often included nonnegligible synchrotron emission from field and cluster-member galaxies, particularly the BCG (e.g., Cooray et al. 1998;Coble et al. 2007). While synchrotron emission is much less problematic at the higher frequencies typical of modern SZ observations (≥ 90 GHz), it can still pose a challenge to higher resolution measurements (e.g., Mason et al. 2010;Romero et al. 2017). However, since the luminosity function of synchrotron sources is relatively flat (e.g., de Zotti et al. 2010), the additional contamination in future, deeper SZ effect observations is expected to be modest. Furthermore, the typical synchrotron spectrum, which decreases as a power-law with frequency, is readily separable from both the tSZ and kSZ effect spectra and their relativistic corrections. Ultimately, the level of contamination must be assessed as new deeper observations are acquired.
Primary CMB Anisotropies: Another astrophysical contaminant comes from the primary anisotropies in the CMB, which become significant ( 1µK) at angular scales larger than 5 (Planck Collaboration et al. 2016a). This is well inside the virial diameter of most clusters, and these CMB anisotropies can therefore impact measurements of large-scale cluster features in SZ images (e.g., Plagge et al. 2010). Fortunately, the CMB spectrum can easily be distinguished from the tSZ effect spectrum, and multi-band observations can effectively separate the two signals (e.g., Planck Collaboration et al. 2016b). The kSZ effect, however, has an identical spectrum to the CMB, and so the only way to separate the two signals is via their spatial templates. At small angular scales where the CMB is relatively smooth, for example when measuring the kinematic SZ signal towards sub-structures in the ICM, this spatial separation can work quite well (e.g., Sayers et al. 2013a;Romero et al. 2017;Adam et al. 2017b;Ruppin et al. 2018). However, separating the full-cluster-scale kSZ effect signal from the CMB is in general significantly more difficult (e.g., Planck Collaboration et al. 2014b, 2018. Dusty Galaxies: Looking towards future SZ effect measurements, thermal dust emission from background galaxies may be the most problematic astrophysical contaminant. Unlike the luminosity function for synchrotron emission, the luminosity function for the thermal dust emission is quite steep (e.g., Glenn et al. 2010), and so the relative amount of contamination grows quickly with deeper SZ effect observations. The average dust temperature in these galaxies is ∼ 40 K (e.g., Chapman et al. 2010), and their redshift distribution peaks near z ∼ 1-3 (e.g., Smith et al. 2017), resulting in a typical observed spectrum with an effective temperature of ∼ 10-20 K. While such a spectrum is highly degenerate with the shape of both the tSZ and kSZ signal at lower frequencies (ν 100 GHz), it can be distinguished with coverage extending to higher frequencies (e.g., Sayers et al. 2013a;Erler et al. 2018). Furthermore, because of the typical sub-arcsecond angular size of these galaxies, they will be unresolved in nearly all SZ effect observations, potentially facilitating spatial separation from any resolved SZ effect features.
While in principal the spatial and spectral properties of these dusty galaxies make them highly distinguishable from the SZ signal, there are many practical challenges to their subtraction from SZ images. Their steep luminosity function results in an effective noise floor due to source confusion that is only a modest function of angular resolution for FWHMs 5 (e.g., Béthermin et al. 2011). Therefore, even relatively high resolution SZ images are likely to contain a non-negligible signal from multiple dusty galaxies per resolution element. Even in the simplest scenario of a single dust temperature per galaxy, this will result in a complex multi-temperature dust spectrum within each resolution element (e.g., Roseboom et al. 2010), and it is unclear how well such a spectrum can be modelled even with broad spectral coverage. Lensing of this galaxy population by the cluster can add complicated spatial variations that present an additional challenge (e.g., Zemcov et al. 2013;Sayers et al. 2018). Further study on this topic is required, but SZ sensitivities significantly better than the confusion limit from dusty galaxies may be extremely difficult to achieve. One potential solution to this problem is ancillary observations from a very high resolution facil-ity such as ALMA, although such observations would be prohibitively expensive for large cluster samples.
Galactic Diffuse Emission: Mainly for clusters located at low galactic latitudes, the diffuse signal from within our galaxy due to synchrotron, free-free, spinning dust, and thermal dust emission can present another possible contaminant (Planck Collaboration et al. 2016d,g). However, even when these signals are non-negligible, they tend to be relatively smooth on the angular scales relevant for clusters. Furthermore, their spectral shapes tend to be highly distinguishable from both the tSZ and kSZ effect signals. As a result, contamination from galactic emission is unlikely to be a major impediment to upcoming SZ effect studies.
SZ confusion: Since the SZ effect signals contain the contributions from all free electrons along the line of sight, with effectively no redshift dependence, the SZ effects produced by the intervening structures (i.e., ionised gas in and around galaxies, groups, clusters, and cosmic web of filaments) in the background and foreground of the cluster will contaminate the SZ signal of the cluster. This SZ confusion is small for massive clusters, but it could become significant for low-mass clusters and groups (White et al. 2002; see also Voit et al. 2001 for the confusion in the X-ray band). Furthermore, the SZ contribution from filaments in the cosmic web is expected at the level of 8%-10% (Hallman et al. 2007). The SZ confusion is also a significant issue for stacking and cross-correlation analyses as well.

Current and future instrumentation
The section provides an overview of current-generation instruments (e.g., MUSTANG-2, NIKA2, ALMA+ACA), from which we provided a science highlights already, and presents a few key aspects of planned and proposed future instruments. Here we caution the reader when comparing measured mapping speeds of current instruments with forecasts for background-limited performance with future instruments. The background-limited mapping speeds have historically been a factor of several times too optimistic in some cases.

Current instrumentation
In this section we summarise the current generation of instrumentation that have measured the tSZ decrement at subarcminute resolution. Table 1 lists the salient aspects of each.

ALMA/ACA
The Atacama Large Millimeter/Submillimeter Array (ALMA) is an observatory consisting of 3 arrays: fifty 12-meter antennas (usually called "ALMA"), twelve 7-meter antennas known as the Morita / Atacama Compact Array (ACA), and four 12-meter antennas known as the 'Total Power Array' (TPA), which share the same design as the 12-meter elements of ALMA.
The best imaging fidelity comes from the 12-meter ALMA main array due to the large collecting area, which achieves low noise, and the superb sampling of the sky in uv-space made possible by having up to 1225 baselines. However, the 12-meter array has the smallest primary beam (and hence FoV, see §5.1), and only recovers scales ∼ 0.5 at its lowest current frequency (see right hand panels of Figure 18).
The ACA observes in a more compact configuration and, due to its smaller 7meter elements, has a FoV 1.7 times larger than that of the main array at the same given frequency. As a trade off, the ACA has less collecting area, lower resolution (∼ 13 in Band 3), and sparser sampling of the sky in uv-space (up to 66 baselines), making imaging more challenging (see middle panels of Figure 18).
And finally, the TPA, which nominally holds the promise of single dish total power measurements, is unfortunately not suitable for SZ observations as it fails to recover continuum signals like the tSZ and kSZ, and -unlike most single dish telescopes in operation -only samples the sky with a single beam per telescope. This is due in large part to the small FoV of the 12-meter elements; the antenna design was not optimised for focal plane arrays. Additionally, the overlap in Fourier modes sampled by a 12-meter is generally insufficient to provide good joint imaging, and several studies have shown a single dish should be at least 3 times the diameter of the interferometric elements with which the data will be combined (see e.g. Frayer 2017). A new, large aperture (> 36 meter), wide FoV telescope could finally provide this capability, improving dramatically the spatial dynamic range of ALMA (e.g. AtLAST, §6.2.8).
Upgrades to ALMA in the next few years will deliver bands covering 35-51 GHz ('ALMA Band 1') and 67-90 GHz ('ALMA Band 2'), with the possibility the latter will be extended to cover 67-116 GHz. Both bands are firmly at frequencies best suited for probing the tSZ (see Figure 2). ALMA Band 1 will improve ALMA's ability to probe larger angular scales by a factor of 2.5. While smaller than the instantaneous FoV bolometer arrays like MUSTANG-2 ( §6.1.2), NIKA2 ( §6.1.3), and TolTEC ( §6.2.1), the larger FoV and the ability of ALMA to jointly image a target in multiple array configurations at the same frequency will further improve imaging fidelity over a broad spatial dynamic range. Additionally, ALMA Band 2 will deliver wider bandwidths (≥ 16 GHz) that anticipate upgrades to the ALMA correlator, and covers frequencies where the tSZ is a factor of ∼ 4× stronger than it is in Band 1. Such capabilities will dramatically improve sensitivity to pressure substructures, shocks, and fluctuations in the tSZ.

MUSTANG-2
The 2nd generation Multiplexed SQUID/TEC Array at Ninety Gigahertz (MUSTANG-2) is a continuum bolometer camera on the 100-meter Green Bank Telescope (GBT) in Green Bank, West Virginia (USA). It consists of 215 feedhorn coupled polarisation insensitive detectors in the atmospheric window centred about 90 GHz (3.3 mm), useful primarily for probing the tSZ decrement. The 100-meter aperture of the GBT yields a 9 resolution, positioning it as one of the highest resolution bolometer cameras to probe the thermal SZ effect. Several examples of MUSTANG-1 results are shown in, e.g, § 3 and § 4. The upgrade from MUSTANG-1 to MUSTANG-2, while bringing newer, more sensitive detectors, was also motivated in large part by the need to recover larger scales.
MUSTANG-2 is now being demonstrated to recover scales comparable to its 4.2 FoV, while MUSTANG-1 was limited to scales similar to its 42 FoV (see discussion of photometric arrays in §5.2).

NIKA2
NIKA2 (Adam et al. 2018b;Catalano et al. 2018) is a millimetre camera, consisting of dichroic arrays of KIDs operating at 150 mK. It was installed at the focus of the IRAM 30-meter telescope in September 2015. NIKA2 observes the sky at 150 and 260 GHz with a wide FoV, 6.5 , a high angular resolution (17.7 and 11.2 at 150 and 260 GHz, respectively), and state-of-the-art sensitivity (8 and 30 mJy.s 1/2 , respectively). It features 616 polarisation-insensitive detectors at 150 GHz, and two arrays of 1140 detectors each at 260 GHz to provide polarisation measurements in this band. The camera's performance at the IRAM 30-meter telescope is described in Adam et al. (2018b); Catalano et al. (2014). The NIKA2 camera is well-suited for high-resolution SZ observations for several reasons. First, it is a dual-band camera operating at frequencies for which the thermal SZ signal is negative and slightly positive respectively. The 260-GHz map may be used for the detection for point sources or as a template of the atmospheric noise. Then, a high angular resolution and a large FoV enable NIKA2 to well match intermediate and high redshift clusters.

A comparison of current instrumentation
The typical scale used in cosmological applications is r 500 , a convention adopted due to the limited ability of state of the art X-ray instruments to probe the full extent of a cluster in reasonable exposure times. One of the new frontiers for cluster astrophysics is to extend the measurements out to and beyond the virial radius ( 3 × r 500 ; see e.g. Walker et al. (2019)) and avoid loss of signal due to spatial filtering, which implies it will be important to build more sensitive, multi-band instruments and telescopes that can access tens of arcminute scales even for z > 0.5. The cluster angular scale corresponding to r 500 is plotted as a function of redshift to highlight the limitations of current instruments in Figure 19. Over-plotted are the radial scales accessible for a few of the instruments discussed in §6.

TolTEC
TolTEC is a new multi-chroic photometric camera being built for the 50-meter Large Millimeter-wave Telescope (LMT) at Sierra Negra, Mexico for operation in 2019. The optical design simultaneously couples the 4 FoV onto 6,300 polarisation-sensitive detectors, divided into three focal planes. Filter passbands will be centred at 150, 220, and 280 GHz, spanning the null in the tSZ effect and peak in the kSZ. The large Table 1 Measured mapping/imaging speeds for current generation subarcminute resolution instruments. The photometric mapping and interferometric imaging speeds are reported for the average central depth for the map size reported. For example, in order for the ALMA 12-meter array in Band 3 to recover a depth of 84 µK over a region as large as a typical NIKA2 or MUSTANG-2 map, one would require ∼55-60 hours of integration time. Note that the map depth generally scales as t −1/2 only when not near the confusion limit. Also, note that the ALMA sensitivity calculator yields sensitivities in Rayleigh-Jeans brightness temperature, rather than ∆T CMB (see §A.4).  Fig. 19 The radial angular extent corresponding to r 500 , θ 500 , of a galaxy cluster as a function of redshift, plotted for a range of masses. The radial scales accessible to a few instruments, namely the ALMA 12meter array in Bands 1 and 3 (see §6.1.1), MUSTANG2 on the 100-meter GBT (see §6.1.2), and NIKA2 on the IRAM 30-meter ( §6.1.3), are plotted for comparison. For nearly all massive clusters at z 1, instruments with ≈ 4.5 FoV are able to recover scales comparable to θ 500 , but provide only a filtered view of scales larger than this radius. Note that TolTEC on the 50-meter LMT ( §6.2.1) has a ∼ 4 FoV, comparable to that of MUSTANG-2. Note also that we have assumed for ALMA that we recover scales as large as the primary beam, while typical scales recovered in the compact configuration are 0.6× that. And finally, note that the scales accessible to the ACA 7-meter array are a factor of 1.7× larger than those accessible to the ALMA 12-meter array at the same frequency. telescope on which this camera will be mounted will result in significant angular resolution gains from previous generation instruments, leading to a 150 GHz FWHM beam size of 9.5 . This is sufficient to measure the full profile of the hot gas on the scale of individual large galaxies. The camera is also expected to acheive mapping speeds of 10 deg 2 /mJy 2 /hour or better at 150 GHz, depending on the coherence scale of atmospheric fluctuations at the site (Bryan et al. 2018).

KISS
The KID Interferometer Spectrometer Surveyor (KISS) is a KID-based MPI spectroscopy camera to be installed at the QUIJOTE 1 telescope at the Teide observatory (at 2390 meters a.s.l. in Tenerife) in November 2018. KISS, which shares a heritage with CONCERTO ( §6.2.3), will operate 600 spatial pixels in the frequency range spanning 70-260 GHz with spectroscopic resolution of ∼3 GHz. With a resolution of ∼ 3 , KISS will target observations of low redshift (z < 0.2) clusters to map their thermodynamic properties via the tSZ effect, and their velocity and mass via the kSZ effect.

CONCERTO
CONCERTO 13 , the CarbON CII line in post-rEionisation and Reionisation epoch (Lagache 2018) is a 5000-element KID-based camera, styled off of NIKA2, coupled to an Martin-Puplet interferometer to provide 1.5 GHz resolution spectroscopy in the atmospheric windows between 125-360 GHz. The instrument will be installed on the APEX telescope, at 5100 meters a.s.l., approximately 100 meters higher than the plateau populated by ALMA. While designed primarily to target [CII] 158 micron emission from high redshift (4 < z < 8) galaxies, it will be capable of targeted tSZ and kSZ observations.

TIME
The TIME instrument (Crites et al. 2014) is designed to measure highly-redshifted 158 micron [CII] (1.9 THz) line emission. However, due to its wide frequency coverage and 150 GHz continuum photometry channel, the instrument will be capable of measuring the tSZ and kSZ in targeted clusters of galaxies. The instrument is a 183-336 GHz spectrometer with a photometer channel centred at 150 GHz. The spectrometer will have a resolution of approximately 100 with 60 spectral elements spanning the 183-336 GHz wavelength range. The instrument will have 16 spatial elements each with a 0.4 beam arranged in a single row. One challenge that normally limits kSZ measurements is effective removal of atmospheric fluctuations. The TIME spectrometer architecture is ideal for removing this source of contamination. Development of the instrument is well underway and is slated to be installed at the Arizona Radio Telescope (ARO) 12-meter prototype ALMA antenna in 2019.

Prime-Cam / CCAT-prime
CCAT-prime will be a 6-m aperture submillimetre telescope located at 5600 meter altitude in Chile and beginning its expected operations in 2021 (Stacey et al. 2018). Prime-Cam is proposed as the first-light bolometer instrument for this telescope, that will cover the frequency range between 190 and 450 GHz with multi-chroic TES detectors, plus host a broadband KIDs array at 860 GHz (Vavagiakis et al. 2018).
The high-altitude location means excellent atmospheric transmission at submillimetre wavelengths, which will place CCAT-prime in a unique position to explore the increment of the tSZ spectrum. The expected sensitivity in the 220 − 400 GHz bands will be an order-of-magnitude or more better than the Planck all-sky data (Table 3, Erler et al. 2018), thanks to the roughly 10,000 detectors in each band hosted inside a wide focal plane (roughly 1.5 • FoV). Combined with mm-wavelength (tSZ decrement) data either from separate optical tubes in the Prime-Cam or from other CMB surveys like advACT and Simons Observatory, CCAT-prime will provide broad coverage of the tSZ spectrum similar to Planck within a ∼ 10, 000 deg 2 sky area, which is expected to facilitate the separation of the tSZ, kSZ, and rSZ contributions in hundreds of high-mass clusters (Mittal et al. 2018, Erler et al. 2018). In the 860 GHz band, the angular resolution will be 15 times better than Planck, providing significant improvements in modelling the CIB and subtracting the contribution from clustercentric dust emission (e.g., Melin et al. 2018a).

Upgrades to current CMB survey instruments
Advanced ACTpol, which is the current generation instrument on ACT, and SPT-3G, the current generation camera on SPT, will also see some upgrades that will allow them to detect on the order of thousands of SZ selected clusters. Both are in the field operating. SPT-3G for example is predicted to find ∼5000 clusters at a signal-to-noise ≥ 4.5 (Benson et al. 2014). The ACT 6-meter will join the Simons Observatory, and both SPT and ACT will likely become part of CMB-S4 ( §6.2.7).

Future CMB survey instruments: Simons Observatory and CMB-S4
The next decade of CMB survey instruments, while not showing any dramatic improvements in spatial resolution, will continue to progress to large detector counts and additional bands from ∼30-300 GHz. First and foremost, the Simons Observatory (S.O.; The Simons Observatory Collaboration et al. 2018) will combine several existing CMB experiments in the Atacama desert, and add a new 6-meter telescope with a similar optical design to CCAT-prime ( §6.2.5). The S.O. large aperture telescope receiver (LATR) will feature ∼60,000 detectors spanning 6 broad receiver bands, centred at 27, 39, 93, 145, 225 and 280 GHz. Combined with CCAT-prime, this will deliver strong constraints on the kSZ effect and on relativistic corrections to the tSZ effect. The broad multi-band capabilities should significantly advance the ability to separate SZ components from each other and contamination by radio synchrotron and dusty sources. It has an anticipated first light in 2021.
Looking further ahead, instrument developments for the S.O. anticipate many of the CMB survey improvements expected for CMB-S4 (Abazajian et al. 2016). The details of the CMB-S4 project are still being determined, but CMB-S4 will likely add up to three 6-meter antennas of similar design as the S.O. and CCAT-prime 6-m, and several more lower resolution 1-meter class antennas. CMB-S4 has an anticipated start date circa 2028. Its dramatic leap forward relies on detector counts 1-2 orders of magnitude larger than S.O. and CCAT-prime ( §6.2.5), finally reaching megapixel counts.

AtLAST
The Atacama Large Aperture Submillimeter/millimeter Telescope (AtLAST) 14 is a broad international community-driven effort to build a 50-meter class single dish observatory with over 1 square degree instantaneous FoV in the Atacama desert in northern Chile. It will likely be located at ≈ 5100 meters a.s.l., close to the ALMA site, which has demonstrated performance over all full ALMA bands (35-950 GHz) and would facilitate accessibility and operational logistics (De Breuck 2018).
Conservative estimates assuming a 1 megapixel photometric array for AtLAST indicate its mapping speed will be well over 10 5 × faster than the full ALMA observatory (c.f. Equations 16 and 17). Currently no submm facilities larger than 12 meters are operational in the southern hemisphere, and AtLAST will serve as a powerful complement to the many lower-resolution primary CMB and SZ survey telescopes located in or planned for the southern hemisphere (see §6.2.7). At the same wavelengths as the LMT 50-meter ( §6.2.1), AtLAST will have the same spatial resolution. However the large (> 1 degree 2 ) FoV, planned higher surface accuracy, and better atmospheric transmission from 5100 meters a.s.l. will allow for simultaneous measurements in the submillimetre (∼ 350 − 950 GHz), to better probe contamination from dusty sources and the increment in the SZ effects. AtLAST will thus build upon the results with CCAT-prime ( §6.2.5) and deliver twice the resolution of the original CCAT 25-meter.

CSST
The Chajnantor Sub/millimeter Survey Telescope (CSST) is a Caltech-led effort to build an inexpensive 30-meter class single-dish survey telescope operating between 90-420 GHz. The concept is based upon a novel telescope design that would allow for a significant reduction in cost in exchange for a somewhat limited range of motion (∼1 radian, Padin 2014). CSST will seek a site on the Chajnantor plateau. The nominal first-light instrument would be a photometric camera with 6 observing bands filling a 1 deg FoV. As a survey telescope, CSST would have access to ∼20,000 deg 2 (half the sky) and would be able to map ∼1,000 deg 2 /year to the confusion limit due to dusty sources at 350 GHz. The actual survey design would likely consist of a "wedding cake" to varying depths. While pointed cluster observations would also be possible, the survey data themselves would provide extremely deep SZ maps of a very large sample of clusters, thus enabling a range of novel SZ science by significantly expanding the spectral coverage available from other large-aperture telescopes like GBT, IRAM, and LMT while providing significantly better angular resolution than any of the planned CMB survey telesopes.

Future space missions
There has also been significant discussion of and planning for SZ science from spaceand balloon-borne platforms. Planck has clearly demonstrated that wide frequency coverage, control of systematics, and coverage of a large sky fraction are some of the key ingredients for obtaining and doing science with large samples of SZ clusters (e.g., Planck Collaboration et al. 2011. Subarcminute angular resolution at mm wavelengths is challenging for space (and near-space), as the size of the primary mirror strongly drives the cost. However, future CMB mission concepts such as CORE (Delabrouille et al. 2018;Remazeilles et al. 2018) and PICO could take a dramatic leap forward in sensitivity and spectral coverage versus that of Planck, enabling an increase in the number of detected SZ clusters to a few × 10 5 (Melin et al. 2018b). This will open the way for detailed SZ studies with unprecedented precision, allow a vast number of cross-correlation studies, and provide additional constraints on cosmological parameters and neutrino mass and number (e.g., see Hill & Pajer 2013;Bolliet et al. 2018;Makiya et al. 2018;Salvati et al. 2018;The Simons Observatory Collaboration et al. 2018).
Higher angular resolution and wide spectral coverage across the mm, submm, and FIR from space could be achieved with the Millimetron Space Observatory (MSO) (Smirnov et al. 2012(Smirnov et al. , 2018. This concept uses a 10-meter folded mirror approach to overcome some of the payload limitations of the launch vehicle. The mirror will be cooled to ∼ 4.5 K to reduce the noise contribution due to thermal emission (which is already low, due to the low intrinsic emissivity of the mirror). MSO will target individual clusters, resolving them at subarcminute resolution with the goal of extracting the rSZ signals and reconstructing cluster profiles.
Finally, even low-resolution all-sky CMB measurements using spectrometer concepts such as PIXIE (Kogut et al. 2011(Kogut et al. , 2016, PRISTINE, and the spectrometer on board of CMB-Bharat could provide a new census of the hot baryons in the Universe by constraining the average all-sky y-parameter and y-weighted temperature . These concepts use a Fourier-Transform-Spectrometer (FTS) approach to sample the absolute CMB sky intensity in hundreds of frequency channels covering a few×(10 GHz − 1 THz). A similar FTS concept was used by COBE/FIRAS to measure the CMB energy spectrum in the 1990's, and by Olimpo (Masi et al. 2003;Schillaci et al. 2014), a balloon-borne SZ experiment targeting clusters that had its first flight in mid-2018.

Conclusions
SZ surveys and a few pioneering instruments have made tremendous progress in the past two decades, and detections and imaging of the SZ effect from massive systems has now become routine. With upcoming advances brought by new instruments, observing across the centimetre/millimetre/submillimetre spectrum (∼ 30 GHz−1 THz) and delivering higher spatial resolution over larger FoV, this progress will continue. It will soon be possible to characterise the intra-group, circumgalactic, and warm/hot intergalactic medium through SZ observations. In order to be able to probe small scales and achieve better spatial fidelity than currently possible, any future SZ facility must not only have high angular resolution, but also a large FoV and broad frequency coverage over most of the mm/submm bands. Fortunately, these are generic requirements to future mm/submm facilities, not endemic to galaxy clusters and SZ measurements.
New facilities should aim for imaging of scales ranging from sub-kpc to several Mpc, and perhaps larger for the warm-hot component of the cosmic web.
In the near term, our SZ window on the warm and hot Universe could potentially benefit by performing a few key, deep ( 100σ) SZ observations, analogous to deep X-ray and HST observations, to address specific topics in astrophysics. For comparison, the most transformative results with Chandra have often come out of very deep (≥ 400 ksec) and 'visionary' (∼ 2 Msec) 15 observations of shocks, cold fronts, AGN outbursts, and nearby clusters (see e.g. McNamara & Nulsen 2007;Markevitch & Vikhlinin 2007;Blanton et al. 2011;Zhuravleva et al. 2014b;Walker et al. 2018). Comparably deep SZ observations should exploit the key strength common to the SZ effects -redshift-independent surface brightness -by e.g. probing the evolution of clusters and proto-clusters in the first half of cosmic history. SZ observations have only begun the scratch the surface, and one can imagine using the tSZ to image the ICM farther into the cluster outskirts, to the accretion shocks and perhaps into the surrounding filaments.
Two decades ago, measurements of mK level temperature fluctuations were state of the art; now, large surveys and deep, targeted observations are achieving RMS values in the few µK range. Our ability to probe the tSZ effect will further advance, and what is now possible through stacking and radial binning will eventually lead to direct imaging of the intra-group, circum-galactic, and warm-hot intergalactic media. Apart from relativistic effects, which provide a handle on ICM temperature, the tSZ has no formal band dependence and can probe gas that is nearly invisible in the soft and medium X-ray bands.
Our ability to probe clusters and other structures through the kSZ will also improve through current and near future instrumentation advances, making kSZ detections and studies of gas motions and cluster rotation more routine. Nonetheless, given the faintness of the kSZ signal in most clusters, such measurements will remain challenging. Astrophysical contaminants such as Galactic cirrus, dust, and radio and submm galaxies ( §5.3) remain difficult to separate and remove fully, due to limited spectral leverage and poor signal to noise of the individual components. Additionally, systems in which the kSZ signal is strong are often merging objects in which temperatures locally boosted by shocks may require non-negligible relativistic corrections. It is clear that to significantly advance the field, instruments going far beyond the state of the art -2-3 colour instruments plus sparse ancillary radio/submm data -will be necessary to constrain all of the SZ contributions as well as the contaminants. 16 Disentangling the velocity of the gas from its optical depth will likely also remain a challenge for kSZ analyses, as it is difficult to model the electron distribution in complex systems. A promising avenue will open when rSZ can be used to make spatially-resolved measurements of ICM temperature. At this point, SZ characterisation of the motions and bulk properties of the ICM and warm/hot baryonic compo- 15 Recall that 100 hours is 360 ksec, and most SZ observations are only a few to tens of hours. 16 We note that much of the kSZ and rSZ work reported here relied on Herschel observations for submm source subtraction, and the number of clusters observed by Herschel during its lifetime is quite limited. On the submm observational side, CCAT-prime ( § 6.2.5) may be the clearest near-term successor to Herschel, with nearly twice the resolution and a longer expected project lifetime, while e.g. ALMA ( § 6.1.1) or AtLAST ( § 6.2.8) could constrain the flux densities of compact sources directly in each band of interest. nents in the Universe may be less directly dependent upon X-ray data. However, we note that direct measurements of the mass-weighted temperature T mw may remain elusive, as rSZ yields a Compton y-weighted temperature unless higher order moments of the temperature-field become accessible (see discussion in § 3.3).
Near future instrumentation will also soon allow us to measure the non-thermal SZ effect ( § 2.4), which has the unique capability of measuring the composition of AGN driven bubbles/cavities ( § 4.6) and, potentially, could even allow one to constrain the nature of dark matter. Both have been subjects of longstanding interest across many branches of astrophysics and cosmology, including large scale simulations and the fields of galaxy evolution and black hole growth. Looking forward to the more distant future, perhaps we can expect stacking and radial averaging to bring more exotic aspects of the SZ, such as the polarised and multiple scattering effects ( § 2.5 and § 2.6), within reach.
practice the under-illumination yields a main beam which, if fitted by single Gaussian would have θ FWHM ≈ 9.5 . Surface imperfections create an error beam, such that the total beam volume is better estimated as Ω bm = 120 square arcseconds.

A.2 Surface brightness
The surface brightness S ν in units of flux density per beam (in Jy/bm) is related to intensity ∆I ν (Equations 2, 7, & 9) by integrating over the beam solid angle: Generically, one can convert between intensity ∆I and a change in the CMB temperature ∆T CMB using the derivative of the blackbody function. The ratio is where the primary CMB intensity normalisation I 0 was defined in Eq. (3).

A.3 CMB survey noise
The RMS noise in maps made using arcminute-resolution CMB instruments, such as those from ACT and SPT, are often compared in units of µK CMB − arcmin, defined as the RMS of the CMB temperature fluctuations ∆T CMB within a map created with pixels that each subtend a solid angle of 1 arcmin 2 . To convert this figure to the RMS of CMB temperature fluctuations within a given instrument's beam, one would divide by the square root of the beam solid angle in square arcminutes (Eq. 19). For example, a map made from an instrument with Ω bm = 120 square arcseconds (0.033 arcmin 2 ) with an RMS noise of 10 µK CMB per beam would correspond to 1.8 µK CMB −arcmin. This conversion assumes the noise properties in the maps are Gaussian and uncorrelated on the scales being binned, which is a simplification that is not generically applicable.

A.4 Rayleigh-Jeans brightness temperature
Many instruments report sensitivities in Rayleigh-Jeans brightness temperature which at ν > 40 GHz can diverge significantly from the temperature decrement in the CMB, ∆T , in units of K CMB . This relation can also be expressed in surface brightness units (e.g. using Eq. 20), as Brightness temperature can be trivially converted to units more directly applicable to CMB and SZ measurements. Following Finkbeiner et al. (1999), this conversion, called the 'Planck correction factor', is where x is defined as it was for Equation 2. In Table 2, we provide computations of the ratio of ∆T CMB /T b for a few representative frequencies.