SEARCHING FOR BULK MOTIONS IN THE INTRACLUSTER MEDIUM OF MASSIVE, MERGING CLUSTERS WITH CHANDRA CCD DATA

, , , and

Published 2016 April 5 © 2016. The American Astronomical Society. All rights reserved.
, , Citation Ang Liu et al 2016 ApJ 821 29 DOI 10.3847/0004-637X/821/1/29

0004-637X/821/1/29

ABSTRACT

We search for bulk motions in the intracluster medium (ICM) of massive clusters showing evidence of an ongoing or recent major merger with spatially resolved spectroscopy in Chandra CCD data. We identify a sample of six merging clusters with >150 ks Chandra exposure in the redshift range 0.1 < z < 0.3. By performing X-ray spectral analysis of projected ICM regions selected according to their surface brightness, we obtain the projected redshift maps for all of these clusters. After performing a robust analysis of the statistical and systematic uncertainties in the measured X-ray redshift zX, we check whether or not the global zX distribution differs from that expected when the ICM is at rest. We find evidence of significant bulk motions at more than 3σ in A2142 and A115, and less than 2σ in A2034 and A520. Focusing on single regions, we identify significant localized velocity differences in all of the merger clusters. We also perform the same analysis on two relaxed clusters with no signatures of recent mergers, finding no signs of bulk motions, as expected. Our results indicate that deep Chandra CCD data enable us to identify the presence of bulk motions at the level of vBM > 1000 km s−1 in the ICM of massive merging clusters at 0.1 < z < 0.3. Although the CCD spectral resolution is not sufficient for a detailed analysis of the ICM dynamics, Chandra CCD data constitute a key diagnostic tool complementing X-ray bolometers on board future X-ray missions.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Clusters of galaxies are the largest virialized systems in the universe, and the strong interplay between their baryonic and dark components is such that they are at the crossroads between astrophysics and cosmology. It is widely accepted that clusters form in a hierarchical manner from a primordial density perturbation field, and evolve through the accretion and merging of other virialized halos. In particular, merging between halos frequently happens at early epochs and still occurs today. Such merging processes strongly affect the dynamic structure of clusters, often producing bulk motions in the hot intracluster medium (hereafter ICM). On one hand, ICM bulk motions could potentially provide an important diagnostic of the formation and evolution of clusters. On the other hand, the ICM velocity field is extremely difficult to measure with our present facilities, and, more importantly, since measurements of the cluster properties rely on the assumption of hydrostatic equilibrium, unnoticed bulk motions may significantly affect observables such as total mass proxies and baryon mass fractions (Nagai et al. 2013; Nelson et al. 2014). Consequently, this will also affect studies of the formation of large-scale structure and constraints on cosmological parameters (see Allen et al. 2011 for a review).

Bremsstrahlung emission from the optically thin, hot ICM is observed in the X-ray band, and it constitutes a powerful diagnostics of the thermodynamics of the ICM itself, allowing one to measure the temperature and density of the emitting gas. In addition, thanks to the almost perfect collisional equilibrium, the measurement of the line emission from highly ionized heavy elements constitutes a direct probe of the metal content of the ICM. Merging processes can introduce several clear signatures in the ICM, such as the irregular distribution of the X-ray surface brightness, the inhomogeneous distribution of temperatures, and fluctuations in the ICM velocity field, as already mentioned. These signatures last for one or more dynamical times until the relaxation processes bring the ICM back to hydrostatic equilibrium. The surface brightness distribution and temperature map of the ICM of massive clusters is routinely obtained from X-ray imaging and spatially resolved spectroscopy. The direct observation of bulk motion instead presents several difficulties. In some rare cases, the bulk motion velocities can be estimated indirectly via the properties of the bow shock resulting from a tangential merging process, as in the well-known cases of the Bullet cluster (Markevitch et al. 2002), A2146 (Russell et al. 2011), A520 (Markevitch et al. 2005), and a few others. The turbulence expected as a result of the stirring of the ICM during the merger process can be investigated by Doppler line broadening in the high-resolution spectra obtained with the Reflection Grating Spectrometer on board XMM-Newton. However, only upper limits in the range 200–600 km s−1 have been obtained, and only for the central regions (Sanders et al. 2010; Pinto et al. 2015). Alternatively, the power spectrum of the velocities can be inferred from the power spectrum of the gas density fluctuations (Zhuravleva et al. 2014).

A straightforward approach is to directly measure the ICM velocity field along the line of sight through the Doppler shift of its emission lines. The most prominent emission line of the ICM X-ray spectra is the ubiquitous K-shell line complex of H-like and He-like iron in the 6.7–6.9 keV rest frame, originally detected by Mitchell et al. (1976). Many other emission lines, including the L-shell from iron and α-element transitions, may also be found, particularly in the soft energy range (0.5–2.0 keV). However, the Kα iron complex is the only one which is detectable in any cluster observed with a total number of net counts as low as ∼1000 in the entire X-ray band (see Yu et al. 2011). In fact, iron line emission has been identified in X-ray clusters as distant as z ∼ 1.5 (Rosati et al. 2004, 2009; Stanford et al. 2005; Tozzi et al. 2013), although the detection may be challenging at z > 1.5 (see, e.g., Tozzi et al. 2015).

Despite this, applying this method to the measurement of ICM bulk motions is extremely challenging. The typical gas velocity of 1000 km s−1 (Nagai et al. 2003) expected in major mergers corresponds only to a ∼20 eV shift of the Kα iron line complex in the 6.7–6.9 keV rest frame. Such a shift is below the CCD energy resolution of XMM-Newton/EPIC, Chandra/ACIS, and Suzaku/XIS (∼100 eV). This constitutes a major obstacle to the direct measurement of the line-of-sight ICM velocity field. Despite these difficulties, there are a few works that have successfully achieved detections of bulk motions in massive, nearby clusters using spatially resolved spectral analysis of CCD data. The first successful detection was made in the Centaurus cluster with ASCA (Dupke & Bregman 2001a). The typical bulk motion velocity was constrained to (2.4 ± 0.1) × 103 km s−1, and was eventually confirmed by Chandra data (Dupke & Bregman 2006). Another case is the Perseus cluster, which was also found to present ICM bulk motions by ASCA (Dupke & Bregman 2001b). The first systematic search for ICM bulk motions was performed by Dupke & Bregman (2005) with ASCA, but only 2 out of 12 low-redshift (z < 0.13) clusters showed reliable signatures of bulk motions. In addition, a significant velocity difference of (5.9 ± 1.6) × 103 km s−1 between two regions of Abell 576 was found by combining Chandra and XMM-Newton data (Dupke et al. 2007).

Thanks to its lower background at high energies, Suzaku/XIS has been used to search for radial bulk motions in the ICM through X-ray redshift measurements. However, only upper limits of the order of 2000 km s−1 are derived for AWM7, Abell 2319, and Coma (Sato et al. 2008, 2011; Sugawara et al. 2009, respectively), and ∼300 km s−1 for Centaurus (Tamura et al. 2014). Only in one case (Abell 2256; Tamura et al. 2011) has a robust detection of a bulk motion of 1500 ± 300 km s−1 been found. Recently, Ota & Yoshida (2015) searched for gas bulk motion in eight nearby clusters with Suzaku, finding signs of large bulk velocity in excess of their uncertainties in only two of them (Abell 2029 and Abell 2255). Despite the low background of Suzaku/XIS, its low angular resolution is a major limiting factor for this investigation. Surprisingly, Chandra ACIS has not been used often for these studies, despite its high angular resolution which allows one to focus on smaller regions and, in principle, to detect more efficiently the effects of bulk motions.

Recently, we proposed a strategy to search for bulk motion in the ICM along the line of sight by performing spatially resolved spectral analysis in Chandra CCD data of the so-called Bullet cluster 1E0657-56 (Liu et al. 2015). This case turned out to be particularly challenging due to the extremely high temperatures, which make the iron emission line complex less prominent, and the fact that the merger is occurring in the plane of the sky. Nevertheless, we were able to find tantalizing ∼2σ evidence for a bulk motion of the significant mass of ICM in two regions near the center. This is consistent with the picture where the ICM is pushed perpendicularly to the bullet trail due to the extremely large pressure reached when crossing the main cluster center. Future missions carrying X-ray bolometers (e.g., Astro-H; Takahashi et al. 2014) may be able to confirm or rule out this result (see the recent simulations by Biffi et al. 2013).

In this work, we intend to apply our technique to a sample of massive, merging clusters imaged by Chandra. Our goal is to investigate whether ICM bulk motions in medium-redshift (z ∼ 0.1–0.3) clusters classified as major mergers can be detected in Chandra ACIS data. We also try to provide robust estimators of the ICM bulk motions whenever possible, such as the global rms velocity difference across the cluster vBM and the maximum velocity difference measured across the cluster Δvmax. We also aim to provide a qualitative description of the dynamical status of the ICM in single clusters through visual inspection of the redshift map, in order to exploit the angular resolution of Chandra.

The paper is organized as follows. In Section 2, we describe the sample selection. In Section 3, we describe the X-ray data reduction and analysis, including a brief description of the strategy already presented in Liu et al. (2015). In Section 4, we present results for the single clusters, while in Section 5 we synthesize the results for the entire sample. In Section 6, we discuss future prospects for improvements in this field. Finally, in Section 7, we summarize our conclusions. Throughout this paper, we adopt the 7 years WMAP cosmology, with Ωm = 0.272, ΩΛ = 0.728, and H0 = 70.4 km s−1 Mpc−1 (Komatsu et al. 2011). Quoted error bars always correspond to a 1σ confidence level.

2. SAMPLE SELECTION

We search the Chandra data archive for clusters classified as mergers based on their X-ray morphology, radio diffuse emission, and/or galaxy dynamics, in the redshift range 0.1 < z < 0.3. Since our goal is to perform spatially resolved spectral analyses with high signal-to-noise ratios (S/Ns), we require long exposure times. We fix the threshold for the minimum total exposure to 150 ks on ACIS-I or ACIS-S. The signatures of an occurring, or impending, major merger, are a clear bimodal or multiple stucture in the redshift distribution of the member galaxies, a strongly disturbed appearance in the X-ray surface brightness, or the presence of a radio halo or radio relic. We found six clusters showing at least one of these prominent properties. The clusters are listed in Table 1 in order of increasing optical redshift, together with Chandra ObsID and total exposure time (after data reduction).

Table 1.  Sample of Merging Clusters and Corresponding Chandra Observations Used in This Work

Cluster Name Optical Redshift ObsID Total Exposure (ks)
A2142 0.090(1) ACIS-S: 15186-16564-16565 155.1
A2034 0.113(2) ACIS-I: 7695-2204-12885 254.5
    13192-12886-13193  
A115 0.197(3) ACIS-I: 13458-13459-15578-15581 310.6
A520 0.203(4) ACIS-I: 4215-9424-9425-9426-9430 516.4
1RXS J0603.3+4214 0.225(5) ACIS-I:15171-15172-15323 235.9
A2146 0.234(6) ACIS-I: 12245-12246-12247-13020 375.3
    13021-13023-13120-13138  

Note. The observing mode is VFAINT for all of the ObsID. The exposure times in the fourth column correspond to the effective values after data reduction. Optical redshift references: (1) Parekh et al. (2015), (2) Muriel & Coenda (2014), (3) Hiroi et al. (2013), (4) Cassano et al. (2013), (5) Van Weeren et al. (2012), (6) Piffaretti et al. (2011).

Download table as:  ASCIITypeset image

A2142 is classified as a merger on the basis of its X-ray morphology (Parekh et al. 2015) and the presence of a radio halo (Cuciti et al. 2015). A2034 is classified as a merger on the basis of a shock visible in the X-ray surface brightness (Owers et al. 2014). Abell 115 has a clear bimodal appearance in the X-ray and has a radio relic (Cuciti et al. 2015). A520 shows a clear radio halo (Cuciti et al. 2015) and a clearly disturbed appearance in the X-ray band. 1RXS J0603.3+4214 has a clear bimodal structure and a bright, peculiar radio relic, known as the "Toothbrush" (Van Weeren et al. 2012). A2146 shows a prominent bow shock in the Chandra X-ray image, indicating a violent merger analogous to the more famous Bullet cluster 1E0657-56, which was the subject of our first investigation of ICM bulk motion (Liu et al. 2015). The properties of each cluster will be discussed in more details in Section 4.

As an immediate check of the robustness of our analysis, we also select two relaxed clusters in the same redshift range and with comparable exposure depth. The two clusters in this control sample are Abell 1689 and Abell 1835, and are listed in Table 2. Clearly, we do not expect to detect any departure from statistical noise in the zX distribution of these relaxed clusters.

Table 2.  Relaxed Clusters and Corresponding Chandra Observations Used in This Work

Cluster Name Optical Redshift ObsID Total Exposure (ks)
A1689 0.183(7) ACIS-I: 5004-6930-7289-7701 151.3
A1835 0.252(8) ACIS-I: 6880-6881-7370 193.7

Note. The observing mode is VFAINT for all of the ObsID. The exposure times in the fourth column correspond to the effective values after data reduction. Optical redshift reference: (7) Hiroi et al. (2013), (8) Girardi et al. (2014).

Download table as:  ASCIITypeset image

3. DATA REDUCTION AND ANALYSIS

3.1. Data Reduction

The Chandra observations used in this work are listed in Tables 1 and 2. All of the observations were taken in VFAINT mode with ACIS-I, except in the case of A2142 which was observed with ACIS-S. The data are reduced using the ciao software (version 4.7) with CALDB 4.6.8. The detailed data reduction procedure has been described in Liu et al. (2015). The area of the clusters used in this analysis is defined as the circular region which provides the maximum S/N in the 0.5–10 keV band once the unresolved sources are removed. Point sources are detected by running the wavdetect task and are then visually inspected, in particular, to remove those deeply embedded in the strong ICM emission. The background files are extracted from a series of circular regions as far as possible from the cluster area, but still on the solid angle defined by the overlap of all of the ObsIDs. To define the cluster subregions, we apply the contour binning technique (Sanders 2006) within the maximum S/N circle. Regions are identified simply on the basis of the surface brightness contours in the 0.5–10 keV band image. To make sure that each region can be used to constrain the X-ray redshift with comparable accuracy, we set a criterion in terms of S/N in the 0.5–10 keV band in each region. Clearly, this criterion does not take into account the expected gradient in the iron abundance, which significantly affects the visibility of the iron emission line complex. We set this threshold to S/N > 140, except in the case of Abell 2146 where we choose a weaker constraint of S/N > 100 due to the relatively lower number of counts. The S/N threshold and the number of regions selected in each cluster are listed in Table 3. For each region, the spectrum is extracted from the merged event file after the removal of unresolved sources, while the calibration files (response matrix files, RMF, and ancillary response files, ARF) are generated independently for each ObsID and then combined by weighting them by the corresponding exposure times. In this way, we keep track of all the differences in the ACIS effective area among the different regions in the detector and among the ObsIDs taken at different epochs.

Table 3.  The Minimum S/N Required in Each Region and the Number of Regions Selected for Spectral Analysis in Each Cluster

Cluster Name Minimum S/N Number of Regions
  (0.5–10 keV)  
A2142 140 52
A2034 140 14
A115 140 18
A520 140 20
1RXS J0603.3+4214 140 9
A2146 100 19
A1689 140 10
A1835 140 11

Download table as:  ASCIITypeset image

3.2. Spectral Analysis

In this work, we use the same analysis strategy presented in Liu et al. (2015) to measure the projected X-ray redshift zX in each region and estimate the total (statistical and systematic) uncertainty. Our analysis is optimized on the basis of the spectral simulations shown in Liu et al. (2015). Here, we simply recall the relevant aspects of the adopted spectral analysis.

The spectra are analyzed using Xspec v12.8.2 (Arnaud 1996). To successfully model the X-ray emission, we use two mekal plasma emission models (Mewe et al. 1985, 1986; Liedahl et al. 1995) that include thermal bremsstrahlung and line emission, with abundances measured relative to the solar values of Asplund et al. (2005) in which Fe/H = 3.6 × 10−5. During the fit, the redshift parameters of the two thermal components are always linked together, while the two temperatures and abundances are left as free parameters.

The main advantage in using a double-temperature thermal spectrum is reducing the possible bias in the measurement of the iron line centroid due to the presence of unnoticed thermal structure along the line of sight. As a matter of fact, different temperature components above 3 keV are not detectable, and a single-temperature thermal component would return an excellent fit even in the presence of widely different temperatures (see Mazzotta et al. 2004). However, in Liu et al. (2015), using spectral simulations we showed that the use of two components provides better results for the measurement of the actual redshift of the ICM, while the temperature structure, which is not our main interest here, is essentially left unconstrained. Clearly, this also introduces a larger systematic uncertainty associated with the temperature structure, which is carefully evaluated during our analysis, as explained below.

Galactic absorption is described by the model tbabs (Wilms et al. 2000). The central values of the Galactic HI column density NHGal are based on Kalberla et al. (2005). We conservatively allow ${\mathrm{NH}}_{\mathrm{Gal}}$ to vary by an interval of ∼10% around the central value. This has a small effects since our best-fit zX are obtained from spectral fits in the hard (2–10 keV) band, which is only marginally affected by Galactic absorption. Cash statistics (Cash 1979) are applied to the unbinned source and background spectra in order to exploit the full spectral resolution of the ACIS-I and ACIS-S CCD. Cash statistics are preferred for faint spectra with respect to the canonical χ2 analysis of binned data (Nousek & Shue 1989). To avoid local minima, we repeat the fit several times before and after running the steppar command on all of the free parameters, particularly on the redshift, with a step of δz = 10−4. Finally, the plots of ΔCstat versus redshift are visually inspected to investigate whether or not there are other possible minima around the best-fit values (see also Yu et al. 2011), which could indicate a noisy spectrum and therefore an unreliable best-fit value for zX. We decide to conservatively exclude all of the results with a secondary minimum closer than ΔC = 6.6, which corresponds to a formal 99% confidence level for one free parameter.

Our reference spectral analysis consists of two steps. As the first step, we fit the spectra in the 2–10 keV energy range to obtain the best-fit redshift.We choose to use only the 2–10 keV energy range to minimize the effects of the presence of low-temperature components on the centroid of the iron line emission complex. This conservative choice is also supported by the spectral simulations we presented in Liu et al. (2015). For each region, we produce a plot showing the ΔCstat value versus the redshift, obtained by varying the redshift parameter and marginalizing the fit with respect to the other parameters. The minimum is generally well defined and roughly symmetric, and therefore provides a robust estimate of the statistical uncertainty.

As a second step, we carefully take into account the systematic uncertainties on redshift due to the unknown thermal structure of the ICM. We conservatively consider all of the possible temperature distributions of the ICM by setting a four-dimensional grid of spectral parameters, namely, temperature and abundance of the two thermal components T1, T2, Z1, and Z2. The steps of the grid are typically 0.5 keV for the temperature and 0.05 for the metal abundance. Clearly, only a subgrid of temperature and abundance values should be taken into account, since some combination of them provides an unacceptable fit. To efficiently select the subgrid, we make use of the full 0.5–10 keV energy range in order to exploit the complete information on the thermal structure and the chemical composition. Therefore, for each region, we measure the best-fit redshift in the full, 0.5–10 keV band leaving all of the parameters free, and collect the absolute minimum Cmin. Then, we select all of the sets of values on the grid that provide a best fit close enough to the absolute best fit. The criterion we adopt is given by ${\rm{\Delta }}{C}_{{\rm{stat}}}\equiv {C}_{{\rm{stat}}}-{C}_{{\rm{min}}}\lt 4.72$, which corresponds to a 1σ confidence level for four free parameters. By applying the criterion ΔCstat < 4.72 for the fits performed in the full 0.5–10 keV band, we select a subgrid of T1, T2, Z1, and Z2 values statistically compatible within 1σ with the spectrum observed in each region. Then, we run the spectral fit on the subgrid, but focusing again in the 2–10 keV band only, and measure a set of best-fit zX values corresponding to the subgrid parameters. The distribution of zX obtained in this manner is used to derive a redshift range defined by the upper and lower 90% percentiles of the zX distribution on the grid, and thus the systematic uncertainties on redshift. Since the systematic uncertainty σsyst is due to the degeneracy of the temperature values, we assume that it is uncorrelated to the statistical error σstat. Therefore, we compute the total 1σ uncertainty as ${\sigma }_{{\rm{tot}}}=\sqrt{{\sigma }_{{\rm{stat}}}^{2}+{\sigma }_{{\rm{syst}}}^{2}}$. We do not consider other relevant sources of systematic errors. Systematics related to calibration issues or the time dependence of the gain are kept under control with a direct check on the fluorescent lines of Au Lα and Ni Kα, which are prominent in the ACIS-I background spectrum at 7.5 and 9.7 keV, respectively. Typically, uncertainties introduced by calibration variation (as a function of position on the CCD or as a function of the observation period) on the merged spectrum of each region correspond to errors ∼5 times smaller than the typical total uncertainty on zX (see Liu et al. 2015). Having assessed that gain-calibration does not significantly affect our conclusions, we do not perform a full treatment of the gain-calibration uncertainties in ACIS-I data.

The complete spectral results for each clusters are presented, where the best-fit redshifts and the corresponding lower and upper statistic and systematic error bars are listed for each region of each cluster (values are rounded to the fourth decimal digit). Results are presented in the form of redshift maps and significance maps. In particular, in the significance maps, we combine the information included in the best-fit redshift and the corresponding total error by showing $({z}_{{\rm{X}}}-\langle {z}_{{\rm{X}}}\rangle )/{\sigma }_{{\rm{tot}}}$, where $\langle {z}_{{\rm{X}}}\rangle $ is the average zX found across the analyzed regions for which we have realiable a spectral fit. Significant (larger than 2σ) deviations in the redshift map are thus visible as bright blue and bright red regions.

3.3. ICM Dynamical Analysis

In this section, we describe how we use the results of the spatially resolved spectral analysis to investigate and possibly characterize the presence of bulk motions in the ICM of each cluster. The first check, analogous to what we did in Liu et al. (2015), is to perform a simple χ2 test on the distribution of zX under the assumption of a unique value for the ICM redshift. Therefore, we compute ${\chi }^{2}\equiv {\rm{\Sigma }}{({z}_{{\rm{X}}}-\langle {z}_{{\rm{X}}}\rangle )}^{2}/{\sigma }_{{\rm{tot}}}^{2}$ and, assuming Nreg − 1 degrees of freedom (dof), where Nreg is the number of selected regions, we evaluate the probability of being inconsistent with the constant redshift hypothesis. To account for the small differences in the lower and upper error bars, we use the corresponding total error for zX values lower or higher than $\langle {z}_{{\rm{X}}}\rangle $.

We also compute the typical bulk motion velocity as the excess of the redshift fluctuations with respect to the statistical and systematic noise. To do that, we simply compute the total root mean square deviation of zX with respect to the mean, defined as ${\sigma }_{{\rm{rms}}}^{2}\equiv {\rm{\Sigma }}{({z}_{{\rm{X}}}-\langle {z}_{{\rm{X}}}\rangle )}^{2}$, and compare it to the average uncertainty $\langle {\sigma }_{{\rm{tot}}}\rangle ={\rm{\Sigma }}{\sigma }_{{\rm{tot}}}/{N}_{{\rm{reg}}}$ (where we use the average value of the lower and upper error bars when they differ). The average excess, which can be ascribed to bulk motions, is computed as ${\sigma }_{{\rm{BM}}}^{2}={\sigma }_{{\rm{rms}}}^{2}-\langle {\sigma }_{{\rm{tot}}}^{2}\rangle $. Therefore, the bulk motion velocity averaged over the region included in the maximum S/N radius is finally estimated as ${v}_{{\rm{BM}}}=c\times {\sigma }_{{\rm{BM}}}/(1+\langle z\rangle )$. The error on vBM is estimated as $\sqrt{2}$ times the total error on the mean. Moreover, in those cases where there is some evidence of bulk motions in the ICM, we also provide the value of the maximum velocity difference Δvmax found across all of the regions with a reliable spectral fit.

We visually inspect the histogram distribution of zX and the redshift and significance maps to infer information on the dynamical status of the ICM. This was done in Liu et al. (2015) for the Bullet cluster, where the two regions with a maximally different redshift are found to be aligned along the bullet trail. We tentatively interpreted them as two regions pushed in opposite directions and perpendicularly to the direction of the bullet at velocities ∼5−6 × 103 km s−1 with respect to the average cluster redshift. Clearly, it is not always possible to draw simple pictures like this one. However, we try to achieve a simple qualitative classification as pre-merger, ongoing, or post-merger according to whether we can clearly identify two well-defined halos with distinct redshift or a few scattered regions with different redshift. Along this same line, clusters with disturbed surface brightness but no signs of bulk motion may be described as post-merger, a phase in which the bulk motions created by the merger already evolved toward a turbulent velocity field with velocity values below the sensitivity achievable in our data. As previously discussed, the regime of evolved merger, when most of the initial kinetic energy due to the merger evolved into a turbulent velocity field of a few ×100 km s−1, can be probed only with the much higher spectral resolution of X-ray bolometers.

4. RESULTS ON INDIVIDUAL CLUSTERS

In this section, we show the results of our analysis. For each of the clusters, we measure the best-fit redshift in all of the selected regions, and hence obtain the redshift map and the significance map as defined in Section 3.2. We also measure the global redshift zXtot, which is defined as the best-fit redshift of the total emission of the cluster within the maximum S/N circle. This value corresponds to the emission-weighted average redshift of the ICM, which in principle may differ significantly from the simple average redshift among the regions $\langle z\rangle $. The statistical significance of the presence of bulk motions is evaluated in each case. Whenever possible, we estimate the typical velocity of the bulk motion vBM and the maximum velocity difference Δvmax. Moreover, for each cluster, we provide a short summary of the available literature concerning their dynamical status.

4.1. Abell 2142

Abell 2142 has already been identified as a merger cluster on the basis of ROSAT observations (Henry & Briel 1996), and more recently has been classified as "non-relaxed" on the basis of the X-ray morphology obtained with Chandra data by Parekh et al. (2015). The two cold fronts, identified in the northwest and southeast regions, clearly separate cooler gas (kT ∼ 7–9 keV) from hotter (kT > 12 keV; Markevitch et al. 2000), making it a promising target to investigate the onset of Kelvin–Helmholtz instabilities at the cold fronts. Evidence of dynamical substructure was also found with optical spectroscopy of the member galaxies (Oegerle et al. 1995). Interestingly, the presence of a narrow-angle tailed radio galaxy may be explained by significant bulk motions in the ICM, as suggested by the statistical analysis of Bliton et al. (1998). Finally, Cuciti et al. (2015) reported A2142 to have a radio halo.

We use the three most recent observation with Chandra ACIS-S for a total of 155.1 ks, as listed in Table 1, while we discard two older observations with ACIS-I (ObsID 5005 and 7692). This choice will minimize the difference in calibration between ObsID. The total number of net counts in the 0.5–10 keV band within a circle of 3 arcmin is about 106. This makes A2142 an ideal target for our analysis, with 52 potentially useful regions for spatially resolved spectral analysis. The best-fit redshift obtained by fitting the total X-ray emission within 3 arcmin is ${z}_{{\rm{Xtot}}}={0.0852}_{-0.0009}^{+0.0012}$. To evaluate the systematic uncertainty on zX, we consider a temperature range from 3 to 16 keV to explore the temperature structure of the ICM, despite the fact that the projected temperature shows only a moderate range from 8 to 12 keV. Metal abundances are mildly scattered around Z = 0.45 Z in the units provided by Asplund et al. (2005), nevertheless we consider a wide range $0.15\lt Z/{Z}_{\odot }\lt 1.5$ in order to account for possible large fluctuations. In 5 regions, we find that the X-ray spectral analysis does not provide a reliable value of zX, and therefore we use only 47 regions. The best-fit zX of the regions used in our analysis with the corresponding statistical and systematic errors are listed in Table 5. The average redshift is $\langle {z}_{{\rm{X}}}\rangle =0.0855$, which is almost coincident with the global, emission-weighted redshift zXtot.

The hypothesis of a constant redshift is rejected at more than 3σ, with a χ2 = 88.02 for 46 dof, as listed in Table 4. However, the histogram distribution of zX (see Figure 1) clearly shows that the evidence of bulk motions comes from three regions (region 2, 8, and 50) maximally distant from the average redshift. The rest of the zX values are distributed according to a Gaussian slightly narrower than the Gaussian whose variance is the average σtot2. This is expected since we assumed a very conservative estimate for σtot. The projected X-ray redshift map and significance map are shown in Figure 2. We measure vBM = 1400 ± 300 km s−1 for the average bulk motion velocity in the ICM of A2142 within a radius of 3 arcmin. Excluding the three extreme zX regions, we are able to establish a hard upper limit of vBM < 1450 km s−1. This result is consistent with the Suzaku X-ray spectral analysis by Ota & Yoshida (2015), which provided Δv < 4200 km s−1. However, thanks to the Chandra angular resolution, we find that region 50 is redshifted with respect to the average velocity by 5000 ± 1500 km s−1, while the contiguous regions 2 and 8 are blueshifted by 7400 ± 2200 km s−1. In particular, we can interpret region 50 as being dominated by a significant amount of ICM mass pushed at high velocity as a result of the ongoing merger. The two blueshifted regions can be tentatively associated with some rotation of the ICM. We note that, given the small field of view of the ACIS-S detector, we repeated the spectral analysis with a different choice of background, finding the same statistical results. Therefore, our conclusions are not affected by uncertainties in the background, despite the fact that the cluster emission covers almost the entire ACIS-S field of view.

Figure 1.

Figure 1. Histogram distribution (solid red line) of the best-fit zX for the 47 regions of Abell 2142 with reliable spectral fits. The green line is the Gaussian with ${\sigma }_{G}=\langle {\sigma }_{{\rm{tot}}}\rangle $ centered on $\langle {z}_{{\rm{X}}}\rangle $.

Standard image High-resolution image
Figure 2.

Figure 2. Top left panel: surface brightness distribution of Abell 2142 in the 0.5–10 keV band; top right panel: region map; bottom left panel: redshift map; bottom right panel: significance map. White regions do not have a reliable redshift measurement.

Standard image High-resolution image

Table 4.  The Results of the χ2 Test on the Distribution of Redshift of the Clusters, and Corresponding Constraints on vBM and Δvmax in km s−1

Cluster Name zXtot $\langle {z}_{{\rm{X}}}\rangle $ Degrees of Freedom χ2 Prob of BM vBM Δvmax Classification
      (Nreg − 1)     (km s−1) (km s−1)  
A2142 ${0.0852}_{-0.0009}^{+0.0012}$ 0.0855 46 88.02 99.98% 1400 ± 300 7400 ± 2200 post-merger
A2034 ${0.1124}_{-0.0044}^{+0.0039}$ 0.1122 10 12.70 76% 1700 ± 840 6200 ± 3300 pre/ongoing merger
A115 ${0.1892}_{-0.0040}^{+0.0030}$ 0.1979 11 30.4 99.9% 4600 ± 1100 8600 ± 3000 pre/ongoing merger
A520 ${0.2082}_{-0.0049}^{+0.0046}$ 0.2083 17 23.4 86% 1800 ± 900 5900 ± 4000 post-merger
1RXS J0603.3+4214 ${0.2316}_{-0.0033}^{+0.0029}$ 0.2303 7 8.6 72% $\lt 2100$ 7200 ± 2800 pre-merger
A2146 ${0.2310}_{-0.0010}^{+0.0022}$ 0.2338 16 19.68 77% 480 ± 490 8200 ± 3600 pre/ongoing merger
A1689 ${0.1814}_{-0.0008}^{+0.0028}$ 0.1814 9 9.3 59% <1600 4600 ± 3000 relaxed
A1835 ${0.2478}_{-0.0012}^{+0.0021}$ 0.2483 10 11.5 38% <1350 4200 ± 2600 relaxed

Note. This information is combined with visual inspection of the redshift maps to obtain the qualitative classification listed in the last column.

Download table as:  ASCIITypeset image

Table 5.  Spectral Analysis Results of Abell 2142

Region z ${\sigma }_{{\mathrm{stat}}_{b}}$ ${\sigma }_{{\mathrm{stat}}_{u}}$ ${\sigma }_{{\mathrm{syst}}_{b}}$ ${\sigma }_{{\mathrm{syst}}_{u}}$ ${\sigma }_{{\mathrm{tot}}_{b}}$ ${\sigma }_{{\mathrm{tot}}_{u}}$
0 0.0839 −0.0030 0.0033 −0.0009 0.0002 −0.0031 0.0034
1 0.0847 −0.0021 0.0026 −0.0007 0.0013 −0.0022 0.0029
2 0.0674 −0.0048 0.0050 −0.0014 0.0005 −0.0051 0.0050
3 0.0881 −0.0013 0.0029 −0.0001 0.0018 −0.0014 0.0034
4 0.0762 −0.0064 0.0061 −0.0040 0.0001 −0.0076 0.0061
5 0.0863 −0.0065 0.0058 −0.0054 0.0026 −0.0084 0.0063
6 0.0917 −0.0022 0.0024 −0.0001 0.0013 −0.0022 0.0027
7 0.0866 −0.0031 0.0032 −0.0001 0.0015 −0.0031 0.0035
8 0.0612 −0.0064 0.0075 −0.0001 0.0455 −0.0064 0.0461
10 0.0808 −0.0039 0.0039 −0.0001 0.0040 −0.0039 0.0055
11 0.0843 −0.0030 0.0032 −0.0001 0.0020 −0.0030 0.0038
12 0.0885 −0.0042 0.0098 −0.0001 0.0010 −0.0042 0.0099
13 0.0946 −0.0034 0.0036 −0.0001 0.0014 −0.0034 0.0039
14 0.0787 −0.0043 0.0036 −0.0007 0.0013 −0.0044 0.0039
15 0.0874 −0.0043 0.0036 −0.0001 0.0016 −0.0043 0.0040
16 0.0860 −0.0036 0.0038 −0.0001 0.0030 −0.0036 0.0049
17 0.0888 −0.0029 0.0031 −0.0020 0.0001 −0.0035 0.0031
18 0.0908 −0.0032 0.0024 −0.0001 0.0013 −0.0032 0.0027
20 0.0844 −0.0029 0.0028 −0.0004 0.0015 −0.0030 0.0031
21 0.0885 −0.0061 0.0083 −0.0015 0.0025 −0.0063 0.0087
22 0.0852 −0.0033 0.0033 −0.0012 0.0020 −0.0036 0.0039
23 0.0841 −0.0037 0.0038 −0.0001 0.0020 −0.0037 0.0043
25 0.0848 −0.0028 0.0027 −0.0001 0.0022 −0.0028 0.0034
26 0.0811 −0.0027 0.0029 −0.0010 0.0010 −0.0028 0.0031
27 0.0882 −0.0040 0.0038 −0.0010 0.0001 −0.0041 0.0038
28 0.0831 −0.0029 0.0024 −0.0010 0.0001 −0.0031 0.0024
29 0.0777 −0.0044 0.0045 −0.0017 0.0043 −0.0047 0.0062
31 0.0914 −0.0036 0.0040 −0.0001 0.0025 −0.0036 0.0047
32 0.0869 −0.0023 0.0022 −0.0010 0.0010 −0.0025 0.0024
33 0.0878 −0.0071 0.0100 −0.0028 0.0132 −0.0076 0.0166
34 0.0847 −0.0027 0.0027 −0.0007 0.0007 −0.0028 0.0028
35 0.0966 −0.0083 0.0058 −0.0087 0.0043 −0.0120 0.0072
36 0.0846 −0.0027 0.0026 −0.0006 0.0024 −0.0028 0.0035
37 0.0821 −0.0053 0.0050 −0.0010 0.0060 −0.0054 0.0078
38 0.0808 −0.0045 0.0046 −0.0013 0.0032 −0.0047 0.0056
39 0.0853 −0.0034 0.0035 −0.0001 0.0016 −0.0034 0.0038
40 0.0832 −0.0033 0.0035 −0.0007 0.0023 −0.0034 0.0042
41 0.0826 −0.0055 0.0041 −0.0001 0.0034 −0.0055 0.0053
42 0.0876 −0.0024 0.0027 −0.0001 0.0030 −0.0024 0.0040
43 0.0820 −0.0049 0.0053 −0.0001 0.0050 −0.0049 0.0073
44 0.0848 −0.0027 0.0023 −0.0008 0.0012 −0.0028 0.0026
45 0.0884 −0.0054 0.0051 −0.0024 0.0036 −0.0059 0.0062
47 0.0823 −0.0037 0.0036 −0.0001 0.0010 −0.0037 0.0037
48 0.0885 −0.0029 0.0031 −0.0001 0.0010 −0.0029 0.0033
49 0.0866 −0.0028 0.0030 −0.0001 0.0014 −0.0028 0.0033
50 0.1124 −0.0052 0.0063 −0.0025 0.0045 −0.0057 0.0077
51 0.0881 −0.0049 0.0044 −0.0001 0.0052 −0.0049 0.0068

Note. Column 1: region number; Column 2: best-fit redshift zX obtained fitting the 2.0–10 keV energy range; Columns 3 and 4: lower and upper 1σ error bars from fit statistics; Columns 5 and 6: lower and upper 1σ error bars from systematics associated to the ICM temperature structure; Columns 7 and 8: total lower and upper 1σ error bars computed as ${\sigma }_{{\rm{tot}}}=\sqrt{{\sigma }_{{\rm{stat}}}^{2}+{\sigma }_{{\rm{syst}}}^{2}}.$

Download table as:  ASCIITypeset image

To summarize, as in the case of the Bullet cluster (Liu et al. 2015), the evidence of bulk motions is coming from a few regions with an extreme redshift difference with respect to the mean after a very conservative selection of reliable spectra and an accurate estimate of the total uncertainty on zX. This may suggest that the cluster is a post-merger or caught in a stage where the shocks have developed since ∼1 dynamical time or more and the velocity field of the ICM still shows a few regions with large velocities with respect to the velocity of the center of mass. The interpretation in terms of ICM dynamics is nevertheless still very hard, and a significantly higher spectral resolution with the same angular resolution is needed before the actual velocity field of the ICM can be accurately measured.

4.2. Abell 2034

Abell 2034, at redshift zopt = 0.113 (Muriel & Coenda 2014), shows many signatures typically associated with merging clusters, from a radio relic near the position of a discontinuity in the X-ray surface brightness, to significant heating of the ICM above its equilibrium temperature, and a significant offset of the cD galaxy from the centroid of the X-ray emission (Kempner & Sarazin 2001; Kempner et al. 2003; Giovannini et al. 2009; Van Weeren et al. 2011). In a recent detailed study by Owers et al. (2014) using a deep Chandra observation of 250 ks, a shock with a Mach number of 1.59 ± 0.06, corresponding to a shock velocity of ∼2000 km s−1, was clearly observed. Combining Chandra data with spectroscopic observations for 328 spectroscopically confirmed member galaxies, they show the presence of a substructure located at the front edge of the shock and that the merger is proceeding along a north–south direction with an inclination angle of ∼23° with respect to the plane of the sky. The cluster is included in the sample of the Merging Cluster Collaboration.5

We use Chandra ACIS-I observations for a total of 254.5 ks. The total number of net counts in the 0.5–10 keV band is 2.66 × 105 within a circle of 3 arcmin in the merged image. We select 14 ICM regions for redshift analysis, but we obtained reliable spectral fits only in 11 out of 14. We find ${z}_{{\rm{Xtot}}}={0.1124}_{-0.0044}^{+0.0039}$, which is in very good agreement with zopt and $\langle {z}_{{\rm{X}}}\rangle =0.1122$. To evaluate the systematic uncertainties on zX, we consider a temperature and abundance range of 6–12 keV and 0.2–0.8 Z, respectively. The best-fit zX of the 11 regions used in our analysis with corresponding errors are listed in Table 6. The histogram distribution of zX shows that most of the best-fit zX values are approximately distributed as a Gaussian, with the exception of two high-redshift values (see Figure 3). However, due to the estimated uncertainties, the probability of a uniform redshift distribution across the ICM is still high at the 24% level (χ2 = 12.70 for 10 dof), which implies no significant detection (less than 2σ) of global bulk motions. Formally, we measure an average bulk motion of vBM = 1700 ± 840 km s−1. In addition, the projected redshift map shows that the two high-redshift regions 4 and 9 are contiguous (see Figure 4, lower panels), which can be interpreted as an enhanced likelihood of a redshift difference between the diffuse emission in the north and the main halo. The average redshift of regions 4 and 9 is 0.1310, while that of the other regions is 0.1081. This implies a velocity difference of Δvmax = 6200 ± 3300 km s−1. The presence of two distinct halos with different redshift may suggest a classification as a pre-merger. Due to the presence of a well-developed shock, this cluster may be caught in the early stages of the merger process, in agreement with the conclusion of Owers et al. (2014), who estimate that the merger is observed only ∼0.3 Gyr after core-passage.

Figure 3.

Figure 3. Histogram distribution of the best-fit zX for the 11 regions of Abell 2034 with reliable spectral fit. Lines are as in Figure 1.

Standard image High-resolution image
Figure 4.

Figure 4. Same as Figure 2, but for Abell 2034.

Standard image High-resolution image

Table 6.  Spectral Analysis Results of Abell 2034

Region z ${\sigma }_{{\mathrm{stat}}_{b}}$ ${\sigma }_{{\mathrm{stat}}_{u}}$ ${\sigma }_{{\mathrm{syst}}_{b}}$ ${\sigma }_{{\mathrm{syst}}_{u}}$ ${\sigma }_{{\mathrm{tot}}_{b}}$ ${\sigma }_{{\mathrm{tot}}_{u}}$
0 0.1056 −0.0057 0.0056 −0.0001 0.0033 −0.0057 0.0065
1 0.1085 −0.006 0.0068 −0.0016 0.0014 −0.0062 0.0069
2 0.1004 −0.0062 0.0088 −0.0001 0.0054 −0.0062 0.0103
3 0.1089 −0.0074 0.008 −0.0017 0.0053 −0.0075 0.0095
4 0.1318 −0.0093 0.0076 −0.0035 0.0001 −0.0099 0.0076
6 0.1116 −0.0038 0.0055 −0.0001 0.0018 −0.0038 0.0057
9 0.1301 −0.0067 0.0058 −0.0049 0.0001 −0.0083 0.0058
10 0.1124 −0.0055 0.0059 −0.0018 0.0001 −0.0057 0.0059
11 0.1099 −0.0077 0.0061 −0.0051 0.0001 −0.0092 0.0061
12 0.1142 −0.0075 0.0102 −0.0035 0.0018 −0.0082 0.0103
13 0.1011 −0.0068 0.0095 −0.0001 0.0040 −0.0068 0.0103

Note. The columns are the same as in Table 5.

Download table as:  ASCIITypeset image

4.3. Abell 115

Abell 115 has been classified as a merger on the basis of its double X-ray morphology (Forman et al. 1981). The ICM thermal structure is best described by two halos that host a cool core with significantly hotter gas between and around them, suggesting that the system is in an advanced stage of merging, despite the fact that no shock front has been found in this region (Shibata et al. 1999; Gutierrez & Krawczynski 2005). A dynamical analysis based on 88 cluster member spectra (Barrena et al. 2007) shows a line-of-sight velocity difference of 1646 km s−1 between the northern and southern subclusters at a projected separation of 0.89 Mpc. In addition, Barrena et al. (2007) interpret the system as a pre-merger based on the location of the BCGs being consistent with the peaks in the X-ray surface brightness distribution. Diffuse radio emission north of the cluster has been found by Govoni et al. (2001) with a 1.4 GHz VLA observation and interpreted as a radio relic that may be associated with the ongoing merger. Abell 115 is also included in the Merging Cluster Collaboration sample.

We used 4 ACIS-I contiguous-in-time observations for a total of 310.6 ks. The total number of counts in the 0.5–10 keV band within a circle of 5.2 arcmin in the merged image is 2.88 × 105. This allows us to obtain 18 useful regions for spectral analysis. Our spatially resolved spectroscopic analysis provides us with reliable zX only in 12 regions out of 18. The average redshift $\langle {z}_{{\rm{X}}}\rangle =0.1979$ is in good agreement with zopt = 0.197 (Hiroi et al. 2013) while the global redshift ${z}_{{\rm{Xtot}}}={0.1892}_{-0.004}^{+0.003}$ is significantly lower. To evaluate systematic uncertainties on zX, we consider temperature and abundance ranges of 2–12 keV and 0.2–0.9 Z, respectively. The best-fit zX of the 12 regions used in our analysis with corresponding errors are listed in Table 7.

Table 7.  Spectral Analysis Results of Abell 115

Region z ${\sigma }_{{\mathrm{stat}}_{b}}$ ${\sigma }_{{\mathrm{stat}}_{u}}$ ${\sigma }_{{\mathrm{syst}}_{b}}$ ${\sigma }_{{\mathrm{syst}}_{u}}$ ${\sigma }_{{\mathrm{tot}}_{b}}$ ${\sigma }_{{\mathrm{tot}}_{u}}$
0 0.1980 −0.0019 0.0036 −0.0000 0.0025 −0.0019 0.0044
1 0.1998 −0.0027 0.0028 −0.0000 0.0002 −0.0027 0.0028
2 0.2026 −0.0036 0.0039 −0.0000 0.0037 −0.0036 0.0053
4 0.1988 −0.0051 0.0042 −0.0000 0.0061 −0.0051 0.0074
6 0.1922 −0.0067 0.0083 −0.0032 0.0018 −0.0074 0.0085
8 0.1897 −0.0124 0.0123 −0.0037 0.0053 −0.0129 0.0134
10 0.2090 −0.0082 0.0086 −0.0040 0.0040 −0.0091 0.0095
11 0.2367 −0.0060 0.0048 −0.0087 0.0000 −0.0106 0.0048
12 0.1679 −0.0094 0.0113 −0.0049 0.0020 −0.0106 0.0115
14 0.1591 −0.0064 0.0094 −0.0000 0.0210 −0.0064 0.0230
15 0.2263 −0.0177 0.0107 −0.0000 0.0000 −0.0177 0.0107
17 0.1943 −0.0069 0.0074 −0.0038 0.0107 −0.0079 0.0130

Note. The columns are the same as in Table 5.

Download table as:  ASCIITypeset image

The hypothesis of a constant redshift across the 12 regions is rejected at more than 3σ (χ2 = 30.4 for 11 dof). The histogram distribution of zX (see Figure 5) shows the presence of two regions significantly above and below the central value. We measure an average bulk motion velocity of vBM = 4600 ± 1100 km s−1.

Figure 5.

Figure 5. Histogram distribution of the best-fit zX for the 12 regions of Abell 115 with reliable spectral fit. Lines are as in Figure 1.

Standard image High-resolution image

Visual inspection of the redshift and significance maps (see Figure 6, lower panels) allows us to identify three main regions with compatible redshift. We first isolate regions belonging to the northern (regions 0, 1, 2, 4, and 10) and southern (3, 6, 11, 17, and 15) clumps. We find a redshift difference that corresponds to about 1σ, and is therefore consistent with having the same redshift. Formally, translated into a velocity difference, this would imply a relative velocity of 3300 ± 3200 km s−1. If we consider instead regions 12 and 14, then we find an average redshift of 0.1634, corresponding to a velocity difference of 8600 ± 3000 km s−1 with respect to the bulk of the ICM. We note that this is also the gas that is found to be significantly hotter (at least by a factor of two, as verified in the projected temperature map not discussed here), and therefore it is likely to be heated and violently displaced by the ongoing merger. This is consistent with a substantial amount of the shocked gas being compressed between the two halos and pushed toward the observer, with the two halos having a minimum in the velocity difference. Therefore, we conclude that A115 is caught in the early phase of the merging process, with the two ICM halos still well separated.

Figure 6.

Figure 6. Same as Figure 2, but for Abell 115.

Standard image High-resolution image

4.4. Abell 520

Abell 520, at zopt = 0.203 (Cassano et al. 2013), was found to host a radio halo elongated in the NE-SW direction, corresponding to the apparent merger axis (Giovannini et al. 1999; Govoni et al. 2001). Its X-ray morphology and the optical spectroscopy suggest a disturbed dynamical state (Proust et al. 2000; Govoni et al. 2001, respectively). Despite recent in-depth studies (Mahdavi et al. 2007; Girardi et al. 2008), the detailed structure of the merger of Abell 520 is still unclear. However, a prominent bow shock with $M={2.1}_{-0.3}^{+0.4}$ has been identified thanks to Chandra data analysis (Markevitch et al. 2005). Abell 115 is also included in the Merging Cluster Collaboration sample.

A520 has the longest Chandra ACIS-I exposure for a total of 516.4 ks. The number of net counts in the 0.5–10 keV band in the merged image is 3.36 × 105 in a circle of 3.2 arcmin radius. The best-fit redshift from the global emission is ${z}_{{\rm{Xtot}}}={0.2082}_{-0.0049}^{+0.0046}$. To estimate the systematic uncertainty on zX, we consider ranges of 3–15 keV and 0.2–1 Z for temperature and abundance, respectively.

We discard only two regions that do not provide reliable redshifts. Therefore, we base our analysis on 18 regions, whose best-fit zX and associated errors are listed in Table 8. The average redshift among those regions is $\langle {z}_{{\rm{X}}}\rangle =0.2083$, in good agreement with zXtot and zopt. The hypothesis of a constant redshift is rejected at an ∼87% confidence level (χ2 = 23.4 for 17 dof), which is therefore less than 2σ. We find vBM ∼ 1800 ± 900 km s−1. This is mostly contributed by the presence of blueshifted regions, mostly in the central 1 arcmin, and the surrounding higher-redshift regions, as can be seen by comparing the histogram distribution (see Figure 7) with the redshift and significance maps (see Figure 8). If we compute the difference between the two sets of regions with comparable redshifts, then we find Δvmax = 5900 ± 4000 km s−1. The overall redshift map suggests a merger with a halo moving toward the observer that was caught shortly after the first pericentric passage.

Figure 7.

Figure 7. Histogram distribution of the best-fit zX for the 18 regions of Abell 520 with a reliable spectral fit. Lines are as in Figure 1.

Standard image High-resolution image
Figure 8.

Figure 8. Same as Figure 2, but for Abell 520.

Standard image High-resolution image

Table 8.  Spectral Analysis Results of Abell 520

Region z ${\sigma }_{{\mathrm{stat}}_{b}}$ ${\sigma }_{{\mathrm{stat}}_{u}}$ ${\sigma }_{{\mathrm{syst}}_{b}}$ ${\sigma }_{{\mathrm{syst}}_{u}}$ ${\sigma }_{{\mathrm{tot}}_{b}}$ ${\sigma }_{{\mathrm{tot}}_{u}}$
0 0.2102 −0.0037 0.0065 −0.0001 0.0021 −0.0037 0.0068
1 0.2069 −0.0047 0.0043 −0.0001 0.0027 −0.0047 0.0050
2 0.1933 −0.0085 0.0096 −0.0001 0.0074 −0.0085 0.0121
3 0.1962 −0.0075 0.0076 −0.0001 0.0024 −0.0075 0.0079
5 0.2276 −0.0155 0.0147 −0.0025 0.0049 −0.0157 0.0154
6 0.2207 −0.0146 0.0079 −0.0122 0.0001 −0.0190 0.0079
8 0.1937 −0.0128 0.0133 −0.0018 0.0056 −0.0129 0.0144
9 0.2048 −0.0087 0.0075 −0.0006 0.0027 −0.0087 0.0079
10 0.2271 −0.0164 0.0138 −0.0039 0.0001 −0.0168 0.0138
11 0.2152 −0.0106 0.0096 −0.0012 0.0058 −0.0106 0.0112
12 0.2138 −0.0046 0.0053 −0.0001 0.0044 −0.0046 0.0068
13 0.2282 −0.0149 0.0126 −0.0175 0.0006 −0.0229 0.0126
14 0.2217 −0.0153 0.0166 −0.0066 0.0001 −0.0166 0.0166
15 0.1921 −0.0088 0.0085 −0.0083 0.0001 −0.0120 0.0085
16 0.195 −0.0124 0.0131 −0.0081 0.0001 −0.0148 0.0131
17 0.1837 −0.0085 0.0086 −0.0001 0.0038 −0.0085 0.0094
18 0.208 −0.0102 0.0068 −0.0057 0.0005 −0.0116 0.0068
19 0.2104 −0.0099 0.0093 −0.0001 0.0114 −0.0099 0.0147

Note. The columns are the same as in Table 5.

Download table as:  ASCIITypeset image

4.5. 1RXS J0603.3+4214

At redshift zopt = 0.225 (Van Weeren et al. 2012), 1RXS J0603.3+4214 (the "Toothbrush" cluster) was discovered in the radio band thanks to its prominent linear radio relic and radio halo plus several additional features which suggest that this is a complex merger. Simulations also show that the merger might be a triple process: in addition to northern and southern subclusters with similar mass, there might be another smaller subcluster infalling in the south (Brüggen et al. 2012). Measurement of the radio spectral index allowed Van Weeren et al. (2012) to constrain the merger Mach number to be in the range 3.3–4.6. However, X-ray observation by XMM-Newton returned a Mach number of <2 (Ogrean et al. 2013), which is inconsistent with the radio constraints. 1RXS J0603.3+4214 is also included in the Merging Cluster Collaboration sample.

We use three Chandra ACIS-I observations for a total of 235.9 ks. The total number of counts in the 0.5–10 keV band within a circle of 3.5 arcmin is 1.74 × 105 in the merged image. To evaluate the systematic uncertainty on zX, we consider ranges of 8–16 keV and 0.3–0.8 Z for the temperature and abundance, respectively. We are able to obtain reliable spectral fits in 8 out of the 9 regions originally considered for analysis. The best-fit zX and corresponding errors are listed in Table 9. The average redshift across the regions is $\langle z\rangle =0.2303$, which is in good agreement with ${z}_{{\rm{Xtot}}}={0.2316}_{-0.0033}^{+0.0029}$ and in reasonable agreement with zopt. The hypothesis of a constant redshift is acceptable (χ2 = 8.7 for 7 dof), as also shown in Figure 9 by the histogram distribution of the best-fit zX values. The upper limit to the average bulk motion velocity is vBM < 2100 km s−1. However, there is a significant uniform gradient in the best-fit redshift from the northern clump (high-z) to the southern clump (low-z; see Figure 10), which corresponds to a maximum difference of ∼0.0295 ± 0.0115 or a maximum velocity difference of Δvmax ∼ 7200 ± 2800 km s−1. On the basis of these results, this cluster can be classified as a pre-merger.

Figure 9.

Figure 9. Histogram distribution of the best-fit zX for the 8 regions of 1RXS J0603.3+4214 with reliable spectral fit. Lines are as in Figure 1.

Standard image High-resolution image
Figure 10.

Figure 10. Same as Figure 2, but for 1RXS J0603.3+4214.

Standard image High-resolution image

Table 9.  Spectral Analysis Results of 1RXS J0603.3+4214

Region z ${\sigma }_{{\mathrm{stat}}_{b}}$ ${\sigma }_{{\mathrm{stat}}_{u}}$ ${\sigma }_{{\mathrm{syst}}_{b}}$ ${\sigma }_{{\mathrm{syst}}_{u}}$ ${\sigma }_{{\mathrm{tot}}_{b}}$ ${\sigma }_{{\mathrm{tot}}_{u}}$
0 0.2317 −0.0062 0.0065 −0.0030 0.0004 −0.0068 0.0065
1 0.2206 −0.0069 0.0072 −0.0021 0.0065 −0.0072 0.0097
2 0.2347 −0.0075 0.0119 −0.0008 0.0071 −0.0075 0.0138
4 0.2493 −0.0067 0.0091 −0.0069 0.0001 −0.0096 0.0091
5 0.2198 −0.0070 0.0063 −0.0045 0.0001 −0.0083 0.0063
6 0.2254 −0.0108 0.0078 −0.0055 0.0001 −0.0121 0.0078
7 0.2294 −0.0094 0.0066 −0.0044 0.0021 −0.0103 0.0069
8 0.2316 −0.0091 0.0075 −0.0001 0.0042 −0.0091 0.0085

Note. The columns are the same as in Table 5.

Download table as:  ASCIITypeset image

4.6. Abell 2146

Despite the lack of diffuse radio emission, A2146 is classified as an extreme binary, head-on cluster merger in Mann & Ebeling (2012). A deep Chandra study reveled the presence of two shock fronts with Mach numbers of M = 2.2 ± 0.8 and M ∼ 1.7 ± 0.3 (Russell et al. 2010). Using the radial velocity of BCGs in the subcluster and the main cluster, Canning et al. (2012) find that the merger axis has an angle of 17° with respect to the plane of the sky. Based on these properties and its appearance in the X-ray, Abell 2146 is often compared to the Bullet cluster.

We use observations with Chandra ACIS-I for a total of 375.3 ks, as listed in Table 1. For simplicity, we do not consider observations with ACIS-S for a total of 45 ks in order to retain uniform calibration throughout our analysis, given the smaller amount of observing time contributed by the ACIS-S exposures. The total number of net counts in the 0.5–10 keV band within a circle of 3.6 arcmin is about 1.99 × 105. The best-fit redshift obtained by fitting the total X-ray emission is ${z}_{\mathrm{Xtot}}={0.2310}_{-0.0010}^{+0.0022}$. Measured temperatures range from the inner ∼3 keV to the outer ∼8 keV in our analysis of the temperature structure (not reported in this work). To evaluate the systematic uncertainty on zX, we consider ranges from 1 to 12 keV and 0.15 to 1.5 Z/Z to explore the temperature and metal abundance throughout the ICM, respectively.

We are able to obtain a reliable spectral fit in 17 out of the 19 regions originally selected in the X-ray image. The best-fit zX of the regions used in our analysis, along with their corresponding statistical and systematic errors, are listed in Table 10. The average redshift across the region is $\langle {z}_{{\rm{X}}}\rangle =0.2338$, which is consistent with zXtot within 1σ. We find χ2 = 19.68 for 16 dof, as listed in Table 4, which is pretty much consistent with a constant redshift across the ICM. The histogram distribution of zX, shown in Figure 11, is consistent with a Gaussian dominated by noise and shows no clear evidence of bulk motions. The average bulk motion is formally vBM = 480 ± 490 km s−1. The projected X-ray redshift map and significance map are shown in Figure 12. The maximum redshift difference is Δvmax = 8200 ± 3600 km s−1 and is obtained between two distinct regions with homogeneous redshift: the "bullet" and its trail, along the direction of the merger, which are blueshifted with respect to the average redshift, and all of the surrounding ICM emission, which is redshifted. Despite the low S/N, the redshift map depicts a clear dynamical situation where the central ICM is dominated by a halo moving toward the observer, while the larger scale emission has, on average, a larger redshift. This provides support for the presence of bulk motions that are still dominated by the infall velocity of the merging halos.

Figure 11.

Figure 11. Histogram distribution of the best-fit zX for the 17 regions of A2146 with reliable spectral fit. Lines are as in Figure 1.

Standard image High-resolution image
Figure 12.

Figure 12. Same as Figure 2, but for Abell 2146.

Standard image High-resolution image

Table 10.  Spectral Analysis Results of Abell 2146

Region z ${\sigma }_{{\mathrm{stat}}_{b}}$ ${\sigma }_{{\mathrm{stat}}_{u}}$ ${\sigma }_{{\mathrm{syst}}_{b}}$ ${\sigma }_{{\mathrm{syst}}_{u}}$ ${\sigma }_{{\mathrm{tot}}_{b}}$ ${\sigma }_{{\mathrm{tot}}_{u}}$
0 0.2292 −0.0045 0.00312787 −0.0000 0.0000 −0.0045 0.0031
1 0.2282 −0.0035 0.00311993 −0.0017 0.0018 −0.0039 0.0036
2 0.2277 −0.0038 0.00571345 −0.0000 0.0083 −0.0038 0.0100
3 0.2282 −0.0051 0.0026654 −0.0027 0.0000 −0.0057 0.0027
4 0.2294 −0.0037 0.00436681 −0.0000 0.0036 −0.0037 0.0057
5 0.2276 −0.0065 0.00615287 −0.0001 0.0044 −0.0065 0.0076
6 0.2325 −0.0074 0.00676923 −0.0055 0.0000 −0.0092 0.0068
7 0.2354 −0.0060 0.00619396 −0.0114 0.0000 −0.0129 0.0062
8 0.2473 −0.0125 0.0069235 −0.0223 0.0000 −0.0255 0.0069
9 0.2436 −0.0072 0.0066689 −0.0036 0.0000 −0.0080 0.0067
10 0.2277 −0.0083 0.00881529 −0.0000 0.0140 −0.0083 0.0165
12 0.2372 −0.0062 0.00682796 −0.0010 0.0010 −0.0063 0.0069
13 0.2326 −0.0088 0.00909447 −0.0025 0.0074 −0.0091 0.0117
14 0.2256 −0.0077 0.00950732 −0.0025 0.0086 −0.0081 0.0128
15 0.2495 −0.0089 0.00691656 −0.0015 0.0000 −0.0090 0.0069
16 0.2226 −0.0114 0.0110643 −0.0060 0.0000 −0.0129 0.0110
18 0.2505 −0.0075 0.00798135 −0.0080 0.0000 −0.0110 0.0080

Note. The columns are the same as in Table 5.

Download table as:  ASCIITypeset image

4.7. Abell 1689

Abell 1689 is a typical relaxed cluster at z ∼ 0.183 (Hiroi et al. 2013). We reduced and analyzed four ACIS-I observations for a total of 151.3 ks (see Table 2). The total number of net counts in the 0.5–10 keV band within a circle of 1 arcmin is about 1.86 × 105 in the ACIS-I merged image. The best-fit redshift obtained fitting the total X-ray emission is ${z}_{{\rm{Xtot}}}={0.1814}_{-0.0008}^{+0.0028}$. To evaluate the systematic uncertainty on zX, we consider ranges from 5 to 15 keV and 0.3 to 0.9 Z/Z to explore the temperature and metal abundance throughout the ICM, respectively.

We are able to obtain a reliable spectral fit in all 10 of the regions originally selected in the X-ray image. The best-fit zX of the regions used in our analysis with the corresponding statistical and systematic errors are listed in Table 11. The average redshift across the region is $\langle {z}_{{\rm{X}}}\rangle =0.1814$. The histogram distribution of zX, shown in Figure 13, is consistent with a Gaussian distribution whose σ is the average statistical error on zX. We find χ2 = 9.31 for 9 dof, as listed in Table 4. Therefore, we do not find any hint of bulk motion, as expected, and find a hard upper limit of vBM < 1600 km s−1.

Figure 13.

Figure 13. Histogram distribution of the best-fit zX for the 10 regions of A1689 with reliable spectral fit. Lines are as in Figure 1.

Standard image High-resolution image
Figure 14.

Figure 14. Same as Figure 2, but for Abell 1689.

Standard image High-resolution image

Table 11.  Spectral Analysis Results of Abell 1689

Region z ${\sigma }_{{\mathrm{stat}}_{b}}$ ${\sigma }_{{\mathrm{stat}}_{u}}$ ${\sigma }_{{\mathrm{syst}}_{b}}$ ${\sigma }_{{\mathrm{syst}}_{u}}$ ${\sigma }_{{\mathrm{tot}}_{b}}$ ${\sigma }_{{\mathrm{tot}}_{u}}$
0 0.1767 −0.0061 0.0048 −0.0001 0.0037 −0.0061 0.0060
1 0.1774 −0.0047 0.0061 −0.0014 0.0001 −0.0049 0.0061
2 0.1864 −0.0103 0.0109 −0.0001 0.0066 −0.0103 0.0127
3 0.1774 −0.0065 0.0049 −0.0019 0.002 −0.0067 0.0052
4 0.1675 −0.0074 0.0066 −0.0019 0.0001 −0.0076 0.0066
5 0.1859 −0.0045 0.0038 −0.0001 0.0016 −0.0045 0.0041
6 0.1895 −0.0084 0.0087 −0.0001 0.0385 −0.0084 0.0394
7 0.1859 −0.006 0.0078 −0.0001 0.0041 −0.0060 0.0088
8 0.1826 −0.0087 0.0073 −0.0001 0.005 −0.0087 0.0088
9 0.1847 −0.0047 0.0071 −0.0001 0.0031 −0.0047 0.0077

Note. The columns are the same as in Table 5.

Download table as:  ASCIITypeset image

4.8. Abell 1835

Abell 1835 is a relaxed cluster at z ∼ 0.234 (Girardi et al. 2014). We reduced and analyzed three ACIS-I observations for a total of 193.7 ks (see Table 2). The total number of net counts in the 0.5–10 keV band within a circle of 1 arcmin is about 2.33 × 105 in the ACIS-I merged data. The best-fit redshift obtained while fitting the total X-ray emission is ${z}_{{\rm{Xtot}}}={0.2478}_{-0.0012}^{+0.0021}$, which is significantly larger than the optical value. To evaluate the systematic uncertainty on zX, we consider ranges from 4 to 14 keV and 0.2 to 0.8 Z/Z to explore the temperature and metal abundance throughout the ICM, respectively.

We are able to obtain a reliable spectral fit in all 11 of the regions originally selected in the X-ray image. The best-fit zX of the regions used in our analysis with the corresponding statistical and systematic errors are listed in Table 12. The average redshift across the region is $\langle {z}_{{\rm{X}}}\rangle =0.2483$. The histogram distribution of zX, shown in Figure 15, is consistent with a Gaussian distribution whose σ is the average statistical error on zX. We find χ2 = 11.5 for 10 dof, as listed in Table 4. Therefore, we do not find any hint of bulk motion, as expected, and find a hard upper limit of vBM < 1350 km s−1.

Figure 15.

Figure 15. Histogram distribution of the best-fit zX for the 11 regions of A1835 with reliable spectral fit. Lines are as in Figure 1.

Standard image High-resolution image
Figure 16.

Figure 16. Same as Figure 2, but for Abell 1835.

Standard image High-resolution image

Table 12.  Spectral Analysis Results of Abell 1835

Region z ${\sigma }_{{\mathrm{stat}}_{b}}$ ${\sigma }_{{\mathrm{stat}}_{u}}$ ${\sigma }_{{\mathrm{syst}}_{b}}$ ${\sigma }_{{\mathrm{syst}}_{u}}$ ${\sigma }_{{\mathrm{tot}}_{b}}$ ${\sigma }_{{\mathrm{tot}}_{u}}$
0 0.2524 −0.0031 0.0041 −0.0024 0.0001 −0.0039 0.0041
1 0.2547 −0.0024 0.003 −0.0001 0.0025 −0.0024 0.0039
2 0.2478 −0.003 0.0052 −0.0001 0.0021 −0.0030 0.0056
3 0.2474 −0.003 0.0027 −0.002 0.0003 −0.0036 0.0027
4 0.2472 −0.0048 0.006 −0.002 0.0018 −0.0052 0.0062
5 0.2408 −0.0028 0.0045 −0.0001 0.0052 −0.0028 0.0068
6 0.2538 −0.01 0.0065 −0.0057 0.0016 −0.0115 0.0066
7 0.2519 −0.0046 0.0064 −0.0028 0.0001 −0.0053 0.0064
8 0.2486 −0.006 0.0052 −0.0001 0.0046 −0.0060 0.0069
9 0.2372 −0.0065 0.0073 −0.0002 0.007 −0.0065 0.0101
10 0.2496 −0.0068 0.0078 −0.0019 0.0023 −0.0070 0.0081

Note. The columns are the same as in Table 5.

Download table as:  ASCIITypeset image

5. STATISTICAL PROPERTIES OF THE CLUSTER SAMPLE

In Figure 17, we summarize our results, showing the χ2 values versus the dof for all of the clusters analyzed in this work and listed in Table 4. Comparison with the lines corresponding to the 1, 2, and 3σ confidence levels shows that, in general, merging clusters always have a larger probability of showing bulk motions in their ICM, but only in two cases do we have a statistical significance larger than 3σ. In two cases, instead, the probability is less than 2σ. Finally, in the last two cases, the global χ2 is not significantly different from that of the two relaxed clusters, which show no signs of bulk motions. Globally, the merging cluster sample always shows a larger probability of hosting bulk motions, but it is hard to reach definitive conclusions in single cases.

Figure 17.

Figure 17. Red squares show the χ2 values vs. the degrees of freedom (equal to Nreg − 1) for merger clusters, while the circled green points show the two relaxed clusters. Dashed lines show the χ2 values as a function of the degrees of freedom corresponding to the 1, 2, and 3σ confidence levels, from the lowest to the highest.

Standard image High-resolution image

However, we note that the analysis based simply on the χ2 including all of the ICM regions tends to dilute the signal, since the majority of the ICM regions are expected to show no significant bulk motions. In fact, the bulk motion signal is created by only a few, large masses of ICM displaced at large velocities during the merger and, in addition, can be partially erased by projection effects. On the bright side, our global evaluation of bulk motion is free from the "look-elsewhere effect" (Lyons 2008). We conclude that, under our most conservative assumptions, only a fraction of the major mergers shows unambiguous signatures of bulk motions in the Chandra CCD data. Specifically, we measure vBM = 1400 ± 300 and 4600 ± 1100 km s−1 in A2142 and A115, respectively. In A2034 and A520 we find vBM = 1700 ± 840 and 1800 ± 900, respectively, although this measurement is not significant (less then 2σ). In the cases of 1RXS J0603.3+4214 and A2146, we find hard (corresponding to 3σ) upper limits of 2100 and ∼1500 km s−1, respectively. These limits are similar to those found for the two relaxed clusters A1689 and A1835, where we find vBM < 1600 and vBM < 1350 km s−1, respectively.

Local velocity differences between ICM regions may achieve much larger values with corresponding larger errors. In the eighth column of Table 4, we list the values of Δvmax found by identifying the ICM regions (or groups of regions) with the largest difference in zX. Clearly, this does not constitute a detection of bulk motion, however, it does provide a rough estimate of the maximum velocity difference we can have across the ICM in our data.

Visual inspection of the redshift and significance maps also allows us to provide a rough classification in terms of pre-merger or ongoing/post-merger. The first case is when two comparable-mass halos are about to begin a head-on merger with large speed along the line of sight, so that in the redshift map the two halos appear to have clearly distinct average redshifts (see, for example, the case of 1RXS J0603.3+4214). The ongoing or post-merger case is when shocks in the ICM are clearly present and the first pericentric passage is recent (i.e., it happened less than one dynamical time ago, see A2034 and A115). Given the uncertainties on the bulk motion velocity, we are not able to explore the regime where the residual kinetic energy of the merger event has already decreased to a more diffuse and chaotic turbulent velocity field of the order of few ×100 km s−1. Our qualitative classification is also provided in the last columns of Table 4.

6. DISCUSSION AND FUTURE PROSPECTS

In this work, we use the same conservative approach that we adopted in Liu et al. (2015), where we explored the presence of bulk motions in the ICM of 1E0657-56. In particular, when estimating the statistical significance of bulk motions, we used five conservative assumptions.

  • 1.  
    We considered only the Fe line emission in the hard band. The use of the soft band would have significantly reduced the statistical error on zX, particulary for the cold regions. However, since, in principle, gas with different temperatures can have different positions along the line of sight and different velocities, this may introduce a source of scatter in zX which is hard to control.
  • 2.  
    We computed the total error on zX in each region, summing in quadrature the statistical error and the 90% percentile of the zX distribution obtained by imposing a varying thermal structure in the ICM. This systematic component is often smaller than, or at best comparable to, the statistical component. The systematic component is likely to be slightly overestimated, as we can verify by comparing the σrms of the zX distribution with the average total error $\langle {\sigma }_{{\rm{tot}}}\rangle $.
  • 3.  
    We discarded all of the regions where the best-fit redshift is not identified clearly, rejecting all of the spectra where the Cstat versus redshift plot shows a secondary minimum at a distance ΔCstat < 6.6 from the absolute minimum.
  • 4.  
    Although the signature of bulk motion is presumably confined to a few, largely discrepant regions, we evaluate the statistical significance with a simple χ2 test on all of the regions. This of course would dilute the effect, but provides a robust evaluation of the actual presence of bulk motions completely free from the "look-elsewhere effect" (Lyons 2008).
  • 5.  
    We use information from the spatial distribution of zX only to provide a qualitative description of the possible dynamical status of the cluster, and not to reinforce the statistical significance of the zX fluctuations. In practice, we do not redefine regions to merge them together according to their zX value, a procedure that would reduce the statistical uncertainty on the final zX maps.

Given our extremely conservative approach, we only hit the tip of the iceberg for signatures of bulk motions in Chandra CCD data.  At the same time, our results show both the difficulty in measuring bulk motions and also the widespread presence of bulk motions in massive, head-on mergers.

Based on these tantalizing but promising results, in the near future we plan to relax some of these constraints. In particular, we will include soft-band spectral information and quantitatively investigate the spatial correlation of the zX fluctuations, possibly merging projected ICM regions that have similar redshifts on the basis of the spectral analysis and surface brightness features. In addition, we can also combine X-ray spectrally resolved analyses with dynamical studies based on the optical redshift distribution of the member galaxies. We are currently applying the caustic method (Diaferio 1999; Diaferio et al. 2005; Serra & Diaferio 2013; Yu et al. 2015) to a few nearby, massive clusters in the Abell cluster catalog (Abell et al. 1989) with the aim of carefully comparing the optical dynamical structure to the X-ray morphology and spatially resolved spectral analysis (H. Yu et al. 2016, in preparation).

Clearly, any improvement that may be achieved with current facilities is limited by the modest spectral resolution. In particular, the inevitable superposition of ICM massive regions with different redshifts along the line of sight heavily blurs the redshift map. Sometimes the presence of more than one minima in the CCD spectra may be interpreted, in principle, as the superposition of two ICM regions moving at different redshift. However, it is impossible to reach this kind of conclusion using the current coarse spectral resolution. In the near future, the X-ray bolometer on board Astro-H (Takahashi et al. 2014) will provide ICM spectra with very high resolution along the line of sight, so that the velocity regime probed in the ICM will reach few ×100 km s−1 as opposed to the current limit of ∼2000 km s−1. However, the lack of adequate spatial resolution will still make the Chandra data a crucial complementary source of information for exploring the ICM dynamics in the coming years.

7. CONCLUSIONS

We search for bulk motions in the ICM in a sample of six clusters at 0.1 < z < 0.3 classified as massive mergers on the basis or their X-ray morphology, diffuse radio emission, or dynamical optical substructure, plus two relaxed clusters. By performing spatially resolved spectral analyses on Chandra CCD data, we obtain the distribution of the ICM redshift zX across their ICM after a careful and extremely conservative evaluation of the total error on each zX value. We accurately investigate the distribution of zX, the zX map, and finally obtain a robust estimate of bulk motions in the ICM. Considering the global distribution of zX across the ICM, our analysis shows that for two of the six merger clusters, there is significant evidence of bulk motions at more than 3σ at the level of vBM = 1400 ± 300 and 4600 ± 1100 km s−1 in A2142 and A115, respectively. We also obtain measurements at lower significance (less than 2σ) of vBM = 1700 ± 840 and 1800 ± 900 km s−1 in A2034 and A520, respectively. Finally, for two clusters (1RXS J0603.3+4214 and A2146), the global analysis of the zX distribution is not significantly different from that of two relaxed clusters (A1689 and A1835), where we do not find any signs of bulk motion, as expected.

Eventually, we consider specific regions by combining the zX distribution with the redshift maps, and we are able to identify regions that show relevant bulk motion velocities in all of the merger clusters. We identify local velocity differences between ICM regions ranging from 5900 to 8600 km s−1, with typical uncertanties of 30%–50%. Whenever possible, we use the spatial distribution of zX to infer a qualitative description of the dynamical status of the clusters, tentatively classifying them as pre-merger, ongoing, or post-merger. Our results are summarized in Table 4.

We conclude that Chandra data can be successfully used to detect ICM bulk motions in massive merging clusters at 0.1 < z < 0.3. Unfortunately, the CCD coarse spectral resolution makes it impossible to provide a meaningful reconstruction of the velocity field of the ICM in single clusters. However, as soon as high-resolution, spatially resolved X-ray spectroscopy is available, the ICM dynamics will become a key diagnostic for understanding the dynamics of whole clusters, the physics of the ICM itself, for achieving accurate measurement of the total mass, and, utimately, for performing accurate cosmological tests based on cluster abundances. We are currently planning a refined strategy to carefully relax some of our assumptions, and may reach a lower velocity regime with current CCD data. The hard limit imposed by the spectral resolution can be partially overcome by near-future X-ray bolometers. However, due to the low angular resolution of X-ray bolometers, Chandra data will be a key complementary probe for ICM dynamics in the coming years.

This work was supported by the Chinese Ministry of Science and Technology National Basic Science Program (Project 973) under grants No.2012CB821804 and 2014CB845806, the Strategic Priority Research Program "The Emergence of Cosmological Structure" of the Chinese Academy of Sciences (No. XDB09000000), the National Natural Science Foundation of China under grants Nos. 11373014, 11073005, and 11403002, and the Fundamental Research Funds for the Central Universities and Scientific Research Foundation of Beijing Normal University. PT is supported by the Recruitment Program of High-end Foreign Experts and he gratefully acknowledges the hospitality of Beijing Normal University. AL is supported by NSF grant AST-1413056.

Footnotes

Please wait… references are loading.
10.3847/0004-637X/821/1/29