Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. II. Propagation of Primary and Secondary nuclei in a magneto-hydrodynamical environment

1. بيان المشكلة

يتم نمذجة انتشار الأشعة الكونية (CR) تقليديًا باستخدام نماذج "الصندوق المتسرب" (leaky-box) الظواهرية أو معادلات الانتشار-الحمل المستقرة (مثل GALPROP) التي تفترض وجود وسط بين نجمي (ISM) ثابت. ومع ذلك، غالبًا ما تتجاهل هذه النماذج التغذية الراجعة الديناميكية للأشعة الكونية على الوسط بين النجمي، مثل الرياح المدفوعة بالأشعة الكونية وإعادة هيكلة المجال المغناطيسي. وبينما توجد نماذج هيدروديناميكية مغناطيسية (MHD) ذاتية الاتساق، إلا أنها غالبًا ما تعامل الأشعة الكونية كـ "سائل رمادي" غير محدد طيفيًا أو تركز فقط على الإلكترونات.

هناك حاجة إلى نموذج ذاتي الاتساق ومحدد طيفيًا يتتبع في آن واحد انتشار أنواع متعددة من نوى الأشعة الكونية الأولية والثانوية مع ربطها ديناميكيًا بالغاز الحراري، والمجالات المغناطيسية، وتغذية المستعرات العظمى (SN) الراجعة. وبشكل محدد، يهدف المؤلفون إلى التحقق من إنتاج وانتشار النوى الثانوية (Li, Be, B) الناتجة عن تفتت (spallation) النوى الأولية (C, N, O) داخل بيئة مجرية ديناميكية وواقعية.

2. المنهجية

الإطار العددي: كود PIERNIK MHD

تستخدم الدراسة كود PIERNIK، وهو كود MHD متعدد السوائل يعتمد على الشبكة. قام المؤلفون بتوسيع خوارزمية CRES (طيف طاقة الأشعة الكونية) الموجودة مسبقًا، والتي كانت تُستخدم للإلكترونات، لتشمل الأنواع الهادرونية للأشعة الكونية.

• نهج العزوم المزدوجة: يحل الكود معادلة فوكر-بلانك لدوال التوزيع باستخدام نهج العزوم المزدوجة (كثافة العدد $n_a$ وكثافة الطاقة الحركية $e_a$) ضمن فئات الزخم.
• تقريب قانون القوة المجزأ: داخل كل فئة زخم، يتم تقريب دالة التوزيع كقانون قوة. يسمح هذا بالتكامل العددي الفعال لمصطلحات كسب/فقد الطاقة دون الحاجة إلى حل دالة التوزيع الكاملة عند كل خطوة.
• الدقة الطيفية: يتتبع النموذج 16 فئة زخم تتراوح من الطاقات غير النسبية إلى النسبية الفائقة ($5 \times 10^{-1}$ إلى $10^5$ من كتلة البروتون).

العمليات الفيزيائية المنفذة

• النوى الأولية: يتم حقن $^{12}\text{C}$ و $^{14}\text{N}$ و $^{16}\text{O}$ عبر بقايا المستعرات العظمية (SNRs). يتم التعامل مع البروتونات كـ "سائل رمادي" (غير محدد طيفيًا) نسبي ليدفع ديناميكيات الوسط بين النجمي.
• النوى الثانوية: يتم توليد الأنواع الثانوية ($^{7}\text{Li}$، $^{9}\text{Be}$، $^{10}\text{Be}$، $^{10}\text{B}$، $^{11}\text{B}$) عبر التفتت (التصادمات الهادرونية) للنوى الأولية ضد بروتونات الوسط بين النجمي الحرارية. يتم أخذ المقاطع العرضية من GALPROP.
• التحلل الإشعاعي: تم تضمين تحلل النظائر غير المستقرة (تحديدًا $^{10}\text{Be}$).
• فقدان الطاقة: يتضمن النموذج خسائر التمدد الأديباتي وتبريد كولوم (باستخدام صيغة بيث-بلوك). تم استبعاد تأثيرات التدفق (streaming) حاليًا.
• آليات النقل: تخضع الأشعة الكونية لعمليات الحمل مع الغاز الحراري والانتشار غير المتماثل المعتمد على الصلابة ($D_{\parallel} \propto R^{\delta}$) على طول خطوط المجال المغناطيسي.

إعداد المحاكاة

• النطاق: صندوق ثلاثي الأبعاد متدرج جاذبيًا ($0.5 \times 1 \times 8$ كيلو بارسيك) يمثل الوسط بين النجمي المحلي بالقرب من النظام الشمسي.
• الشروط الأولية: توازن هيدروستاتيكي مع $n_{\text{ISM}} \approx 1 \text{ cm}^{-3}$ و $T \approx 7000$ كلفن.
• التغذية الراجعة: تقوم المستعرات العظمية الموزعة عشوائيًا بحقن بروتونات الأشعة الكونية (10% من $10^{51}$ إرج لكل مستعر أعظم) لدفع النظام.
• دراسة المعلمات: أجرى المؤلفون محاكاة بتغيير ثلاث معلمات رئيسية:
  1. معدل المستعر الأعظم (SN Rate): 20، 60، و 80 كيلو بارسيك/ميجا سنة.
  2. معامل الانتشار ($D_{\parallel}$): $3 \times 10^{28}$ و $6 \times 10^{28}$ سم$^2$/ثانية.
  3. مؤشر الصلابة ($\delta$): 0.3 و 0.5.

3. المساهمات الرئيسية

1. توسيع CRESP ليشمل الأنواع الهادرونية: يعد هذا أول تنفيذ لخوارزمية CRESP في PIERNIK للتعامل مع أنواع متعددة من النوى الأولية والثانوية المحددة طيفيًا في آن واحد، بما في ذلك مصطلحات مصدر التفتت والتحلل الإشعاعي.
2. اقتران ذاتي الاتساق بين الأشعة الكونية والوسط بين النجمي: على عكس النماذج المستقرة، توضح هذه الدراسة كيف تشكل الأشعة الكونية ديناميكيًا الوسط بين النجمي (عبر الطفو والرياح) وكيف يقوم هذا الوسط الذي أعيد تشكيله بدوره بتغيير نقل الأشعة الكونية.
3. اكتشاف جديد حول الانتشار ونسبة البورون إلى الكربون (B/C): يحدد البحث علاقة غير متوقعة بين معامل الانتشار ونسبة البورون إلى الكربون (B/C)، مدفوعة بإعادة هيكلة المجال المغناطيسي.

4. النتائج الرئيسية

التطور الديناميكي

• التدفقات المدفوعة بالأشعة الكونية: تؤدي حقن الأشعة الكونية إلى دفع تدفقات غازية رأسية وإنشاء هياكل عدم استقرار باركر (حلقات مغناطيسية). تؤدي معدلات المستعر الأعظم الأعلى إلى رياح أقوى وسرعات غاز رأسية أسرع.
• إعادة هيكلة المجال المغناطيسي: يخلق طفو الأشعة الكونية مكونات مجال مغناطيسي رأسية. تعتمد قوة هذه المجالات الرأسية على كفاءة الانتشار.

الخصائص الطيفية

• انحرافات الميل: تظهر أطياف الأشعة الكونية المتطورة انحرافات عن ميل الحقن الأولي ($q \approx 4.1$).
  • عند الطاقات النسبية، يطابق الفرق في ميل الطيف بين الأولي والثانوي ($\Delta a_R$) مؤشر الانتشار $\delta$.
  • عند الطاقات غير النسبية، تكون الميول أكثر حدة بسبب تبريد كولوم والتفتت، حيث يظهر القرص أطيافًا أكثر صلابة من الهالة.
• قمم كثافة الطاقة: تظهر الأطياف قممًا حول 10–30 جيجا إرج/نوكليون بسبب التفاعل بين الحقن، والتبريد، والانتشار.

الملاحظات: نسب B/C و $^{10}\text{Be}/^{9}\text{Be}$

• تأثير معدل المستعر الأعظم:
  • معدل مستعر أعظم أعلى $\leftarrow$ نسبة B/C أقل. تزيد الرياح القوية من عملية الحمل، مما يقلل من زمن بقاء النوى الأولية في القرص الكثيف حيث يحدث التفتت.
  • تؤدي معدلات المستعر الأعظم الأعلى أيضًا إلى نسب $^{10}\text{Be}/^{9}\text{Be}$ أقل بسبب الهروب والتحلل الأسرع.
• تأثير معامل الانتشار ($D_{\parallel}$):
  • اكتشاف جديد: خلافًا لتوقعات نموذج الصندوق المتسرب القياسي (حيث يعني $D$ الأعلى B/C أقل)، يظهر هذا النموذج أن $D_{\parallel}$ الأعلى يؤدي إلى B/C أعلى.
  • الآلية: يوزع الانتشار الأعلى الأشعة الكونية بشكل أكثر توازنًا على طول خطوط المجال الأفقي، مما يقلل من تدرجات ضغط الأشعة الكونية المحلية التي تدفع الطفو المغناطيسي. ينتج عن ذلك مجالات مغناطيسية رأسية أضعف ($\langle |B_z|/|\vec{B}| \rangle$ ينخفض). تعمل المجالات الرأسية الأضعف على حبس النوى الأولية لفترة أطول في القرص، مما يزيد من إنتاج النوى الثانوية (البورون).
• تأثير مؤشر الصلابة ($\delta$):
  • يؤدي زيادة $\delta$ من 0.3 إلى 0.5 إلى جعل نسبة B/C أكثر حدة عند الطاقات النسبية، مما يجعل نتائج المحاكاة في توافق أفضل مع بيانات رصد AMS-02 (التي تفضل $\delta \approx 1/3$).

5. الأهمية والقيود

الأهمية:

• تثبت الدراسة أن محاكاة MHD ذاتية الاتساق ضرورية لتفسير ملاحظات الأشعة الكونية بدقة. إن "الانتشار الفعال" الذي تختبره الأشعة الكونية ليس مجرد دالة لمعلمة $D_{\parallel}$، بل يتم تعديله بواسطة هندسة المجال المغناطيسي، والتي هي نفسها نتاج تغذية الأشعة الكونية الراجعة.
• توفر إطارًا عدديًا قويًا لمقارنة نماذج انتشار الأشعة الكونية المعقدة مع البيانات عالية الدقة من AMS-02 و Voyager.

القيود والعمل المستقبلي:

• الدقة: دقة الشبكة الحالية (62.5 بارسيك) أخشن من عمليات المحاكاة الحديثة (1–4 بارسيك)، مما قد يؤدي لتفويت الاضطرابات صغيرة النطاق.
• الفيزياء المبسطة: يتجاهل النموذج تدفق الأشعة الكونية، والبنية الحرارية متعددة الأطوار للوسط بين النجمي، وتغذية المستعر الأعظم الحرارية/الحركية التفصيلية.
• الهندسة: تستخدم المحاكاة بقعة محلية مع حدود دورية، مما قد يؤدي إلى تحيزات مقارنة بنموذج هالة مجرية كامل.

الخلاصة:
نجح المؤلفون في إثبات أن كود PIERNIK MHD، الممتد بخوارزمية CRESP، يمكنه نمذجة إنتاج وانتشار نوى الأشعة الكونية بشكل ذاتي الاتساق. تسلط النتائج الضوء على أن ملاحظات نقل الأشعة الكونية (مثل B/C) حساسة للاقتران الديناميكي بين الأشعة الكونية، والغاز، والمجالات المغناطيسية، مما يتحدى افتراضات النماذج الظواهرية الثابتة.