arxiv: v1 [astro-ph.he] 29 Aug 2012

Size: px

Start display at page:

Download "arxiv: v1 [astro-ph.he] 29 Aug 2012"

Chloe Collins
5 years ago
Views:

1 Astronomy & Astrophysics manuscript no. G draft c ESO 01 August 30, 01 The Supernova Remnant G in X-rays (Research Note) T. Prinz 1 and W. Becker 1, 1 Max Planck Institute for extraterrestrial Physics, Giessenbachstrasse, 8741 Garching, Germany tprinz@mpe.mpg.de Max-Planck Institut für Radioastronomie, Auf dem Hügel 69, 311 Bonn, Germany arxiv:108.98v1 [astro-ph.he] 9 Aug 01 Received date / Accepted date ABSTRACT Aims. We present a detailed study of the supernova remnant (SNR) G in the 0.-1 kev X-ray band. Methods. Using data from XMM-Newton we performed a spectro-imaging analysis of G in order to deduce the basic parameters of the remnant and to search for evidence of a young neutron star associated with it. Results. In X-rays the remnant is characterized by a bright arc located in the south-west direction. Its X-ray spectrum can best be described by an absorbed non-equilibrium collisional plasma model with a hydrogen density of n H = cm and a plasma temperature of million Kelvin. The analysis revealed a remnant age of 800 to 7600 years and a distance of kpc. The latter suggests a spatial connection with a close-by HII region. We did not find evidence for a young neutron star associated with the remnant. Key words. ISM: supernova remnants - ISM: individual objects: G Stars: neutron 1. Introduction Every year, several hundred supernova (SN) events are observed, representing the end state of stellar evolution in which a corecollapse of a massive star or the thermonuclear disruption of a white dwarf takes place. Most of these events are observed in the optical band by ground based observatories, some in comparing new and old photographic plates. On the other hand, SNe are the most energetic events which we can observe in the universe, making them often as bright as a whole galaxy. Especially the latter property supports observations of SN events up to distances of Gpc (Rodney et al. 011). At this distance, though, the only information obtainable from them is the characteristic raising and fading of their light as a function of time, i.e. their photometric light-curves, and its spectral evolution, which both allow their classification (Sako et al. 008). In our own Galaxy, the rate of observed SNe is small and on average only two per century (Keane & Kramer 008). Many SNe, though, remain unobserved due to optical light extinction. The last recorded SN in our own Galaxy, the Kepler SN, was observed in AD Not more than six other SN have been detected in the Galaxy in the past two thousand years (Green & Stephenson 003). The optical light from the two youngest supernovae known in our Galaxy, G and Cas A, were not observed about 100 and 30 years ago (Reynolds et al. 008; Green & Stephenson 003). However, even if the direct light form a SN event was missed, for several nearby ones we have the chance to study at least the light of their remnant, which remains visible in various spectral bands for about 10 years. Although the light echo from the SN CAS A has been found recently (Rest et al. 008), it is one of the few cases so far where the SN light could be studied a few hundred years after the SN event. The most promising way to learn about the evolution of a SN shock front and the feedback on the evolution of their host galaxy is therefore to study the diffuse emission of supernova remnants. In the last years several new supernova remnants were detected thanks to the increasing sensitivity of modern X-ray observatories. One of those remnants is G , which was first detected in X-rays by Schaudel et al. (00) and later identified to be a SNR by Schaudel (003) and Robbins et al. (01). It is a shell-like SNR from which X-ray and radio emission has been detected. In addition, filaments in the infrared and H α band were detected in the near proximity of the source. A first detailed study of the SNR G was presented by Robbins et al. (01) using ROSAT PSPC data. As this detector had roughly five independent energy channels in the range 0.1- kev, the spectral results deduced in their analysis were very limited and suggested that the X-rays were emitted from a thermal plasma. Robbins et al. (01) therefore were not able to put any constraints with high confidence on the derived spectral parameters. The ROSAT data did not allow to deduce parameters like the age, the expansion velocity of the remnant or its distance. In this letter we report on an XMM-Newton observation which was targeted to SNR G The results of the spatial and spectral analysis of the X-ray data is presented in Section. A discussion is given in Section 3 in which we use the inferred spectral parameters of G to derive an estimate for the age, the radius, the expansion velocity and the distance. Section 4 provides the concluding summary.. X-ray observation and data reduction G was observed for 13.6 ksec on 8 June 011 with the EPIC cameras (Strüder et al. 001) on board the X-ray observatory XMM-Newton (ObsID ). The two MOS and the PN cameras were operated in Full-Frame mode using the medium filter. Part of the observation was taken with the filter 1

wheel in closed position because the observation was strongly affected by particle background radiation. Therefore, the performed duration was only 4.4 ksec and.

2 wheel in closed position because the observation was strongly affected by particle background radiation. Therefore, the performed duration was only 4.4 ksec and. ksec for the PN and MOS1/ detectors, respectively. We used the XMM SAS version to reprocess and reduce the data. Times of high background activity were identified by inspecting the light curves of the MOS1/ and PN data at energies above 10 kev. After rejecting these times the effective exposures of the PN, MOS1 and MOS cameras were 3.7 ksec,.1 ksec, and.3 ksec, respectively. We extracted images and exposure maps in the five standard bands of XMM-Newton using all three EPIC instruments. Single and double events were selected from the PN data and single to quadruple events from the MOS data sets. For the spectral fitting we used the X-Ray Spectral Fitting Package (XSPEC) version 1.7.0u. For the spectral analysis of the extended source emission we restricted the energy range to kev because the count rate detected at higher energies was too sparse for a meaningful spectral analysis. Below 0.4 kev the detector and telescope response is not well established. All given uncertainties in this letter represent the 1σ confidence range for one parameter of interest, unless stated otherwise..1. Spatial analysis Figure 1 shows the X-ray color image of SNR G A bright arc-like structure is clearly detectable in the south-east which is elongated to a elliptically shaped sphere. The center of this ellipse is at RA c = 11 h m ṣ 3±0 ṣ 7, DEC c = ± 3 with a semi major axis of and a semi minor axis of Declination (J000) -6:: :00:00.0 0: :00.0 1:00.0 0: :00.0 7: :6:00.0 :00.0 4:00.0 Right Ascension (J000) Fig XMM-Newton MOS1/ color image of G (red kev, green kev and blue 1.3- kev). The superimposed images are binned with 6 per pixel and smoothed by a Gaussian kernel ofσ=30. The photons detected in the annular sector shown in white were used for spectral analysis of the supernova remnant. We searched for point sources with a signal-to-noise ratio S/N of at least 3σ in the five standard bands using a sliding box source detection algorithm (SAS tool edetect chain). Five point 3 4 sources were detected (see Fig. 1), which all have a S/N 7σ G (σ G is defined asσ G = 1+ c bg + 0.7, where c bg are the background counts). Their position, positional error and signal-tonoise ratio are listed in Table 1. We cross-correlated the source position with the help of the VizieR online tool 1 by searching for optical counterparts within a square of side length 6 in all available catalogs. This is equal to 3σ, wereσis the error in the source position including the XMM-Newton absolute astrometric accuracy of (r.m.s.). We found a match for every detected source. These counterparts with their positional displacement from the X-ray source are listed in Table 1. No source has been detected close to the geometrical center of the supernova remnant. The 3σ upper limit on the count rate of a central compact object is cts/s in the energy range 0. to 1 kev using the merged MOS1/ data... Spectral analysis The energy spectrum of G was extracted from an elliptical annular sector with the center at RA c and DEC c, semi major axes of and. 9, semi minor axes of 3. 7 and. and an opening angle of 140. The background spectrum was derived from a nearby region placed on the same MOS1/ chip with the same size as the source region. The background contribution was found to be 0% in the two MOS cameras and 71% in the PN data. After subtracting the background 13, 136, and 3077 source counts remained in the MOS1, MOS, and PN data. For the spectral analysis of the remnant these counts were binned to have at least 7 counts per bin in the case of the MOS1/ observations and 10 counts per bin for the PN data. We checked whether the spectral fitting results will change, if we model the instrumental background according to the suggestions made by Kuntz & Snowden (008). They proposed to add a Gaussian at the fluorescent line of Al Kα at E = 1.49 kev with zero width and a power law model for modeling the remaining soft proton contamination. This power law model is not convolved with the instrumental response. As it turns out from this analysis the model fits were not significantly better than without adding these components using the standard F-test. Thus, the following results are without modeling the instrumental background separately. We fitted the X-ray spectrum of G with various models: A hot diffuse gas model (MEKAL,χ = for 9 d.o.f.), a model for a collisionally ionized diffuse gas (APEC, χ = for 9 d.o.f.), a power law model (χ = 3.3/9 d.o.f.), a thermal bremsstrahlung model (χ = 43.3 for 9 d.o.f.), a Raymond-Smith diffuse gas model (χ = 11.7/9 d.o.f.), a non-equilibrium collisional plasma that allows a temperature evolution (GNEI,χ = 88.6 for 93 d.o.f.) or has a constant temperature (NEI,χ = 89.3 for 94 d.o.f.), a plane-parallel shock plasma model (PSHOCK,χ = 87. for 94 d.o.f.), an ionization equilibrium collisional plasma model (EQUIL,χ = 96.8 for 9 d.o.f.), and a Sedov model with separate ion and electron temperature (χ = 90.0/93 d.o.f.). For spectral models with χ red < 1.4, all fit parameters are listed in Table. The NEI model fit is shown in Figure. In the following, we will use the best-fit results of the NEI model to derive related SNR parameters. For all models with a non-equilibrium equation of state we measure an ionization timescale valueτ=t 0 N e < 10 1 s cm 3 (t 0 is the age of the remnant and N e the post-shock electron number density), which is 1 xmm.vilspa.esa.es/docs/documents/cal-tn-0018.pdf.

3 Table 1. Detected sources in XMM-Newton observation The sources are denoted as in Fig. 1. Source RA (J000) DEC (J000) δ RA/DEC S/N optical counterpart (positional discrepancy) h:m:s d:m:s arcsec σ G 1 11:7: :04: MASS J (. ) 11:6: :0: USNO-B (. 3) 3 11:: :03: MASS J (4. 9) 4 11:4: :13: MASS J (. 7) 11:4: :01: HD (4. 0) Table. Spectral parameters of the best-fit models for SNR G model n H [10 cm ] k B T [kev] τ α [10 11 s/cm 3 ] k B T e β / kt γ [kev] Norm δ [10 cm ] χ /d.o.f. f X ǫ [10 11 erg cm s 1 ] NEI / / PSHOCK / / /93 8±3 SEDOV < 0.4/ GNEI 1.7± / / EQUIL 1.10± ± / /9. ± APEC 1.3± / / Notes. (α) ionization timescale (β) electron temperature immediately behind the shock front (γ) ionization timescale-averaged temperature kt = t0 T (t)n t e (t)dt/τ (δ) Norm= s 4π[D A (1+z)] Ne N H dv, where D A is the angular diameter distance to the source in cm, N e and N H are the post-shock electron and hydrogen densities in cm 3, respectively. (ǫ) X-ray flux in the energy range 0. to 4.0 kev. smaller than the expected timescale at which collisional ionization equilibrium is reached (Borkowski et al. 001). Therefore, the APEC and EQUIL model can be neglected. For the GNEI model the temperature and the ionization timescale-averaged temperature are almost the same, which is no improvement over the NEI model. Additionally, we will not use the results of the SEDOV model any further, because the determination of the postshock temperature T e in that model is only possible at energies above 3 to 4 kev (Borkowski et al. 001) where the spectrum of G is not well constrained. The value for n H in the NEI model fit is lower than the average integrated hydrogen column density toward SNR G , which is n LAB H = cm (Kalberla et al. 00). This value is based on HI emission line measurements at a radio frequency of 1 cm and refers to the entire hydrogen column density in the line of sight. The temperature of the plasma is Kelvin and no significant deviation of the best-fit statistic was obtained when allowing the abundances to differ from the solar values. Using the derived parameters of the NEI model fit the flux in the 0. to 4 kev band is erg/cm /s. In addition we investigated the spatial variation of the spectral parameters by extracting all photons in the northern and southern part of the remnant. However, no difference can be assessed within the derived errors. 3. Discussion 3.1. Comparison with the ROSAT results Fig.. Spectrum and fitted model (NEI) of the X-ray emission of SNR G In black the MOS1, red the MOS and green the PN data. Schaudel (003) and Robbins et al. (01) analyzed two pointed ROSAT PSPCB observations independently of each other. These data were taken between 1st and 8th of February 1993 and 19th and 1st of February Both authors found a temperature of the X-ray emitting gas and a hydrogen column density n H toward the source that is lower by at least a factor of two when compared with the results deduced in our work. Only Schaudel (003) found an n H that is comparable with our value listed in Table by fitting the Raymond-Smith model, though his value had a much higher uncertainty. As already mentioned in the introduction, the PSPCB had about five independent energy channels, which limits the conclusions drawn from their spectral fittings, especially when 70 spectral bins were used as in the work of Robbins et al. (01). This corresponds to an oversampling of about 1 times the spectral resolution of the detector! 3.. Supernova remnant Using the deduced spectral parameters from the NEI fit we can derive basic properties of the remnant, such as the distance d, 3

4 post-shock hydrogen density N H, sweep-up mass M, the age of the remnant t, the radius in pc R s, and the shock velocity v s. We used the following set of equations, which is described in detail in Prinz & Becker (01, and references therein). d Sedov d Reddening ( ) / ( ) 1/ = 740 θ 3/ E1 f [kpc] T s Norm = n H / ( n H /A V A V /E B V E B V /kpc ) [kpc] V emit = f θ 3 (d Sedov [kpc]) 3 [cm 3 ] Norm N H = 70.7 d Sedov f θ 3 [cm 3 ] M = 1.4 N H m H V [kg] ( ) 1 t = E1 3 T 6 s [yr] R s = 0.34 ( E1 N H v s = R s [km/s]. t N H ) 1 t [pc] Hereinθ = major minor axis is the reduced radius of the remnant in arcmin, E 1 the explosion energy in units of 10 1 erg, T s the fitted plasma temperature, f the filling factor to correct for the morphology of the SNR, Norm the normalization of the spectral fit, n H the fitted hydrogen column density, V emit is the X-ray emitting volume, and m H the mass of a hydrogen atom. From the distribution of mean color excess E B V per kiloparsec derived by Lucke (1978) we found E B V /kpc=0.±0.10 in the direction of the remnant. In addition, we used the relation between n H and the visual extinction A V of Predehl & Schmitt (199) n H /A V = (1.79±0.03) 10 1 cm and the relation between A V and the color excess A V /E B V = 3.±0. (Zombeck 007). This leads to a distance of d Reddening = 9±4 kpc. The Sedov-analysis based distance is d Sedov = kpc. In the following, we give all important quantities in units of d 9.8 = d/9.8 kpc as d Sedov has a smaller uncertainty than other distance estimates. Furthermore, the mean color excess E B V which we used to calculate d Reddening was derived for a reddening layer up to kpc and thus is just a rough estimate. From Figure 1 we infer that only f= 1% of the remnant is bright enough to be used for spectral analysis. Therefore, the post-shock hydrogen density N H is and the sweep-up mass is M= 9M d / 9.8. Assuming that the explosion energy E is equal to the canonical value of 10 1 erg, the age of the remnant is 800 to 7600 years and the radius is R s = d1/6 9.8 pc. We derive a shock velocity v s of km/s. Using the flux values deduced in Section we compute its luminosity to be L 0. 4 X = d9.8 erg/s Compact central object To obtain a rough estimate of the flux upper limit for a hypothetical compact source in the center of the remnant, we assumed that the source is CCO-like. No compact source nor a radio pulsar has been detected in G so far (Robbins et al. 01). Like other CCOs, e.g. the one in the SNR Puppis A, we assume that the spectrum is dominated by blackbody emission with a temperature of Kelvin and a luminosity in the kev band of at least 10 3 erg/s (see Becker 009, for a review). This corresponds to a flux in this energy range of erg cm s 1. Using the WebPIMMS tool with the fitted n H of the SNR the count rate in the 0.-1 kev range is> cts/s for the merged MOS1 and MOS observations, an order of magnitude lower than the derived 3σ upper limit for a point-like source at the center of the remnant. 4. Conclusion and summary The remnant is characterized by a bright arc in the south-west direction and by diffuse emission with low surface brightness in its western part. We showed that the X-ray emission of G is in agreement with coming from a collisionally heated plasma which has not reached equilibrium yet. The Sedov analysis lead to the conclusion that the SNR is about 6600 years old and is expanding with a velocity on the order of 70 km/s. The close-by HII region G has a velocity of +±1 km s 1 based on measured hydrogen recombination lines (Caswell & Haynes 1987). With the standard IAU parameter for the distance to the center of our Galaxy R 0 = 8. kpc and the solar orbit velocity of V 0 = 0 km/s derived by Kerr & Lynden-Bell (1986) and the Galactic rotation model of Fich et al. (1989) we deduce a distance to the HII region of 9.3±0.6 kpc. For the error estimate we assumed an uncertainty in the velocity-to-distance conversion of 7 km/s (e.g., Clifton et al. 1988). The deduced distance d Sedov is in perfect agreement with the distance of the close-by HII region G This is a strong indicator for a spatial connection between the SNR and the HII region as already suggested by Robbins et al. (01). The used XMM-Newton observation was strongly affected by radiation, which lead to a net observation time that was shorter by a factor of three than the approved exposure time. Therefore, only limited statements can be made about the existence of a compact source located near the center of the supernova remnant. Deeper observations are needed to clarify this question in more detail. Acknowledgements. We acknowledge the use of the XMM-Newton data archive. T.P. acknowledges support from and participation in the International Max-Planck Research School on Astrophysics at the Ludwig-Maximilians University. References Becker, W. 009, in Astrophysics and Space Science Library, Vol. 37, Astrophysics and Space Science Library, ed. W. Becker, 91 Borkowski, K. J., Lyerly, W. J., & Reynolds, S. P. 001, ApJ, 48, 80 Caswell, J. L. & Haynes, R. F. 1987, A&A, 171, 61 Clifton, T. R., Frail, D. A., Kulkarni, S. R., & Weisberg, J. M. 1988, ApJ, 333, 33 Fich, M., Blitz, L., & Stark, A. A. 1989, ApJ, 34, 7 Green, D. A. & Stephenson, F. R. 003, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 98, Supernovae and Gamma-Ray Bursters, ed. K. Weiler, 7 19 Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 00, A&A, 440, 77 Keane, E. F. & Kramer, M. 008, MNRAS, 391, 009 Kerr, F. J. & Lynden-Bell, D. 1986, MNRAS, 1, 103 Kuntz, K. D. & Snowden, S. L. 008, A&A, 478, 7 Lucke, P. B. 1978, A&A, 64, 367 Predehl, P. & Schmitt, J. H. M. M. 199, A&A, 93, 889 Prinz, T. & Becker, W. 01, A&A, 44, A7 Rest, A., Welch, D. L., Suntzeff, N. B., et al. 008, ApJ, 681, L81 Reynolds, S. P., Borkowski, K. J., Green, D. A., et al. 008, ApJ, 680, L41 Robbins, W. J., Gaensler, B. M., Murphy, T., Reeves, S., & Green, A. J. 01, MNRAS, 419, 63 Rodney, S. A., Riess, A. G., Strolger, L., et al. 011, in American Astronomical Society Meeting Abstracts #18, #19.01 Sako, M., Bassett, B., Becker, A., et al. 008, AJ, 13, 348 Schaudel, D. 003, PhD thesis, LMU, Germany 4

5 Schaudel, D., Becker, W., Voges, W., et al. 00, in Astronomical Society of the Pacific Conference Series, Vol. 71, Neutron Stars in Supernova Remnants, ed. P. O. Slane & B. M. Gaensler, 391 Strüder, L., Briel, U., Dennerl, K., et al. 001, A&A, 36, L18 Zombeck, M. 007, Handbook of Space Astronomy and Astrophysics: Third Edition, ed. Zombeck, M. (Cambridge University Press)

Radio Observations of TeV and GeV emitting Supernova Remnants

Radio Observations of TeV and GeV emitting Supernova Remnants Denis Leahy University of Calgary, Calgary, Alberta, Canada (collaborator Wenwu Tian, National Astronomical Observatories of China) outline