بقايا المستعر الأعظم Kes 17: مسرّع فعال للأشعة الكونية داخل سحابة جزيئية

Joseph D. Gelfand NYU Abu Dhabi PO Box 903, New York, NY 10276 jg168@cosmo.nyu.edu Daniel Castro Kavli Institute for Astrophysics & Space Research, Massachusetts Institute of Technology Patrick O. Slane Harvard-Smithsonian Center for Astrophysics 60 Garden Street, Cambridge, MA 02138 Tea Temim Observational Cosmology Lab, Code 665 NASA Goddard Space Flight Center Greenbelt, MD 20771, USA John P. Hughes Department of Physics and Astronomy Rutgers University 136 Frelinghuysen Rd, Piscataway, NJ 08854, USA Cara Rakowski United States Patent and Trademark Office 600 Dulany Street, Alexandria, VA cara.rakowski@gmail.com
الملخص

إن بقايا المستعر الأعظم Kes 17 (SNR G304.6+0.1) هي واحدة من عدد قليل، لكنه آخذ في الازدياد، من البقايا المرصودة عبر الطيف الكهرومغناطيسي. في هذه الورقة، نحلل أرصادا حديثة راديوية وبالأشعة السينية وبأشعة $\gamma$ لهذا الجرم، ونبين أن تسريعا فعالا للأشعة الكونية مطلوب لتفسير طيفه العريض غير الحراري. وتشير هذه الأرصاد أيضا إلى أن Kes 17 تتمدد داخل سحابة جزيئية، مع أن تحديدنا لعمرها يعتمد على ما إذا كان التوصيل الحراري أم تبخر الكتل هو المسؤول الرئيس عن مورفولوجيتها الحرارية المملوءة مركزيا في الأشعة السينية. إن الدليل على تسريع فعال للأشعة الكونية في Kes 17 يدعم عملا نظريا حديثا مفاده أن المجال المغناطيسي القوي والاضطراب والطبيعة المتكتلة للسحب الجزيئية تعزز إنتاج الأشعة الكونية في بقايا المستعرات العظمى. ومع أن ثمة حاجة إلى أرصاد إضافية لتأكيد هذا التفسير، فإن مزيدا من دراسة Kes 17 مهم لفهم كيفية تسريع الأشعة الكونية في بقايا المستعرات العظمى.

Subject headings:
ISM: الأشعة الكونية, ISM: أجرام منفردة: Kes 17, ISM: بقايا المستعرات العظمى, أشعة غاما: ISM, الأشعة السينية: جرم منفرد: Kes 17
11affiliation: Affiliate Member, Center for Cosmology and Particle Physics, New York University22affiliation: Oak Ridge Associated Universities, Oak Ridge, TN 37831, USA

1. المقدمة

يعتقد أن بقايا المستعرات العظمى (SNRs) تؤدي دورا مهما في تكوين الوسط بين النجمي متعدد الأطوار (ISM) الموجود داخل المجرات المكونة للنجوم وفي تنظيمه (مثلا، McKee & Ostriker 1977): فهي توزع المعادن المنتجة في انفجار السلف عبر المجرة المضيفة، وتنتج الغبار (مثلا، Salpeter 1977)، وتسرع الأشعة الكونية حتى طاقات $E\sim10^{15.5}~{\rm eV}$ (مثلا، Arnett & Schramm 1973). غير أن الأدلة الرصدية المباشرة الداعمة لهذا الادعاء الأخير نادرة. ومع أن عدة بقايا مستعرات عظمى في درب التبانة حددت على أنها منتجة للأشعة الكونية، فإن واحدة فقط (بقايا مستعر Tycho الأعظم) تظهر دليلا على تسريع البروتونات حتى “الركبة” في طيف الأشعة الكونية، وهي التي يعتقد أنها تفصل الأشعة الكونية المجرية عن خارج المجرية (Eriksen et al., 2011).

يتطلب تحديد ما إذا كانت بقايا المستعرات العظمى مسؤولة عن جماعة الأشعة الكونية المرصودة دراسة بقايا منفردة لتحديد ما إذا كانت تسرع الأشعة الكونية وكيف تفعل ذلك. فبقايا المستعرات العظمى أجرام شديدة التعقيد تتألف من مادة ساخنة من الوسط بين النجمي صدمتها موجة الانفجار المستعرية المتمددة (“الصدمة الأمامية”)، ومن مقذوفات المستعر الأعظم المسخنة بالموجة الصدمية التي يدفعها الوسط بين النجمي المصدوم إلى داخل البقايا (“الصدمة العكسية”)، ومن مقذوفات باردة غير مصدومة، ومن إلكترونات وأيونات نسبية تسرعت عند الصدمة الأمامية و/أو العكسية. إضافة إلى ذلك، يعتمد التطور الديناميكي لبقايا المستعر الأعظم اعتمادا قويا على محيطها (مثلا، Lozinskaya 1992). ويتطلب تحديد خواص الجسيمات النسبية المسرعة داخل بقايا مستعر أعظم أولا قياس الخواص الفيزيائية لهذه المكونات المختلفة.

ويتطلب ذلك تحليل انبعاث بقايا المستعر الأعظم عبر كامل الطيف الكهرومغناطيسي. يتتبع الانبعاث الراديوي لبقايا المستعر الأعظم إلكترونات بطاقة GeV تسرعت في البقايا. وينتج انبعاثها تحت الأحمر (IR) من الغبار والغاز الذري والجزيئي داخل البقايا وخارجها بعد تسخينهما بالصدمات وبالانبعاثات الأعلى طاقة. ويتتبع الانبعاث الحراري في الأشعة السينية لبقايا المستعر كلا من مادة الوسط بين النجمي التي صدمتها الصدمة الأمامية والمقذوفات التي صدمتها الصدمة العكسية. وأخيرا، يتتبع انبعاث أشعة $\gamma$ الإلكترونات النسبية، وربما الهادرونات (الأشعة الكونية)، المسرعة في بقايا المستعر الأعظم (مثلا، Ackermann et al. 2013). وبفضل قدرات رصدية جديدة عند أطوال موجية تحت حمراء (مثلا، Spitzer وAKARI وHerschel) وعند طاقات أشعة $\gamma$ (مثلا، Fermi وH.E.S.S وV.E.R.I.T.A.S وMAGIC)، فإن عدد بقايا المستعرات العظمى المكتشفة في النطاقات الموجية الأربعة جميعا يزداد بسرعة. ومن هذه البقايا Kes 17 (SNR G304.6+0.1). في هذه الورقة، نحلل أرصادا حديثة راديوية وتحت حمراء وبالأشعة السينية وبأشعة $\gamma$ لهذه البقايا (§2)، ونستخدم هذه النتائج لتحديد الخواص الفيزيائية لهذه البقايا ومحيطها (§3). وأخيرا، نلخص نتائجنا ونناقش تبعاتها على التفاعل بين بقايا المستعرات العظمى وبيئاتها (§4).

2. الأرصاد وتحليل البيانات

في هذا القسم، نحلل أرصادا حديثة راديوية (§2.1)، وتحت حمراء (§2.2)، وبالأشعة السينية (§2.3)، وبأشعة $\gamma$2.4) لهذا المصدر.

Refer to caption
Figure 1.— صورة عند 1.4 GHz لـ Kes 17. أعدت صورة الشدة الكلية من البيانات المدمجة الموصوفة في §2.1 باستخدام ترجيح منتظم وتركيب متعدد الترددات وإزالة التفاف بالإنتروبيا العظمى. لهذه الصورة ضجيج rms مقداره 0.72 mJy beam-1 ودقة مقدارها $23\farcs3\times18\farcs4$ (يمثل القطع الناقص في الزاوية اليسرى السفلى حجم الحزمة واتجاهها.) تبين الكنتورات البيضاء مستويات السطوع السطحي 10، 25، 50، 75، و100 $\sigma$. ويشير الخط المتقطع إلى المنطقة المستخدمة لصنع منحنى السطوع المبين في الشكل 2.

2.1. الراديو

رصدت Kes 17 أول مرة عند 408 MHz و5 GHz بواسطة Shaver & Goss (1970) الذي صنف هذا المصدر، استنادا إلى الطيف غير الحراري الذي توحي به تدفقاته عند هذين الترددين، على أنه بقايا مستعر أعظم. وقد دعم هذا التحديد كشف انبعاث مستقطب عند 5 GHz (Milne & Dickel, 1975)، والمورفولوجيا غير المنتظمة الشبيهة بالقشرة التي كشفها تحليل البيانات المأخوذة أثناء مسح مستوى المجرة بتلسكوب Molonglo Observatory Synthesis Telescope (MOST) (Whiteoak & Green, 1996). يشير تحليل طيف H I باستخدام طريقة امتصاص H I القياسية باتجاه Kes 17 إلى مسافة $d>9.7~{\rm kpc}$ (Caswell et al., 1975). وأخيرا وليس آخرا، يتطلب كشف انبعاث ميزر OH (1720 MHz) (Frail et al., 1996) وجود مادة جزيئية مصدومة (Elitzur 1976، وانظر Wardle & McDonnell 2012 لمراجعة حديثة).

Frequency [GHz] Flux Density [Jy] Reference
0.408 29.8 Shaver & Goss (1970)
0.843 18 Whiteoak & Green (1996)
1.4 10.9±0.14 This work
5.0 6.7 Shaver & Goss (1970)
Table 1كثافات التدفق لـ Kes 17 عند عدة ترددات راديوية.

رصدت مصفوفة Australia Telescope Compact Array، وهي في تهيئة 1.5A، هذه البقايا في 2004 مارس 14 عند كل من 1.4 و 2.4 GHz. استخدم هذا الرصد إعداد المقرن ذي أكبر عرض نطاق متاح (حزم 128 MHz عبر 13 قناة)، مع حزمة متمركزة عند 1384 MHz وأخرى عند 2368 MHz. سجل هذا الرصد أنماط الاستقطاب الخطي الأربعة كلها (XX وYY وXY وYX). استخدمنا حزمة البرمجيات MIRIAD (Sault et al., 1995) لمعايرة كثافة التدفق باستخدام رصد لـ PKS B1934-638، ولمعايرة الطور باستخدام بيانات من أرصاد منتظمة لـ PKS 1329665، ولتصوير بيانات Kes 17. ولتحسين حساسيتنا للانبعاث المنتشر عند 1.4 GHz في الحقل، جمعنا مرئيات 1.4 GHz لـ Kes 17 مع بيانات متصلة من Southern Galactic Plane Survey (McClure-Griffiths et al., 2005).

Refer to caption
Figure 2.— منحنى الشدة عند 1.4 GHz لـ Kes 17 على امتداد خط ذي ميل ثابت $\delta=-62^\circ 41^\prime 37\farcs77$ (J2000). عند هذا الميل، تمتلك Kes 17 أكبر امتداد زاوي.

كما هو مبين في الشكل 1، تمتلك هذه البقايا عند 1.4 GHz مورفولوجيا قشرة جزئية بقطر $\sim7\farcm5$ تهيمن عليها حافتان في المنطقتين الجنوبية والشمالية الغربية، تصل بينهما سمة شبيهة بـ“شق” في الجنوب الغربي. ومع أن المنطقة الشمالية الشرقية لها سطوع سطحي $\sim6\times$ أخفض من الحافتين الجنوبية والشمالية الغربية، فإن هناك هبوطا حادا في كثافة التدفق يحدد حافة هذه البقايا (الشكل 2). كما يكشف انبعاث راديوي منتشر داخل القشرة (الشكلان 1 و2). تبلغ كثافة التدفق الكلية (Stokes I) عند 1.4 GHz لـ Kes 17 هي $10.9\pm0.14$ Jy ولم يكشف أي انبعاث مستقطب. وقد حال نقص البيانات ذات تباعد u–v قصير عند 2.4 GHz دون صنع صورة وقياس كثافة تدفق بجودة مماثلة عند هذا التردد الأعلى.

2.2. تحت الأحمر

تحقق أول كشف تحت أحمر لـ Kes 17 بواسطة Infrared Astronomical Satellite (IRAS)، الذي كشف انبعاثا شبيها بالقشرة من بقايا المستعر هذه (Arendt, 1989). وقد كشفت أرصاد أحدث بكاميرا Infrared Array Camera (IRAC) على متن Spitzer من مسح GLIMPSE انبعاثا شديد السطوع في حزم 3.6-8.0 $\micron$ (Lee, 2005; Reach et al., 2006) من قشرة بقايا مستعر أكثر انتشارا. وتكون البنية الخيطية على امتداد الحافة الشمالية الغربية ساطعة بوجه خاص عند 4.5 $\micron$، مما يوحي بأن الانبعاث صادر من H2 مصدوم. واستنادا إلى الألوان وأوجه التشابه المورفولوجي في صور IRAC، خلص Reach et al. (2006) إلى أن هذا الانبعاث تنتجه صدمات جزيئية.

كشف تحليل أرصاد متابعة طيفية عند 5–38 $\micron$ بواسطة Infrared Spectrograph (IRS) على متن Spitzer عن خطوط دورانية نقية ساطعة لـ H2، ما يرجح وجود تفاعل بين بقايا المستعر ومادة جزيئية كثيفة (Hewitt et al., 2009). وتشير نماذج الصدمات إلى أن إثارة خطوط H2 المرصودة تتطلب مكونين صدميين؛ صدمة C أبطأ مقدارها 10 kms1 عبر كتل أكثف ذات n0=106cm3، وصدمة C أسرع مقدارها 40 kms1 تعبر وسطا أقل كثافة ذا n0=104cm3 (Hewitt et al., 2009). ويكشف تحليل الأطياف أيضا خطوط بنية دقيقة ذرية لـ Fe II وNe II و Ne III وS III وS I وS II، تؤدي تدفقات خطوطها النسبية إلى كثافات في مجال 100 – 1000 ${\rm
  cm}^{-3}$ وسرعات صدمية قدرها 150 – 200 kms1 (Hewitt et al., 2009).

ومؤخرا، كشف الانبعاث العريض النطاق في منتصف وتحت الأحمر البعيد من Kes 17 بواسطة Multiband Imaging Photometer for Spitzer (MIPS) عند 24 $\micron$ (Lee et al., 2011; Pinheiro Gonçalves et al., 2011) وبواسطة القمر AKARI عند 15 و24 و65 و90 و140 و160 $\micron$ (Lee et al., 2011). يتركز الانبعاث عند هذه الأطوال الموجية في القشرتين الغربية والجنوبية، ويتداخل جزئيا مع الحافة الراديوية الغربية. ويناسب توزيع الطاقة الطيفي العريض في تحت الأحمر جيدا جسمان أسودان معدلان مع خليط من تراكيب حبيبات غبار كربونية وسيليكاتية. درجات الحرارة الأفضل ملاءمة هي 79±6K و27±3K مع كتل غبار مقدارها 6.2×104M للمكون الساخن و 6.7M للمكون البارد على مسافة 8 kpc (Lee et al., 2011). ومع أن هذه المسافة لا تتسق إلا بدرجة متواضعة مع تلك المقدرة من طيف امتصاص H I لهذه البقايا (§2.1)، فإن هذا الاختلاف لا يغير هذه الكتل تغيرا مهما.

2.3. الأشعة السينية

كشف انبعاث الأشعة السينية من Kes 17 أول مرة في رصد غير منشور مدته $\sim11$ ks بواسطة ASCA (ObsID 57013000) في 1999 فبراير 12. ثم رصد XMM-Newton هذه البقايا في 2005 أغسطس 12 (ObsID 0303100201) لمدة $\sim20$ ks بعد إزالة توهجات الخلفية (Combi et al., 2010). وقد كشف تحليل هذا الرصد أن انبعاث الأشعة السينية كان منتشرا وممتدا، وأكثر سطوعا داخل القشرة الراديوية، وأشار إلى وجود مكونين من الانبعاث السيني غير الحراري والحراري وإلى تغيرات مكانية في طيفه السيني (Combi et al., 2010).

وفي وقت أحدث، رصدت Kes 17 في 2010 سبتمبر 3 بواسطة مرصد Suzaku لمدة $\sim100$ ks (ObsID 505074010). قمنا أولا بإعادة معالجة بيانات Suzaku باستخدام مهمة aepipeline الجديدة في حزمة البرمجيات ftools v6.1111 1 متاحة على http://heasarc.gsfc.nasa.gov/ftools/ (Blackburn, 1995)، ثم استخدمنا منطقة دائرية نصف قطرها $4\farcm7$ متمركزة على Kes 17 لاستخراج طيفها، وحددنا طيف الخلفية باستخدام بيانات من حلقة بين $4\farcm7$ و$7\farcm4$ في نصف القطر ومتمركزة أيضا على Kes 17. وللاتساق، استخدمنا منطقتي المصدر والخلفية نفسيهما للكواشف الثلاثة كلها. أنشئت الأطياف باستخدام xselect، وكذلك RMF وARF لطيف المصدر. ثم لائمنا الطيف المرصود المطروح منه الخلفية في المجال 0.5–10 keV مع نماذج انبعاث مختلفة باستخدام حزمة البرمجيات Sherpa (Freeman et al., 2001)، وتحققنا من نتائجنا باستخدام XSpec (Arnaud, 1996) لأن هاتين الحزمتين تستخدمان خوارزميات مختلفة لتحديد أفضل ملاءمة وأخطاء معاملات النموذج.

Refer to caption
Figure 3.— طيف الأشعة السينية لـ Kes 17 بواسطة Suzaku، وعليه تنبؤ نموذج phabs $\times$ vnei المعطى في الجدول 2.
Parameters vnei vray vgnei vsedov vpshock vnpshock
NH[1022cm2] 3.790.14+0.07 3.690.13+0.13 3.720.09+0.09 3.970.15+0.73 3.780.05+0.05 4.220.10+0.04
kTe [keV] 0.760.01+0.03 0.760.02+0.02 0.760.03+0.03 0.54+0.03 0.760.01+0.03 1.370.03+0.17
kTshock [keV] 0.58+0.16 0.510.07+1.15
kT [keV] 1.10.12+0.26
[Mg] 1.690.29+0.23 1.710.28+0.31 1 1.7 1.680.21+0.23 1
[S] 0.580.07+0.07 0.620.07+0.07 0.580.07+0.07 0.680.37+0.28 0.560.06+0.08 0.690.08+0.15
τ [cm-3 s] >2×1012 2.70.3+0.4×1011 7.71.8×1011 5×1013 6.30.6+23×1010
K 0.0380.003+0.002 0.0320.002+0.003 0.0360.003+0.003 0.0500.012+0.009 0.0380.004+0.002 0.0710.005+0.001
Foreground Source
kT [keV] 0.650.15+0.13 0.660.15+0.12 0.660.14+0.12 0.670.66+14 0.650.16+0.14 0.830.12
K [×105] 1.110.29+0.41 1.120.29+0.29 1.130.29+0.29 1.130.50+1.43 1.060.22+0.40 1.810.62+0.60
χ2 1290.86 1242.07 1273.04 1813.52 1291.76 1848.89
d.o.f 1628 1629 1628 1628 1628 1628
Table 2نتائج الملاءمة المشتركة لطيف الأشعة السينية لـ Kes 17 المقاس بالكواشف الثلاثة كلها في Suzaku. يبين الصف الأول النموذج المستخدم لملاءمة الانبعاث من بقايا المستعر الأعظم Kes 17؛ وفي كل حالة ضرب في نموذج phabs لاحتساب الامتصاص الكهروضوئي. في جميع الملاءمات نمذج الانبعاث الأمامي المذكور في §2.3 باستخدام بلازما Raymond-Smith بدرجة $N_H\equiv0$ وبوفورات شمسية. تشير $\tau$ إلى مقياس زمن التأين (المعرف في §3.2)، وتشير $K$ إلى التطبيع (المعرف في §3.2؛ Arnaud 1996). في نموذج vgnei، تمثل $\langle kT
  \rangle$ درجة حرارة البلازما المتوسطة على مقياس زمن التأين. في نموذجي vpshock وvnpshock، تمثل $\tau$ أعلى مقياس زمن تأين في البلازما – وثبتت أدنى قيمة عند $\tau_l\equiv 0$. ثبتت وفرة Mg في نموذج vgnei عند القيمة الشمسية لأنها كانت القيمة المفضلة عندما سمح لها بالتغير. في نموذجي vnpshock وvsedov، فضلت الملاءمات وفورات Mg غير فيزيائية، لذلك ثبتت عند القيم المقتبسة. حيثما ترد الأخطاء، فإنها تدل على مجال ثقة 90%، وإلا فإن النموذج لم يستطع تقييد قيمة هذا المعامل. وأخيرا، تعني “d.o.f” درجات الحرية.

بسبب الدقة الزاوية على مستوى الدقيقة القوسية في Suzaku، فإن أطيافنا ملوثة بأجرام غير مرتبطة في منطقة المصدر. يشير تحليلنا لرصد XMM السابق إلى وجود ثلاثة مصادر ساطعة وناعمة للأشعة السينية، على الأرجح نجوم أمامية، ضمن منطقة مصدر Suzaku – أحدها $\sim30\times$ أكثر سطوعا من الاثنين الآخرين مجتمعين. ويجد تحليلنا الطيفي لانبعاثها المجمع أنه يعاد إنتاجه جيدا ببلازما Raymond–Smith ذات $N_H\equiv0$ ووفورات شمسية (Anders & Grevesse, 1989). وبما أن الانبعاث من هذه المصادر موجود في أطياف Suzaku، أدرجنا مثل هذا المكون في جميع الملاءمات الطيفية الموصوفة أدناه، مع جعل درجة الحرارة والتطبيع معاملين حرين. في كل ملاءمة، تتسق خواص هذا المكون مع ما قيس لأكثر النجوم سطوعا في بيانات XMM – مما يؤكد الأصل الأمامي لهذا الانبعاث.

كما هو مبين في الشكل 3، توجد عدة خطوط بارزة في طيف الأشعة السينية لـ Kes 17، أبرزها Mg xi وربما Mg xii عند $\sim1.5$ keV، وSi xiii وSi xiv عند $\sim2$ keV، وS xv عند $\sim2.5$ keV، وهي دالة على انبعاث حراري. في معظم النماذج الحرارية، يؤدي افتراض أن البلازما المشعة ذات وفورات شمسية إلى ملاءمات تقلل تقدير تدفق خطوط Mg هذه وتزيد تقدير تدفق خطوط S هذه. وللتحقق مما إذا كانت هذه الفروق تنجم عن معايرة غير مؤكدة لكاشف xis1 حول خط Si K، أعدنا ملاءمة البيانات مع استبعاد هذه القنوات وبافتراض وفورات شمسية. ومع ذلك، رصدت السمة نفسها في بقايا بيانات xis0 وxis3. لذلك سمحنا لوفرتي Mg وS بالتغير في ملاءماتنا.

كما هو مبين في الجدول 2، فإن نمذجة الانبعاث السيني الحراري بنموذج تأين غير متزن (مثلا، نموذج vnei في XSpec وSherpa) تنتج ملاءمة ممتازة ($\chi^2 \approx 0.8$ مخفض) لدرجة حرارة إلكترونية $kT_{\rm
  e} \sim 0.7$ keV، ووفرة دون شمسية لـ S لكن وفرة فوق شمسية لـ Mg، ومقياس زمن تأين كبير جدا $\tau$ (τ2×1012 cm-3 s مع ثقة 90%)، قريب من شرط التأين المتزن (Smith & Hughes, 2010). وقد أبلغ أيضا عن قيمة كبيرة لـ $\tau$ بواسطة Gok & Sezer (2012) في تحليلهم المستقل لهذا الرصد بواسطة Suzaku. وليس مفاجئا أن نمذجة الطيف المرصود بنموذج Raymond & Smith (1977) لبلازما ساخنة منتشرة، الذي يفترض اتزان التأين، تقدم أيضا ملاءمة جيدة جدا للبيانات لدرجة حرارة إلكترونية ووفورات مماثلة (الجدول 2). غير أن هذه النماذج تفترض أن للبلازما درجة حرارة ثابتة وموحدة وحالة تأين واحدة (أي إن البلازما كلها صدمت في الوقت نفسه إلى درجة الحرارة نفسها)، وهذا أمر غير مرجح في بقايا مستعر أعظم. ونتيجة لذلك، لائمنا أيضا طيف الأشعة السينية المرصود لـ Kes 17 بنماذج أكثر تحفيزا فيزيائيا تسمح بمدى من درجات الحرارة لكن بحالة تأين واحدة (vgnei وvsedov) أو بدرجة حرارة واحدة لكن بمدى من مقاييس زمن التأين (vpshock وvnpshock). ومن بين هذه النماذج الأكثر تحفيزا فيزيائيا، كان vgnei وحده قادرا على إعادة إنتاج الطيف المرصود (الجدول 2). ولا يتطلب هذا النموذج وفرة فوق شمسية لـ Mg، ويفضل مقياس زمن تأين $\tau \sim 3\times10^{11}$ cm-3 s، أخفض بكثير من المطلوب في النماذج الأخرى (الجدول 2).

إن نجاح النماذج الحرارية الصرفة في إعادة إنتاج طيف الأشعة السينية المرصود يتعارض مع التحليلات السابقة لـ Gok & Sezer (2012) و Combi et al. (2010) التي تتطلب انبعاثا غير حراري كبيرا. ففي تحليلهم لبيانات Suzaku نفسها، يعيد Gok & Sezer (2012) إنتاج درجة حرارة الإلكترونات ووفرتي Mg وS الواردة في الجدول 2، لكنهم يتطلبون مكون قانون قوة ذي معامل فوتوني $\Gamma \sim 1.4$. (وهم يستخدمون نموذجا حراريا يفترض اتزان التأين، متسقا مع نموذج Raymond-Smith الموصوف أعلاه.) ويحفز مكون قانون القوة فائض انبعاث دون 1.5 keV ناتج عن ملاءمة الطيف المرصود بنموذج حراري واحد (الشكل 3 في Gok & Sezer 2012). في ملاءماتنا الطيفية، يهيمن على هذا المجال الطاقي المكون الأمامي الموصوف أعلاه. إن التشابه بين خواص مكوننا الأمامي والخواص الطيفية السينية لهذه النجوم المقاسة بتحليلنا لرصد XMM يقيم حجة قوية لوصف حراري صرف لانبعاث الأشعة السينية بواسطة Suzaku من Kes 17.

ولا يفسر ذلك الانبعاث السيني غير الحراري الذي ادعاه Combi et al. (2010) في تحليلهم لبيانات XMM. فقد قسم هؤلاء المؤلفون Kes 17 إلى ثلاث مناطق مكانية (لم تتضمن أي منها النجوم الأمامية المذكورة أعلاه) واحتاجوا إلى مكون قانون قوة مهم ذي $\Gamma \sim 1-3$ لإعادة إنتاج الانبعاث المرصود $>4$ keV (الشكل 2 في Combi et al. 2010) في كل منطقة. ويؤكد تحليلنا لبيانات XMM باتباع إجراءاتهم هذه النتيجة. وحتى إذا لم يقسم المرء انبعاث الأشعة السينية المرصود من Kes 17 إلى ثلاث مناطق مكانية، يبقى مكون غير حراري لازما لتفسير طيف الأشعة السينية المقاس بواسطة XMM – فإن ملاءمة الطيف السيني المركب لـ Kes 17 كما قيس بـ XMM ببلازما Raymond-Smith ممتصة واحدة ذات وفرتي Mg وS غير شمسيتين22 2 أفضل معاملات الملاءمة هي $N_H=3.3^{+0.1}_{-0.1}\times10^{22}$ cm-2 و $kT=0.80_{-0.03}^{+0.03}$ keV (تدل الأخطاء على مجال ثقة 90%)، وكانت لهذه الملاءمة قيمة $\chi^2=1674.91$ في 1607 درجات حرية. تقلل بانتظام من تقدير التدفق عند $>4$ keV. وتحسن إضافة مكون قانون قوة إلى هذا النموذج الملاءمة، مخفضة $\chi^2$ إلى 1665.22 مع 1605 درجة حرية. لهذا المكون ذي قانون القوة تطبيع $K_{\rm PL} = 2.0_{-1.9}^{+6.3} \times
10^{-2}$ photons cm-2 s-1 keV-1 عند 1 keV (تدل الأخطاء على مجال ثقة 90%) ومعامل فوتوني $\Gamma=7.5_{-3.6}^{+2.3}$. وهذه القيم متسقة هامشيا مع تحليل Combi et al. (2010)، الذي يذكر مجموعا قدره $K_{\rm PL} \sim
10^{-3}$ photons cm-2 s-1 keV-1 عند 1 keV ومعامل فوتوني $\Gamma=3.1\pm0.3$ في المنطقة الشمالية التي يزعمون أنها تهيمن على انبعاث الأشعة السينية غير الحراري. وهذا المعامل الفوتوني أنعم بكثير من الانبعاث السيني غير الحراري المكتشف من بقايا مستعرات عظمى أخرى (مثلا، Reynolds 2008). ومع ذلك، وفقا لاختبار f هذا، فإن التحسن في $\chi^2$ بإضافة مكون قانون قوة له احتمال $\sim1\%$ أن يكون ناشئا عن المصادفة، ولذلك فله دلالة $\lesssim 3\sigma$.

لتحديد ما إذا كان انبعاث الأشعة السينية غير الحراري الذي أبلغ عنه Combi et al. (2010) متسقا مع تحليلنا لبيانات Suzaku الأعمق بكثير، استخدمنا XSpec لمحاكاة الطيف المتوقع للمكون الأمامي ومكون Raymond-Smith ممتص زائد قانون قوة ذي معامل فوتوني معطى $\Gamma$ وتطبيع $K_{\rm PL}$، ثم لائمناه باستخدام نموذج Raymond-Smith + المكون الأمامي الموصوف أعلاه. والحد الأعلى الناتج على $K_{\rm PL}$ هو أعلى قيمة لـ $K_{\rm PL}$ استطاع نموذجنا الحراري الصرف عندها ملاءمة الطيف المحاكى بقيمة $\chi^2 < 2$ مخفضة. وبما أن Combi et al. (2010) يزعمون أن المعامل الفوتوني $\Gamma$ لانبعاث قانون القوة يتغير بين $\Gamma \sim 2
- 3$ في مناطق مختلفة، حددنا الحد الأعلى على $K_{\rm
  PL}$ لكل من $\Gamma=2$ و$\Gamma=3$. عند $\Gamma=2$، نطلب $K_{\rm PL}<1.5\times10^{-4}$، أما عند $\Gamma=3$ فـ $K_{\rm
  PL}<5\times10^{-4}$. وكلا الحدين الأعلىين غير متسق مع نتائج Combi et al. (2010)، إذ تطلبت ملاءماتهم للمناطق الجنوبية والوسطى والشمالية من Kes 17 قيمة مجمعة $K_{\rm PL} \sim 10^{-3}$. وبسبب الدلالة الإحصائية المنخفضة نسبيا للمكون غير الحراري في طيف بقايا المستعر المركب المستخرج من بيانات XMM، وتدفقه الأعلى بمقدار $\sim2-5\times$ مما يسمح به الطيف المستخرج من رصد Suzaku الذي كشف $\sim3\times$ فوتونا أكثر من Kes 17 مقارنة بـ XMM، نخلص إلى أنه لا يوجد انبعاث سيني غير حراري مهم مكتشف من Kes 17.

مع أن Suzaku لا يملك الدقة الزاوية اللازمة لكشف تغيرات مكانية مباشرة في انبعاث الأشعة السينية من Kes 17، فإن من الممكن استخدام مجموعة البيانات هذه لاختبار مثل هذه الادعاءات (Combi et al., 2010). فإذا كانت صحيحة، فإن نمذجة طيف Suzaku المرصود بثلاثة نماذج pshock ممتصة + قانون قوة التي استخدمها Combi et al. (2010) بالإضافة إلى المكون الأمامي المذكور أعلاه يجب أن تنتج ملاءمة أفضل من النماذج الحرارية المفردة المستخدمة أعلاه. ولم يكن هذا هو الحال. ومع ذلك، وعلى الرغم من امتلاكها درجات حرية أقل، فقد كانت للملاءمة الناتجة قيمة $\chi^2$ أسوأ من الملاءمات الطيفية الواردة في الجدول 2، حتى عندما ثبتنا قيم $N_H$ و$kT$ والوفرة و$\tau$ و$\Gamma$ لكل مكون pshock + قانون قوة على القيم التي أوردها Combi et al. (2010) (وسمحنا للتطبيعات بالتغير لاحتساب اختلاف أحجام مناطق الاستخراج). لذلك نخلص إلى أن التغيرات المكانية في طيف الأشعة السينية لـ Kes 17 التي أبلغ عنها Combi et al. (2010) غير متسقة مع الطيف الكلي لهذه البقايا، وهي على الأرجح نتيجة مزيج من اختلاف الدقات المكانية والطيفية بين XMM وSuzaku ومن الارتياب المنهجي والإحصائي في خلفية XMM.

Refer to caption
Figure 4.— خريطة TS ملساء من Fermi LAT للأحداث المحولة في القسم front في المجال من 2 إلى 200 GeV لمنطقة $0\fdg6 \times
0\fdg6$، متمركزة على SNR Kes 17. تقسيم البكسلات هو $0\fdg01$، وملست الخرائط بغوسيات عرضها $0\fdg2$. تمثل الكنتورات الخضراء الانبعاث الراديوي (0.843 GHz) من أرصاد MOST. وتعرض إحصاءات الاختبار كنتورات بيضاء (81-100-121).

2.4. أشعة غاما

تكشف Kes 17 أيضا عند طاقات أشعة $\gamma$ في نطاق GeV (Wu et al., 2011). ولتحديد خواصها، حللنا بيانات 39 شهرا (من 2008 أغسطس حتى 2012 فبراير) جمعها Large Area Telescope (Fermi-LAT) على متن Fermi Gamma-ray Space Telescope. لا ندرج في هذا التحليل إلا الأحداث المنتمية إلى صنف Pass 7 V6 Source، الذي يخفض معدل الخلفية المتبقية (Abdo et al. in prep). ونستخدم أيضا دوال استجابة الجهاز المحدثة (“Pass7 version 6”؛ Rando & for the Fermi LAT Collaboration 2009، Abdo et al. in prep) (IRFs)، ونخفض مساهمة أشعة $\gamma$ الأرضية المرتدة بتعيين زاوية سمت عظمى للفوتونات الواردة تساوي 100 (Abdo et al., 2009b). استخدمنا Fermi Science Tools v9r23p133 3 توزع حزمة Science Tools والوثائق المتعلقة بها بواسطة Fermi Science Support Center على http://fermi.gsfc.nasa.gov/ssc، وطبقنا تقنية ملاءمة الاحتمالية العظمى لتحليل الخصائص المورفولوجية والطيفية لمصدر أشعة $\gamma$ (Mattox et al., 1996). وننمذج انبعاث الخلفية المنتشر في gtlike بمكون مجري ينتج من تفاعلات الأشعة الكونية مع كل من ISM والفوتونات، وبمكونات متناحية تمثل الخلفيات خارج المجرية والمتبقية. يستخدم ملف mapcube gal_2yearp7v6_v0.fits لوصف انبعاث أشعة $\gamma$ من درب التبانة، وينمذج المكون المتناحي باستخدام جدول iso_p7v6source.txt.

Refer to caption
Figure 5.— توزيع الطاقة الطيفي لمصدر أشعة $\gamma$ المتوافق موضعيا مع Kes 17. تمثل الأسهم حدود ثقة 95% عند هذه الطاقات. وتمثل أشرطة الخطأ السوداء الارتيابات الإحصائية (تقديرات 1$\sigma$ مبنية على هسيان عكسي عند القيمة المثلى لسطح لوغاريتم الاحتمالية). وتمثل أشرطة الخطأ الحمراء الارتيابات المنهجية، وهي مجموع تربيعي للارتياب المرتبط بدوال استجابة الجهاز (IRFS)، الذي نحصل عليه من فريق الجهاز (كما ذكر في §2.4) والارتياب المرتبط بتغيرات شدة الخلفية المجرية المنتشرة، المستمد بتغيير (وتثبيت) تطبيع مكون الخلفية المجرية في مكتبة المصادر للملاءمة إلى 94% و106% من أفضل قيمة ملاءمة تم الحصول عليها من ملاءمة البيانات في خانة طاقة معينة.

درست الخصائص المكانية لانبعاث أشعة $\gamma$ في حقل Kes 17 باستخدام فوتونات بين 2 و200 GeV حُولت في القسم front من LAT. ولهذه العينة الفرعية من بيانات أشعة $\gamma$، تبلغ زاوية نصف قطر الاحتواء 68% للفوتونات ذات السقوط العمودي $\leq 0\fdg3$. أنشأنا خرائط إحصاء الاختبار44 4 إحصاء الاختبار هو النسبة اللوغاريتمية بين احتمالية وجود مصدر نقطي في موضع معين ضمن شبكة $L_{\rm ps}$، واحتمالية النموذج دون المصدر الإضافي $L_{\rm null}$، $2\rm{log}(L_{\rm{ps}}/L_{\rm{null}})$. (TS) مع احتساب الخلفيتين المجرية والمتناحية باستخدام gttsmap واستخدمنا هذه الخريطة لتحديد الدلالة الإحصائية وموضع المصدر وامتداده المحتمل. وكما هو مبين في الشكل 4، فإن مصدرا بدلالة $\sim 11.1\sigma$ (ذروة TS $\approx124$) غير محلول (نصف قطر ثقة $95\%$ مقداره $=4\farcm2$) لأشعة $\gamma$ وبمركز (α2000,δ2000=13h06m05s,624254′′) يتوافق مع الانبعاث الراديوي. ولا تظهر خريطة TS المتبقية، التي بنيت بنمذجة مصدر نقطي عند مركز الانبعاث الأفضل ملاءمة، أي دليل على أن المصدر ممتد مكانيا، إذ إن قيم TS المتبقية هي $<3\sigma$ ضمن 1 من المركز.

حددنا توزيع الطاقة الطيفي (SED) لمصدر أشعة $\gamma$ المرتبط بـ Kes 17 باستخدام بيانات من فوتونات ذات طاقة بين 0.2 و204.8 GeV حُولت في قسمي front و back. استبعدنا الفوتونات دون 200 MeV لأنه، في هذا المجال الطاقي، تتغير المساحة الفعالة للجهاز بسرعة وتوجد ارتيابات كبيرة مرتبطة بالنموذج المجري المنتشر. استخدمنا gtlike لنمذجة التدفق في كل خانة طاقة، وقدّرنا معاملات أفضل ملاءمة من خلال تقنية الاحتمالية العظمى. ولنمذجة الخلفية في ملاءمات الاحتمالية ندرج مصادر من كتالوج Fermi LAT Second Source Catalog ذي 24 شهر (Nolan et al., 2012)55 5 تتاح بيانات المصادر 1873 في Fermi LAT Second Source Catalog بواسطة Fermi Science Support Center على http://fermi.gsfc.nasa.gov/ssc/data/access/lat/2yr_catalog/. تمتلك دوال IRFs “Pass7 version 6” التي استخدمناها ارتيابات منهجية معتمدة على الطاقة في المساحة الفعالة: 10% عند 100 MeV، تنخفض إلى 5% عند 560 MeV، وترتفع إلى 20% عند 10 GeV (Abdo et al. 2009a، Abdo et al. in prep والمراجع الواردة هناك). كما قربنا أثر عدم اليقين في المستوى المجري المنتشر الأساسي عبر تغيير تطبيع الخلفية المجرية اصطناعيا بمقدار $\pm6\%$ من أفضل قيمة ملاءمة في كل خانة طاقة، على نحو مشابه لتحليل Castro & Slane (2010). وكما هو مبين في الشكل 5، بالنسبة إلى الطاقات $<800$ MeV و$>51.2$ GeV تحدد البيانات حدودا عليا للتدفق فقط. ويصف توزيع الطاقة الطيفي الناتج جيدا قانون قوة بمعامل طيفي $\Gamma = 2.0 \pm 0.3$ وتدفق فوتوني متكامل فوق 100 MeV قدره $F_{>100\rm{ MeV}}\approx
1.6\times10^{-8}$ photons/cm2/s – مماثل لما قاسه Wu et al. (2011).

3. التفسير

كما وصف في §1، يمكن من خلال دراسة الانبعاث عريض النطاق من بقايا مستعر أعظم تحديد خواص كل من المادة داخل البقايا وفي الوسط بين النجمي المحيط. نحلل أولا الانبعاث غير الحراري المرصود من Kes 17 لتحديد الأصل الفيزيائي لانبعاث أشعة $\gamma$ منها (§3.1)، ثم نستخدم تلك النتائج لتقدير عمر هذه البقايا وطبيعة بيئتها (§3.2). وبما أن المورفولوجيا الراديوية الشبيهة بالقشرة لـ Kes 17 توحي بأن هذا الانبعاث ينشأ عند الصدمة الأمامية أو قربها (§2.1، الشكل 1)، نفترض أن لهذه البقايا نصف قطر $R_{\rm snr} \approx 10d_{10}~{\rm pc}$ عند مسافة $d=10d_{10}~{\rm kpc}$.

Refer to caption
Refer to caption
Refer to caption
Figure 6.— ملاءمات عريضة النطاق لأرصاد الراديو (مربعات مفتوحة؛ Shaver & Goss 1970; Whiteoak & Green 1996) وFermi-LAT (دوائر حمراء) لـ Kes 17 باستخدام نماذج A1 (أعلى)، وA2 (وسط)، وA3 (أسفل). وتعرض الأطياف المنمذجة من الانبعاث السنكروتروني (أسود)، وانبعاث كومبتون العكسي (أخضر)، واضمحلال $\pi^0$ (أحمر)، والفرملة غير الحرارية (أزرق). يشير الخط البنفسجي المتقطع إلى الحد الأعلى لانبعاث الأشعة السينية غير الحراري المحدد من أرصاد Suzaku، باستخدام نموذج قانون قوة بمعامل $\Gamma=2$.

3.1. منشأ انبعاث أشعة $\gamma$ من Kes 17

يتطلب تحديد خواص الإلكترونات والبروتونات التي سرّعتها هذه البقايا نمذجة الخصائص الطيفية عريضة النطاق لانبعاثها غير الحراري. نفترض أن الانبعاث الراديوي غير الحراري هو إشعاع سنكروتروني إلكتروني، في حين أن انبعاث أشعة $\gamma$ في نطاق GeV هو مزيج من تشتت كومبتون العكسي (IC) للفوتونات المحيطة بواسطة إلكترونات نشطة، ومن الفرملة غير الحرارية (NT)، ومن اضمحلال $\pi^0$s المنتجة في تصادمات بين هادرونات عالية الطاقة (بصورة رئيسية بروتونات) وبروتونات أدنى طاقة.

تشير ملاءمة كثافات التدفق المرصودة $S_\nu$ عند ترددات راديوية مختلفة (الجدول 1) بقانون قوة ($S_\nu
\propto \nu^\alpha$) إلى معامل طيفي راديوي $\alpha \approx
-0.6$. ونقيد ملاءماتنا أكثر باستخدام الحدود العليا لانبعاث الأشعة السينية غير الحراري المستمدة في §2.3. ننظر في ثلاثة سيناريوهات لمنشأ أشعة $\gamma$ المرصودة في نطاق GeV، لكل منها آلية انبعاث مهيمنة مختلفة: انبعاث IC أو الفرملة NT أو اضمحلال $\pi^0$. نفترض أن التوزيع الطيفي $dN_{p,e}/dE$ للجسيمات المسرعة في Kes 17 هو:

dNp,edE = ap,eEΓp,eexp[EE0p,e], (1)

حيث تمثل $E_{0\,p,e}$ طاقة القطع للبروتونات/الإلكترونات (مثلا، Reynolds 2008; Castro & Slane 2010). وتعطى نسبة الإلكترونات إلى البروتونات عند الطاقات النسبية بمعاملات تطبيع توزيعات هذه الجسيمات، $K_{ep}\equiv a_e/a_p$. وتحصل هذه المعاملات بجعل الطاقة المتكاملة الكلية في الجسيمات المسرعة داخل قشرة بقايا المستعر مساوية لـ $E_{\rm
  cr}=\eta_{\rm cr} E_{\rm sn}$، حيث إن $\eta_{\rm cr}$ هو متوسط كفاءة الصدمة في إيداع الطاقة في بروتونات الأشعة الكونية و$E_{\rm sn}$ هي الطاقة الحركية الابتدائية لمقذوفات المستعر الأعظم. في جميع النماذج نفترض $E_{\rm sn}=10^{51}$ erg، و$\Gamma_ p=
\Gamma_e = 2.0$ (وهو كل من المعامل الذي تتنبأ به آلية فيرمي الأساسية للتسريع والقيمة المستمدة من ملاءمة طيف أشعة $\gamma$ المرصود بنموذج قانون قوة؛ §2.4)، وكثافة عددية للإلكترونات $n_e=1.23 \bar{n}$ (وهي تقابل مادة ذات وفورات شمسية)، حيث $\bar{n}$ هي الكثافة العددية المتوسطة حجميا للبروتونات في ISM المحيط بـ Kes 17.

Kep n¯ B2 ηcr E0e E0p FIC Fbrem Fπ
Model [cm-3] [μG] [TeV] [TeV] [109 ph cm-2 s-1]
IC Dominated (A1) 0.1 1.0 10 0.6 2.2 2.2 6.2 1.5 2.6
NT Bremss. Dominated (A2) 0.5 12 15 0.08 1.5 1.5 2.9 8.9 3.2
π0 Decay Dominated (A3) 0.01 9.0 70 0.4 0.9 30 0.3 0.8 13
Table 3نتائج ملاءمة الطيف عريض النطاق لـ Kes 17 عند افتراض آليات مهيمنة مختلفة لانبعاث أشعة $\gamma$ في نطاق GeV. $K_{ep}$ هي نسبة الإلكترونات إلى البروتونات عند الطاقات النسبية، و$\bar{n}$ هي الكثافة المتوسطة لـ ISM المحيط، و$\eta_{\rm cr}$ هي كفاءة تسريع الأشعة الكونية، و$B_2$ هي شدة المجال المغناطيسي مباشرة خلف الصدمة، و$E_{\rm 0e}$ هي طاقة القطع للإلكترون المسرع، و$F_{\rm IC}$ و$F_{\rm brem}$ و$F_{\pi}$ هي على الترتيب تدفق انبعاث كومبتون العكسي والفرملة غير الحرارية واضمحلال $\pi^0$ عند $>100$ MeV.

ننمذج الانبعاث من اضمحلال $\pi^0$ باستخدام عمل Kamae et al. (2006)، الذي يتضمن عامل تحجيم 1.85 للهيليوم والنوى الأثقل (Mori, 2009) كما وصفه Castro & Slane (2010). وتتبع مكونات الانبعاث السنكروتروني وIC النماذج التي قدمها Baring et al. (1999) والمراجع الواردة هناك، ونمذج انبعاث الفرملة NT باستخدام الوصفة التي قدمها Bykov et al. (2000). ونفترض أن حقل الفوتونات المهيمن لتشتت IC هو الخلفية الكونية الميكروية (CMB؛ $kT_{\rm{CMB}} = 2.725$ K). ونفترض أيضا أن Kes 17 في طور Sedov من تطورها، حيث تنضغط المادة المصدومة بعامل 4. وإذا كانت المادة المجروفة في Kes 17 تبرد إشعاعيا، فإنها ستنضغط بعامل $\gg4$. وسيؤدي ذلك إلى زيادة ملحوظة في كثافة بروتونات الأشعة الكونية المحيطة المجروفة بواسطة بقايا المستعر المتمددة، معززا انبعاث أشعة $\gamma$ منها (Chevalier, 1999)، وكذلك كثافة إلكترونات الأشعة الكونية المحيطة المجروفة، وربما معززا انبعاث أشعة $\gamma$ منها أيضا. ومع ذلك، فإن الأعمار المقدرة في §3.2 تشير إلى أن الأمر ليس كذلك.

بنينا سيناريوهات تهيمن فيها كل آلية ممكنة لانبعاث أشعة $\gamma$ (IC والفرملة NT وانبعاث اضمحلال $\pi^0$) وذلك بتعديل قيم $K_{ep}$ و$\bar{n}$ و$E_{0e}$ وشدة المجال المغناطيسي بعد الصدمة $B_2$، ثم ملاءمة الطيف العريض المرصود. وكما هو مبين في الشكل 6، تستطيع آليات انبعاث أشعة $\gamma$ الثلاث المهيمنة في نطاق GeV جميعا إعادة إنتاج الطيف العريض لهذه البقايا، شريطة أن تكون معاملات النماذج التمثيلية مدرجة في الجدول 3. غير أن نموذجنا يتطلب Kep0.1 لإعادة إنتاج كثافة التدفق الراديوية المرصودة إذا كانت إشعاعات IC (A1) أو الفرملة NT (A2) تهيمن على انبعاث أشعة $\gamma$ – وهذا لا يتسق مع القيمة المقاسة محليا للأشعة الكونية $K_{ep} \sim 0.01$ (Yoshida, 2008). لذلك نخلص إلى أن اضمحلال $\pi^0$ هو المسؤول الرئيس عن انبعاث أشعة $\gamma$ في نطاق GeV المكتشف من Kes 17. لاحظ أننا نفترض $K_{ep}=0.01$ في سيناريو اضمحلال $\pi^0$؛ وتحصل نتائج مماثلة لقيم أخفض من $\eta_{\rm cr}$ و$K_{ep}$ وقيم أعلى من $\bar{n}$ لأن $\eta_{\rm cr} \propto \bar{n}^{-1}$ في نموذج انبعاث اضمحلال $\pi^0$ هذا (Drury et al., 1994). ونتيجة لذلك، تتطلب نمذجتنا أن تكون n¯9 cm-3 وηcr0.4 لكي يهيمن اضمحلال $\pi^0$ على انبعاث أشعة $\gamma$ من Kes 17.

اضمحلال $\pi^0$ هو آلية انبعاث أشعة $\gamma$ المهيمنة حتى لو كانت ادعاءات وجود انبعاث سيني غير حراري من Combi et al. (2010) أو Gok & Sezer (2012) صحيحة. فإعادة نمذجة الأطياف العريضة المرصودة لتدفقات الأشعة السينية غير الحرارية المذكورة في هاتين الورقتين تتطلب مرة أخرى قيما لـ $K_{ep}$ أعلى بكثير مما يقاس محليا إذا هيمن انبعاث IC و/أو الفرملة NT عند طاقات GeV. وتتطلب هذه التدفقات أيضا قيمة أعلى لـ $E_{0e}$، مما يغير لاحقا $\bar{n}$ و$B_2$ و$\eta_{\rm cr}$، لكن بسبب التغايرات بين هذه المعاملات، يكون الأثر الصافي غير مؤكد.

يمكن أن يؤدي إهمال حقول الفوتونات الخلفية غير CMB (كثافة الطاقة $u_{\rm CMB}=0.260$ eV cm-3) إلى مبالغتنا في تقدير $K_{\rm
  ep}$ عند النظر في انبعاث أشعة $\gamma$ المهيمن عليه إشعاع IC. يوحي موقع Kes 17 في مستوى المجرة بوجود حقول فوتونية أكثر طاقة، ومن شأن إدراجها أن يخفض القيمة المطلوبة لـ $K_{\rm ep}$. ومع أن انبعاثا تحت أحمر يكشف من Kes 17 نفسها (§2.2؛ Lee et al. 2011)، فإن كثافة طاقته 66 6 نحسب $u\sim8\times10^{-4}$ eV cm-3 للجسم الأسود المعدل عند $T
  \approx 79$ K و$L\approx 1500~L_\odot$ و$u\sim6.5\times10^{-3}$ eV cm-3 للجسم الأسود المعدل عند $T\approx 27$ K و$L\approx
  12500~L_\odot$. منخفضة جدا بحيث لا تعدل قيمة $K_{\rm ep}$ تعديلا مهما.

حقول الفوتونات الأخرى الممكنة هي ضوء النجوم المحيط ($T\sim 5000$ K) والانبعاث من الغبار الدافئ ($T \sim 25$ K). ولتحديد ما إذا كان بإمكانها السماح لـ IC بالهيمنة على أشعة $\gamma$ المرصودة، نقدر كثافة طاقة الفوتونات اللازمة ليكون تشتت IC للإلكترونات النسبية في البقايا على هذه الفوتونات مسؤولا رئيسيا عن معظم انبعاث أشعة $\gamma$ إذا كان $K_{\rm ep}=0.01$. وبما أن كثافة طاقة الفوتونات المطلوبة تنخفض عند طاقة أكبر في الإلكترونات النسبية $E_{\rm elec}$، فإننا نستخدم أقصى طاقة إلكترونية تسمح بها البيانات. وتقترح كفاءة تسريع أشعة كونية $\eta=0.4$ (كما يوحي نموذج اضمحلال $\pi^0$) أن $E_{\rm
  cr} = 0.4\times10^{51}$ ergs. ولـ $K_{\rm ep}=0.01$ وطاقات قطع الإلكترونات المعطاة في الجدول 3، تكون الطاقة الكلية في الإلكترونات النسبية $E_{\rm elec} =
3.4\times10^{48}$ ergs. وباختيار كفاءة تسريع أشعة كونية عالية إلى هذا الحد، من المرجح أننا نبالغ في تقدير الطاقة الحقيقية في الإلكترونات النسبية، ولذلك فإن تحليلنا يقلل من كثافة الطاقة المطلوبة.

إذا هيمن على الفوتونات الخلفية الانبعاث من الغبار الدافئ، فيجب أن تكون كثافة طاقتها udust40uCMB. وإذا هيمن ضوء النجوم على الفوتونات الخلفية، فيجب أن تكون كثافة طاقتها ustarlight500uCMB. وبما أن مزيجا من الاثنين محتمل، فقد لائمنا أيضا كثافة طاقة ضوء النجوم المطلوبة بافتراض $u_{\rm dust} = 20u_{\rm CMB}$. في هذه الحالة تكون القيمة المطلوبة ustarlight200uCMB. في كل سيناريو، تكون كثافة الطاقة المطلوبة لكل من انبعاث الغبار وضوء النجوم أعلى بكثير من القيم المقدرة من نمذجة انبعاث أشعة $\gamma$ المنتشر المرصود في درب التبانة (Strong et al., 2000). ونتيجة لذلك، فمن غير المرجح أن تمتلك الفوتونات من ضوء النجوم والغبار الدافئ الساقطة على Kes 17 كثافات الطاقة العالية اللازمة لكي يكون انبعاث IC هو آلية انبعاث أشعة $\gamma$ المهيمنة. لذلك نخلص إلى أن اضمحلال $\pi^0$ هو غالبا المسؤول عن معظم انبعاث أشعة $\gamma$ المرصود من Kes 17.

إذا كان هذا صحيحا، فإن Kes 17 هي واحدة من بقايا مستعرات عظمى قليلة (10) ذات دليل رصدي مباشر على تسريع البروتونات إلى طاقات عالية. في الجدول 4، نقارن الخواص الفيزيائية لـ Kes 17 بتلك البقايا الأخرى. في كثير من هذه البقايا، تفترض النمذجة الطيفية عريضة النطاق أيضا أن $K_{ep} \sim 0.01$. وبسبب غياب المعلومات الطيفية عند طاقات TeV، لا نستطيع إلا تقييد $E_{0p} > 500$ GeV لـ Kes 17 (الجدول 3) – وهي أعلى من القيمة المرصودة لـ SNR IC 443، لكنها متسقة مع القيمة المقاسة لبقايا مستعرات عظمى أخرى (الجدول 4). كما أن شدة المجال المغناطيسي بعد الصدمة، $B_{\rm snr}$، تقع ضمن المجال الذي تمتده بقايا المستعرات العظمى المنتجة للأشعة الكونية – وإن كانت أقرب إلى القيمة المستنتجة للبقايا الأقدم (tage4000 سنة) (مثل IC 443 وW51C وW28) منها إلى البقايا الأصغر عمرا ($t_{\rm age} < 1000$ سنة) (مثل Cas A وبقايا مستعر Tycho الأعظم) ضمن هذه المجموعة. أما الكثافة المتوسطة المسموح بها $\bar{n}$ فتقع ضمن المجال الذي تمتده هذه المجموعة، لكن متوسط كفاءة تسريع الأشعة الكونية $\eta_{\rm cr}$ المطلوب لأدنى قيمة مسموحة من $\bar{n}$ مرتفع جدا: 5× أعلى من كفاءة البقايا الأصغر عمرا و$\sim2\times$ أعلى من أي من بقايا المستعرات الأقدم المنتجة للأشعة الكونية. ولا يمكننا تحديد ما إذا كانت هذه البقايا منتجا فعالا بوجه خاص للأشعة الكونية إلا بتحديد الخواص الفيزيائية للبيئة المحيطة بـ Kes 17، كما نفعل في §3.2.

3.2. بيئة Kes 17

كما نوقش في §3.1، فإن منشأ هادروني في المقام الأول لأشعة $\gamma$ في نطاق GeV المكتشفة من Kes 17 يتطلب أن تتمدد هذه البقايا في ISM ذي كثافة متوسطة حجميا مقدارها n¯9(ηcr/0.4)1 cm-33.1، الجدول 3). في هذا القسم، نرغب في تحديد ما إذا كانت هذه البيئة متسقة مع نصف القطر (§2.1)، وكتلة الغبار (§2.2)، ودرجة حرارة الإلكترونات (§2.3) المرصودة من هذه البقايا. ويتيح لنا هذا التحليل أيضا استنتاج الخواص الأساسية، مثل عمرها $t_{\rm age}$ وسرعة تمددها الحالية $v_{\rm snr}$، اللازمة لفهم آلية تسريع الجسيمات الكامنة (مثلا، Reynolds & Keohane 1999; Reynolds 2008).

إن كمية الغبار 311 M في قشرة بقايا المستعر المستنتجة من أرصاد تحت الأحمر (§2.2؛ Lee et al. 2011) يغلب أن تهيمن عليها حبيبات غبار موجودة سلفا جرفتها المقذوفات المتمددة أو غبار تكون داخل البقايا. إذا كانت Kes 17 تتمدد في وسط ذي كثافة متوسطة حجميا $\bar{n} = 10\bar{n}_{10}$ cm-3 (حيث n¯101، §3.1)، فإن كتلة المادة المجروفة بواسطة مقذوفات المستعر الأعظم المتمددة $M_{\rm sw}$ هي:

Msw = 43πRsnr3n¯mp (2)
1300n¯10d103M. (3)

إذا كانت لهذا الوسط نسبة كتلة غبار إلى غاز نموذجية مقدارها $\sim0.1\% -
0.5\%$ (مثلا، Pei 1992)، فإن Kes 17 قد جرفت $\sim (1 -
5)\bar{n}_{10}~M_{\odot}$ من الغبار بين النجمي، وهي كمية قابلة للمقارنة مع الكتلة المقدرة من الأرصاد. وإذا كانت هناك كتلة أكبر بكثير عند درجات حرارة أدنى (Lee et al., 2011)، فإن ذلك لا يتطلب أن يكون الغبار الإضافي قد تكون داخل بقايا المستعر، بل يرجح أنه يدل على $\bar{n}_{10} > 1$ و/أو نسبة كتلة غبار إلى غاز أعلى، وهو أمر ممكن إذا كانت Kes 17 تتمدد في سحابة جزيئية. وفي الواقع، إن كشف ميزر OH (1720 MHz) (§2.1) وانبعاث الصدمات الجزيئية (§2.2) يشير إلى أن هذه البقايا تتمدد داخل سحابة جزيئية (Yusef-Zadeh et al., 2003). وهذا يعني أن الكتل المرصودة في تحت الأحمر هي على الأرجح نتيجة مادة كثيفة داخل السحابة جرفتها وصدمتها المقذوفات المتمددة (§2.2).

على الرغم من أن للسحب الجزيئية بنية كثافة معقدة جدا (مثلا، Williams et al. 1995)، يمكن تقريب هذه البيئة بوصفها مجموعة كتل ذات كثافة متوسطة $\bar{n}_{\rm clump}$ وعامل ملء حجمي $f_{\rm clump}$ مغمورة في وسط موحد بين الكتل ذو كثافة $n_{\rm ic}$ وعامل ملء حجمي $1-f_{\rm
  clump}$ (Chevalier, 1999). في هذا النموذج، تكون $\bar{n}$:

n¯ = n¯clumpfclump+nic(1fclump). (4)

كما لوحظ في §3.1، فإن كفاءة التسريع $\eta_{\rm
  cr} \propto \bar{n}^{-1}$، و$\eta_{\rm cr} \sim 10\%$ (القيمة المستنتجة لبقايا أخرى؛ الجدول 4) تتطلب $\bar{n} \sim 90$ cm-3. وقد قيست هذه المعاملات ($\bar{n}$ و$n_{\rm
  ic}$ و$\bar{n}_{\rm clump}$ و$f_{\rm clump}$) لعدد قليل من السحب الجزيئية، ووجدت أن القيم النموذجية عادة هي nic10cm3 و$f_{\rm clump}\sim10\%$ و$\bar{n} \sim 20~{\rm
  cm}^{-3}$ و$\bar{n}_{\rm clump} \sim 200 - 1000$ cm-3 مع تباين كبير بين السحب (مثلا، Blitz 1993; Williams et al. 1995). ومن تحليل الطيف الحراري للأشعة السينية لـ Kes 17 المقدم في §3.2.1 و§3.2.2، نقدر $n_{\rm ic} \lesssim 0.4~{\rm cm}^{-3}$. إذا كانت $\bar{n} \sim
9$ cm-3، و$\bar{n}_{\rm clump}=200$ cm-3، و$n_{\rm
  ic}\sim0.4$ cm-3، فإن $f_{\rm clump} \sim 4\%$ – وهو متسق مع القيم المرصودة. غير أنه إذا كانت $\bar{n} \sim 90$ cm-3، فإن $f_{\rm clump}$ قيمة عالية للغاية $\sim 45\%$ لهذه القيم من $\bar{n}_{\rm clump}$ و$n_{\rm ic}$. لكن إذا كانت $\bar{n}_{\rm
  clump} \sim 1000$ cm-3، كما قيس حول منطقة Hii NGC 2244 الموجودة داخل سحابة جزيئية (Williams et al., 1995)، فإن $f_{\rm clump}\sim 10\%$ لقيم $\bar{n} \sim 90$ cm-3 و$n_{\rm
  ic}\sim0.4$ cm-3. ويشير مدى كثافات الكتل (من $n_{\rm clump}
\sim 100 - 1000$ cm-3 حتى $n\sim10^4 - 10^6$ cm-3) المستنتج من تحليل طيف تحت الأحمر لـ Kes 17 (§2.2) إلى أن $\bar{n}_{\rm clump} \sim 1000$ cm-3 قيمة معقولة.

يتطلب تفسير نصف قطر Kes 17 ودرجة حرارتها السينية فهم تطورها الديناميكي. يقذف المستعر الأعظم مادة ذات كتلة $M_{\rm ej}$ وطاقة حركية ابتدائية $E_{\rm sn}$ إلى محيطه. في البداية، تتمدد المقذوفات بسرعة فوق صوتية بالنسبة لبيئتها، دافعة صدمة تسمى “الصدمة الأمامية” في محيطها. عند هذه الصدمة، تتسارع المادة المحيطة المجروفة وتنضغط وتسخن إلى ضغط أعلى بكثير من ضغط المقذوفات المتمددة. ونتيجة لذلك، تدفع المادة المحيطة المصدومة موجة صدمية، تسمى “الصدمة العكسية”، إلى داخل المقذوفات المتمددة، فتبطئها وتضغطها وتسخنها. في النموذج التطوري القياسي لبقايا المستعرات العظمى (مثلا، Chevalier 1977; Truelove & McKee 1999 والمراجع الواردة هناك)، عند الأزمنة المبكرة تتمدد الصدمة الأمامية بسرعة ثابتة تقريبا، بحيث يكون نصف قطر الصدمة الأمامية $R_{\rm snr} \propto t$. وبما أنها تتمدد بسرعة ثابتة، فإن المقذوفات تفقد قليلا من الطاقة الحركية خلال هذه المرحلة. وتنتهي مرحلة “التمدد الحر” هذه عندما تكون الصدمة العكسية قد مرت خلال جميع المقذوفات، وذلك تقريبا عندما تكون الكتلة المجروفة بواسطة الصدمة الأمامية $M_{\rm sw} \approx M_{\rm
  ej}$. خلال هذه المرحلة، التي يشار إليها عادة باسم “طور Sedov-Taylor”، يتوقع أن تهيمن الفواقد الأديباتية على تطور طاقة المقذوفات، وتتمدد بقايا المستعر وفق $R_{\rm snr} \propto
t^{2/5}$ (مثلا، Chevalier 1977; Truelove & McKee 1999 والمراجع الواردة هناك). وعندما يصبح زمن التبريد الإشعاعي للغاز المصدوم قابلا للمقارنة مع عمر بقايا المستعر، تهيمن الفواقد الإشعاعية، فتغير تطورها الديناميكي تغيرا كبيرا (Blondin et al., 1998).

ومع ذلك، سيكون تطور Kes 17 مختلفا اختلافا كبيرا بسبب تمددها في سحابة جزيئية متكتلة (مثلا، Chevalier 1999) وبسبب تسريعها الفعال للأشعة الكونية (مثلا، Ellison et al. 2007, 2010; Ferrand et al. 2010; Castro et al. 2011). والدليل الرصدي على أن Kes 17 تطورت بشكل مختلف هو طبيعتها ذات “المورفولوجيا المختلطة” (Combi et al., 2010)، المعرفة بالمزيج المرصود من قشرة راديوية حادة الطيف وانبعاث سيني حراري داخلي (Rho & Petre, 1998). ويمكن لملف الكثافة ودرجة الحرارة غير السيدوفي الذي توحي به مورفولوجيتها الحرارية السينية المملوءة مركزيا أن يعدل نمو بقايا المستعر (مثلا، White & Long 1991). وفي الوقت الراهن، فإن التفسيرين الفيزيائيين الرئيسيين لبقايا المستعرات ذات المورفولوجيا المختلطة هما أن التوصيل الحراري يدفع الغاز المسخن عند الصدمة الأمامية إلى المركز (مثلا، Cui & Cox 1992; Chevalier 1999) أو أن كتلا كثيفة تتبخر داخل البقايا (مثلا، White & Long 1991). ومع أن أيا من النموذجين لا يعيد بدقة إنتاج ملفات درجة الحرارة والسطوع السطحي السيني المرصودة لجميع بقايا المستعرات ذات المورفولوجيا المختلطة (مثلا، Slane et al. 2002)، فسنفسر نصف قطر ودرجة حرارة الإلكترونات في Kes 17 باستخدام هذه النماذج لتقدير عمرها وبيئتها تقديرا تقريبيا.

3.2.1 التوصيل الحراري

إذا كان التوصيل الحراري هو المسؤول الرئيس عن الأشعة السينية الحرارية المرصودة في مركز Kes 17، فمن المرجح أن يكون تطورها مشابها للتسلسل “القياسي” المبين أعلاه. وكما من قبل، نقرب بيئة السحابة الجزيئية بوصفها مجموعة من كتل منفصلة صغيرة وعالية الكثافة مغمورة في غاز منخفض الكثافة وثابت بين الكتل (Chevalier, 1999). ستصدم المقذوفات المتمددة كلا من الغاز بين الكتل والكتل، لكن بالنسبة إلى كثافات الكتل المستنتجة من أرصاد تحت الأحمر ($n_{\rm clump} \sim 10^2 -
10^6$ cm-3؛ §2.2)، تكون الصدمة المنقولة بطيئة جدا بحيث لا تسخن هذه المادة إلى درجات حرارة باعثة للأشعة السينية. لذلك ينبغي ألا تتجاوز كتلة الغاز الباعث للأشعة السينية في Kes 17 $M_{X}$ كتلة المادة بين الكتل المجروفة $M_{\rm ic}$ $(M_{X} \leq M_{\rm
  ic})$.

يمكن تقدير $M_{X}$ من ملاءمات الطيف الحراري للأشعة السينية المقدمة في §2.3. كتلة الغاز الباعث للأشعة السينية تساوي:

MX = 43πRsnr3fXnH,Xmp, (5)

حيث $m_p$ هي كتلة البروتون، و$n_{H,X}$ هي كثافة الغاز الباعث للأشعة السينية، و$f_X$ هي نسبة حجم بقايا المستعر المملوء بالبلازما الباعثة للأشعة السينية. ويمكننا تقدير $n_{H,X}$ من التطبيع $K$ لنماذج الانبعاث الحراري السيني المستخدمة في §2.3، إذ إن (Arnaud, 1996):

K = 10144πd2ne,XnH,X𝑑V (6)
= 0.3fXμ1nH,X2dθsnr3×1014 (7)

حيث $\mu$ هي النسبة العددية للإلكترونات إلى البروتونات في البلازما ($n_{H,X} = n_{e,X} / \mu$؛ $\mu=1.23$ للوفرة الشمسية)، و$d$ هي المسافة إلى المصدر بوحدة cm، و$\theta_{\rm snr}$ هو نصف القطر الزاوي لـ Kes 17 بالراديان. للقيم المقاسة من $K$ (الجدول 2، §2.3)، يكون $n_{H,X}$ تقريبا:

nH,X 0.4fX12d1012cm3. (8)

في بقايا مستعر أعظم قياسية، $f_X \approx 1/12$، لكن بقايا ذات مورفولوجيا مختلطة يرجح أن تمتلك $f_X$ أكبر من هذه القيمة. وبما أن $1/12 \leq
f_X \leq 1$، فإننا نقدر:

nH,X (0.41.4)d1012cm3. (9)

يتطلب ربط $M_X$ بـ $M_{\rm ic}$ تقدير زمن تبريد الغاز الباعث للأشعة السينية $t_{\rm cool}$. فإذا كانت Kes 17 أكبر عمرا بكثير من زمن التبريد، فإن $M_{X} \ll M_{\rm ic}$، أما إذا كانت Kes 17 أصغر عمرا من زمن التبريد، فإن $M_{X} \approx M_{\rm ic}$. ويمكن تقريب زمن التبريد على أنه:

tcool EthermalLthermal, (10)

حيث $E_{\rm thermal}$ هي الطاقة الحرارية لهذا الغاز و$L_{\rm
  thermal}$ هي لمعانه الحراري. وتكون الطاقة الحرارية تقريبا:

Ethermal = 43πRsnr3nH,XfXkBTX, (11)

حيث $k_B$ هو ثابت Boltzmann و$T_X$ هي درجة حرارة الأشعة السينية، في حين أن اللمعان الحراري $L_{\rm thermal}$ هو:

Lthermal = 43πRsnr3Λ (12)

حيث، بالنسبة إلى وفورات شمسية ودرجة حرارة $kT_X \sim 1~{\rm keV}$ (متوسط درجة حرارة الإلكترونات والتركيب الكيميائي التقريبي اللذان توحي بهما نمذجتنا؛ §2.3)، تكون $\Lambda$ (Raymond et al., 1976):

Λ 5ne,XnH,X×1023ergss1cm3 (13)
5μnH,X2×1023ergss1cm3. (14)

وهذا يشير إلى أن:

tcool kBTxfX5μnH,x×1023 (15)
2fX12×106years (16)
(0.62)×106years (17)

للمدى المذكور أعلاه من $f_X$. وهذا أعلى بكثير من عمر بقايا المستعرات الأخرى المحددة بوصفها مسرعات فعالة للأشعة الكونية (الجدول 4). لذلك، في هذا السيناريو، من المرجح أن $M_X \approx M_{\rm ic}$.

كتلة الغاز بين الكتل المجروفة بواسطة المقذوفات المتمددة هي:

Mic = 43πRsnr3(1fclump)nicmp. (18)

ولـ $M_X \approx M_{\rm ic}$، نحصل على:

nic fX1fclumpnH,X (19)
0.4fX12(1fclump)1d1012cm3 (20)

إذا كانت $f_{\rm clump}\sim10\%$ كما تشير الأرصاد (مثلا، Blitz 1993، Williams et al. 1995)، فإن المدى المسموح لـ $f_X$ يفضل $n_{\rm ic} \sim 0.1 - 0.4~{\rm cm}^{-3}$.

كما ذكر أعلاه، ينبغي في هذا السيناريو أن يكون تطور Kes 17 مشابها لتطور بقايا المستعر “القياسية”. لذلك فإن نصف قطرها المرصود $R_{\rm snr}$ يوحي بعمر (مثلا، Lozinskaya 1992):

tage = (Rsnr1.15)52(nicμmpEsn)12 (21)
4200d1052(nic0.4)12E5112years, (22)

حيث $E_{\rm sn}=10^{51}E_{51}$ ergs. للقيم $n_{\rm ic}$ المقدرة أعلاه، يكون عمر Kes 17 فقط $t_{\rm age} \sim 2000 - 4200$ سنة. وهذا أخفض بكثير من مقياس زمن التبريد $t_{\rm
  cool}$ المحسوب أعلاه، ومتسق مع افتراضنا أن $M_X
\approx M_{\rm ic}$. ومع أن تقدير العمر هذا يهمل أثر تسريع الأشعة الكونية، تشير الدراسات العددية إلى أن هذا التحليل يقلل العمر بمقدار $\sim20\%$ في حالة تسريع جسيمات فعال للغاية ($\eta_{\rm cr} \sim 40\%$؛ مثلا، Castro et al. 2011). لذلك نقدر في هذا السيناريو أن عمر Kes 17 هو $t_{\rm age} \sim 2000 - 5200$ سنة.

إذا كان ذلك صحيحا، فإن Kes 17 تتمدد حاليا بسرعة $v_{\rm snr}$ (مثلا، Lozinskaya 1992; Truelove & McKee 1999):

vsnr 0.43(Esnnicmp)12Rsnr32 (23)
570(niccm3)12kms. (24)

إذا كان $n_{\rm ic} \sim 0.1-0.4~{\rm cm}^{-3}$ كما استنتج أعلاه، فعندئذ $v_{\rm snr} \sim 900 - 1800~{\rm km~s}^{-1}$. ويتوقع من مثل هذه الصدمة أن تسخن الإلكترونات إلى درجة حرارة $T_e$:

kTe 316mevsnr2 (25)
0.93eV, (26)

وهي أدنى بكثير من $kT_{e} \sim 0.8~{\rm keV}$ المستنتجة من نمذجتنا لطيف الأشعة السينية المرصود (§2.3). غير أن الأيونات تسخن إلى درجة حرارة $T_i$:

kTi 316mpvsnr2 (27)
1.56keV, (28)

وهي أعلى من درجة حرارة الإلكترونات المرصودة. لكن آليات عديدة لتسخين الإلكترونات تعمل داخل بقايا المستعر الأعظم. وتشير الأرصاد إلى أنه، عند صدمة أمامية تتمدد بسرعة $v_{\rm snr} \sim 900 - 1800~{\rm
  km~s}^{-1}$، ستسخن الإلكترونات إلى درجة حرارة $T_e \sim
(0.1-0.8)T_i$ عند الصدمة الأمامية لهذا المدى من السرعات الصدمية (Ghavamian et al., 2007)، وربما تعزز هذه العملية بتسريع جسيمات فعال عند الصدمة الأمامية (Castro et al., 2011). إضافة إلى ذلك، تسخن الأيونات الإلكترونات داخل بقايا المستعر عبر تصادمات Coulomb. ويوحي مقياس زمن التأين العالي المستنتج من نمذجتنا لطيف الأشعة السينية الحراري المرصود بوجود اتزان حراري تقريبي على الأقل بين الإلكترونات والأيونات في هذه البقايا. لذلك يتسق هذا السيناريو مع درجة حرارة الإلكترونات المرصودة.

خلاصة القول، إذا كان التوصيل الحراري هو الآلية المهيمنة المسؤولة عن الطبيعة ذات المورفولوجيا المختلطة لـ Kes 17، فإن هذه البقايا عمرها $\sim2000-5000$ سنة وتتمدد في وسط متكتل بكثافة بين-كتلية مقدارها $\sim0.1-0.4~{\rm cm}^{-3}$.

3.2.2 تبخر الكتل

بدلا من ذلك، قد تنجم الطبيعة ذات المورفولوجيا المختلطة لـ Kes 17 عن تبخر الكتل الكثيفة التي جرفتها المقذوفات المتمددة داخل البقايا ثم تسخينها إلى درجات حرارة باعثة للأشعة السينية بواسطة الغاز الساخن بين الكتل الذي صدم عند الصدمة الأمامية. في هذه الحالة، نتوقع $M_X > M_{\rm ic}$. وبتكرار تحليل §3.2.1، يتطلب ذلك أن:

nic 0.4fX12(1fclump)1d1012cm3. (29)

وبما أن $f_X < 1$، ففي هذا السيناريو $n_{\rm ic} < 0.4~{\rm cm}^{-3}$.

يمكن لتبخر الكتل داخل بقايا المستعر أن يؤثر تأثيرا كبيرا في ديناميكيات الصدمة الأمامية، التي يتوقع الآن أن تتمدد وفق:

Rsnr = [25(γ+1)κEsn16πnicmp]15t25, (30)

حيث $\gamma=5/3$ هو معامل أديباتية المادة المحيطة و$\kappa$ هو نسبة الطاقة الحرارية إلى الحركية في بقايا المستعر (White & Long, 1991). لذلك يكون عمر Kes 17:

tage = [25(γ+1)16πmp]12(κEsnnic)12Rsnr52 (31)
6000κ12E5112nic12years (32)

لـ $d_{10}=1$، حيث يمكن أن تتغير $\kappa$ بين 0.01 و1 (White & Long, 1991). ويشير اختيار $\kappa = 0.01$ و$n_{\rm ic} = 0.4~{\rm
  cm}^{-3}$ إلى أن عمرها $t_{\rm age} < 40000$ سنة. ومع أن هذا التحليل يهمل أثر تسريع الجسيمات الفعال في تطور الصدمة الأمامية (White & Long, 1991)، فإن الخطأ الناتج $\sim20\%$ الذي توحي به المحاكاة (مثلا، Castro et al. 2011) أقل بكثير من الارتياب الناتج عن القيم المجهولة لـ $\kappa$ و$n_{\rm ic}$.

في هذا السيناريو، تعطى $v_{\rm snr}$ بـ:

vsnr = 25[25(γ+1)κEsn16πnicmp]15tage35, (33)
120,000(κE51nic)15(tage1year)35kms. (34)

وبالنظر إلى حالة “العمر الأقصى” أعلاه ($\kappa=0.01$، و$n_{\rm
  ic}=0.4$ cm-3، و$t_{\rm age}=40000$ سنة)، نجد أن حاليا $v_{\rm snr} \geq 100$ km s-1. ومن المعادلتين 25 و27، تكون هذه السرعة الدنيا بطيئة جدا بحيث لا تسخن حاليا الإلكترونات أو الأيونات إلى درجة حرارة الإلكترونات المقاسة $kT_{e}\sim0.8~{\rm keV}$2.3). ويمكن تصحيح ذلك إما بأن تكون Kes 17 أصغر عمرا من العمر الأقصى المقدر أعلاه، مما يوحي بـ $\kappa > 0.01$ و/أو $n_{\rm ic} <
0.4~{\rm cm}^{-3}$، أو بأن المادة المشعة قد سخنت في وقت سابق عندما كانت البقايا تتمدد أسرع. وهذا معقول لأن الزمن اللازم لكي يبرد الغاز الباعث للأشعة السينية ($t_{\rm cool}
\gtrsim 6\times10^5~{\rm years}$؛ §3.2.1، المعادلة 15) أطول من العمر الأقصى لـ Kes 17 في هذا السيناريو. لذلك فمن المعقول أن المادة التي سخنت في أزمنة أبكر ما زالت تشع اليوم.

خلاصة القول، إذا كان تبخر الكتل هو السبب المهيمن للطبيعة ذات المورفولوجيا المختلطة لـ Kes 17، فإنها تتمدد في وسط ذي كثافة بين-كتلية $n_{\rm ic} < 0.4~{\rm cm}^{-3}$ وعمرها $<40000$ سنة.

4. الاستنتاجات

في هذه الورقة، نحلل ونفسر أرصادا حديثة لبقايا المستعر الأعظم Kes 17 عبر الطيف الكهرومغناطيسي (§2). يشير تحليلنا إلى أن لهذه البقايا قشرة راديوية جزئية بقطر $\sim7\farcm5$، وهو ما يترجم إلى نصف قطر فيزيائي مقداره $R\sim10d_{10}$ pc عند مسافة $d=10d_{10}$ kpc. ويوحي كشف انبعاث ميزر OH عند 1720 MHz وطيف تحت الأحمر لـ Kes 17 بأن هذه البقايا تتمدد في سحابة جزيئية، كما يوحي تحليلنا لرصد حديث لهذه البقايا بواسطة Suzaku بأن الانبعاث السيني المرصود حراري في معظمه، تصدره بلازما ذات كثافة $n_{H,X} \sim 0.4$ cm-3، وبوفورات قريبة من الشمسية باستثناء نقص في وفرة S، وتتألف من غاز ذي درجة حرارة إلكترونية $kT_{\rm e} \approx 0.8$ keV في اتزان حراري تقريبي (§2.3 و§3.2؛ الجدول 2). وأخيرا، يكشف تحليلنا لأرصاد Fermi لهذا الحقل بقوة انبعاث أشعة $\gamma$ في نطاق GeV متوافقا موضعيا مع هذه البقايا ومن شبه المؤكد أنه صادر من قشرة بقايا المستعر (§2.4).

من خلال نمذجة الانبعاث غير الحراري عريض النطاق لـ Kes 17، حددنا أن أشعة $\gamma$ في نطاق GeV تنتج في الغالب من بروتونات الأشعة الكونية المسرعة عند الصدمة الأمامية لبقايا المستعر وهي تصطدم بمادة مجروفة داخل البقايا، منتجة $\pi^0$s تضمحل إلى أشعة $\gamma$3.1). ويتطلب هذا التفسير أن تتمدد هذه البقايا في وسط بكثافة متوسطة n¯9 cm-33.1)، وهو ما يتسق مع بيئة السحابة الجزيئية التي يدل عليها انبعاث ميزر OH وكتلة الغبار الكبيرة المستنتجة من طيفها تحت الأحمر (§3.2). ويعتمد عمر Kes 17 وكثافة الوسط بين الكتل داخل السحابة $n_{\rm ic}$ على ما إذا كان التوصيل الحراري أو تبخر الكتل الكثيفة هو المسؤول الرئيس عن الطبيعة ذات المورفولوجيا المختلطة لهذه البقايا. فإذا كان التوصيل الحراري هو المسؤول، فإن Kes 17 تتمدد في بيئة ذات $n_{\rm ic} \sim 0.1-0.4~{\rm cm}^{-3}$ وعمرها ليس إلا $\sim2000-5000$ سنة (§3.2.1). أما إذا كان تبخر الكتل الكثيفة هو المسؤول الرئيس، فإن Kes 17 تتمدد داخل بيئة ذات $n_{\rm ic} < 0.4~{\rm cm}^{-3}$ وقد يبلغ عمرها $\sim40000$ سنة (§3.2.2). وإذا كانت كفاءة الأشعة الكونية $\eta_{\rm cr}$ لـ Kes 17 مشابهة لتلك في بقايا مستعرات أخرى يعتقد أنها تسرع الأشعة الكونية (الجدول 4)، فإنها تتمدد في بيئة ذات $\bar{n}\sim90$ cm-33.2). ويتطلب ذلك أن تكون للكتل المحيطة كثافة متوسطة مقدارها $\bar{n}_{\rm
  clump}\sim1000$ cm-3 من أجل كسر كتلة معقول للكتل $f_{\rm clump} \sim 10\%$3.2). وقد رصدت مثل هذه الكثافة المتوسطة للكتل في بعض السحب الجزيئية (Blitz, 1993)، وكذلك حول نجوم ضخمة كونت فقاعة رياح نجمية أو منطقة Hii داخل سحابة جزيئية (Williams et al., 1995).

إن كفاءة تسريع الأشعة الكونية المرتفعة المحتملة المستنتجة في Kes 17 مثيرة جدا للاهتمام – ولا سيما بالنظر إلى تمددها المرجح في وسط متكتل. وما زال كثير من آلية تسريع الجسيمات داخل بقايا المستعرات العظمى، ولا سيما كيفية اعتماد كفاءتها على محيطها، غير معروف. وتشير أعمال نظرية حديثة إلى أن التمدد في بيئة مضطربة ومتكتلة وقوية المغنطة يعزز تسريع الأشعة الكونية (مثلا، Bykov et al. 2000; Zhang et al. 2009)، وستختبر دراسة إضافية لـ Kes 17 هذه النتائج. ويتطلب ذلك فهما أفضل لبيئة Kes 17، وتحديدا مدى تكتل محيطها. كما توجد حاجة إلى أرصاد إضافية لقياس خواص الأشعة الكونية المسرعة من أجل اختبار نماذج آلية التسريع. ويمكن إنجاز ذلك بطرائق متنوعة؛ فمثلا، سيسمح قياس طيف أشعة $\gamma$ في نطاق TeV بتحديد الطاقة العظمى للأشعة الكونية المسرعة في بقايا المستعر $E_{0p}$، كما سيتيح رسم خرائط انبعاث CO حول Kes 17 قياس الكثافة المتوسطة $\bar{n}$ لبيئتها (مثلا، Williams et al. 1995)، وسيتيح رصد أعمق بالأشعة السينية لنا استخدام الانبعاث الحراري السيني المرصود (مثلا، Ellison et al. 2010) لتقييد كفاءة تسريع الأشعة الكونية $\eta_{\rm cr}$، وتقييد أفضل لطاقة أي إلكترونات أشعة كونية مسرعة في هذه البقايا، وتحديد منشأ مظهرها ذي المورفولوجيا المختلطة. وعلى أي حال، فإن الدراسة المستقبلية لـ Kes 17 مهمة للغاية لفهم كيفية تسريع الجسيمات النشطة في بقايا المستعرات العظمى.

يشكر المؤلفون الحكم المجهول على تعليقاته واقتراحاته المفيدة. يقر JDG بدعم منحة NASA NNX10AR51G، وزمالة NSF لما بعد الدكتوراه في علم الفلك والفيزياء الفلكية AST-0702957، وكذلك بضيافة Center for Cosmology and Particle Physics في New York University حيث أنجز قدر كبير من هذا العمل. ويقر POS بالدعم من عقد NASA NAS8-03060. يمول Commonwealth of Australia تلسكوب Australia Telescope لتشغيله مرفقا وطنيا تديره CSIRO. وقد استخدم هذا البحث بيانات حصل عليها من القمر Suzaku، وهي بعثة مشتركة بين وكالتي الفضاء في اليابان (JAXA) و الولايات المتحدة الأمريكية (NASA).

References

  • Abdo et al. (2010a) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010a, ApJ, 710, L92
  • Abdo et al. (2010b) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010b, ApJ, 718, 348
  • Abdo et al. (2009a) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009a, Physical Review Letters, 103, 251101
  • Abdo et al. (2009b) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009b, Phys. Rev. D, 80, 122004
  • Abdo et al. (2010c) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010c, Science, 327, 1103
  • Abdo et al. (2010d) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010d, ApJ, 712, 459
  • Ackermann et al. (2013) Ackermann, M., Ajello, M., Allafort, A., et al. 2013, Science, 339, 807
  • Aleksić et al. (2012) Aleksić, J., Alvarez, E. A., Antonelli, L. A., et al. 2012, A&A, 541, A13
  • Anders & Grevesse (1989) Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197
  • Arendt (1989) Arendt, R. G. 1989, ApJS, 70, 181
  • Arnaud (1996) Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17
  • Arnett & Schramm (1973) Arnett, W. D. & Schramm, D. N. 1973, ApJ, 184, L47+
  • Baring et al. (1999) Baring, M. G., Ellison, D. C., Reynolds, S. P., et al. 1999, ApJ, 513, 311
  • Berezhko et al. (2003) Berezhko, E. G., Pühlhofer, G., & Völk, H. J. 2003, A&A, 400, 971
  • Berezhko & Völk (2004) Berezhko, E. G. & Völk, H. J. 2004, A&A, 419, L27
  • Blackburn (1995) Blackburn, J. K. 1995, in Astronomical Society of the Pacific Conference Series, Vol. 77, Astronomical Data Analysis Software and Systems IV, ed. R. A. Shaw, H. E. Payne, & J. J. E. Hayes, 367
  • Blitz (1993) Blitz, L. 1993, in Protostars and Planets III, ed. E. H. Levy & J. I. Lunine, 125–161
  • Blondin et al. (1998) Blondin, J. M., Wright, E. B., Borkowski, K. J., et al. 1998, ApJ, 500, 342
  • Bykov et al. (2000) Bykov, A. M., Chevalier, R. A., Ellison, D. C., et al. 2000, ApJ, 538, 203
  • Cassam-Chenaï et al. (2007) Cassam-Chenaï, G., Hughes, J. P., Ballet, J., et al. 2007, ApJ, 665, 315
  • Castro & Slane (2010) Castro, D. & Slane, P. 2010, ApJ, 717, 372
  • Castro et al. (2011) Castro, D., Slane, P., Patnaude, D. J., et al. 2011, ApJ, 734, 85
  • Caswell et al. (1975) Caswell, J. L., Murray, J. D., Roger, R. S., et al. 1975, A&A, 45, 239
  • Chevalier (1977) Chevalier, R. A. 1977, ARA&A, 15, 175
  • Chevalier (1999) —. 1999, ApJ, 511, 798
  • Combi et al. (2010) Combi, J. A., Albacete Colombo, J. F., Sánchez-Ayaso, E., et al. 2010, A&A, 523, A76+
  • Cui & Cox (1992) Cui, W. & Cox, D. P. 1992, ApJ, 401, 206
  • Drury et al. (1994) Drury, L. O., Aharonian, F. A., & Voelk, H. J. 1994, A&A, 287, 959
  • Elitzur (1976) Elitzur, M. 1976, ApJ, 203, 124
  • Ellison et al. (2007) Ellison, D. C., Patnaude, D. J., Slane, P., et al. 2007, ApJ, 661, 879
  • Ellison et al. (2010) Ellison, D. C., Patnaude, D. J., Slane, P., et al. 2010, ApJ, 712, 287
  • Eriksen et al. (2011) Eriksen, K. A., Hughes, J. P., Badenes, C., et al. 2011, ApJ, 728, L28
  • Ferrand et al. (2010) Ferrand, G., Decourchelle, A., Ballet, J., Teyssier, R., & Fraschetti, F. 2010, A&A, 509, L10
  • Frail et al. (1996) Frail, D. A., Goss, W. M., Reynoso, E. M., et al. 1996, AJ, 111, 1651
  • Freeman et al. (2001) Freeman, P., Doe, S., & Siemiginowska, A. 2001, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 4477, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. J.-L. Starck & F. D. Murtagh, 76–87
  • Ghavamian et al. (2007) Ghavamian, P., Laming, J. M., & Rakowski, C. E. 2007, ApJ, 654, L69
  • Giordano et al. (2012) Giordano, F., Naumann-Godo, M., Ballet, J., et al. 2012, ApJ, 744, L2
  • Giuliani et al. (2011) Giuliani, A., Cardillo, M., Tavani, M., et al. 2011, ApJ, 742, L30
  • Giuliani et al. (2010) Giuliani, A., Tavani, M., Bulgarelli, A., et al. 2010, A&A, 516, L11
  • Gok & Sezer (2012) Gok, F. & Sezer, A. 2012, ArXiv e-prints
  • Hewitt et al. (2009) Hewitt, J. W., Rho, J., Andersen, M., et al. 2009, ApJ, 694, 1266
  • Kamae et al. (2006) Kamae, T., Karlsson, N., Mizuno, T., et al. 2006, ApJ, 647, 692
  • Kim et al. (2008) Kim, Y., Rieke, G. H., Krause, O., et al. 2008, ApJ, 678, 287
  • Koo et al. (2010) Koo, B.-C., Heiles, C., Stanimirović, S., et al. 2010, AJ, 140, 262
  • Koo et al. (1995) Koo, B.-C., Kim, K.-T., & Seward, F. D. 1995, ApJ, 447, 211
  • Lee (2005) Lee, H.-G. 2005, Journal of Korean Astronomical Society, 38, 385
  • Lee et al. (2011) Lee, H.-G., Moon, D.-S., Koo, B.-C., et al. 2011, ApJ, 740, 31
  • Lozinskaya (1992) Lozinskaya, T. A. 1992, Supernovae and stellar wind in the interstellar medium, ed. Lozinskaya, T. A.
  • Mattox et al. (1996) Mattox, J. R., Bertsch, D. L., Chiang, J., et al. 1996, ApJ, 461, 396
  • McClure-Griffiths et al. (2005) McClure-Griffiths, N. M., Dickey, J. M., Gaensler, B. M., et al. 2005, ApJS, 158, 178
  • McKee & Ostriker (1977) McKee, C. F. & Ostriker, J. P. 1977, ApJ, 218, 148
  • Milne & Dickel (1975) Milne, D. K. & Dickel, J. R. 1975, Australian Journal of Physics, 28, 209
  • Mori (2009) Mori, M. 2009, Astroparticle Physics, 31, 341
  • Morlino & Caprioli (2012) Morlino, G. & Caprioli, D. 2012, A&A, 538, A81
  • Nolan et al. (2012) Nolan, P. L., Abdo, A. A., Ackermann, M., et al. 2012, ApJS, 199, 31
  • Pei (1992) Pei, Y. C. 1992, ApJ, 395, 130
  • Pinheiro Gonçalves et al. (2011) Pinheiro Gonçalves, D., Noriega-Crespo, A., Paladini, R., et al. 2011, AJ, 142, 47
  • Rando & for the Fermi LAT Collaboration (2009) Rando, R. & for the Fermi LAT Collaboration. 2009, ArXiv e-prints
  • Raymond et al. (1976) Raymond, J. C., Cox, D. P., & Smith, B. W. 1976, ApJ, 204, 290
  • Raymond & Smith (1977) Raymond, J. C. & Smith, B. W. 1977, ApJS, 35, 419
  • Reach et al. (2005) Reach, W. T., Rho, J., & Jarrett, T. H. 2005, ApJ, 618, 297
  • Reach et al. (2006) Reach, W. T., Rho, J., Tappe, A., et al. 2006, AJ, 131, 1479
  • Reynolds (2008) Reynolds, S. P. 2008, ARA&A, 46, 89
  • Reynolds & Keohane (1999) Reynolds, S. P. & Keohane, J. W. 1999, ApJ, 525, 368
  • Rho & Petre (1998) Rho, J. & Petre, R. 1998, ApJ, 503, L167+
  • Salpeter (1977) Salpeter, E. E. 1977, ARA&A, 15, 267
  • Sault et al. (1995) Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Society of the Pacific Conference Series, Vol. 77, Astronomical Data Analysis Software and Systems IV, ed. R. A. Shaw, H. E. Payne, & J. J. E. Hayes, 433–+
  • Shaver & Goss (1970) Shaver, P. A. & Goss, W. M. 1970, Australian Journal of Physics Astrophysical Supplement, 14, 133
  • Slane et al. (2002) Slane, P., Smith, R. K., Hughes, J. P., et al. 2002, ApJ, 564, 284
  • Smith & Hughes (2010) Smith, R. K. & Hughes, J. P. 2010, ApJ, 718, 583
  • Strong et al. (2000) Strong, A. W., Moskalenko, I. V., & Reimer, O. 2000, ApJ, 537, 763
  • Tang et al. (2011) Tang, Y. Y., Fang, J., & Zhang, L. 2011, ApJ, 739, 11
  • Tavani et al. (2010) Tavani, M., Giuliani, A., Chen, A. W., et al. 2010, ApJ, 710, L151
  • Troja et al. (2008) Troja, E., Bocchino, F., Miceli, M., & Reale, F. 2008, A&A, 485, 777
  • Truelove & McKee (1999) Truelove, J. K. & McKee, C. F. 1999, ApJS, 120, 299
  • Uchiyama et al. (2010) Uchiyama, Y., Blandford, R. D., Funk, S., et al. 2010, ApJ, 723, L122
  • Uchiyama et al. (2012) Uchiyama, Y., Funk, S., Katagiri, H., et al. 2012, ApJ, 749, L35
  • Wardle & McDonnell (2012) Wardle, M. & McDonnell, K. 2012, ArXiv e-prints
  • White & Long (1991) White, R. L. & Long, K. S. 1991, ApJ, 373, 543
  • Whiteoak & Green (1996) Whiteoak, J. B. Z. & Green, A. J. 1996, A&AS, 118, 329
  • Williams et al. (1995) Williams, J. P., Blitz, L., & Stark, A. A. 1995, ApJ, 451, 252
  • Wu et al. (2011) Wu, J. H. K., Wu, E. M. H., Hui, C. Y., et al. 2011, ApJ, 740, L12
  • Yoshida (2008) Yoshida, K. 2008, Advances in Space Research, 42, 477
  • Yusef-Zadeh et al. (2003) Yusef-Zadeh, F., Wardle, M., Rho, J., et al. 2003, ApJ, 585, 319
  • Zhang et al. (2009) Zhang, W., MacFadyen, A., & Wang, P. 2009, ApJ, 692, L40
SNR Age n¯ ISM Bsnr ηcr Kep E0,p Citations
Kes 17 2,00040,000 >9 cm-3 Clumpy 35 μG <0.4 0.02 >500 GeV
Cas A 330 years 30 cm-3 Clumpy 0.51 mG 0.0050.02 0.0040.02 1030 TeV aaBerezhko et al. (2003), bbBerezhko & Völk (2004), ccAbdo et al. (2010a), vvKim et al. (2008)
Tycho 440 years 0.3 cm-3 Uniform 200300μG 0.060.075 0.0016 >470 TeV hhCassam-Chenaï et al. (2007), iiEriksen et al. (2011), jjGiordano et al. (2012), kkMorlino & Caprioli (2012)
IC 443 4,000 years 250 cm-3 Clumpy 10 μG 0.0060.02 0.010.03 100200 GeV ddTroja et al. (2008), eeAbdo et al. (2010d), ffTavani et al. (2010), ggTang et al. (2011)
W44 20,000 years 100 cm-3 Clumpy 40800μG 0.030.15 0.010.05 nnReach et al. (2005), ooAbdo et al. (2010c), ppUchiyama et al. (2010), qqGiuliani et al. (2011), ggTang et al. (2011), rrUchiyama et al. (2012)
W51C 30,000 years 10 cm-3 Clumpy <150μG 0.16 0.0125 120 TeV ssKoo et al. (1995), ttKoo et al. (2010), ggTang et al. (2011), uuAleksić et al. (2012)
W28 40,000 years 100 cm-3 Clumpy 40 – 160 μG 0.01 – 0.03 0.01 llAbdo et al. (2010b), mmGiuliani et al. (2010), ggTang et al. (2011)
Table 4الخواص الفيزيائية لبقايا مستعرات عظمى غير Kes 17 تمتلك دليلا رصديا مباشرا على تسريع البروتونات. $\bar{n}$ هي الكثافة المتوسطة لـ ISM المحيط، و$B_{\rm
    snr}$ هي شدة المجال المغناطيسي داخل بقايا المستعر، و$\eta_{\rm cr}$ هي النسبة بين $E_{\rm cosmic ray}$ والطاقة الحركية الابتدائية للمستعر الأعظم السلف (المفترض في حالات كثيرة أنها $10^{51}$ ergs)، و$K_{ep}$ هي التطبيع النسبي بين الإلكترونات والبوزيترونات المسرعة، و$E_{0,p}$ هي طاقة القطع في طيف البروتونات المسرعة. وتعكس المجالات المعطاة لقيم مختلفة اختلافات في الأدبيات، وفي بعض بقايا المستعرات تكون القيم المبلغ عنها افتراضات مستخدمة في النمذجة لا قيما ملائمة. لا تعطى قيمة لـ $E_{0,p}$ لبقايا المستعرات W28 وW44 لأن طيف حقن الأشعة الكونية بقانون قوة مكسور، بدلا من قانون قوة بقطع أسي، مطلوب لإعادة إنتاج أطيافها غير الحرارية.