Issue 
A&A
Volume 617, September 2018



Article Number  A64  
Number of page(s)  15  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201732458  
Published online  18 September 2018 
Resolving the hydrostatic mass profiles of galaxy clusters at z ∼ 1 with XMMNewton and Chandra
^{1}
IRFU, CEA, Université ParisSaclay, 91191 GifSurYvette, France
email: iacopo.bartalucci@cea.fr
^{2}
Université Paris Diderot, AIM, Sorbonne Paris Cité, CEA, CNRS, 91191 GifsurYvette, France
Received:
13
December
2017
Accepted:
7
May
2018
We present a detailed study of the integrated total hydrostatic mass profiles of the five most massive M_{500}^{SZ} < 5 × 10^{14} M_{⊙} galaxy clusters selected at z ∼ 1 via the Sunyaev–Zel’dovich effect. These objects represent an ideal laboratory to test structure formation models where the primary driver is gravity. Optimally exploiting spatiallyresolved spectroscopic information from XMMNewton and Chandra observations, we used both parametric (forward, backward) and nonparametric methods to recover the mass profiles, finding that the results are extremely robust when density and temperature measurements are both available. Our Xray masses at R_{500} are higher than the weak lensing masses obtained from the Hubble Space Telescope (HST), with a mean ratio of 1.39_{−0.35}^{+0.47}. This offset goes in the opposite direction to that expected in a scenario where the hydrostatic method yields a biased, underestimated, mass. We investigated halo shape parameters such as sparsity and concentration, and compared to local Xray selected clusters, finding hints for evolution in the central regions (or for selection effects). The total baryonic content is in agreement with the cosmic value at R_{500}. Comparison with numerical simulations shows that the mass distribution and concentration are in line with expectations. These results illustrate the power of Xray observations to probe the statistical properties of the gas and total mass profiles in this high mass, highredshift regime.
Key words: Xrays: / galaxies: clusters – dark matter
© ESO 2018
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
In the current ΛCDM paradigm, structure formation in the Universe is driven by the gravitational collapse of the dark matter component. In this context, the form of the dark matter density profile is a sensitive test not only of the structure formation scenario, but also of the nature of the dark matter itself. In addition, it is impossible to fully comprehend the baryonic physics without first achieving a full understanding of the dominant dark matter component.
Cosmological numerical simulations uniformly predict a quasiuniversal cusped dark matter density profile, whose form only depends on mass and redshift. Perhaps the bestknown parameterisation of dark matter density profiles is the NavarroFrenkWhite (NFW) profile suggested by Navarro et al. (1997).
This profile is flexible; in scaled coordinates (i.e. radius scaled to the virial radius) its shape is characterised by a single parameter, the concentration c, the ratio of the scale radius to the virial radius, r_{s}/R_{Δ}^{1}. Its normalisation, for a given concentration, is proportional to the mass. The concentration is known to exhibit a weak dependence on mass and redshift (typically a decrease of a factor 1.5 at z = 1, e.g. Duffy et al. 2008), although the exact dependence is a matter of some debate in the literature (e.g. Diemer & Kravtsov 2015).
In the local (z ≲ 0.3) Universe, there is now strong observational evidence for NFWtype dark and total matter density profiles with typical concentrations in line with expectations from simulations. Such evidence comes both from Xray observations (e.g. Pointecouteau et al. 2005; Vikhlinin et al. 2006; Buote et al. 2007), and more recently, from gravitational lensing studies (e.g. Merten et al. 2015; Okabe & Smith 2016). While encouraging, more work is needed to make the different observations converge, and observational biases and selection effects are still an issue (e.g. Groener et al. 2016).
In contrast, constraints on distant systems, and the evolution to the present, are sparse. The recent compilation of weak and strong lensing observations of 31 clusters at z > 0.8 by Sereno & Covone (2013) illustrates the difficulty of obtaining firm constraints on cluster mass profiles in this redshift regime with lensing (their Fig. 1). Stacking the velocity data of ten clusters in the redshift range 0.87 < z < 1.34, Biviano et al. (2016) derived a concentration , in agreement with theoretical expectations. Perhaps the strongest constraints come from the Xray observations of Schmidt & Allen (2007, 0.06 < z < 0.7) and Amodeo et al. (2016, 0.4 < z < 1.2). The evolution factor of these c–M relations, expressed as (1+z)^{α}, is consistent with theoretical expectations, but with large uncertainties (α = 0.71 ± 0.52, and α = 0.12 ± 0.61, respectively).
Fig. 1 3D temperature profiles of all the clusters of our sample. Radii are scaled by . Section 2.3.2 describes the temperature profile calculation. For each panel: the black points represent the nonparametriclike 3D temperature profiles measured using Chandra and XMMNewton, with round and squared points, respectively. The grey shaded area represents the bestfitting 3D parametric model (Vikhlinin et al. 2006). The blue and red areas represent the result of the backward fit (BP approach, see Sect. 2.4.1), assuming hydrostatic equilibrium and an NFW or an Einasto mass profile, respectively. The parametric models were estimated only in the radial range covered by the density profile. The shaded regions correspond to the 68% confidence level regions. 
The poor constraints at high redshift are due in part to the difficulty in detecting objects at these distances. Surveys using the Sunyaev–Zel’dovich (SZ) effect have the advantage of the redshift independent nature of the signal and the tight relation between the signal and the underlying total mass (da Silva et al. 2004). The advent of such surveys (Planck Collaboration XXIX 2014; Hasselfield et al. 2013; Bleem et al. 2015; Planck Collaboration XXVII 2016; Hilton et al. 2018) has transformed the quest for highredshift clusters. Samples taken from such surveys are thus ideal for testing the theory of the dark matter collapse and its evolution. In this context, Xray observations, while not the most accurate for measuring the mass because of the need for the assumption of hydrostatic equilibrium (HE), can give more precise results than other methods because of their good spatial resolution and signaltonoise ratios. A combination with theoretical modelling can give crucial insights into both the dark matter collapse and the coeval evolution of the baryons in the potential well.
Here we present a pilot study of the Xray hydrostatic mass profiles of the five most massive SZdetected clusters at z ∼ 1, where the mass is M_{500} > 5 × 10^{14} M_{⊙} as estimated from their SZ signal. Initial results, obtained by optimally combining spatially and spectrally resolved XMMNewton and Chandra observations, concerned the evolution of gas properties, and were described in Bartalucci et al. (2017; hereafter B17). Here we used the same observations to probe the total mass and its spatial distribution. We discuss the various Xray mass estimation methods used in Sect. 2, and the robustness of the recovered mass distribution in Sect. 3. Results are compared with local systems to probe evolution in Sect. 4 and cosmological numerical simulations in Sect. 5. We discuss our conclusions in Sect. 6.
We adopt a flat Λcold dark matter cosmology with Ω_{m} = 0.3, Ω_{Λ} = 0.7, H_{0} = 70 km Mpc s^{−1}, and h(z) = (Ω_{m}(1 + z)^{3} + Ω_{Λ})^{1/2} throughout. Uncertainties are given at the 68% confidence level (1σ). All fits were performed via χ^{2} minimisation.
2. Data sample and analysis
2.1. Sample
A detailed description of the sample used here, including the data reduction, is given in B17. Briefly, the sample is drawn from the South Pole Telescope (SPT) and Planck SZ catalogues (Bleem et al. 2015; Planck Collaboration XXVI 2011), and consists of the five galaxy clusters with the highest SZ mass proxy value^{2} () at z > 0.9 (see Fig. 1 of B17). All five objects were detected in the SPT survey; PLCK G266.6+27.3 was also independently detected in the Planck SZ survey. All five have been observed by both XMMNewton and Chandra, using the European Photon Imaging Camera (EPIC, Turner et al. 2001 and Strüder et al. 2001) and the Advanced CCD Imaging Spectrometer (ACIS, Garmire et al. 2003), respectively. Four objects were the subject of an XMMNewton Large Programme, for which the exposure times were tuned so as to enable extraction of temperature profiles up to R_{500}. Shorter archival Chandra observations were also used. The fifth object, PLCK G266.6+27.3, was initially the subject of a snapshot XMMNewton observation (Planck Collaboration XXVI 2011), and was then subsequently observed in a deep Chandra exposure.
Dedicated pipelines, described in full in B17, were used to produce cleaned and reprocessed data products for both observatories. These pipelines apply identical background subtraction and effective area correction techniques to prepare both XMMNewton and Chandra data for subsequent analysis. The definition of surface brightness and temperature profile extraction regions was also identical, and point source lists were combined.
2.2. Preliminaries
Under the assumptions of spherical symmetry and HE, the integrated mass profile of a cluster is given by (1) where μ = 0.6 is the mean molecular weight in a.m.u^{3}, m_{H} is the hydrogen atom mass, and T(r) and n_{e}(r) are the 3D temperature and density radial profiles, respectively. The key observational inputs needed for this calculation are thus the radial density and temperature profiles, plus their local gradients. A complication is that these quantities are observed in projection on the sky, and thus the binaveraged 2D annular (projected) measurements must be converted to the corresponding measurements in the 3D shell (deprojected) quantities.
A number of approaches exist in the literature for the specific case of cluster mass modelling (for a review, see e.g. Ettori et al. 2013, and references therein). Generally speaking, one can either model the mass distribution and fit the projected (2D) quantities (backwardfitting), or deproject the observable quantities to obtain the 3D profiles and calculate the resulting integrated mass profile (forwardfitting). This deprojection in turn can either be performed either by using parametric functions or be undertaken nonparametrically.
In the following, we chose to calculate all deprojected quantities at the emissionweighted effective radius, r_{w}, assigned to each projected annulus, i, defined as in Lewis et al. (2003): (2) Formally, r_{w} should be calculated iteratively from the density profile, but McLaughlin (1999) has shown that the above equation is an excellent approximation for a wide range of density profile slopes.
2.3. Density and temperature profiles
2.3.1. Density
We used the combined XMMNewton–Chandra density profiles detailed in B17, which were derived from the [0.3–2] keV band surface brightness profiles using the regularised nonparametric deprojection technique described in Croston et al. (2006). As shown in B17, the resulting 3D (deprojected) density distributions from XMMNewton and Chandra agree remarkably well.
We then fitted these profiles simultaneously with a parametric model based on that described in Vikhlinin et al. (2006, see Appendix A), allowing us to obtain for each object a combined density profile that fully exploits the high angular resolution of Chandra in the core and the large effective area of XMMNewton in the outskirts. The resulting 3D density distribution is technically a parametric profile. However, in view of the much better statistical quality of the density profiles (compared to that of the temperature profiles), this last parametric step does not overconstrain the resulting mass distribution. As in B17, to avoid extrapolation, the minimum and maximum radii for the parametric models were set to match those of the measured deprojected profiles.
SPT–CL J0546–5345 presents a clear substructure in its southwest sector which was not masked in B17. Since here our focus is on the measurement of integrated mass profiles, such substructures should generally be excluded from the analysis. We thus computed a new combined density profile with the substructure masked for this system. The new profile we use in this work is described in Appendix B and is shown in Fig. B.1.
2.3.2. Temperature
We base our 3D (deprojected) temperature profiles on those published in B17. In a first step, we extracted spectra from concentric annuli centred on the Xray peak and determined the 2D (projected) temperature profile by measuring the temperature in each bin. We iteratively modified the annular binning scheme defined in B17, to ensure that the fell within the outermost radius of the final annulus of each profile.
We then employed two methods to obtain the 3D temperature profile:
Parametric: we fitted a model similar to that proposed by Vikhlinin et al. (2006), reducing the number of free parameters when necessary, to the 2D profiles. This model was convolved with a response matrix to take into account projection and (for XMMNewton) PSF redistribution; during this convolution, the weighting scheme proposed by Vikhlinin (2006, see also Mazzotta et al. 2004) was used to correct for the bias introduced by fitting isothermal models to a multitemperature plasma. Uncertainties were computed via 1000 Monte Carlo simulations of the projected temperature profiles.
Nonparametriclike: analytical models such as those described above tend to be overconstrained, and do not reflect the fact that the temperature distribution is measured only at the points corresponding to the limited number of annuli within which spectra are extracted. To overcome these limitations we define the nonparametericlike temperature profile by estimating the parametric model temperature at the weighted radii corresponding to the 2D annular binning scheme, and imposing the uncertainty on the annular spectral fit as a lower limit to the uncertainty in the 3D bin.
The resulting profiles are shown in Fig. 1, where the smooth grey envelope represents the parametric 3D temperature distribution, and the black points with errors represent the 3D nonparametriclike temperature profile. 3D density and temperature profiles are available in electronic format.
2.4. Mass profiles
2.4.1. Mass profile calculation
Total mass profiles were determined following Eq. (1). To examine the robustness of the recovered profiles, we used both forwardfitting and backwardfitting methods, as we describe below.

Forward nonparametriclike (FNPL): this is our baseline mass measurement. It uses the combined 3D density profile (Sect. 2.3.1) and the nonparametriclike 3D temperature profile as input, and produces a mass profile estimate at each weighted radius, r_{w}. The mass measurement and its uncertainty were calculated using a similar scheme to that first presented in Pratt & Arnaud (2003) and further developed in Démoclès et al. (2010). In this procedure, a random temperature was generated at each r_{w}, and a cubic spline was used to compute the derivative. One thousand Monte Carlo simulations of this type were performed; the final mass profile and its uncertainties were then derived from the median and associated 68% confidence region. The mass profiles derived from these realisations were constrained to respect the monotonic condition (i.e. M(r + dr) > M(r)) and to be convectively stable (i.e. d ln T/d ln n_{e} < 2/3).
The resulting mass profiles are shown with their corresponding error bars in Fig. 2. The relative errors are of the order of 30% in the inner core, and (somewhat counterintuitively) decrease to ∼10 −15% at large radii. This effect is an intrinsic property of the typical amplitude and uncertainty on the logarithmic density and temperature gradients, and is quantified in more detail in Appendix C.

Forward parametric (FP): here the fully parametric 3D density and temperature profiles were used to compute the total mass distribution on the radial grid of the combined density profile. Uncertainties were calculated using 1000 Monte Carlo realisations, and we did not impose any condition on the resulting mass profiles. The grey shaded areas in Fig. 2 correspond to the 68% dispersion envelopes.
This method may lead to nonphysical results, as can be seen at large radii in SPT–CL J2146–4633 and SPT–CL J0546–5345, where the cumulative total mass profiles start decreasing. For this reason, we do not compute a median profile and we do not use these results to perform quantitative analyses. However, these profiles retain the maximum amount of information on the intrinsic dispersion, allowing us to explore the dispersion related to density and temperature measurement errors. Additionally, these mass profiles are estimated on the finer radial grid and wider radial range of the density profiles and so they can be used to qualitatively investigate the behaviour in the cluster core and outskirt regions. We note that we did not extrapolate the parametric model of the density profiles, i.e. we did not attempt to estimate masses in regions where there are no observational constraints.
Backward parametric (BP): here we assumed that the total mass distribution could be described by an NFW (Navarro et al. 1997) or Einasto (Einasto 1965; Navarro et al. 2004) distribution, and inverted Eq. (1), taking into account the 3D density profile, to obtain the corresponding 3D temperature profile. This was then projected and convolved with the instrument response and PSF, and fitted to the 2D temperature profile. Uncertainties were estimated through a Monte Carlo randomisation procedure using 1000 realisations. The resulting temperature and mass profiles are shown in Figs. 1 and 2. The analysis was again restricted to the radial range covered by the density profile.
Direct fit (DF): we also directly fitted the FNPL mass profiles using the NFW and Einasto functional forms. The resulting best fits, computed on the combined density profile radial grid, are shown in Fig. 2. The corresponding uncertainties were estimated by repeating the fitting procedure on 1000 Monte Carlo realisations of the FNPL mass profile. The NFW fit concentrations at R_{500}, c_{500} ≡ R_{500}/r_{s} where r_{s} is the scale radius, are given in Table 1.
Fig. 2 Scaled mass profiles of all the clusters of our sample, derived assuming hydrostatic equilibrium (HE). All calculations used the combined XMMNewtonChandra density profile, and full details of the mass calculation methods are given in Sect. 2.4.1. Various methods give very consistent results in the radial range with temperature information, but may diverge at small and large radius in spite of the density information. For each panel: the black points represent the mass profiles obtained from the forward nonparametriclike method, using the HE equation and the nonparametriclike temperature profiles shown as black points in Fig. 1. The blue and red solid lines represent the fit of these forward nonparametriclike profiles using a NFW and an Einasto model, respectively. The grey area is the mass profile computed assuming hydrostatic equilibrium and using the parametric temperature profiles shown with a grey area in Fig. 1. The blue and red envelopes represent the mass profile computed using the backward method, i.e. fitting the observed temperature profile with a model derived from the HE equation and assuming a NFW and an Einasto profile, respectively, for the underlying total mass distribution. The parametric mass profiles are estimated in the wider radial range covered by the density profile. 
Relevant quantities computed at fixed radii and overdensities.
2.4.2. Determination of mass at fixed radius and density contrast
The value of (and consequently ) was determined iteratively using the M_{500}–Y_{X} relation, as calibrated in Arnaud et al. (2010), assuming selfsimilar evolution. Here Y_{X} is defined as the product of the gas mass computed at and the temperature measured in the [0.15−0.75] region (Kravtsov et al. 2006). As the radial density bin widths used here differ from those used in B17, as described above in Sect. 2.3, the gas mass profiles and the quantities based on M_{500}−Y_{X} were updated. For this reason, the values in Table 1 differ slightly (∼1%) from those published in Table 2 of B17.
We determined the FNPL masses at density contrasts Δ = [2500,500], namely and , at radii and , respectively. We also interpolated all the mass profiles (except FP) described in Sect. 2.4, at . These are referred to as in the following text and figures. Radii and the corresponding masses are given in Table 1.
The of SPT–CL J2106–5844 reaches a nonphysical value of ∼40 × 10^{14} M_{⊙}. The is at the outer edge of the last temperature bin, so extrapolation is required. As the mass profile of this object is very steep, the radius at which Δ = 500 is boosted, and the corresponding mass reaches nonrealistic values. The resulting and estimates are provided in Table 1, although they are not used for any quantitative analysis. The DF NFW yields more reasonable M_{500} estimates, although they are poorly constrained. The bestfit c_{500} value is equal to the minimum value allowed by the fit (c_{500} = 0.01), corresponding to the quasipower law behaviour of the mass profile, and yields an M_{500} that is significantly greater than . A more conservative lower value of c_{500} = 1 forces the curve to be higher in the core and the fit is then driven by the third point (R ∼ R_{2500}) because of its small relative error. This analysis yielded a ∼5 times higher χ^{2} and a value of M_{500} = 7.6 ± 2.1 × 10^{14} M_{⊙}, now in agreement with . This result must simply be considered as an NFW extrapolation, with priors on c_{500}, of the welldetermined mass at R_{2500}.
Two objects from our sample, SPT–CL J0546–5345 and SPT–CL J2106–5844, were also analysed by Amodeo et al. (2016) using Chandra only datasets. The authors estimated M_{200} and c_{200} using the BP approach and the NFW functional form. Using the concentration and mass values published in their Table 2 to compute M_{500} yields M_{500} = 4.0 ± 2.9 × 10^{14} M_{⊙} and M_{500} = 6.5 ± 3.9 × 10^{14} M_{⊙} for SPT–CL J0546–5345 and SPT–CL J2106–5844, respectively. These are perfectly consistent with the present BPNFW estimates; however, our deeper observations and extended radial coverage allowed us to better constrain the measurements, the relative errors being ∼5 times smaller.
3. Robustness of Xray mass
In this section, we first examine the robustness of the HE mass estimate to the Xray analysis method. As the HE assumption is a known source of systematics, through the HE bias, we then compare the HE mass to lensing mass estimates, which do not rely on this assumption.
3.1. Mass profile shape
Figure 2 shows the mass profiles resulting from the different mass estimation methods discussed above. The BP results indicate that while the NFW model is a good description in the case of relaxed objects (e.g. PLCK G266.6+27.3) and some perturbed systems (e.g. SPT–CL J2341–5119), the Einasto model is generally a better fit for our sample (as is evident from the figure, and from the χ^{2} value) and is more able to fit a wider range of dynamical states. This is unsurprising given the larger number of parameters in the Einasto model. Forward and backward methods also give extremely consistent results. The limitations of the NFW model can be seen in SPT–CL J0546–5345, where this form is clearly a poor description of the data, leading to the BP NFW masses being somewhat different to those from other methods.
Overall, all the mass estimation methods yield remarkably robust and consistent results within the radial range covered by the spectroscopic data, i.e. within the minimum and maximum effective radii of the temperature profile bins, except in cases where the underlying model is insufficiently flexible. Mass profile uncertainties are quite different between methods, however, with the FP method yielding the smallest and the FNPL method yielding the largest (or most conservative). This simply reflects the restrictions each method places on the possible shape of the profile.
Outside the radial range covered by the spectroscopic data, the results are most robust and agnostic to the mass estimation method when the profiles are regular and can be described by a simple model (e.g. NFW). However, when the radial sampling is poor (the profiles have few points) or when the profile is irregular (e.g. SPT–CL J0546–5345), estimation of the mass outside the radial range probed by the spectroscopic data is less robust and will depend strongly on the method used to measure the mass. In addition, outside the region covered by the spectroscopic data, the uncertainties rapidly increase with the distance from effective radius of the final temperature measurement, in spite of the density information.
3.2. Mass within
We now turn to the robustness of the mass determined within a fixed radius, calculated as described in Sect. 2.4.2. Figure 3 shows the ratio between the mass obtained employing the different methods, with as a reference mass. For SPT–CL J2146–4633 and PLCK G266.6+27.3 the HE mass measurements are in excellent agreement, the difference being within a small percent. SPT–CL J2341–5119 and SPT–CL J0546–5345 present larger differences (∼10%) according to the mass estimation method between mass estimates. Interestingly, the BP masses of SPT–CL J2341–5119 are closest of all the objects to its .
Fig. 3 Comparison of the hydrostatic mass computed at fixed radius, , using the different methods, in units of . There is excellent agreement, with differences of less than 10%, when the radius is enclosed in the radial range covered by the spectroscopic data. 
SPT–CL J2106–5844 is the only cluster for which all the methods yield masses greater than , by a factor of ∼40%, except if we further restrict the possible range of concentration parameters. The difference between mass estimates is also noticeably larger than for the other objects, due to the limited radial coverage. Even if the masses are estimated at a fixed radius, , this radius falls barely within the outermost temperature radial bin. We conclude that that in order to perform robust measurements, the radius at which the mass is to be estimated should lie within the weighted radial range covered by the spectroscopic data.
3.3. Comparison to weak lensing
Weak lensing mass measurements represent an additional and independent method of investigating the robustness of our mass determinations; furthermore, understanding the systematic differences between weak lensing and Xray masses at z ∼ 1 is crucial for any future cosmological or physical exploitation of such samples. We compared our results with the weak lensing masses published in Schrabback et al. (2018), who determined M_{500} for 13 SPT clusters observed with the Hubble Space Telescope. Four of their objects are in common with our sample.
Schrabback et al. (2018) give different weak lensing M_{500} estimates, depending on the choice of centre (Xray peak and SZ peak). The top left panel of Fig. 4 show the comparison between M_{500} measured using the Xray peak as centre, , and , as listed in Table 1. Formally, there is good agreement, with the masses for each individual cluster being consistent at 1σ. However, there is a clear systematic offset in the sense that all Xray masses are higher than the WL masses, with an errorweighted mean ratio of . The right panel shows the comparison with the HE masses, computed at (instead of ) to avoid an artificial increase in differences due to different apertures. The difference is similar to . We found the same results by comparing the Xray masses with the weak lensing masses centred on the SZ peak, , as shown in the bottom panels of Fig. 4.
Fig. 4 Comparison between our Xray masses and the weak lensing masses published in Schrabback et al. (2018). All estimates for a given cluster are consistent within the statistical errors. However, there is a general trend of smaller lensing mass than the HE mass, contrary to expectation. Note also the higher statistical precision of the Xray masses. Top left panel: comparison between and weak lensing masses estimated at R_{500}, centred on the Xray peak. The grey area for the weak lensing represents the statistical errors. The black solid bars represent the sum in quadrature of systematic and statistical errors. The blue and red lines represent the bias, (1 − b), between the Xray hydrostatic and weak lensing mass as measured by Weighting the Giants (WtG, von der Linden et al. 2014) and by the Canadian Cluster Comparison Project (CCCP, Hoekstra et al. 2015), respectively. To better visualise the points, we crop the lower values of SPTCLJ23415119 and SPTCLJ05465345, which are of the order of ∼10^{−1} × 10^{14} M_{⊙}. Top right panel: same as the top left panel, except showing the comparison between the hydrostatic mass computed at , , and the weak lensing masses. Bottom left and right panels: same as the top panels except that weak lensing masses are computed using the SunyavezZeldovich (SZ) peak as the centre. The error–weighted mean ratio and corresponding errors are reported in each panel. 
This result is unexpected. The socalled “hydrostatic bias”, owing to the assumption of HE, is believed to result in a net underestimate of the total mass in Xray measurements, while lensing observations, although slightly biased, are expected to yield results that are closer to the true value. Indeed, such a trend has been found, for example in the Weighing the Giants (WtG, von der Linden et al. 2014) project and by the Canadian Cluster Comparison Project (CCCP, Hoekstra et al. 2015), where the Xray hydrostatic masses are ∼30% and 20% lower than the WL values^{4}, respectively. While there are only four objects in our sample, we find the opposite trend here. With a HEtoWL mass ratio of , our results are marginally consistent at 1σ with the Schrabback et al. (2018) results, and inconsistent with WtG at the 2σ level.
This comparison underlines the capability and complementarity of Xray observations with respect to optical observations, especially at these redshifts. The Xray statistical errors are significantly smaller than the weak lensing uncertainties; furthermore, the Xray results are remarkably robust, as we demonstrate in the previous sections. The results we find here show that while Xray observations at high redshift are expensive and challenging, they offer a robust and precise tool which can efficiently complement measurements in other wavelengths.
4. Evolution of cluster properties
4.1. – relation and evolution of the ratio
The M_{500}−Y_{X} relation we use in this work was calibrated using hydrostatic masses derived from the relaxed subsample of 12 REXCESS objects (Arnaud et al. 2010; Pratt et al. 2010), plus eight additional relaxed systems from Arnaud et al. (2007). This relation was derived from local objects and we assumed selfsimilar evolution. The present observations offer the opportunity to investigate the robustness of this relation when applied to a highredshift sample dominated by disturbed objects. Figure 5 shows the resulting comparison of with and , in the left and right panels, respectively.
Fig. 5 Left panel: comparison between computed iteratively through the M_{500}−Y_{X} relation and the hydrostatic mass, . The estimates are consistent within the statistical errors. The of SPT‒CL J2106‒5844 is ∼3.5 greater than (see Sect. 2.4.2). For this reason, the point is off the scale and its is instead shown with the black arrow. The black dotted line is the 1 : 1 relation. Right panel: same as the left panel, except showing the comparison between and . Error–weighted mean ratio and corresponding errors are reported in each panel. 
In both cases there is excellent agreement between individual measurements. The only exception is the of SPT–CL J2106–5844, which is subject to the systematic uncertainty discussed above. The errorweighed mean ratios are and , consistent with unity. This suggets that the relation is robust, even when applied to such an extreme sample. The good agreement is consistent either with no evolution of the ratio between the two quantities, or with an evolution of the ratio where the evolution is counterbalanced by some other effects. However the latter explanation is unlikely given that the evolution is perfectly compensated, such that the agreement between and is excellent as a function of redshift. We note that this result does not necessarily imply that there is no evolution of the bias between the hydrostatic mass and the true mass. However, our comparison with weak lensing above would suggest that the bias cannot be dramatic.
4.2. Scaled mass and total density profiles
We calculated the total density profiles for our sample using the bestfitting DF Einasto model. The resulting profiles are shown compared to those from REXCESS in the left panel of Fig. 6. The right panel shows the corresponding cumulative total mass profiles.
Fig. 6 Left panel: scaled total density profiles computed using the mass profiled derived from the DF Einasto model. For each cluster, the total radial range is that of the combined density profile; estimates beyond the radial range covered by temperature measurements are marked with dotted lines. The grey lines represent the scaled total density profiles derived from the REXCESS sample. Right panel: scaled mass profiles. The colour scheme is the same as in the left panel. The black error bars in both panels represent the 68% dispersion of the REXCESS profiles at 0.1 and 0.5 . 
A consistent picture emerges from these comparisons. At all the mass profiles are in excellent agreement: the dispersion is similar, and interestingly is centred around unity (i.e. the is comparable to , consistent with our findings in the previous section). Apart from the profile of SPT–CL J2106–5844, which is affected by poor radial coverage especially at large radius, all the density profiles of the z ∼ 1 sample lie within the envelope of the REXCESS profiles at high radii (>0.5R_{500}). However, the profiles tend to be shallower on average in the central regions. The profiles of SPT–CLJ2146–4633 and SPT–CL J2106–5844 are even shallower in the core (<0.3R_{500}) than the leastpeaked REXCESS profile. The other three systems lie within the 1σ dispersion of REXCESS profiles, but tend to trace the lower envelope of the distribution in total density and total mass, especially towards the most central parts.
Unfortunately, due to the small size of the sample and the poor quality of SPT–CL J210–5844, we cannot quantify whether there is a significant difference compared to REXCESS in median mass profile shape and/or an increase in the intrinsic scatter around it. If these differences are confirmed, this behaviour can be interpreted either as evolution in the core regions or as being due to a difference between Xray and SZ selection. Comparison with an Xray selected sample at similar redshifts or comparison to a similar SZselected sample at lower redshift, would help to clarify this point.
4.3. Sparsity
The halo sparsity was introduced by Balmès et al. (2014) to characterise the form of the mass distribution in a way that is independent of any parametric model. It is defined as the ratio of masses integrated within two fixed overdensities, (3) where Δ_{1,2} represent the overdensities at which the masses are calculated, with Δ1 < Δ2. As is discussed in Balmès et al. (2014), the properties of the sparsity are independent of the choice of the Δ as long as the definition of the halo is not ambiguous (Δ_{1} not too small), and that dynamical interaction between baryons and dark matter can be neglected (Δ_{2} not too large).
We chose to measure the sparsity within overdensity of Δ_{1} = 500 and Δ_{2} = 2500 with respect to the critical density. These overdensities are well matched to the sensitivity of the Xray observations discussed here, and are sufficiently distant to properly sample the form of the mass profile.
Figure 7 shows the resulting sparsity measurements for our z ∼ 1 sample. These data are compared to those from REXCESS, which exhibit a peaked distribution in a narrow range, 1 < S < 3.
Fig. 7 Number of clusters as a function of their sparsity. The blue and gold shaded bins represent the sparsity distribution of the five highz clusters and the REXCESS sample, respectively. Individual objects are identified by symbols overplotted on the blue bins. The grey arrows represent the lower limit of the sparcity of SPT–CL J2146–4633 and SPT–CL J2106–5844. 
Three of the clusters in our z ∼ 1 sample have sparsity values that lie well within the REXCESS distribution. SPT–CL J2146–4633 and SPT–CL J2106–5844 lie outside this distribution, their sparsity being ∼4–5 times the mean value (∼2) compared to REXCESS. This result reflects what we already found for the mass profiles. This study and the recent parallel study of Corasaniti et al. (2018) represent the first applications of this quantity to a large sample of objects. The narrow distributions in Fig. 7 show its effectiveness in tracing the population characteristics.
4.1. Baryon fraction
The baryon fraction determined at the radius R is defined as (4) where M_{star} is the total stellar mass, M_{gas} is the gas mass, M_{tot} is the halo total mass, and all quantities are integrated within R. In B17 we presented the baryon fraction derived using for the total mass estimate. Here we extended this analysis by deriving the baryon fraction using for the M_{tot} term. This is fundamental to understand possible systematics related to the fact that the gas mass profiles and measurements are correlated (i.e. the Y_{X} is based on the gas mass). Figure 8 shows the baryon fraction as a function of mass computed for this work, the results from B17, and the mean derived from REXCESS (Pratt et al. 2009). For M_{star} at R_{500} we used the stellar masses published in Chiu et al. (2016, stellar mass for SPT–CL J2146–4633 is not available).
Fig. 8 Baryon fractions computed at as a function of mass. The baryon fraction does not show any dependence with respect to the mass at this high z. Black points represent the baryon fraction computed using the hydrostatic mass profiles at . The grey points represent the baryon fractions published in Bartalucci et al. (2017), computed using the . Gas masses were computed using the gas mass profiles derived from the combined density profiles. We used the stellar masses published in Chiu et al. (2016); stellar mass for SPT–CL J2146–4633 is not available. The yellow shaded area represents the baryon fraction published in Planck Collaboration XIII (2016). 
The baryon fractions for PLCK G266.6+27.3, SPT–CL J2341–5119 and SPT–CL J0546–5345 are in exce llent agreement with the previous results published in B17. SPT–CL J2106–5844 presents a larger deviation, but the hydrostatic mass computation for this object is affected by the lack of radial coverage. The use of hydrostatic mass measurements here confirms and consolidates what we found in B17: in this redshift regime, the baryon fraction does not show any dependence with respect to the mass. The density enclosed within a certain radius is higher hence more energy is required to expel the gas. We also confirm the good agreement between the baryon fraction of our sample with the fraction derived by the Planck collaboration (Planck Collaboration XIII 2016).
5. Comparison with simulations
5.1. Mass profiles
We now turn to a comparison with cosmological numerical simulations. We use the same simulated sample of five z = 1 galaxy clusters in the [4–6] × 10^{14} × M_{⊙} mass range described in Sect. 6 of B17, selected from the AGN 8.0 model of the suite of hydrodynamical cosmological simulations cosmoOWLS (Le Brun et al. 2014). These simulations include baryonic physics, and represent an extension to larger volumes of the OverWhelmingly Large Simulations project (Schaye et al. 2010).
From the simulated datasets we extracted and fitted the pressure profiles using a generalised NFW model. We then derived the simulated mass profiles by applying the hydrostatic assumption to the gNFW pressure profile in combination with the density profile (see e.g. Pratt et al. 2016). Figure 9 shows the comparison between the observed FNPL and simulated mass profiles, scaled by and , respectively, where is defined as the sum of all the particles within .
Fig. 9 Scaled hydrostatic mass profiles derived in this work and from the suite of cosmological simulations published in Le Brun et al. (2014), shown with coloured and grey solid lines. The black error bars represent the 68% dispersion of the simulated profiles computed at 0.1 and 0.5 . Simulated and Xray mass profiles were scaled by their and , respectively. Our sample and simulated cluster radial profiles were scaled by their and , respectively. The common scaled radius is indicated with R_{500}. 
The agreement over the full radial range is remarkably good. The shape, normalisation, and scatter of the simulated profiles seem to reproduce well the observations, four of the five observed profiles lie within the 68% dispersion of the theoretical profiles computed at 0.1 and 0.5 . Interestingly, there is also excellent agreement of the profiles at R_{500}, hinting that represents a robust estimate of the true mass in this mass and redshift regime. Unfortunately, we were not able to investigate the behaviour of the profiles in the core regions below 0.1R_{500}. Furthermore, as the five simulated clusters discussed here are the only objects in the cosmoOWLS cosmological box that fulfil the mass and redshift criteria, the qualitative agreement might be coincidental. A larger number of higher resolution simulations and better sampling of the Xray profiles are needed in order to make progress on this front.
5.2. Concentration
The NFW concentration is known to evolve with redshift and mass (e.g. Dutton & Macciò 2014), although at the highest masses there is surprisingly little evolution (Le Brun et al. 2018). While the mass dependence in the local Universe has been confirmed in a number of works (e.g. Pratt & Arnaud 2005; Pointecouteau et al. 2005; Vikhlinin et al. 2006; Voigt & Fabian 2006; Gastaldello et al. 2007; Buote et al. 2007; Ettori et al. 2010), evolution has received less attention (Sereno & Covone 2013; Schmidt & Allen 2007; Amodeo et al. 2016). The constraints are especially poor in the highz regime; the typical uncertainties on concentration parameters of the five clusters at z > 0.9 studied by Amodeo et al. (2016) are of the order of ±[60–80]% (their Table 2).
The very precise measurements afforded by the present observations allow us to further investigate the c–M relation and its evolution. Figure 10 shows the concentrations for four of the clusters in our sample compared to the theoretical predictions derived from the simulations of Dutton & Macciò (2014). The predictions were computed for a set of clusters in the local and distant universe, at z = 0 and z = 1, respectively, at Δ = 200. We translated their results at Δ = 500 using the NFW profile. At z = 1, we considered the concentrations plotted in their Fig. 10 rather than their power law fit, the c–M relation flattening to c_{200} = 4. (c_{500} = 2.6) in the present highmass range. Bhattacharya et al. (2013) found that the dispersion for the c–M relation is ∼30%. We used this result to roughly estimate the typical dispersion for the c–M at z = 1.
Fig. 10 Concentrationmass relation. The c_{500} is derived from the DF NFW model. Blue and red solid lines represent the theoretical relations from the suite of cosmological simulations of Dutton & Macciò (2014). Dotted lines represent the 30% scatter (Bhattacharya et al. 2013) for the z = 1 relation. The concentration of SPT–CL J2106–5844 is not reported because the data radial coverage does not allow a robust determination of c_{500}. 
Two clusters are within the 1σ dispersion of the mean expected relation, while two are [1.5–2]σ away. We iteratively computed the mean concentration, taking into account statistical errors and intrinsic scatter. Our results agrees, at the 1σ level, with the expectations: the mean concentration is ⟨c_{500}⟩ = 2.06 ± 0.67, with an estimated intrinsic dispersion of 1.2 ± 0.6. Although the sample size is small, this is the first test of the c–M relation at these redshifts with precise individual concentration measurements (i.e. errors smaller than the expected scatter).
6. Discussion and conclusions
We have presented the individual hydrostatic mass profiles of the five most distant (z ∼ 1) and massive galaxy clusters from the SPT and Planck cluster catalogues, measuring for the first time the profiles up to R_{500}. The combination of Chandra and XMMNewton, following the technique developed in B17, allowed us to overcome cosmological dimming and to derive robust measurements from the core regions out to R_{500}. The temperature profiles cover a typical radial range of [0.08–1] R_{500}, while the combined XMMNewton/Chandra density profiles are typically in the range [0.01–1.7] R_{500}. We considered both parametric (forward and backward) and nonparametric approaches to measuring the mass profiles. The main results regarding the robustness of the Xray profiles are the following:
Xray hydrostatic mass measurements at this redshift regime are remarkably robust and methodindependent. All the profiles are consistent within the uncertainties as long as they are determined in the radial range where there are density and temperature measurements. This robustness is also reflected in the determination of mass at fixed radius or at a particular density contrast.
In the very core region R < 0.08 R_{500}, where only density information is available, parametric models are necessary. The density information brings a certain constraint to the shape of the mass profile, but with an uncertainty that increases with decreasing radius.
At R_{500}, it is essential to have a temperature measurement to anchor the total mass at this radius and constrain the shape of the mass profile. Robust M_{500} estimates are only possible when this condition is fulfilled. In the absence of this constraint, model extrapolation can be rapidly divergent and can yield unphysical results.
Generally, when the radial sampling is poor (the profiles have few points) or when the profile is irregular, estimation of the mass outside the radial range probed by the temperature data is less robust and will depend strongly on the method used to measure the mass. On the other hand, if the shape of the profile is well reproduced using an NFW or Einastotype model, the resulting mass estimate outside the range with measured temperatures is more robust.
We compared and for four clusters of our sample with weak lensing (WL) mass measurements from HST observations, finding that the Xray and WL mass measurements are in agreement within the uncertainties. There is, however, an offset on average, in the sense that the Xray masses appear to be systematically higher by a factor of than the WL masses. This offset goes in the opposite direction to what has been found in previous works (e.g. WtG at the 2σ level), and is contrary to the expectations for a “hydrostatic bias”.
The above results confirm the power of combining XMMNewton and Chandra for measuring the mass profile distribution and for estimating the hydrostatic M_{500} up to z ∼ 1. We expect these results to be even less sensitive to systematic effects such as background estimation, contamination by background or foreground point sources, and the absolute temperature calibration than for local clusters. This is because object angular size is much smaller than the field of view, yielding a better constraint on the background, and the spectrum is redshifted to lower energies, where the effective area calibration of Xray telescopes is more robust. In parallel, the Chandra observations allow robust point source detection and density measurement very deep into the core regions. In contrast, WL measurements become increasingly challenging at these redshifts. The statistical quality of the WL mass data is much poorer than is reachable with Xrays, even with HST, and control of systematic effects (in particular the measure of the redshift distribution of background sources or the removal of contamination by cluster members) becomes more demanding. The fact that we find a positive HE bias is probably linked to these effects.
We then investigated the evolution by comparison with local data and with expectations from numerical simulations. The main results were:
The agreement between the hydrostatic masses at or at Δ = 500 and is remarkably good, suggesting that the M_{500}–Y_{X} relation is robust and that it can be extended to samples of disturbed and distant objects. It also suggests that there is no significant evolution between and with redshift. The comparison with WL masses would further suggest that there is no dramatic increase in the bias between the hydrostatic mass and the true mass. However, it is clear that better WL data are needed to settle this point.
We compared the scaled mass and total density profiles to those of the Xray selected local sample REXCESS. This comparison shows that on average there is excellent agreement with REXCESS at large radii. The clusters of our sample exhibit a larger dispersion in shape over the full radial range, and systematically trace the lower envelope of the REXCESS distribution in the core region. These results suggest either the presence of evolution, or an Xray/SZ selection effect.
We computed the sparsity for a large sample of clusters (the five highz objects plus REXCESS) studied with Xray observations. The sparsity enables efficient characterisation of the mass distribution in the cluster halo, and the comparison with REXCESS confirms the above.
We extended and strengthened the baryon fraction results found in B17. Using the hydrostatic mass measurements we confirmed our previous finding indicating that the baryon fraction at this redshift does not depend significantly on the halo mass, and agrees with the value from Planck Collaboration XIII (2016).
A comparison with the cosmoOWLS simulations (Le Brun et al. 2014) showed that there is excellent agreement between observed and simulated profiles, the latter derived by imitating an Xray approach. The scatter of our sample is also well reproduced by the simulations over the full radial range. We also studied the concentrationmass relation for the first time at high precision in this mass and redshift range, and found good agreement with the evolution predicted by Dutton & Macciò (2014).
This work represents the first full application of the method developed in B17, confirming that the combination of Chandra and XMMNewton is crucial in order to study highredshift objects, and allowing us to investigate the statistical properties of the mass profiles of cluster haloes in the highmass, highredshift z ∼ 1 range. Despite the small sample size, we were able to obtain a first insight into the statistical properties of these cluster haloes, suggesting profiles that are slightly less peaked than in local systems, in line with the expected theoretical evolution. However, a robust lowredshift SZselected anchor for the radial mass distribution is badly needed, especially taking into account the now wellknown issue of Xray versus SZ selection effects (Lovisari et al. 2017; AndradeSantos et al. 2017; Rossetti et al. 2017). Larger sample sizes are needed to better consolidate the average behaviour and its dispersion. In parallel, higher resolution numerical simulations of larger volumes (e.g. Le Brun et al. 2018) are needed to provide the theoretical counterparts to the type of objects we have studied here.
Published SPT masses are estimated “true” mass from the SZ signal significance, as detailed in Bleem et al. (2015). Masses in the Planck catalogue are derived iteratively from the Y_{SZ}–M_{500} relation calibrated using hydrostatic masses from XMMNewton. They are not corrected for hydrostatic bias and are on average 0.8 times smaller. In Fig. 1 of B17, and in this work, the SPT masses were renormalised by a factor of 0.8 to the Planck standard.
Any variation of the mean molecular weight with metallicity is negligible. The typical radial or redshift dependence of metallicity in clusters (Mantz et al. 2017) yields less than 0.5% variations on μ.
Acknowledgments
The results reported in this article are based on data obtained from the Chandra Data Archive and observations obtained with XMMNewton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. This work was supported by CNES. The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP720072013) ERC grant agreement no 340519.
References
 Amodeo, S., Ettori, S., Capasso, R., & Sereno, M. 2016, A&A, 590, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 AndradeSantos, F., Jones, C., Forman, W. R., et al. 2017, ApJ, 843, 76 [NASA ADS] [CrossRef] [Google Scholar]
 Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2007, A&A, 474, L37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Balmès, I., Rasera, Y., Corasaniti, P.S., & Alimi, J.M. 2014, MNRAS, 437, 2328 [NASA ADS] [CrossRef] [Google Scholar]
 Bartalucci, I., Arnaud, M., Pratt, G. W., et al. 2017, A&A, 598, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bhattacharya, S., Habib, S., Heitmann, K., & Vikhlinin, A. 2013, ApJ, 766, 32 [NASA ADS] [CrossRef] [Google Scholar]
 Biviano, A., van der Burg, R. F. J., Muzzin, A., et al. 2016, A&A, 594, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bleem, L. E., Stalder, B., de Haan, T., et al. 2015, ApJS, 216, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Buote, D. A., Gastaldello, F., Humphrey, P. J., et al. 2007, ApJ, 664, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Chiu, I., Mohr, J., McDonald, M., et al. 2016, MNRAS, 455, 258 [Google Scholar]
 Corasaniti, P. S., Ettori, S., Rasera, Y., et al. 2018, ApJ, 862, 40 [NASA ADS] [CrossRef] [Google Scholar]
 Croston, J. H., Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2006, A&A, 459, 1007 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Croston, J. H., Pratt, G. W., Böhringer, H., et al. 2008, A&A, 487, 431 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 da Silva, A. C., Kay, S. T., Liddle, A. R., & Thomas, P. A. 2004, MNRAS, 348, 1401 [NASA ADS] [CrossRef] [Google Scholar]
 Démoclès, J., Pratt, G. W., Pierini, D., et al. 2010, A&A, 517, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Diemer, B., & Kravtsov, A. V. 2015, ApJ, 799, 108 [NASA ADS] [CrossRef] [Google Scholar]
 Duffy, A. R., Schaye, J., Kay, S. T., & Dalla Vecchia, C. 2008, MNRAS, 390, L64 [NASA ADS] [CrossRef] [Google Scholar]
 Dutton, A. A., & Macciò, A. V. 2014, MNRAS, 441, 3359 [NASA ADS] [CrossRef] [Google Scholar]
 Einasto, J. 1965, Trudy Astrofizicheskogo Instituta AlmaAta, 5, 87 [NASA ADS] [Google Scholar]
 Ettori, S., Gastaldello, F., Leccardi, A., et al. 2010, A&A, 524, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ettori, S., Donnarumma, A., Pointecouteau, E., et al. 2013, Space Sci. Rev., 177, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Garmire, G. P., Bautz, M. W., Ford, P. G., Nousek, J. A., & Ricker, G. R. Jr., 2003, in XRay and GammaRay Telescopes and Instruments for Astronomy, eds. J. E. Truemper, & H. D. Tananbaum Proc. SPIE, 4851, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Gastaldello, F., Buote, D. A., Humphrey, P. J., et al. 2007, ApJ, 669, 158 [NASA ADS] [CrossRef] [Google Scholar]
 Groener, A. M., Goldberg, D. M., & Sereno, M. 2016, MNRAS, 455, 892 [NASA ADS] [CrossRef] [Google Scholar]
 Hasselfield, M., Hilton, M., Marriage, T. A., et al. 2013, J. Cosmol. Astropart. Phys., 7, 008 [NASA ADS] [CrossRef] [Google Scholar]
 Hilton, M., Hasselfield, M., Sifón, C., et al. 2018, ApJS, 235, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Hoekstra, H., Herbonnet, R., Muzzin, A., et al. 2015, MNRAS, 449, 685 [Google Scholar]
 Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 [NASA ADS] [CrossRef] [Google Scholar]
 Le Brun, A. M. C., McCarthy, I. G., Schaye, J., & Ponman, T. J. 2014, MNRAS, 441, 1270 [NASA ADS] [CrossRef] [Google Scholar]
 Le Brun, A. M. C., Arnaud, M., Pratt, G. W., & Teyssier, R. 2018, MNRAS, 473, L69 [NASA ADS] [CrossRef] [Google Scholar]
 Lewis, A. D., Buote, D. A., & Stocke, J. T. 2003, ApJ, 586, 135 [NASA ADS] [CrossRef] [Google Scholar]
 Lovisari, L., Forman, W. R., Jones, C., et al. 2017, ApJ, 846, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Mantz, A. B., Allen, S. W., Morris, R. G., et al. 2017, MNRAS, 472, 2877 [NASA ADS] [CrossRef] [Google Scholar]
 Mazzotta, P., Rasia, E., Moscardini, L., & Tormen, G. 2004, MNRAS, 354, 10 [NASA ADS] [CrossRef] [Google Scholar]
 McLaughlin, D. E. 1999, AJ, 117, 2398 [NASA ADS] [CrossRef] [Google Scholar]
 Merten, J., Meneghetti, M., Postman, M., et al. 2015, ApJ, 806, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Hayashi, E., Power, C., et al. 2004, MNRAS, 349, 1039 [NASA ADS] [CrossRef] [Google Scholar]
 Okabe, N., & Smith, G. P. 2016, MNRAS, 461, 3794 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XXVI. 2011, A&A, 536, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIX. 2014, A&A, 571, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVII. 2016, A&A, 594, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pointecouteau, E., Arnaud, M., & Pratt, G. W. 2005, A&A, 435, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pratt, G. W., & Arnaud, M. 2003, A&A, 408, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pratt, G. W., & Arnaud, M. 2005, A&A, 429, 791 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pratt, G. W., Böhringer, H., Croston, J. H., et al. 2007, A&A, 461, 71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pratt, G. W., Croston, J. H., Arnaud, M., & Böhringer, H. 2009, A&A, 498, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pratt, G. W., Arnaud, M., Piffaretti, R., et al. 2010, A&A, 511, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pratt, G. W., Pointecouteau, E., Arnaud, M., & van der Burg, R. F. J. 2016, A&A, 590, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rossetti, M., Gastaldello, F., Eckert, D., et al. 2017, MNRAS, 468, 1917 [NASA ADS] [CrossRef] [Google Scholar]
 Schaye, J., Dalla Vecchia, C., Booth, C. M., et al. 2010, MNRAS, 402, 1536 [Google Scholar]
 Schmidt, R. W., & Allen, S. W. 2007, MNRAS, 379, 209 [NASA ADS] [CrossRef] [Google Scholar]
 Schrabback, T., Applegate, D., Dietrich, J. P., et al. 2018, MNRAS, 474, 2635 [NASA ADS] [CrossRef] [Google Scholar]
 Sereno, M., & Covone, G. 2013, MNRAS, 434, 878 [NASA ADS] [CrossRef] [Google Scholar]
 Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vikhlinin, A. 2006, ApJ, 640, 710 [NASA ADS] [CrossRef] [Google Scholar]
 Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 691 [NASA ADS] [CrossRef] [Google Scholar]
 Voigt, L. M., & Fabian, A. C. 2006, MNRAS, 368, 518 [NASA ADS] [CrossRef] [Google Scholar]
 von der Linden, A., Allen, M. T., Applegate, D. E., et al. 2014, MNRAS, 439, 2 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Parametric models used
In this section we report the parametric models based on Vikhlinin et al. (2006) we used in this work. We fitted the combined density profiles with (A.1) where n_{01}, r_{c1}, α, β_{1}, r_{s}, ϵ, n_{02}, r_{c2}, and β_{2} were free parameters. Deprojected 3D temperature profiles were fitted with: (A.2) where T_{0}, T_{min}, r_{t}, b, c, and r_{cool} were free parameters. We fitted the temperature profile of SPT–CL J0546–5345 fixing both a_{cl} and b to 2. The temperature profile of SPT–CL J2106–5844 was fit using a thirddegree polynomial.
Appendix B: Density profile of SPTCLJ05465345
This work is focused on the extraction of mass profiles under the assumption of HE. For this reason, we masked the substructure in the southwest sector of SPT–CL J0546–5345, highlighted with the blue dotted circle in Fig. A.1 in B17, and derived the density and temperature profiles centred on the Xray peak. The details of profiles extraction are given in B17. Figure B.1 shows the deprojected density profiles of SPT–CL J0546–5345 using Chandra and XMMNewton datasets. Given the excellent agreement between the two, the profiles were simultaneously fitted using the parametric model of Vikhlinin et al. (2006). Its uncertainties were calculated using a Monte Carlo procedure.
Fig. B1 Normalised, scaled, and deprojected density profile of SPTCLJ05465345 measured by Chandra and XMMNewton with red and blue polygons, respectively. The black solid line and the dotted lines represent the simultaneous fit with the density parametric model of Vikhlinin et al. 2006 and its 1σ error, respectively. 
Appendix C: Mass profile errors
The relative errors of the hydrostatic mass profiles derived using the FNPL and the FP method shown in Fig. 2 exhibit a radial dependency. In the core the relative errors are larger than in the outskirts of the profiles. This is counterintuitive because the mass profiles were derived from the density profiles for which the relative error is negligible (∼1–2%) and from temperature profiles which were defined to have a the same signaltonoise ratio in each radial bin (see Sect. 3.4 of B17). The relative error of the mass profile M(<r) derived using Eq. (1) and neglecting the error on the density profile is: (C.1) where β is defined as the logarithmic derivative of temperature and density, namely β_{T}(R) and β_{ne}(R), respectively, with respect to the logarithmic derivative of the radius, β ≡ d log x/d log r. Equation (C.1) shows that the relative error is proportional to the sum of the term A = ΔT/T, and of the term B = Δβ_{T} /(β_{T} + β_{ne}).
The behaviour of the two terms as a function of the radius can be studied using simple models for the density and temperature. We employed the density and temperature parametric models in Eqs. (A.1) and (A.2), respectively, to generate “toy model” temperature and density profiles. The parameters are reported in Table C.1. We neglected the error on the density profiles and assumed a constant relative error on the temperature profile of 6%. These assumptions are a good approximation of a realistic case, where density profiles are well constrained and temperature profiles are tailored to have constant signaltonoise ratio. We generated three density profiles with different inner slopes and three temperature profiles with different shapes, shown in panels a and b of Fig. C.1, respectively. These profiles are representative of what is generally found in large samples of clusters. The density profile in the core strongly varies from cluster to cluster (coolcore, dynamically disturbed, relaxed, etc.). In the outskirts the behaviour is selfsimilar and all the profiles present a steep gradient. The temperature profile strongly depends on the cluster characteristics. The shape ranges from the “bell” shape of the coolcore clusters (T model X) to being almost flat (T model Z). For a gallery of individual density and temperature profiles, see e.g. Vikhlinin et al. (2006), Pratt et al. (2007), and Croston et al. (2008).
Fig. C1 Panel a: density profiles generated using the parametric model of Vikhlinin et al. (2006). The errors are set to 0. Panel b: temperature profiles generated using the parametric model of Vikhlinin et al. (2006). The relative error on each bin is fixed to 6%. For clarity, bins and corresponding errors are reported only for “T model X”. We report only the shape form of “T model Y” and “T model Z” with solid grey and magenta lines, respectively. Panel c: hydrostatic mass profiles derived from the three density profiles and from the temperature profile “T model X” shown in panel a and b, respectively, using the FP method. Panel d: β is defined as the logarithmic derivative of the density and temperature, namely β(n_{e}) and β(kT), with respect to the logarithm of the radius, β≡∂ ln x/β ln r. The errors of β(kT) are derived via a MonteCarlo procedure. Panel e: temperature and mass profile relative errors as a function of the radius. The dotted line represents the 6% relative error. The mass profiles are computed using the three density profiles and the “T model X”. Panel f: ratio between the error of β(kT) and the distance D between β(n_{e}) and β(kT), the distance being defined as D ≡ β(kT) – β(n_{e}). Panel g: mass profile relative errors as a function of the same quantity shown on the y axis of panel f. For this plot, we also included the results using all the temperature profiles. Points are colourcoded in order to clearly identify the inner core (blue points, upper right part of the plot) and the outskirts (red points, lower left part of the plot). For all the panels except panel g: points represent the quantity measured at fixed radial bins. The solid line is showed to guide the eye. 
Parameters of the toy models.
We took as reference the temperature profile “T model X”, and computed the hydrostatic mass profiles for the three density toy models. The results are shown in panel c of Fig. C.1. We observe the following:
the logarithmic slope, β, of the density and temperature profiles strongly depends on the radius. Panel d shows that β_{T} is small at all radii and smaller than β_{ne}. The difference between the two increases with radius;
the relative error of the mass profile ΔM/M, shown in panel e, reflects this different behaviour of the βs. In the outskirts, β_{ne} is much larger than β_{T} so that the Bterm in Eq. (C.1) becomes negligible compared to the Aterm. For this reason, the mass profile relative errors in the outskirts tend to coincide with the temperature relative error, i.e. the Aterm. This is not the case in the core where the difference between β_{T} and β_{ne} is less important and the Bterm is no longer negligible. For this reason, in the core the mass relative errors are larger than the Aterm only. The relative importance of the Bterm is related to the distance between the two βs;
the behaviour of the Bterm as a function of radius can be visualised studying the ratio between Δβ_{T} and the distance D between the β_{ne} and β_{T}, the distance being defined as D ≡ β_{T} – β_{ne}. Panel f in Fig. C.1 shows that this ratio decreases with radius as D increases. The mass measurements in the final radial bin slightly deviate from this behaviour because of the larger error of the β_{T} in the last bin. We derived β_{T} and its error Δβ_{T} by estimating the median and the 68% deviation within 1000 realisations of the temperature profile. For each realisation, we estimated the gradient for the nth bin determining the slope using the [n – 1,n,n + 1] bins. For the boundary bins we determined the slope selecting the first and last three radial bins, respectively. This is less constraining for the gradient so that the dispersion within all the realisations is greater and the resulting error Δβ_{T} is larger;
there is a clear correlation in the loglog space between the relative error on the mass ΔM/M and this ratio Δβ_{T/D}, as seen in panel g. Results using the other two temperature profiles are also shown and present the same behaviour.
The behaviour of the relative error on the mass profile is thus an intrinsic property of the hydrostatic equation, Eq. (1), and does not depend on the temperature or density profile shape. In particular, this effect is tightly linked to the general behaviour of the temperature and density gradients in galaxy clusters.
All Tables
All Figures
Fig. 1 3D temperature profiles of all the clusters of our sample. Radii are scaled by . Section 2.3.2 describes the temperature profile calculation. For each panel: the black points represent the nonparametriclike 3D temperature profiles measured using Chandra and XMMNewton, with round and squared points, respectively. The grey shaded area represents the bestfitting 3D parametric model (Vikhlinin et al. 2006). The blue and red areas represent the result of the backward fit (BP approach, see Sect. 2.4.1), assuming hydrostatic equilibrium and an NFW or an Einasto mass profile, respectively. The parametric models were estimated only in the radial range covered by the density profile. The shaded regions correspond to the 68% confidence level regions. 

In the text 
Fig. 2 Scaled mass profiles of all the clusters of our sample, derived assuming hydrostatic equilibrium (HE). All calculations used the combined XMMNewtonChandra density profile, and full details of the mass calculation methods are given in Sect. 2.4.1. Various methods give very consistent results in the radial range with temperature information, but may diverge at small and large radius in spite of the density information. For each panel: the black points represent the mass profiles obtained from the forward nonparametriclike method, using the HE equation and the nonparametriclike temperature profiles shown as black points in Fig. 1. The blue and red solid lines represent the fit of these forward nonparametriclike profiles using a NFW and an Einasto model, respectively. The grey area is the mass profile computed assuming hydrostatic equilibrium and using the parametric temperature profiles shown with a grey area in Fig. 1. The blue and red envelopes represent the mass profile computed using the backward method, i.e. fitting the observed temperature profile with a model derived from the HE equation and assuming a NFW and an Einasto profile, respectively, for the underlying total mass distribution. The parametric mass profiles are estimated in the wider radial range covered by the density profile. 

In the text 
Fig. 3 Comparison of the hydrostatic mass computed at fixed radius, , using the different methods, in units of . There is excellent agreement, with differences of less than 10%, when the radius is enclosed in the radial range covered by the spectroscopic data. 

In the text 
Fig. 4 Comparison between our Xray masses and the weak lensing masses published in Schrabback et al. (2018). All estimates for a given cluster are consistent within the statistical errors. However, there is a general trend of smaller lensing mass than the HE mass, contrary to expectation. Note also the higher statistical precision of the Xray masses. Top left panel: comparison between and weak lensing masses estimated at R_{500}, centred on the Xray peak. The grey area for the weak lensing represents the statistical errors. The black solid bars represent the sum in quadrature of systematic and statistical errors. The blue and red lines represent the bias, (1 − b), between the Xray hydrostatic and weak lensing mass as measured by Weighting the Giants (WtG, von der Linden et al. 2014) and by the Canadian Cluster Comparison Project (CCCP, Hoekstra et al. 2015), respectively. To better visualise the points, we crop the lower values of SPTCLJ23415119 and SPTCLJ05465345, which are of the order of ∼10^{−1} × 10^{14} M_{⊙}. Top right panel: same as the top left panel, except showing the comparison between the hydrostatic mass computed at , , and the weak lensing masses. Bottom left and right panels: same as the top panels except that weak lensing masses are computed using the SunyavezZeldovich (SZ) peak as the centre. The error–weighted mean ratio and corresponding errors are reported in each panel. 

In the text 
Fig. 5 Left panel: comparison between computed iteratively through the M_{500}−Y_{X} relation and the hydrostatic mass, . The estimates are consistent within the statistical errors. The of SPT‒CL J2106‒5844 is ∼3.5 greater than (see Sect. 2.4.2). For this reason, the point is off the scale and its is instead shown with the black arrow. The black dotted line is the 1 : 1 relation. Right panel: same as the left panel, except showing the comparison between and . Error–weighted mean ratio and corresponding errors are reported in each panel. 

In the text 
Fig. 6 Left panel: scaled total density profiles computed using the mass profiled derived from the DF Einasto model. For each cluster, the total radial range is that of the combined density profile; estimates beyond the radial range covered by temperature measurements are marked with dotted lines. The grey lines represent the scaled total density profiles derived from the REXCESS sample. Right panel: scaled mass profiles. The colour scheme is the same as in the left panel. The black error bars in both panels represent the 68% dispersion of the REXCESS profiles at 0.1 and 0.5 . 

In the text 
Fig. 7 Number of clusters as a function of their sparsity. The blue and gold shaded bins represent the sparsity distribution of the five highz clusters and the REXCESS sample, respectively. Individual objects are identified by symbols overplotted on the blue bins. The grey arrows represent the lower limit of the sparcity of SPT–CL J2146–4633 and SPT–CL J2106–5844. 

In the text 
Fig. 8 Baryon fractions computed at as a function of mass. The baryon fraction does not show any dependence with respect to the mass at this high z. Black points represent the baryon fraction computed using the hydrostatic mass profiles at . The grey points represent the baryon fractions published in Bartalucci et al. (2017), computed using the . Gas masses were computed using the gas mass profiles derived from the combined density profiles. We used the stellar masses published in Chiu et al. (2016); stellar mass for SPT–CL J2146–4633 is not available. The yellow shaded area represents the baryon fraction published in Planck Collaboration XIII (2016). 

In the text 
Fig. 9 Scaled hydrostatic mass profiles derived in this work and from the suite of cosmological simulations published in Le Brun et al. (2014), shown with coloured and grey solid lines. The black error bars represent the 68% dispersion of the simulated profiles computed at 0.1 and 0.5 . Simulated and Xray mass profiles were scaled by their and , respectively. Our sample and simulated cluster radial profiles were scaled by their and , respectively. The common scaled radius is indicated with R_{500}. 

In the text 
Fig. 10 Concentrationmass relation. The c_{500} is derived from the DF NFW model. Blue and red solid lines represent the theoretical relations from the suite of cosmological simulations of Dutton & Macciò (2014). Dotted lines represent the 30% scatter (Bhattacharya et al. 2013) for the z = 1 relation. The concentration of SPT–CL J2106–5844 is not reported because the data radial coverage does not allow a robust determination of c_{500}. 

In the text 
Fig. B1 Normalised, scaled, and deprojected density profile of SPTCLJ05465345 measured by Chandra and XMMNewton with red and blue polygons, respectively. The black solid line and the dotted lines represent the simultaneous fit with the density parametric model of Vikhlinin et al. 2006 and its 1σ error, respectively. 

In the text 
Fig. C1 Panel a: density profiles generated using the parametric model of Vikhlinin et al. (2006). The errors are set to 0. Panel b: temperature profiles generated using the parametric model of Vikhlinin et al. (2006). The relative error on each bin is fixed to 6%. For clarity, bins and corresponding errors are reported only for “T model X”. We report only the shape form of “T model Y” and “T model Z” with solid grey and magenta lines, respectively. Panel c: hydrostatic mass profiles derived from the three density profiles and from the temperature profile “T model X” shown in panel a and b, respectively, using the FP method. Panel d: β is defined as the logarithmic derivative of the density and temperature, namely β(n_{e}) and β(kT), with respect to the logarithm of the radius, β≡∂ ln x/β ln r. The errors of β(kT) are derived via a MonteCarlo procedure. Panel e: temperature and mass profile relative errors as a function of the radius. The dotted line represents the 6% relative error. The mass profiles are computed using the three density profiles and the “T model X”. Panel f: ratio between the error of β(kT) and the distance D between β(n_{e}) and β(kT), the distance being defined as D ≡ β(kT) – β(n_{e}). Panel g: mass profile relative errors as a function of the same quantity shown on the y axis of panel f. For this plot, we also included the results using all the temperature profiles. Points are colourcoded in order to clearly identify the inner core (blue points, upper right part of the plot) and the outskirts (red points, lower left part of the plot). For all the panels except panel g: points represent the quantity measured at fixed radial bins. The solid line is showed to guide the eye. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.