EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 599 (March 2017) A&A, 599 (2017) A138 Full HTML
Free Access
Issue
A&A
Volume 599, March 2017
Article Number A138
Number of page(s) 15
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201628971
Published online 15 March 2017

© ESO, 2017

1. Introduction

Galaxy clusters represent the place where astrophysics and cosmology meet: while their overall internal dynamics is dominated by gravity, the astrophysical processes taking place on galactic scales leave observable imprints on the diffuse hot gas trapped within their potential wells (Giacconi et al. 2009). Galaxy clusters (e.g., Vikhlinin et al. 2009a, b; Mantz et al. 2010), in combination with supernovae (e.g., Riess et al. 2011), cosmic microwave background (e.g., Bennett et al. 2013; Hinshaw et al. 2013; Planck Collaboration XVI 2016), baryon acoustic oscillations (e.g., Dawson et al. 2013) and cosmological weak lensing (e.g., Schrabback et al. 2010; Laureijs et al. 2011; Marian et al. 2011; Heymans et al. 2013) can constrain the dark energy equation-of-state parameter, w, at both late and early cosmological epochs. Upcoming experiments such as eROSITA will improve the statistical power of cluster cosmology by a few orders of magnitude (e.g., Predehl et al. 2010; Merloni et al. 2012). However, there are two major challenges in upcoming surveys for high-precision cluster cosmology: accurate mass calibrations and efficient follow-up, such as optical identifications and redshift measurements.

The cluster cosmological applications are degenerate with the mass calibrations of galaxy clusters (e.g., Zhang & Wu 2003; Stanek et al. 2010; Merloni et al. 2012). Large X-ray surveys select galaxy clusters by their observables (e.g., Ebeling et al. 2000; Böhringer et al. 2004; Clerc et al. 2012; Hilton et al. 2012; Takey et al. 2013), particularly the X-ray luminosity (L), and not by their masses (M). The L−M relation is required to recover the selection function in terms of cluster masses and predict the cluster masses, hence the cluster mass function. The relations of bolometric X-ray luminosity versus mass calibrated with different samples differ significantly (e.g., Chen et al. 2007; Maughan 2007; Pacaud et al. 2007; Pratt et al. 2009; Mantz et al. 2010; Reichert et al. 2011), in which the slope varies from ~1.5 to ~2.0. The predicted numbers of detected clusters in the eROSITA survey based on different calibrations of the L−M scaling relations differ by up to a factor of two (e.g., Pillepich et al. 2012). In addition to the luminosity segregation between cool-core (CC) and non-cool-core (NCC) clusters, merging affects the properties of both the intracluster medium (ICM) and cluster galaxies (e.g., Poole et al. 2006; Evrard et al. 2008), but on different timescales (e.g., Roettiger et al. 1997). Mass estimates from optical spectroscopy and gravitational lensing, for instance, which are independent of the X-ray luminosity measurement, provide an X-ray blind reference of the cluster mass to calibrate the L−M relation (e.g., Kellogg et al. 1990; Wu et al. 1998, 1999; Zhang et al. 2008; Leauthaud et al. 2010; Israel et al. 2014, 2015; von der Linden et al. 2014; Sereno 2015; Sereno & Ettori 2015; Andreon et al. 2016). However, the dynamical mass estimates are sensitive to the cluster galaxy selection and member statistic and overestimate the mass of a merging cluster when the merger axis is along the line of sight (e.g., Biviano et al. 2006; Gifford & Miller 2013; Old et al. 2013; Saro et al. 2013; Wu & Huterer 2013; Rabitz et al. 2017).

The HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS, Reiprich & Böhringer 2002) is a flux-limited sample of 64 clusters selected from the ROSAT All-Sky Survey (RASS; Ebeling et al. 2000; Böhringer et al. 2004). We analyze high-quality XMM-Newton and ROSAT pointed observations as well as optical spectroscopic redshifts of 13 650 cluster galaxies for all 64 HIFLUGCS clusters. Excluding 2A 0335+096 with only three redshifts, we present the relation of X-ray luminosity versus dynamical mass for the 63 HIFLUGCS clusters and quantify the effect of mergers as well as CC systems on the L−M relation. A simultaneous mass calibration and cosmological application procedure (e.g., Bocquet et al. 2015; Mantz et al. 2015) can break the degeneracy between them. This approach is promising for the use of the combined X-ray and optical surveys in the near future because a vast number of telescopes are dedicated to carry out the optical spectroscopic surveys of galaxy clusters, such as eBOSS/SPIDERS1 (Clerc et al. 2016) and 4MOST2. This method, however, relies not only on independent measurements of at least two cluster quantities, for example, X-ray luminosity and dynamical mass, but also on well-understood knowledge of the uncertainties and potential biases in the two measured quantities (e.g., Applegate et al. 2014; von der Linden et al. 2014; Planck Collaboration XIII 2016).

Furthermore, efficient follow-up identifications of a large number of galaxy clusters become one of the main tasks in upcoming surveys for high-precision cluster cosmology (e.g., Merloni et al. 2012; Nandra et al. 2013; Pointecouteau et al. 2013). To optimize the optical spectroscopy follow-up, it is invaluable to investigate in which redshift range existing and upcoming multiwavelength surveys are suitable to identify groups and clusters of galaxies and to measure their redshifts. For further applications of the optical spectroscopic follow-up, it is worth testing in which redshift range these multiwavelength surveys can provide accurate dynamical masses that are sufficient for the L−M calibration.

In this paper, we simulate the optical spectroscopic follow-up of clusters by Monte Carlo (MC) resampling of the HIFLUGCS cluster galaxy redshifts at hand according to eight optical spectroscopic setups. We calibrate the redshift and dynamical mass estimates and quantify in which redshift range the tested optical spectroscopic setups are reliable for measuring the cluster redshifts and dynamical masses. We organize the paper as follows. In Sect. 2 we briefly describe the sample, data, and analysis method. The results on the L−M relation of the HIFLUGCS sample are presented in Sect. 3. The results on the redshift and dynamical mass calibrations using simulations of the re-sampled cluster galaxy redshifts are presented in Sect. 4. We summarize our conclusions in Sect. 5. Throughout the paper, we assume Ωm = 0.3, ΩΛ = 0.7, σ8 = 0.81, and H0 = 70 km s-1 Mpc-1. We quote a 68% confidence level. Unless explicitly stated otherwise, we apply the BCES bisector regression method (Akritas & Bershady 1996), taking into account measurement errors on both variables to determine the parameters and their errors in the fit.

Table 1

Cluster properties.

2. Sample, data, and analysis method

The HIFLUGCS is a sample of galaxy clusters with an X-ray flux limit (0.1–2.4 keV) of 2 × 10-11 erg/s/cm2 and Galactic latitude | b | ≥ 20.0 degrees covering two-thirds of the sky (Reiprich & Böhringer 2002; Hudson et al. 2010), selected from the ROSAT All-Sky Survey (RASS; Ebeling et al. 2000; Böhringer et al. 2004).

Table 2

Ratios of X-ray bolometric luminosity measurements in different annuli for the whole HIFLUGCS sample of 64 clusters.

2.1. X-ray data and analysis

We analyzed nearly 4 Ms of raw data from XMM-Newton for 63 clusters (Zhang et al. 2009, 2011, 2012). After cleaning and selecting the longest observation closest to the cluster center, we obtained ~ 1.3 Ms of XMM-Newton data for 59 clusters. We measured the X-ray properties for all 64 clusters in the HIFLUGCS by combining XMM-Newton and ROSAT data. The X-ray flux-weighted centroids are listed in Zhang et al. (2011, 2012) with the method described in Sect. 2.3 of Zhang et al. (2010).

2.1.1. X-ray luminosity

The X-ray luminosity was measured within the cluster radius, r500. This radius is defined as that within which the matter overdensity is 500 times the critical density of the Universe. It was derived from the dynamical mass in Sect. 2.2.2 centered on the X-ray flux-weighted centroid except for 2A 0335+096, which has three galaxy redshifts in total and of which the cluster radius r500 was derived from the scaling relation of mass versus gas mass. We note that Zhang et al. (2011) measured the X-ray luminosity within the r500 derived from the gas mass using the scaling relation of total mass versus gas mass derived in Pratt et al. (2009). To suppress the scatter of the L−M relation caused by cool cores, we used the core-corrected X-ray luminosity (Lcor hereafter), which is derived from the integration of the surface brightness assuming a constant value of the surface brightness within 0.2r500 equal to the value at 0.2r500, SX(R< 0.2r500) = SX(0.2r500), following Zhang et al. (2007). Here R is the projected cluster-centric distance. We note that this correction was only applied to determine the X-ray luminosity. We use the core-corrected bolometric X-ray luminosity (Table 1) for the relation of luminosity versus mass of the HIFLUGCS throughout the paper.

Some past work used a radius larger than r500 (e.g., Reiprich & Böhringer 2002) to derive the X-ray luminosity, and some used different methods to account for the cool cores (e.g., Pratt et al. 2009). Different methods in deriving the X-ray luminosity may lead to different normalization values of the L−M relation and also to different slopes. We therefore tabulate the ratios of the X-ray luminosity derived within different annuli for the whole HIFLUGCS of 64 clusters in Table 2. These values help to ensure a fair comparison between our results and those in other papers. Since the X-ray luminosity varies very little with the truncation radius, the correlation introduced by using the dynamical mass-determined r500 in deriving the X-ray luminosity to the Lcor−M relation is negligible.

2.2. Optical data and analysis

The positions of the brightest cluster galaxies (BCGs) for 63 clusters are listed in Table 1 in Zhang et al. (2011) and for RXC J1504.1−0248 in Zhang et al. (2012). We obtained the velocity of the cluster galaxies from the literature (updated until April 2013, including the compilation in Andernach et al. 2005; and Zhang et al. 2011, 2012). Since the individual error estimates are inhomogeneous, we decided not to weight the calculation of the average velocity. When there was more than one velocity per galaxy available, we calculated an average of the measurements, excluding discordant values and those with large errors.

We checked recent redshift surveys and found significant (> 90) numbers of new redshifts for Abell clusters: A2029 from Tyler et al. (2013), A2142 from Owers et al. (2011), A2255 from Tyler et al. (2014), and A2256 from Rines et al. (2016). For 2A 0335+96, NED now offers more redshifts from Huchra et al. (2012), resulting in 14 galaxies between 8400 and 13 200 km s-1, with a velocity dispersion of about 720 km s-1. For NGC 1550, we found 42 galaxies with a velocity dispersion of about 700 km s-1.

This work focuses on calibrating the uncertainties of the redshift and dynamical mass estimates when only reduced numbers of galaxy redshifts per cluster are used. The effect on our study by adding the new redshifts mentioned above is negligible. We therefore decided to stay with our compilation without updating it. In our upcoming study, we aim to follow up clusters with low numbers of spectroscopic members to ensure more than about 45 members for all clusters in the sample so that we can derive reliable mass estimates for all individual clusters and calibrate the relation of luminosity versus mass by taking the sample selection effect into account, that is, the Malmquist bias (Malmquist 1922) that arises from working with a flux-limited sample. We shall therefore revise the dynamical mass estimates, including the new redshifts mentioned above.

2.2.1. Cluster galaxy selection

In hierarchical structure-formation scenarios, the spherical infall model predicts a trumpet-shaped region in the diagram of line-of-sight velocity versus projected distance, the so-called caustic (e.g., Kaiser 1987). The boundary of the caustic defines galaxies inside as cluster galaxies and those outside as fore- and background galaxies. For each cluster, we plotted the line-of-sight velocity of the selected galaxies as a function of their projected distance from the BCG and located the caustic that efficiently excludes interlopers (e.g., Diaferio 1999; Katgert et al. 2004; Rines & Diaferio 2006). We considered only the galaxies inside the caustic as cluster galaxies and excluded the others from the subsequent analysis. We list the number of cluster galaxies (ngal) for the HIFLUGCS in Table 1. We gathered redshifts of a total of 13 650 cluster galaxies after the caustic member selection. Since 2A 0355+096 has only three redshifts, we excluded it from the study involving cluster redshift and dynamical mass estimates. The median number of spectroscopic members for the remaining sample is 188 per cluster.

There are another six clusters that have fewer than 45 cluster galaxies with spectroscopic redshifts in the HIFLUGCS (i.e., 13 galaxies for A0478, 22 for NGC1550, 42 for EXO 0422−086, 37 for Hydra A, 20 for S1101, and 44 for A2597), which were maintained in the redshift and dynamical mass study. All six systems with <45 cluster galaxy redshifts are CC clusters. The reason for having <45 cluster galaxy redshifts is that these clusters are not covered by large spectroscopic surveys such as SDSS or 2dF/6dF. We aim to observe them with ground-based telescopes. Recently, we obtained new redshifts from VLT/VIMOS for Abell S1101 (Rabitz et al. 2017).

2.2.2. Redshift, dynamical mass, and cluster radius

We applied the bi-weight estimator (e.g., Beers et al. 1990) to the cluster galaxies to measure the cluster redshift (z) and velocity dispersion (σ). The errors were estimated through 1000 bootstrap simulations.

The caustic method mentioned in Sect. 2.2.1 can provide an estimate of the escape velocity and thus the matter distribution for the halo according to the distribution of the caustic amplitude (e.g., Diaferio 1999). However, only with a large number of spectroscopic members, for instance, ~200 members, are we likely to recover a rather complete sampling of the caustic, which yields a reliable measurement of the underlying dark matter distribution of the cluster, as shown, e.g., by the simulations in Serra & Diaferio (2013), for instance. Their simulations also showed that fewer members tend to provide a slightly reduced amplitude of the caustic because the caustic is undersampled, which causes underestimation of the total mass.

To optimize the optical spectroscopic follow-up for galaxy clusters, particularly at high redshifts where the cluster galaxy member statistic is low, we focused on the method for deriving the dynamical mass based on a small number of cluster galaxies, which differs from the caustic method. There are a number of such methods that are calibrated by samples of simulated clusters (e.g., Biviano et al. 2006; Munari et al. 2013). We note that the difference between the dynamical masses derived with these methods is rather small, which is within the uncertainties of our dynamical mass estimates. In addition, we focused only on the relative variation in the dynamical mass measurements that is normalized to the input dynamical mass in the investigations of the systematic bias and uncertainties of redshift and dynamical mass estimates based on MC resampling (Sect. 4).

Our method of mass estimation is based entirely on the velocity dispersion, as explained in Sect. 3 of Biviano et al. (2006, Mσ therein). We followed the Navarro et al. (NFW, 1997) model used in Biviano et al. (2006) to derive M500 from Mσ. The dimensionless Hubble parameter, E(z)=Ωm(1+z)3+ΩΛ\hbox{$E(z)=\sqrt{\Omega_{\rm m}(1+z)^3 +\Omega_{\Lambda}}$}, was added to their formulation since not all our clusters are at z ~ 0. We used the dynamical mass at the radius within which the overdensity was 500 times the critical density as the cluster mass (M500, Table 1). We derived the cluster radius (r500) from the dynamical mass (M500), different from the method used by Zhang et al. (2011).

2.3. Cluster morphology

Although the cluster properties are driven by gravity, predominantly from dark matter, the baryon physics on small scales modifies the scaling relation from the self-similar prediction through cooling, merging, and feedback from star formation and AGN activities. We divided the HIFLUGCS into undisturbed versus disturbed subsamples as well as CC and NCC subsamples to study the effect of mergers and CC clusters, respectively, on our results.

thumbnail Fig. 1

Left and middle panels: histograms of the luminosity concentration (left) and the offset between the BCG and X-ray flux-weighted centroid (middle) of all 64 HIFLUGCS clusters with their Gaussian fits shown as solid curves. The mean of the Gaussian fits and their 1σ clipping are shown as dashed and solid vertical lines, respectively. Right panel: cL versus ΔR/r500 with its best fit, log 10cL = (0.593 ± 0.053) + (0.125 ± 0.028)log 10(ΔR/r500), as dash-dotted line. The vertical and horizontal solid lines denote the 1σ clipping below which the clusters are considered as undisturbed ones (open squares). The disturbed clusters are shown as open triangles.

2.3.1. Undisturbed versus disturbed clusters

Many studies found that a large offset between the BCG and X-ray peak/centroid probably indicates that more merging has occurred in the past (e.g., Katayama et al. 2003; Hudson et al. 2010; Zhang et al. 2010, 2011). We converted the angular separation between the X-ray flux-weighted center and BCG position into the physical separation at the cluster redshift for all 64 clusters in the HIFLUGCS. As shown in Zhang et al. (2011), the offset in projection between the X-ray flux-weighted center and BCG position scaled by the cluster radius, ΔR/r500, serves very well to separate the disturbed from the undisturbed clusters.

The luminosity concentration is also one of the parameters suitable to separate the clusters (e.g., Santos et al. 2008). Here we set this parameter to be cL = L(R ≤ r500) /L(R ≤ 0.2r500), and we integrated the surface brightness within a projected cluster-centric distance of r500 and 0.2r500 to derive the bolometric X-ray luminosity of L(R ≤ r500) and L(R ≤ 0.2r500), respectively. All HIFLUGCS clusters also have Chandra data, but mainly in the field of the cluster cores, not even out to 0.2r500 for a number of cases. We therefore preferred to rely on XMM-Newton data to compute the cL parameter.

In the following we combine the cuts of cL and ΔR/r500 to effectively separate the sample into undisturbed and disturbed clusters. As shown in Fig. 1 and Table 3, the histograms of both log 10cL and log 10(ΔR/r500) of the 64 clusters follow approximately Gaussian distributions. The best power-law fit between the two parameters for the 64 clusters is log 10cL = (0.593 ± 0.053) + (0.125 ± 0.028)log 10(ΔR/r500). We divided the whole sample into disturbed and undisturbed clusters using the 1σ-clipping of the Gaussian-mean values of these two histograms. Those with low values of both cL and ΔR/r500 are considered as undisturbed clusters (see the right panel of Fig. 1). The cuts divide the sample into 39 undisturbed clusters and 24 disturbed clusters.

Table 3

Gaussian fit parameters of the histograms of the logarithmic values of cL and ΔR/r500 for the full HIFLUGCS sample of all 64 clusters.

2.3.2. CC versus NCC clusters

The relation of gas mass versus total mass in recent work (e.g., Arnaud et al. 2005; Zhang et al. 2008) shows a slope shallower than unity, the self-similar prediction, which may be accounted for by the mass dependence of the gas mass fraction in CC and NCC systems (e.g., Eckert et al. 2012). Therefore, we also investigated our results for the subsamples of CC and NCC clusters, respectively.

The central cooling time can be accurately estimated from Chandra data because of its smaller point-spread function (PSF). We used the central cooling time calculated at 0.004r500 from Eq. (15) in Hudson et al. (2010) and divided the sample of the 64 clusters into 28 CC clusters and 36 NCC clusters as listed in Table 2 in Zhang et al. (2011).

thumbnail Fig. 2

a) Core-corrected bolometric X-ray luminosity versus dynamical mass for the 63 clusters with the best power-law fits for the 63 and 57 (ngal ≥ 45) clusters as black solid and dashed lines, respectively. b) Core-corrected X-ray bolometric luminosity versus dynamical mass for the 57 (ngal ≥ 45) clusters with the best power-law fits for the NCC, CC, disturbed, and undisturbed clusters as red dashed, blue dashed, black dot-dot-dot-dashed, and black dotted lines, respectively.

3. X-ray luminosity versus dynamical mass

We carried out the X-ray and optical analyses independently, in addition to taking the dynamical mass-determined r500 in computing the X-ray bolometric luminosity. Since the luminosity values measured within r500 and 2.5r500 differ only by ~15% on average, as shown in Table 2, using the cluster radius, r500, derived from the dynamical mass to compute the X-ray bolometric luminosity is not expected to cause any significant bias in our result.

We fit the Lcor−M relation by a power-law, log10􏽥L=Alog10􏽥M+B,\begin{equation} \log_{10}\widetilde{L}=A\log_{10}\widetilde{M}+B , \label{e:lm} \end{equation}(1)in which 􏽥L=LcorE(z)1044erg/s\hbox{$\widetilde{L}=\frac{L^{\rm cor}}{E(z)10^{44}{\rm erg/s}}$} and 􏽥M=M E(z)1014M⊙\hbox{$\widetilde{M}=\frac{M\,E(z)}{10^{14}M_{\odot}}$}. We note that the HIFLUGCS is a flux-limited sample, which consists of clusters covering a broad range of redshifts (i.e., from z = 0.0037 to z = 0.2153). Therefore, we included the E(z) factor in the scaling relation to account for the evolution of the geometry of the Universe in this redshift range. We assumed a constant standard deviation of the luminosity from this relation in logarithmic scale and used σlog 10L for the intrinsic scatter of the luminosity in logarithmic scale. We chose the relatively simple algorithm given in Akritas & Bershady (1996) to perform the fit to the relation y = Ax + B, in which χ2=∑i(yi−Axi−B)2σy,i2+A2σx,i2·\begin{equation} \chi^2= \sum_{i} \frac{(y_i - A x_i -B)^2}{\sigma^2_{y,i}+A^2 \sigma^2_{x,i}} \cdot \label{e:chi2} \end{equation}(2)The scatter of y consists of two components, the statistical and intrinsic scatter, which satisfy σy,i2=σy,i,sta2+σy,i,int2\hbox{$\sigma^2_{y,i}=\sigma^2_{y,i,{\rm sta}}+\sigma^2_{y,i, {\rm int}}$} (e.g., Weiner et al. 2006). In Table 4 we list the best-fit parameters and also include the normalization values by setting the slope to the self-similar prediction, namely 4/3. The data still agree with a slope of 4/3 within 2σ uncertainties of the best-fit slope parameters.

Table 4

Best power-law fits, log10􏽥L=Alog10􏽥M+B\hbox{$\log_{10}\widetilde{L}=A\log_{10}\widetilde{M}+B$}, in which 􏽥L=LcorE(z)1044erg/s\hbox{$\widetilde{L}=\frac{L^{\rm cor}}{E(z)10^{44}{\rm erg/s}}$} and 􏽥M=M E(z)1014M⊙\hbox{$\widetilde{M}=\frac{M\,E(z)}{10^{14}M_{\odot}}$}, of the Lcor−M relations.

3.1. Lcor – M relation for the ngal ≥ 45 clusters

The dynamical masses tend to be underestimated for clusters with too few redshift measurements according to the clusters in simulations (e.g., Biviano et al. 2006; Zhang et al. 2011). The effect becomes negligible (i.e., on a level of a few per cent) when clusters with a sufficient number of spectroscopic cluster galaxies are considered.

Zhang et al. (2011) found that the velocity dispersion measurements, and thus the dynamical masses, may be significantly underestimated for clusters with fewer than 45 spectroscopic member galaxies in the sample. To avoid any bias due to extreme outliers with too few cluster galaxies with spectroscopic redshifts, we show the Lcor−M relations for the 63 clusters (Fig. 2a) and for the 57 clusters that have at least 45 cluster galaxies with spectroscopic redshifts (Fig. 2b), respectively. For the 63 clusters, the slope, 1.29 ± 0.09, agrees with the self-similar prediction. However, most recent work shows steeper slopes, such as 1.77 ± 0.07 in Pratt et al. (2009). Taking a close look, the systems with very few redshifts appear to cause the shallow slope. In contrast to that of the 63 clusters, the slope of the 57 clusters is steeper. The selection effect, that is, the Malmquist bias that should be corrected for a flux-limited sample, should lead to an even steeper slope for the 57 clusters, which tends to agree better with recent X-ray calibrated L−M relations (e.g., Pratt et al. 2009). At the same time, the fact that the six discarded clusters are all CC clusters may also have an effect on the total best-fit slope, although there appears little difference between the slopes of the L−M relations between CC and NCC clusters.

thumbnail Fig. 3

X-ray core-to-r500 luminosity fraction versus dynamical mass for the CC (blue) and NCC (red) clusters. The curves show the regression fits for all (black), CC (blue) and NCC (red) clusters, respectively. There is no obvious dependence considering the observed scatter.

3.2. Lcor – M relations for the CC/NCC and undisturbed/disturbed subsamples

We list the best-fit Lcor−M relations for the different classes of clusters in Table 4. As shown in Figs. 2b and 3, the NCC clusters are undersampled in the low-mass regime since the HIFLUGCS is a flux-limited sample. The slope parameters agree within 1σ error between the undisturbed and disturbed clusters as well as between the CC and NCC clusters. There is no obvious mass dependence in the scatter of the core-to-total luminosity ratio (1 /cL) for the CC and NCC subsamples, respectively, as shown in Fig. 3. The Lcor−M relation for the disturbed clusters has smaller intrinsic scatter than that for the undisturbed clusters. The intrinsic scatter of the NCC clusters is also smaller than that of the CC clusters even when the core-corrected luminosity is used. This partially depends on how broad the cut is that was used to classify the corresponding cluster subsamples. The normalization of the disturbed (NCC) clusters agrees better with the X-ray calibrated relation (e.g., Zhang et al. 2008, see Table 3 therein; Pratt et al. 2009). Assuming a fixed slope of 4/3, the normalization agrees between different cluster types within the uncertainties.

A comparison of the best fits for the undisturbed (CC) clusters in the samples of the 57 and 63 clusters in Fig. 2 clearly shows that the scatter for the undisturbed (CC) clusters in the sample of 57 clusters with Ngal> 45 is lower than that for the sample of 63 clusters. As shown in Fig. 2b for the 57 clusters, the normalization on average agrees within 17% between the undisturbed and disturbed clusters, and within 32% between the CC and NCC clusters. The slopes still agree among the subsamples. As the scatter is dominated by the CC clusters, the scatter of the Lcor−M relation may mainly be due to the integrated effect of the structure formation history instead of recent extreme merging.

4. Systematic bias and uncertainties of redshift and dynamical mass estimates based on MC resampling

The results from our sample show that clusters with only few spectroscopic members tend to have underestimated dynamical masses. It is therefore of prime importance to quantify any systematic bias and uncertainties of the dynamical mass estimate and their implications on mass calibrations using upcoming optical spectroscopic surveys.

Both large optical spectroscopic surveys and individual pointed observations will be used to follow up galaxy clusters detected in upcoming large X-ray surveys such as eROSITA. The huge amount of optical spectroscopic follow-up data might be combined with the X-ray data for mass calibrations, in particular, in different redshift ranges. This allows constraining the evolution of the mass scaling relations, which becomes important when using high-redshift clusters for cosmology. In the following, we calculate the cluster redshift and estimate the dynamical mass using subsamples of member galaxies that we resampled according to eight optical spectroscopic setups, which demonstrates their ability for studying mass calibrations using the follow-up of the X-ray detected clusters. We note that the resampled galaxies were selected from the catalog of the caustic-bounded galaxies. Therefore, the probability of contamination from interlopers is low. The systematic uncertainties are mainly caused by statistical fluctuations in the resampling. Our study is also useful to assign the mass uncertainty in the stacking analysis using cluster samples with a rather clean member selection but a small number of members (e.g., Clerc et al. 2016).

4.1. Optical spectroscopic survey setups

The extended Baryon Oscillation Spectroscopic Survey (eBOSS; e.g., Schlegel et al. 2011) is a redshift survey covering a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000−4800. This survey targets objects up to redshift 2. With the eBOSS setup, the SPectroscopic IDentification of ERosita Sources (SPIDERS; e.g., Merloni et al. 2012) survey has been designed to follow up X-ray selected active galactic nuclei (AGN) and clusters over ~7500 degree2 area. Before eROSITA data will be available, SPIDERS is targeting ROSAT and XMM-Newton sources (see Clerc et al. 2016; Dwelly et al., in prep.). In the later years of operations, SPIDERS will target eROSITA sources detected in the first two years of observations (see Merloni et al. 2012). The 4-m Multi-Object Spectroscopic Telescope (4MOST, e.g., de Jong et al. 2012) is designed to obtain >1 million redshifts for >50 000 clusters in the southern sky in the eROSITA survey. The Euclid mission can follow up high-redshift galaxies through their emission lines (e.g., Laureijs et al. 2011; Pointecouteau et al. 2013). In this study, however, we do not discuss the optimization of the follow-up of the Euclid survey since Euclid provides slitless spectroscopy, which has no fiber constraint. Euclid will instead focus on how to improve the situation due to overlapping spectra in densely populated regions for optimizing its strategy, for example.

In practice, only several cluster galaxies per cluster can be targeted because of the fiber spacing constraint for multi-object spectroscopic surveys; this is not the case for slitless spectroscopic surveys like Euclid. Moreover, toward the high-redshift regime, a bright limiting magnitude cut can strongly limit the number of spectroscopic members. To separate the two effects, we performed the investigations in two steps: (i) for the setups using the closest fiber separation alone and (ii) for the setups including constraints from both the closest fiber separation and limiting magnitude cuts. SPIDERS uses a minimum fiber separation of 65′′. To demonstrate the effect of the fiber separation, we used a fiber spacing constraint of 55′′ and 65′′, respectively, in our MC simulations for SPIDERS. The limiting magnitude is the fiber magnitude. Toward the very high-redshift regime, optical spectroscopy of several of the bright galaxies already requires high-sensitivity optical spectrographs, for instance, on the Very Large Telescope (VLT). We also investigated the setups with a limited number of cluster redshifts. Below are the eight setups in which the BCG was always first selected for targeting.

  • (I)

    SPIDERS_55: it has a closest fiber separation of 55′′. Each position can be pointed only once.

  • (II)

    SPIDERS_65: it has a closest fiber separation of 65′′. Each position can be pointed only once.

  • (III)

    4MOST: it has a closest fiber separation of 20′′. Each position can be pointed four times.

  • (IV)

    SPIDERS_55m: it has a closest fiber separation of 55′′. Each position can be pointed only once. The limiting magnitude is iAB = 21.

  • (V)

    SPIDERS_65 m: it has a closest fiber separation of 65′′. Each position can be pointed only once. The limiting magnitude is iAB = 21.

  • (VI)

    4MOSTm: it has a closest fiber separation of 20′′. Each position can be pointed four times. The limiting magnitude is rAB = 22.

  • (VII)

    10zs: it samples the BCG and nine randomly selected remaining member galaxies per cluster.

  • (VIII)

    05zs: it samples the BCG and four randomly selected remaining member galaxies per cluster.

4.2. Input cluster redshifts, zin, in the I–VI setups

Follow-up of high-redshift systems is particularly important to constrain the evolution of mass calibrations, which is still poorly understood. Technically, nearby clusters can tolerate high values of the closest fiber separation in the optical spectroscopy. The angular size of a cluster decreases toward high redshift. The number of possible spectroscopic targets decreases with redshift given a fixed value of the closest fiber separation. Therefore, we tested the I–VI setups assuming input cluster redshifts at the HIFLUGCS cluster redshifts (z in Table 1) as well as placing the clusters at the redshift values of 0.2, 0.4, 0.6, and 0.8, respectively. We note that the relative rest-frame line-of-sight velocity of any member galaxy to its host cluster was not modified in this process.

4.3. Limiting magnitude cuts in the IV–VI setups

We included the limiting magnitude cuts in resampling the member galaxies in the IV–VI setups as follows.

4.3.1. Shape of the galaxy luminosity function (GLF)

We assumed that the observed spectroscopic member galaxies are those at the bright end of the GLF. We adopted the Schechter function fit of the GLF of the bright galaxies within a projected cluster-centric distance of 2 Mpc in the i and r bands (Popesso et al. 2005, Table 2, local background) for the SPIDERS and 4MOST surveys, respectively. Based on the results in Hansen et al. (2009), we assumed that the shape of the GLF of the bright red members does not change with cluster-centric radius for simplicity.

4.3.2. Normalization of the GLF

Several richness versus mass calibrations have been reported (e.g., Hansen et al. 2005, 2009; Johnston et al. 2007; Reyes et al. 2008). We took the richness versus weak-lensing mass calibration given by Eq. (10) in Hansen et al. (2009), in which the masses are measured within the radius where the mass overdensity is 200 relative to the critical density, and calculated the richness from the cluster dynamical mass for each HIFLUGCS cluster.

For simplicity, we assumed that the red member galaxies brighter than (i∗ + 1) are similar to those brighter than (r∗ + 1). Therefore, the number of the red member galaxies brighter than (i∗ + 1) is similar to the number of the red members brighter than (r∗ + 1). We obtained the richness by integrating the GLF down to (i∗ + 1) magnitude for the SPIDERS survey and (r∗ + 1) magnitude for the 4MOST survey, respectively. Assuming the richness obtained from the dynamical mass equal to the richness derived by integrating the GLF, we obtained the normalization of the GLF.

4.3.3. Galactic extinction and K-corrections

We corrected apparent magnitudes for Galactic extinction using the maps of Schlegel et al. (1998). We assumed that photometric errors at bright magnitudes are 0.05. We used typical colors and their color uncertainties of luminous red galaxies of Maraston et al. (2009), and applied K-correction and evolutionary correction using the luminous red galaxy template in kcorrect (v4_2, Blanton & Roweis 2007).

4.3.4. Member selection

The maximum number of member galaxies, nmax, that can be spectroscopically followed-up was computed by integrating the GLF down to the limiting magnitude. In Sect. 4.3.1 we assumed that the observed spectroscopic member galaxies are those filling in the bright end of the GLF. Therefore, if the number of cluster members for the HIFLUGCS clusters (ngal in Table 1) is larger than nmax, then the fraction of the ngal spectroscopic members that are brighter than the limiting magnitude is fobs = nmax/ngal; otherwise, all existing spectroscopic members are brighter than the limiting magnitude, nmax = ngal.

thumbnail Fig. 4

Bias (solid lines) and dispersion (dashed lines) of the redshift (left panels) and dynamical mass (right panels) measurements in the SPIDERS_55 (upper panels), SPIDERS_65 (middle panels) and 4MOST (lower panels) setups from the input cluster values as a function of the resampled number of redshifts in use. The curves are smoothed with a boxcar average of the specified width of 11 to avoid spikes due to undersampling for a few clusters.

thumbnail Fig. 5

Bias (solid curves) and dispersion (dashed curves) of the redshift (left panel) and dynamical mass (right panel) measurements in the 10zs setup from the input cluster values versus input number of cluster galaxies. The curves are smoothed with a boxcar average of the specified width of 11 to avoid spikes due to under-sampling for a few clusters. Note that ngal is the input number of galaxies and the output number of galaxies is always 10.

Many cluster members in the observational sample have no i and r magnitudes from the NED. As a preliminary selection, we first took the BCG and then randomly took fobsngal−1 additional spectroscopic members using the MC method. The initially selected members were filtered further according to the values of the closest fiber separation.

4.3.5. Effect of flux loss of the aperture magnitude

The angular sizes of very bright cluster galaxies in the nearby Universe (i.e., z< 0.4) are usually larger than the fiber aperture. The flux-loss correction may become important in deriving the apparent magnitudes from the measured flux in the fiber aperture when considering the limiting magnitude cut. Additionally, the PSF affects the measured apparent magnitudes due to convolution of galaxy images caused by the seeing. As shown by our model of the flux loss in targeting galaxies with spectroscopic fibers in Appendix A, this effect is negligible and skipped in our further analyses.

4.4. Redshift and dynamical mass measurements

To estimate the average of the systematic uncertainties, we carried out 500 resampling runs per cluster per setup per input cluster redshift (zin = z,0.2,0.4,0.6,0.8 for I–VI setups and zin = z for VII–VIII setups). In each run, we measured the cluster redshift, velocity dispersion, and dynamical mass following the procedure in Sect. 2.2.2 when the number of resampled members reached ten, ngal,i ≥ 10. Otherwise (ngal,i< 10), we computed the cluster redshift, which is the mean of the member galaxy redshifts, but not the velocity dispersion and dynamical mass.

The Gaussian mean and bi-weight (see Sect. 2.2.2) mean values of the histograms of the redshift and dynamical mass estimates from 500 runs based on the resampled members agree well in all cases except for the 10zs setup. The bi-weight mean of the histograms is in general closer to the input value than the Gaussian mean for the 10zs setup. Therefore, we used the bi-weight mean to calculate the results. The number of member redshifts used per cluster per setup per input cluster redshift, ngal,MC, was computed as the bi-weight mean of ngal,i with i = 0,...,499. As ngal,i< 10 may be obtained for some resampling runs, ngal,MC can be lower than ten for some cases involving the dynamical mass estimates, in which the dynamical mass per cluster per setup per input cluster redshift, M500,MC, is computed as the bi-weight mean of M500,j with j = i if ngal,i ≥ 10. Simulated clusters with fewer than ten members were excluded from the mass computation.

4.5. Summary of results

4.5.1. Fraction of selected members

When we used the collisional distance constraint alone (setups I−III), almost all members were selected in the 4MOST setup, and more than 95% of the members were selected in the SPIDERS_55 and SPIDERS_65 setups. In the case of ten or five redshifts per cluster, only a small fraction of the members were selected. The fraction of selected members including the limiting magnitude cut (setups IV–VI) decreases faster toward the high-redshift end than that in the resampling assuming only the collisional distance constraint.

As shown in the right panels of Figs. 4–6, the scatter in the dynamical mass estimates for the clusters with at least ten members is within 50%, which is the minimum required precision for a stacking analysis for the mass calibration given the known intrinsic scatter in the mass-observable relation (e.g., Biviano et al. 2006). We therefore inspected the fraction of the clusters with at least ten members. All of the 63 clusters in the SPIDERS and 4MOST cases that assume no magnitude cut have more than ten galaxy members per cluster. In the setups with the magnitude limit instead, the fraction of the clusters with at least ten members per cluster decreases with redshift as shown in Fig. 7. In the redshift bins of 0.2, 0.4, 0.6, and 0.8 in the SPIDERS_55m setup, there are 1, 3, 8, and 18 clusters out of 63 clusters that have fewer than ten redshifts per cluster. The case for the SPIDERS_65m setup, which uses a larger fiber separation (65′′), is slightly worse. In the redshift bins of 0.4, 0.6 and 0.8 using the 4MOSTm setup, there are 1, 2, and 4 clusters out of 63 clusters that have fewer than ten redshifts per cluster.

thumbnail Fig. 6

Bias (solid curves) and dispersion (dashed curves) of the redshift (left panels) and dynamical mass (right panels) measurements in the SPIDERS_55m (upper panels), SPIDERS_65m (middle panels) and 4MOSTm (lower panels) setups from the input cluster values as a function of the re-sampled number of redshifts in use. The curves are smoothed with a boxcar average of the specified width of 11 to avoid spikes due to under-sampling for a few clusters.

4.5.2. Redshift estimates

We considered all clusters with ngal,MC> 1 in the redshift estimates. Here we did not consider ngal,MC = 1 clusters because the BCG is always taken according to the selection method (Sect. 4.3.4), which leads to zero dispersion based on 500 MC realizations. The cluster redshifts are well recovered to a level of a few per cent in all cases. There is no obvious systematic dependence on the cluster redshift. The left panels of Figs. 4–6 and Fig. 8 show the bias of the redshifts obtained in the resampled procedure from the input cluster redshifts and its dispersion based on the 500 MC realizations. We note that the catalog of the galaxy redshifts of the HIFLUGCS is a rather clean member galaxy input catalog, probably almost free of interlopers. Uncertainties derived from the scatter in the down-sampling therefore do not account for the effect of interlopers. This artificially restricts the uncertainties of the redshift estimates to less than 0.01. In Figs. 5 and 8 the scatter increases slightly with the number of the input galaxies, which may indicate that the scatter of redshift with a fixed number of spectroscopic members depends on the richness. Richness-dependent values of minimum numbers of galaxy redshifts for different richness populations may thus be required for a stacking analysis to ensure an equivalent scatter introduced by individual systems in the redshift estimates.

4.5.3. Dynamical mass estimates

There is no redshift or mass dependence in the systematic uncertainties of the dynamical mass estimates. The right panels of Figs. 4–6 show the bias of the dynamical mass estimates obtained in the resampled procedure from the cluster total mass and its dispersion based on the 500 MC realizations. Here we did not consider ngal,i< 10 realizations, as noted in Sect. 4.4. Therefore, toward the low-ngal,MC end, the bias and dispersion is slightly underestimated because not all 500 MC realizations can be used.

The dynamical mass is well recovered with a bias lower than ~20% on average. The dispersion in Figs. 4–6 demonstrates the precision that the corresponding setup can reach as a function of redshift. This information can be used for the sample selection corresponding to the required precision for the purpose of mass calibration in upcoming surveys. Given the same required accuracy, the 4MOST survey can provide mass calibration up to a higher redshift than the SPIDERS survey. With ten redshifts per cluster, the dynamical mass measurement can easily be biased by a factor of two for individual clusters, which would not be suitable for a reliable mass calibration unless the clusteres were stacked. Similar to that shown in the redshift estimates, the scatter again slightly increases with the number of the input galaxies, this may indicate a dependence on richness. Richness-dependent cuts of minimum numbers of galaxy redshifts in different richness bins are also be required for a stacking analysis to ensure an equivalent scatter introduced by individual systems in the dynamical mass estimates.

4.5.4. Limitations

The current paper describes an ideal situation of using an input member catalog almost free of contamination from interlopers for the spectroscopic follow-up; this does not represent the reality in general. According to real surveys, we need to construct our input galaxy candidate sample for the spectroscopic follow-up, which cannot avoid contamination from interlopers, for example. Additionally, we need to select a number of certain member galaxies from a large catalog for the follow-up to determine the redshift and velocity dispersion, which requires knowledge of the selection effects for recovering the mass estimates in a reliable way. The input member galaxy sample for the simulation in the current paper, however, probably contains no interlopers. Neither can we quantify the selection effect because the sample is a collection of cluster galaxies from the available literature and not a spectroscopic follow-up of a well-defined photometric galaxy member sample, for instance. A further study based on mock data in simulations with known selection functions as well as sufficient fore- and background galaxies in the galaxy sample will help to address these questions. Another limitation is that clusters may evolve from z = 0.8 to the present epoch, an effect that is not fully accounted for by simply shifting the HIFLUGCS clusters to higher redshifts. It is therefore not obvious that the HIFLUGCS sample is the best test-case for forecasting the 4MOST and SPIDERS surveys at higher redshifts.

thumbnail Fig. 7

Fraction of the clusters with at least ten redshifts per cluster after the resampling using the SPIDERS_55m (black triangles), SPIDERS_65m (red diamonds, with 0.03 offset) and 4MOSTm (blue squares, with 0.06 offset) setups as a function of the assigned input cluster redshift (zin) in the re-sampling. zin = 0 case refers to the resampled clusters at the cluster redshifts z.

5. Conclusions

We analyzed ~ 1.3 Ms of clean XMM-Newton data and ROSAT pointed observations as well as optical spectroscopic redshifts of 13 650 cluster member galaxies for all 64 HIFLUGCS clusters. Excluding 2A 0335+096 with only three redshifts, we presented the Lcor−M relation for 63 nearby clusters of galaxies in the HIFLUGCS. For the optimal use of the optical spectroscopic surveys for high-redshift galaxy clusters and groups observed in upcoming X-ray surveys for mass calibrations, we carried out MC resampling of the galaxy members with spectroscopic redshifts and calibrated the systematic uncertainties in the redshift and dynamical mass estimates. We predicted the redshift and dynamical mass estimates assuming the SPIDERS and 4MOST optical spectroscopic survey setups, respectively, using MC resampling of the 63 HIFLUGCS clusters by placing the sample at the actual cluster redshifts as well as at redshifts of 0.2, 0.4, 0.6, and 0.8. Aiming for high-redshift cluster or group systems, we also predicted the redshift estimates based on five and ten spectroscopic members per cluster, respectively, and the mass estimates based on ten spectroscopic members per cluster. We list the conclusions in detail as follows.

  • Given sufficient numbers (i.e., ≥45) of member galaxies in computing the dynamical masses, the Lcor−M relations agree between the disturbed and undisturbed clusters.

  • The CC clusters still dominate the scatter in the Lcor−M relation even when the core-corrected X-ray luminosity is used. This indicates that the scatter of the Lcor−M scaling relation mainly reflects the structure formation history of the clusters.

  • The dynamical masses can be measured within 10% uncertainty using the SPIDERS_55, SPIDERS_65 and 4MOST setups, which are independent of cluster redshift and mass. With ten or more redshifts per cluster, the dynamical masses can be recovered with a bias lower than 20% on average, in which the dynamical masses are underestimated for most systems.

  • The bias of the cluster dynamical mass estimates increases toward the high-redshift end. The underestimation of the cluster masses on average using the SPIDERS_65m setup, in which both the collisional distance of 65′′ and magnitude cut are considered, is better than 19%, 28%, 34%, and 37% in the redshift bins of 0.2, 0.4, 0.6, and 0.8. The dynamical mass is recovered with an underestimation lower than 20% up to redshift 0.6 using the 4MOSTm setup, in which both the collisional distance and magnitude cut are considered. In the redshift bin of 0.8, the underestimation of the dynamical masses on average, according to the 4MOSTm setup, is still lower than 24%. Assuming the SPIDERS_55m (4MOSTm) setup, the dynamical masses can be used as an independent reference blind to the X-ray observables to calibrate the cluster mass with an underestimation lower than 20% up to redshift 0.2 (0.6) with 2% (3%) catastrophic outliers (i.e., fewer than ten members per cluster) in upcoming X-ray surveys. Empirically, this bias can be corrected using a complete sample from observations or a mock sample in simulations, which includes the causes of the bias, such as contamination from fore- and background galaxies.

thumbnail Fig. 8

Bias (solid curves) and dispersion (dashed curves) of the redshift measurements in the 05zs setup from the input cluster values versus input number of cluster galaxies. The curves are smoothed with a boxcar average of the specified width of 11 to avoid spikes due to undersampling for a few clusters. ngal is the input number of galaxies, and the output number of galaxies is always 5.

Acknowledgments

The XMM-Newton project is an ESA Science Mission with instruments and contributions directly funded by ESA Member States and the USA (NASA). The XMM-Newton project is supported by the Bundesministerium für Wirtschaft und Technologie/Deutsches Zentrum für Luft- und Raumfahrt (BMWi/DLR, FKZ 50 OX 0001) and the Max-Planck Society. 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. We thank Tom Dwelly for helpful discussion and an anonymous referee for the insight and expertise that improved the work. Y.Y.Z. acknowledges support by the German BMWi through the Verbundforschung under grant 50 OR 1506. T.H.R. acknowledges support from the DFG through the Heisenberg research grant RE 1462/5 and grant RE 1462/6.

References

  1. Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 [NASA ADS] [CrossRef] [Google Scholar]
  2. Andernach, H., Tago, E., Einasto, M., Einasto, J., & Jaaniste, J. 2005, Nearby Large-Scale Structures and the Zone of Avoidance, eds. A. P. Fairall, & P. Woudt, ASP Conf. Ser., 329, 283 [Google Scholar]
  3. Andreon, S., Serra, A. L., Moretti, A., & Trinchieri, G. 2016, A&A, 585, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Applegate, D. E., von der Linden, A., Kelly, P. L., et al. 2014, MNRAS, 439, 48 [NASA ADS] [CrossRef] [Google Scholar]
  5. Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 893 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Beers, T. C., Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bennett, C. L., Larson, D., Weiland, J. L., et al. 2013, ApJS, 208, 20 [Google Scholar]
  8. Bernardi, M., Hyde, J. B., Sheth, R. K., Miller, C. J., & Nichol, R. C. 2007, AJ, 133, 1741 [NASA ADS] [CrossRef] [Google Scholar]
  9. Biviano, A., Murante, G., Borgani, S., et al. 2006, A&A, 456, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bocquet, S., Saro, A., Mohr, J. J., et al. 2015, ApJ, 799, 214 [NASA ADS] [CrossRef] [Google Scholar]
  12. Böhringer, H., Schuecker, P., Guzzo, L., et al. 2004, A&A, 425, 367 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Chen, Y., Reiprich, T. H., Böhringer, H., Ikebe, Y., & Zhang, Y.-Y. 2007, A&A, 466, 805 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Clerc, N., Sadibekova, T., Pierre, M., et al. 2012, MNRAS, 423, 3561 [NASA ADS] [CrossRef] [Google Scholar]
  15. Clerc, N., Merloni, A., Zhang, Y.-Y., et al. 2016, MNRAS, 463, 4490 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dawson, K. S., Schlegel, D. J., Ahn, C. P., et al. 2013, AJ, 145, 10 [Google Scholar]
  17. de Jong, R. S., Bellido-Tirado, O., Chiappini, C., et al., 2012, Proc. SPIE, 8446, 84460T [Google Scholar]
  18. Diaferio, A. 1999, MNRAS, 309, 610 [NASA ADS] [CrossRef] [Google Scholar]
  19. Ebeling, H., Edge, A.C., Allen, S. W., et al. 2000, MNRAS, 365, 1021 [Google Scholar]
  20. Eckert, D., Vazza, F., Ettori, S., et al. 2012, A&A, 541, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Evrard, A. E., Bialek, J., Busha, M., et al. 2008, ApJ, 672, 122 [NASA ADS] [CrossRef] [Google Scholar]
  22. Giacconi, R., Borgani, S., Rosati, P., et al. 2009, Astro2010: The Astronomy and Astrophysics Decadal Survey. Science white Papers, 90 [Google Scholar]
  23. Gifford, D., & Miller, C. J. 2013, ApJ, 768, L32 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hansen, S. M., McKay, T. A., Wechsler, R. H., et al. 2005, ApJ, 633, 122 [NASA ADS] [CrossRef] [Google Scholar]
  25. Hansen, S. M., Sheldon, E. S., Wechsler, R. H., & Koester, B. P. 2009, ApJ, 699, 1333 [NASA ADS] [CrossRef] [Google Scholar]
  26. Heymans, C., Grocutt, E., Heavens, A., et al. 2013, MNRAS, 432, 2433 [NASA ADS] [CrossRef] [Google Scholar]
  27. Hilton, M., Romer, A. K., Kay, S. T., et al. 2012, MNRAS, 424, 2086 [NASA ADS] [CrossRef] [Google Scholar]
  28. Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [NASA ADS] [CrossRef] [Google Scholar]
  29. Huchra, J. P., Macri, L. M., Masters, K. L., et al. 2012, ApJS, 199, 26 [NASA ADS] [CrossRef] [Google Scholar]
  30. Hudson, D. S., Mittal, R., Reiprich, T. H., et al. 2010, A&A, 513, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Hyde, J. B., & Bernardi, M. 2009, MNRAS, 394, 1978 [NASA ADS] [CrossRef] [Google Scholar]
  32. Israel, H., Reiprich, T. H., Erben, T., et al. 2014, A&A, 564, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Israel, H., Schellenberger, G., Nevalainen, J., Massey, R., & Reiprich, T. H. 2015, MNRAS, 448, 814 [NASA ADS] [CrossRef] [Google Scholar]
  34. Johnston, D. E., Sheldon, E. S., Wechsler, R. H., et al. 2007, ArXiv e-prints [arXiv:0709.1159] [Google Scholar]
  35. Kaiser, N. 1987, MNRAS, 227, 1 [NASA ADS] [CrossRef] [Google Scholar]
  36. Katayama, H., Hayashida, K., Takahara, F., & Fujita, Y. 2003, ApJ, 585, 687 [NASA ADS] [CrossRef] [Google Scholar]
  37. Katgert, P., Biviano, A., & Mazure, A. 2004, ApJ, 600, 657 [NASA ADS] [CrossRef] [Google Scholar]
  38. Kellogg, E., Falco, E., Forman, W., Jones, C., & Slane, P. 1990, in Gravitational Lensing, ASP Conf. Ser., 360, 141 [Google Scholar]
  39. Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, ArXiv e-prints [arXiv:1110.3193] [Google Scholar]
  40. Leauthaud, A., Finoguenov, A., Kneib, J.-P., et al. 2010, ApJ, 709, 97 [NASA ADS] [CrossRef] [Google Scholar]
  41. Malmquist, K. G. 1922, Meddelanden fran Lunds Astronomika Observatorium, Ser. I, No. 100, [Google Scholar]
  42. Mantz, A., Allen, S. W., Rapetti, D., & Ebeling, H. 2010, MNRAS, 406, 1759 [NASA ADS] [Google Scholar]
  43. Mantz, A. B., von der Linden, A., Allen, S. W., et al. 2015, MNRAS, 446, 2205 [NASA ADS] [CrossRef] [Google Scholar]
  44. Maraston, C., Strömbäck, G., Thomas, D., Wake, D. A., & Nichol, R. C. 2009, MNRAS, 394, L107 [NASA ADS] [CrossRef] [Google Scholar]
  45. Marian, L., Hilbert, S., Smith, R. E., Schneider, P., & Desjacques, V. 2011, ApJ, 728, L13 [NASA ADS] [CrossRef] [Google Scholar]
  46. Maughan, B. J. 2007, ApJ, 668, 772 [NASA ADS] [CrossRef] [Google Scholar]
  47. Merloni, A., Predehl, P., Becker, W., et al. 2012, ArXiv e-prints [arXiv:1209.3114] [Google Scholar]
  48. Munari, E., Biviano, A., Borgani, S., Murante, G., & Fabjan, D. 2013, MNRAS, 430, 2638 [NASA ADS] [CrossRef] [Google Scholar]
  49. Nandra, K., Barret, D., Barcons, X., et al. 2013, ArXiv e-prints [arXiv:1306.2307] [Google Scholar]
  50. Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
  51. Old, L., Gray, M. E., & Pearce, F. R. 2013, MNRAS, 434, 2606 [NASA ADS] [CrossRef] [Google Scholar]
  52. Owers, M. S., Nulsen, P. E. J., & Couch, W. J. 2011, ApJ, 741, 122 [NASA ADS] [CrossRef] [Google Scholar]
  53. Pacaud, F., Pierre, M., Adami, C., et al. 2007, MNRAS, 382, 1289 [NASA ADS] [CrossRef] [Google Scholar]
  54. Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 124, 266 [NASA ADS] [CrossRef] [Google Scholar]
  55. Pillepich, A., Porciani, C., & Reiprich, T. H. 2012, MNRAS, 422, 44 [NASA ADS] [CrossRef] [Google Scholar]
  56. Planck Collaboration XVI. 2016, 594, A16 [Google Scholar]
  57. Planck Collaboration XIII. 2016, 594, A13 [Google Scholar]
  58. Pointecouteau, E., Reiprich, T. H., Adami, C., et al. 2013, ArXiv e-prints [arXiv:1306.2319] [Google Scholar]
  59. Poole, G. B., Fardal, M. A., Babul, A., et al. 2006, MNRAS, 373, 881 [NASA ADS] [CrossRef] [Google Scholar]
  60. Popesso, P., Böhringer, H., Romaniello, M., & Voges, W. 2005, A&A, 433, 415 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Pratt, G. W., Croston, J. H., Arnaud, M., & Böhringer, H., 2009, A&A, 498, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Predehl, P., Böhringer, H., Brunner, H., et al. 2010, Proc. SPIE, 7732 [Google Scholar]
  63. Rabitz, A., Zhang, Y.-Y., Schwope, A., et al. 2017, A&A, 597, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Reiprich, T. H., & Böhringer, H. 2002, ApJ, 567, 716 [NASA ADS] [CrossRef] [Google Scholar]
  65. Reichert, A., Böhringer, H., Fassbender, R., & Mühlegger, M. 2011, A&A, 535, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Reyes, R., Mandelbaum, R., Hirata, C., Bahcall, N., & Seljak, U. 2008, MNRAS, 390, 1157 [NASA ADS] [CrossRef] [Google Scholar]
  67. Riess, A. G., Macri, L., Casertano, S., et al. 2011, ApJ, 730, 119 [NASA ADS] [CrossRef] [Google Scholar]
  68. Rines, K., & Diaferio, A. 2006, ApJ, 132, 1275 [Google Scholar]
  69. Rines, K. J., Geller, M. J., Diaferio, A., & Hwang, H. S. 2016, ApJ, 819, 63 [NASA ADS] [CrossRef] [Google Scholar]
  70. Roettiger, K., Loken, C., & Burns, J. O. 1997, ApJS, 109, 307 [NASA ADS] [CrossRef] [Google Scholar]
  71. Santos, J. S., Rosati, P., Tozzi, P., et al. 2008, A&A, 483, 35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Saro, A., Mohr, J. J., Bazin, G., & Dolag, K. 2013, ApJ, 772, 47 [CrossRef] [Google Scholar]
  73. Sereno, M. 2015, MNRAS, 450, 3665 [NASA ADS] [CrossRef] [Google Scholar]
  74. Sereno, M., & Ettori, S. 2015, MNRAS, 450, 3675 [NASA ADS] [CrossRef] [Google Scholar]
  75. Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [NASA ADS] [CrossRef] [Google Scholar]
  76. Schlegel, D., Abdalla, F., Abraham, T., et al. 2011, ArXiv e-prints [arXiv:1106.1706] [Google Scholar]
  77. Schrabback, T., Hartlap, J., Joachimi, B., et al. 2010, A&A, 516, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Serra, A. L., & Diaferio, A. 2013, ApJ, 768, 116 [NASA ADS] [CrossRef] [Google Scholar]
  79. Stanek, R., Evrard, A. E., Böhringer, H., Schuecker, P., & Nord, B. 2006, ApJ, 648, 956 [NASA ADS] [CrossRef] [Google Scholar]
  80. Stanek, R., Rasia, E., Evrard, A. E., Pearce, F., & Gazzola, L. 2010, ApJ, 715, 1508 [NASA ADS] [CrossRef] [Google Scholar]
  81. Takey, A., Schwope, A., & Lamer, G. 2013, A&A, 558, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Tyler, K. D., Rieke, G. H., & Bai, L. 2013, ApJ, 773, 86 [NASA ADS] [CrossRef] [Google Scholar]
  83. Tyler, K. D., Bai, L., & Rieke, G. H. 2014, ApJ, 794, 31 [NASA ADS] [CrossRef] [Google Scholar]
  84. Vikhlinin, A., Burenin, R. A., Ebeling, H., et al. 2009a, ApJ, 692, 1033 [NASA ADS] [CrossRef] [Google Scholar]
  85. Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009b, ApJ, 692, 1060 [NASA ADS] [CrossRef] [Google Scholar]
  86. von der Linden, A., Mantz, A., Allen, S. W., et al. 2014, MNRAS, 443, 1973 [NASA ADS] [CrossRef] [Google Scholar]
  87. Weiner, B. J., Willmer, C. N. A., Faber, S. M., et al. 2006, ApJ, 653, 1049 [NASA ADS] [CrossRef] [Google Scholar]
  88. Wu, H.-Y., & Huterer, D. 2013, MNRAS, 434, 2556 [NASA ADS] [CrossRef] [Google Scholar]
  89. Wu, X.-P., Fang, L.-Z., & Xu, W. 1998, A&A, 338, 813 [NASA ADS] [Google Scholar]
  90. Wu, X.-P., Xue, Y.-J., & Fang, L.-Z. 1999, ApJ, 524, 22 [NASA ADS] [CrossRef] [Google Scholar]
  91. Zhang, Y.-Y., & Wu, X.-P., 2003, ApJ, 583, 529 [NASA ADS] [CrossRef] [Google Scholar]
  92. Zhang, Y.-Y., Finoguenov, A., Böhringer, H., et al. 2007, A&A, 467, 437 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Zhang, Y.-Y., Finoguenov, A., Böhringer, H., et al. 2008, A&A, 482, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Zhang, Y.-Y., Reiprich, T. H., Finoguenov, A., Hudson, D. S., & Sarazin, C. L. 2009, ApJ, 699, 1178 [NASA ADS] [CrossRef] [Google Scholar]
  95. Zhang, Y.-Y., Okabe, N., Finoguenov, A., et al. 2010, ApJ, 711, 1033 [NASA ADS] [CrossRef] [Google Scholar]
  96. Zhang, Y.-Y., Andernach, H., Caretta, C. A., et al. 2011, A&A, 526, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  97. Zhang, Y.-Y., Verdugo, M., Klein, M., & Schneider, P. 2012, A&A, 542, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Effect of flux loss of the aperture magnitude

Appendix A.1: Fraction of recovered flux

We used a Sérsic profile for the two-dimensional (2D) surface brightness profiles of galaxies (e.g., Peng et al. 2002), Σ(x,y)∝exp[−κrellRe)(1/n−1]\hbox{$ \Sigma(x,y) \propto \exp \left[-\kappa \left(\frac{r_{\rm ell}}{R_{\rm e}} \right)^{1/n}-1 \right]$} with rell=[|x|(c+2)+|yq|(c+2)]1/(c+2)\hbox{$ r_{\rm ell} = \left[\left\vert x \right\vert ^{(c+2)}+\left\vert \frac{y}{q} \right\vert ^{(c+2)} \right]^{1/(c+2)}$}. The parameter q is the ellipticity, and we used its distribution observed in SDSS (Fig. 3 in Hyde & Bernardi 2009, peaked at ~0.8). For a de Vaucouleurs profile, c = 0, n = 4, κ = 7.67 for n = 4, and Σ(x,y)∝exp[−7.67rell(/Re)1/4]\hbox{$\Sigma(x,y) \propto \exp \left[ -7.67 \left(r_{\rm ell}/R_{\rm e}\right)^{1/4}\right]$}. Following Bernardi et al. (2007), we computed the effective radius in units of kpc from the r-band absolute magnitude for a passive, elliptical galaxy, ⟨ log Re | Mr ⟩ = −0.681(Mr + 21)/2.5 + 0.343, using the spectral energy distribution template in Maraston et al. (2009).

We assumed \hbox{$1\farcs4$} seeing for SPIDERS and 1″ seeing for 4MOST (private communication with T. Dwelly), in which the seeing is better than the given value for 90% of the time. The diameter of the aperture is ?= 2″ for SPIDERS and \hbox{$1\farcs5$} for 4MOST. The resulting flux was computed by integrating the surface brightness profile convolved with the seeing within the aperture. The fractions of recovered flux fsynth in the fibers of the SPIDERS and 4MOST surveys are shown in Fig. A.1. Our results are consistent with that of the objects classified as “LRG” or “GALAXY”, which have spectroscopic redshifts from the BOSS in SDSS DR11.

Appendix A.2: Effect on our sampling

For galaxies with an intrinsic magnitude only slightly brighter than the limiting magnitude cut, its aperture magnitude may become fainter than the limiting magnitude when considering flux loss. We inspected this effect on our sampling by comparing the aperture magnitudes of the member galaxies used in the resampling with the limiting magnitude cut. Although there is quite a dramatic difference (~1–1.5 mag) between the total magnitude and the 2″ aperture magnitude, we found that the aperture magnitudes of all member galaxies are well above the limiting magnitude cut for the clusters at zin< 0.4 and zin = 0.8. This is because the galaxies are much brighter than the limiting magnitude at the low redshift of zin = 0.2, and the angular sizes of the galaxies are well within the aperture size at the high redshift of zin = 0.8 such that flux loss plays no role. At redshifts of zin = 0.4 and 0.6, the aperture magnitude considering flux loss leads to a number of member galaxies fainter than the limiting magnitude in a couple of clusters among 63 clusters in the resampling. No more than three out of 63 clusters are affected by the flux loss. We therefore consider this effect negligible. Flux loss plays a main role in degrading the spectral quality, which is not part of this study.

thumbnail Fig. A.1

Fractions of recovered flux fsynth in the fibers of the SPIDERS (left panel) and 4MOST (right panel) configurations. The x-axis and y-axis are the galaxy redshift and r-band absolute magnitude, respectively.

All Tables

Table 1

Cluster properties.

In the text
Table 2

Ratios of X-ray bolometric luminosity measurements in different annuli for the whole HIFLUGCS sample of 64 clusters.

In the text
Table 3

Gaussian fit parameters of the histograms of the logarithmic values of cL and ΔR/r500 for the full HIFLUGCS sample of all 64 clusters.

In the text
Table 4

Best power-law fits, log10􏽥L=Alog10􏽥M+B\hbox{$\log_{10}\widetilde{L}=A\log_{10}\widetilde{M}+B$}, in which 􏽥L=LcorE(z)1044erg/s\hbox{$\widetilde{L}=\frac{L^{\rm cor}}{E(z)10^{44}{\rm erg/s}}$} and 􏽥M=M E(z)1014M⊙\hbox{$\widetilde{M}=\frac{M\,E(z)}{10^{14}M_{\odot}}$}, of the Lcor−M relations.

In the text

All Figures

thumbnail Fig. 1

Left and middle panels: histograms of the luminosity concentration (left) and the offset between the BCG and X-ray flux-weighted centroid (middle) of all 64 HIFLUGCS clusters with their Gaussian fits shown as solid curves. The mean of the Gaussian fits and their 1σ clipping are shown as dashed and solid vertical lines, respectively. Right panel: cL versus ΔR/r500 with its best fit, log 10cL = (0.593 ± 0.053) + (0.125 ± 0.028)log 10(ΔR/r500), as dash-dotted line. The vertical and horizontal solid lines denote the 1σ clipping below which the clusters are considered as undisturbed ones (open squares). The disturbed clusters are shown as open triangles.
In the text
thumbnail Fig. 2

a) Core-corrected bolometric X-ray luminosity versus dynamical mass for the 63 clusters with the best power-law fits for the 63 and 57 (ngal ≥ 45) clusters as black solid and dashed lines, respectively. b) Core-corrected X-ray bolometric luminosity versus dynamical mass for the 57 (ngal ≥ 45) clusters with the best power-law fits for the NCC, CC, disturbed, and undisturbed clusters as red dashed, blue dashed, black dot-dot-dot-dashed, and black dotted lines, respectively.
In the text
thumbnail Fig. 3

X-ray core-to-r500 luminosity fraction versus dynamical mass for the CC (blue) and NCC (red) clusters. The curves show the regression fits for all (black), CC (blue) and NCC (red) clusters, respectively. There is no obvious dependence considering the observed scatter.
In the text
thumbnail Fig. 4

Bias (solid lines) and dispersion (dashed lines) of the redshift (left panels) and dynamical mass (right panels) measurements in the SPIDERS_55 (upper panels), SPIDERS_65 (middle panels) and 4MOST (lower panels) setups from the input cluster values as a function of the resampled number of redshifts in use. The curves are smoothed with a boxcar average of the specified width of 11 to avoid spikes due to undersampling for a few clusters.
In the text
thumbnail Fig. 5

Bias (solid curves) and dispersion (dashed curves) of the redshift (left panel) and dynamical mass (right panel) measurements in the 10zs setup from the input cluster values versus input number of cluster galaxies. The curves are smoothed with a boxcar average of the specified width of 11 to avoid spikes due to under-sampling for a few clusters. Note that ngal is the input number of galaxies and the output number of galaxies is always 10.

In the text
thumbnail Fig. 6

Bias (solid curves) and dispersion (dashed curves) of the redshift (left panels) and dynamical mass (right panels) measurements in the SPIDERS_55m (upper panels), SPIDERS_65m (middle panels) and 4MOSTm (lower panels) setups from the input cluster values as a function of the re-sampled number of redshifts in use. The curves are smoothed with a boxcar average of the specified width of 11 to avoid spikes due to under-sampling for a few clusters.

In the text
thumbnail Fig. 7

Fraction of the clusters with at least ten redshifts per cluster after the resampling using the SPIDERS_55m (black triangles), SPIDERS_65m (red diamonds, with 0.03 offset) and 4MOSTm (blue squares, with 0.06 offset) setups as a function of the assigned input cluster redshift (zin) in the re-sampling. zin = 0 case refers to the resampled clusters at the cluster redshifts z.
In the text
thumbnail Fig. 8

Bias (solid curves) and dispersion (dashed curves) of the redshift measurements in the 05zs setup from the input cluster values versus input number of cluster galaxies. The curves are smoothed with a boxcar average of the specified width of 11 to avoid spikes due to undersampling for a few clusters. ngal is the input number of galaxies, and the output number of galaxies is always 5.

In the text
thumbnail Fig. A.1

Fractions of recovered flux fsynth in the fibers of the SPIDERS (left panel) and 4MOST (right panel) configurations. The x-axis and y-axis are the galaxy redshift and r-band absolute magnitude, respectively.
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.