Issue |
A&A
Volume 514, May 2010
|
|
---|---|---|
Article Number | L2 | |
Number of page(s) | 4 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/200913713 | |
Published online | 18 May 2010 |
LETTER TO THE EDITOR
Spatial identification of the overionized plasma in W49B
M. Miceli1,2 - F. Bocchino2 - A. Decourchelle3 - J. Ballet3 - F. Reale1,2
1 - Dipartimento di Scienze Fisiche ed Astronomiche, Sezione di
Astronomia, Università di Palermo, Piazza del Parlamento 1, 90134 Palermo,
Italy
2 - INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1,
90134 Palermo, Italy
3 - DSM/DAPNIA/Service d'Astrophysique, AIM-UMR 7158, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France
Received 20 November 2009 / Accepted 9 April 2010
Abstract
Context. Recent Suzaku X-ray observations of the
ejecta-dominated supernova remnant W49B have shown that there is a
clear indication in the global spectrum for overionized plasma, whose
physical origin is still under debate.
Aims. In order to ascertain the physical origin of this rapidly
cooling plasma, we focus on its spatial localization within the X-ray
emitting ejecta.
Methods. We confirm a saw-edged excess (interpreted as a strong
radiative recombination continuum) in the global spectrum above
8 keV, which emerges above the ionization-equilibrium model. We
produce a hardness-ratio map to determine where the plasma is
overionized and perform a spectral analysis of the regions with and
without strong overionization.
Results. We find that the overionized plasma is localized in the
center of the remnant and in its western jet, while it is not detected
in the bright eastern jet, where the expansion of the ejecta is
hampered by their interaction with a dense interstellar cloud.
Conclusions. The location of overionized plasma suggests that
the inner ejecta are rapidly cooling by expansion, unlike the outer
ejecta, for which expansion is hampered by interstellar clouds seen in H2.
Key words: X-rays: ISM - ISM: supernova remnants - ISM: individual object: W49B
1 Introduction
X-ray observations of young supernova remnants (SNRs) allow us to probe the physical and chemical conditions of the shock-heated ejecta and to study the physical processes involved in their interaction with the circumstellar ambient medium.
The source W49B is one of the brightest ejecta-dominated SNRs observed
in X-rays, where it shows a jet-like morphology. Its global X-ray
spectrum, characterized by intense emission lines from He-like and
H-like ions of overabundant metals (Si, S, Ar, Ca, Fe, and also Cr
and Mn) indicates a plasma at collisional ionization equilibrium
(CIE, see for example Fujimoto et al. 1995; Hwang et al. 2000).
These results are confirmed by spatially resolved spectroscopy with XMM-Newton (Miceli et al. 2006, hereafter M06), which reveals a significant Ni overabundance (
)
in the center. Miceli et al. (2006)
found similar chemical composition and temperatures of the ejecta in
the central and eastern regions, while in the western part the
abundances (in particular for Fe) and the temperatures are lower.
Both the comparison of the observed ejecta abundances with the yields
predicted by explosive nucleosynthesis models (performed by M06) and
the multipole expansion spatial analysis (performed by Lopez et al. 2009a) concur in indicating a core-collapse origin for W49B.
The eastern border of the remnant shows bright radio emission (Lacey et al. 2001) and is spatially coincident with a shocked molecular H2
cloud that hampers the expansion of the ejecta in the east direction
(being about three order of magnitude denser than the ejecta, see Keohane et al. 2007
and M06), thus distorting the eastern jet southward. Molecular hydrogen
and [Fe II] and radio emission have been observed in the
south-western edge of W49B, while in the center the [Fe II] emission
reveals a barrel-like structure ``surrounding'' the jet (Keohane et al. 2007). Miceli et al. (2008) have shown that W49B seems to be the result of an aspherical jet-like supernova explosion with an explosion energy of
erg and a mass of the shocked ejecta of
6
.
Kawasaki et al. (2005) analyzed the global spectrum of W49B and measured the intensity ratio of the H-like to the He-like K
lines of Ar and Ca, finding that the ionization temperature Tz is higher than the electron temperature,
,
thus claiming overionization in W49B.
However, M06 found
for Ca and no evidence of overionization for Ar and also pointed out
that the estimate of the electron temperature derived from the global
spectrum is not reliable, given that it is significantly lower than
that obtained from the spatially resolved spectral analysis of
homogeneous regions. Moreover, M06 performed the same analysis on the
uniform central region of W49B. Finding no evidence for overionization
in the spectrum, they concluded that the temperatures of the Ar and Ca
ions are consistent with
.
Recently, Ozawa et al. (2009) (hereafter O09) analyzed the Suzaku
spectrum extracted from the whole remnant and revealed a saw-edged bump
above 8 keV. They demonstrated that the bump is associated with a
strong radiative recombination continuum of iron and produced an
accurate spectral model to describe both the recombination continuum
and the lines and to obtain diagnostics for the overionized ejecta
(similar radiative recombination features have been also observed by Yamaguchi et al. 2009 in IC 443, but for Si and S). O09 derived that for Fe the ionization temperature is
keV, while
is
1.5 keV.
Moreover, O09 estimated that the overionized plasma has the same
emission measure as the bremsstrahlung-emitting plasma, thus suggesting
a common origin for the two components. These results strongly indicate
that the Fe ejecta in W49B are cooling so fast that they are
overionized. Nevertheless, there are still a few important open issues.
First of all the electron temperature derived from the global
spectrum by O09 is systematically lower than that observed in the
uniform spectral regions discussed in M06, where
ranges between
1.8 keV and
3 keV, and in the accurately selected regions presented in Lopez et al. (2009b) (based on Chandra data), where
keV. This discrepancy in
may arise because the global spectrum originates from physically
non-uniform regions of W49B with different temperatures and abundances.
Secondly, it is important to verify if all the ejecta in W49B are
overionized. While M06 did not find overionization effects for Ar
and Ca, O09 detected overionization for Fe. Notice that Lopez et al. (2009b)
have shown that the iron morphology is very distinct from that of other
elements, being more asymmetric and more segregated (and localized in
the central and eastern parts of the remnant), while Ar and Ca appear
well mixed and more isotropic. So in principle it is possible that the
physical conditions in the Fe-rich ejecta are different from those of
the other ejecta.
Finally, a spatial localization of the overionized plasma may
contribute to ascertain the physical origin of its rapid cooling, which
is still not understood.
Here we present the analysis of archive XMM-Newton observations of W49B, specifically devoted to constrain the spatial distribution of the overionized plasma. The data analysis procedure is shown in Sect. 2 and our results are presented and discussed in Sects. 3 and 4, respectively.
2 Data processing
The XMM-Newton data analyzed here consist of the two observations 0084100401 and 0084100501 (PI A. Decourchelle) already presented in detail in M06 (see their Sect. 2 and Table 1). We re-processed the data using the Science Analysis System and followed the same procedure and criteria described in M06. The images are a superposition of the pn images of the two observations and are background-subtracted, vignetting-corrected, and adaptively smoothed following the recipe described in M06. In the images, the background is subtracted by normalizing blank sky observations to the level of flux of the science observation in the 10-12 keV (12-14 keV) band for the MOS (PN) cameras.
Table 1: Best-fit parameters for the spectra extracted in region Y (black ellipse in Fig. 2) and region N (i.e. the entire remnant minus region Y).
For spectral analysis we used the same background event files used
in M06, and the background spectrum was extracted from the same region
positions on the CCD and subtracted from the source spectrum. The
background flux is always a few percents of the signal in the
8-10 keV band; as an example, for the pn camera in the spectral
region Y (see Sect. 3), it is 4% of the total signal in the 8-10 keV band and
9%
of the total flux in the 8-12 keV band. We first analyzed the
spectra extracted from the two observations independently, and verified
that they are compatible. Then we summed the spectra extracted from the
two observations for each camera. Spectra were rebinned to achieve a
signal-to-noise ratio per bin >
and ancillary response files were produced with the SAS ARFGEN task. The spectral fits were performed simultaneously on the MOS1, MOS2, and pn spectra using XSPEC (Arnaud 1996).
We adopt the best-fit model described in M06. The model consists of two
absorbed thermal components with oversolar abundances (Si, S, Ar, Ca,
Fe, and Ni) in ionization equilibrium, plus two narrow Gaussians to
model the Cr and Mn lines. We are interested in the hard X-ray
emission, so we fitted the spectra in the 4.4-12 keV band. Miceli et al. (2006)
have found small scale variations in the best-fit parameters of the
cold component and, because we will focus here on relatively extended
regions, we expect that these inhomogeneities may be an issue in the
spectral fittings. However, in the hard 4.4-12 keV band the cold
component contribution is negligible. Because it is not possible to
constrain the parameters of the cold component and the Si, S, Ar, and
Ca abundances (whose emission lines are all below 4.4 keV), we
derived the best-fit temperature and emission measure of the cold
component, the
,
and the Si, S, Ar, and Ca abundances by fitting the spectra in the
1-9 keV band and then froze them to their best-fit values when we
performed the fits in the 4.4-12 keV band. All the reported errors
are at 90% confidence.
3 Results
![]() |
Figure 1: pn (upper) and MOS (lower) global spectra of W49B in the 4.4-12 keV band, together with their best-fit model and residuals. The best-fit model consists of two absorbed thermal components with oversolar abundances in collisional ionization equilibrium. |
Open with DEXTER |
We first extracted the spectrum of the entire remnant and we modeled it
in the 4.4-12 keV band by adopting the best-fit model described in
Sect. 2. We found a best-fit electron temperature
kT2=1.72+0.03-0.02 keV, quite similar to that reported by O09.
The spectrum in Fig. 1 shows quite clearly a saw-edged excess above 8.3 keV, first associated with a radiative recombination continuum by O09. We also got a poor fit in the Fe (at
6.7 keV) and Fe+Ni (around 8 keV) lines in qualitative agreement with O09. Even though the excess is weaker in the XMM-Newton spectra than in the Suzaku data (where there are several spectral bins at more than 5 sigmas from the model, see O09), it is still significant.
![]() |
Figure 2:
Upper panel: vignetting-corrected and background-subtracted pn
map of the ratio of the count-rate in the hard (8.3-12 keV) and
bremsstrahlung (4.4-6.2 keV) bands (i.e. H-S map). The image is
adaptively smoothed to a signal-to-noise ratio 10 and the bin-size
is 12''. The white contour lines mark the EPIC count-rate (at |
Open with DEXTER |
To investigate the spatial distribution of the overionized ejecta, we
produced a hardness ratio map, i.e. a map of the ratio of the pn
count-rate in the high-energy (8.3-12 keV) band to that in the
bremsstrahlung (4.4-6.2 keV) band (upper panel of Fig. 2,
H-S map). A high ratio indicates radiative recombination continuum
and/or a hardening of the bremsstrahlung. We remark that if the
overionization were not present we would expect high values of the H-S
map in the eastern end of the remnant, where the ejecta temperature is
higher (see M06 and Lopez et al. 2009b).
Instead, in our H-S map the low ratios are to the East. To quantify the
deviations from CIE accounting for the temperature inhomogeneities, we
calculated the H-S values expected assuming CIE in the seven regions
discussed in M06, (
,
(i=1,...,7), computed using the kT2,i reported in Table 3 of M06) and compared them to those obtained in our H-S map,
.
The lower panel of Fig. 2 shows the ratios
.
If the H-S variations were due to inhomogeneities in the plasma
temperature, these ratios would be consistent with one, while values
>1 indicate that the hard flux is higher than that expected in
equilibrium of ionization. In regions 1 and 2 (at the eastern
edge of the remnant) the hardness ratio agrees with that expected in
CIE, while in regions 3-6 (at the central-western part
of W49B) it is significantly higher. In region 7 (to the
South-West) the difference is lower, but still significant. Figure 2
therefore clearly shows that the hard radiative recombination continuum
is stronger in the center and in the western region of W49B, while at
the eastern edge the H-S values are consistent with CIE. As a
cross-check we also produced an ``excess image'' (not shown here) by
subtracting from the pn image in the 8.3-12 keV band the
underlying bremsstrahlung contribution (estimated by fitting the global
spectrum in the 4.4-6.2 keV band) and obtained the same results.
The image analysis therefore indicates that the ejecta are overionized in the center and in the western region of W49B, while in the bright eastern edge of the remnant, where the ejecta are interacting with a dense molecular cloud, there is no significant excess. There are weak indications of overionization also in the south-western part of the remnant, characterized by strong radio and IR emission and faint X-rays.
As a further check we extracted a spectrum from the black ellipse shown in Fig. 2,
where we localized strong overionization effects (we call this region
``Y'') and compared it with the spectrum extracted from the
complementary part of the remnant (i. e. the entire remnant minus the overionized region, hereafter region ``N''). The two spectra are shown in Fig. 3
together with their best-fit models and residuals. To estimate the
significance of the differences in the two spectra we imposed the N
spectrum to have the same spectral bins as the lower-statistics Y
spectrum (where the signal-to-noise ratio per bin is >5). The best-fit parameters are shown in Table 1.
Figure 3 clearly
confirms that the hard X-ray excess is dominant above 8.3 keV in
region Y, while it is much less significant in region N. The
last row of Table 1 also
confirms that the excess is very large only in region Y. We remark
that the low significance of the excess in region N is not due to
statistic effects, because in this region the number of counts in the
4.4-12 keV band is much higher than in region Y (52 000 counts in region N vs.
36 000 counts
in region Y). Both the image and spectral analysis therefore
indicate that the overionization effects are mainly localized in the
central-western part of the remnant.
![]() |
Figure 3:
Upper panel: pn (green) and MOS (black and red) spectra of region Y (black ellipse in Fig. 2) together with their best-fit model and residuals. Lower panel:
same as upper panel for region N (the complementary part of the
remnant). We imposed the N spectrum to have the same spectral bins as
the lower-statistics Y spectrum (where the signal-to-noise ratio per
bin is > |
Open with DEXTER |
4 Discussion
Our analysis confirms the radiative recombination continuum in the high-energy X-ray spectrum of W49B. As explained by O09, this indicates that the ejecta (at least the Fe ions) are overionized and are rapidly cooling to reach the ionization equilibrium. Nevertheless, we have shown that the bright eastern ejecta of W49B do not present tracers of overionization. It is important to stress that the eastern jet of ejecta cannot expand freely because a large and dense molecular cloud hampers the expansion in the east direction (Keohane et al. 2007). The western jet, instead, is expanding undisturbed, because it does not interact with any dense structure.
The overionization due to a rapid cooling of the ejecta is therefore
observed only where the ejecta can expand freely. This rapid cooling is
then likely associated with the adiabatic expansion of the
central-western jet. A similar scenario has been invoked by Yamaguchi et al. (2009)
to explain the overionization effects in the supernova remnant
IC 443.
We expect that the overionized Fe ejecta experienced a strong heating
in the early phases of the SNR evolution. This is quite reasonable
considering that the explosion occurred in a very complex environment,
likely shaped by winds from the progenitor stars, and with the
circumstellar material showing a barrel-like morphology (Keohane et al. 2007).
The interaction of the remnant with this environment can produce strong
reverse and reflected shocks that can efficiently heat and ionize the
ejecta. The subsequent rapid expansion of the ejecta in the center and
to the West then produced their rapid cooling. We have estimated that a
temperature
K and
s cm-3
would benecessary to have an equal ionization (to H-like Fe) and
electron-heating time. Nevertheless, an accurate analysis requires a
modeling of the hydrodynamic evolution of the plasma including the
effects of non-equilibrium of ionization (e.g. Reale & Orlando 2008).
Not all the ejecta of W49B show signatures of overionization for the Fe ions, and Fig. 2 shows that in regions with similar (and large) Fe abundances, like the eastern jet and the central regions (see M06 and Lopez et al. 2009b),
there are different states of ionization (equilibrium to the East and
overionization in the center), while the central and western regions
(that have different Fe abundances, see M06 and Lopez et al. 2009b)
are both overionized. We conclude that while we derive overionization
effects for the Fe ions, not all the Fe-rich ejecta are overionized. On
the other hand we remind that in the center of W49B no overionization
effects are visible for the Ar and Ca ions (M06). We expect therefore
different ionization temperatures and different ratios
for different species and, at least for Fe, for different regions of
the remnant. In principle it is also possible that the temperature of
the electrons responsible for the Fe recombination is different from
that of the electrons responsible for the bremsstrahlung emission in
the 4.4-6.2 keV band. However, from the global spectrum of W49B,
O09 found that both the radiative recombination and the bremsstrahlung
continua can be associated with electrons at the same temperature.
Nevertheless, a spatially resolved study appears necessary in future
works.
Another important issue is related to the determination of the value of the ion and electron temperatures in the ejecta. Ozawa et al. (2009)
have modeled both the radiative recombination continuum and line
emission by deriving the ionization temperature from the global
spectrum, finding
keV, which is higher than the electron temperature
keV. In our large regions we also found similar values of
,
as shown in Table 1.
Nevertheless, these values are significantly lower than those derived
from smaller regions by M06. For example, our region Y is approximately
the union of regions 3, 4, 5, and 6 in M06, where they find
keV, which is significantly higher than the value
keV reported in Table 1.
It is then important to analyze spectra extracted from relatively small
regions to accurately constrain the ion and electron temperature
distribution. Unfortunately, the available statistics do not allow us
such a detailed analysis, and we therefore stress the importance of
obtaining deeper observations of W49B to address this point.
References
- Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, ASP Conf. Ser., 101, 17 [Google Scholar]
- Fujimoto, R., Tanaka, Y., Inoue, H., et al. 1995, PASJ, 47, L31 [NASA ADS] [Google Scholar]
- Hwang, U., Petre, R., & Hughes, J. P. 2000, ApJ, 532, 970 [NASA ADS] [CrossRef] [Google Scholar]
- Kawasaki, M., Ozaki, M., Nagase, F., Inoue, H., & Petre, R. 2005, ApJ, 631, 935 [NASA ADS] [CrossRef] [Google Scholar]
- Keohane, J. W., Reach, W. T., Rho, J., et al. 2007, ApJ, 654, 938 [NASA ADS] [CrossRef] [Google Scholar]
- Lacey, C. K., Lazio, T. J. W., Kassim, N. E., et al. 2001, ApJ, 559, 954 [NASA ADS] [CrossRef] [Google Scholar]
- Lopez, L. A., Ramirez-Ruiz, E., Badenes, C., et al. 2009a, ApJ, 706, L106 [NASA ADS] [CrossRef] [Google Scholar]
- Lopez, L. A., Ramirez-Ruiz, E., Pooley, D. A., et al. 2009b, ApJ, 691, 875 [NASA ADS] [CrossRef] [Google Scholar]
- Miceli, M., Decourchelle, A., Ballet, J., et al. 2008, Adv. Space Res., 41, 390 [NASA ADS] [CrossRef] [Google Scholar]
- Miceli, M., Decourchelle, A., Ballet, J., et al. 2006, A&A, 453, 567 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ozawa, M., Koyama, K., Yamaguchi, H., Masai, K., & Tamagawa, T. 2009, ApJ, 706, L71 [NASA ADS] [CrossRef] [Google Scholar]
- Reale, F., & Orlando, S. 2008, ApJ, 684, 715 [NASA ADS] [CrossRef] [Google Scholar]
- Yamaguchi, H., Ozawa, M., Koyama, K., et al. 2009, ApJ, 705, L6 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Best-fit parameters for the spectra extracted in region Y (black ellipse in Fig. 2) and region N (i.e. the entire remnant minus region Y).
All Figures
![]() |
Figure 1: pn (upper) and MOS (lower) global spectra of W49B in the 4.4-12 keV band, together with their best-fit model and residuals. The best-fit model consists of two absorbed thermal components with oversolar abundances in collisional ionization equilibrium. |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Upper panel: vignetting-corrected and background-subtracted pn
map of the ratio of the count-rate in the hard (8.3-12 keV) and
bremsstrahlung (4.4-6.2 keV) bands (i.e. H-S map). The image is
adaptively smoothed to a signal-to-noise ratio 10 and the bin-size
is 12''. The white contour lines mark the EPIC count-rate (at |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Upper panel: pn (green) and MOS (black and red) spectra of region Y (black ellipse in Fig. 2) together with their best-fit model and residuals. Lower panel:
same as upper panel for region N (the complementary part of the
remnant). We imposed the N spectrum to have the same spectral bins as
the lower-statistics Y spectrum (where the signal-to-noise ratio per
bin is > |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.