Abstract

We measure the evolution of the velocity dispersion–temperature (σvTX) relation up to z = 1 using a sample of 38 galaxy clusters drawn from the XMM Cluster Survey. This work improves upon previous studies by the use of a homogeneous cluster sample and in terms of the number of high-redshift clusters included. We present here new redshift and velocity dispersion measurements for 12 z > 0.5 clusters observed with the Gemini Multi Object Spectographs instruments on the Gemini telescopes. Using an orthogonal regression method, we find that the slope of the relation is steeper than that expected if clusters were self-similar, and that the evolution of the normalization is slightly negative, but not significantly different from zero (σv ∝ T0.86±0.14E(z)−0.37±0.33). We verify our results by applying our methods to cosmological hydrodynamical simulations. The lack of evolution seen in our data is consistent with simulations that include both feedback and radiative cooling.

INTRODUCTION

Clusters of galaxies are the largest coherent gravitationally bound objects in our Universe. By studying galaxy clusters, information can be gained about the formation of galaxies, and the effect of ongoing processes such as merging and AGN feedback. They can also be used as a probe of cosmology by studying the evolution of their number density with mass and redshift (e.g. Vikhlinin et al. 2009; Hasselfield et al. 2013; Reichardt et al. 2013; Planck Collaboration XXIV 2015). However, the mass of galaxy clusters is not a quantity that can be directly measured, and therefore it needs to be determined using observable mass tracers such as X-ray properties (e.g. luminosity and temperature), the Sunyaev-Zel'dovich (SZ) effect signal, and optical properties, such as richness, line-of-sight velocity dispersion of member galaxies, and shear due to gravitational lensing (e.g. Ortiz-Gil et al. 2004; Vikhlinin et al. 2006; Rozo et al. 2009; Sifón et al. 2013; Nastasi et al. 2014; von der Linden et al. 2014; Hoekstra et al. 2015). Scaling relations are power laws between galaxy cluster properties that have the potential to allow us to measure the mass of clusters using easily observable properties such as the X-ray luminosity, X-ray temperature and the velocity dispersion. These power laws can be predicted if we assume clusters are formed in the manner described by Kaiser (1986). In this model, known as the self-similar model, all galaxy clusters and groups are essentially identical objects which have been scaled up or down (Maughan et al. 2012). Strong self-similarity refers to when galaxy clusters have been scaled by mass and weak self-similarity refers to a scaling due to the changing density of the Universe with redshift (Bower 1997). This model makes some key assumptions, as described by Kravtsov & Borgani (2012) and Maughan et al. (2012). The first assumption is that we are in an Einstein-de-Sitter Universe, Ωm = 1, so the clusters form via a single gravitational collapse at the observed redshift. Secondly, gravitational energy as a result of the collapse is the only source of energy to the intracluster medium (ICM). By introducing these assumptions we greatly simplify the problem so that properties of the density field depend on only two control parameters, the slope of the power spectrum of the initial perturbations and its normalization. The strong self-similarity determines the slope and is not expected to evolve with redshift while the weak self-similarity is responsible for the evolution of the normalization since in this simplified model it is due only to a change in density with redshift (Bryan & Norman 1998). The most commonly studied scaling relation is the luminosity–temperature relation (LXT), however there is still no consensus on how it evolves with redshift and if self-similarity holds. Some studies have found that the evolution of the normalization of this relation is consistent with self-similarity (e.g. Vikhlinin et al. 2002; Lumb et al. 2004; Maughan et al. 2006), while other studies have found zero or negative evolution (e.g. Ettori et al. 2004; Branchesi et al. 2007; Clerc et al. 2012, 2014; Hilton et al. 2012). Maughan et al. (2012) also found that the evolution of the LXT relation was not self-similar, but concluded that this could plausibly be explained by selection effects. In this paper we focus on the lesser studied relationship between the velocity dispersion of member galaxies (σv) and the X-ray temperature (TX) of the ICM. Since the velocity dispersion is a measure of the kinetic energy of the galaxies in the cluster, and temperature is related to the kinetic energy of the gas, both the gas and galaxies are tracers of the gravitational potential. One would expect a self-similar relationship of the form σv ∝ T0.5, if clusters were formed purely due to the action of gravity (Quintana & Melnick 1982; Kaiser 1986; Voit 2005). However, almost all previous studies of the relation have found a steeper power-law slope than this (see Table 1). The relation is also not expected to evolve with redshift. To date this has been tested only by Wu, Fang & Xu (1998) and Nastasi et al. (2014). Even then, all but four clusters in the Wu et al. (1998) sample are at z < 0.5. Nastasi et al. (2014) made a measurement of the relation at 0.6 < z < 1.5 using a sample of 12 clusters, obtaining results consistent with previous studies at low redshift. One may expect evolution in cluster scaling relations due to the increase of star formation and AGN activity at high redshift (e.g. Silverman et al. 2005; Magnelli et al. 2009), or due to the increase in frequency of galaxy cluster mergers with increasing redshift (e.g. Cohn & White 2005; Kay et al. 2007; Mann & Ebeling 2012). Galaxy cluster mergers are among the most energetic events in the Universe, and simulations have shown that these could result in the boosting of cluster X-ray temperatures (e.g. Randall, Sarazin & Ricker 2002; Ritchie & Thomas 2002; Poole et al. 2007). Fig. 9 in Ritchie & Thomas (2002) shows how the temperature is boosted when two equal mass systems have a head on collision with varying initial distances between their centres. All of these processes add energy into the ICM, and so we might expect to see an overall increase in the average temperatures of galaxy clusters above that expected from the self-similar case at a given redshift.

Table 1.

Previous measurements of the velocity dispersion–temperature relation. Here the relation is in the form σv = 10ATB, where σv is measured in km s−1 and T is measured in keV.

PaperNumber of clustersABRedshift rangeFitting method
Edge & Stewart (1991)232.60 ± 0.080.46 ± 0.12z < 0.1Least squares
Lubin & Bahcall (1993)412.52 ± 0.070.60 ± 0.11z < 0.2χ2
Bird, Mushotzky & Metzler (1995)222.50 ± 0.090.61 ± 0.13z < 0.1Bisector
Girardi et al. (1996)372.53 ± 0.040.61 ± 0.05z < 0.2Bisector
Ponman et al. (1996)272.54 ± 0.040.55 ± 0.05z < 0.15Bisector
White, Jones & Forman (1997)352.53 ± 0.080.60 ± 0.10z < 0.2Orthogonal
Wu et al. (1998)942.47 ± 0.060.67 ± 0.09z < 0.9Orthogonal
Wu et al. (1998)1102.57 ± 0.030.49 ± 0.05z < 0.1Orthogonal
Wu et al. (1998)392.57 ± 0.080.56 ± 0.090.1 < z < 0.9Orthogonal
Wu, Xue & Fang (1999)922.49 ± 0.030.64 ± 0.02z < 0.45Orthogonal
Xue & Wu (2000)1092.53 ± 0.030.58 ± 0.05z < 0.2Orthogonal
Nastasi et al. (2014)122.47 ± 0.190.64 ± 0.340.64 < z < 1.46Bisector
PaperNumber of clustersABRedshift rangeFitting method
Edge & Stewart (1991)232.60 ± 0.080.46 ± 0.12z < 0.1Least squares
Lubin & Bahcall (1993)412.52 ± 0.070.60 ± 0.11z < 0.2χ2
Bird, Mushotzky & Metzler (1995)222.50 ± 0.090.61 ± 0.13z < 0.1Bisector
Girardi et al. (1996)372.53 ± 0.040.61 ± 0.05z < 0.2Bisector
Ponman et al. (1996)272.54 ± 0.040.55 ± 0.05z < 0.15Bisector
White, Jones & Forman (1997)352.53 ± 0.080.60 ± 0.10z < 0.2Orthogonal
Wu et al. (1998)942.47 ± 0.060.67 ± 0.09z < 0.9Orthogonal
Wu et al. (1998)1102.57 ± 0.030.49 ± 0.05z < 0.1Orthogonal
Wu et al. (1998)392.57 ± 0.080.56 ± 0.090.1 < z < 0.9Orthogonal
Wu, Xue & Fang (1999)922.49 ± 0.030.64 ± 0.02z < 0.45Orthogonal
Xue & Wu (2000)1092.53 ± 0.030.58 ± 0.05z < 0.2Orthogonal
Nastasi et al. (2014)122.47 ± 0.190.64 ± 0.340.64 < z < 1.46Bisector
Table 1.

Previous measurements of the velocity dispersion–temperature relation. Here the relation is in the form σv = 10ATB, where σv is measured in km s−1 and T is measured in keV.

PaperNumber of clustersABRedshift rangeFitting method
Edge & Stewart (1991)232.60 ± 0.080.46 ± 0.12z < 0.1Least squares
Lubin & Bahcall (1993)412.52 ± 0.070.60 ± 0.11z < 0.2χ2
Bird, Mushotzky & Metzler (1995)222.50 ± 0.090.61 ± 0.13z < 0.1Bisector
Girardi et al. (1996)372.53 ± 0.040.61 ± 0.05z < 0.2Bisector
Ponman et al. (1996)272.54 ± 0.040.55 ± 0.05z < 0.15Bisector
White, Jones & Forman (1997)352.53 ± 0.080.60 ± 0.10z < 0.2Orthogonal
Wu et al. (1998)942.47 ± 0.060.67 ± 0.09z < 0.9Orthogonal
Wu et al. (1998)1102.57 ± 0.030.49 ± 0.05z < 0.1Orthogonal
Wu et al. (1998)392.57 ± 0.080.56 ± 0.090.1 < z < 0.9Orthogonal
Wu, Xue & Fang (1999)922.49 ± 0.030.64 ± 0.02z < 0.45Orthogonal
Xue & Wu (2000)1092.53 ± 0.030.58 ± 0.05z < 0.2Orthogonal
Nastasi et al. (2014)122.47 ± 0.190.64 ± 0.340.64 < z < 1.46Bisector
PaperNumber of clustersABRedshift rangeFitting method
Edge & Stewart (1991)232.60 ± 0.080.46 ± 0.12z < 0.1Least squares
Lubin & Bahcall (1993)412.52 ± 0.070.60 ± 0.11z < 0.2χ2
Bird, Mushotzky & Metzler (1995)222.50 ± 0.090.61 ± 0.13z < 0.1Bisector
Girardi et al. (1996)372.53 ± 0.040.61 ± 0.05z < 0.2Bisector
Ponman et al. (1996)272.54 ± 0.040.55 ± 0.05z < 0.15Bisector
White, Jones & Forman (1997)352.53 ± 0.080.60 ± 0.10z < 0.2Orthogonal
Wu et al. (1998)942.47 ± 0.060.67 ± 0.09z < 0.9Orthogonal
Wu et al. (1998)1102.57 ± 0.030.49 ± 0.05z < 0.1Orthogonal
Wu et al. (1998)392.57 ± 0.080.56 ± 0.090.1 < z < 0.9Orthogonal
Wu, Xue & Fang (1999)922.49 ± 0.030.64 ± 0.02z < 0.45Orthogonal
Xue & Wu (2000)1092.53 ± 0.030.58 ± 0.05z < 0.2Orthogonal
Nastasi et al. (2014)122.47 ± 0.190.64 ± 0.340.64 < z < 1.46Bisector

In this paper, we study a sample of 38 z < 1.0 galaxy clusters drawn from the XMM Cluster Survey (XCS; Mehrtens et al. 2012). We divide the sample into two groups: a low-redshift sample (0.0 < z < 0.5), and a high-redshift sample (0.5 < z < 1.0), such that each group has an equal number of clusters in each, and then proceed to test for evolution in the σvT relation. We describe the sample and processing of the optical and X-ray data in Section 2. Section 3 discusses the method used to determine cluster membership and for measuring the velocity dispersion and describes the methods used for fitting the σvT relation, and we present our results in Section 4. We discuss our findings in Section 5 and present our conclusions in Section 6. We assume a cosmology with Ωm = 0.27, ΩΛ = 0.73 and H0 = 70 km s−1 Mpc−1 throughout.

SAMPLE AND DATA REDUCTION

The cluster sample for this work is drawn from XCS, a serendipitous X-ray cluster survey being conducted using archival XMM–Newton data. Data Release 1 (DR1) of the XCS is described in Mehrtens et al. (2012). The overall aims of the XCS project are to measure cosmological parameters through the evolution of the cluster mass function with redshift (Sahlén et al. 2009), study the evolution of galaxies in clusters (Collins et al. 2009; Hilton et al. 2009, 2010; Stott et al. 2010) and investigate the X-ray scaling relations as a way to study the evolution of the cluster gas with redshift (Hilton et al. 2012). The XCS Automated Pipeline Algorithm (XAPA) described in Lloyd-Davies et al. (2011) was used to search the XMM archive for cluster candidates. Mehrtens et al. (2012) describes confirmation of a subset of these candidates as clusters using the combination of data from the literature and optical follow-up observations. This left a final sample of 503 X-ray confirmed galaxy clusters, 255 which were previously unknown and 356 of which were new X-ray detections. Of these, 464 have redshift estimates, and 402 have temperature measurements. For XCS-DR1 the cluster-averaged X-ray temperatures (TX) were measured using an automated pipeline described in detail in Lloyd-Davies et al. (2011). In summary, this pipeline operates as follows: spectra were generated in the 0.3–7.9 keV band using photons in the XAPA source ellipse, (where the XAPA ellipse corresponds 0.08–0.56 of R500, with a median value of to 0.36 R500, where R500 is calculated using equation 2 and table 2 from Arnaud et al. 2005); an in-field background subtraction method was used; and model fitting was done inside xspec (Schafer 1991) using an absorbed mekal (Mewe & Schrijver 1986) model and Cash statistics (Cash 1979). In the fit, the hydrogen column density was fixed to the Dickey (1990) value and the metal abundance to 0.3 times the solar value. For this paper, we have updated the TX values compared to Mehrtens et al. (2012). The pipeline is very similar to that described in Lloyd-Davies et al. (2011), but using updated versions of the XMM calibration and xspec. The median X-ray count for all the clusters in our final sample was 1919 with a minimum count of 220. We note that for only one of the clusters in our sample are the X-ray counts used for the spectral analysis less than 300. This is the minimum threshold defined by Lloyd-Davies et al. (2011) for reliable, i.e. with a fractional error of <0.4, Tx measurements at Tx > 5keV (see fig. 16 from Lloyd-Davies et al. 2011). The remaining cluster was fit using 220 counts, but has a temperature of 3.5 keV (so still has an expected fractional error of 0.4) For this paper both the samples were constructed from XCS DR1, except for one of the high-redshift sample clusters (XMMXCS J113602.9−032943.2) which is a previously unreported XCS cluster detection. Fig. 1 shows the redshift and temperature distributions of the two samples. The high-redshift sample contains more high-temperature clusters than the low-redshift sample, which may be due to selection effects which result in higher luminosity and hence higher temperature clusters being chosen at higher redshift.

The redshift and temperature distributions of the low- and high-redshift cluster samples used in this work. The solid grey marks the low-redshift sample (z < 0.5) and the shading with diagonal lines marks the high-redshift sample (z > 0.5). Note that the high-redshift sample (Tmedian = 4.5 keV) contains more high-temperature clusters than the low-redshift sample (Tmedian = 3.0 keV).
Figure 1.

The redshift and temperature distributions of the low- and high-redshift cluster samples used in this work. The solid grey marks the low-redshift sample (z < 0.5) and the shading with diagonal lines marks the high-redshift sample (z > 0.5). Note that the high-redshift sample (Tmedian = 4.5 keV) contains more high-temperature clusters than the low-redshift sample (Tmedian = 3.0 keV).

Low-redshift sample

The low-redshift sample contains 19 clusters whose properties can be found in Table 3. In order to obtain this sample, we excluded all clusters from the DR1 sample that did not have temperatures or which had a redshift z > 0.5, leaving us with a sample of 320 clusters. We performed a search in NED1 for galaxies surrounding each cluster. We included only clusters which had spectroscopic redshifts resulting in our sample size being decreased from 320 to 296. Since NED collects data from many different sources the reliability of the redshifts cannot be guaranteed. Hence, where possible we use only one source of redshifts per cluster to ensure homogeneity. These redshifts are specified to the fourth decimal place but unfortunately for most an uncertainty is not included in the original sample and therefore we assumed an accuracy of 1 per cent. We excluded galaxies located at a projected radial distance >R200 (the radius within which the mean density is 200 times the critical density of the Universe at the cluster redshift) as such galaxies are unlikely to be cluster members. To ensure we did not exclude possible members, for this initial step R200 was calculated using a fiducial velocity dispersion of 2000 km s−1 following Finn et al. (2005),
(1)
Here, σv is the line-of-sight velocity dispersion (see Section 3.1.3) and z is the redshift of the cluster. Equation (1) assumes that the galaxy velocity distribution follows an isothermal sphere dark matter profile. The fiducial R200 values span the range 2–4 Mpc. Section 3.1.2 below describes how this initial cluster membership selection was refined to give the final cluster members. We then excluded all clusters which had less than 10 galaxies as this would provide us with two few members for accurate velocity dispersion calculation leaving us with a sample size of 19 clusters.

High-redshift sample

The high-redshift sample is made up of 19 clusters whose properties can be found in Table 4. Member redshifts were determined from observations using the Gemini telescopes for 12 of these clusters (see Section 2.2.1). The other seven clusters used data obtained from Nastasi et al. (2014). They drew both on new observations and on existing data. For example, the observations of three of the Nastasi et al. (2014) clusters we have used in this paper (XMMXCS J105659.5−033728.0, XMMXCS J113602.9−032943.2 and XMMXCS J182132.9+682755.0 in Table 4) were presented in Tran et al. (1999), respectively. The observations of the other four clusters we have used in this paper were presented for the first time in Nastasi et al. (2014). These four were discovered independently (to XCS) by the XMM Newton Distant Cluster Project (XDCP; Fassbender et al. 2011). Nastasi et al. (2014) also presented galaxy redshift data for another six XDCP clusters, however we have not used those in this paper because there are insufficient galaxies to derive an accurate velocity dispersion.2 For the seven clusters that relied on Nastasi et al. (2014) data, new temperatures were obtained using XCS pipelines and the velocity dispersion was recalculated using the Nastasi et al. (2014) cluster redshift together with the method described in Section 3.

Observations

Observations of 12 z > 0.5 clusters were obtained using the Gemini Multi Object Spectographs (GMOS) on both the Gemini telescopes from 2010 to 2012. The nod-and-shuffle mode (Glazebrook & Bland-Hawthorn 2001) was used to allow better sky subtraction and shorter slit lengths when compared to conventional techniques. For all observations the R400 grating and OG515 order blocking filter were used, giving wavelength coverage of 5400–9700 Å. The GMOS field-of-view samples out to R200 at the redshifts of our sample Sifón et al. (2013). A total of 30 masks were observed with a varying number of target slitlets. Each slitlet had length 3 arcsec and width 1 arcsec. Target galaxies were selected to be fainter than the brightest cluster galaxy (which was also targeted in the slit masks), on the basis of i-band pre-imaging obtained from Gemini. We also used colour or photo-z information, where available, to maximize our efficiency in targeting cluster members. For five clusters which had r, z-band photometry from the National Optical Astronomy Observatory–XMM Cluster Survey (NXS; described in Mehrtens et al. 2012), we preferentially selected galaxies with r − z colours expected for passively evolving galaxies at the cluster redshift (see Mehrtens et al. 2012, for details). For four clusters, we used photometric redshifts for galaxies from SDSS DR7 (Abazajian et al. 2009). For XMMXCS J113602.9−032943.2, we used galaxy photo-zs that were measured from our own riz photometry obtained at the William Herschel Telescope on 2011 May 5. Observations at three different central wavelengths (7500, 7550 and 7600 Å) were used to obtain coverage over the gaps between the GMOS CCDs. For all observations an 85 percentile image quality and 50 percentile sky transparency were requested. The details of the individual observations are given in Table A1.

Spectroscopic data reduction

The data were reduced in a similar manner to Hilton et al. (2010), using pyraf and the Gemini iraf3 package. We used the tools from this package to subtract bias frames; make flat-fields; apply flat-field corrections and create mosaic images. We then applied nod-and-shuffle sky subtraction using the gnsskysub task. Wavelength calibration was determined from arc frames taken between the science frames, using standard iraf tasks. All data were then combined using a median, rejecting bad pixels using a mask constructed from the nod-and-shuffle dark frames. Finally, we combined the pairs of spectra corresponding to each nod position, and extracted one-dimensional spectra using a simple boxcar algorithm.

Galaxy redshift measurements

We measured galaxy redshifts from the spectra by cross-correlation with SDSS spectral templates4 using the RVSAO/XCSAO package for iraf (Kurtz & Mink 1998). XCSAO implements the method described by Tonry & Davis (1979). The spectra were compared to six different templates over varying redshifts with the final redshift measurement being determined after visual inspection. Redshifts were assigned a quality flag according to the following scheme: Q = 3 corresponds to two or more strongly detected features; Q = 2 refers to one strongly detected or two weakly detected features; Q = 1 one weakly detected feature and Q = 0 when no features could be identified. The features used were spectral lines, with the most commonly identified being [O ii] 3727 Å, H, K, H β and the [O iii] 4959, 5007 Å lines. Only galaxies with a quality rating of Q ≥ 2 were used in this study because these have reasonably secure redshifts. Fig. 2 shows spectra of some member galaxies of the cluster XMMXCS J025006.4−310400.8 as an example. Tables of redshifts for galaxies in each cluster field as well as histograms depicting the included/excluded members and the best-fitting Gaussian can be found in Appendices B and C, respectively. Table 2 is shown, as an example, below for cluster XMMXCS J025006.4−310400.8.

The z = 0.91 cluster XMMXCS J025006.4 − 310400.8. The left-hand panel shows the Gemini i-band image overlayed with the X-ray contours in blue. The red squares represent possible galaxy cluster members. Each possible member is labelled Mx.y, where x is the mask number and y is the object ID. The right-hand panel shows the Gemini spectra (black lines) for a subset of these galaxies. The grey bands indicate regions affected by telluric absorption lines. The red line is the best-fitting SDSS template. The green dotted vertical lines show the positions of the H and K lines at the galaxy redshift.
Figure 2.

The z = 0.91 cluster XMMXCS J025006.4 − 310400.8. The left-hand panel shows the Gemini i-band image overlayed with the X-ray contours in blue. The red squares represent possible galaxy cluster members. Each possible member is labelled Mx.y, where x is the mask number and y is the object ID. The right-hand panel shows the Gemini spectra (black lines) for a subset of these galaxies. The grey bands indicate regions affected by telluric absorption lines. The red line is the best-fitting SDSS template. The green dotted vertical lines show the positions of the H and K lines at the galaxy redshift.

Table 2.

We depict the galaxy redshifts for the cluster XMMXCS J025006.4−310400.8. Column 1 gives an arbitrary ID for each galaxy, columns 2 and 3 give the right ascension and declination, respectively, and column 4 gives the redshift of the galaxy. Column 5 gives the quality flag as explained in Section 2.2.3 and column 6 shows whether or not the galaxy was included as a member for the determination of the velocity dispersion.

IDMaskRA(J2000)Dec(J2000)zQualityMember
1GS-2010B-Q-46-0602h50m22.s92−31°03′53.00.83373
8GS-2010B-Q-46-0602h50m15.s25−31°03′33.50.72633
9GS-2010B-Q-46-0602h50m13.s27−31°03′35.70.61683
10GS-2010B-Q-46-0602h50m12.s05−31°03′08.00.71463
14GS-2010B-Q-46-0602h50m06.s63−31°03′13.70.84963
15GS-2010B-Q-46-0602h50m08.s70−31°03′49.70.90523
16GS-2010B-Q-46-0602h50m03.s78−31°03′51.50.35333
17GS-2010B-Q-46-0602h50m06.s89−31°03′51.50.92173
18GS-2010B-Q-46-0602h50m05.s48−31°03′53.00.69723
19GS-2010B-Q-46-0602h50m04.s50−31°03′51.50.91493
21GS-2010B-Q-46-0602h50m02.s79−31°04′04.90.98313
22GS-2010B-Q-46-0602h50m06.s48−31°03′56.90.90693
23GS-2010B-Q-46-0602h50m09.s04−31°04′06.30.75673
25GS-2010B-Q-46-0602h50m03.s83−31°04′34.00.89883
26GS-2010B-Q-46-0602h50m04.s24−31°04′50.60.90953
28GS-2010B-Q-46-0602h50m04.s59−31°05′41.70.61973
34GS-2010B-Q-46-0602h50m04.s26−31°07′05.20.12613
1GS-2012B-Q-011-0902h50m22.s92−31°03′53.00.59242
2GS-2012B-Q-011-0902h50m18.s11−31°03′10.50.93262
4GS-2012B-Q-011-0902h50m16.s05−31°03′23.11.00772
5GS-2012B-Q-011-0902h50m14.s89−31°03′32.10.72453
6GS-2012B-Q-011-0902h50m14.s94−31°02′56.80.99272
9GS-2012B-Q-011-0902h50m08.s98−31°03′01.10.90863
10GS-2012B-Q-011-0902h50m10.s24−31°03′27.80.90563
11GS-2012B-Q-011-0902h50m07.s01−31°01′00.90.52043
13GS-2012B-Q-011-0902h50m06.s54−31°03′44.70.91263
17GS-2012B-Q-011-0902h50m03.s99−31°03′53.00.91763
18GS-2012B-Q-011-0902h50m07.s33−31°04′10.60.91063
19GS-2012B-Q-011-0902h50m02.s67−31°03′26.00.97973
20GS-2012B-Q-011-0902h50m07.s38−31°05′28.70.64942
21GS-2012B-Q-011-0902h49m58.s10−31°03′39.60.86963
22GS-2012B-Q-011-0902h50m05.s48−31°04′40.50.90262
27GS-2012B-Q-011-0902h49m58.s68−31°05′25.80.62743
28GS-2012B-Q-011-0902h49m59.s99−31°05′07.80.88162
29GS-2012B-Q-011-0902h49m58.s60−31°04′59.20.92163
30GS-2012B-Q-011-0902h50m00.s51−31°04′44.50.56543
31GS-2012B-Q-011-0902h50m04.s26−31°07′05.20.98273
IDMaskRA(J2000)Dec(J2000)zQualityMember
1GS-2010B-Q-46-0602h50m22.s92−31°03′53.00.83373
8GS-2010B-Q-46-0602h50m15.s25−31°03′33.50.72633
9GS-2010B-Q-46-0602h50m13.s27−31°03′35.70.61683
10GS-2010B-Q-46-0602h50m12.s05−31°03′08.00.71463
14GS-2010B-Q-46-0602h50m06.s63−31°03′13.70.84963
15GS-2010B-Q-46-0602h50m08.s70−31°03′49.70.90523
16GS-2010B-Q-46-0602h50m03.s78−31°03′51.50.35333
17GS-2010B-Q-46-0602h50m06.s89−31°03′51.50.92173
18GS-2010B-Q-46-0602h50m05.s48−31°03′53.00.69723
19GS-2010B-Q-46-0602h50m04.s50−31°03′51.50.91493
21GS-2010B-Q-46-0602h50m02.s79−31°04′04.90.98313
22GS-2010B-Q-46-0602h50m06.s48−31°03′56.90.90693
23GS-2010B-Q-46-0602h50m09.s04−31°04′06.30.75673
25GS-2010B-Q-46-0602h50m03.s83−31°04′34.00.89883
26GS-2010B-Q-46-0602h50m04.s24−31°04′50.60.90953
28GS-2010B-Q-46-0602h50m04.s59−31°05′41.70.61973
34GS-2010B-Q-46-0602h50m04.s26−31°07′05.20.12613
1GS-2012B-Q-011-0902h50m22.s92−31°03′53.00.59242
2GS-2012B-Q-011-0902h50m18.s11−31°03′10.50.93262
4GS-2012B-Q-011-0902h50m16.s05−31°03′23.11.00772
5GS-2012B-Q-011-0902h50m14.s89−31°03′32.10.72453
6GS-2012B-Q-011-0902h50m14.s94−31°02′56.80.99272
9GS-2012B-Q-011-0902h50m08.s98−31°03′01.10.90863
10GS-2012B-Q-011-0902h50m10.s24−31°03′27.80.90563
11GS-2012B-Q-011-0902h50m07.s01−31°01′00.90.52043
13GS-2012B-Q-011-0902h50m06.s54−31°03′44.70.91263
17GS-2012B-Q-011-0902h50m03.s99−31°03′53.00.91763
18GS-2012B-Q-011-0902h50m07.s33−31°04′10.60.91063
19GS-2012B-Q-011-0902h50m02.s67−31°03′26.00.97973
20GS-2012B-Q-011-0902h50m07.s38−31°05′28.70.64942
21GS-2012B-Q-011-0902h49m58.s10−31°03′39.60.86963
22GS-2012B-Q-011-0902h50m05.s48−31°04′40.50.90262
27GS-2012B-Q-011-0902h49m58.s68−31°05′25.80.62743
28GS-2012B-Q-011-0902h49m59.s99−31°05′07.80.88162
29GS-2012B-Q-011-0902h49m58.s60−31°04′59.20.92163
30GS-2012B-Q-011-0902h50m00.s51−31°04′44.50.56543
31GS-2012B-Q-011-0902h50m04.s26−31°07′05.20.98273
Table 2.

We depict the galaxy redshifts for the cluster XMMXCS J025006.4−310400.8. Column 1 gives an arbitrary ID for each galaxy, columns 2 and 3 give the right ascension and declination, respectively, and column 4 gives the redshift of the galaxy. Column 5 gives the quality flag as explained in Section 2.2.3 and column 6 shows whether or not the galaxy was included as a member for the determination of the velocity dispersion.

IDMaskRA(J2000)Dec(J2000)zQualityMember
1GS-2010B-Q-46-0602h50m22.s92−31°03′53.00.83373
8GS-2010B-Q-46-0602h50m15.s25−31°03′33.50.72633
9GS-2010B-Q-46-0602h50m13.s27−31°03′35.70.61683
10GS-2010B-Q-46-0602h50m12.s05−31°03′08.00.71463
14GS-2010B-Q-46-0602h50m06.s63−31°03′13.70.84963
15GS-2010B-Q-46-0602h50m08.s70−31°03′49.70.90523
16GS-2010B-Q-46-0602h50m03.s78−31°03′51.50.35333
17GS-2010B-Q-46-0602h50m06.s89−31°03′51.50.92173
18GS-2010B-Q-46-0602h50m05.s48−31°03′53.00.69723
19GS-2010B-Q-46-0602h50m04.s50−31°03′51.50.91493
21GS-2010B-Q-46-0602h50m02.s79−31°04′04.90.98313
22GS-2010B-Q-46-0602h50m06.s48−31°03′56.90.90693
23GS-2010B-Q-46-0602h50m09.s04−31°04′06.30.75673
25GS-2010B-Q-46-0602h50m03.s83−31°04′34.00.89883
26GS-2010B-Q-46-0602h50m04.s24−31°04′50.60.90953
28GS-2010B-Q-46-0602h50m04.s59−31°05′41.70.61973
34GS-2010B-Q-46-0602h50m04.s26−31°07′05.20.12613
1GS-2012B-Q-011-0902h50m22.s92−31°03′53.00.59242
2GS-2012B-Q-011-0902h50m18.s11−31°03′10.50.93262
4GS-2012B-Q-011-0902h50m16.s05−31°03′23.11.00772
5GS-2012B-Q-011-0902h50m14.s89−31°03′32.10.72453
6GS-2012B-Q-011-0902h50m14.s94−31°02′56.80.99272
9GS-2012B-Q-011-0902h50m08.s98−31°03′01.10.90863
10GS-2012B-Q-011-0902h50m10.s24−31°03′27.80.90563
11GS-2012B-Q-011-0902h50m07.s01−31°01′00.90.52043
13GS-2012B-Q-011-0902h50m06.s54−31°03′44.70.91263
17GS-2012B-Q-011-0902h50m03.s99−31°03′53.00.91763
18GS-2012B-Q-011-0902h50m07.s33−31°04′10.60.91063
19GS-2012B-Q-011-0902h50m02.s67−31°03′26.00.97973
20GS-2012B-Q-011-0902h50m07.s38−31°05′28.70.64942
21GS-2012B-Q-011-0902h49m58.s10−31°03′39.60.86963
22GS-2012B-Q-011-0902h50m05.s48−31°04′40.50.90262
27GS-2012B-Q-011-0902h49m58.s68−31°05′25.80.62743
28GS-2012B-Q-011-0902h49m59.s99−31°05′07.80.88162
29GS-2012B-Q-011-0902h49m58.s60−31°04′59.20.92163
30GS-2012B-Q-011-0902h50m00.s51−31°04′44.50.56543
31GS-2012B-Q-011-0902h50m04.s26−31°07′05.20.98273
IDMaskRA(J2000)Dec(J2000)zQualityMember
1GS-2010B-Q-46-0602h50m22.s92−31°03′53.00.83373
8GS-2010B-Q-46-0602h50m15.s25−31°03′33.50.72633
9GS-2010B-Q-46-0602h50m13.s27−31°03′35.70.61683
10GS-2010B-Q-46-0602h50m12.s05−31°03′08.00.71463
14GS-2010B-Q-46-0602h50m06.s63−31°03′13.70.84963
15GS-2010B-Q-46-0602h50m08.s70−31°03′49.70.90523
16GS-2010B-Q-46-0602h50m03.s78−31°03′51.50.35333
17GS-2010B-Q-46-0602h50m06.s89−31°03′51.50.92173
18GS-2010B-Q-46-0602h50m05.s48−31°03′53.00.69723
19GS-2010B-Q-46-0602h50m04.s50−31°03′51.50.91493
21GS-2010B-Q-46-0602h50m02.s79−31°04′04.90.98313
22GS-2010B-Q-46-0602h50m06.s48−31°03′56.90.90693
23GS-2010B-Q-46-0602h50m09.s04−31°04′06.30.75673
25GS-2010B-Q-46-0602h50m03.s83−31°04′34.00.89883
26GS-2010B-Q-46-0602h50m04.s24−31°04′50.60.90953
28GS-2010B-Q-46-0602h50m04.s59−31°05′41.70.61973
34GS-2010B-Q-46-0602h50m04.s26−31°07′05.20.12613
1GS-2012B-Q-011-0902h50m22.s92−31°03′53.00.59242
2GS-2012B-Q-011-0902h50m18.s11−31°03′10.50.93262
4GS-2012B-Q-011-0902h50m16.s05−31°03′23.11.00772
5GS-2012B-Q-011-0902h50m14.s89−31°03′32.10.72453
6GS-2012B-Q-011-0902h50m14.s94−31°02′56.80.99272
9GS-2012B-Q-011-0902h50m08.s98−31°03′01.10.90863
10GS-2012B-Q-011-0902h50m10.s24−31°03′27.80.90563
11GS-2012B-Q-011-0902h50m07.s01−31°01′00.90.52043
13GS-2012B-Q-011-0902h50m06.s54−31°03′44.70.91263
17GS-2012B-Q-011-0902h50m03.s99−31°03′53.00.91763
18GS-2012B-Q-011-0902h50m07.s33−31°04′10.60.91063
19GS-2012B-Q-011-0902h50m02.s67−31°03′26.00.97973
20GS-2012B-Q-011-0902h50m07.s38−31°05′28.70.64942
21GS-2012B-Q-011-0902h49m58.s10−31°03′39.60.86963
22GS-2012B-Q-011-0902h50m05.s48−31°04′40.50.90262
27GS-2012B-Q-011-0902h49m58.s68−31°05′25.80.62743
28GS-2012B-Q-011-0902h49m59.s99−31°05′07.80.88162
29GS-2012B-Q-011-0902h49m58.s60−31°04′59.20.92163
30GS-2012B-Q-011-0902h50m00.s51−31°04′44.50.56543
31GS-2012B-Q-011-0902h50m04.s26−31°07′05.20.98273

ANALYSIS

Membership determination and velocity dispersion measurements

In this section we describe the methodology used to determine cluster membership and calculate the velocity dispersion of each cluster.

Cluster redshifts

For all of the clusters an estimate of the redshift is known either from the literature or from previous observations and this is used as a starting point. The peculiar velocity of each of the galaxies is calculated relative to this redshift estimate using
(2)
where vi is the peculiar velocity of the ith galaxy, zi is the redshift of the ith galaxy, z¯ is the redshift of the cluster and c is the speed of light. Extreme foreground and background sources were removed by applying a 3000 km s−1 cut with respect to the cluster redshift and then the redshift was recalculated using the biweight location method described by Beers, Flynn & Gebhardt (1990). This process was iterated until the redshift converged.

Cluster membership

A fixed gapper method, similar to that of Fadda et al. (1996) and Crawford, Wirth & Bershady (2014), was applied to determine which galaxies are cluster members. The reasoning behind this method is that by studying a histogram of the redshifts of possible members there should be a clear distinction between the cluster and the fore/background galaxies. Therefore we can exclude interlopers by finding the velocity difference between adjacent galaxies and setting a fixed gap that should not be exceeded. De Propris et al. (2002) found this optimum gap to be 1000 km s−1, which avoids the merging of subclusters but also prevents the breaking up of real systems into smaller groups. Therefore all our galaxies were sorted by peculiar velocity and the difference between all adjacent pairs was calculated. Any galaxies which had a difference between adjacent galaxies of greater than 1000 km s−1 were considered interlopers and were removed. This process was iterated until the number of galaxies converged.

Velocity dispersion

We used our confirmed galaxy cluster members to calculate an initial estimate of the velocity dispersion of each cluster using the biweight scale method described in Beers et al. (1990). We then calculated R200 using equation (1), and excluded all galaxies located at projected cluster-centric radial distances outside R200. The velocity dispersion of each cluster was then recalculated. This final radial cut did not remove more than two galaxies from the final sample for each cluster. Tables 3 and 4 list the final redshifts, velocity dispersions and R200 values for the low- and high-redshift samples, respectively.

Table 3.

Low-redshift sample (0.0 < z < 0.5): column 1 gives the name of the XCS cluster, columns 2 and 3 give its J2000 right ascension and declination. Column 4 gives the redshift, with the uncertainty found using bootstrapping and column 5 gives the redshift from literature. Column 6 gives the temperature with its positive and negative 1σ uncertainty. Column 7 gives the number of confirmed members and columns 8 and 9 give the calculated velocity dispersion and R200, respectively. The references for the redshifts are as follows: 1. Mehrtens et al. (2012), 2. Cappi, Held & Marano (1998), 3. Yoon et al. (2008), 4. Mulchaey et al. (2006), 5. Vikhlinin et al. (1998), 6.Takey, Schwope & Lamer (2013), 7. Finoguenov et al. (2007), 8. Takey, Schwope & Lamer (2011), 9. Hao et al. (2010), 10. Hennawi et al. (2008), 11. Struble & Rood (1999), 12. Burenin et al. (2007), 13. Koester et al. (2007), 14. Mullis et al. (2003) and 15. Sakelliou & Merrifield (1998).

NameRADeczzlitTMembersσvR200
(J2000)(J2000)(keV)(km s−1)(Mpc)
XMMXCS J000013.9−251052.100h00m13.s9−25°10′52.10.0845 ± 0.00040.0811.80 0.20+0.4019410  ± 0801.11
XMMXCS J003430.1−431905.600h34m30.s1−43°19′05.60.3958 ± 0.00100.4013.50 0.20+0.2022920 ± 1501.96
XMMXCS J005603.0−373248.000h56m03.s0−37°32′48.00.1659 ± 0.00090.1625.20 0.20+0.3022900 ± 1402.06
XMMXCS J015315.0+010214.201h53m15.s0+01°02′14.20.0593 ± 0.00020.0631.08 0.02+0.0212240  ± 0800.55
XMMXCS J072054.3+710900.507h20m54.s3+71°09′00.50.2309 ± 0.00050.2342.90 0.40+0.5029550  ± 0601.20
XMMXCS J081918.6+705457.508h19m18.s6+70°54′57.50.2298 ± 0.00050.2353.00 0.60+0.8019410  ± 0700.83
XMMXCS J094358.2+164120.709h43m58.s2+16°41′20.70.2539 ± 0.00050.2511.50 0.20+0.4027590  ± 0901.54
XMMXCS J095957.6+251629.009h59m57.s6+25°16′29.00.0523 ± 0.00050.0861.40 0.05+0.0515510 ± 2201.79
XMMXCS J100047.4+013926.910h00m47.s4+01°39′26.90.2202 ± 0.00060.2273.30 0.20+0.2016560 ± 1401.41
XMMXCS J100141.7+022539.810h01m41.s7+02°25′39.80.1233 ± 0.00050.1281.43 0.03+0.0626590 ± 1301.05
XMMXCS J104044.4+395710.410h40m44.s4+39°57′10.40.1389 ± 0.00070.1693.54 0.03+0.0317860 ± 1502.12
XMMXCS J111515.6+531949.511h15m15.s6+53°19′49.50.4663 ± 0.00100.47105.40 0.90+1.5016910 ± 3101.75
XMMXCS J115112.0+550655.511h51m12.s0+55°06′55.50.0791 ± 0.00030.08111.66 0.04+0.0416330 ± 1001.50
XMMXCS J123144.4+413732.012h31m44.s4+41°37′32.00.1735 ± 0.00090.18122.70 0.40+0.6010480 ± 1001.26
XMMXCS J151618.6+000531.315h16m18.s6+00°05′31.30.1200 ± 0.00050.13135.40 0.10+0.1035870 ± 2202.01
XMMXCS J161132.7+541628.316h11m32.s7+54°16′28.30.3372 ± 0.00130.3384.60 0.80+1.2012790 ± 1501.69
XMMXCS J163015.6+243423.216h30m15.s6+24°34′23.20.0625 ± 0.00030.07143.50 0.40+0.6062710 ± 1302.20
XMMXCS J223939.3−054327.422h39m39.s3−05°43′27.40.2451 ± 0.00030.24142.80 0.20+0.2068560  ± 0701.32
XMMXCS J233757.0+271121.023h37m57.s0+27°11′21.00.1237 ± 0.00070.12153.40 0.40+0.6012460 ± 1101.49
NameRADeczzlitTMembersσvR200
(J2000)(J2000)(keV)(km s−1)(Mpc)
XMMXCS J000013.9−251052.100h00m13.s9−25°10′52.10.0845 ± 0.00040.0811.80 0.20+0.4019410  ± 0801.11
XMMXCS J003430.1−431905.600h34m30.s1−43°19′05.60.3958 ± 0.00100.4013.50 0.20+0.2022920 ± 1501.96
XMMXCS J005603.0−373248.000h56m03.s0−37°32′48.00.1659 ± 0.00090.1625.20 0.20+0.3022900 ± 1402.06
XMMXCS J015315.0+010214.201h53m15.s0+01°02′14.20.0593 ± 0.00020.0631.08 0.02+0.0212240  ± 0800.55
XMMXCS J072054.3+710900.507h20m54.s3+71°09′00.50.2309 ± 0.00050.2342.90 0.40+0.5029550  ± 0601.20
XMMXCS J081918.6+705457.508h19m18.s6+70°54′57.50.2298 ± 0.00050.2353.00 0.60+0.8019410  ± 0700.83
XMMXCS J094358.2+164120.709h43m58.s2+16°41′20.70.2539 ± 0.00050.2511.50 0.20+0.4027590  ± 0901.54
XMMXCS J095957.6+251629.009h59m57.s6+25°16′29.00.0523 ± 0.00050.0861.40 0.05+0.0515510 ± 2201.79
XMMXCS J100047.4+013926.910h00m47.s4+01°39′26.90.2202 ± 0.00060.2273.30 0.20+0.2016560 ± 1401.41
XMMXCS J100141.7+022539.810h01m41.s7+02°25′39.80.1233 ± 0.00050.1281.43 0.03+0.0626590 ± 1301.05
XMMXCS J104044.4+395710.410h40m44.s4+39°57′10.40.1389 ± 0.00070.1693.54 0.03+0.0317860 ± 1502.12
XMMXCS J111515.6+531949.511h15m15.s6+53°19′49.50.4663 ± 0.00100.47105.40 0.90+1.5016910 ± 3101.75
XMMXCS J115112.0+550655.511h51m12.s0+55°06′55.50.0791 ± 0.00030.08111.66 0.04+0.0416330 ± 1001.50
XMMXCS J123144.4+413732.012h31m44.s4+41°37′32.00.1735 ± 0.00090.18122.70 0.40+0.6010480 ± 1001.26
XMMXCS J151618.6+000531.315h16m18.s6+00°05′31.30.1200 ± 0.00050.13135.40 0.10+0.1035870 ± 2202.01
XMMXCS J161132.7+541628.316h11m32.s7+54°16′28.30.3372 ± 0.00130.3384.60 0.80+1.2012790 ± 1501.69
XMMXCS J163015.6+243423.216h30m15.s6+24°34′23.20.0625 ± 0.00030.07143.50 0.40+0.6062710 ± 1302.20
XMMXCS J223939.3−054327.422h39m39.s3−05°43′27.40.2451 ± 0.00030.24142.80 0.20+0.2068560  ± 0701.32
XMMXCS J233757.0+271121.023h37m57.s0+27°11′21.00.1237 ± 0.00070.12153.40 0.40+0.6012460 ± 1101.49
Table 3.

Low-redshift sample (0.0 < z < 0.5): column 1 gives the name of the XCS cluster, columns 2 and 3 give its J2000 right ascension and declination. Column 4 gives the redshift, with the uncertainty found using bootstrapping and column 5 gives the redshift from literature. Column 6 gives the temperature with its positive and negative 1σ uncertainty. Column 7 gives the number of confirmed members and columns 8 and 9 give the calculated velocity dispersion and R200, respectively. The references for the redshifts are as follows: 1. Mehrtens et al. (2012), 2. Cappi, Held & Marano (1998), 3. Yoon et al. (2008), 4. Mulchaey et al. (2006), 5. Vikhlinin et al. (1998), 6.Takey, Schwope & Lamer (2013), 7. Finoguenov et al. (2007), 8. Takey, Schwope & Lamer (2011), 9. Hao et al. (2010), 10. Hennawi et al. (2008), 11. Struble & Rood (1999), 12. Burenin et al. (2007), 13. Koester et al. (2007), 14. Mullis et al. (2003) and 15. Sakelliou & Merrifield (1998).

NameRADeczzlitTMembersσvR200
(J2000)(J2000)(keV)(km s−1)(Mpc)
XMMXCS J000013.9−251052.100h00m13.s9−25°10′52.10.0845 ± 0.00040.0811.80 0.20+0.4019410  ± 0801.11
XMMXCS J003430.1−431905.600h34m30.s1−43°19′05.60.3958 ± 0.00100.4013.50 0.20+0.2022920 ± 1501.96
XMMXCS J005603.0−373248.000h56m03.s0−37°32′48.00.1659 ± 0.00090.1625.20 0.20+0.3022900 ± 1402.06
XMMXCS J015315.0+010214.201h53m15.s0+01°02′14.20.0593 ± 0.00020.0631.08 0.02+0.0212240  ± 0800.55
XMMXCS J072054.3+710900.507h20m54.s3+71°09′00.50.2309 ± 0.00050.2342.90 0.40+0.5029550  ± 0601.20
XMMXCS J081918.6+705457.508h19m18.s6+70°54′57.50.2298 ± 0.00050.2353.00 0.60+0.8019410  ± 0700.83
XMMXCS J094358.2+164120.709h43m58.s2+16°41′20.70.2539 ± 0.00050.2511.50 0.20+0.4027590  ± 0901.54
XMMXCS J095957.6+251629.009h59m57.s6+25°16′29.00.0523 ± 0.00050.0861.40 0.05+0.0515510 ± 2201.79
XMMXCS J100047.4+013926.910h00m47.s4+01°39′26.90.2202 ± 0.00060.2273.30 0.20+0.2016560 ± 1401.41
XMMXCS J100141.7+022539.810h01m41.s7+02°25′39.80.1233 ± 0.00050.1281.43 0.03+0.0626590 ± 1301.05
XMMXCS J104044.4+395710.410h40m44.s4+39°57′10.40.1389 ± 0.00070.1693.54 0.03+0.0317860 ± 1502.12
XMMXCS J111515.6+531949.511h15m15.s6+53°19′49.50.4663 ± 0.00100.47105.40 0.90+1.5016910 ± 3101.75
XMMXCS J115112.0+550655.511h51m12.s0+55°06′55.50.0791 ± 0.00030.08111.66 0.04+0.0416330 ± 1001.50
XMMXCS J123144.4+413732.012h31m44.s4+41°37′32.00.1735 ± 0.00090.18122.70 0.40+0.6010480 ± 1001.26
XMMXCS J151618.6+000531.315h16m18.s6+00°05′31.30.1200 ± 0.00050.13135.40 0.10+0.1035870 ± 2202.01
XMMXCS J161132.7+541628.316h11m32.s7+54°16′28.30.3372 ± 0.00130.3384.60 0.80+1.2012790 ± 1501.69
XMMXCS J163015.6+243423.216h30m15.s6+24°34′23.20.0625 ± 0.00030.07143.50 0.40+0.6062710 ± 1302.20
XMMXCS J223939.3−054327.422h39m39.s3−05°43′27.40.2451 ± 0.00030.24142.80 0.20+0.2068560  ± 0701.32
XMMXCS J233757.0+271121.023h37m57.s0+27°11′21.00.1237 ± 0.00070.12153.40 0.40+0.6012460 ± 1101.49
NameRADeczzlitTMembersσvR200
(J2000)(J2000)(keV)(km s−1)(Mpc)
XMMXCS J000013.9−251052.100h00m13.s9−25°10′52.10.0845 ± 0.00040.0811.80 0.20+0.4019410  ± 0801.11
XMMXCS J003430.1−431905.600h34m30.s1−43°19′05.60.3958 ± 0.00100.4013.50 0.20+0.2022920 ± 1501.96
XMMXCS J005603.0−373248.000h56m03.s0−37°32′48.00.1659 ± 0.00090.1625.20 0.20+0.3022900 ± 1402.06
XMMXCS J015315.0+010214.201h53m15.s0+01°02′14.20.0593 ± 0.00020.0631.08 0.02+0.0212240  ± 0800.55
XMMXCS J072054.3+710900.507h20m54.s3+71°09′00.50.2309 ± 0.00050.2342.90 0.40+0.5029550  ± 0601.20
XMMXCS J081918.6+705457.508h19m18.s6+70°54′57.50.2298 ± 0.00050.2353.00 0.60+0.8019410  ± 0700.83
XMMXCS J094358.2+164120.709h43m58.s2+16°41′20.70.2539 ± 0.00050.2511.50 0.20+0.4027590  ± 0901.54
XMMXCS J095957.6+251629.009h59m57.s6+25°16′29.00.0523 ± 0.00050.0861.40 0.05+0.0515510 ± 2201.79
XMMXCS J100047.4+013926.910h00m47.s4+01°39′26.90.2202 ± 0.00060.2273.30 0.20+0.2016560 ± 1401.41
XMMXCS J100141.7+022539.810h01m41.s7+02°25′39.80.1233 ± 0.00050.1281.43 0.03+0.0626590 ± 1301.05
XMMXCS J104044.4+395710.410h40m44.s4+39°57′10.40.1389 ± 0.00070.1693.54 0.03+0.0317860 ± 1502.12
XMMXCS J111515.6+531949.511h15m15.s6+53°19′49.50.4663 ± 0.00100.47105.40 0.90+1.5016910 ± 3101.75
XMMXCS J115112.0+550655.511h51m12.s0+55°06′55.50.0791 ± 0.00030.08111.66 0.04+0.0416330 ± 1001.50
XMMXCS J123144.4+413732.012h31m44.s4+41°37′32.00.1735 ± 0.00090.18122.70 0.40+0.6010480 ± 1001.26
XMMXCS J151618.6+000531.315h16m18.s6+00°05′31.30.1200 ± 0.00050.13135.40 0.10+0.1035870 ± 2202.01
XMMXCS J161132.7+541628.316h11m32.s7+54°16′28.30.3372 ± 0.00130.3384.60 0.80+1.2012790 ± 1501.69
XMMXCS J163015.6+243423.216h30m15.s6+24°34′23.20.0625 ± 0.00030.07143.50 0.40+0.6062710 ± 1302.20
XMMXCS J223939.3−054327.422h39m39.s3−05°43′27.40.2451 ± 0.00030.24142.80 0.20+0.2068560  ± 0701.32
XMMXCS J233757.0+271121.023h37m57.s0+27°11′21.00.1237 ± 0.00070.12153.40 0.40+0.6012460 ± 1101.49
Table 4.

High-redshift sample (0.5 < z < 1.0): all columns are as explained in Table 3. The superscripts in column one indicate the origin of redshift data when it did not come from our own observations; aare all clusters taken from the XDCP survey and were observed with the VLT–FORS2 spectograph (Nastasi et al. 2014), balso known as MS1054-03 was observed with Keck for 8.6 h (Tran et al. 1999), cwas observed with Keck (Donahue et al. 1999), dalso known as RXJ1821.6+6827 was observed with CFHT, Keck and the 2.2 m telescope at the University of Hawaii (Gioia et al. 2004). The references for the redshifts are as follows: 1. Nastasi et al. (2014), 2. Scharf et al. (1997), 3. Mehrtens et al. (2012), 4. Adami et al. (2011), 5. Šuhada et al. (2011), 6. Bellagamba et al. (2011), 7. Gioia & Luppino (1994), 8. Basilakos et al. (2004), 9. Gioia et al. (2004) and 10. Perlman et al. (2002).

NameRADeczzlitTMembersσvR200
(J2000)(J2000)(keV)(km s−1)(Mpc)
aXMMXCS J000216.1−355633.800h02m16.s1−35°56′33.80.7709 ± 0.00210.7714.83 0.76+1.015131100 ± 1901.77
XMMXCS J005656.6−274031.900h56m56.s6−27°40′31.90.5601 ± 0.00070.5623.30 0.63+0.9415380 ± 600.66
XMMXCS J015241.1−133855.901h52m41.s1−13°38′55.90.8268 ± 0.00100.8233.23 0.31+0.3829840 ± 1501.33
XMMXCS J021734.7−051326.902h17m34.s7−05°13′26.90.6467 ± 0.00120.6542.23 0.44+0.9012620 ± 2101.11
XMMXCS J025006.4−310400.802h50m06.s4−31°04′00.80.9100 ± 0.00240.9024.50 0.88+1.33131120 ± 2601.66
XMMXCS J030205.1−000003.603h02m05.s1−00°00′03.60.6450 ± 0.00070.6555.82 1.32+2.0916610 ± 1801.04
aXMMXCS J095417.1+173805.909h54m17.s117°38′05.90.8272 ± 0.00170.8213.65 0.51+0.6210940 ± 3101.42
XMMXCS J095940.7+023113.409h59m40.s7+02°31′13.40.7291 ± 0.00050.7265.02 0.55+0.6825470 ± 900.88
bXMMXCS J105659.5−033728.010h56m59.s5−03°37′28.00.8336 ± 0.00130.8277.57 0.40+0.43291010 ± 1201.57
XMMXCS J112349.4+052955.111h23m49.s4+05°29′55.10.6550 ± 0.00070.6534.62 0.95+1.5517600 ± 2101.05
XMMXCS J113602.9−032943.211h36m02.s9−03°29′43.20.8297 ± 0.00113.32 0.78+1.2021700 ± 1101.06
cXMMXCS J114023.0+660819.011h40m23.s9+66°08′19.00.7855 ± 0.00150.7877.47 0.77+0.9222950 ± 1001.51
aXMMXCS J124312.2−131307.212h43m12.s2−13°13′07.20.7910 ± 0.00140.8014.92 1.54+2.9311790 ± 4601.19
XMMXCS J134305.1−000056.813h43m05.s1−00°00′56.80.6894 ± 0.00110.6784.49 0.57+0.7223920 ± 1701.72
XMMXCS J145009.3+090428.814h50m09.s3+09°04′28.80.6412 ± 0.00070.6033.84 0.55+0.6622630 ± 901.07
dXMMXCS J182132.9+682755.018h21m32.s9+68°27′55.00.8166 ± 0.00110.8294.49 0.56+0.7919860 ± 1301.34
XMMXCS J215221.0−273022.621h52m21.s0−27°30′22.60.8276 ± 0.00110.8232.18 0.45+0.6715530 ± 1500.86
XMMXCS J230247.7+084355.923h02m47.s7+08°43′55.90.7187 ± 0.00140.72105.29 0.50+0.59221010 ± 1301.60
aXMMXCS J235616.4−344144.323h56m16.s4−34°41′44.30.9391 ± 0.00120.9414.57 0.41+0.4810670 ± 2600.91
NameRADeczzlitTMembersσvR200
(J2000)(J2000)(keV)(km s−1)(Mpc)
aXMMXCS J000216.1−355633.800h02m16.s1−35°56′33.80.7709 ± 0.00210.7714.83 0.76+1.015131100 ± 1901.77
XMMXCS J005656.6−274031.900h56m56.s6−27°40′31.90.5601 ± 0.00070.5623.30 0.63+0.9415380 ± 600.66
XMMXCS J015241.1−133855.901h52m41.s1−13°38′55.90.8268 ± 0.00100.8233.23 0.31+0.3829840 ± 1501.33
XMMXCS J021734.7−051326.902h17m34.s7−05°13′26.90.6467 ± 0.00120.6542.23 0.44+0.9012620 ± 2101.11
XMMXCS J025006.4−310400.802h50m06.s4−31°04′00.80.9100 ± 0.00240.9024.50 0.88+1.33131120 ± 2601.66
XMMXCS J030205.1−000003.603h02m05.s1−00°00′03.60.6450 ± 0.00070.6555.82 1.32+2.0916610 ± 1801.04
aXMMXCS J095417.1+173805.909h54m17.s117°38′05.90.8272 ± 0.00170.8213.65 0.51+0.6210940 ± 3101.42
XMMXCS J095940.7+023113.409h59m40.s7+02°31′13.40.7291 ± 0.00050.7265.02 0.55+0.6825470 ± 900.88
bXMMXCS J105659.5−033728.010h56m59.s5−03°37′28.00.8336 ± 0.00130.8277.57 0.40+0.43291010 ± 1201.57
XMMXCS J112349.4+052955.111h23m49.s4+05°29′55.10.6550 ± 0.00070.6534.62 0.95+1.5517600 ± 2101.05
XMMXCS J113602.9−032943.211h36m02.s9−03°29′43.20.8297 ± 0.00113.32 0.78+1.2021700 ± 1101.06
cXMMXCS J114023.0+660819.011h40m23.s9+66°08′19.00.7855 ± 0.00150.7877.47 0.77+0.9222950 ± 1001.51
aXMMXCS J124312.2−131307.212h43m12.s2−13°13′07.20.7910 ± 0.00140.8014.92 1.54+2.9311790 ± 4601.19
XMMXCS J134305.1−000056.813h43m05.s1−00°00′56.80.6894 ± 0.00110.6784.49 0.57+0.7223920 ± 1701.72
XMMXCS J145009.3+090428.814h50m09.s3+09°04′28.80.6412 ± 0.00070.6033.84 0.55+0.6622630 ± 901.07
dXMMXCS J182132.9+682755.018h21m32.s9+68°27′55.00.8166 ± 0.00110.8294.49 0.56+0.7919860 ± 1301.34
XMMXCS J215221.0−273022.621h52m21.s0−27°30′22.60.8276 ± 0.00110.8232.18 0.45+0.6715530 ± 1500.86
XMMXCS J230247.7+084355.923h02m47.s7+08°43′55.90.7187 ± 0.00140.72105.29 0.50+0.59221010 ± 1301.60
aXMMXCS J235616.4−344144.323h56m16.s4−34°41′44.30.9391 ± 0.00120.9414.57 0.41+0.4810670 ± 2600.91
Table 4.

High-redshift sample (0.5 < z < 1.0): all columns are as explained in Table 3. The superscripts in column one indicate the origin of redshift data when it did not come from our own observations; aare all clusters taken from the XDCP survey and were observed with the VLT–FORS2 spectograph (Nastasi et al. 2014), balso known as MS1054-03 was observed with Keck for 8.6 h (Tran et al. 1999), cwas observed with Keck (Donahue et al. 1999), dalso known as RXJ1821.6+6827 was observed with CFHT, Keck and the 2.2 m telescope at the University of Hawaii (Gioia et al. 2004). The references for the redshifts are as follows: 1. Nastasi et al. (2014), 2. Scharf et al. (1997), 3. Mehrtens et al. (2012), 4. Adami et al. (2011), 5. Šuhada et al. (2011), 6. Bellagamba et al. (2011), 7. Gioia & Luppino (1994), 8. Basilakos et al. (2004), 9. Gioia et al. (2004) and 10. Perlman et al. (2002).

NameRADeczzlitTMembersσvR200
(J2000)(J2000)(keV)(km s−1)(Mpc)
aXMMXCS J000216.1−355633.800h02m16.s1−35°56′33.80.7709 ± 0.00210.7714.83 0.76+1.015131100 ± 1901.77
XMMXCS J005656.6−274031.900h56m56.s6−27°40′31.90.5601 ± 0.00070.5623.30 0.63+0.9415380 ± 600.66
XMMXCS J015241.1−133855.901h52m41.s1−13°38′55.90.8268 ± 0.00100.8233.23 0.31+0.3829840 ± 1501.33
XMMXCS J021734.7−051326.902h17m34.s7−05°13′26.90.6467 ± 0.00120.6542.23 0.44+0.9012620 ± 2101.11
XMMXCS J025006.4−310400.802h50m06.s4−31°04′00.80.9100 ± 0.00240.9024.50 0.88+1.33131120 ± 2601.66
XMMXCS J030205.1−000003.603h02m05.s1−00°00′03.60.6450 ± 0.00070.6555.82 1.32+2.0916610 ± 1801.04
aXMMXCS J095417.1+173805.909h54m17.s117°38′05.90.8272 ± 0.00170.8213.65 0.51+0.6210940 ± 3101.42
XMMXCS J095940.7+023113.409h59m40.s7+02°31′13.40.7291 ± 0.00050.7265.02 0.55+0.6825470 ± 900.88
bXMMXCS J105659.5−033728.010h56m59.s5−03°37′28.00.8336 ± 0.00130.8277.57 0.40+0.43291010 ± 1201.57
XMMXCS J112349.4+052955.111h23m49.s4+05°29′55.10.6550 ± 0.00070.6534.62 0.95+1.5517600 ± 2101.05
XMMXCS J113602.9−032943.211h36m02.s9−03°29′43.20.8297 ± 0.00113.32 0.78+1.2021700 ± 1101.06
cXMMXCS J114023.0+660819.011h40m23.s9+66°08′19.00.7855 ± 0.00150.7877.47 0.77+0.9222950 ± 1001.51
aXMMXCS J124312.2−131307.212h43m12.s2−13°13′07.20.7910 ± 0.00140.8014.92 1.54+2.9311790 ± 4601.19
XMMXCS J134305.1−000056.813h43m05.s1−00°00′56.80.6894 ± 0.00110.6784.49 0.57+0.7223920 ± 1701.72
XMMXCS J145009.3+090428.814h50m09.s3+09°04′28.80.6412 ± 0.00070.6033.84 0.55+0.6622630 ± 901.07
dXMMXCS J182132.9+682755.018h21m32.s9+68°27′55.00.8166 ± 0.00110.8294.49 0.56+0.7919860 ± 1301.34
XMMXCS J215221.0−273022.621h52m21.s0−27°30′22.60.8276 ± 0.00110.8232.18 0.45+0.6715530 ± 1500.86
XMMXCS J230247.7+084355.923h02m47.s7+08°43′55.90.7187 ± 0.00140.72105.29 0.50+0.59221010 ± 1301.60
aXMMXCS J235616.4−344144.323h56m16.s4−34°41′44.30.9391 ± 0.00120.9414.57 0.41+0.4810670 ± 2600.91
NameRADeczzlitTMembersσvR200
(J2000)(J2000)(keV)(km s−1)(Mpc)
aXMMXCS J000216.1−355633.800h02m16.s1−35°56′33.80.7709 ± 0.00210.7714.83 0.76+1.015131100 ± 1901.77
XMMXCS J005656.6−274031.900h56m56.s6−27°40′31.90.5601 ± 0.00070.5623.30 0.63+0.9415380 ± 600.66
XMMXCS J015241.1−133855.901h52m41.s1−13°38′55.90.8268 ± 0.00100.8233.23 0.31+0.3829840 ± 1501.33
XMMXCS J021734.7−051326.902h17m34.s7−05°13′26.90.6467 ± 0.00120.6542.23 0.44+0.9012620 ± 2101.11
XMMXCS J025006.4−310400.802h50m06.s4−31°04′00.80.9100 ± 0.00240.9024.50 0.88+1.33131120 ± 2601.66
XMMXCS J030205.1−000003.603h02m05.s1−00°00′03.60.6450 ± 0.00070.6555.82 1.32+2.0916610 ± 1801.04
aXMMXCS J095417.1+173805.909h54m17.s117°38′05.90.8272 ± 0.00170.8213.65 0.51+0.6210940 ± 3101.42
XMMXCS J095940.7+023113.409h59m40.s7+02°31′13.40.7291 ± 0.00050.7265.02 0.55+0.6825470 ± 900.88
bXMMXCS J105659.5−033728.010h56m59.s5−03°37′28.00.8336 ± 0.00130.8277.57 0.40+0.43291010 ± 1201.57
XMMXCS J112349.4+052955.111h23m49.s4+05°29′55.10.6550 ± 0.00070.6534.62 0.95+1.5517600 ± 2101.05
XMMXCS J113602.9−032943.211h36m02.s9−03°29′43.20.8297 ± 0.00113.32 0.78+1.2021700 ± 1101.06
cXMMXCS J114023.0+660819.011h40m23.s9+66°08′19.00.7855 ± 0.00150.7877.47 0.77+0.9222950 ± 1001.51
aXMMXCS J124312.2−131307.212h43m12.s2−13°13′07.20.7910 ± 0.00140.8014.92 1.54+2.9311790 ± 4601.19
XMMXCS J134305.1−000056.813h43m05.s1−00°00′56.80.6894 ± 0.00110.6784.49 0.57+0.7223920 ± 1701.72
XMMXCS J145009.3+090428.814h50m09.s3+09°04′28.80.6412 ± 0.00070.6033.84 0.55+0.6622630 ± 901.07
dXMMXCS J182132.9+682755.018h21m32.s9+68°27′55.00.8166 ± 0.00110.8294.49 0.56+0.7919860 ± 1301.34
XMMXCS J215221.0−273022.621h52m21.s0−27°30′22.60.8276 ± 0.00110.8232.18 0.45+0.6715530 ± 1500.86
XMMXCS J230247.7+084355.923h02m47.s7+08°43′55.90.7187 ± 0.00140.72105.29 0.50+0.59221010 ± 1301.60
aXMMXCS J235616.4−344144.323h56m16.s4−34°41′44.30.9391 ± 0.00120.9414.57 0.41+0.4810670 ± 2600.91

Fitting the velocity dispersion–temperature relation

To determine the scaling relation between the velocity dispersion and temperature, we fitted a power law of the form
(3)
Here, 5 keV and 1000 km s−1 are the pivot temperature and velocity dispersion, respectively, for our fit. These were chosen to reduce the covariance between the normalization A and the slope B, and for ease of comparison to previous studies. In the above, evolution of the normalization is parametrized as E(z)C, where E(z)=Ωm(1+z)3+ΩΛ describes the redshift evolution of the Hubble parameter. For the self-similar case, B = 0.5 and C = 0 are expected. Similarly to Hilton et al. (2012), the best-fitting values for these parameters were found using Markov Chain Monte Carlo (MCMC) with the Metropolis algorithm. Both orthogonal and bisector regression methods were used. For the orthogonal method, the probability for a given cluster to be drawn from the model scaling relation is
(4)
where r − rmodel is the orthogonal distance of the cluster from the model relation, Δr is the error on the orthogonal distance and S is the intrinsic scatter orthogonal to the model relation. Δr is calculated from the projection in the direction orthogonal to the model line of the ellipse defined by the errors on log σv and log T, chosen according to the position of a given point relative to the model fit line. For the bisector method, the intrinsic scatter and measurement errors are treated independently for each axis. Therefore in the equation for Pmodel, rmodel is replaced by
(5)
and
(6)
where r and Δr are replaced by x, Δx or y, Δy as appropriate. The intrinsic scatter S is replaced by two parameters Sx and Sy. For both methods, the likelihood L of a given model is simply the product of Pmodel for each cluster in the sample, i.e. in the orthogonal case
(7)
where we assume generous, uniform priors on each parameter, as listed in Table 5.
Table 5.

Priors on σvT relation fit parameters.

ParameterUniform priorNotes
A(−5.0,5.0)
B(0.0,2.0)
C(−1.0,1.0)
S(0.01,1.0)Orthogonal method only
Sx(0.01,1.0)Bisector method only
Sy(0.01,1.0)Bisector method only
ParameterUniform priorNotes
A(−5.0,5.0)
B(0.0,2.0)
C(−1.0,1.0)
S(0.01,1.0)Orthogonal method only
Sx(0.01,1.0)Bisector method only
Sy(0.01,1.0)Bisector method only
Table 5.

Priors on σvT relation fit parameters.

ParameterUniform priorNotes
A(−5.0,5.0)
B(0.0,2.0)
C(−1.0,1.0)
S(0.01,1.0)Orthogonal method only
Sx(0.01,1.0)Bisector method only
Sy(0.01,1.0)Bisector method only
ParameterUniform priorNotes
A(−5.0,5.0)
B(0.0,2.0)
C(−1.0,1.0)
S(0.01,1.0)Orthogonal method only
Sx(0.01,1.0)Bisector method only
Sy(0.01,1.0)Bisector method only

RESULTS

Evolution of the slope and intrinsic scatter

For the model given in equation (3), it is assumed that the slope (parameter B) is not evolving with redshift. To test this, the σvT relation was fitted with C = 0 in two redshift bins, 0.0 < z < 0.5 and 0.5 < z < 0.9, with 19 clusters in each bin. The parameters A, B and S were obtained using the MCMC method described above for the high- and low-redshift samples individually. The results for this are shown in Figs 3 and 4. Using the orthogonal method we found B = 1.12 ± 0.41 for the high-redshift sample and B = 0.89 ± 0.16 for the low-redshift sample. However, we found that the slope of the relation for the high-redshift sample is unconstrained if the prior on B is relaxed further. We assume for the remainder of this paper that the slope does not evolve with redshift, though clearly either a larger sample or more accurate measurements of individual clusters are needed to confirm that this is true. The intrinsic scatter is S = 0.05 ± 0.02 for the low-redshift sample and S = 0.08 ± 0.04 for the high-redshift sample. Therefore there is no evidence that the intrinsic scatter varies with redshift.

The σv–T relation assuming no evolution, i.e. C = 0 in equation (3), for low (left – 0.0 < z < 0.5) and high (right – 0.5 < z < 0.9) redshift samples. The solid blue line shows an orthogonal regression fit to the data with the dashed line representing the 95 per cent confidence interval. The dot–dashed line shows a bisector regression fit to the data (see Section 3.2). A model of the form seen in equation (3) was used in the Metropolis algorithm to determine a line of best fit (see Section 4.1). It is interesting to note that the two best-measured systems (XMMXCS J105659.5−033728.0 and XMMXCS J114023.0+660819.0) in the high-redshift subsample are relatively far off the best-fitting relation, with a higher than predicted temperature. Our current observations do not provide good enough spatial resolution or deep enough multicolour photometry to determine the exact reason for this and require further study and re-observations.
Figure 3.

The σvT relation assuming no evolution, i.e. C = 0 in equation (3), for low (left – 0.0 < z < 0.5) and high (right – 0.5 < z < 0.9) redshift samples. The solid blue line shows an orthogonal regression fit to the data with the dashed line representing the 95 per cent confidence interval. The dot–dashed line shows a bisector regression fit to the data (see Section 3.2). A model of the form seen in equation (3) was used in the Metropolis algorithm to determine a line of best fit (see Section 4.1). It is interesting to note that the two best-measured systems (XMMXCS J105659.5−033728.0 and XMMXCS J114023.0+660819.0) in the high-redshift subsample are relatively far off the best-fitting relation, with a higher than predicted temperature. Our current observations do not provide good enough spatial resolution or deep enough multicolour photometry to determine the exact reason for this and require further study and re-observations.

Corner plots for the low (left – 0.0 < z < 0.5) and high (right – 0.5 < z < 0.9) redshift sample showing all the one and two-dimensional projections of the posterior probability distributions of the three parameters when using the orthogonal method. The values of each parameter are given in the top centre. The histograms show the one-dimensional marginalized distribution for each parameter and the other plots show the two-dimensional version, where the contours show 1σ, 2σ and 3σ.
Figure 4.

Corner plots for the low (left – 0.0 < z < 0.5) and high (right – 0.5 < z < 0.9) redshift sample showing all the one and two-dimensional projections of the posterior probability distributions of the three parameters when using the orthogonal method. The values of each parameter are given in the top centre. The histograms show the one-dimensional marginalized distribution for each parameter and the other plots show the two-dimensional version, where the contours show 1σ, 2σ and 3σ.

Evolution of the normalization

To test for the evolution of the normalization (parameter A in equation 3), the low- and high-redshift samples were combined and C was allowed to vary in the MCMC analysis. The results obtained are shown in Figs 6 and 7. We found C = −0.53 ± 0.27, meaning that for a given σv, a higher T is obtained at higher redshift. However, the no evolution relation falls within the 95 per cent confidence interval and therefore we conclude that there is no significant evidence in favour of evolution.

We also applied a statistical test known as the Akaike information criterion (AIC) to determine whether the model with or without evolution (Fig. 5) was preferred. The AIC estimates the quality of each model relative to each other and is therefore a means of model selection. It is defined by (Burnham & Anderson 2002) as
(8)
where L is the maximized likelihood function (equation 7) and k is the number of free parameters. The AIC includes a penalty for using extra parameters as a way to discourage overfitting and rewards goodness of fit based on the likelihood function. Therefore the lower the value of the AIC, the better the fit. For the combined sample with the no evolution model the AIC value was −64.6 and when the fourth parameter for evolution (C) was included this increased to −62.1. Therefore, combining this with the results from the σvT relation fit, it can be concluded that the preferred model is the one with no evolution of the normalization of the scaling relation.
The σv–T relation assuming no evolution, i.e. C = 0 in equation (3), for the combined sample. All lines are as explained in Fig. 3.
Figure 5.

The σvT relation assuming no evolution, i.e. C = 0 in equation (3), for the combined sample. All lines are as explained in Fig. 3.

DISCUSSION

Comparison with previous results

Table 6 and Figs 3–6 present the results of applying the orthogonal and bisector fitting methods to the low redshift, high redshift and combined samples. We see that the bisector and orthogonal method give very similar results especially for our total sample without evolution. Hogg, Bovy & Lang (2010) suggest that the bisector method should be avoided, as by simply finding the difference between a forward and reverse fitting method large systematic errors will be introduced, but it has been widely used for scaling relation measurements in the past and is therefore included for completeness. Results from previous studies of the σvT relation are collected in Table 1. All of these studies, except for Edge & Stewart (1991) and the low-redshift sample of Wu et al. (1998), obtained a slope steeper than the expected self-similar slope of σv ∝ T0.5. We measured B = 0.72 ± 0.12 using the orthogonal fitting method and B = 0.77 ± 0.08 using the bisector fitting method for our combined sample. Therefore both the orthogonal and bisector slopes are in agreement with each other and the previous values in the literature, except for the result obtained by Edge & Stewart (1991) which is only consistent with the orthogonal result. Except for work done by Wu et al. (1998) and Nastasi et al. (2014), all the previous results were obtained for low-redshift samples and no test for evolution was performed. Wu et al. (1998) divided their sample into two groups, z < 0.1 and z ≥ 0.1, and found no significant evolution, however their sample included only four clusters in the redshift range 0.5 < z < 1.0. Nastasi et al. (2014) had a sample of 12 galaxy clusters and found a very large error of more than 50 per cent on their slope. They concluded that their sample size was too small to accurately measure evolution. We conclude that the data presented in this paper – a homogeneous cluster sample that is larger than those used in previous studies at z > 0.5 – are consistent with previous results.

Plots showing the σv–T relation for the combined sample with varying evolution, i.e. C is a free parameter in equation (3). The velocity dispersion is scaled to take into the account the evolution by multiplying by E(z)−C. All lines are as explained in Fig. 3.
Figure 6.

Plots showing the σvT relation for the combined sample with varying evolution, i.e. C is a free parameter in equation (3). The velocity dispersion is scaled to take into the account the evolution by multiplying by E(z)C. All lines are as explained in Fig. 3.

Plot showing the evolution of the normalization of the σv−T relation obtained for the data with the 95 per cent confidence intervals in the dashed lines as compared to the self-similar relation which predicts no evolution shown as the solid line. The black points show our sample data.
Figure 7.

Plot showing the evolution of the normalization of the σvT relation obtained for the data with the 95 per cent confidence intervals in the dashed lines as compared to the self-similar relation which predicts no evolution shown as the solid line. The black points show our sample data.

Table 6.

Best-fitting σvT scaling relation parameters using both the orthogonal and bisector regression methods (see Section 3.2).

MethodParameterLow redshiftHigh redshiftCombinedCombined
(no evolution)(with evolution)
OrthogonalA0.02 ± 0.06− 0.08 ± 0.05− 0.05 ± 0.040.02 ± 0.05
B0.89 ± 0.161.12 ± 0.410.72 ± 0.120.86 ± 0.14
S0.05 ± 0.030.08 ± 0.040.06 ± 0.030.06 ± 0.03
C000− 0.37 ± 0.33
BisectorA0.02 ± 0.04− 0.07 ± 0.03− 0.04 ± 0.020.02 ± 0.04
B0.85 ± 0.111.01 ± 0.170.77 ± 0.080.86 ± 0.09
Sx0.15 ± 0.030.19 ± 0.030.15 ± 0.020.14 ± 0.02
Sy0.07 ± 0.030.12 ± 0.040.09 ± 0.020.09 ± 0.02
C000− 0.49 ± 0.25
MethodParameterLow redshiftHigh redshiftCombinedCombined
(no evolution)(with evolution)
OrthogonalA0.02 ± 0.06− 0.08 ± 0.05− 0.05 ± 0.040.02 ± 0.05
B0.89 ± 0.161.12 ± 0.410.72 ± 0.120.86 ± 0.14
S0.05 ± 0.030.08 ± 0.040.06 ± 0.030.06 ± 0.03
C000− 0.37 ± 0.33
BisectorA0.02 ± 0.04− 0.07 ± 0.03− 0.04 ± 0.020.02 ± 0.04
B0.85 ± 0.111.01 ± 0.170.77 ± 0.080.86 ± 0.09
Sx0.15 ± 0.030.19 ± 0.030.15 ± 0.020.14 ± 0.02
Sy0.07 ± 0.030.12 ± 0.040.09 ± 0.020.09 ± 0.02
C000− 0.49 ± 0.25
Table 6.

Best-fitting σvT scaling relation parameters using both the orthogonal and bisector regression methods (see Section 3.2).

MethodParameterLow redshiftHigh redshiftCombinedCombined
(no evolution)(with evolution)
OrthogonalA0.02 ± 0.06− 0.08 ± 0.05− 0.05 ± 0.040.02 ± 0.05
B0.89 ± 0.161.12 ± 0.410.72 ± 0.120.86 ± 0.14
S0.05 ± 0.030.08 ± 0.040.06 ± 0.030.06 ± 0.03
C000− 0.37 ± 0.33
BisectorA0.02 ± 0.04− 0.07 ± 0.03− 0.04 ± 0.020.02 ± 0.04
B0.85 ± 0.111.01 ± 0.170.77 ± 0.080.86 ± 0.09
Sx0.15 ± 0.030.19 ± 0.030.15 ± 0.020.14 ± 0.02
Sy0.07 ± 0.030.12 ± 0.040.09 ± 0.020.09 ± 0.02
C000− 0.49 ± 0.25
MethodParameterLow redshiftHigh redshiftCombinedCombined
(no evolution)(with evolution)
OrthogonalA0.02 ± 0.06− 0.08 ± 0.05− 0.05 ± 0.040.02 ± 0.05
B0.89 ± 0.161.12 ± 0.410.72 ± 0.120.86 ± 0.14
S0.05 ± 0.030.08 ± 0.040.06 ± 0.030.06 ± 0.03
C000− 0.37 ± 0.33
BisectorA0.02 ± 0.04− 0.07 ± 0.03− 0.04 ± 0.020.02 ± 0.04
B0.85 ± 0.111.01 ± 0.170.77 ± 0.080.86 ± 0.09
Sx0.15 ± 0.030.19 ± 0.030.15 ± 0.020.14 ± 0.02
Sy0.07 ± 0.030.12 ± 0.040.09 ± 0.020.09 ± 0.02
C000− 0.49 ± 0.25

Description of simulations

Comparison to simulations are important for two main reasons. First, we can determine if there is any bias due to sample selection as the simulations provide both a bigger temperature and redshift range. It also allows us to compare different simulation models and learn about the nature of the non-gravitational physics through their effect on the gas temperature. The Millennium Gas Project is a set of hydrodynamical simulations described in Short et al. (2010) which uses the same initial perturbations as the Millennium Simulation (Springel et al. 2005). These simulations include a variety of models, including gravity only; energy injection with radiative cooling; and feedback only. For comparison to the data presented in this paper, the feedback only model (FO) in a volume of 250 h−1 Mpc3 was used. This model includes supernova and AGN feedback using a semi-analytic galaxy formation model. Heating due to supernovae and AGN and the star formation rate are obtained using the model of De Lucia & Blaizot (2007). The AGN feedback model used is described in Bower, McCarthy & Benson (2008), which is dependent on the matter accreted by the central black hole and the efficiency with which the matter is converted to energy near the event horizon, with the upper limit being at 2 per cent of the Eddington rate. As a comparison to the velocity dispersion of the cluster, two proxies were considered, the velocity dispersion of the stars (σStars) and that obtained from the dark matter particles (σDM). The temperatures used from the simulation were spectroscopic-like temperatures (Tsl; Mazzotta et al. 2004). To ensure that only clusters similar to those in our sample were included, we excluded all groups from the simulation with a mass less than 1014 M. We also included a temperature cut, 2 < T(keV) < 11, and a redshift cut, 0 < z < 1, to match our sample. We also compared to the results of the BAHAMAS hydrodynamical simulation (McCarthy et al. 2016 and Caldwell et al. in preparation). Here, a 400 h−1 Mpc3 box is used, with initial conditions based on Planck 2013 cosmological parameters (Planck Collaboration XVI 2014), and both AGN and supernovae feedback models as described by Le Brun et al. (2014). A galaxy mass lower limit of 5 × 109 M and a cluster mass lower limit of 1014 M were implemented. This simulation reproduces a large number of X-ray, SZ and optical scaling relations of groups and clusters. However, unlike previous simulations, the new simulation also reproduces the observed galaxy stellar mass function remarkably well over a wide range of stellar masses. The velocity dispersion is traced by galaxies and is calculated using the gapper technique described by Beers et al. (1990). The temperatures used from the simulation were spectroscopic (TS). We note that spectroscopic-like temperatures, as used in both simulations, are most robust at T > 2 keV, where the bremsstrahlung mechanism dominates (Short et al. 2010). Therefore, while we have applied mass, temperature and redshift cuts to the simulated cluster catalogues that are a reasonable match to our observed sample, the correspondence is not exact, as six of the observed clusters have T < 2 keV. This is a compromise aimed at limiting the potential impact of low-mass clusters with less reliable temperature measurements in the simulations. Matching the velocity dispersions to the simulations, however, is not as straightforward, and is outside the scope of this paper, but for completeness we briefly discuss causes of bias identified in previous studies. Old, Gray & Pearce (2013) studied the recovery of velocity dispersions from simulation data and explored how sample selection can impact the measurements and cause a bias. They introduced I-band magnitude limits and found that the velocity dispersion recovered from the haloes was systematically higher than that from the galaxies. When this sample was further limited to just the brightest galaxies, this discrepancy was enhanced. They suggest that the reason for this is that dynamical friction greatly affects the velocity of the galaxies, and therefore to reduce this bias a strictly magnitude-limited sample should be avoided. Old et al. (2013) also calculated the velocity dispersion over different radial distances to see how it varied as a function of distance from the cluster centre. They found that the velocity dispersion was sensitive to the radius in which it was measured with a difference of 10 per cent in measurements being found between 0.5 and 1 R200. Sifón et al. (2016) also studied the impact of the choice of radius on velocity dispersion measurements, and found that the bias is negligible for measurements that sample beyond 0.7 R200 (see their fig. 4). Since our observations sample out to at least 0.7 R200 for all clusters, and beyond R200 for more than half the sample, we expect our velocity dispersion measurements to be unaffected by this source of bias. We used these two sets of simulations to test our orthogonal fitting methods both with and without evolution.

Comparison with simulations - fitting with no evolution

The orthogonal fitting method described in Section 3.2 was applied to both sets of simulations with C = 0. The parameters A, B and S for both the Millennium Gas Project and BAHAMAS simulations are shown in Table 7. The σvT relation for the Millennium Gas Project with the two different σv proxies are shown in Fig. 8. The slope is slightly steeper for the stars (B = 0.62 ± 0.01) than for the dark matter (B = 0.55 ± 0.08) but both are consistent with previous studies of the σvT relation and the results obtained from our data. Fig. 9 shows that the orthogonal fit to the full BAHAMAS sample systematically overestimates the average velocity dispersion at T > 5 keV. This may be due, in part, to the model not being a complete description of the data, as the scatter appears to vary with temperature. This is not captured in our orthogonal regression model (equation 4), i.e. S is assumed to be constant with both T and z. A comparison was made to a fit performed by Caldwell et al. (in preparation) to the BAHAMAS data. In this method the σ2kT relation was derived by first parametrically determining the mean functions and redshift evolution of velocity dispersion and temperature, separately, with respect to mass (M500 critical). The velocity dispersion measurements are averaged in 0.25 log10(M500c) mass bins, to avoid biases from high cluster counts and the large scatter seen at low T (Fig. 9). These mean values of velocity dispersion and mass are fit with a power law, and the temperature relation is derived with the same method. This was converted into the σvT relation (with slope B = 0.545) that is plotted in Fig. 9. This method provides a better fit to BAHAMAS data over the full temperature range than the orthogonal method. We use the slope obtained using the Caldwell et al. method for further studies and comparisons. We examine the effect of bias in the recovered slope on our results in Section 5.4. We note that there is no single method which gives the underlying ‘true’ scaling relation in the presence of errors on both variables and intrinsic scatter: the recovered slope and normalization depend upon the details of the method used. The fitting of a scaling relation is also affected by the selection processes used when determining your sample and for this study this has not been corrected for. Since our clusters are selected on X-ray luminosity rather than temperature or velocity dispersion we believe that our values are not biased and therefore selection effects will not have as big an impact on our results.

The σv–T relation for the Millennium Gas Project simulations using dark matter and stars as proxies for the velocity dispersion. The blue dots are the data obtained from the simulation and the solid black line shows the fit using the orthogonal regression method. The slope is slightly steeper for the stars (B = 0.62 ± 0.01) than for the dark matter (B = 0.55 ± 0.08) but both are consistent with previous studies of the σv–T relation and the results obtained from our data.
Figure 8.

The σvT relation for the Millennium Gas Project simulations using dark matter and stars as proxies for the velocity dispersion. The blue dots are the data obtained from the simulation and the solid black line shows the fit using the orthogonal regression method. The slope is slightly steeper for the stars (B = 0.62 ± 0.01) than for the dark matter (B = 0.55 ± 0.08) but both are consistent with previous studies of the σvT relation and the results obtained from our data.

Table 7.

Best-fitting values for the parameters in equation (4) (slope, intercept and scatter) for the various models obtained from simulations without evolution. For the Millennium Gas Project, we use dark matter (DM) and stars as the tracers for the velocity dispersion. The BAHAMAS simulation uses galaxies. The Millennium Gas simulations use spectroscopic-like temperatures (Tsl) and the BAHAMAS simulation use spectroscopic temperatures (Ts). Caldwell et al. (in preparation) present a different method for determining the σvT relation as discussed in Section 5.2, the results of which are also shown below.

SimulationσtracerTmodelABS
Millennium GasDMTsl−0.011 ± 0.0020.553 ± 0.0080.028 ± 0.001
Millennium GasStarsTsl−0.034 ± 0.0030.621 ± 0.0100.034 ± 0.001
BAHAMASGalaxiesTs−0.055 ± 0.0030.848 ± 0.0120.055 ± 0.001
BAHAMAS (Caldwell et al. in preparation)GalaxiesTs−0.1330.545
SimulationσtracerTmodelABS
Millennium GasDMTsl−0.011 ± 0.0020.553 ± 0.0080.028 ± 0.001
Millennium GasStarsTsl−0.034 ± 0.0030.621 ± 0.0100.034 ± 0.001
BAHAMASGalaxiesTs−0.055 ± 0.0030.848 ± 0.0120.055 ± 0.001
BAHAMAS (Caldwell et al. in preparation)GalaxiesTs−0.1330.545
Table 7.

Best-fitting values for the parameters in equation (4) (slope, intercept and scatter) for the various models obtained from simulations without evolution. For the Millennium Gas Project, we use dark matter (DM) and stars as the tracers for the velocity dispersion. The BAHAMAS simulation uses galaxies. The Millennium Gas simulations use spectroscopic-like temperatures (Tsl) and the BAHAMAS simulation use spectroscopic temperatures (Ts). Caldwell et al. (in preparation) present a different method for determining the σvT relation as discussed in Section 5.2, the results of which are also shown below.

SimulationσtracerTmodelABS
Millennium GasDMTsl−0.011 ± 0.0020.553 ± 0.0080.028 ± 0.001
Millennium GasStarsTsl−0.034 ± 0.0030.621 ± 0.0100.034 ± 0.001
BAHAMASGalaxiesTs−0.055 ± 0.0030.848 ± 0.0120.055 ± 0.001
BAHAMAS (Caldwell et al. in preparation)GalaxiesTs−0.1330.545
SimulationσtracerTmodelABS
Millennium GasDMTsl−0.011 ± 0.0020.553 ± 0.0080.028 ± 0.001
Millennium GasStarsTsl−0.034 ± 0.0030.621 ± 0.0100.034 ± 0.001
BAHAMASGalaxiesTs−0.055 ± 0.0030.848 ± 0.0120.055 ± 0.001
BAHAMAS (Caldwell et al. in preparation)GalaxiesTs−0.1330.545
Plot showing the σv–T relation for the BAHAMAS simulation data. The blue circles represent the data points from the simulation. The solid line is the fit obtained using the orthogonal method and the dashed line is the fit obtained using the method described in Section 5.2 and performed by Caldwell et al. (in preparation). From this it can be seen that the orthogonal fit overestimates the velocity dispersion at TX > 5 keV (see Section 5.2).
Figure 9.

Plot showing the σvT relation for the BAHAMAS simulation data. The blue circles represent the data points from the simulation. The solid line is the fit obtained using the orthogonal method and the dashed line is the fit obtained using the method described in Section 5.2 and performed by Caldwell et al. (in preparation). From this it can be seen that the orthogonal fit overestimates the velocity dispersion at TX > 5 keV (see Section 5.2).

The effect of biased slope measurements on the evolution of the normalization

Having seen, using the BAHAMAS simulation, that the slope recovered using the orthogonal regression method may be biased high, we now discuss the potential impact of a biased slope measurement on our conclusions regarding the observed cluster sample in Section 5.1.

To investigate this, we generated 1000 mock samples (each containing 38 clusters) from the BAHAMAS simulation with the same temperature distribution as the observed sample, and applied the orthogonal regression method. Fig. 10 shows the distribution of recovered slope values. The average is B = 0.69 ± 0.13, which is 2σ higher than the slope obtained from the fit performed by Caldwell et al. (in preparation, Section 5.2). Therefore, if the BAHAMAS sample is representative of the real cluster population, then we would conclude that the slope we have measured for the observed cluster sample is biased high. To check if a biased slope estimate affects our conclusions regarding the lack of significant evidence for evolution of the normalization of the relation (Section 4.2), we fixed the slope to B = 0.545 and re-ran the orthogonal fit for the observed cluster sample. We found C = 0.15 ± 0.28, which is consistent with no evolution (Fig. 11). Therefore, even if the slope value of B = 0.86 ± 0.14 that we measured was biased high for any reason, this does not affect our conclusion that we do not see significant evidence in favour of evolution.

This histogram shows the probability of getting a specific value for the slope of the best-fitting σ–T relation given a set T distribution. We chose various subsamples from the BAHAMAS simulation which had the same T distribution as our sample and calculated the slope for each. The mean slope obtained is B = 0.69 ± 0.13, which is within 2σ of the value obtained from Caldwell et al (in preparation), so there is a slight bias from the distribution of the sample.
Figure 10.

This histogram shows the probability of getting a specific value for the slope of the best-fitting σ–T relation given a set T distribution. We chose various subsamples from the BAHAMAS simulation which had the same T distribution as our sample and calculated the slope for each. The mean slope obtained is B = 0.69 ± 0.13, which is within 2σ of the value obtained from Caldwell et al (in preparation), so there is a slight bias from the distribution of the sample.

Plot showing the evolution of the normalization of the σv–T relation, with B = 0.545, obtained for the observed cluster sample with the 95 per cent confidence intervals in the dashed lines, as compared to the self-similar relation which predicts no evolution shown as the solid line. The black points show the measurements for the clusters in our sample.
Figure 11.

Plot showing the evolution of the normalization of the σvT relation, with B = 0.545, obtained for the observed cluster sample with the 95 per cent confidence intervals in the dashed lines, as compared to the self-similar relation which predicts no evolution shown as the solid line. The black points show the measurements for the clusters in our sample.

Comparison with simulations – fitting with evolution

We now investigate evolution in the normalization of the σvT relation in the simulations by fitting for the value of C, as we did for the observed sample (see Section 4.2). The results are shown in Table 8 and graphically in Fig. 12. The BAHAMAS simulation was tested both with a slope that was allowed to vary and a fixed slope. Although, in Section 5.4 we showed that the overestimated slope does not affect evolution, we included a fixed slope for further comparison. Both were found to be consistent with zero evolution, further proving that the biased slope does not affect evolution. However, the simulations from the Millennium Gas Project show small but significant positive evolution (C = 0.273 ± 0.013 for σStars and Tsl). To see the reason for this, we can re-write the σvT relation in terms of the σvM and TM relations, where M is the cluster mass (see e.g. Maughan 2014). We define
(9)
where
(10)
Here, A, B and C have the same meaning as before, and the subscripts indicate the corresponding relation (e.g. BTM indicates the slope of the TM relation). If we set CσvM=1/3, CTM = 2/3 and BσvT=1/2 as predicted by the self-similar relation, then we obtain CσvT=0 as expected. We performed fits to determine the values of CσvM, CTM and BσvT in the Millennium Gas simulation at z = 0 and 0.5. We found that CσvM=1/3 when using either σStars or σDM as the measure of σv, and that BσvT varied from 0.55 to 0.6 (depending on whether spectroscopic-like or mass-weighted temperature estimates were used), which is slightly higher than the self-similar value, but not by enough to explain the positive evolution measured in the σvT relation. This leads to the conclusion that the evolution is driven by the value of CTM, and it was found that the measured value for the dark matter was CTM = 2/3 as expected, but that this decreased to values between 0 and 0.2 for the gas. Therefore in the Millennium Gas simulation, the lack of redshift evolution in the T − M relation drives the positive evolution in the σvT relation. The most likely explanation for the lack of redshift evolution in the T − M relation in the Millennium Gas simulation is the absence of radiative cooling. When both cooling and feedback are included in simulations (as in BAHAMAS), the feedback acts as a regulation mechanism, heating the surrounding dense gas and expelling it from the cluster core. This in turn leads to higher entropy gas flowing inwards. In the Millennium Gas simulation, the feedback model heats the gas and directly increases its entropy, which is eventually distributed throughout the cluster. This builds up over time as more and more energy is pumped into the gas from the growing black holes, and has the effect of slowing down the evolution of the T − M relation (compared to the evolution expected due to the decreasing background density with redshift). This in turn leads to the positive evolution of the σvT relation. It is likely that the more sophisticated feedback model used in BAHAMAS, where the entropy evolution is driven by radiative cooling, is the more realistic of the two.
We compared the evolution of the normalization of σv–T relation of the Millennium Gas and BAHAMAS simulations with the self-similar relation and that found from our data using a fixed slope. The solid line shows the line representing the self-similar relation i.e. C = 0, the dot–dashed line represents the BAHAMAS simulation results with a fixed B = 0.545 to avoid bias and the vertically dashed line represents the result from the Millennium Gas simulation. The blue dashed line and black points are our orthogonal fit and observed sample, respectively.
Figure 12.

We compared the evolution of the normalization of σvT relation of the Millennium Gas and BAHAMAS simulations with the self-similar relation and that found from our data using a fixed slope. The solid line shows the line representing the self-similar relation i.e. C = 0, the dot–dashed line represents the BAHAMAS simulation results with a fixed B = 0.545 to avoid bias and the vertically dashed line represents the result from the Millennium Gas simulation. The blue dashed line and black points are our orthogonal fit and observed sample, respectively.

Table 8.

Best-fitting values for the parameters in equation (4) (slope, intercept, scatter and evolution) for the various models obtained from simulations. All abbreviations are as in Table 7.

SimulationσtracerTmodelABSC
Millennium GasDMTsl−0.031 ± 0.0020.551 ± 0.0060.0220 ± 0.00100.371 ± 0.014
Millennium GasStarsTsl−0.056 ± 0.0020.619 ± 0.0090.0295 ± 0.00100.397 ± 0.019
BAHAMAS (fixed B)GalaxiesTs−0.135 ± 0.0020.5450.0390 ± 0.00100.046 ± 0.016
BAHAMAS (varying B)GalaxiesTs−0.071 ± 0.0050.779 ± 0.0140.0570 ± 0.0010−0.029 ± 0.024
SimulationσtracerTmodelABSC
Millennium GasDMTsl−0.031 ± 0.0020.551 ± 0.0060.0220 ± 0.00100.371 ± 0.014
Millennium GasStarsTsl−0.056 ± 0.0020.619 ± 0.0090.0295 ± 0.00100.397 ± 0.019
BAHAMAS (fixed B)GalaxiesTs−0.135 ± 0.0020.5450.0390 ± 0.00100.046 ± 0.016
BAHAMAS (varying B)GalaxiesTs−0.071 ± 0.0050.779 ± 0.0140.0570 ± 0.0010−0.029 ± 0.024
Table 8.

Best-fitting values for the parameters in equation (4) (slope, intercept, scatter and evolution) for the various models obtained from simulations. All abbreviations are as in Table 7.

SimulationσtracerTmodelABSC
Millennium GasDMTsl−0.031 ± 0.0020.551 ± 0.0060.0220 ± 0.00100.371 ± 0.014
Millennium GasStarsTsl−0.056 ± 0.0020.619 ± 0.0090.0295 ± 0.00100.397 ± 0.019
BAHAMAS (fixed B)GalaxiesTs−0.135 ± 0.0020.5450.0390 ± 0.00100.046 ± 0.016
BAHAMAS (varying B)GalaxiesTs−0.071 ± 0.0050.779 ± 0.0140.0570 ± 0.0010−0.029 ± 0.024
SimulationσtracerTmodelABSC
Millennium GasDMTsl−0.031 ± 0.0020.551 ± 0.0060.0220 ± 0.00100.371 ± 0.014
Millennium GasStarsTsl−0.056 ± 0.0020.619 ± 0.0090.0295 ± 0.00100.397 ± 0.019
BAHAMAS (fixed B)GalaxiesTs−0.135 ± 0.0020.5450.0390 ± 0.00100.046 ± 0.016
BAHAMAS (varying B)GalaxiesTs−0.071 ± 0.0050.779 ± 0.0140.0570 ± 0.0010−0.029 ± 0.024

CONCLUSIONS

We have studied the evolution of the velocity dispersion–temperature (σvT) relation using a cluster sample spanning the range 0.0 < z < 1.0 drawn from XCS. This work improves upon previous studies in terms of the use of a homogeneous cluster sample and the number of z > 0.5 clusters included. We present new redshift and velocity dispersion measurements based on Gemini data for 12 such z > 0.5 XCS clusters. We used an orthogonal regression method to measure the normalization, slope and intrinsic scatter of the σvT relation for two subsamples: 19 clusters at z < 0.5, and 19 clusters with z > 0.5. In both cases, we found the slope of the relation to be consistent with the findings of previous studies, i.e. slightly steeper than expected from self-similarity. Under the assumption that the slope of the relation does not evolve with redshift, we measured the evolution of the normalization of the relation using the complete sample of 38 clusters. We found this to be slightly negative but not significantly different from the self-similar solution (σv ∝ T0.86±0.14E(z)−0.37±0.33). Moreover, a no evolution model is the preferred choice when considering the AIC. We applied the same scaling relation analysis methods to the BAHAMAS and Millennium Gas cosmological hydrodynamical simulations. The σvT relation does not evolve in BAHAMAS, in agreement with our findings for the observed cluster sample. However, positive evolution is seen in the Millennium Gas simulation. The difference is most likely due to the inclusion of self-consistent modelling of radiative cooling in BAHAMAS, which is absent in the Millennium Gas simulation. This leads to a very slowly evolving TM relation in the Millennium Gas simulation, which in turn drives the positive evolution of the σvT relation. While this work has improved upon previous studies in terms of the number of high-redshift clusters included, we note that the uncertainties on the scaling relation parameters are still rather large, and a combination of better measurements of individual cluster properties and a larger sample are required to make further progress. Future studies will look at implementing the Bayesian method described by Kelly (2007) to account for intrinsic scatter and measurement errors and looking at possible selection effects.

SW and MH acknowledge financial support from the National Research Foundation and SKA South Africa. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. Based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), Ministério da Ciência, Tecnologia e Inovação (Brazil) and Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina). PAT acknowledges support from the Science and Technology Facilities Council (grant number ST/L000652/1). JPS gratefully acknowledges support from a Hintze Research Fellowship. STK acknowledges support from STFC (grant number ST/L000768/1).

1

This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

2

The methodology described in Section 3 was applied to these six clusters before they were excluded from our study.

3

iraf is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.

REFERENCES

Abazajian
K. N.
et al.
2009
ApJS
182
543

Adami
C.
et al.
2011
A&A
526
A18

Arnaud
M.
Pointecouteau
E.
Pratt
G. W.
2005
A&A
441
893

Basilakos
S.
Plionis
M.
Georgakakis
A.
Georgantopoulos
I.
Gaga
T.
Kolokotronis
V.
Stewart
G. C.
2004
MNRAS
351
989

Beers
T. C.
Flynn
K.
Gebhardt
K.
1990
AJ
100
32

Bellagamba
F.
Maturi
M.
Hamana
T.
Meneghetti
M.
Miyazaki
S.
Moscardini
L.
2011
MNRAS
413
1145

Bird
C. M.
Mushotzky
R. F.
Metzler
C. A.
1995
ApJ
453
40

Bower
R. G.
1997
MNRAS
288
355

Bower
R. G.
McCarthy
I. G.
Benson
A. J.
2008
MNRAS
390
1399

Branchesi
M.
Gioia
I. M.
Fanti
C.
Fanti
R.
2007
A&A
472
739

Bryan
G. L.
Norman
M. L.
1998
ApJ
495
80

Burenin
R. A.
Vikhlinin
A.
Hornstrup
A.
Ebeling
H.
Quintana
H.
Mescheryakov
A.
2007
ApJS
172
561

Burnham
K. P.
Anderson
D. R.
2002
Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach
2 edn
Springer-Verlag
New York

Cappi
A.
Held
E. V.
Marano
B.
1998
A&AS
129
31

Cash
W.
1979
ApJ
228
939

Clerc
N.
Sadibekova
T.
Pierre
M.
Pacaud
F.
Le Fèvre
J.-P.
Adami
C.
Altieri
B.
Valtchanov
I.
2012
MNRAS
423
3561

Clerc
N.
et al.
2014
MNRAS
444
2723

Cohn
J. D.
White
M.
2005
Astropart. Phys.
24
316

Collins
C. A.
et al.
2009
Nature
458
603

Crawford
S. M.
Wirth
G. D.
Bershady
M. A.
2014
ApJ
786
30

De Lucia
G.
Blaizot
J.
2007
MNRAS
375
2

De Propris
R.
et al.
2002
MNRAS
329
87

Dickey
J. M.
1990
Thronson
H. A.
Jr
Shull
J. M.
Astrophysics and Space Science Library Vol. 161, The Interstellar Medium in Galaxies
Kluwer
Dordrecht
473

Donahue
M.
Voit
G. M.
Scharf
C. A.
Gioia
I. M.
Mullis
C. R.
Hughes
J. P.
Stocke
J. T.
1999
ApJ
527
525

Edge
A. C.
Stewart
G. C.
1991
MNRAS
252
428

Ettori
S.
Tozzi
P.
Borgani
S.
Rosati
P.
2004
A&A
417
13

Fadda
D.
Girardi
M.
Giuricin
G.
Mardirossian
F.
Mezzetti
M.
1996
ApJ
473
670

Fassbender
R.
et al.
2011
New J. Phys.
13
125014

Finn
R. A.
et al.
2005
ApJ
630
206

Finoguenov
A.
et al.
2007
ApJS
172
182

Gioia
I. M.
Luppino
G. A.
1994
ApJS
94
583

Gioia
I. M.
Wolter
A.
Mullis
C. R.
Henry
J. P.
Böhringer
H.
Briel
U. G.
2004
A&A
428
867

Girardi
M.
Fadda
D.
Giuricin
G.
Mardirossian
F.
Mezzetti
M.
Biviano
A.
1996
ApJ
457
61

Glazebrook
K.
Bland-Hawthorn
J.
2001
PASP
113
197

Hao
J.
et al.
2010
ApJS
191
254

Hasselfield
M.
et al.
2013
J. Cosmol. Astropart. Phys.
7
8

Hennawi
J. F.
et al.
2008
AJ
135
664

Hilton
M.
et al.
2009
ApJ
697
436

Hilton
M.
et al.
2010
ApJ
718
133

Hilton
M.
et al.
2012
MNRAS
424
2086

Hoekstra
H.
Herbonnet
R.
Muzzin
A.
Babul
A.
Mahdavi
A.
Viola
M.
Cacciato
M.
2015
MNRAS
449
685

Hogg
D. W.
Bovy
J.
Lang
D.
2010
preprint (arXiv:1008.4686)

Kaiser
N.
1986
MNRAS
222
323

Kay
S. T.
da Silva
A. C.
Aghanim
N.
Blanchard
A.
Liddle
A. R.
Puget
J.-L.
Sadat
R.
Thomas
P. A.
2007
MNRAS
377
317

Kelly
B. C.
2007
ApJ
665
1489

Koester
B. P.
et al.
2007
ApJ
660
239

Kravtsov
A. V.
Borgani
S.
2012
ARA&A
50
353

Kurtz
M. J.
Mink
D. J.
1998
PASP
110
934

Le Brun
A. M. C.
McCarthy
I. G.
Schaye
J.
Ponman
T. J.
2014
MNRAS
441
1270

Lloyd-Davies
E. J.
et al.
2011
MNRAS
418
14

Lubin
L. M.
Bahcall
N. A.
1993
ApJ
415
L17

Lumb
D. H.
et al.
2004
A&A
420
853

Magnelli
B.
Elbaz
D.
Chary
R. R.
Dickinson
M.
Le Borgne
D.
Frayer
D. T.
Willmer
C. N. A.
2009
A&A
496
57

Mann
A. W.
Ebeling
H.
2012
MNRAS
420
2120

Maughan
B. J.
2014
MNRAS
437
1171

Maughan
B. J.
Jones
L. R.
Ebeling
H.
Scharf
C.
2006
MNRAS
365
509

Maughan
B. J.
Giles
P. A.
Randall
S. W.
Jones
C.
Forman
W. R.
2012
MNRAS
421
1583

Mazzotta
P.
Rasia
E.
Moscardini
L.
Tormen
G.
2004
MNRAS
354
10

McCarthy
I. G.
Schaye
J.
Bird
S.
Le Brun
A. M. C.
2016
MNRAS
preprint (arXiv:1603.02702)

Mehrtens
N.
et al.
2012
MNRAS
423
1024

Mewe
R.
Schrijver
C. J.
1986
A&A
169
178

Mulchaey
J. S.
Lubin
L. M.
Fassnacht
C.
Rosati
P.
Jeltema
T. E.
2006
ApJ
646
133

Mullis
C. R.
et al.
2003
ApJ
594
154

Nastasi
A.
et al.
2014
A&A
564
A17

Old
L.
Gray
M. E.
Pearce
F. R.
2013
MNRAS
434
2606

Ortiz-Gil
A.
Guzzo
L.
Schuecker
P.
Böhringer
H.
Collins
C. A.
2004
MNRAS
348
325

Perlman
E. S.
Horner
D. J.
Jones
L. R.
Scharf
C. A.
Ebeling
H.
Wegner
G.
Malkan
M.
2002
ApJS
140
265

Planck Collaboration XVI
2014
A&A
571
A16

Planck Collaboration XXIV
2015
preprint (arXiv:1502.01597)

Ponman
T. J.
Bourner
P. D. J.
Ebeling
H.
Böhringer
H.
1996
MNRAS
283
690

Poole
G. B.
Babul
A.
McCarthy
I. G.
Fardal
M. A.
Bildfell
C. J.
Quinn
T.
Mahdavi
A.
2007
MNRAS
380
437

Quintana
H.
Melnick
J.
1982
AJ
87
972

Randall
S. W.
Sarazin
C. L.
Ricker
P. M.
2002
ApJ
577
579

Reichardt
C. L.
et al.
2013
ApJ
763
127

Ritchie
B. W.
Thomas
P. A.
2002
MNRAS
329
675

Rozo
E.
et al.
2009
ApJ
699
768

Sahlén
M.
et al.
2009
MNRAS
397
577

Sakelliou
I.
Merrifield
M. R.
1998
MNRAS
293
489

Schafer
R. A.
1991
XSPEC, An X-ray Spectral Fitting Package: Version 2 of the User's Guide
European Space Agency
Paris

Scharf
C. A.
Ebeling
H.
Perlman
E.
Malkan
M.
Wegner
G.
1997
ApJ
477
79

Short
C. J.
Thomas
P. A.
Young
O. E.
Pearce
F. R.
Jenkins
A.
Muanwong
O.
2010
MNRAS
408
2213

Sifón
C.
et al.
2013
ApJ
772
25

Sifón
C.
et al.
2016
MNRAS
461
248

Silverman
J. D.
et al.
2005
ApJ
618
123

Springel
V.
et al.
2005
Nature
435
629

Stott
J. P.
et al.
2010
ApJ
718
23

Struble
M. F.
Rood
H. J.
1999
ApJS
125
35

Šuhada
R.
et al.
2011
A&A
530
A110

Takey
A.
Schwope
A.
Lamer
G.
2011
A&A
534
A120

Takey
A.
Schwope
A.
Lamer
G.
2013
A&A
558
A75

Tonry
J.
Davis
M.
1979
AJ
84
1511

Tran
K.-V. H.
Kelson
D. D.
van Dokkum
P.
Franx
M.
Illingworth
G. D.
Magee
D.
1999
ApJ
522
39

Vikhlinin
A.
McNamara
B. R.
Forman
W.
Jones
C.
Quintana
H.
Hornstrup
A.
1998
ApJ
502
558

Vikhlinin
A.
van Speybroeck
L.
Markevitch
M.
Forman
W. R.
Grego
L.
2002
ApJ
578
L107

Vikhlinin
A.
Kravtsov
A.
Forman
W.
Jones
C.
Markevitch
M.
Murray
S. S.
Van Speybroeck
L.
2006
ApJ
640
691

Vikhlinin
A.
et al.
2009
ApJ
692
1060

Voit
G. M.
2005
Rev. Mod. Phys.
77
207

von der Linden
A.
et al.
2014
MNRAS
439
2

White
D. A.
Jones
C.
Forman
W.
1997
MNRAS
292
419

Wu
X.-P.
Fang
L.-Z.
Xu
W.
1998
A&A
338
813

Wu
X.-P.
Xue
Y.-J.
Fang
L.-Z.
1999
ApJ
524
22

Xue
Y.-J.
Wu
X.-P.
2000
ApJ
538
65

Yoon
J. H.
Schawinski
K.
Sheen
Y.-K.
Ree
C. H.
Yi
S. K.
2008
ApJS
176
414

SUPPORTING INFORMATION

Additional Supporting Information may be found in the online version of this article:

APPENDIX B. REDSHIFT CATALOGUE.

APPENDIX C. VELOCITY HISTOGRAMS.

(Supplementary Data).

Please note: Oxford University Press is not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.

APPENDIX A: OBSERVATIONS LOG

Table A1.

Spectroscopic observations log. For all observations the R400 grating and the OG515 filter was used.

Cluster NameMaskSlitsAirmass rangeObservation dateFrames(s)Seeing (arcsec)
XMMXCS J005656.6−274031.9GS-2012B-Q-011-03331.222012-10-191 × 1830
1.48–1.762012-10-162 × 1830
1.01–1.042012-10-153 × 18300.76–0.80
GS-2012B-Q-011-04351.05–1.352012-11-144 × 1830
XMMXCS J015241.1−133855.9GS-2011B-Q-050-01331.05–1.212011-12-026 × 1830
GS-2011B-Q-050-02341.05–1.652011-12-036 × 1830
XMMXCS J021734.7−051326.9GS-2012B-Q-011-06341.14–1.482012-12-054 × 1830
XMMXCS J025006.4−310400.8GS-2012B-Q-011-09321.11–1.202012-11-242 × 1830
1.01–1.192012-11-214 × 1830
GS-2010B-Q-046-06351.06–1.122010-11-142 × 18300.50–0.60
1.06–1.442010-11-135 × 1830, 1 × 7621
XMMXCS J030205.1−000003.6GS-2011B-Q-050-03321.172011-12-011 × 1830
1.182011-11-201 × 1098
1.17–1.452011-11-184 × 1830
GS-2011B-Q-050-04321.322011-12-311 × 1098
1.23–1.742011-12-304 × 18300.85–1.40
GS-2011B-Q-050-05331.27–1.572012-01-172 × 18300.7
XMMXCS J095940.7+023113.4GS-2010B-Q-046-02351.19–1.232011-01-093 × 1830
1.19–1.252011-01-084 × 1830
GS-2012A-Q-46-01351.19–1.292012-03-184 × 1830
1.19–1.232012-03-022 × 1830
GS-2012A-Q-46-02341.20–1.352012-03-273 × 1830
1.23–1.462012-03-233 × 18300.8
GS-2012A-Q-46-03341.21–1.542012-03-226 × 18300.65–0.70
XMMXCS J112349.3+052956.8GS-2012A-Q-46-05331.23–1.332012-04-225 × 1830
1.472012-04-211 × 1830
GS-2012A-Q-46-06321.25–1.652012-05-154 × 18300.63–0.76
1.45–1.662012-04-222 × 1830
GS-2010B-Q-046-03331.262011-01-311 × 1830
1.35–1.642011-01-292 × 1525, 1 × 975
1.23–1.242011-01-272 × 1830
XMMXCS J113602.9–032943.2GS-2012A-Q-46-07361.142012-05-241× 1830
1.122012-05-231 × 1830
1.12–1.162012-05-203 × 1830
1.122012-05-191 × 1830
GS-2012A-Q-46-08331.48–1.762012-07-152 × 1830
1.41–1.802012-07-113 × 18300.50–0.70
1.52012-07-101 × 1830
XMMXCS J134305.1−000056.8GS-2012A-Q-46-10361.16–1.232012-05-244 × 1830
1.242012-05-211 × 1830
GS-2012A-Q-46-11341.252012-07-101 × 1830
1.16–1.192012-07-092 × 1830
1.22012-07-061 × 1830
1.54–1.842012-06-222 × 1830
XMMXCS J145009.3+090428.8GN2012A-Q-070-05321.02–1.052012-07-092 × 18001.15
1.11–1.622012-06-264 × 18000.84–0.98
GN2012A-Q-070-06341.02–1.042012-07-072 × 1800
1.09–1.172012-07-062 × 1800
1.48–1.792012-06-272 × 1800
GN2012A-Q-070-07331.22–1.592012-07-223 × 1800
1.04–1.162012-07-083 × 18001
XMMXCS J215221.0–273022.6GS-2010B-Q-046-04361.14–1.242010-11-122 × 1830
1.02–1.212010-09-144 × 1830
GS-2011B-Q-050-06341.07–1.152011-10-052 × 1830
1.12–1.562011-09-184 × 18300.60–1.00
GS-2011B-Q-050-07341.00–1.102011-10-244 × 1830
1.05–1.122011-10-162 × 1830
XMMXCS J230247.7+084355.9GN2012A-Q-070-10341.372012-08-081 × 18000.60–0.68
1.02–1.112012-07-305 × 18000.43–0.86
GN2012A-Q-070-11331.18–1.312012-08-132 × 1800
1.02–1.082012-08-093 × 18000.60–0.68
1.192012-08-081 × 18001
Cluster NameMaskSlitsAirmass rangeObservation dateFrames(s)Seeing (arcsec)
XMMXCS J005656.6−274031.9GS-2012B-Q-011-03331.222012-10-191 × 1830
1.48–1.762012-10-162 × 1830
1.01–1.042012-10-153 × 18300.76–0.80
GS-2012B-Q-011-04351.05–1.352012-11-144 × 1830
XMMXCS J015241.1−133855.9GS-2011B-Q-050-01331.05–1.212011-12-026 × 1830
GS-2011B-Q-050-02341.05–1.652011-12-036 × 1830
XMMXCS J021734.7−051326.9GS-2012B-Q-011-06341.14–1.482012-12-054 × 1830
XMMXCS J025006.4−310400.8GS-2012B-Q-011-09321.11–1.202012-11-242 × 1830
1.01–1.192012-11-214 × 1830
GS-2010B-Q-046-06351.06–1.122010-11-142 × 18300.50–0.60
1.06–1.442010-11-135 × 1830, 1 × 7621
XMMXCS J030205.1−000003.6GS-2011B-Q-050-03321.172011-12-011 × 1830
1.182011-11-201 × 1098
1.17–1.452011-11-184 × 1830
GS-2011B-Q-050-04321.322011-12-311 × 1098
1.23–1.742011-12-304 × 18300.85–1.40
GS-2011B-Q-050-05331.27–1.572012-01-172 × 18300.7
XMMXCS J095940.7+023113.4GS-2010B-Q-046-02351.19–1.232011-01-093 × 1830
1.19–1.252011-01-084 × 1830
GS-2012A-Q-46-01351.19–1.292012-03-184 × 1830
1.19–1.232012-03-022 × 1830
GS-2012A-Q-46-02341.20–1.352012-03-273 × 1830
1.23–1.462012-03-233 × 18300.8
GS-2012A-Q-46-03341.21–1.542012-03-226 × 18300.65–0.70
XMMXCS J112349.3+052956.8GS-2012A-Q-46-05331.23–1.332012-04-225 × 1830
1.472012-04-211 × 1830
GS-2012A-Q-46-06321.25–1.652012-05-154 × 18300.63–0.76
1.45–1.662012-04-222 × 1830
GS-2010B-Q-046-03331.262011-01-311 × 1830
1.35–1.642011-01-292 × 1525, 1 × 975
1.23–1.242011-01-272 × 1830
XMMXCS J113602.9–032943.2GS-2012A-Q-46-07361.142012-05-241× 1830
1.122012-05-231 × 1830
1.12–1.162012-05-203 × 1830
1.122012-05-191 × 1830
GS-2012A-Q-46-08331.48–1.762012-07-152 × 1830
1.41–1.802012-07-113 × 18300.50–0.70
1.52012-07-101 × 1830
XMMXCS J134305.1−000056.8GS-2012A-Q-46-10361.16–1.232012-05-244 × 1830
1.242012-05-211 × 1830
GS-2012A-Q-46-11341.252012-07-101 × 1830
1.16–1.192012-07-092 × 1830
1.22012-07-061 × 1830
1.54–1.842012-06-222 × 1830
XMMXCS J145009.3+090428.8GN2012A-Q-070-05321.02–1.052012-07-092 × 18001.15
1.11–1.622012-06-264 × 18000.84–0.98
GN2012A-Q-070-06341.02–1.042012-07-072 × 1800
1.09–1.172012-07-062 × 1800
1.48–1.792012-06-272 × 1800
GN2012A-Q-070-07331.22–1.592012-07-223 × 1800
1.04–1.162012-07-083 × 18001
XMMXCS J215221.0–273022.6GS-2010B-Q-046-04361.14–1.242010-11-122 × 1830
1.02–1.212010-09-144 × 1830
GS-2011B-Q-050-06341.07–1.152011-10-052 × 1830
1.12–1.562011-09-184 × 18300.60–1.00
GS-2011B-Q-050-07341.00–1.102011-10-244 × 1830
1.05–1.122011-10-162 × 1830
XMMXCS J230247.7+084355.9GN2012A-Q-070-10341.372012-08-081 × 18000.60–0.68
1.02–1.112012-07-305 × 18000.43–0.86
GN2012A-Q-070-11331.18–1.312012-08-132 × 1800
1.02–1.082012-08-093 × 18000.60–0.68
1.192012-08-081 × 18001
Table A1.

Spectroscopic observations log. For all observations the R400 grating and the OG515 filter was used.

Cluster NameMaskSlitsAirmass rangeObservation dateFrames(s)Seeing (arcsec)
XMMXCS J005656.6−274031.9GS-2012B-Q-011-03331.222012-10-191 × 1830
1.48–1.762012-10-162 × 1830
1.01–1.042012-10-153 × 18300.76–0.80
GS-2012B-Q-011-04351.05–1.352012-11-144 × 1830
XMMXCS J015241.1−133855.9GS-2011B-Q-050-01331.05–1.212011-12-026 × 1830
GS-2011B-Q-050-02341.05–1.652011-12-036 × 1830
XMMXCS J021734.7−051326.9GS-2012B-Q-011-06341.14–1.482012-12-054 × 1830
XMMXCS J025006.4−310400.8GS-2012B-Q-011-09321.11–1.202012-11-242 × 1830
1.01–1.192012-11-214 × 1830
GS-2010B-Q-046-06351.06–1.122010-11-142 × 18300.50–0.60
1.06–1.442010-11-135 × 1830, 1 × 7621
XMMXCS J030205.1−000003.6GS-2011B-Q-050-03321.172011-12-011 × 1830
1.182011-11-201 × 1098
1.17–1.452011-11-184 × 1830
GS-2011B-Q-050-04321.322011-12-311 × 1098
1.23–1.742011-12-304 × 18300.85–1.40
GS-2011B-Q-050-05331.27–1.572012-01-172 × 18300.7
XMMXCS J095940.7+023113.4GS-2010B-Q-046-02351.19–1.232011-01-093 × 1830
1.19–1.252011-01-084 × 1830
GS-2012A-Q-46-01351.19–1.292012-03-184 × 1830
1.19–1.232012-03-022 × 1830
GS-2012A-Q-46-02341.20–1.352012-03-273 × 1830
1.23–1.462012-03-233 × 18300.8
GS-2012A-Q-46-03341.21–1.542012-03-226 × 18300.65–0.70
XMMXCS J112349.3+052956.8GS-2012A-Q-46-05331.23–1.332012-04-225 × 1830
1.472012-04-211 × 1830
GS-2012A-Q-46-06321.25–1.652012-05-154 × 18300.63–0.76
1.45–1.662012-04-222 × 1830
GS-2010B-Q-046-03331.262011-01-311 × 1830
1.35–1.642011-01-292 × 1525, 1 × 975
1.23–1.242011-01-272 × 1830
XMMXCS J113602.9–032943.2GS-2012A-Q-46-07361.142012-05-241× 1830
1.122012-05-231 × 1830
1.12–1.162012-05-203 × 1830
1.122012-05-191 × 1830
GS-2012A-Q-46-08331.48–1.762012-07-152 × 1830
1.41–1.802012-07-113 × 18300.50–0.70
1.52012-07-101 × 1830
XMMXCS J134305.1−000056.8GS-2012A-Q-46-10361.16–1.232012-05-244 × 1830
1.242012-05-211 × 1830
GS-2012A-Q-46-11341.252012-07-101 × 1830
1.16–1.192012-07-092 × 1830
1.22012-07-061 × 1830
1.54–1.842012-06-222 × 1830
XMMXCS J145009.3+090428.8GN2012A-Q-070-05321.02–1.052012-07-092 × 18001.15
1.11–1.622012-06-264 × 18000.84–0.98
GN2012A-Q-070-06341.02–1.042012-07-072 × 1800
1.09–1.172012-07-062 × 1800
1.48–1.792012-06-272 × 1800
GN2012A-Q-070-07331.22–1.592012-07-223 × 1800
1.04–1.162012-07-083 × 18001
XMMXCS J215221.0–273022.6GS-2010B-Q-046-04361.14–1.242010-11-122 × 1830
1.02–1.212010-09-144 × 1830
GS-2011B-Q-050-06341.07–1.152011-10-052 × 1830
1.12–1.562011-09-184 × 18300.60–1.00
GS-2011B-Q-050-07341.00–1.102011-10-244 × 1830
1.05–1.122011-10-162 × 1830
XMMXCS J230247.7+084355.9GN2012A-Q-070-10341.372012-08-081 × 18000.60–0.68
1.02–1.112012-07-305 × 18000.43–0.86
GN2012A-Q-070-11331.18–1.312012-08-132 × 1800
1.02–1.082012-08-093 × 18000.60–0.68
1.192012-08-081 × 18001
Cluster NameMaskSlitsAirmass rangeObservation dateFrames(s)Seeing (arcsec)
XMMXCS J005656.6−274031.9GS-2012B-Q-011-03331.222012-10-191 × 1830
1.48–1.762012-10-162 × 1830
1.01–1.042012-10-153 × 18300.76–0.80
GS-2012B-Q-011-04351.05–1.352012-11-144 × 1830
XMMXCS J015241.1−133855.9GS-2011B-Q-050-01331.05–1.212011-12-026 × 1830
GS-2011B-Q-050-02341.05–1.652011-12-036 × 1830
XMMXCS J021734.7−051326.9GS-2012B-Q-011-06341.14–1.482012-12-054 × 1830
XMMXCS J025006.4−310400.8GS-2012B-Q-011-09321.11–1.202012-11-242 × 1830
1.01–1.192012-11-214 × 1830
GS-2010B-Q-046-06351.06–1.122010-11-142 × 18300.50–0.60
1.06–1.442010-11-135 × 1830, 1 × 7621
XMMXCS J030205.1−000003.6GS-2011B-Q-050-03321.172011-12-011 × 1830
1.182011-11-201 × 1098
1.17–1.452011-11-184 × 1830
GS-2011B-Q-050-04321.322011-12-311 × 1098
1.23–1.742011-12-304 × 18300.85–1.40
GS-2011B-Q-050-05331.27–1.572012-01-172 × 18300.7
XMMXCS J095940.7+023113.4GS-2010B-Q-046-02351.19–1.232011-01-093 × 1830
1.19–1.252011-01-084 × 1830
GS-2012A-Q-46-01351.19–1.292012-03-184 × 1830
1.19–1.232012-03-022 × 1830
GS-2012A-Q-46-02341.20–1.352012-03-273 × 1830
1.23–1.462012-03-233 × 18300.8
GS-2012A-Q-46-03341.21–1.542012-03-226 × 18300.65–0.70
XMMXCS J112349.3+052956.8GS-2012A-Q-46-05331.23–1.332012-04-225 × 1830
1.472012-04-211 × 1830
GS-2012A-Q-46-06321.25–1.652012-05-154 × 18300.63–0.76
1.45–1.662012-04-222 × 1830
GS-2010B-Q-046-03331.262011-01-311 × 1830
1.35–1.642011-01-292 × 1525, 1 × 975
1.23–1.242011-01-272 × 1830
XMMXCS J113602.9–032943.2GS-2012A-Q-46-07361.142012-05-241× 1830
1.122012-05-231 × 1830
1.12–1.162012-05-203 × 1830
1.122012-05-191 × 1830
GS-2012A-Q-46-08331.48–1.762012-07-152 × 1830
1.41–1.802012-07-113 × 18300.50–0.70
1.52012-07-101 × 1830
XMMXCS J134305.1−000056.8GS-2012A-Q-46-10361.16–1.232012-05-244 × 1830
1.242012-05-211 × 1830
GS-2012A-Q-46-11341.252012-07-101 × 1830
1.16–1.192012-07-092 × 1830
1.22012-07-061 × 1830
1.54–1.842012-06-222 × 1830
XMMXCS J145009.3+090428.8GN2012A-Q-070-05321.02–1.052012-07-092 × 18001.15
1.11–1.622012-06-264 × 18000.84–0.98
GN2012A-Q-070-06341.02–1.042012-07-072 × 1800
1.09–1.172012-07-062 × 1800
1.48–1.792012-06-272 × 1800
GN2012A-Q-070-07331.22–1.592012-07-223 × 1800
1.04–1.162012-07-083 × 18001
XMMXCS J215221.0–273022.6GS-2010B-Q-046-04361.14–1.242010-11-122 × 1830
1.02–1.212010-09-144 × 1830
GS-2011B-Q-050-06341.07–1.152011-10-052 × 1830
1.12–1.562011-09-184 × 18300.60–1.00
GS-2011B-Q-050-07341.00–1.102011-10-244 × 1830
1.05–1.122011-10-162 × 1830
XMMXCS J230247.7+084355.9GN2012A-Q-070-10341.372012-08-081 × 18000.60–0.68
1.02–1.112012-07-305 × 18000.43–0.86
GN2012A-Q-070-11331.18–1.312012-08-132 × 1800
1.02–1.082012-08-093 × 18000.60–0.68
1.192012-08-081 × 18001

Supplementary data