Issue 
A&A
Volume 604, August 2017



Article Number  A100  
Number of page(s)  14  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201630209  
Published online  18 August 2017 
On the evolution of the entropy and pressure profiles in Xray luminous galaxy clusters at z > 0.4
^{1} Dipartimento di Fisica e Astronomia Università di Bologna, via Piero Gobetti, 93/2, 40129 Bologna, Italy
email: vittorio.ghirardini2@unibo.it
^{2} INAF, Osservatorio Astronomico di Bologna, via Piero Gobetti, 93/3, 40129 Bologna, Italy
^{3} INFN, Sezione di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy
^{4} LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06, 75014 Paris, France
^{5} Department of Physics, LudwigMaximiliansUniversitaet, Scheinerstr. 1, 81679 Muenchen, Germany
^{6} Excellence Cluster Universe, Boltzmannstr. 2, 85748 Garching, Germany
Received: 7 December 2016
Accepted: 4 April 2017
Context. Galaxy clusters are the most recent products of hierarchical accretion over cosmological scales. The gas accreted from the cosmic field is thermalized inside the cluster halo. Gas entropy and pressure are expected to have a selfsimilar behaviour with their radial distribution following a power law and a generalized NavarroFrenkWhite profile, respectively. This has also been shown in many different hydrodynamical simulations.
Aims. We derive the spatially resolved thermodynamical properties of 47 Xray galaxy clusters observed with Chandra in the redshift range 0.4 <z< 1.2, which is one of the largest samples investigated so far with Xray spectroscopy and masses reconstructed via the hydrostatic equilibrium equation, with particular care taken to reconstruct the gas entropy and pressure radial profiles. We search for deviation from the selfsimilar behaviour and look for possible evolution with redshift.
Methods. Under the assumption of a spherically symmetric distribution of the intracluster plasma, we combine the deprojected gas density and deprojected spectral temperature profiles via the hydrostatic equilibrium equation to constrain the concentration and scale radius, which are the parameters that describe a NavarroFrenkWhite profile for each of the clusters in our sample. The temperature profile, which combined with the observed gas density profile reproduces the bestfit mass model, is then used to reconstruct the profiles of entropy and pressure. These profiles cover a median radial interval of [0.04 R_{500}–0.76 R_{500}]. After interpolating on the same radial grid and partially extrapolating up to R_{500}, these profiles are then stacked in three independent redshift bins to increase the precision of the analysis. The gas mass fraction is then used to improve the selfsimilar behaviour of the profiles by reducing the scatter among the profiles by a factor 3.
Results. The entropy and pressure profiles lie very close to the baseline prediction from gravitational structure formation. We show that these profiles deviate from the baseline prediction as function of redshift, in particular at z> 0.75, where, in the central regions, we observe higher values of the entropy (by a factor of ~2.2) and systematically lower estimates (by a factor of ~2.5) of the pressure. The effective polytropic index, which retains information about the thermal distribution of the gas, shows a slight linear positive evolution with the redshift and concentration of the dark matter distribution. A prevalence of noncool core, disturbed systems, as we observe at higher redshifts, can explain such behaviours.
Key words: galaxies: clusters: intracluster medium / galaxies: clusters: general / Xrays: galaxies: clusters / intergalactic medium
© ESO, 2017
Fig. 1 Redshift (left) and total mass (right) distribution of the clusters in our sample. 
1. Introduction
Cosmic structures evolve hierarchically from primordial density fluctuations, growing to form larger and larger systems under the action of gravity. Clusters of galaxies are the biggest virialized structures in the Universe and aggregate as bound objects at a relatively late time (z< 3). The cosmic baryons fall into the gravitational potential of the cold dark matter (CDM) halo and under the action of the collapse and subsequent shocks, adiabatic compression, and turbulence heat up to the virial temperature of few keV and form a fully ionized Xray emitting intracluster medium (ICM; Tozzi & Norman 2001; Voit et al. 2005; Zhuravleva et al. 2014).
Modelling the ICM emission by thermal brehmsstrahlung allows observations in the Xray band to provide a direct probe of the gas (electron) density, n_{e}, and gas temperature, T. Assuming the perfect gas law and a monoatomic gas, the pressure is recovered as P = Tn_{e} and the specific entropy as .
Entropy is a fundamental quantity to track the thermal history of a cluster, since it always rises when heat is produced (see Voit et al. 2005). In the presence of nonradiative processes only, lowentropy gas would sink to the centre of the cluster while highentropy gas would expand, producing the typical selfsimilar radial distribution that follows a power law with a characteristic slope of 1.1 (e.g. Tozzi & Norman 2001). Deviations from this predicted behaviour have been observed in the central region of clusters (see Ponman et al. 1999). Simulations show that nongravitational cooling and heating processes, such as radiative cooling and subsequent AGN and supernovae feedback, break selfsimilarity in the inner region of galaxy clusters (e.g. Voit et al. 2005; McCarthy et al. 2017). To justify the observed deviations from selfsimilarity, Tozzi & Norman (2001) propose a model in which an entropy floor is present before the gas is accreted by the dark matter halo.
Some dedicated papers have studied the evolution of the entropy profiles with redshift only in the last few years. McDonald et al. (2013), McDonald et al. (2014), and more recently Bartalucci et al. (2017) have concluded that both cool core clusters (CC) and noncool core clusters (NCC) in samples selected through their SunyaevZeldovich (Sunyaev & Zeldovich 1980, hereafter SZ) signal have similar entropy profiles with no relevant changes with the cosmic time.
On the other hand, pressure is the quantity the least affected by the dynamical history and nongravitational physics (Arnaud et al. 2010). Using the model introduced by Nagai et al. (2007), Arnaud et al. (2010) show that pressure has a very regular radial profile with only small deviations about the mean (less than 30 per cent outside the core regions). This is confirmed by recent observations of SZselected clusters (e.g. McDonald et al. 2014; Bartalucci et al. 2017) showing no significant deviation from the “universal profile” within R_{500} up to redshift z = 0.6.
In this work, we use the analysis of the mass distribution of 47 galaxy clusters in the redshift range 0.405–1.235 presented in Amodeo et al. (2016) to study the radial shape of pressure and entropy at different redshifts, looking for deviations from the selfsimilar behaviour and for evolution with cosmic time.
In the present study, we assume a flat ΛCDM cosmology with matter density parameter Ω_{M} = 0.3 and an Hubble constant of H_{0} = 70 km s^{1} Mpc^{1}.
Therefore the critical density of the Universe is written as (1)where E(z) ≡ H(z) /H_{0} = [Ω_{M}(1 + z)^{3} + Ω_{Λ}]^{1 / 2}.
In our study, we also consider a rescaling dependent on the mass of the halo. To apply this, we measure the quantities of interest over the regions of the cluster defined by an overdensity Δ, which is defined as a region for which the mean mass density is Δ times the critical density of the Universe.
In the following analysis, we choose Δ = 500, considering that our profiles have a radial extent of the same order of magnitude as R_{500} (see Fig. C.3). By definition, M_{500} is then equal to .
All the quoted statistical uncertainties are at 1σ level of confidence.
All the fitting processes are, unless otherwise stated, executed using the IDL routine MPFIT (Markwardt 2009), which performs a LevenbergMarquardt leastsquares fit weighting with both the errors on x and yaxes.
The paper is organized as follows: in Sect. 2 we present the sample, in Sect. 3 we present the method we apply to reconstruct the entropy and pressure profiles. The data analysis is detailed in Sect. 4 with an exhaustive discussion of our results presented in Sect. 5. We draw our conclusions in Sect. 6.
2. Sample properties
The sample selection and Xray analysis is thoroughly explained in Amodeo et al. (2016). In this work, however, we recall some of the most important properties. Because they are massive objects, the relatively long Chandra exposure time (t_{exp}> 20 ks) considered for each cluster permits us to extract at least three independent spectra over the cluster’s emission; this allows us to carry out a complete Xray analysis as is usually performed for low redshift clusters. (see Amodeo et al. 2016, for details on the sample selection and Xray analysis). We point out that clusters that are undergoing a major merger are excluded from the analysis since they would strongly break the hydrostatic equilibrium assumption, which is essential in recovering the mass profile. A completeness analysis of the sample has been made (see Amodeo et al. 2016): the selection criteria have effectively chosen the very massive high end of the cluster halo function in the investigated redshift range.
In Fig. 1, we show the redshift and mass distributions of the objects in our sample, with nine systems at z> 0.8 and eight with an estimated M_{500} larger than 2 × 10^{15}M_{⊙}. We use the centroid shift w and the concentration c_{SB} (both resulting from the analysis made in Amodeo et al. 2016) to distinguish clusters with a cool core and clusters without it. We follow the results of Cassano et al. (2010), choosing clusters with w< 0.012 and c_{SB}> 0.2 to be the classical relaxed CC. These objects are located in the upperleft quadrant of Fig. 2. The main goal of this paper is to study the evolution of entropy and pressure with cosmic time, therefore the sample is divided into redshift bins. Nevertheless a parallel analysis of the cluster sample has been made, in which clusters are divided in CC and NCC based on their morphological properties, using the criteria found by Cassano et al. (2010) mentioned above.
Fig. 2 Morphological parameters distribution in the plane of w vs. concentration c_{SB}. 
In our sample, the gas density profiles, obtained from the geometrical deprojection of the observed surface brightness profiles, cover the median radial range of [0.04 R_{500}–0.76 R_{500}] with a mean relative error of 21%. In Fig. C.3, we show the observed distribution of the minimum and maximum radius in the gas density profiles. We point out that half of the clusters have a radial coverage that extends above 0.77 R_{500}.
3. Method to reconstruct K(r) and P(r)
Amodeo et al. (2016) presented the method applied to constrain the mass distribution of the galaxy clusters in our sample under the assumption that a spherically symmetric ICM is in hydrostatic equilibrium with the underlying dark matter potential. The backward method adopted (see Ettori et al. 2013) allows us to constrain the parameters of a mass model, i.e. the concentration and scale radius for a NFW model (Navarro et al. 1997), using both the gas density profile, obtained from the geometrical deprojection of the Xray surface brightness profile, and the spatially resolved spectroscopic measurements of the gas temperature. As a byproduct of the bestfit mass model, we obtain the 3D temperature profile associated with the gas density measured in each radial bin. In other words, we obtain an estimate of the ICM temperature in each volume’s shell where a gas density is measured from the geometrical deprojection of the Xray surface brightness profile, in such a way that the bestfit mass model is reproduced by inserting the temperature and density profiles into the hydrostatic equilibrium equation (Binney & Tremaine 1987). From the combination of these profiles, thermodynamical properties such as pressure and entropy are recovered.
Following Voit (2005), temperature, pressure, and entropy associated with this halo’s overdensity are, respectively, where Eq. (5) is the Voit et al. (2005) prediction for which the radial dependence of K is rescaled from Δ = 200 to 500 with the ratio as predicted from a NFW mass model with c_{200} = 4, which is typical for massive systems.
The importance of using a proper rescaling, which includes the gas mass fraction to reach selfsimilarity in the entropy profiles, is pointed out in the work of Pratt et al. (2010). The gas fraction, f_{gas}( <r) = M_{gas}( <r) /M_{tot}( <r), is here defined as the ratio between the gas mass obtained from the integration of the gas density over the cluster’s volume and hydrostatic mass and is thus reconstructed directly from our data. We show the radial profiles of the gas fraction in Fig. 3. Throughout this paper we consider the universal baryon fraction to be f_{b} = 0.15 (Planck Collaboration XIII 2016). In Table A.1 we list additional information that is not presented in Amodeo et al. (2016), such as references on the redshifts and gas mass fraction at R_{500}.
4. Data analysis
Fig. 3 Gas fraction profile for all our clusters. The black dashed line represents the universal baryon fraction from Planck (Planck Collaboration XIII 2016). 
We combine the profiles because of the poor statistics of each profile. In order to have at least one point in each radial bin, chosen logarithmically with total number of bins equal to 30^{1}, over the range [ 0.015 R_{500},R_{500} ], we first extrapolated the data using the bestfit power law plus a constant for the entropy profile and the functional form introduced by Nagai et al. (2007) for the pressure profile, (6)where x = R/R_{500} and γ, α, and β are the central slope, intermediate slope, and outer slope, respectively, defined by a scale parameter r_{s} = R_{500}/c_{500} (R ≪ r_{s}, R ~ r_{s} and R ≫ r_{s}, respectively). All the parameters of this function are left free except for the outer slope β and the inner slope γ, which is fixed to the “universal” values (5.49 and 0.308; Arnaud et al. 2010).
We calculated the value of the thermodynamical quantities in the extrapolated radial points with the functional forms indicated above. The error associated with the thermodynamic quantities in the extrapolated radial points is the sum in quadrature of the propagated bestfitting parameters and the median uncertainty estimated in the last five radial bins of the raw profiles.
Furthermore, to investigate the average behaviour of these profiles as function of redshift, we divided the dataset into three redshift bins chosen to have approximatively the same number of clusters so that the resulting profiles have an approximatively constant signaltonoise ratio: 15 with z ∈ [ 0.4,0.52 ]; 20 with z ∈ [ 0.52,0.77 ]; and 12 with z> 0.77. In each redshift bin, the profiles are stacked in logarithmic space using the inverse of the 1σ error as weights, meaning that at each radial point the weighted mean is :where x is the logarithm of the considered thermodynamic quantities (pressure or entropy).
This “stacking” procedure ensures a higher statistical significance of the measured “mean” quantities, which propagates into a lower uncertainties in constraining the bestfit parameters of the functional forms adopted.
Fig. 4 Complete sample of our data is plotted and colour coded with respect to redshift. We can see the effect of rescaling from top to bottom, from largely scattered data to coherent data. No error bars are drawn for sake of clarity. 
4.1. Fitting procedure for the entropy profile
In Fig. C.5, we show the reconstructed entropy profiles, rescaled as described below using K_{500}, of the 47 clusters in our sample prior to the application of the extrapolation and stacking procedure.
The radial behaviour of the entropy distribution is commonly described with a power law plus a constant term accounting for the combined action of cooling and heating feedbacks, which affect the central regions (see Cavagnolo et al. 2009), (7)This functional form has an underlaying physical sense when we rescale with 100 kpc, since typically deviations from nonradiative simulations are seen below this radius, where cool core clusters and noncool core clusters actually differ. The parameter K_{0} is the central entropy, which has been used in several works (i.e. Cavagnolo et al. 2009; Voit et al. 2005; McDonald et al. 2013) to discriminate between relaxed CC clusters, with K_{0} ~ 30 keV cm^{2}, and disturbed NCC, with K_{0}> 70 keV cm^{2}.
We also consider a functional form where the scaling is carried out with respect to 0.15 R_{500} to take into account the dimension of the core in systems at different mass and redshift as follows: (8)As shown in Pratt et al. (2010), we expect the scatter among entropy profiles of the clusters to be suppressed even more when the renormalization includes both the global and radial dependence on the gas mass fraction, (9)In Fig. 4 we show how our refinement in the rescaling procedure of the entropy profiles from “no rescaling” at all to the inclusion of the dependence upon the gas fraction (see Eq. (9)), reduces the scatter among the profiles, improving the agreement with the selfsimilar prediction (Voit et al. 2005).
4.2. No rescaling
In the first panel of Fig. 4 we plot all the entropy profiles without applying any rescaling, i.e. the actual physical size in kpc versus entropy in keV cm^{2}. We name these profiles “raw data”, because the physical values are reported without any rescaling. We can deduce that our clusters have very different thermodynamic histories and in fact at each radial point the profiles spans more than one order of magnitude, which excludes selfsimilarity.
We fit all the entropy profiles via Eq. (7) to look at the occupation of the parameter space in our sample. For each cluster we obtain a value for each one of the parameters K_{0}, K_{100}, and α and we plot all those values in Fig. C.1. We fitted each one of the three histograms with a lognormal distribution. We dedicated special attention to the parameter K_{0}, since the work of Cavagnolo et al. (2009) showed that central entropy may have a bimodal distribution. We obtained ~10 clusters with K_{0} = 0 and ~3 with α = 0.
In the central regions of galaxy clusters, no significant emission from gas at very low temperature is observed (e.g. Peterson & Fabian 2006), limiting the central value of the gas entropy to be positive (K_{0}> 0). For this reason these points are excluded from Fig. C.1.
The parameter K_{0} (top panel) does not have a clear bimodal distribution. The data exhibit a peak between 150 and 200 keV cm^{2} and a significative tail for low values of K_{0}. The best fit obtained with two lognormal distributions shows one peak at 130 keV cm^{2} and one at 5 keV cm^{2}, with , while the unimodal fits yield a peak at 100 keV cm^{2} with . We used the Bayesian information criterion (BIC; Schwarz 1978) to distinguish if there are statistical differences between the two fitted models: a ΔBIC between 2 and 6 indicates positive evidence against the model with higher BIC, while at values greater than 6 the evidence is strong. We obtained ΔBIC = 0.02 and therefore we cannot discriminate between a unimodal and bimodal distribution. Even though the bestfit centers are close to the values found by Cavagnolo et al. (2009).
Fig. 5 Stacked entropy profiles and fits in the 3 redshift bins. The shaded areas represents the errors on the fit. The colours of the points are representative of the redshift bin considered. The black line represents the nonradiative prediction (Voit et al. 2005). 
The distribution of the parameter K_{100} (middle panel of Fig. C.1) is well fitted by a lognormal distribution with median value corresponding to ^{2} keV cm^{2} and . The distribution of this parameter is smooth without any significative tail; this indicates that fitting considering 100 kpc as the typical core dimension generates a welldefined distribution for K_{100} resembling a lognormal distribution.
The distribution of the power law index α (bottom panel of Fig. C.1) is, as K_{100}, well fitted by a lognormal distribution with median value ^{2} and . The visible peak has a higher value than the nonradiative prediction (Voit et al. 2005), but nevertheless it is compatible within 1σ.
Fig. 6 Fits in different redshift bins showing the errors on the fits using the shaded area. The blue profile (z> 0.77) is compatible with the low redshift profiles for intermediate radii (R> 0.07 R_{500}) while at low radii there is a difference, where the high redshift entropy profile has a flat profile. The black line represents the nonradiative prediction (Voit et al. 2005). 
4.3. Rescaling using K_{500}
The profile of a thermodynamical quantity, including entropy, should have a unique shape for galaxy clusters, after adequate rescaling (Voit 2005). Nonradiative simulations (see Voit et al. 2005) have shown that we should rescale with quantities defined with respect to the critical density in order to achieve this. As it is stated above, we use an overdensity of 500 and rescale the entropy profile using K_{500}, which is defined in Eq. (4). We can observe the effect of the scaling in the third panel of Fig. 4. The profiles we get are less scattered than the raw data, even though the dispersion about the mean is still quite high, about one order of magnitude. Nevertheless we observe that above 0.15 R_{500} most of the clusters have a selfsimilar behaviour. This is because nongravitational processes are less relevant in the outskirts of galaxy clusters (Voit 2005).
As described in Sect. 4, we radially stack the data and we fit using Eq. (8). We show the regrouped data and the fitting results with their errors bars in Fig. 5 and in top part of Table 1. As we can observe from the ratio between the data and the predicted profile (bottom panel of Fig. 5), selfsimilarity is present below 0.6 R_{500} for the two low redshift bins while it is reached only between 0.15 R_{500} and 0.6 R_{500} at high redshift (z> 0.77). Moreover the high redshift stacked profile is slightly flatter than the others with a slope of 1.0 ± 0.1 and with its best fit requiring a constant term different from 0 to reproduce the data.
In the top part of Table 1 we show the best fit dividing the sample in CC and NCC. These parameters refers to Fig. B.1. As we can see both from the best fit and the plot, CC clusters are compatible with nonradiative prediction (Voit et al. 2005) with both slope and normalization. On the other hand, the NCC subset of clusters have a rather flat entropy profile.
4.4. Rescaling using the gas fraction
Fig. 7 Dimensionless entropy K/K_{500} as a function of mass M_{500} at different radii. Here the radii of 0.15 R_{500}, R_{2500}, and R_{1000} are used. The black dashed curve is the expectation from the selfsimilar model, and the magenta line is the best fit using a power law with index α on all the data points. Red, green, and blue points indicate data points from z ∈ [ 0.4 − 0.52 ], z ∈ [ 0.52 − 0.77 ] , and z> 0.77 respectively. 
The entropy distribution depends on baryon fraction with a mass (or equivalently, temperature) dependence. Consequently entropy has both a radial and global dependence on the gas fraction (Pratt et al. 2010). When corrected by this effect, data become compatible with the nonradiative prediction (Voit et al. 2005) and the dispersion drops dramatically (last panel of Fig. 4).
A practical way to quantify the deviation from the selfsimilar prediction is shown in Fig. 7, where we show the behaviour of the dimensionless entropy profiles at some specific radii (0.15 R_{500}, R_{2500}, and R_{1000})^{2}, by interpolating the surrounding data points, with respect to the mass. When the mass decreases the deviation from the selfsimilar prediction increases. By modelling this dependance with a simple power law, we obtain a slope value which becomes smaller (in modulus) at larger radii. In particular, we point out that at the highest radius considered (R_{1000}), the profile is compatible with a constant, even though the influence of just 2 or 3 points makes the bestfit slope slightly negative.
We renormalized the entropy profiles multiplying them by gas fraction profiles (K → K × (E(z)f_{gas}(R) /f_{b})^{2 / 3}). The resulting profiles are visible in the last panel of Fig. 4. The selfsimilarity of the entropy profile is now finally clear in our dataset of clusters at different redshift. At 0.15 R_{500}, for instance, the scatter is reduced by a factor ~3 when the rescaling by the gas fraction is applied.
We then stack our profiles and we show the resulting entropy radial distribution in Fig. 6. We notice that the stacked profiles are compatible with the nonradiative prediction (Voit et al. 2005) in the radial range [0.05 R_{500}–0.7 R_{500}], with the two low redshift stacked profile pushing this compatibility down to the lower limit of our analysis. Moreover at large radii, R> 0.7 R_{500}, all the stacked profiles are slightly below the prediction.
We fit the stacked entropy profiles with a power law plus a constant and we show the results in Fig. 6 and in the bottom part of Table 1. The goodness of the fit has improved with respect to the rescaling without gas fraction and the parameters we get are closer to the nonradiative prediction.
The slope values we obtain are slightly larger than 1.1, indicating profiles that are steeper than the simulated profile, which agrees with several recent works (Voit et al. 2005; McDonald et al. 2013; Morandi & Ettori 2007; Cavagnolo et al. 2009). The situation is different in the case of z> 0.77, where the value of the central entropy is significantly different from zero. This indicates a high average central entropy for clusters at high redshifts.
In the bottom part of Table 1 we show the best fit dividing the sample into CC and NCC. This parameters refers to Fig. B.2, which is the plot analogous to Fig. 6. Similar to what happens to the split in redshift, the entropy profile of CC gets steeper than the prediction. The NCC are still flatter than the prediction and need a central entropy that is different from zero to fit the data.
5. Discussion
In the recent past, the dichotomy between cool core clusters, with a steep density profile and a drop of the temperature in the center, and noncool core clusters, with a rather flat density profile and a flat temperature profile in the center, has been studied in numerous works. The particular shape of the density profile reflects in the behaviour of the entropy profile, where CC have a low entropy floor in the center while NCC have a higher entropy floor (see Figs. B.1 and B.2).
5.1. Selfsimilarity
Nonradiative simulations predict that the thermodynamic properties of clusters of galaxies should be selfsimilar once rescaled to specific physical quantities. In the previous section we have shown that by using a proper rescaling we reach selfsimilarity (see Fig. 4). Selfsimilarity is observed for all the redshift ranges. However, only a proper rescaling using the gas fraction makes the agreement with the prediction (Voit et al. 2005) within 20% above 0.05 R_{500}. This is shown in the bottom panel of Figs. 5 and 6, where the ratios between data rescaled as indicated in Eqs. (8) and (9), respectively, and the nonradiative prediction (Voit et al. 2005) is shown. We observe that rescaling by the gas fraction is only needed in the high redshift bin to recover the self similar behaviour that is reached between 0.05 and 0.7 R_{500} within 20% from the theoretical value.
The deviations present in the inner part of the profiles may be interpreted as some form of residual energy (Morandi & Ettori 2007), which may be due to some nongravitational physics processes.
Even though the agreement with simulations is remarkable, we get slightly steeper entropy profiles, i.e. 1.13 ± 0.01 for the low redshift bins, while we get a flatter profile in the high redshift bin, with a power law index that is compatible with the prediction, 1.16 ± 0.08, but with a nonzero central entropy.
The agreement with the results of Cavagnolo et al. (2009) is excellent in all three redshift bins since we find all slopes between 1.1 and 1.2. For the high redshift bins adding a constant to reproduce the inner part of the profile is required to make our results on the slope compatible with both nonradiative predictions (Voit et al. 2005) and observations (Cavagnolo et al. 2009). Otherwise a simple power law would yield a flat profile with slope 0.90 ± 0.03.
The excess with respect to the Voit et al. (2005) selfsimilar prediction is present at low radii where most of our data are above the prediction. This extra entropy is more pronounced in low mass systems, as shown in Fig. 7 and consistent with the results obtained by Pratt et al. (2010).
5.2. Angular resolution effect
In a thermalized system, low entropy gas sinks in the center while high entropy gas floats out in the outskirts, producing an entropy profile that increases monotonously with the radius. The net effect is that the larger the central bin considered for the analysis, the higher the measured value of the entropy. In the top panel of Fig. 8, we plot the value of K_{0} for each cluster versus the radius of the innermost data point and measure an evident positive correlation (Pearson’s rank correlation of 0.76, corresponding to a significance of 8.6 × 10^{8} of its deviation from zero that is associated with the case with no correlation). This correlation becomes even more significative if we consider the same points rescaled by the properties of the halos. This is shown in the bottom panel of Fig. 8, where we measure a very significant Pearson’s rank correlation of 0.83. A similar result is shown by Panagoulia et al. (2014): the smaller the innermost radial bin, the smaller the central entropy we measure.
However, as we show in Fig. 8, the correlation between the rescaled central entropy and the innermost radial bin shows a slope value of 1.78 ± 0.04, which is about 8σ away from the expected value, and therefore does not reproduce the predicted radial dependence from Voit et al. (2005; see Eq. (5)), suggesting that the flattening has a different origin form the lack of spatial resolution.
Fig. 8 Top: central entropy vs. the innermost data point. A clear postiive correlation is measured. Bottom: central entropy versus innermost data point rescaled, an even tighter correlation is present. 
5.3. Sample completeness
In Amodeo et al. (2016), an extended study on the completeness of the sample analysed here is performed. In their Fig. 13, the completeness functions are compared to the distribution of the objects in the M_{200} − z plane, showing that the applied selection criteria effectively selected the very massive high end of the cluster halo function in the investigated redshift range.
The observed trend to detect lower masses at higher redshifts is intrinsic to the halo mass function predicted for the hierarchical structure formation in the CDM dominated Universe.
Therefore, we conclude that, within the limits of our sample selection and statistics, the observed redshift and mass dependences are a reasonable representation of the behaviour of the Xray luminous cluster population in the high mass end between z = 0.4 and z = 1.2.
5.4. Bimodality
Cavagnolo et al. (2009) have shown that the distribution of the values of the cluster central entropy reflects the dichotomy between CC and NCC clusters, finding two distinct populations peaking at 15 keV cm^{2} and 150 keV cm^{2}. This bimodality has not been confirmed in later work (e.g. Santos et al. 2010; Pratt et al. 2010).
In Sect. 3.2 we have shown that using the BIC there are no statistical differences between a unimodal and bimodal distribution. We consider whether this lack of statistical evidence is because of the poor statistics of the cluster sample. We build a bootstrap analysis of our distribution by selecting 10 000 samples of our objects (also allowing for repetitions). This approach permits us to determine whether a random sample extracted from our data show an evident bimodal distribution. The results of this analysis are shown in Fig. C.4. We observe a median value of and conclude that the distribution of this BIC is compatible with no differences between a unimodal and bimodal distribution.
In Fig. C.2 the effect of redshift on the distribution of the central entropy is shown. As the number of clusters in each redshift bin is too small, we are not able to accurately prove any redshift evolution of the central entropy, however we have indications that an evolution may be present.
At high redshift we have an important peak at entropy higher than 100 keV cm^{2}, which represents the NCC systems, and only one cluster with central entropy lower than 10 keV cm^{2}.
We observe that from high redshift to low redshift the peak at high entropy becomes less prominent and a larger fraction of clusters get a smaller value of the central entropy, so that in the lower redshift bins the majority of our clusters have a central entropy below 100 keV cm^{2}.
Nevertheless there is an indication of evolution from many NCC systems at high redshift towards mostly relaxed CC clusters at low redshift. Owing to the poor statistics of the sample we are not able to prove this scenario using statistical tests.
5.5. Evolution with redshift
We observe, within the central region, an evolution with redshift of the entropy (r< 0.1 R_{500}). It suggests that the entropy profiles are flatter at high redshifts in massive objects, or, alternatively, there is a nonzero central entropy at high redshift. This resembles an evolution in the entropy profiles. However, as we show in Fig. C.2, the fraction of NCC clusters is much higher at high redshifts than at low redshifts. This would imply that at z> 0.77 the clusters in our dataset are not able to develop a cool core similar to clusters at low redshift, meaning that cool cores are less common in the past, flattening the entropy profiles and creating a nonzero value for the central entropy. Several studies (e.g. Vikhlinin et al. 2007; Santos et al. 2008, 2010; McDonald 2011) that have investigated the evolution of the cool coreness of clusters, support this scenario, in particular on the lower relative abundance of the strongest CC at high redshift. In other works (e.g. McDonald et al. 2013, 2017), the cuspiness of the gas density is shown to decrease with increasing redshift, as a consequence of a nonevolving core that is embedded in an ambient ICM, which evolves selfsimilarly.
If we define the CC clusters as those with central entropy that is lower than 100 keV cm^{2} then from Fig. C.2 we derive that in the low redshift bins 67% of clusters are CC, which is a percentage that reduces to 50% and 40% in the intermediate and high redshift bins, respectively.
This result slightly deviates, but not in a significant way, from what has been presented in previous studies (e.g. Vikhlinin et al. 2007; McDonald et al. 2013). For example, McDonald et al. (2013; see their Fig. 12) obtain that CCs represent 40% (10−40%) of the cluster population at low (high) redshift. Their sample is SZ selected and is, therefore, less biased towards CC clusters with respect to an Xray selected sample.
This result is consistent with the hierarchical scenario of the growth of structure, given that at high redshift (z ~ 1) clusters are in the middle of their formation history and cool cores could have easily been destroyed by one of the many merger events, or not even built, if they did not have enough time to relax in the center. Moreover, it is remarkable that even in our sample, which is selected in order to have specific Xray properties and thus prone to include Xray bright centrally peaked cool core objects, we observe a relatively larger fraction of NCC at higher redshift.
Fig. 9 Top: pressure profiles for all our clusters rescaled using an overdensity of 500, colour coded with respect to the redshift bin to which each curve belongs, compared with the bestfit results in Arnaud et al. (2010; black solid line). Bottom: stacked pressure profiles (in logarithmic space) compared with the best fit of Arnaud et al. (2010); the panel at the bottom shows the ratio with respect to the “universal” pressure profile. 
5.6. Pressure
Fig. 10 Left: polytropic index for the objects in our sample vs. concentration: γ increases with larger values of c_{NFW}. Black points are the observed data points; the green and blue lines are the results of Komatsu & Seljak (2001) and Ascasibar et al. (2006), respectively; the red line with pink shaded area is the best fit using classical least squares minimization with 1σ confidence region; and the magenta line with lily shaded area is the best fit using an MCMC algorithm. Right: polytropic index for all clusters as a function of redshift. The data are better fitted by a linear relation: the effective polytropic index grows with redshift. Black points indicate the observed data points and the red line with pink shaded area indicates the best fit with 1σ dispersion using MCMC algorithm. The empty circles indicate the value of adiabatic index using the mean result for the clusters in each bin; the empty squares indicate the polytropic index of the stacked profiles in each redshift bin; the brown dashed line and dashdotted magenta line indicate the best linear fit on the “mean” and “stacked” adiabatic index. The parameters describing these lines are listed in Table 3. 
For the same dataset, we study the behaviour of the electronic pressure profile P(r) = n_{e}(r) T(r).
A generalized NFW profile, as introduced by Nagai et al. (2007) (see Eq. (6)), has been widely used to study the radial rescaled pressure profile. The bestfitting parameters [ P_{0},c_{500},γ,α,β ] = [ 8.403,1.177,0.3081,1.0510,5.49 ] obtained from Arnaud et al. (2010) represent the socalled “universal pressure profile” for galaxy clusters.
Pressure is the quantity less affected by the thermal history of the cluster (Arnaud et al. 2010). Indeed McDonald et al. (2014) found no significant evolution of the pressure profile in the analysis of SPT SZselected clusters, just a mild flattening of the profile below 0.1 R_{500}. Battaglia et al. (2012), however, suggested from the analysis of cosmological hydrodynamical simulations that a significant evolution of the pressure profile should occur beyond z = 0.7, and only outside R_{500}, as consequence of the increasing nonthermal support towards the outskirts of galaxy clusters.
In the functional form shown in Eq. (6), we fix the parameters β and γ to the fiducial value found in the work of Arnaud et al. (2010). These parameters are degenerate and therefore fixing at least one of the slopes is advised for tighter parameter distributions and a better comparison (Arnaud et al. 2010).
The pressure profiles are plotted in the top panel of Fig. 9, together with the curve of the universal profile, in order to make a comparison. All the profiles show a very similar shape, surrounding the “universal” one from both sides with an apparent discrepancy only in the normalization of the profiles. At each radius the scatter is of about one order of magnitude.
We applied the same procedure described in Sect. 4 and used in the analysis of the entropy profiles. We interpolated over the same radial grid, extrapolated up to R_{500} using the bestfit model, and stacked these pressure profiles. We obtained that the stacked curves in the two low redshift bins are compatible with the “universal” pressure profile at radii below 0.1 R_{500}, while above this value the data points are slightly below. Nevertheless the best fit on these stacked data points include the Arnaud et al. (2010) result at all radii within the error bars.
On the other hand, we observe a distinctive flattening below 0.1 R_{500} for the high redshift scaled profile for the inner part of these profiles with values that are about 30% of the “universal” pressure profile. Above this radius the data points are very close to the other redshift bins, meaning slightly below the Arnaud et al. (2010) result.
Pressure is indeed a thermodynamic property that is hardly affected by the thermodynamic history of clusters. The observed flattening at high redshift at low radii is also observed in the work of McDonald et al. (2014) and can be easily explained by the minimal presence of CC clusters at high redshift.
In Table 2 we show the best fit of Eq. (6), where we have kept β and γ fixed to the “universal” pressure profile best fit. We obtain a very good fit in all the three redshift bins; only the normalization of this functional in the high redshift bin shows a distinct discrepancy with the results of Arnaud et al. (2010), greater than 5σ. In fact the bottom panel of Fig. 9 shows that the low redshift bins points are almost on top of the “universal” profile, while the high redshift points are compatible only above 0.1 R_{500} and below this threshold the discrepancy grows to be factor 3 at 0.01 R_{500}. For completeness in the same table we show the fitting performed for the subsets of CC and NCC, referring to the data in Fig. B.3. We observe that NCC clusters bestfit results resembles the highz subsample, while CC clusters resemble lowz objects.
5.7. Polytropic index
The polytropic index γ, equal to the ratio of specific heats C_{P}/C_{V} for an ideal gas, is a common proxy when evaluating the physical state and thermal distribution of the gas. It is defined as (10)with values of γ expected to be in the range [1.1, 1.2] (Komatsu & Seljak 2001; Ascasibar et al. 2006; Shi 2016), i.e. between 1, the value describing an isothermal gas, and 5:3, the value of an isoentropic gas, when the gas is well mixed and the gas entropy per atom is constant.
Studying the evolution of the polytropic index with redshift and its relation with the concentration c of the dark matter distribution can provide a more consistent picture of the processes that regulate the hierarchical structure formation. A correlation between γ and c is expected if the radial structure of the ICM and the host halo depend on the halo mass. Komatsu & Seljak (2001) require a linear relation between concentration and polytropic index by assuming in this theoretical work that the gas traces the dark matter distribution outside the core. Ascasibar et al. (2006) have shown that c and γ conspire to produce the observed scaling relations, matching the selfsimilar slope at many overdensities.
We estimate an effective polytropic index γ by fitting the pressure with a power law as a function of the gas density. As a first step, γ is calculated for each single cluster. Then, we evaluate the weighted mean in each bin. We also calculate γ for the stacked profiles.
We look for correlations between the polytropic index and dark matter concentration as recovered from the best fit with a NFW model in Amodeo et al. (2016). Ascasibar et al. (2006) and Komatsu & Seljak (2001) have shown that between concentration and polytropic index there is a linear relation with slope of 0.005 and 0.01, respectively. Using the Markov chain Monte Carlo (MCMC) code emcee (ForemanMackey et al. 2013), we obtain a slope that is much stepper than the theoretical predictions, although with a relative uncertainty of about 70% (see bestfit values in Table 3), which makes it compatible with previous results (Komatsu & Seljak 2001; Ascasibar et al. 2006) within 1σ, and even compatible with 0 at 2σ. Moreover, the intercept we get is much smaller than what has been previously calculated.
In the right panel of Fig. 10, we show the polytropic index as a function of redshift. We measure a positive evolution with redshift, with larger values of γ (by more than 2σ) at higher redshift (see Table 3).
Value of the polytropic index, dependence on concentration, and evolution with redshift.
6. Conclusions
From the sample described in Amodeo et al. (2016), which contains one of the largest collections of clusters at z> 0.8 that are homogeneously analysed in their Xray spectral properties, we have extracted the entropy and pressure profiles of 47 clusters observed with Chandra in a redshift range from 0.4 to 1.24.
We observe higher values of the gas entropy in the central region at higher redshift, which cannot be explained as an effect owing to spatial resolution. A plausible explanation of this result is the fact that at high redshift we observe a lack of cool core clusters with respect to the low redshift sample.
Moreover at intermediate radii, between 0.1 R_{500} and 0.7 R_{500}, the selfsimilarity is recovered when we use entropy dependence on both the redshift and gas fraction and the scatter between the profiles is reduced by a factor ~3. The best fit of the stacked profiles is very similar to the Voit et al. (2005) prediction from nonradiative simulations. We also show that the pressure profiles flatten at high redshift at radii below 0.1 R_{500} with lower values by about 50% than those observed at z ≲ 0.5.
Overall, these results agrees with a scenario in which galaxy clusters are the last gravitationally bound structures to form according to the hierarchical evolution. They start forming at z ≈ 3, and at z ~ 1 − 1.5 they are still in the middle of their formation. At this epoch, cool cores could either be destroyed by merger events or have not formed yet, thereby reducing their relative number at earlier epoch. Moreover, the merging processes ongoing at high redshift would imply that objects at z ~ 1 are mostly unrelaxed, with a flatter entropy profile, which produces a clear excess in the inner parts and a deficit in the outskirts. As we show in Figs. 5 and 6, high redshift clusters indeed have a rather flat stacked entropy profile, supporting the evidence that the floating and sinking of the gas entropy has not been completed yet. The thermodynamical disturbed condition of the high redshift systems is further supported by the observed flattening of thermal pressure in the inner part of the stacked pressure profile (see Fig. 9).
Moreover, we measure a slightly significant evolution of the effective polytropic index of the ICM, we measure it by estimating dlog P_{e}/ dlog n_{e}, using the dark matter concentration and redshift with an 87% significance (1.5σ) for concentration and 73% (1.1σ) for redshift, indicating that the gas possesses a slightly larger polytropic index in systems that have a more concentrated mass distribution at higher redshift. This result supports the observational evidence that at high redshift we recover more isentropic (i.e. more flat) entropy profiles.
References
 Allen, S. W., Schmidt, R. W., Ebeling, H., Fabian, A. C., & van Speybroeck, L. 2004, MNRAS, 353, 457 [NASA ADS] [CrossRef] [Google Scholar]
 Amodeo, S., Ettori, S., Capasso, R., & Sereno, M. 2016, A&A, 590, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92, 11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ascasibar, Y., Sevilla, R., Yepes, G., Müller, V., & Gottlöber, S. 2006, MNRAS, 371, 193 [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]
 Battaglia, N., Bond, J. R., Pfrommer, C., & Sievers, J. L. 2012, ApJ, 758, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Binney, J., & Tremaine, S. 1987, Galactic dynamics (Princeton, NJ; Princeton University Press) [Google Scholar]
 Böhringer, H., Voges, W., Huchra, J. P., et al. 2000, ApJS, 129, 435 [NASA ADS] [CrossRef] [Google Scholar]
 Bonamente, M., Joy, M. K., LaRoque, S. J., et al. 2006, ApJ, 647, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Borgani, S., Governato, F., Wadsley, J., et al. 2001, ApJ, 559, L71 [NASA ADS] [CrossRef] [Google Scholar]
 Bryan, G. L. 2000, ApJ, 544, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Cassano, R., Ettori, S., Giacintucci, S., et al. 2010, ApJ, 721, L82 [NASA ADS] [CrossRef] [Google Scholar]
 Cavagnolo, K. W., Donahue, M., Voit, G. M., & Sun, M. 2009, ApJS, 182, 12 [NASA ADS] [CrossRef] [Google Scholar]
 Coble, K., Bonamente, M., Carlstrom, J. E., et al. 2007, AJ, 134, 897 [NASA ADS] [CrossRef] [Google Scholar]
 De Propris, R., Stanford, S. A., Eisenhardt, P. R., Holden, B. P., & Rosati, P. 2007, AJ, 133, 2209 [NASA ADS] [CrossRef] [Google Scholar]
 Donahue, M., Voit, G. M., O’Dea, C. P., Baum, S. A., & Sparks, W. B. 2005, ApJ, 630, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Ebeling, H., Jones, L. R., Fairley, B. W., et al. 2001, ApJ, 548, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Ebeling, H., Barrett, E., Donovan, D., et al. 2007, ApJ, 661, L33 [NASA ADS] [CrossRef] [Google Scholar]
 Edge, A. C., Ebeling, H., Bremer, M., et al. 2003, MNRAS, 339, 913 [NASA ADS] [CrossRef] [Google Scholar]
 Ellis, S. C., & Jones, L. R. 2004, MNRAS, 348, 165 [NASA ADS] [CrossRef] [Google Scholar]
 Ettori, S., Donnarumma, A., Pointecouteau, E., et al. 2013, Space Sci. Rev., 177, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Fassbender, R., Böhringer, H., Santos, J. S., et al. 2011, A&A, 527, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Foley, R. J., Andersson, K., Bazin, G., et al. 2011, ApJ, 731, 86 [NASA ADS] [CrossRef] [Google Scholar]
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [CrossRef] [Google Scholar]
 Gioia, I. M., & Luppino, G. A. 1994, ApJS, 94, 583 [NASA ADS] [CrossRef] [Google Scholar]
 Gioia, I. M., Braito, V., Branchesi, M., et al. 2004, A&A, 419, 517 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gonzalez, A. H., Zaritsky, D., Dalcanton, J. J., & Nelson, A. 2001, ApJS, 137, 117 [NASA ADS] [CrossRef] [Google Scholar]
 Gralla, M. B., Sharon, K., Gladders, M. D., et al. 2011, ApJ, 737, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Henry, J. P., Gioia, I. M., Mullis, C. R., et al. 1997, AJ, 114, 1293 [NASA ADS] [CrossRef] [Google Scholar]
 Hicks, A. K., Ellingson, E., Bautz, M., et al. 2008, ApJ, 680, 1022 [NASA ADS] [CrossRef] [Google Scholar]
 High, F. W., Stalder, B., Song, J., et al. 2010, ApJ, 723, 1736 [NASA ADS] [CrossRef] [Google Scholar]
 Holden, B. P., Stanford, S. A., Squires, G. K., et al. 2002, AJ, 124, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., & Seljak, U. 2001, MNRAS, 327, 1353 [NASA ADS] [CrossRef] [Google Scholar]
 Kotov, O., & Vikhlinin, A. 2006, ApJ, 641, 752 [NASA ADS] [CrossRef] [Google Scholar]
 Mantz, A. B., Allen, S. W., Morris, R. G., et al. 2014, MNRAS, 440, 2077 [NASA ADS] [CrossRef] [Google Scholar]
 Markwardt, C. B. 2009, Astronomical Data Analysis Software and Systems XVIII, ASP Conf. Ser., 411, 251 [NASA ADS] [Google Scholar]
 McCarthy, I. G., Schaye, J., Bird, S., & Le Brun, A. M. C. 2017, MNRAS, 465, 2936 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, M. 2011, ApJ, 742, L35 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, M., Benson, B. A., Vikhlinin, A., et al. 2013, ApJ, 774, 23 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, M., Benson, B. A., Vikhlinin, A., et al. 2014, ApJ, 794, 67 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, M., Allen, S. W., Bayliss, M., et al. 2017, ApJ, 843, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Morandi, A., & Ettori, S. 2007, MNRAS, 380, 1521 [NASA ADS] [CrossRef] [Google Scholar]
 Morandi, A., & Cui, W. 2014, MNRAS, 437, 1909 [NASA ADS] [CrossRef] [Google Scholar]
 Mullis, C. R., McNamara, B. R., Quintana, H., et al. 2003, ApJ, 594, 154 [NASA ADS] [CrossRef] [Google Scholar]
 MuñozCuartas, J. C., Macciò, A. V., Gottlöber, S., & Dutton, A. A. 2011, MNRAS, 411, 584 [NASA ADS] [CrossRef] [Google Scholar]
 Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007, ApJ, 668, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Panagoulia, E. K., Fabian, A. C., & Sanders, J. S. 2014, MNRAS, 438, 2341 [NASA ADS] [CrossRef] [Google Scholar]
 Perlman, E. S., Horner, D. J., Jones, L. R., et al. 2002, ApJS, 140, 265 [NASA ADS] [CrossRef] [Google Scholar]
 Peterson, J. R., & Fabian, A. C. 2006, Phys. Rep., 427, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ponman, T. J., Cannon, D. B., & Navarro, J. F. 1999, Nature, 397, 135 [NASA ADS] [CrossRef] [Google Scholar]
 Pratt, G. W., Croston, J. H., Arnaud, M., Bohringer, 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]
 Romer, A. K., Nichol, R. C., Holden, B. P., et al. 2000, ApJS, 126, 209 [NASA ADS] [CrossRef] [Google Scholar]
 Rosati, P., Tozzi, P., Ettori, S., et al. 2004, AJ, 127, 230 [NASA ADS] [CrossRef] [Google Scholar]
 Santos, J. S., Rosati, P., Tozzi, P., et al. 2008, A&A, 483, 35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Santos, J. S., Tozzi, P., Rosati, P., Bohringer, H. 2010, A&A, 521, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shi, X. 2016, MNRAS, 461, 1804 [NASA ADS] [CrossRef] [Google Scholar]
 Schindler, S., Guzzo, L., Ebeling, H., et al. 1995, A&A, 299, L9 [NASA ADS] [Google Scholar]
 Schwarz G. 1978, Ann. Statist., 5, 461 [Google Scholar]
 Sereno, M., & Ettori, S. 2015, MNRAS, 450, 3675 [NASA ADS] [CrossRef] [Google Scholar]
 Song, J., Zenteno, A., Stalder, B., et al. 2012, ApJ, 761, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Stocke, J. T., Morris, S. L., Gioia, I. M., et al. 1991, ApJS, 76, 813 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, I. B. 1980, ARA&A, 18, 537 [NASA ADS] [CrossRef] [Google Scholar]
 Takey, A., Schwope, A., & Lamer, G. 2011, A&A, 534, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tozzi, P., & Norman, C. 2001, ApJ, 546, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Tozzi, P., Scharf, C., & Norman, C. 2000, ApJ, 542, 106 [NASA ADS] [CrossRef] [Google Scholar]
 Vanderlinde, K., Crawford, T. M., de Haan, T., et al. 2010, ApJ, 722, 1180 [NASA ADS] [CrossRef] [Google Scholar]
 Vikhlinin, A., McNamara, B. R., Forman, W., et al. 1998, ApJ, 502, 558 [NASA ADS] [CrossRef] [Google Scholar]
 Vikhlinin, A., Burenin, R., Forman, W. R., et al. 2007, in Heating versus Cooling in Galaxies and Clusters of Galaxies (Berlin Heiderlberg: SpringerVerlag), ESO Astrophys. Symp., 48 [Google Scholar]
 Voit, G. M. 2005, Rev. Mod. Phys., 77, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Voit, G. M., Kay, S. T., Bryan G. L. 2005, MNRAS, 364, 909 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, Q. D., & Walker, S. 2014, MNRAS, 439, 1796 [NASA ADS] [CrossRef] [Google Scholar]
 Wen, Z. L., Han, J. L., & Liu, F. S. 2009, ApJS, 183, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Zhuravleva, I., Churazov, E., Schekochihin, A. A., et al. 2014, Nature, 515, 85 [Google Scholar]
Appendix A: Additional information
Addition to Table 2 in Amodeo et al. (2016): information on the sizes and gas content of these clusters.
Appendix B: CC versus NCC
Fig. B.1 Staked entropy profiles for the subsets of CC and NCC based on Cassano et al. (2010) criteria. We observe that CC profile is almost on top of the Voit et al. (2005) prediction. 
Fig. B.2 Stacked entropy profiles corrected using the gas fraction for the subsets of CC and NCC. The CC entropy profile still lies very close to the prediction, while the NCC profile has gotten closer but still flatter than the prediction. 
Fig. B.3 Pressure profiles for the subsets of CC and NCC. These profiles are very similar to those in Fig. 9, however the NCC profiles within 0.1 R_{500} are about a factor of 2 higher than the high redshift one. 
Appendix C: Plots
Fig. C.1 Top: distribution of the parameter K_{0}. The red line represents the kernel density plot with a smoothing width of 0.1. The black and blue line represens the best fit obtained using one and two lognormal distributions, respectively. The shaded grey and blue regions enclose the 68.3 % probability region around the best fit due to parameter uncertainties. Middle and bottom: distribution of the parameters K_{100} and α, respectively. The blue line with shaded area represents the best fit with a single lognormal distribution and the 1σ probability regions around it. The median value of the fit is shown directly on the graph. 
Fig. C.2 Evolution of the central entropy distribution with redshift. We can clearly see an evolution with evidence of bimodality at high redshift. The red lines represents the kernel density estimation of the distributions. The clusters that are best fitted by K_{0} = 0 are added as if they have the value of 1 kev cm^{2}. 
Fig. C.3 Distribution of the rescaled innermost (top) and outermost (bottom) radial spatial bin colour coded with redshift. 
Fig. C.4 Distribution of the bootstraps results. The black solid line represents the position of the median, while the dashed lines represent the region comprehending 68% of the distribution. The 1σ region is compatible with ΔBIC = 0. 
Fig. C.5 Rescaled entropy profiles of clusters in three redshift bins. Each colour represents data from a single cluster. 
All Tables
Value of the polytropic index, dependence on concentration, and evolution with redshift.
Addition to Table 2 in Amodeo et al. (2016): information on the sizes and gas content of these clusters.
All Figures
Fig. 1 Redshift (left) and total mass (right) distribution of the clusters in our sample. 

In the text 
Fig. 2 Morphological parameters distribution in the plane of w vs. concentration c_{SB}. 

In the text 
Fig. 3 Gas fraction profile for all our clusters. The black dashed line represents the universal baryon fraction from Planck (Planck Collaboration XIII 2016). 

In the text 
Fig. 4 Complete sample of our data is plotted and colour coded with respect to redshift. We can see the effect of rescaling from top to bottom, from largely scattered data to coherent data. No error bars are drawn for sake of clarity. 

In the text 
Fig. 5 Stacked entropy profiles and fits in the 3 redshift bins. The shaded areas represents the errors on the fit. The colours of the points are representative of the redshift bin considered. The black line represents the nonradiative prediction (Voit et al. 2005). 

In the text 
Fig. 6 Fits in different redshift bins showing the errors on the fits using the shaded area. The blue profile (z> 0.77) is compatible with the low redshift profiles for intermediate radii (R> 0.07 R_{500}) while at low radii there is a difference, where the high redshift entropy profile has a flat profile. The black line represents the nonradiative prediction (Voit et al. 2005). 

In the text 
Fig. 7 Dimensionless entropy K/K_{500} as a function of mass M_{500} at different radii. Here the radii of 0.15 R_{500}, R_{2500}, and R_{1000} are used. The black dashed curve is the expectation from the selfsimilar model, and the magenta line is the best fit using a power law with index α on all the data points. Red, green, and blue points indicate data points from z ∈ [ 0.4 − 0.52 ], z ∈ [ 0.52 − 0.77 ] , and z> 0.77 respectively. 

In the text 
Fig. 8 Top: central entropy vs. the innermost data point. A clear postiive correlation is measured. Bottom: central entropy versus innermost data point rescaled, an even tighter correlation is present. 

In the text 
Fig. 9 Top: pressure profiles for all our clusters rescaled using an overdensity of 500, colour coded with respect to the redshift bin to which each curve belongs, compared with the bestfit results in Arnaud et al. (2010; black solid line). Bottom: stacked pressure profiles (in logarithmic space) compared with the best fit of Arnaud et al. (2010); the panel at the bottom shows the ratio with respect to the “universal” pressure profile. 

In the text 
Fig. 10 Left: polytropic index for the objects in our sample vs. concentration: γ increases with larger values of c_{NFW}. Black points are the observed data points; the green and blue lines are the results of Komatsu & Seljak (2001) and Ascasibar et al. (2006), respectively; the red line with pink shaded area is the best fit using classical least squares minimization with 1σ confidence region; and the magenta line with lily shaded area is the best fit using an MCMC algorithm. Right: polytropic index for all clusters as a function of redshift. The data are better fitted by a linear relation: the effective polytropic index grows with redshift. Black points indicate the observed data points and the red line with pink shaded area indicates the best fit with 1σ dispersion using MCMC algorithm. The empty circles indicate the value of adiabatic index using the mean result for the clusters in each bin; the empty squares indicate the polytropic index of the stacked profiles in each redshift bin; the brown dashed line and dashdotted magenta line indicate the best linear fit on the “mean” and “stacked” adiabatic index. The parameters describing these lines are listed in Table 3. 

In the text 
Fig. B.1 Staked entropy profiles for the subsets of CC and NCC based on Cassano et al. (2010) criteria. We observe that CC profile is almost on top of the Voit et al. (2005) prediction. 

In the text 
Fig. B.2 Stacked entropy profiles corrected using the gas fraction for the subsets of CC and NCC. The CC entropy profile still lies very close to the prediction, while the NCC profile has gotten closer but still flatter than the prediction. 

In the text 
Fig. B.3 Pressure profiles for the subsets of CC and NCC. These profiles are very similar to those in Fig. 9, however the NCC profiles within 0.1 R_{500} are about a factor of 2 higher than the high redshift one. 

In the text 
Fig. C.1 Top: distribution of the parameter K_{0}. The red line represents the kernel density plot with a smoothing width of 0.1. The black and blue line represens the best fit obtained using one and two lognormal distributions, respectively. The shaded grey and blue regions enclose the 68.3 % probability region around the best fit due to parameter uncertainties. Middle and bottom: distribution of the parameters K_{100} and α, respectively. The blue line with shaded area represents the best fit with a single lognormal distribution and the 1σ probability regions around it. The median value of the fit is shown directly on the graph. 

In the text 
Fig. C.2 Evolution of the central entropy distribution with redshift. We can clearly see an evolution with evidence of bimodality at high redshift. The red lines represents the kernel density estimation of the distributions. The clusters that are best fitted by K_{0} = 0 are added as if they have the value of 1 kev cm^{2}. 

In the text 
Fig. C.3 Distribution of the rescaled innermost (top) and outermost (bottom) radial spatial bin colour coded with redshift. 

In the text 
Fig. C.4 Distribution of the bootstraps results. The black solid line represents the position of the median, while the dashed lines represent the region comprehending 68% of the distribution. The 1σ region is compatible with ΔBIC = 0. 

In the text 
Fig. C.5 Rescaled entropy profiles of clusters in three redshift bins. Each colour represents data from a single cluster. 

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.