Issue |
A&A
Volume 502, Number 1, July IV 2009
|
|
---|---|---|
Page(s) | 131 - 137 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200809735 | |
Published online | 27 May 2009 |
Constraints on the kinematics of the 44Ti ejecta of Cassiopeia A from INTEGRAL/SPI
P. Martin1,2 - J. Knödlseder1,2 - J. Vink3 - A. Decourchelle4 - M. Renaud5
1 - Centre d'Etude Spatiale des Rayonnements (CESR), 9 avenue colonel Roche, BP44346, 31028 Toulouse Cedex 4, France
2 - Université de Toulouse (UPS), Centre National de la Recherche Scientifique (CNRS), UMR 5187, France
3 - Astronomical Institute, University Utrecht, PO Box 80000, 3508 TA Utrecht, The Netherlands
4 - Laboratoire AIM, CEA/DSM-CNRS-Université Paris Diderot, DAPNIA/Sap, CEA-Saclay, 91191 Gif-sur-Yvette, France
5 - Laboratoire APC, CNRS-UMR 7164-Université Paris Diderot, 10 rue A. Domon et L. Duquet, 75205 Paris Cedex 13, France
Received 6 March 2008 / Accepted 3 April 2009
Abstract
Context. The medium-lived 44Ti isotope is synthesised by explosive Si-burning in core-collapse supernovae. It is extremely sensitive to the dynamics of the explosion and therefore can be used to indirectly probe the explosion scenario. The young supernova remnant Cassiopeia A is to date the only source of gamma-ray lines from 44Ti decay. The emission flux has been measured by CGRO/COMPTEL, BeppoSAX/PDS and INTEGRAL/IBIS.
Aims. The high-resolution spectrometer SPI on-board the INTEGRAL satellite can provide spectrometric information about the emission. The line profiles reflect the kinematics of the 44Ti in Cassiopeia A and can thus place constraints on its nucleosynthesis and potentially on the associated explosion process.
Methods. Using 4 years of data from INTEGRAL/SPI, we have searched for the gamma-ray signatures from the decay of the 44Ti isotope. The overwhelming instrumental background noise required an accurate modelling and a solid assessment of the systematic errors in the analysis.
Results. Due to the strong variability of the instrumental background noise, it has not been possible to extract the two lines at 67.9 and 78.4 keV. Regarding the high-energy line at 1157.0 keV, no significant signal is seen in the 1140-1170 keV band, thereby suggesting that the line signal from Cassiopeia A is broadened by the Doppler effect. From our spectrum, we derive a 500 km s-1 lower limit at 2
on the expansion velocity of the 44Ti ejecta.
Conclusions. Our result does not allow us to constrain the location of 44Ti since the velocities involved throughout the remnant, derived from optical and X-ray studies, are all far above our lower limit.
Key words: gamma rays: observations - ISM: supernova remnants - ISM: individual objects: Cassiopeia A - nuclear reactions, nucleosynthesis, abundances
1 Introduction
Until very recently, Cassiopeia A (hereafter Cas A) occupied a specific place among galactic supernova remnants as being the result of the most recent galactic supernova known to date (the youngest one might now be G1.9+0.3, with an estimated age of about a century, see Reynolds et al. 2008). Based on an optical study of outer high-velocity knots, Thorstensen et al. (2001) derived an explosion date of 1671 or slightly later for Cas A. Combined with the 3.4 kpc distance estimated by Reed et al. (1995), the explosion that eventually lead to Cas A might well be the transient 6th magnitude star noticed by Flamsteed in 1680. From the light-echo scattered by interstellar dust, the event has now been identified as an SNIIb, which indicates that the progenitor was a red supergiant that had lost most of its hydrogen envelope at the moment of collapse (Krause et al. 2008).
Cas A retains interesting features that could be valuable clues in the long-standing effort to understand the supernova explosion mechanism. Indeed, growing evidence has arisen that the explosion and/or the expansion of Cas A were asymmetric and highly turbulent. X-ray images taken with Chandra have revealed a jet (and to a lesser extent a counter-jet) protruding in the north-east direction far beyond the nearly spherical forward shock (Vink 2004; Hwang et al. 2004). Another indication for an asymmetric process is the substantial mixing suggested by the Chandra X-ray images in which Fe-rich material appears, in at least two extended regions, at a greater projected radius than Si-rich material (Hwang et al. 2000; Hughes et al. 2000). On a smaller scale, the wide variety in the so-called optical Fast-Moving Knots (FMKs) spectra regarding the relative strengths of O and S lines, together with the absence of a correlation between these spectral properties and the velocity of the FMKs, also imply some turbulence or uneven burning (Chevalier & Kirshner 1979).
In addition to these peculiarities, Cas A is known as the only source of gamma-ray line emission from 44Ti decay. The unstable 44Ti isotope is a by-product of explosive Si-burning resulting from alpha-rich freeze-out (The et al. 1998; Woosley et al. 1973). The 44Ti is expected to lie deep in the ejecta and its yield is very sensitive to the early dynamics of the explosion and to the position of the mass cut (the boundary that separates the ejected material from the material that falls back onto the compact object). As such, 44Ti provides indirect evidence about the flow pattern that mediated the explosion energy and can therefore be used to probe the explosion scenario.
After a mean lifetime of 85 yr (Ahmad et al. 2006), 44Ti decays by electron capture to 44Sc and then to 44Ca, thereby emitting three de-excitation photons at 67.9, 78.4 and 1157.0 keV. Given its youth and relative proximity, Cas A constitutes an ideal target to search for the characteristic decay emission of the medium-lived 44Ti. The first detection of the 44Ti emission from Cas A was reported by Iyudin et al. (1994) from CGRO/COMPTEL observations at 1157.0 keV, with a downward revision of the flux by Dupraz et al. (1997). Subsequent detections of the low-energy lines by BeppoSax/PDS (Vink et al. 2001) and recently by INTEGRAL/IBIS/ISGRI (Renaud et al. 2006) have allowed investigators to tighten the flux value for the 44Ti emission from Cas A to
(
)
10-5 ph cm-2 s-1 (Renaud et al. 2006). With the age and distance mentioned above for Cas A and assuming that the 44Ti is not strongly ionized (which could delay its decay by electron capture, see Mochizuki et al. 1999), this flux translates into an initial 44Ti mass of about
.
While the total 44Ti flux, and hence the total 44Ti mass, is rather well constrained today, little is known about the kinematics of the 44Ti ejecta as none of the above-mentioned instruments had the capabilities required for fine spectrometry. The high-resolution spectrometer SPI on-board the INTEGRAL satellite opened the way for a determination of line profiles and thus for an estimate of the 44Ti velocity. In particular, the high-energy line at 1157.0 keV is expected to reflect the motion of the 44Ti ejecta. Such a measurement is essential to achieve a complete understanding of the nucleosynthesis of 44Ti in Cas A and is relevant to fundamental questions such as how core-collapse supernovae explode and why Cas A is the only source of 44Ti emission in the sky (The et al. 2006).
In this paper we present the INTEGRAL/SPI observations of the 44Ti decay line at 1157.0 keV and their implications in terms of velocity of the 44Ti ejecta. From this constraint, the site of the nucleosynthesis of 44Ti is then discussed in the context of all other observational data available on Cas A.
2 Data analysis
2.1 Instrument characteristics
The SPI instrument onboard the INTEGRAL gamma-ray space observatory is a high-resolution spectrometer with imaging capabilities (Roques et al. 2003; Vedrenne et al. 2003). The spectrometry of gamma radiation with energies between 20 keV and 8 MeV is performed by an array of 19 high-purity germanium detectors, with a spectral resolution of about 2 keV FWHM at 1 MeV. All gamma-ray sources, diffuse or point-like, emitting in this energy range can also be imaged indirectly thanks to a coded mask system with an angular resolution
of 2.8
and a fully-coded field of view of 16
.
For a given pointing, a coded-mask achieves a spatial modulation of the celestial signal and the reconstruction of the sky intensity distribution exploits the correlations between the recorded image and the mask pattern or a decoding array derived from it (see Fenimore & Cannon 1978). The instrument response of SPI, however, is quite complex and inversion by correlation methods is not possible. More information should be accumulated and the observation of an astrophysical object or region with SPI therefore consists of a series of pointings around the target, following a so-called dithering pattern (thus adding a temporal modulation of the astrophysical signal). The sky signal is then reconstructed through a model-fitting approach.
The photon interactions recorded by SPI are divided into two classes: single events, be they from astrophysical or instrumental origin, depositing all of their energy in a single detector, and multiple events, in which case the interactions involve several adjacent detectors (which therefore form a so-called pseudo-detector). In the work presented here, only single events (hereafter SE) and double events (hereafter ME2) were considered. Detector 2 and 17 failed during orbit 142 and 210 respectively. Apart from a reduction in effective area, this led to an increase of the background in the SE recorded by the detectors adjacent to the failed ones.
2.2 Data set
An observation group (hereafter OG) has been built from all pointings
within 20
around Cas A achieved between revolution 7 and revolution 539, plus all pointings outside the galactic plane (
)
during the same period. The former are mainly dedicated observations of the Cas A remnant complemented by some pointings from the periodic Galactic Plane Scan, while the latter are empty fields (at the energy considered) that will help to constrain the background level. All SE and ME2 data with energies between 1130 and 1200 keV were considered (ME2 account for about a quarter of all counts at these energies). This initial selection of pointings has been screened for abnormally high counting rates such as those occurring as the satellite crosses the radiation belts or experiences solar flares. The eventual OG is a three-dimensional space with the following characteristics: 8328 pointings (time axis), 19 single detectors and 42 double detectors (space axis) and 70 one-keV bins ranging from 1130 to 1200 keV (energy axis). This amounts to 6.65 Ms of effective observation time towards Cas A and 12.5 Ms of effective observation time of empty regions, spread over a 4.2 year interval between December 2002 and March 2007.
The corresponding raw spectrum and count rate are presented in Fig. 1. Most of these reflect the instrumental background noise, the origin and handling of which are detailed in the following section.
![]() |
Figure 1: Raw SPI spectrum with arrow indicating the nominal energy of the high-energy 44Ti line ( top). Mean count rate in the 1140-1170 keV band as a function of pointing: black crosses are SPI data and the red line is the sum of background and source models ( middle). Residual count rate as a function of pointing after subtraction of the fitted background and source models to the data, pointings have been rebinned in groups of 20 for clarity ( bottom). |
Open with DEXTER |
2.3 Instrumental background noise
The constant bombardment of the satellite by solar and cosmic-ray particles generates a very strong background noise through the constant activation or excitation of the on-board material. This background noise manifests in the form of fluorescence and nuclear lines on top of a continuum (Jean et al. 2003). It overwhelms the astrophysical signal by a factor of 100 and exhibits time, spectral and even spatial variations within the detector plane. The strong time variability of the background noise is illustrated in the middle plot of Fig. 1, where the increasing trend is due to the evolution of the solar activity that modifies the cosmic-ray flux at the satellite position, while the discontinuities are the result of the selected pointings of the OG being non-contiguous in time. In addition, some of the instrumental background lines correspond exactly to the transitions searched for by the nuclear gamma-ray astronomers. In particular, a background line exists at 1157.0 keV (see upper plot of Fig. 1), and corresponds to the production and subsequent nuclear de-excitation of 44Ca in the instrument. Care should therefore be taken when analysing the data to ensure complete background line subtraction.
Standard methods for analysing the data of a coded-mask instrument make use of detector patterns, which give the relative background count rates in each detector at a given energy and are calibrated from empty field observations. For a given observation, the detector pattern is scaled to the data and the excess counts are then accounted for by fitting a source model or exploited by some image reconstruction algorithm. In the specific case of astrophysical gamma-ray lines, however, the signals are so weak that this approach becomes prohibitive. The variations of the background noise level from pointing to pointing are of the order of the searched signals and the relative count rates in the detectors are subject to some long-term evolution. The search for a weak source using a detector pattern would therefore require a frequent scaling of the pattern (for each pointing in fact) and a few readjustments of it to cope with long-term variations. With such a large number of degrees of freedom, the resulting sensitivity is far too low for most astrophysical gamma-ray lines. Bringing the signals to the fore thus necessitates alternative methods involving fewer parameters so as to improve the sensitivity of the instrument. The true challenge here is the reproduction of the time evolution of the background noise, at each energy and for each detector.
In the approach we have used here, the background is modelled as a linear combination of various templates, which are time-series intended to account for the different sources and temporal trends of the background noise. It should be emphasized, however, that despite all efforts (see for example Weidenspointner et al. 2003), the true nature of the background remains mostly unknown, and as a consequence our approach is empirical. The central element of this stage is the so-called GEDSAT activity tracer, which represents the count rates of all events that saturate the detectors, that is all events with incoming energies above 8 MeV. It gives an estimate of the flux of high-energy cosmic-ray (and solar) particles, and therefore traces the instant activation or excitation rate of the on-board material. The use of the GEDSAT in our modelling is an hypothesis on the time variability of the real background and, as such, it could introduce some bias in the results. Yet, many analyses of SPI gamma-ray line observations have been performed under this prescription (see Knödlseder et al. 2005; Diehl et al. 2006; Wang et al. 2007, for the annihilation, 60Fe and 26Al lines respectively) and produced successful results, so that adopting a similar strategy for the study of the 44Ti lines seems justified; in addition we will ensure that the results are not affected by systematic effects through a thorough analysis of the residuals.
For the low-energy lines at 67.9 and 78.4 keV, the background noise turned out to be so high and so variable that it was impossible to build a satisfactory model. For the analysis of the 1157.0 keV line, on the contrary, the background proved to be quite well represented by a simple three-component background model of the following form:
![]() |
= | ![]() |
|
![]() |
|||
![]() |
(1) |
In the above formula, BGM is the resulting background model and the indices p, d and e indicate the pointing, the detector and the energy bin or band considered. A, B and C are the scaling coefficients of each component, and the functions j, k and l define the parameter set for each component. Tp and T0 are respectively the mean date of pointing p and the date of the beginning of the mission, and GEDSATp,d and LIVEp,d are the rate of saturating events and the livetime of detector d during pointing p. The first component (controlled by factors Aj) is a constant component, supposedly accounting for the part of the background that originates from the long-lived natural radioactivity of the elements making up the satellite. The second component (controlled by factors Bk) is a linear growth from the day the satellite was put into space, and accounts for any build-up of a specific medium-lived radioisotope due to a continuous activation by cosmic-ray bombardment. The third component (controlled by factors Cl) models the prompt component of the cosmic-ray induced background, that is, all background radiation arising from short-lived nuclei or excitation states (short compared to the typical duration of a pointing). The Aj, Bk and Cl factors cannot be estimated a priori with sufficient accuracy and should therefore be determined through fitting to the data.
2.4 Model fitting
In this work, we want to obtain spectral information on a fiducial source of 44Ti emission whose position is well established. The 44Ti signal is therefore extracted through model fitting (as opposed to image reconstruction). Cas A is considered as a point source at (l, b) = (111.73,
-2.13
), since the size of the remnant (roughly 5' on the Chandra images from Hwang et al. 2004) is far below the SPI angular resolution. This source model is mapped into the data-space through a convolution with the instrument response function and then forms the sky model. All background and sky components are then fitted simultaneously to the data using a Maximum Likelihood criterion for Poisson statistics. The source model is fitted independently for each energy bin to extract the source spectra, while each component of the background model is fitted independently for each detector, energy bin and for 3 time periods delimited by the failures of detectors 2 and 17 (because the failures have changed the background noise in the detectors surrounding the dead ones qualitatively and quantitatively).
2.5 Assessment of systematic effects
Before proceeding with the interpretation of the result, it is important to check that our analysis of the SPI data is satisfactory. We want to be sure that our modelling of the background has not introduced systematic errors that might have prevented the detection of a signal.
The residuals left from the subtraction of all fitted models (background and sky) to the data are projected along the time and energy dimensions of the data-space. It is checked using chi-square tests that the residuals are statistically distributed and do not suspiciously group over certain periods or energy ranges. The temporal residuals for the present analysis are shown in the lower plot of
Fig. 1, where they have been rebinned for clarity. The
reduced
when residuals are determined on a pointing basis is 1.0. Yet, when computing the residuals over increasing time-scales, the reduced
increases up to values of about 2.0, indicating that long-term systematic effects are present. Some long-term trends may not be properly accounted for by our background model, which can be explained by the fact that our main ingredient, the GEDSAT rate, traces instant activity only and not build-up effects. However, the propagation of these systematics from the data-space to the final results is not straightforward. Firstly, the instrument response is quite complex and image reconstruction is based on specific correlations between detectors. Secondly, when working with large OG including several thousand pointings and spanning a long time range, it may well be that systematic errors average out.
In order to assess the impact of the systematic errors on the result, we have studied the distribution of sky residuals. Sky residuals correspond to a back-projection of the residuals from the detector plane to the sky, through the coded mask. It is the sky space, rather than the data space, that we are interested in, so that this analysis tells us more than the temporal residuals discussed previously. Formally, sky residuals result from the following calculation:
where
is the image of sky residuals, expressed in number
of
.
and
are respectively the data counts and model (background and sky) counts back-projected on the plane of the sky from the np,d and
,
which are respectively the data and model counts recorded by detector d during pointing p, for a given energy range.
is the instrument response function (IRF) that gives for each direction in the sky the expected distribution of events over the detector plane (derived from Sturner et al. 2003). A normalisation is necessary to ensure that the number of back-projected counts is identical to the number of counts in the data space.
To allow qualitative and quantitative assessment of our sky residuals, simulations were performed for the same OG, using as the true background function our fitted background model and as the true source function a point-source at the position of Cas A with a flux of
ph cm-2 s-1. The sky residuals obtained in the 1150-1165 keV band are presented in Fig. 2, together with a typical simulation. The images are qualitatively similar, suggesting that our residuals are close to what is expected from statistical fluctuations only. This is better illustrated in Fig. 3, where the distribution of our residuals is compared to the purely statistical distribution obtained from 10 simulations. Both distributions match, thereby indicating that the systematic errors revealed by the temporal residuals do not convert into any bias on the sky properties. We can therefore consider that the results of our 1157 keV line analysis are robust.
![]() |
Figure 2: Images of the sky residuals from our analysis ( top) and from a typical simulation ( bottom). |
Open with DEXTER |
![]() |
Figure 3: Distribution of the sky residuals obtained from our analysis (black crosses, only non-zero points are shown) compared to the expected statistical distribution obtained from simulations (red curve). |
Open with DEXTER |
3 Results and discussion
3.1 The spectrum of Cassiopeia A
![]() |
Figure 4: Cassiopeia A spectrum at 1157.0 keV combining SE and ME2; the red solid curve is the fit of a Gaussian shape. |
Open with DEXTER |
The model-fitting procedure described in 2.4 yields the source spectrum presented in Fig. 4, where SE and ME2 have been combined. The spectrum has been fitted by a Gaussian shape with a flat continuum underneath and the best-fit profile corresponds to a line flux of
ph cm-2 s-1 at a position of
keV. This is definitely not significant (although consistent with the expected flux and position within the error bars) and it suggests that the 1157.0 keV line signal from Cas A may be broadened by the Doppler effect, with its flux spread over a too wide energy range and consequently lost in the statistical fluctuations. From this nearly flat spectrum we can derive a statistical lower limit on the line width.
Assuming a total line flux of
ph cm-2 s-1 and no position shift, the 2
lower limit on the measured line width is 3.8 keV FWHM (see the
-curve shown in Fig. 5). The spectral shape eventually appearing in the source spectrum results from the convolution of the intrinsic astrophysical 1157.0 keV line profile with the spectral response of the instrument. Under the assumption that both profiles can be described by Gaussian functions, the measured line width is a quadratic sum of the instrumental and astrophysical line widths. Since the spectral resolution of SPI at 1 MeV is about 2 keV FWHM, our lower limit on the measured line width translates into a lower limit on the astrophysical width of about 3.2 keV FWHM, which in turn implies an
500 km s-1 Doppler velocity at 1157.0 keV. This value is an order of magnitude lower limit on the expansion velocity of the 44Ti ejecta, since our spectrum clearly does not allow us to discriminate between different sophisticated velocity patterns.
3.2 The iron and velocity distribution in Cassiopeia A
From SPI and IBIS we have no detailed positional information about the 44Ti emission but the constraint on the expansion velocity derived above, together with some nucleosynthesis considerations, can be interpreted in the context of all observational data currently available. Since 44Ti is synthesized in explosive Si-burning, a strong spatial correlation between 56Fe and 44Ti is expected. So a good starting point for an interpretation of our 44Ti result is no doubt an inventory of all the iron recorded so far in Cas A.
From XMM-Newton data, Willingale et al. (2003) performed an extensive spectral analysis of Cas A and determined a shocked ejecta mass of 2.2
and a swept-up CSM mass of 7.9
.
The abundances they obtained from spectral modelling indicate a total X-ray emitting iron mass of 0.058
,
which is a sizeable fraction of the canonical 0.1
of iron supposedly ejected in core-collapse events. Accounting for pre-existing iron in the progenitor reduces the value to about 0.04
for the initial solar metallicity, which still constitutes a substantial amount. If we consider that the reverse-shocked ejecta may contain some fraction of cold material (as suggested by Hwang & Laming 2003), presently undetectable in X-rays because of an earlier reverse-shock crossing and subsequent cooling or because of a lower density, the work of Willingale et al. (2003) seems to indicate that a fair amount of the iron produced by Cas A currently lies between the forward and reverse shocks, with rms radial velocities of about 1700-1800 km s-1.
The total iron yield of the Cas A event remains however uncertain. The low brightness of the event recorded (or not, see Stephenson & Green 2005) by Flamsteed in the 1680 s suggests a low 56Ni yield, but the visual extinction towards Cas A might range from 4 to 8 mag. Taking this into account, Young et al. (2006) estimated that the maximum 56Ni mass that could have been produced in the event is 0.2 .
The recent identification of Cas A as the remnant of a SNIIb explosion has allowed us to tighten the constraints on the total iron yield; it very likely lies between 0.07 and 0.15
,
as indicated by the canonical SNIIb SN1993J (see Krause et al. 2008, and references therein). So in the most optimistic scenario, Cas A harbours about 0.09-0.011
in addition to the 0.04-0.06
estimated by Willingale et al. (2003).
From nucleosynthesis considerations, a substantial fraction of these 0.09-0.11
of iron can be expected to lie inside the reverse shock. Indeed, in a spherical supernova explosion 56Fe is thought to form the deepest ejecta. According to the analytic hydrodynamic model used by Laming & Hwang (2003) and Hwang & Laming (2003), 0.2
of material are still freely expanding, so the most of iron could still be there, possibly accompanied by 44Ti. The free expansion velocity at the reverse-shock position was estimated to around 5000 km s-1 by Morse et al. (2004) from nearly undecelerated optical clumps. Therefore, the expansion velocities of the unshocked 56Fe and 44Ti ejecta should range between 0 and 5000 km s-1.
![]() |
Figure 5:
|
Open with DEXTER |
An alternative has been put forward in the particular context of 44Ti production. Maeda & Nomoto (2003) showed that a bipolar supernova explosion caused by a pair of opposite hydrodynamic jets fed by accretion brings the material on the jet axis to higher temperatures than that reached in spherical models, and hence causes a more rapid expansion and hence a stronger alpha-rich freeze-out for that material. This results primarily in an augmented 44Ti yield, but the hydrodynamics associated with that scenario also leads to an inversion of the velocities along the jet where 44Ti and 56Ni have the highest velocities around 20 000 km s-1, followed by 28Si and 16O. The simulation of Maeda & Nomoto (2003) has been invoked to explain the high 44Ti mass
and 44Ti/56Ni ratio inferred for Cas A (Prantzos 2004) and also provides an explanation for the outward increasing iron abundance suggested by Laming et al. (2006) and the outward increasing sulfur-to-oxygen ratio reported by van den Bergh (1971) and Fesen (2001) for the optical FMKs of the jet. Moreover, Laming & Hwang (2003) showed from a comparison of Chandra spectra of individual knots with self-similar hydrodynamical simulations that more energy was injected in the direction of the north-east jet. In the specific case of Cas A, the jet lies within 20
of the plane of the sky (Fesen 2001), so the 20 000 km s-1 of the 44Ti ejecta in this jet-induced nucleosynthesis scenario translates into maximum velocities along the line of sight of about
6400 km s-1.
The 44Ti ejecta can therefore be present in three distinct locations: the freely-expanding ejecta, the reverse-shocked ejecta and the jet. All three sites are associated with expansion velocities that are far higher than our
500 km s-1 lower limit and so the only conclusion we can draw from our INTEGRAL/SPI observations is that the 44Ti does not lie strictly at the mass cut (the asymptotic velocity of which is theoretically zero), contrary to what is obtained in spherically symmetric models.
4 Conclusion
We have searched for the three 44Ti decay lines in Cassiopeia A with the high-resolution spectrometer SPI on board INTEGRAL. Due to the very high instrumental background noise with strong time variability, it was not possible to extract the two lines at 67.9 and 78.4 keV. For the high-energy line at 1157.0 keV, no significant signal is seen in the 1140-1170 keV band, suggesting that the 1157.0 keV line signal from Cassiopeia A is broadened by the Doppler effect. From our spectrum, we derive a 500 km s-1 lower limit at 2
on the expansion velocity of the 44Ti ejecta.
When interpreted in the context of all observational data currently available, especially in optical and X-rays, this result does not allow us to constrain the nucleosynthesis site of 44Ti or the explosion process since the velocities involved throughout the remnant are all far above our lower limit. Further progress on this issue may follow from renewed efforts to understand and model the background noise of the SPI instrument at low energies. In the medium term, interesting results could also arise from the search for the 44Sc K
fluorescence emission that follows the electron capture of 44Ti. Investigations with Chandra as reported by Theiling & Leising (2006) have been fruitless,
but line detection for specific regions may be achieved through extended observations.
Acknowledgements
The SPI project has been completed under the responsibility and leadership of CNES. We are grateful to ASI, CEA, CNES, DLR, ESA, INTA, NASA and OSTC for their support.
References
- Ahmad, I., Greene, J. P., Moore, E. F., et al. 2006, Phys. Rev. C, 74, 065803 [NASA ADS] [CrossRef] (In the text)
- Chevalier, R. A., & Kirshner, R. P. 1979, ApJ, 233, 154 [NASA ADS] [CrossRef] (In the text)
- Diehl, R., Halloin, H., Kretschmer, K., et al. 2006, Nature, 439, 45 [NASA ADS] [CrossRef]
- Dupraz, C., Bloemen, H., Bennett, K., et al. 1997, A&A, 324, 683 [NASA ADS] (In the text)
- Fenimore, E. E., & Cannon, T. M. 1978, Appl. Opt., 17, 337 [NASA ADS] [CrossRef] (In the text)
- Fesen, R. A. 2001, ApJS, 133, 161 [NASA ADS] [CrossRef] (In the text)
- Hughes, J. P., Rakowski, C. E., Burrows, D. N., & Slane, P. O. 2000, ApJ, 528, L109 [NASA ADS] [CrossRef]
- Hwang, U., Holt, S. S., & Petre, R. 2000, ApJ, 537, L119 [NASA ADS] [CrossRef]
- Hwang, U., & Laming, J. M. 2003, ApJ, 597, 362 [NASA ADS] [CrossRef] (In the text)
- Hwang, U., Laming, J. M., Badenes, C., et al. 2004, ApJ, 615, L117 [NASA ADS] [CrossRef]
- Iyudin, A. F., Diehl, R., Bloemen, H., et al. 1994, A&A, 284, L1 [NASA ADS] (In the text)
- Jean, P., Vedrenne, G., Roques, J. P., et al. 2003, A&A, 411, L107 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Knödlseder, J., Jean, P., Lonjou, V., et al. 2005, A&A, 441, 513 [NASA ADS] [CrossRef] [EDP Sciences]
- Krause, O., Birkmann, S. M., Usuda, T., et al. 2008, Science, 320, 1195 [NASA ADS] [CrossRef] (In the text)
- Laming, J. M., & Hwang, U. 2003, ApJ, 597, 347 [NASA ADS] [CrossRef] (In the text)
- Laming, J. M., Hwang, U., Radics, B., Lekli, G., & Takács, E. 2006, ApJ, 644, 260 [NASA ADS] [CrossRef] (In the text)
- Maeda, K., & Nomoto, K. 2003, ApJ, 598, 1163 [NASA ADS] [CrossRef] (In the text)
- Mochizuki, Y., Takahashi, K., Janka, H.-T., Hillebrandt, W., & Diehl, R. 1999, A&A, 346, 831 [NASA ADS] (In the text)
- Morse, J. A., Fesen, R. A., Chevalier, R. A., et al. 2004, ApJ, 614, 727 [NASA ADS] [CrossRef] (In the text)
- Prantzos, N. 2004, in 5th INTEGRAL Workshop on the INTEGRAL Universe, ed. V. Schoenfelder, G. Lichti, & C. Winkler, ESA SP, 552, 15 (In the text)
- Reed, J. E., Hester, J. J., Fabian, A. C., & Winkler, P. F. 1995, ApJ, 440, 706 [NASA ADS] [CrossRef] (In the text)
- Renaud, M., Vink, J., Decourchelle, A., et al. 2006, ApJ, 647, L41 [NASA ADS] [CrossRef] (In the text)
- Reynolds, S. P., Borkowski, K. J., Green, D. A., et al. 2008, ArXiv e-prints, 803 (In the text)
- Roques, J. P., Schanne, S., von Kienlin, A., et al. 2003, A&A, 411, L91 [NASA ADS] [CrossRef] [EDP Sciences]
- Stephenson, F. R., & Green, D. A. 2005, Journal for the History of Astronomy, 36, 217 [NASA ADS] (In the text)
- Sturner, S. J., Shrader, C. R., Weidenspointner, G., et al. 2003, A&A, 411, L81 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- The, L.-S., Clayton, D. D., Jin, L., & Meyer, B. S. 1998, ApJ, 504, 500 [NASA ADS] [CrossRef]
- The, L.-S., Clayton, D. D., Diehl, R., et al. 2006, A&A, 450, 1037 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Theiling, M. F., & Leising, M. D. 2006, New Astron. Rev., 50, 544 [NASA ADS] [CrossRef] (In the text)
- Thorstensen, J. R., Fesen, R. A., & van den Bergh, S. 2001, AJ, 122, 297 [NASA ADS] [CrossRef] (In the text)
- van den Bergh, S. 1971, ApJ, 165, 457 [NASA ADS] [CrossRef] (In the text)
- Vedrenne, G., Roques, J.-P., Schönfelder, V., et al. 2003, A&A, 411, L63 [NASA ADS] [CrossRef] [EDP Sciences]
- Vink, J. 2004, New Astron. Rev., 48, 61 [NASA ADS] [CrossRef]
- Vink, J., Laming, J. M., Kaastra, J. S., et al. 2001, ApJ, 560, L79 [NASA ADS] [CrossRef] (In the text)
- Wang, W., Harris, M. J., Diehl, R., et al. 2007, A&A, 469, 1005 [NASA ADS] [CrossRef] [EDP Sciences]
- Weidenspointner, G., Kiener, J., Gros, M., et al. 2003, A&A, 411, L113 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Willingale, R., Bleeker, J. A. M., van der Heyden, K. J., & Kaastra, J. S. 2003, A&A, 398, 1021 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Woosley, S. E., Arnett, W. D., & Clayton, D. D. 1973, ApJS, 26, 231 [NASA ADS] [CrossRef]
- Young, P. A., Fryer, C. L., Hungerford, A., et al. 2006, ApJ, 640, 891 [NASA ADS] [CrossRef] (In the text)
All Figures
![]() |
Figure 1: Raw SPI spectrum with arrow indicating the nominal energy of the high-energy 44Ti line ( top). Mean count rate in the 1140-1170 keV band as a function of pointing: black crosses are SPI data and the red line is the sum of background and source models ( middle). Residual count rate as a function of pointing after subtraction of the fitted background and source models to the data, pointings have been rebinned in groups of 20 for clarity ( bottom). |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Images of the sky residuals from our analysis ( top) and from a typical simulation ( bottom). |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Distribution of the sky residuals obtained from our analysis (black crosses, only non-zero points are shown) compared to the expected statistical distribution obtained from simulations (red curve). |
Open with DEXTER | |
In the text |
![]() |
Figure 4: Cassiopeia A spectrum at 1157.0 keV combining SE and ME2; the red solid curve is the fit of a Gaussian shape. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
|
Open with DEXTER | |
In the text |
Copyright ESO 2009
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.