Characterization of diffuse X-ray emission from IGR~J17448-3232: an implication of a line of sight merging activity
Shoko Watanabe, Shigeo Yamauchi, Kumiko K. Nobukawa, and Hiroki, Akamatsu

TL;DR
This study analyzes the diffuse X-ray emission from galaxy cluster IGR J17448-3232, revealing non-ionization equilibrium conditions likely caused by a line-of-sight merging event, with detailed spectral modeling of the intracluster medium.
Contribution
It provides the first detailed spectral analysis indicating NEI plasma in IGR J17448-3232 and links this to a merging activity along the line of sight.
Findings
ICM in the core has high electron temperature (~13-15 keV)
Outer region shows lower temperature (5-8 keV)
Evidence suggests a line-of-sight merging event
Abstract
Results of the spectral analysis for the galaxy cluster IGR J17448-3232 are presented. The intracluster medium (ICM) in the central region (r<300", 320 kpc) has a high electron temperature plasma of kTe~13-15 keV and an ionization temperature estimated from an intensity ratio of Fe XXVI Ly alpha/Fe XXV He alpha lines is lower than the electron temperature, which suggests that the ICM is in the non-ionization equilibrium (NEI) state. The spectrum in the central region can be also fitted with a two-component model: a two-temperature plasma model in a collisional ionization equilibrium (CIE) with temperatures of 7.9 keV and >34 keV or a CIE+power law model with a temperature of 9.4 keV and a photon index of 1.1. The two component models can represent the intensity ratio of Fe XXVI Ly alpha/Fe XXV He alpha lines. On the other hand, the spectrum in the outer region (r>300'') can be explained…
| Component | Parameter | Value |
| Abs1 | (cm-2) | 5.6 (fixed) |
| FE | (keV) | 0.09 (fixed) |
| Normalization∗ | (2.0)10-2 | |
| (keV) | 0.59 (fixed) | |
| Normalization∗ | (2.9) | |
| Abundance † (solar) | 0.05 (fixed) | |
| Abs2 | (cm-2) | (1.3) |
| LP | (keV) | 0.88 |
| Normalization∗ | (1.20.3)10-5 | |
| Abundance† (solar) | 1 (fixed) | |
| HP | (keV) | 4.90.9 |
| Normalization∗ | (7.31.6) | |
| Abundance† (solar) | =LP | |
| RC | 2.13 (fixed) | |
| (1.60.6)10-8 | ||
| Equivalent width (eV) | 457 (fixed) | |
| Abs3 | (cm-2) | =2 |
| CXB | 1.412 (fixed) | |
| Normalization§ | 8.1710-7 (fixed) | |
| /d.o.f. | 296/215 |
| Parameter | Value | ||||
|---|---|---|---|---|---|
| 0′′–100′′ | 100′′–200′′ | 200′′–300′′ | 300′′–500′′ | 500′′–700′′ | |
| ( cm-2) | 2.06 | 2.060.04 | 1.93 | 2.1 | 2.5 |
| (keV) | 13.3 | 13.6 | 15.0 | 8.2 | 5.0 |
| Energy (keV) | 6.330.03 | 6.34 | 6.340.03 | 6.31 | 6.33 (fixed) |
| Intensity (10-5 photons s-1 cm-2) | 0.9 | 1.8 | 1.30.4 | 2.90.8 | 1.7 |
| Energy (keV) | 6.60 | 6.62 | 6.63 | 6.60 | 6.61 (fixed) |
| Intensity (10-5 photons s-1 cm-2) | 0.60.3 | 0.9 | 0.4 | 1.20.7 | 1.5 |
| /d.o.f. | 538/490 | 1050/970 | 468/450 | 186/211 | 197/146 |
| Parameter | Value | |||
|---|---|---|---|---|
| Model | CIE | NEI | 2 CIE | CIEpower law |
| ( cm-2) | 2.06 | 1.960.05 | 1.97 | 1.97 |
| Normalization∗ | 3.130.06 | 3.170.06 | — | — |
| (keV) | 11.8 | 14.1 | — | — |
| (keV) | — | 5.0 (fixed) | — | — |
| (1011 cm-3 s) | — | 2.5 | — | — |
| Normalization | — | — | 1.90.3 | 2.5 |
| (keV) | — | — | 7.9 | 9.4 |
| Normalization | — | — | 1.6 | — |
| (keV) | — | — | 34 | — |
| Redshift | 0.055 (fixed) | 0.055 (fixed) | 0.055 (fixed) | 0.055 (fixed) |
| Abundance† | 0.330.04 | 0.190.04 | 0.360.05 | 0.33 |
| Normalization | — | — | — | 1.010-3 |
| — | — | — | 1.1 | |
| /d.o.f. | 1125/1120 | 1083/1119 | 1098/1118 | 1102/1118 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
\draft\Received\Accepted
\KeyWords
galaxies: clusters: individual (IGR J174483232) — galaxies: clusters: intracluster medium — X-rays: galaxies: clusters
Characterization of diffuse X-ray emission from IGR J174483232: an implication of a line of sight merging activity
Shoko Watanabe11affiliation: Faculty of Science, Nara Women’s University, Kitauoyanishimachi, Nara 630-8506
Shigeo Yamauchi11affiliationmark:
Kumiko K. Nobukawa11affiliationmark:
and Hiroki Akamatsu22affiliation: SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands
Abstract
Results of the spectral analysis for the galaxy cluster IGR J174483232 are presented. The intracluster medium (ICM) in the central region (300*′′, 320 kpc) has a high electron temperature plasma of kT_{\rm e}$$\sim13–15 keV and an ionization temperature estimated from an intensity ratio of Fe\emissiontypeXXVI Ly/Fe\emissiontypeXXV He lines is lower than the electron temperature, which suggests that the ICM is in the non-ionization equilibrium (NEI) state. The spectrum in the central region can be also fitted with a two-component model: a two-temperature plasma model in a collisional ionization equilibrium (CIE) with temperatures of 7.9 keV and 34 keV or a CIEpower law model with a temperature of 9.4 keV and a photon index of 1.1. The two component models can represent the intensity ratio of Fe\emissiontypeXXVI Ly/Fe\emissiontypeXXV He lines. On the other hand, the spectrum in the outer region (300′′*) can be explained by a single CIE plasma model with a temperature of 5–8 keV. Based on the spectral feature and its circular structure, we propose that the NEI plasma was produced by merging along the line-of-sight direction.
1 Introduction
Galaxy clusters, consisting of galaxies, thin hot plasma, and dark matter, are the largest virialized objects in the Universe. The galaxy clusters have grown primarily via merging surrounding subclusters and accretion flow from large-scale structures. Some evidence of the past merging is found in their spatial and temperature structures: the shock wave produced by the merging process heats intracluster medium (ICM), and hence makes a very hot plasma (e.g., 15 keV for 1E 065756: [Markevitch et al. (2002)]), while the ICM exhibits a non-circular structure. From dynamical (equilibration) time scale (109 yr) and their low-density (– cm*-3*) nature in the central region, the ICM would be in a collisional ionization equilibrium (CIE). In fact, the CIE model well represents spectra of the ICM.
Theoretical researches predict that the ICM is in a non-equilibrium ionization (NEI) state after a merging process (e.g., [Takizawa (1999), Akahori & Yoshikawa (2010)]). Fujita et al. (2008) examined an ionization temperature using an intensity ratio of Fe\emissiontypeXXVI Ly/Fe\emissiontypeXXV He lines of the Ophiuchus galaxy cluster, one of the hottest known galaxy clusters, but found that the ICM has reached the ionization equilibrium state. Recently, Inoue et al. (2016) reported that an ionization temperature derived from the intensity ratio is below an electron temperature in the hottest part of A754, which suggests that the ICM is in the NEI state. There are also some arguments that no significant evidence of the NEI plasma is found even in merging clusters (e.g., Russell et al. (2012)). Thus, it is now debatable whether galaxy clusters contain NEI plasmas or not.
IGR J174483232 is a source located near to the Galactic plane [(, )(\timeform356.D8, \timeform-1.D7)], discovered by INTEGRAL (Bird et al. 2007, 2010). Chandra observed the field including IGR J174483232 and found an extended source, CXOU J174453.4323254, which is a counterpart of IGR J174483232 (Tomsick et al., 2009). Follow up observations with XMM-Newton revealed that IGR J174483232 is a galaxy cluster hidden behind the Galactic bulge (Barrière et al., 2015). The redshift was estimated to be =0.0550.001 from the X-ray spectroscopy. The X-ray emission exhibits roughly a circular structure and the radial profile is well represented by an isothermal model with a scale radius of \timeform3’.27 and of 0.56 (Barrière et al., 2015). The temperature is quite high (10 keV) in the inner parts, and gradually decrease to =4 keV at the 10*′–13′* annulus region as is apart from the center (Barrière et al., 2015).
The existence of the high temperature plasma may suggest that IGR J174483232 is in the post merger phase. As mentioned above, IGR J174483232 has several striking features including a large core radius, and spherical symmetric X-ray distribution, which is non-trivial in the merging clusters. To understand the orgin of these features, further X-ray analysis is required. Thus, we analyzed the XMM-Newton data and found a sign that an ionization temperature derived from an intensity ratio of Fe\emissiontypeXXVI Ly/Fe\emissiontypeXXV He lines does not match with the electron temperature. In this paper, we report results of spectral analysis. We utilize cosmological parameters of =0.3, =0.7, and H0=70 km s*-1* Mpc*-1*. Adopting a redshift of =0.055 (Barrière et al., 2015), 1*′*=64 kpc. The quoted errors are at the 90% confidence level unless otherwise mentioned.
2 Observation and Data Reduction
Observation of the galaxy cluster IGR J174483232 was carried out with the European Photon Imaging Camera (EPIC) on board XMM-Newton (OBSID 0672260101). The EPIC system is composed of two different detectors; two MOS cameras (MOS 1 and MOS 2: Turner et al. (2001)) and one pn camera (Strder et al., 2001). The MOS and pn cameras were operated in the full-frame mode. The spatial resolution of the XMM-Newton X-ray telescopes is 15*′′*.
Data reduction and analysis were made using the ESAS (Extended Source Analysis Software) package (Snowden et al., 2008) version 15.0.0 and HEASOFT version 6.24. We excluded the data including background flare events. The resultant exposure times are 16.8, 25.2, and 26.5 ks for pn, MOS 1, and MOS 2, respectively. In the present observation, MOS 1 CCD6 and MOS 2 CCD5 are out-of-function. Within 300*′′* region, pn, MOS 1 and MOS 2 are available, but for the analysis of the spectra in the outer region (300*′′*), only the pn data are used (see figure 1).
3 Analysis and Results
3.1 Spatical distribuion
Figure 1 is an X-ray image of IGR J174483232 in the 0.5–10.0 keV energy band. Contributions of point sources which have been reported in Barrière et al. (2015) are excluded. In order to maximize photon statistics, the data of MOS 1, MOS 2, and pn are co-added. It clearly shows a single peak structure with no substructures. The peak position is estimated to be (, )J2000.0=(\timeform266.D20888, \timeform-32.D56344).
We made a radial profile and confirmed that it is well represented by an isothermal model with a core radius, , of \timeform3.’0\timeform0.’5 and of 0.550.09. The resultant is consistent with those of relaxed clusters, on the other hand, shows somewhat higher value than those of relaxed clusters (e.g., Akahori & Masai (2005)).
3.2 Spectral analysis
3.2.1 Estimation of the sky background
Since IGR J174483232 is located near to the Galactic plane, we take account of the Galactic diffuse X-ray emission (GDXE) together with the Cosmic X-ray background (CXB) as a sky background. Thus, we at first estimate the sky background spectrum using data obtained on (, )=(\timeform356.D4, \timeform-1.D5) (the separation angle from IGR J174483232 of \timeform0.D4) with the X-ray Imaging Spectrometer (XIS: Koyama et al. (2007)) onboard Suzaku (Mitsuda et al., 2007) (Obs. ID 505084010).
Uchiyama et al. (2013) reported that the GDXE spectra are well represented by a model consisting of a foreground emission (FE), low- and high-temperature plasma components (LP and HP, respectively), and a reflection component (RC). The FE is composed of two components, 0.09 keV and 0.59 keV plasmas. The RC is composed of an absorbed power-law function (PL) and K (6.400 keV) and K (7.058 keV) lines from neutral Fe. The sky background model is expressed as
[TABLE]
where TP and ABS show a thermal plasma model (vapec in XSPEC) and photoelectric absorption (phabs in XSPEC), respectively. The K line intensity was fixed to be 0.125 K line intensity and the equivalent width of K line was fixed to be 457 eV (Uchiyama et al., 2013). The metal abundances, electron temperatures, and the value for FE were fixed to those in Uchiyama et al. (2013) and the spectral parameters of the CXB were fixed to the values in Kushino et al. (2002). The value of ABS3 is assumed to be twice of that of ABS2. The abundance tables were taken from Anders & Grevesse (1989). The cross sections of photoelectric absorption were taken from Balucinska-Church & McCammon (1992).
The best-fit spectral parameters are listed in table 1, while the best-fit model is shown in figure 2. This model gave a reduced (/d.o.f.) value of 296/214=1.38. There are possible systematic errors, such as an energy scale, atomic data of Fe L-lines, and small calibration errors near the Si\emissiontypeI K-edge energy. Taking account of these systematic errors, we use this model as a reasonable approximation.
3.2.2 Spectrum of IGR J174483232
In order to examine the spatial variation, spectra are extracted from 5 annular regions of 0*′′–100′′, 100′′–200′′, 200′′–300′′, 300′′–500′′, and 500′′–700′′. Point sources with 210-14* erg s*-1* cm*-2* in the 0.4–2.3 keV energy band (reported in Barrière et al. (2015)) are excluded. The source regions are displayed in figure 1.
The spectra are fitted with a model consisting of the non-X-ray background (NXB), GDXE, CXB, and emission from IGR J174483232. According to Snowden et al. (2008), the NXB is modeled as a PL function representing the soft proton contamination (a time variable flare component and a quiescent continuum component) plus Gaussians representing instrumental fluorescence lines. The normalization and index of the PL function and normalizations of the Gaussians were set to be free. Since the intensities of the LP, HP, and RC components of the GDXE depend on the Galactic latitude, the intensities are adjusted to those at the position of IGR J174483232 according to Yamauchi et al. (2016). Taking account of the scale heights of the 6.4 and 6.7 keV lines, we scaled down the intensities of RC and thermal components (LP and HP) by 0.80 and 0.89, respectively. On the other hand, since the spatial distribution of the FE component of the GDXE has not been well studied, there is large ambiguity. Thus, we used data above 1.2 keV where the contribution of the FE component is low. The spectral parameters of the CXB were fixed to the values in Kushino et al. (2002). The contribution of the GDXE and the CXB was in the range from 3 % (0*′′–300′′* region) to 27 % (500*′′–700′′* region) of the total count rates.
The spectrum of IGR J174483232 clearly exhibits Fe\emissiontypeXXV He and Fe\emissiontypeXXVI Ly lines, as is shown in figure 3. In order to estimate an ionization temperature, we simultaneously fitted the MOS and pn spectra in the 1.2–11 keV band with a model of thermal bremsstrahlung (zbremss in XSPEC) plus Gaussian lines with fixing the redshift of zbremss to 0.055 (Barrière et al., 2015). Free parameters are normalization and electron temperature of zbremss, a center energy and an intensity of Gaussian lines, and a column density. The best-fit parameters are listed in table 2 and the best-fit model is plotted in figure 3. Some residuals due to the instrumental line were found at 8 keV for the pn, which does not affect the results. The other instrumental lines of the MOS and the pn were well reproduced.
Figure 4 displays a correlation plot between the electron temperature and the intensity ratio of Fe\emissiontypeXXVI Ly line/Fe\emissiontypeXXV He lines. The ratios in the inner regions (300*′′) are 0.3–0.7, which corresponds to that of 8–10 keV plasma in the CIE state, while the temperature estimated from the continuum shape is 13–15 keV. On the other hand, the ratios in the outer region (300′′*) are consistent with that in the CIE state. In order to check influence of the background estimation on the intensity ratio of the lines and the electron temperature, we varied the intensity of the sky background and the detector background by 10 % and then obtained the consistent results within the statistical errors.
As shown in figure 4, the spectra within 300*′′* exhibit the NEI feature. We fitted the spectrum extracted from a circle with a radius of 300*′′* with a single temperature CIE model and found that the CIE model does not represent the ratio of Fe\emissiontypeXXVI Ly line/Fe\emissiontypeXXV He line well. The spectrum and the best-fit model is shown in figure 5a, while the best-fit parameters are listed in table 3 (/d.o.f.=1125/1120=1.004). Thus, we next examined an NEI plasma model (vrnei in XSPEC) for the spectrum within 300*′′. The initial electron temperature is assumed to be 5 keV that is in the outermost region (see table 2). The NEI model well represented the intensity ratio of Fe\emissiontypeXXVI Ly/Fe\emissiontypeXXV He lines. The best-fit parameters are listed in table 3 (/d.o.f.=1083/1119=0.968) and the best-fit model is plotted in figure 5b. The value is estimated to be (2.5)1011* cm*-3* s.
Although the single-temperature NEI model represents the spectrum well, a multi-component model might be another possibility. We also carried out a spectral fitting using a two-temperature CIE model (2CIE model). The abundances are linked between the low and the high temperature components. The model gave an acceptable fit (/d.o.f.=1098/1118=0.982). The best-fit parameters are listed in table 3 and the best-fit model is plotted in figure 5c. The obtained temperatures are 7.9 keV and 34 keV. In the case that the abundances are unlinked, we found that the abundance of the high temperature component cannot be constrained ( solar) and the improvement is low ( of 2).
We also applied a CIEPL model and found that the model also gave an acceptable fit (/d.o.f.=1102/1118=0.986). The best-fit parameters are listed in table 3 and the best-fit model is plotted in figure 5d. The best-fit temperature and photon index are 9.4 keV and 1.1, respectively.
We carried out the same model fit using the other abundance table in Lodders (2003) and obtained consistent results within given statistical errors.
4 Discussion
IGR J174483232 roughly has a circular structure with no substructure (see figure 1). When a single-temperature plasma is assumed, the ICM exhibits a high electron temperature of 13–15 keV in the central region and a low electron temperature of 5–8 keV in the outer region and has no evidence of a cool core that found in Ophiuchus cluster (Fujita et al., 2008). We found a sign that the intensity ratio of Fe\emissiontypeXXVI Ly/Fe\emissiontypeXXV He lines in the central region is lower than that of the CIE plasma with a temperature of 13–15 keV. The spectrum of 300*′′* can be fitted with not only the single-temperature NEI model but also the two-component model (see table 3).
When the 2CIE model is applied, the electron temperatures are estimated to be 7.9 keV and 34 keV (table 3). In the case of the CIEPL model, the electron temperature of the thermal component and the photon index of the PL model are 9.4 keV and 1.1, respectively. The temperatures of 8–10 keV are consistent with those derived from the Fe\emissiontypeXXVI Ly/Fe\emissiontypeXXV He line ratio (see figures 5c and 5d).
Using the spatial distribution derived from the isothermal model fit, the bolometric luminosities within a cluster radius, (a radius with a mean overdensity of 500 times the critical density of the Universe), of the 8 keV (2CIE) or 9 keV (CIEPL) component are estimated to be 81044 erg s*-1* or 11045 erg s*-1*, respectively, which are lower than those expected from the – relation (e.g., Reichert et al. (2011)). As for the hard components, the electron temperature of the hard component from the 2CIE model fit is higher than typical values of the relaxed (e.g., Reichert et al. (2011)) and merging galaxy clusters (e.g., 1E 065756: Markevitch et al. (2002); RX J1347.51145: Ota et al. (2008)), while the photon index of 1.1 from the CIEPL model fit is quite hard. In addition, the 2–10 keV band flux of the hard component corresponds to about a half (2CIE) or about one-fourth (CIEPL) of the total flux. Such a high flux of the high-temperature or non-thermal component relative to the thermal one has not been reported so far. Furthermore, the NEI model gives a better fit than the 2CIE or CIEPL model although all the models are statistically acceptable. Based on these results, we argue that the NEI plasma model in the central region is likely rather than the 2CIE or the CIEPL model.
The NEI plasma in the central region would be produced by shock-heating via a merger in the past. Using the obtained plasma parameters, we estimate an elapsed time from the shock-heating. A volume emission measure of 300*′′* is calculated to be 21067 cm*-3*. Assuming a cylindrical structure with a radius of 300*′′* and a height of 2, the electron density, , is estimated to be 10*-3* cm*-3*. Hence, we can calculate an elapsed time to be t_{\rm elapse}$$\sim$$n_{\rm e}t/n_{\rm e}$$\sim8106 yr, which suggests the merging event would be occurred 8106 yr ago. The crossing time of the cluster, , by a sound velocity, , is estimated to be =2/c_{\rm s}$$\sim109 yr. Since is much longer than , it is unlikely that the merger event was occurred on the celestial plane. IGR J174483232 has a circular structure without substructures and only the inner part has the NEI plasma with a high electron temperature. Based on these facts, we speculate that the merging event of IGR J174483232 is occurred along a line-of-sight direction. In order to check the scenario, further observational and theoretical studies are encouraged.
5 Conclusion
We analyzed the XMM-Newton data of the galaxy cluster IGR J174483232. The results are summarized as follows.
- •
The X-ray emission has a circular structure with no substructure.
- •
The ICM in the central region exhibits high electron temperature of 13–15 keV, while that in the outer region the temperature is 5–8 keV.
- •
An ionization temperature is estimated from an intensity ratio of Fe\emissiontypeXXVI Ly/Fe\emissiontypeXXV He lines. The value in the inner region is lower than 13 keV, while that in the outer region is consistent with those in the CIE state. These suggest that the ICM in the inner region is in the NEI state.
- •
The spectrum in the central region can be also represented by a two-component model: a two-temperature CIE model with temperatures of 7.9 keV and 34 keV or a CIEpower law model with a temperature of 9.4 keV and a photon index of 1.1. The high-temperature component or the power law component is quite hard.
- •
The obtained results support an idea that IGR J174483232 is a merging galaxy cluster. We speculate that the merging event of IGR J174483232 is occurred along the line-of-sight direction.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Akahori & Masai (2005) Akahori, T., & Masai, K. 2005, PASJ, 57, 419
- 2Akahori & Yoshikawa (2010) Akahori, T., & Yoshikawa, K. 2010, PASJ, 62, 335
- 3Anders & Grevesse (1989) Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197
- 4Balucinska-Church & Mc Cammon (1992) Balucinska-Church, M., & Mc Cammon, D. 1992, Ap J, 400, 699
- 5Barrière et al. (2015) Barrière, N. M., Tomsick, J. A., Wik, D. R., Chaty, S., & Rodoriguez, J. 2015, Ap J, 799, 24
- 6Bird et al. (2007) Bird, A. J., et al. 2007, Ap JS, 170, 175
- 7Bird et al. (2010) Bird, A. J., et al. 2010, Ap JS, 186, 1
- 8Fujita et al. (2008) Fujita, Y., et al. 2008, PASJ, 60, 1133
