A&A 376, 248-253 (2001)
DOI: 10.1051/0004-6361:20010882
F. Bocchino1,2 - A. M. Bykov3
1 -
Astrophysics Division, Space Science Dept. ESA-ESTEC,
Postbus 299, 2200 AG Noordwijk, The Netherlands
2 -
Osservatorio Astronomico di Palermo, Piazza del Parlamento 1,
90134 Palermo, Italy
3 -
A. F. Ioffe Institute for Physics and Technology,
194021 St. Petersburg, Russia
Received 11 April 2001 / Accepted 15 June 2001
Abstract
XMM-Newton observations of the X-ray feature 1SAX J0617.1+2221 in
the IC 443 supernova remnant are reported. We resolve the structure of the
nebula into a compact core with a hard spectrum of photon index
in the 2-10 keV energy range. The nebula also has
an extended (
)
X-ray halo, much
larger than the radio emission extension. The photon index softens,
following a linear scaling with distance from the centroid, similar to
other known X-ray plerions. The index range is compatible with
synchrotron burn-off models. All the observational evidence points
toward a confirmation of the plerionic nature of the nebula, as
recently suggested by a Chandra observation, but with characteristics
more similar to "non Crab-like" plerions. We discuss the implications
on the synchrotron nebula magnetic field if the >100 MeV emission
reported by CGRO EGRET is produced by the synchrotron emission.
We also constrain the thermal emission of the central object, arguing that
the surface temperature should be around 0.1 keV, although other
possible fits cannot be excluded on the base of the XMM-Newton data.
Key words: pulsars: general - shock waves - ISM: supernova remnants - ISM: individual object: IC 443
The processes of conversion of pulsar spin-down energy into high-energy emission in pulsar wind nebula are of great physical interest (e.g., Kennel & Coroniti 1984; Chevalier 2000). Modern arcsec resolution instruments such as on XMM-Newton and Chandra that are sensitive to photons up to 10 keV provide unique possibilities to study pulsar nebulae, because the non-thermal emission of the nebulae may be easily detected and studied in this band (e.g., Vela-X, Helfand et al. 2001; Pavlov et al. 2001; G21.5-0.9, Slane et al. 2000; Warwick et al. 2001; and 3C 58, Bocchino et al. 2001), even when immersed in the soft X-ray emission of the companion shell, as is often the case.
Indeed, the shell SNR IC 443, once though to be mostly thermal in the
X-ray band (Petre et al. 1988; Asaoka & Aschenbach 1994) has been discovered to emit hard
X-ray emission by Wang et al. (1992). ASCA Gas Imaging Scintillator (GIS)
observations by
Keohane et al. (1997) discovered the localized character of the hard X-ray
emission and its non-thermal nature. They concluded that most of the
2-10 keV photons came from an isolated emitting feature and
from the South East elongated ridge of hard emission. Preite-Martinez et al. (1999) and
Bocchino & Bykov (2000) reported a hard component detected with the
Phoswich Detector System (PDS) on BeppoSAX
and two compact X-ray sources corresponding to the ASCA sources
detected with the BeppoSAX Medium-Energy Concentrator Spectrometer
(MECS) (1SAX J0617.1+2221 and 1SAX J0618.0+2227).
Very recently, 1SAX J0617.1+2221 has been observed by Chandra as reported
by Olbert et al. (2001), who also show a VLA observation at 1.46, 4.86 and
8.46 GHz, and a polarization measurement. They argue that the hard
radio spectral index, the amount of polarization and the overall X-ray and
radio morphology strongly suggest that the source is a plerion nebula
with a point source in it, whose characteristic cometary shape is due
to supersonic motion of the neutron star. However, the limited
counting statistics of the 10 ks Chandra observation do not allow a
detailed spectral study of the nebula, which is required to compare
this new plerion with current theoretical models. Moreover, IC 443 is also a
possible candidate for the CGRO EGRET -ray source 3EG
J0617+2238 (Hartman et al. 1999) having a flux above 100 MeV of
ph s-1 cm-2 with a photon index of
.
More detailed measurements of the nebula spectral properties are needed
to study the relation between 3EG J0617+2238 and the nebula.
In this paper an XMM-Newton study of the recently discovered plerion nebula in IC 433 is presented. In particular, we use the large effective area of the EPIC instrument to address the synchrotron burn-off effect in the nebula, to resolve its structure and measure the flux, and to constrain the thermal radiation of the central object in the nebula.
IC 443 was observed as part of the Cal/PV phase of the XMM-Newton Observatory (Jansen et al. 2001). Here, we have used the two
observations centered on
(J2000) performed on 2000 September 27.
Data from the two MOS (Turner et al. 2001) cameras and the PN (Strueder et al. 2001)
camera are used. The MOS and PN cameras are CCD arrays which collect
X-ray photons between 0.1 and 15 keV and have a field of view of
diameter. The pixel size is
and
for MOS
and PN, respectively, while the mirror Point Spread Function is 6
Full-width at half maximum.
The data were acquired with the medium
filter and in full image mode, and therefore the temporal resolution is
2.5 s and 73 ms for the MOS and PN, respectively. The poorer spatial
resolution of the PN is compensated for by its greater sensitivity, on the
average 20-30% more than the combined MOS cameras.
The Standard Analysis System (SAS) software used (version 5.0.1, xmmsas-20001215) takes cares of most of the required events screening. However, we have further screened the data to eliminate some residual hot pixels and occasional background enhancement during intervals of intense incident flux of soft protons. In particular, we extracted the background lightcurve at energies >10 keV and identified time intervals of unusually high count rates (typically more than 1 s-1) and removed them from subsequent analysis. We merged the the event files of the subpointings before continuing with the analysis. The total screened and unscreened exposure times are 24 and 32 ks, respectively.
We extracted PN images with a 4
pixel size and exposure images
at the same resolution. We defined two energy bands in which to
extract images, namely 0.5-2 keV (hereafter the soft band) and 3-10 keV
(hard band). Figure 1, top panel, shows the 0.5-2 keV PN
image of the South Eastern section of IC 443. Most of the emission in this
band is of thermal origin, as established by Asaoka & Aschenbach (1994) based on
ROSAT data.
Figure 1, bottom panel, shows the 3-10 keV PN image of the
same field with a logarithmic color scaling and a heavier smoothing
chosen in order to emphasize the faint diffuse emission. Most of the
thermal emission associated with IC 443 is not present in this band, and
therefore it is well suited to the study of any hard compact
X-ray sources in the vicinity of this SNR. Indeed, this image shows
several point sources, besides the plerion nebula itself. We see that
the nebula is much more extended then reported by Olbert et al. (2001) and can
be represented by an ellipse of
.
By fitting a two dimensional Gaussian to the image, we find that the
location of the plerion nebula centroid is
(J2000, therefore XMMU J061705.5+222130), with an estimated
accuracy of 5
,
due to both statistical uncertainties (at
)
and systematic errors introduced by attitude reconstruction,
and therefore in agreement with the position of the Chandra source CXOU
J061705.3+222127. Figure 2 shows the combined MOS 3-10 keV image
of the nebula core. The morphology of the faint rim matches somewhat
that of the bright core, which may suggest a common origin, in
agreement with the plerionic interpretation. Figure 2 also
shows that the radio halo at 8.46 GHz observed by Olbert et al. (2001) has
somewhat different shape, and is smaller than the X-ray full extent of
the nebula.
We also note that the source 1SAX J0618.0+2227 is resolved with XMM-Newton into two closeby sources (see Fig. 1). Moreover, there are
other hard X-ray sources in the field of view. Some of these sources
might be related to the SNR. The location of these sources is
consistent with a scenario concerning the interaction of the SNR with
molecular cloud, as proposed by Bocchino & Bykov (2000).
![]() |
Figure 2:
MOS1+2 3-10 keV image closeup of the core of the nebula, smoothed with a Gaussian of 4
![]() ![]() |
The spectral analysis of the plerion is rather complicated because of
the presence of diffuse thermal emission on large spatial scales which is
unassociated with the plerion nebula. Figure 1 shows that
the plerion is partially superimposed on a large region of soft
emission, running from the North East to the
South West, which links the plerion region with
the very bright soft region around
.
In order to cope with these problems, we adopted the following procedure. First, we identified a "background area" (Bkg in Fig. 1) which is outside the thermal emission of IC 443. This background have been corrected for vignetting and area effects prior to subtraction of the source spectra. The background correctly includes emission from G189.6+3.3, the "companion" SNR identified by Asaoka & Aschenbach (1994), which encompasses all the field of view of the XMM-Newton observation. Then, we defined a sample of the thermal region in which the plerion is partially immersed (Th2 in Fig. 1). Unfortunately, the Th2 region cannot be located outside the plerion nebula, because the spectral properties may be different in a region so far from the nebula. The thermal region was fitted with the two-temperature optically thin plasma model of Mewe et al. (1985) and Liedahl et al. (1995) with interstellar absorption of Morrison & McCammon (1983), using the appropriate PN response matrix for off-axis sources. An additional power-law component was used to model the non-thermal emission of the plerion falling inside the Th2 extraction region. The above 2T model has been successfully used by Asaoka & Aschenbach (1994) to describe the ROSAT data of this part of IC 443. A proper study of the remnant thermal emission including more appropriate models (e.g., Non-Equilibrium of Ionization) is beyond the scope of our present work.
The results obtained from the thermal regions show that the two
temperature thermal model used by Asaoka & Aschenbach (1994) provides a good
description of the data with best-fit parameters in agreement with
the findings of ROSAT (
kT1=0.18+0.10-0.03,
kT2=0.90+0.10-0.05,
,
).
We analyzed the spectrum of the compact source XMMU J061705.5+222130, extracting its
spectrum from a 6
radius region centered on the source
centroid. In order to take into account the residual thermal emission
in the fittings, we used a power-law model together with the 2T
model used in the Th2 region to fit the data. The absorption,
temperatures and normalization ratio of the two thermal component were
kept fixed to the values found in Th2. The best-fit results are
reported in Table 1. Since the best-fit normalization of
the thermal component is a factor of 2 higher than that expected by
scaling the value found in Th2, we have tested for the presence of an
additional black-body component in the spectrum of this central part of
the nebula, by adding this component to the power-law instead of the 2T
component. There is no a significant improvement in
,
but a
best-fit is found with kT=0.13 keV and a neutron star (NS) radius of 3.3 km. The corresponding values of the black-body temperature and radius
as observed at infinite distance allowed by our data are shown in
Fig. 3. At the temperature derived by Olbert et al. (2001) with Chandra, kT=0.7 keV, the upper-limit to the 1-5 keV unabsorbed flux of the BB
component is
erg cm-2 s-1, (consistent
with the Chandra detection of a flux of
erg cm-2 s-1), but
the high temperature and very low emitting area (even for heated polar
cap models) points against this region of the parameter space.
We have also investigated the presence of periodic signal by performing
an FFT of the PN events collected in the NS source region, i.e. a
circular region
centered on the centroid with a 6
radius. We did not find a
periodic signal at the 99% confidence level in the
10-4-6.5 Hz
range, with an upper limit, at the same confidence level, of 16.8% of
sinusoidal pulsated fraction.
In order to test the plerionic interpretation, we have tried to detect
the softening of the spectrum toward the outer regions of the nebula, a
well known effect which has been observed in other plerions (e.g. 3C 58,
Torii et al. 2000; Bocchino et al. 2001; G21.5-0.9 by Slane et al. 2000; Warwick et al. 2001)
and which is caused by the short
lifetime of high energy electrons compared with those with lower
energies. Given the peculiar shape of this nebula, we
extracted spectra from four elliptical annuli, shown in
Fig. 2, and fit them separately. We chose the
ellipse sizes in order to minimize the point spread function
(PSF) contamination effects. For
instance, the first region is a 12
radius circle, in which
of the PSF integral is included.
Because of the contamination by thermal emission, we have followed the
same approach used in fitting the NS region, i.e. a fixed parameter 2T
model and a power-law. Whilst this procedure does not take into account
the uncertainties associated in the Th2 fit (which however are not
large), it provides an accurate way of not rejecting useful data,
thus reducing the overall uncertainties in the spectral parameter of
the nebula. For comparison, we also show the results obtained with
fits to the hard band spectrum (E>3 keV) and a single power-law.
Name |
![]() |
![]() |
![]() ![]() |
![]() |
erg cm-2 s-1 | ||||
NS | - | 1.63+0.11-0.17 | ![]() |
27/24 |
PlCore | 1.71+0.44-0.23 | 1.63+0.11-0.10 | ![]() |
40/48 |
Pl1 | 1.95+0.35-0.18 | 1.81+0.08-0.06 |
![]() |
119/108 |
Pl2 | 1.88+0.31-0.16 | 1.88+0.07-0.06 |
![]() |
149/149 |
Pl3 | 2.70+0.31-0.26 | 2.30+0.09-0.09 |
![]() |
451/435 |
a Fit to power-law only, 3-10 keV,
cm-2.
b Fit to power-law + fixed parameters 2T components
(
,
kT1=0.2, kT2=0.9,
norm1/norm2=10.
Full bandwidth (0.5-10 keV).
c Power-law unabsorbed flux in the
2-10 keV band.
The results of the spectral analysis are shown in Table 1 and the
spectrum of the central core is shown in Fig. 4. All the fits
provide a good description of the data. In Fig. 5 we
present the spectral softening in the nebula going from the inner to
the outer regions, as well as the softening observed in other plerions,
for comparison.
Letting the absorption value vary does not yield a better
fit, therefore we argue that there is no evidence for
any absorption variation across the plerion nebula.
![]() |
Figure 5:
![]() ![]() ![]() ![]() |
The XMM-Newton observations allow us to resolve the structure of the hard X-ray
source XMMU J061705.5+222130. The observed power-law index steepening (Fig. 5)
is consistent with the synchrotron burn-off effect
(Kennel & Coroniti 1984) and remarkably similar to the effect observed in other
plerions, thus confirming the pulsar wind nebula interpretation
suggested by Olbert et al. (2001). Note that the frequency index steepening
is consistent with the standard value of 0.5 from the one-zone model by
Chevalier (2000), and that the hard spectral index of the core suggests a
low X-ray efficiency, similar to the Vela-X nebula and CTB 80 as
reported by Chevalier (2000). Following the one-zone Chevalier (2000),
away from cooling regime the
ratio is dependent on the
fraction
of energy density of the emitting region going into
magnetic fields (
), and this
quantity is rather uncertain. However, assuming
,
like in the bow-shock nebula CTB 80, and using the sum of the flux in
all the elliptical regions (
erg s-1 in the
0.2-4 keV band at 1.5 kpc), we obtain
erg s-1 (
erg s-1 with the empirical relation of
Seward & Wang 1988), similar to the value adopted by Olbert et al. (2001).
Another important characteristic of this plerion nebula is that its X-ray extent appears larger then its radio extension, an effect also observed in G21.5-0.9 (Warwick et al. 2001) and 3C 58 (Bocchino et al. 2001). This is unexpected if both X-ray and radio are just the synchrotron emission of the nebula. However, the inverse Compton and non-thermal bremsstrahlung components of the X-ray emission, as well as other thermal components related to a possible shell, may become important in the outer halo of the nebula. Unfortunately, there is no way to verify if, apart from the contamination of the IC 443 thermal emission, the spectrum of the outer nebula is purely non-thermal (like in G21.5-0.9, Warwick et al. 2001), or contains an additional thermal component due to the expansion of the nebula into the ejecta (as in 3C 58, Bocchino et al. 2001). On the other hand, such a small radio nebula may be the result of a low sensitivity of the VLA observations reported by Olbert et al. (2001). The X-ray to radio luminosity ratio, defined for instance as in Woltjer et al. (1997), is only a factor of 2 less then the Crab's and 100 more than 3C 58. But if deeper radio observations reveal a nebula comparable with, or larger then, the X-ray nebula the real value of the ratio may be reduced.
By extrapolating back the sum of the spectra of the different nebula regions in the radio regime, and comparing them with the radio spectrum reported by Olbert et al. (2001), we have found a spectral break at 100 GHz, instead of 10000 GHz as reported by Olbert et al. (2001). The difference is due to the contribution of the hard inner nebula regions, which were not properly modeled using ASCA data, on which the Olbert et al. (2001) estimation relies upon. Such a low frequency break is also found in 3C 58 and G21.5-0.9, as well as in other examples of other "non Crab-like" plerion reported by Woltjer et al. (1997).
It should be noted, that there are some differences between the IC 443 plerion and 3C 58 (and G21.5-0.9). In particular, the spectral index steepening in 3C 58 and G21.5-0.9 is higher than 0.5, the value predicted by the one-zone model of Chevalier (2000). Another key difference is the central source in the center, which is very prominent in case of the IC 443 plerion. However, on the basis of the published Chandra data on G21.5-0.9 (Slane et al. 2000), the presence of central compact sources in G21.5-0.9 cannot be excluded.
As for the NS, a better interpretation of the its temperature
constraint requires a careful analysis of the NS atmosphere effect
(e.g., Pavlov et al. 1995) and general relativitistic corrections. The effect of
the hydrogen atmosphere could reduce the effective temperature
typically by a factor 1.3 to 2 and thus could make the effective NS
radius consistent with the standard NS radius of 10 km. It is worth
noting that the NS temperature of
keV is nicely
consistent with the IC 443 age of 30000 years (Chevalier 1999), if the
standard NS cooling curve of superfluid NS is used (e.g., Yakovlev et al. 1999,
- see their Fig. 7), but a lower age cannot be excluded. The
discrepancy between our best guess for the NS surface temperature and
the temperature derived by Olbert et al. (2001) by fitting the point source
CXOU J061705.3+222127 may be due to the low value of
assumed by
them (
cm-2). We verified that this
value
is rejected at 99% confidence level by all our fits. A deeper
Chandra observation will better constrain the physical properties of the
compact source.
If the hard power-law spectrum of the nebula core region is
extrapolated up to GeV regime, it would provide a flux of 2.0
(0.4-15.2)
ph cm-2 s-1, which is consistent
with the EGRET flux of 3EG J0617+2238 (
ph cm-2 s-1). The nebula is placed outside the 95% uncertainty
circle of the 3EG J0617+2238, but the systematic uncertainties of the
EGRET circle allows us to speculate about possible association of
the two sources. If, however, the GeV regime EGRET photons indeed originate from synchrotron emission of
relativistic electrons, this would seriously constrain the nebula model
of Kennel & Coroniti (1984) where e
are accelerated by a relativistic wind
termination shock and then radiate in the downstream region. The
maximum photon energy
(measured in GeV) produced by an
accelerated electron of maximal Lorentz factor
in the
magnetic field of the emitting region
(measured in G) can be
obtained from
G. The e
synchrotron loss time,
s, must exceed the e
acceleration time by relativistic shock,
(where
is the
magnetic field in the acceleration region). Since
cannot be less
then the particle gyroperiod, one may obtain
G. Thus, the magnetic field amplification in the shock
downstream region
is required. The
amplification effect can be studied with high spatial resolution X-ray
imaging. The EGRET luminosity of
erg s-1 would not be a
negligible fraction of the spin-down power estimated from the
Seward & Wang (1988) relation. An alternative explanation for the EGRET
emission is the relativistic bremsstrahlung of radio emitting electrons
accelerated by the supernova radiative shock (Bykov et al. 2000), and future
INTEGRAL and GLAST observations will resolve that alternative, as well
as providing an accurate position for 3EG J0617+2238.
We report XMM-Newton Cal/PV observations of the hard X-ray nebula inside
the supernova remnant IC 443, recently discovered by Chandra (Olbert et al. 2001). The longer exposure time and the high XMM-Newton effective
area have allowed us to trace and to study the diffuse emission of the
nebula far beyond the extension reported by Olbert et al. (2001). We confirm
the plerionic nature of the nebula, as suggested by Olbert et al. (2001), by
virtue of the observation of the synchrotron burn-off, which is
in agreement with plerion model expectations. The
relatively hard power-law index of the core, the flat radio spectrum, the low frequency break,
and the X-ray extension much larger then the radio extension at 8.46
GHz make this plerion more similar to 3C 58 and G21.5-0.9 (two well
known examples of "non Crab-like" or "second kind" plerions,
Woltjer et al. 1997), rather than to the Crab nebula. Our data allows us to
place only upper-limit on the presence of a thermal compact X-ray
source in the nebula core. However, we argue that its temperature
is likely be
keV, significantly lower that the
temperature suggested by Olbert et al. (2001). We also argue that, if the
nebula is the counterpart of the EGRET source 3EG J0617+2238, as not
excluded by the extrapolation of its X-ray spectrum, a strong constraint
on magnetic field amplification in the nebula may be posed.
Acknowledgements
F. Bocchino acknowledges an ESA Research Fellowship. The work of A.M.B was supported by the INTAS-ESA 99-1627 grant. We thank A. N. Parmar and the anonymous referee for their useful comments on the manuscript. We also thank the ESTEC XMM-Newton staff for their suggestions about data analysis.