{'fraction_non_alphanumeric': 0.04661877378291214, 'fraction_numerical': 0.03402556183130556, 'mean_word_length': 4.329574390687522, 'pattern_counts': {'":': 0, '<': 47, '<?xml version=': 0, '>': 24, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Redshift Evolution of the Mean Temperature, Pressure, and Entropy Profiles in 80 SPT-Selected Galaxy Clusters

(Abridged) We present the results of an X-ray analysis of 80 galaxy clusters selected in the 2500 deg^2 South Pole Telescope survey and observed with the Chandra X-ray Observatory. We divide the full sample into subsamples of ~20 clusters based on redshift and central density, performing an X-ray fit to all clusters in a subsample simultaneously, assuming self-similarity of the temperature profile. This approach allows us to constrain the shape of the temperature profile over 0<r<1.5R500, which would be impossible on a per-cluster basis, since the observations of individual clusters have, on average, 2000 X-ray counts. The results presented here represent the first constraints on the evolution of the average temperature profile from z=0 to z=1.2. We find that high-z (0.6<z<1.2) clusters are slightly (~40%) cooler both in the inner (r<0.1R500) and outer (r>R500) regions than their low-z (0.3<z<0.6) counterparts. Combining the average temperature profile with measured gas density profiles from our earlier work, we infer the average pressure and entropy profiles for each subsample. Overall, our observed pressure profiles agree well with earlier lower-redshift measurements, suggesting minimal redshift evolution in the pressure profile outside of the core. We find no measurable redshift evolution in the entropy profile at r<0.7R500. We observe a slight flattening of the entropy profile at r>R500 in our high-z subsample. This flattening is consistent with a temperature bias due to the enhanced (~3x) rate at which group-mass (~2 keV) halos, which would go undetected at our survey depth, are accreting onto the cluster at z~1. This work demonstrates a powerful method for inferring spatially-resolved cluster properties in the case where individual cluster signal-to-noise is low, but the number of observed clusters is high.

1. INTRODUCTION Galaxy clusters, despite what the name implies, consist primarily of matter that is not in galaxies. A typical cluster is well-modeled by a central dark matter halo (∼85% by mass) and a diffuse, optically-thin plasma (∼15% by mass). The response of this hot ( 10 7 K) plasma, known as the intracluster medium (ICM), to the evolving gravitational potential is one of our best probes of the current state and evolution of galaxy clusters. X-ray imaging and spectroscopy of the ICM allow estimates of the cluster mass profile via the spectroscopic temperature and gas density (e.g., Forman & Jones 1982;Nevalainen et al. 2000;Sanderson et al. 2003;Arnaud et al. 2005;Kravtsov et al. 2006;Vikhlinin et al. 2006;Arnaud et al. 2007), the enrichment history of the cluster via the ICM metallicity (e.g., De Young 1978;Matteucci & Greggio 1986;de Plaa et al. 2007;Bregman et al. 2010;Bulbul et al. 2012), the cooling history via the cooling time or entropy (e.g., White et al. 1997;Peres et al. 1998;Cavagnolo et al. 2008;McDonald et al. 2013), the feedback history via the presence of X-ray bubbles (e.g., Rafferty et al. redshift. The construction of a "Universal" pressure profile, for example, can allow comparisons to simulated galaxy clusters, as well as provide a functional form for matched-filtering algorithms, such as those that are used to select galaxy clusters using the Sunyaev-Zel'dovich (SZ;Sunyaev & Zeldovich 1972) effect. Much effort has been made to quantify the average temperature (e.g., Loken et al. 2002;Vikhlinin et al. 2006;Pratt et al. 2007;Leccardi & Molendi 2008b;Baldi et al. 2012), entropy (e.g., Voit et al. 2005;Piffaretti et al. 2005;Cavagnolo et al. 2009;Pratt et al. 2010), and pressure (e.g. Arnaud et al. 2010;Sun et al. 2011;Bonamente et al. 2012;Planck Collaboration et al. 2013b) profiles for low-redshift galaxy groups and clusters based on both X-ray and SZ selection. In all cases, the average profiles have a substantial amount of scatter at r 0.2R 500 , due to the presence (or lack) of a cool, dense core (e.g., Vikhlinin et al. 2006;Cavagnolo et al. 2009;Arnaud et al. 2010), but collapse onto the selfsimilar expectation at larger radii. This suggests that non-gravitational processes (e.g., cooling, AGN feedback) are important in the central region of the cluster while gravity is the dominant force in the outer region.
While the aforementioned studies have made significant progress in quantifying the average temperature, entropy, and pressure profiles of galaxy groups and clusters, they have focused almost entirely on low-redshift (z 0.2) systems. This is in part due to the relative ease with which one can measure the temperature profile in nearby systems, but also due to the fact that, until recently, large, well-selected samples of galaxy clusters at high redshift did not exist. This has changed, with the recent success of large SZ surveys from the Atacama Cosmology Telescope (ACT; Marriage et al. 2011;Hasselfield et al. 2013), Planck (Planck Collaboration et al. 2011, 2013a, and the South Pole Telescope (SPT; Vanderlinde et al. 2010;Reichardt et al. 2013). These surveys have discovered hundreds of new galaxy clusters at z > 0.3, allowing the study of galaxy cluster evolution for the first time out to z > 1 using large, homogeneous data sets.
In this paper, we present a joint-fit spectroscopic analysis of 80 SPT-selected galaxy clusters in the SPT-XVP sample (McDonald et al. 2013, Benson et al. in prep). Utilizing uniform-depth X-ray observations of these clusters we can, for the first time, constrain the redshift evolution of the average ICM temperature, entropy, and pressure profiles. We present the details of this analysis in §2, including the resulting projected and deprojected temperature profiles in §2.2 and §2.3, respectively. In §3 we infer the average pressure ( §3.1) and entropy ( §3.2) profiles. In §4 we discuss the implications of the observed evolution, specifically in the inner ∼100 kpc and outskirts (r R 500 ) of the mean pressure and entropy profiles, and assess any potential biases in our analysis. Finally, we summarize these results in §5 before suggesting future applications of these data. Throughout this work, we assume H 0 =70 km s −1 Mpc −1 , Ω M = 0.27, and Ω Λ = 0.73.

DATA AND ANALYSIS
2.1. Sample The majority of the observations for this program were obtained as part of a Chandra X-ray Visionary Project to observe the 80 most massive SPT-selected galaxy clusters at 0.4 < z < 1.2 (PI: B. Benson). This survey is described in more detail in McDonald et al. (2013) and Benson et al. (in prep). We begin our sample definition by identifying 91 galaxy clusters that are detected in the SPT 2500 deg 2 survey and have been observed with Chandra. We first exclude four clusters which are detected with the SPT at S/N < 6 (SPT-CLJ0236-4938, SPT-CLJ0522-5818, SPT-CLJ2011-5725, SPT-CLJ2332-5053), which gives our sample a uniform SZ selection at S/N ∼ 6.5. We exclude two additional clusters (SPT-CLJ0330-5228, SPT-CLJ0551-5709) which suffer from severe projection effects due to extended foreground sources (i.e., nearby, low-mass groups). In the remaining 85 systems, we identify a tight distribution of X-ray photon counts at r > 0.2R 500 (see Figure  1), and exclude four clusters with exceptionally high signal-to-noise (SPT-CLJ0658-5556 (the Bullet Cluster), SPT-CLJ2248-4431, SPT-CLJ0232-4421, SPT-CLJ0102-4915), which could dominate the stacking analysis, and one cluster with very low counts (SPT-CLJ0037-5047), which will not contribute meaningful signal. What remains is a sample of 80 clusters that occupy a tight sequence in signal-to-noise at large radii. These 80 clusters define our sample of massive (M 500 3 × 10 14 E(z) −1 M ⊙ ) galaxy clusters with uniform-depth X-ray imaging, spanning the redshift range 0.3 z 1.2. We note that, in the outermost annulus, the combined Galactic and extragalactic background is roughly an order of magnitude brighter than the source emission. However by combining ∼20 spectra ( √ 20 improvement) and joint-fitting the background in the on-and off-source regions ( √ 2 improvement), we improve the signal-to-noise from ∼3 per spectrum to ∼20 -enough to constrain the spectroscopic temperature to within ∼10%.
This sample of 80 clusters was divided into six subsamples, based on individual cluster redshift and the presence (or lack) of a cool core, in order to probe the redshift evolution and cooling-dependence of the universal temperature, pressure, and entropy profiles. For simplicity, and so that subsamples are of equivalent signal-to-noise, we divide the low-redshift and high-redshift subsamples in half, with the 50% "cuspiest" clusters, where cuspiness (α) is defined as the slope of the gas density profile at 0.04R 500 , making up the cool core (CC) subsample and the 50% least cuspy clusters defining the non-cool core (NCC) subsample. This yields the six subsamples summarized in Table 1. The choice of z = 0.6 as a dividing line was motivated by the desire to have an equal number of clusters in both the high-z and low-z bins. The mean redshift for the two redshift bins are z low = 0.46 and z high = 0.82, which provides a broad baseline for comparison to previous studies at z ∼ 0.1 (see §1).
Global cluster properties (e.g., M 500 , kT 500 ) used in this work are derived in Benson et al. (in prep), following closely the procedures described in Andersson et al. (2011). Briefly, we estimate R 500 , the radius within which the average enclosed density is 500 times the critical density, by iteratively adjusting R 500 until the measured Y X (≡ M g × kT ) satisfies the Y X,500 −M 500 relation (Vikhlinin et al. 2009), which assumes a purely self-similar evolution, M 500 ∝ Y X,500 E(z) −2/5 . Once -Number of X-ray counts per annulus in the energy range 0.7-7.0 keV for 85 SPT-selected galaxy clusters observed with Chandra. Blue and red curves correspond to clusters classified as cool core and non-cool core, respectively (see §2.1). The shaded gray band represents the measured Galactic and extragalactic background in each annulus. Five clusters which are either too high or too low signal-to-noise to be used in the stacking analysis are shown in black. Overall, the remaining 80 clusters have a tight distribution of counts at r > 0.2R 500 . We stress that this has not been scaled for exposure time -we have nearly uniform depth (∼2000 counts per cluster) over this full sample. the radius converges, kT 500 is measured in a core-excised annulus from 0.15R 500 < r < R 500 . Further details on the derivation of global X-ray properties can be found in Andersson et al. (2011). We point out that we have also examined the effects of using an M gas -derived temperature to normalize the temperature profiles, following Vikhlinin et al. (2006), and confirm that the results we will present below are independent of our choice of normalization.
In Figure 2 we show the distributions of redshift (z), mass (M 500 ), and total X-ray counts for the full sample of 80 clusters, as well as the low-and high-redshift subsamples. Here M 500 is computed assuming the Y X -M relation from Vikhlinin et al. (2009), as described in McDonald et al. (2013) and Benson et al. (in prep). This figure demonstrates that the low-z and high-z subsamples are comprised of similar clusters, in terms of their total mass, and have been observed to similar depths, allowing a fair comparison.  (z), mass (M 500 ) and total X-ray counts for the 80 clusters used in this work. Here M 500 is computed assuming the Y X -M relation from Vikhlinin et al. (2009), as described in McDonald et al. (2013) and Benson et al. (in prep). We show separately the distributions for low-z (blue; z < 0.6) and high-z (red; z > 0.6) clusters, demonstrating the similarity in mass and data quality between these two subsamples.

Joint Spectral Fitting
For each cluster, we extract spectra in 12 annuli spanning 0 < r < 1.5R 500 . At r > 0.3R 500 , we use logarithmically-spaced bins, while interior to 0.3R 500 , we choose a binning scheme which achieves fine sampling while also maintaining suitable signal-to-noise per annulus. The bin edges are tabulated in Table A.1. This two-part binning scheme yields roughly equal signal-tonoise in all bins, without requiring overly narrow/wide bins at any radius. The total number of annuli was chosen, through trial and error, so that the cluster-to-cluster scatter in kT /kT 500 within a given annulus is similar to the fit uncertainty (i.e., there is negligible improvement in scatter from widening the bins). The number of counts per bin, for each cluster, is shown in Figure 1. This figure demonstrates that, even in our outermost bin (1.2R 500 < r < 1.5R 500 ), we have ∼100 X-ray counts per cluster spectrum.
For each cluster, we extract spectra in each of these annuli, using ciao v4.6 and caldb v4.6.1.1, along with accompanying background (both blank-sky and adjacent to source) and response files. Blank-sky background spectra are rescaled based on the observed 9.5-12 keV flux and subtracted from both on-and off-source spectra, which are then simultaneously modeled in the procedures outlined below. All spectra are binned in energy with wide bins at low ( 1 keV) and high ( 4 keV) energy where the signal-to-noise is low. The average binning was ∆(log 10 E) = 0.1. Given a subsample of N clusters (e.g., 40 low-z clusters), we randomly draw N clusters, allowing repeats (i.e., bootstrap analysis). This means that, in some cluster realizations, the contributions from cluster i may be double-counted, while cluster j will be excluded. The spectra for these N clusters are simultaneously fit over the energy range 0.5-10.0 keV with xspec (v12.8.0; Arnaud 1996), using a combination of a single-temperature plasma (apec; Smith et al. 2001), a soft X-ray background contribution Note. -1 : Constrained to within 15% of the Galactic value from Kalberla et al. (2005). Other spectra of the same field share the same value of N H . 2 : Constrained to within 5% of the optical redshift (see Song et al. 2012;Ruel et al. 2013). Other spectra of the same cluster share the same value of z. 3 : Limited range to 0 < Z < 2Z ⊙ . 4 : 10 −14 4π(D A (1+z)) 2 nen H dV (apec, kT = 0.18 keV), a hard X-ray background contribution (bremss, kT = 40 keV), and a Galactic absorption model (phabs) 1 . These additional soft (Galactic foreground) and hard (unresolved CXB) emission models account for any residual emission after the blank-sky backgrounds are subtracted. The various free parameters of the plasma model are constrained as described in Table 2.
This method, which has 2N + 1 free parameters per annulus (see Table 2), makes the assumption that, in the subsample of N clusters, there is a universal temperature profile of the form kT /kT 500 = f (r/R 500 ). We do not make any further assumptions about the form of this function. We simultaneously fit spectra for all N clusters, including multiple observations where available, along with off-source background spectra (at r 3R 500 ) unique to each cluster. Goodness-of-fit was determined using the χ 2 parameter, with weighting based on Churazov et al. (1996), which has been shown to yield In all panels, the colored regions represent the data, with dark and light regions corresponding to 1σ and 90% confidence, respectively, while the solid line corresponds to the median value. These projected temperature profiles are derived from joint spectral fits, as described in §2.2. The best fit projected profile is shown in solid black, with the 1σ uncertainty range in dotted black (see §2.3). For comparison, average projected temperature profiles from previous studies Baldi et al. 2012) are shown where applicable. In general, the agreement with previous works at similar redshifts is qualitatively good. Note. -All uncertainties are 1σ. The procedures for producing these profiles are described in detail in §2.1-2.2. unbiased parameter estimates for spectra containing as few as ∼50 total counts. We repeat this process 100 times for each annulus in order to assess the cluster-tocluster variation in kT /kT 500 . We note that, while the metallicity (Z i ) is left as a free parameter, we do not have sufficiently deep data to put meaningful constraints on the shape of Z(r/R 500 ). We leave this parameter free so that, in the few clusters with a strong Fe xxv emission line, the temperature is not skewed high in order to improve the fit. This method was found to yield the smallest residuals in temperature for simulated clusters with low counts and unknown metallicity 2 .
In Figure 3 and Table 3 we show the result of this jointfit analysis for each of the 6 subsamples described in §2.1. This figure demonstrates our ability to constrain the projected temperature to within ∼10% over 100 realizations, despite the fact that the average individual cluster contributes only ∼100-200 X-ray counts to each of these bins ( Figure 1). The uncertainty range shown in Figure 3 represents the cluster-to-cluster scatter in the kT /kT 500 profile, which dominates over the statistical uncertainty in Colored curves correspond to the best-fit deprojected model and the 1σ uncertainty in this model. In the right-most panels, we overplot the cool core and non-cool core profiles for comparison. In all panels we show the average deprojected temperature profile from Vikhlinin et al. (2006) for comparison. This figure highlights the lower-temperature cores and outskirts in high-z clusters, relative to their low-z counterparts. The best-fit parameters which describe these curves are provided in Table 4. the joint fit. Overall the general shapes are as expected -centrally-concentrated clusters have cool cores, while their counterparts do not. The agreement with previous average projected temperature profiles (e.g., Pratt et al. 2007;Baldi et al. 2012) is qualitatively good. In particular, when the results from Baldi et al. (2012) are scaled to the same normalization (r/R 500 , kT /kT 500 ), they agree at the 1σ level with our low-z subsample. In general, the ICM in the outskirts of high-z clusters appears to have a steeper temperature gradient than the ICM in their low-z counterparts. We note, however, that the temperature profiles shown in Figure 3 are projected -in order to determine the true radial distribution of ICM temperature, entropy, and pressure, we must first perform a three-dimensional deprojection.

Deprojecting Mean Temperature Profiles
To determine the three-dimensional average temperature profile, we project an analytic function onto two dimensions and fit this projected model to the data. Our three-dimensional temperature model is inspired by Eq. 6 from Vikhlinin et al. (2006): (1) where x = r/R 500 . This equation models the temperature profile in two parts: 1) the core region, which has a temperature decline parametrized by a minimum temperature (T min ), scale radius (r c ), and shape (a cool ) (Allen et al. 2001) and 2) a broken power law with a characteristic inner slope (a), transition steepness (b), and outer slope (c), and a transition radius (r t ). Since we only sample the cool core (r 0.1R 500 ) with 2-3 data points ( Figure 3), we have removed two degrees of freedom from the more general parametrization presented in Vikhlinin et al. (2006), fixing a cool = 2 and a = 0, as otherwise the fit would be completely unconstrained in the inner parts. To project this function onto two dimensions, we follow the procedures described in detail by Vikhlinin (2006), which require two additional ingredients, aside from the analytic temperature profile: the three-dimensional electron density and metallicity profiles. For the latter, we assume the average metallicity profile from Leccardi & Molendi (2008a), but confirm that radically different metallicity profiles (i.e., flat, inverted) result in 5% differences in the deprojected temperature, and then only at low temperatures ( 2 keV). For the electron density profile, we utilize the deprojected density profiles for each cluster from McDonald et al. (2013). Since each bootstrapped temperature profile is actually the weighted average of N clusters, we compute the appropriate three-dimensional gas density profile as follows: where C i (r) is the number of X-ray counts for cluster i at radius r. This produces a mean gas density profile, weighted in approximately the same way as the mean temperature profile. For each bootstrap realization, we use the same N clusters in the calculation of both the mean temperature and density profiles. From these profiles, we project the temperature profile along a given line of sight, fol- Note.
-The functional form for these fits, which are shown in Figure 4, is provided in Eq. 2. This parametrization is adopted from Vikhlinin et al. (2006), holding a cool = 2 and a = 0 fixed. Both rc and rt are expressed in units of R 500 .
lowing Vikhlinin (2006). This procedure accounts for different contributions from continuum and line emission along the line integral, providing an accurate estimate of the projected single-temperature model. To correctly factor in the temperature-sensitive detector response, we convert to absolute temperature units (from kT /kT 500 ) using the average kT 500 for the subsample of N clusters. For each radial bin, the projected temperature was computed by numerically integrating along the line of sight over −4R 500 < r < 4R 500 as well as along the bin in the radial direction. Once complete, this procedure yielded a projected temperature profile which was fit to the data using a least-squares χ 2 minimization routine. This process was repeated for each realization of the projected temperature profile (Figure 3), allowing an estimate of the uncertainty in the deprojected model.
The resulting deprojected temperature profiles, along with the 1σ uncertainty regions, are shown in Figure 4 and Tables 3 and A.1. These results suggest that highredshift clusters tend to be, on average, ∼30-40% cooler at large radii (r ∼ 1.5R 500 ) than their low-redshift counterparts. However, when non-cool cores are considered on their own, this trend is reversed (albeit at low significance). This may be a result of, on average, more significant "clumping" in the outskirts of high-z clusters (Nagai & Lau 2011). Infalling subhalos will tend to have lower temperature and higher density than the surrounding ICM, leading to a bias towards low temperatures in the emission measure-weighted spectra if they are unresolved. We will discuss this possibility further in §4, in the context of the average entropy profile. Similar to previous works, we find that the peak temperature for relaxed (cool core) systems is reached at ∼0.2-0.3R 500 . The cores of high-z cool core clusters appear to be cooler than their low-z counterparts at the ∼30% level, an effect still pronounced in the the full ensemble. Given that our typical 1-σ uncertainty in the central bin is ∼15%, this is only marginally significant.

MEAN PRESSURE AND ENTROPY PROFILES
Previous studies (see §1) have measured the "Universal" pressure and/or entropy profiles by simply taking the mean or median of a large number of individual profiles, which are individually calculated as follows: This approach is unfeasible for this sample, given that each individual cluster only has ∼2000 X-ray counts and, thus, the individual kT (r) profiles are unconstrained. However, if we define P 500 and K 500 as follows (from Nagai et al. 2007): where n g,500 = (µ e /µ)n e,500 = 500f b ρ crit /(µm p ), and we assume that µ = 0.59 is the mean molecular weight and (Gonzalez et al. 2013) is the universal baryon fraction, then we can express the scaled pressure and entropy profiles as follows: These expressions allow us to measure P/P 500 and K/K 500 as a function of radius, without actually measuring P (r) or K(r) for any given cluster, as is usually done. Instead, we simply combine the average ρ(r)/ρ crit profiles  with the average kT (r)/kT 500 profile from our joint-fit analysis (see §2). The average pressure and entropy profiles, computed using Eqs. 6 and 7, are provided in Table A.1. 3.1. Pressure In Figure 5 we show the average pressure profile, P (r)/P 500 , from this work. We compare our joint-fit results to previous studies by Vikhlinin et al. (2006) (2013b), in order to account for small differences in mass range between samples. We apply an additional scaling to profiles derived based on XMM-Newton data, which accounts for a 16% difference in temperature normalization between Chandra and XMM-Newton for massive clusters (Schellenberger et al. 2014). Both A10 and P13 define P 500 ∝ M 2/3 500 ∝ kT 0.37 500 , so a 16% normalization error in temperature corresponds to a 10% normalization error in P/P 500 . Likewise, since R 500 ∝ M 1/3 500 ∝ kT 0.19 500 , we apply a 3% renormalization in R 500 to account for differences between the two X-ray telescopes. Figure 5 demonstrates our ability to reach larger radii than Arnaud et al. (2010), who used data from XMM-Newton, due to the fact that the Chandra ACIS-I field of view represents a larger physical size at high redshift. While we are unable to probe as deep into the outskirts as Planck, the unmatched angular resolution of Chandra allows us to sample a factor of ∼10 finer than the native Planck resolution in the cluster core. Thus, this work bridges the gap between A10, who primarily samples the inner pressure profile, and P13, who primarily samples the cluster outskirts.
For our low-z subsample, there is good overall agreement between our results and previous work. Specifi- In all panels, the colored regions represent the data, with dark and light regions corresponding to 1σ and 90% confidence, respectively, while the solid horizontal line corresponds to the median value. The best fit GNFW profile (Eq. 8) is shown in solid black, while comparisons to V06, A10, and P13 are shown in dotted and dashed lines. The relative spatial coverage of A10 and P13 are shown at the top, for reference. For the low-z subsample, the agreement between all three samples is excellent. For high-z clusters, there appears to be slightly less pressure both in the core (r < 0.1R 500 ) and outskirts (r > R 500 ).
cally, at 0.2R 500 < r < 1.5R 500 , our results are fully consistent at the ∼ 1σ level with both A10 and P13. At r < 0.2R 500 , the cool cores from A10 appear to have higher pressure. This may be due to their different cool core classification, based on central density, which also varies as a function of total mass . Alternatively, the difference could be due to a difference in centering -we use the large-area centroid, while A10 use the X-ray peak -but these should both agree for relaxed, centrally-concentrated systems. The difference in the core pressure is less dramatic when the full samples are considered.
Non-cool core clusters (both high-and low-z) appear to have a deficit of pressure in their cores relative to A10 and P13, both at the 1-2σ level. Again, this may be a result of different centering choices: A10 and, presumably, P13 choose the X-ray peak as the cluster center, while we have measured the cluster center at larger radii in order to minimize the effects of core sloshing (see e.g., ZuHone et al. 2010). Our choice of centering allows for the X-ray peak to be offset from the chosen center, which would lead to lower central pressure. An example of the difference between the X-ray peak and the "large-radius centroid" is shown in Figure 6 for SPT-CLJ0014-4952. At large radii, the cluster is clearly circular on the sky. However, the X-ray peak, presumably representing a sloshing cool core, is located at a distance of ∼ 0.3R 500 from the large-radius centroid. It is most likely these off-center emission peaks that are leading to the small differences between our stacked pressure profile and those of A10 and P13, specifically at r 0.1R 500 and r R 500 . We would argue that, for non-cool core clusters, our choice of the cluster center is more representative of the true dark matter distribution.
We find that high-z cool core clusters have lower cen- The X-ray peak and large-annulus centroid are depicted with blue and red crosses, respectively. The two red dashed circles correspond to radii of 0.3R 500 and R 500 . Bottom: Electron density profile, derived using the two different centers. For "disturbed" clusters the X-ray peak can be significantly offset from the largescale centroid, leading to a deficit of core pressure (∼50%) and slight excess pressure at large radii (∼15%) if the latter is chosen as the cluster center.

1.043
Note. -Generalized NFW fit to our average pressure profile (see Figure 5). For comparison, fits from Arnaud et al. (2010) and Planck Collaboration et al. (2013b) are shown. For clarity, rows that are bolded are original to this work. Values with a † mark have been held fixed during the fitting procedure. Quoted uncertainties represent the 1σ confidence intervals for each parameter. We note that there are strong covariances between parameters which are not encapsulated in these uncertainty ranges. For example, if all other parameters were held fixed, the uncertainty on P 0 would be of order ∼5% for all fits.
tral pressure than those at intermediate-and low-z, by factors of ∼3 and ∼6, respectively, consistent with our earlier work which showed a rapidly-evolving central gas density between z ∼ 0 and z ∼ 1, with the central, low-entropy core being less massive at z 0.6 (McDonald et al. 2013). This drop in central pressure with increasing redshift is a result of both lower central temperature (see Figure 4) and lower central density ) in relaxed, high-z clusters. We note that, under the assumption of hydrostatic equilibrium, increasing the gas mass of the core ought to decrease the central temperature (assuming dark matter dominates the mass budget). The fact that we observe the opposite implies that some additional form of heating, either gravitational (e.g., increased dark matter density, adiabatic compression) or feedback-related, is raising the central temperature of low-z clusters. The observed change in central pressure would seem to suggest that high-z X-ray cavities, inflated by radio-mode AGN feedback, should be a factor of 3 1/3 larger in radius in order to maintain the same energetics (PdV) as low-z cavities. This factor of ∼40% in size evolution is currently smaller than the typical scatter in cavity size for high-z clusters (Hlavacek-Larrondo et al. 2012, Hlavacek-Larrondo et al. in prep).
To allow a more direct comparison to A10 and P13, we fit the stacked pressure profiles with a generalized NFW (GNFW) profile following Nagai et al. (2007): where x = r/R 500 and f (M ) = (M 500 /3 × 10 14 h −1 70 M ⊙ ) 0.12 . This generalized version of the NFW profile (Navarro et al. 1997) is a broken power law with a characteristic inner slope (γ), outer slope (β), curvature (α), and scale radius (c 500 ). The results of this fit, which is constrained over the radial range 0.01R 500 r 1.5R 500 , are provided in Table 5. At low-z, these fits are consistent with earlier work by A10 and P13. At high-z, these represent the first constraints on the shape of the average pressure profile, allowing a comparison to simulations over 8 Gyr of cosmic time (see §4).

Entropy
In Figure 7 we show the evolution of the stacked entropy profile. Here, we compare our data to the nonradiative simulations of Voit et al. (2005), rescaled to ∆ = 500 by Pratt et al. (2010). This represents the "base" entropy due to purely gravitational processes. The solid purple line corresponds to the median entropy profile for nearby REXCESS clusters (0.05 < z < 0.18; Pratt et al. 2010), rescaled to include the crosscalibration differences between XMM-Newton and Chandra (see §3.1; Schellenberger et al. 2014). We are unable to compare to Cavagnolo et al. (2009) due to the lack of a normalized (r/R 500 , K/K 500 ) median profile, or the appropriate quantities to do such a scaling ourselves.
The agreement between our low-z average entropy profile and that of Pratt et al. (2010) is excellent, suggesting that there has been little evolution in the entropy structure of clusters since z ∼ 0.6. At large radii (r > R 500 ), the low-z average entropy profiles are consistent with the gravity-only simulations (Voit et al. 2005), with the exception of the very last NCC data point, suggesting that the process responsible for injecting excess entropy at r < R 500 is unimportant at larger radii. On average, the NCC clusters have ∼ 2.5× higher entropy in the central bin, corresponding to a factor of ∼4 increase in central cooling time (Cavagnolo et al. 2009). The average CC central entropy of ∼0.085K 500 corresponds to K ∼ 85 keV cm 2 (assuming kT 500 = 6.5 and z = 0.45), or a cooling time of ∼ 3 × 10 9 yrs (Cavagnolo et al. 2009). We note that our CC/NCC In all panels, the colored regions represent the data, with dark and light regions corresponding to 1σ and 90% confidence, respectively, while the solid horizontal line corresponds to the median value. In all panels we show the baseline gravity-only entropy profile from simulations (Voit et al. 2005) for comparison. In the right-most panels, we show results from Pratt et al. (2010) for comparison. In the low-z subsample, the agreement with Pratt et al. (2010) is excellent, suggesting very little evolution from z ∼ 0.1 to z ∼ 0.5. At high-z, the entropy flattens out at large radii (r R 500 ) for the cool core subsample, as well as CC+NCC sample, which may be indicative of increased clumping in these more-distant systems. In general, the entropy profile appears to be unevolving, or evolving very slowly, at r 0.7R 500 , suggesting a remarkably stable balance between cooling and feedback. separation is such that the CC class will contain both "strong" and "weak" cool cores, leading to a higher average cooling time (see e.g., Hudson et al. 2010).
Consistent with McDonald et al. (2013), we find no significant evolution in the core entropy of CC clusters over the full redshift range studied here, despite the fact that both the central density ) and temperature ( Figure 4) have evolved significantly (≥ 2σ) over the same timescale. The universality of the average CC entropy profile at r 0.7R 500 , shown more clearly in Figure 8, may be indicative of a long-standing balance between cooling and feedback processes on large scales (see also Hlavacek-Larrondo et al. in prep).
At z > 0.6, we measure a distinct entropy decrement at r > 0.5R 500 in high-z clusters when compared to their low-z counterparts. This flattening of the entropy profile is apparent in the high-z CC subsample as well as the combined sample (at > 90% significance in all cases), indicating that it is not being driven by one or two extreme clusters. Figure 8 shows that this flattening becomes significant at r 0.7R 500 . Such a flattening could plausibly be due to clumping of the intracluster medium (e.g., Nagai & Lau 2011) -we will discuss this possibility and the inferred clumping required to explain the measurements in the next section.
In summary, we find no significant evolution in the average entropy profile from z ∼ 0.1  to z ∼ 1 within r 0.7R 500 , suggesting that the balance between cooling and feedback is exceptionally wellregulated over long periods of time (∼ 8 Gyr). Fig. 8.-This plot compares the low-and high-z entropy profiles for three different subsamples. Here, the y-axis is the ratio of the measured entropy for the high-z subsample to the low-z subsample at a given radius. Error bars represent 1σ uncertainty. All three profiles are consistent with negligible evolution at r < 0.7R 500 and significant evolution at r > 0.7R 500 .

DISCUSSION
Below, we discuss the results of §2 and §3, comparing these to previous observations and simulations in order to aid in their interpretation. In addition, we investigate  (Schellenberger et al. 2014) (see discussion in §3.1). We show, for comparison, average profiles for simulated clusters with M 500 > 3 × 10 14 M ⊙ and at similar redshifts, taken from Battaglia et al. (2012), and Dolag et al. (in prep). In the lower panels, we show the ratio of the simulated profiles to the data, with the 1σ scatter depicted with dashed lines. Our best-fit GNFW profiles (see Table 5 and Figure 5) are shown with solid black lines. At all redshifts, the data agree well with the simulated clusters for r 0.1R 500 . In their cores (r < 0.1R 500 ), simulated high-z clusters from both Battaglia et al. (2012) and Dolag et al. (in prep) have significantly higher pressure than the data. There is mild evidence for a steepening of the pressure profile at large radii for z > 0.6, although this may be driven by clumping in the intracluster medium near the virial radius. potential systematic errors and/or biases which may conspire to influence our conclusions.

Comparison to Simulations
In Figure 9 we compare our average pressure profiles to those from simulations, as presented in Battaglia et al. (2012) and Dolag et al. (in prep). For the latter simulations, pressure profiles are computed and presented in Liu et al. (in prep). For each simulation, we show three redshift slices, similar to our low-and high-z subsamples as well as a z ∼ 0 subsample for comparison to P13, and have made mass cuts similar to the SPT 2500 deg 2 survey selection function. All samples have been normalized by f (M ) (see §3.1). We do not show comparisons for CC and NCC subsamples here, since i) we do not have cuspiness measurements for the simulated clusters, and ii) it is unlikely that the simulated clusters will span the full range of properties from cool core to non-cool core clusters. In general, simulations struggle to get the complicated balance between cooling and feedback right in the cores of clusters (r < 0.1R 500 ), but perform well outside of the core where gravity dominates.
As shown in Figure 9, our measured pressure profiles and the pressure profiles from both sets of simulations agree reasonably well at r > 0.1R 500 for all redshift slices. At r < 0.1R 500 , the difference between the simulations and the data becomes worse with increasing redshift. The simulated clusters appear to have massive cool cores in place already at z ∼ 1, while the observed clusters are becoming more centrally concentrated over the past ∼8 Gyr (McDonald et al. 2013). At large radii, the best-fit profile is consistent with Dolag et al. (in prep), and slightly steeper than that predicted by Battaglia et al. (2012), but all profiles are consistent at the 1σ level with the data. We stress that any steepening of the pressure profile may be artificial, indicative of a bias due to clumping of the ICM at higher redshift, a point we will address below.
We can also compare our observations of an unevolving entropy core (Figure 7) to simulations, this time by Gaspari et al. (2012) who focus on the delicate balance between AGN feedback and cooling in the cores of simulated galaxy clusters. These simulations demonstrate that, while at the very center (<10 kpc) of the cluster the entropy can fluctuate significantly (factors of ∼2-3) on short (Myr) timescales, the entropy at 20 kpc is relatively stable of ∼5 Gyr timescales. These simulations, which reproduce realistic condensation rates of cool gas from the ICM, suggest that a gentle, nearlycontinuous injection of mechanical energy from the central AGN is sufficient both to offset the majority of the cooling (preventing the cooling catastrophe) and to effectively "freeze" the entropy profile in place.
Overall, the agreement between observations and simulations is encouraging. The primary difference between the two occurs at r 0.1R 500 , with excess pressure in the simulated cores. As these are the radii where the complicated interplay between ICM cooling, bulk ICM motions, and AGN feedback is most important, it is perhaps unsurprising that the deviations between data and model are most severe in this regime.

Cluster Outskirts: Halo Accretion?
In recent years, a number of different studies have observed a flattening of the entropy profile for a number of different galaxy clusters at the virial radius (e.g., Bautz et al. 2009;Walker et al. 2013;Reiprich et al. 2013;Urban et al. 2014). This flattening, while not observed in all clusters (e.g., Eckert et al. 2013), has been attributed to clumping in the intracluster medium (see e.g., Simionescu et al. 2011;Nagai & Lau 2011;Urban et al. 2014). If a substantial fraction of the ICM beyond the virial radius is in small, overdense clumps, the measured electron density (n e ) over a large annulus will be biased high, due to surface brightness being proportional to n 2 e . These clumps are thought to be the halos of infalling galaxies or small groups. Due to their low mass, they ought to be cooler than the ambient ICM, which could also lead to the measured temperature being biased low.
Given that we measure, on average, lower temperatures (and entropies) at large radii (r R 500 ) in high-z clusters, it is worth discussing whether this result could be driven by clumping and, specifically, how massive these clumps could be. If we assume that an extended source with <20 X-ray counts would go undetected against the diffuse cluster emission, we can estimate a limiting X-ray luminosity at z = 0.8 of L X ∼ 2 × 10 43 erg s −1 , corresponding to a halo mass of M 500 ∼ 8 × 10 13 M ⊙ and temperature of ∼2 keV (Vikhlinin et al. 2009). Thus, it is quite possible that the measured temperature in the outskirts of clusters at z > 0.6 is biased low due to our inability to detect and mask group-sized halos which are in the process of accreting onto the massive cluster.
The entropy flattening that we measure in Figure 7 is driven primarily by the evolution in the temperature profile (Figures 3-4), with only a small, insignificant evolution measured in the outer part of the gas density profile . This makes sense, if the temperature profile is in fact biased by infalling > 10 13 M ⊙ groups at ∼R 500 . Figure 10 illustrates this scenario, showing the density and temperature profiles for a typical SPT-selected cluster (M 500 = 6 × 10 14 M ⊙ , kT 500 = 6.5 keV), and an infalling group-sized system (M 500 = 6 × 10 13 M ⊙ , kT 500 = 1.5 keV). For simplicity, we assume that the infalling group is isothermal and constant density, with ρ g = M g,500 / 4 3 πR 3 500 , where both M g,500 and R 500 can be derived from the group mass, assuming a gas fraction of 0.12. This simple test shows that, at r 1.7R 500 , group-sized halos will significantly bias the measured density high, while at r 1.7R 500 they will bias the measured temperature low. At ∼R 500 , where we measure a flattening of the entropy profile, the density of the infalling group and the ambient ICM are roughly equal, with a factor of ∼3 difference in temperature. This temperature contrast would result in an artificial steepening of the temperature profile, as we observe in Figures 3-4). Following (Vikhlinin 2006), we estimate that the group-sized halos would need to contribute ∼30-40% of the total X-ray counts in the outer annuli to bias the temperature low by the observed 40%, with the exact fraction depending on the relative temperature of the cluster and group.
Simulations suggest that at z 1 there is significantly more massive substructure in the outskirts of galaxy clus- Fig. 10.-Idealized depiction of a group-sized (M 500 = 6 × 10 13 M ⊙ ; blue lines) halo falling into a massive (M 500 = 6 × 10 14 M ⊙ ; red lines) galaxy cluster. The infalling group is assumed to be isothermal and constant density, with the density equal to ρg = 0.12M 500 / 4 3 πR 3 500 and temperature taken from the M-T X relation (Vikhlinin et al. 2009). This figure demonstrates that, as a group-sized halo falls into a massive cluster, it will first significantly bias the density high at r 1.7R 500 (right of dashed vertical line), and then bias the temperature low at r 1.7R 500 (left of dashed vertical line). The latter effect may be driving the steep temperature profile (Figures 3-4) and entropy flattening (Figure 7) that we observe in high-z clusters.
ters. For example, Tillson et al. (2011) find that the accretion rate onto massive halos evolves by a factor of ∼3.5 from z ∼ 1.5 to z ∼ 0, while Fakhouri & Ma (2010) find that 10 14 M ⊙ halos are accreting 10 13 M ⊙ subhalos at a rate ∼3 times higher at z ∼ 1 than at z ∼ 0. These results suggest that the entropy flattening which we measure (Figure 7) is consistent with a temperature bias due to our inability to detect (and mask) large substructures in the outskirts of SPT-selected clusters. We stress that this "superclumping" is qualitatively different than the "clumping" inferred in nearby clusters (e.g., Simionescu et al. 2011;Nagai & Lau 2011;Urban et al. 2014), which is commonly interpreted as large numbers of small subhalos raining onto clusters at the virial radius, where group-sized halos would be detected and masked.

Cool Core Evolution
In an earlier analysis of this dataset , we saw evidence for evolution in the central gas density of cool cores over the past 8 Gyr but no evidence that the minimum ICM entropy in the central ∼10 kpc had evolved since z ∼ 1, maintaining a floor at ∼10 keV cm 2 . Now, with a more rigorous joint-fit analysis to constrain the central temperature, providing a more accurate estimate of the central entropy, we revisit this result.
From Figure 7, we see no measureable evolution in the central entropy bin (0 < r < 0.04R 500 ), from K/K 500 = 0.095 +0.02 −0.02 at low-z to K/K 500 = 0.102 +0.02 −0.01 at low-z. Indeed, the average cool core entropy profile shows no evidence for evolution interior to r < 0.7R 500 since z ∼ 1 (Figure 8). In the absence of feedback or redistribution of entropy, one would expect the average entropy to drop rapidly in the cores of these clusters, on Gyr or shorter timescales. Given the 5 Gyr spanned by this sample, and the consistency with the z ∼ 0 work by Pratt et al. (2010), we can argue that some form of feedback is precisely offsetting cooling between z ∼ 1 and z ∼ 0. Specifically, as the central gas density increases, the core temperature also increases. This trend is contrary to what one would expect from simple hydrostatic equilibrium in a dark matter-dominated halo, but is consistent with the expectation for adiabatic compression of the gas. A likely culprit for this heat injection is radio-mode feedback (e.g., Churazov et al. 2001;Fabian 2012;McNamara & Nulsen 2012), which has been shown to be operating steadily over similar timescales (Hlavacek-Larrondo et al. 2012). Indeed, Gaspari et al. (2011) demonstrate that the immediate result of a burst of AGN feedback is to increase the core temperature of the gas, while leave the large-scale (r 0.1R 500 ) distribution of temperatures unchanged. We finish by stressing that this work and that of McDonald et al. (2013) refer to the entropy in the inner ∼40 kpc as the "central entropy". This annulus, which contains all of the lowest entropy gas falling onto the central cluster galaxy, is limited in size by our relatively shallow exposures. Indeed, Panagoulia et al. (2013) show that with improved angular resolution the entropy continues to drop toward the central AGN. Our discussion of an "entropy floor" is always referring to a fixed radius, within which the mean entropy is not evolving.

Systematic Biases/Uncertainties
Below we briefly address three potential issues with our data analysis: whether the low signal-to-noise in cluster outskirts is driving the entropy flattening, whether joint spectral fitting yields the same results as averaging individual fits, and whether the average temperature profile is mass-dependent.

X-ray Spectrum Signal-to-Noise
While our observing program was designed to obtain 2000 X-ray counts per cluster, a variety of effects conspired to create the scatter in the observed number of counts per cluster (see Figure 2). These factors include uncertainties in the ξ-L X relation, uncertainties on early redshift measurements, and the presence or lack of a cool core. Here, we investigate how strongly the measured average entropy profile depends on the S/N of the included observations.
In Figure 11 we have divided the low-z and highz subsamples by the S/N in the three outermost bins (r > 0.75R 500 ), specifically so that we can test whether the observed entropy flattening is a function of S/N. For the low-S/N subsamples, there are a total of ∼1000 Xray counts in each of the three outermost bins and ∼2700 Fig. 11.-Joint-fit entropy profile for both the low-and high-z subsamples (see also Figure 7). The red and blue points correspond to the joint-fit profiles for low-and high-S/N subsamples, respectively, as described in §4.4.1. We find that, at large radii, the flattening of the entropy profile correlates with both increasing redshift and decreasing S/N. The most significant flattening is present in the high-z, low-S/N subsample, which contains 7 of the counts per radial bin over the full radial range, compared to ∼2800 (outer) and ∼4600 (full radial range) per bin for the high-S/N subsamples. For the low-z clusters, the measured entropy profile appears to be independent of the S/N -the difference of a factor of ∼2 in the total number of X-ray counts used in the spectral modeling does not appear to have a significant affect on the resulting entropy profile. For the high-z clusters, the low-and high-S/N profiles are identical at r < 0.6R 500 , with more flattening at larger radii in the low-S/N clusters. Since the low-S/N clusters also tend to be higher redshift (the high-z, low-S/N subsample contains 7 of the 8 clusters at z > 1 and all 4 clusters at z > 1.1), it is not clear which effect is most responsible for the flattening. In general, there is a trend of more flattening going to both higher redshift and lower S/N. We do not expect a significant bias from low cluster counts, due to our background modeling on an observation-by-observation basis ( §2.2), but we can not rule out this possibility. Given that the lowz, low-S/N clusters have equally low S/N to the high-z, low-S/N clusters, we suggest that the flattening is more significantly driven by redshift evolution.

Joint-Fitting Versus Profile Averaging
To test whether our joint-fitting technique is introducing a systematic bias, we compute individual temperature profiles for our low-z subsample ( Figure 12). Given that each annulus has on the order of ∼ 100 X-ray counts, these individual fits are poorly constrained. However by averaging ∼40 profiles (unweighted), we can constrain the average temperature profile for this subsample. For comparison, we show the results of our joint-fit analysis for the same clusters. We find that the joint-fit method and the averaging method yield consistent results. Since the uncertainty on the joint-fit analysis is really the scatter in the mean for a number of realizations (black points), we have shown the standard error on the mean (standard deviation divided by √ N ) in the average profile (red points) in order to make a fair comparison.
This simple test confirms that our method of jointfitting multiple spectra is largely unbiased with respect to the true average profile. Naively, one might expect a joint-fit analysis to be biased towards the highest signalto-noise spectra, since each cluster is essentially weighted by its total X-ray counts, while each cluster is weighted equally in the averaging method. However, this test shows that any bias that would be imparted by jointfitting spectra of varying signal-to-noise is offset by randomly drawing and fitting subsamples of spectra. Vikhlinin et al. (2006) show that, for a sample of relaxed, low-redshift clusters, low-mass systems tend to have higher central temperatures than their high-mass counterparts. We explore this idea in Figure 13 by dividing our high-z subsample by total gas mass, M g,500 , rather than by cuspiness (the following results hold for the low-z subsample as well). This figure confirms that the temperature profiles of galaxy clusters are not selfsimilar at r 0.3R 500 . We find that low mass systems have temperatures ∼20-30% higher in their cores, consistent with work by Vikhlinin et al. (2006) which covered a larger mass range. At r > 0.3R 500 there appear to be no deviations from self-similarity, suggesting that non-gravitational processes are most likely driving the differences in the core. This figure demonstrates how important a well-selected sample is for such a joint-fit analysis to be successful and yield results representative of the true population. We expect that, given the similar mass distribution of our low-and high-z subsamples (see Figure 2), this mass bias is not driving any of the trends discussed in §3.

Mass Bias
It is also possible that our use of the Y X,500 -M 500 relation to infer R 500 could impart a bias in these results, if the assumed evolution on this relation is incorrect. To investigate this potential bias, we re-extracted spectra using R 500 estimates based on the M gas -M 500 relation, and repeated the analysis described in §2.2. The resulting temperature profiles were consistent with what we have presented here, suggesting that our assumed evolution on the Y X,500 -M 500 relation is appropriate out to z ∼ 1.2. 2). Individual cluster profiles are shown as red dashes, while the average of these profiles is shown as thick red points. The uncertainty shown for the average profile is the standard error on the mean (standard deviation divided by √ N ) to allow a better comparison to the joint-fit uncertainties, which are measuring the scatter in the mean temperature for a large number of realizations. The joint-fit result, which is fully consistent with the average profile, is shown in black. This figure demonstrates that our joint-fit analysis is not strongly affected by combining spectra of varying signal-to-noise. Fig. 13.-Average temperature profiles for high-z clusters. We show the combined fits in grey, low-mass systems in blue, and highmass systems in red. This figure demonstrates that the deviation from self-similarity interior to 0.3R 500 , consistent with earlier work by Vikhlinin et al. (2006), is present out to z ∼ 1. Beyond 0.3R 500 , there is no evidence for a mass bias.

SUMMARY
We have presented a joint-fit analysis of X-ray spectra for 80 SPT-selected galaxy clusters spanning redshifts 0.3 < z < 1.2. These spectra, which individually only contain ∼2000 X-ray counts, are divided into subsamples of ∼20 clusters each, and the spectra in each subsample are simultaneously modeled assuming a self-similar temperature profile. This allows us to constrain the redshift evolution of the temperature, pressure, and entropy pro-files for massive clusters. Our major results are summarized below: • We are able to constrain the average temperature profile out to ∼ 1.5R 500 for both low-z (0.2 < z < 0.6) and high-z (0.6 < z < 1.2) clusters. The temperature profile for our low-z subsample is consistent with earlier works by Vikhlinin et al. (2006), Pratt et al. (2007), and Baldi et al. (2012). Combined with density profiles from McDonald et al. (2013), we constrain the pressure and entropy profiles over 0.01R 500 < r < 1.5R 500 , providing the first constraints on the redshift evolution of the Universal pressure profile.
• The cores of high-z cool core galaxy clusters appear to be marginally (∼2σ) cooler than those of their low-z counterparts by ∼30%. This is precisely what is needed to maintain constant central entropy since z ∼ 1, given the observed evolution in the central electron density, as reported by McDonald et al. (2013).
• The average temperature profile in the outskirts of high-z cool core clusters is steeper than in the outskirts of low-z cool core clusters. This results in a steepening of the outer pressure profile and a flattening of the outer entropy profile. These data are consistent with an increase in the number of group-mass (∼1.5 keV) halos falling into the cluster at R 500 which our relatively shallow exposures are unable to detect. This "superclumping" should be a factor of ∼3 times more common at z ∼ 1 than it is today. Failure to mask these massive subhalos can bias the temperature at R 500 low by the observed amount (∼40%).
• The cores of low-z clusters have significantly higher pressure than those of high-z clusters, increasing by a factor of ∼10 between z ∼ 1 and z = 0. This is driven primarily by the increase in central density with decreasing redshift , but is also boosted by the increasing central temperature with decreasing redshift.
• We find good overall agreement between our low-z average pressure profile and those of Arnaud et al. (2010) andPlanck Collaboration et al. (2013b).
• Simulated clusters from Battaglia et al. (2012) and Dolag et al. (in prep) reproduce the evolution of the observed pressure profile at r 0.1R 500 . The growth of cool cores, resulting in a factor of ∼10 increase in the central pressure over the past ∼8 Gyr is not reproduced in simulations.
• We measure no significant redshift evolution in the entropy profile for cool cores at r 0.7R 500 , suggesting that the average entropy profile for massive clusters is stable on long timescales and over a large range of radii. This may be a result of a long-standing balance between ICM cooling and AGN feedback.
This work demonstrates that a joint-spectral-fit X-ray analysis of low signal-to-noise cluster observations can be used to constrain the average temperature, pressure, and entropy profile to large radii. This has proven to be a powerful method for analyzing high-redshift clusters, where obtaining >10,000 X-ray counts per cluster is unfeasible. These techniques will add additional power to future surveys by, for example, eRosita, or serendipitous surveys like ChaMP (Barkhouse et al. 2006), XCS (Mehrtens et al. 2012), and XXL (Pierre et al. 2011), where the number of clusters with data is high, but the data quality per cluster is low.
A. AVERAGE GAS DENSITY, TEMPERATURE, ENTROPY, AND PRESSURE PROFILES Below, we provide all of the three-dimensional deprojected quantities used in this work. Details on the derivation of these quantities, which are based on the two-dimensional temperature profiles (Table 3) and the three-dimensional gas density profiles from McDonald et al. (2013), can be found in §2 and §3. All uncertainties are 1σ. Deprojected average temperature, entropy, and pressure profiles for our sample of 80 SPT-selected clusters.