Multi-messengers from the radioactive decay of $r$-process nuclei

ملخص تقني: الرسائل المتعددة من التحلل الإشعاعي لنوى عملية الـ r (الالتقاط السريع للنيوترونات)

بيان المشكلة
تعد عملية الالتقاط السريع للنيوترونات (r-process) مسؤولة عن تخليق العناصر الثقيلة، وهي تحدث في بيئات فيزيائية فلكية قاسية مثل اندماج النجوم النيوترونية، ومن المحتمل أيضاً في أكينات (cocoons) انفجارات أشعة غاما. وبينما تُغذى ظواهر "الكيلونوفا" (kilonova) الناتزلة عن طريق التحلل الإشعاعي لنوى عملية الـ r، فإن النمذجة الحالية للمنحنيات الضوئية للكيلونوفا تعتمد على فرضيات مبسطة فيما يتعلق بترسيب الطاقة من نواتج التحلل. وتحديداً، غالباً ما تحسب النماذج الحالية معدلات توليد الطاقة عبر قيم Q للتحلل وتفترض توزيعات حرارية بسيطة للجسيمات المنبعثة (الإلكترونات، والنيوترينوهات، وأشعة غاما، والنيوترونات). ويتجاهل هذا النهج السمات الطيفية غير الحرارية التفصيلية لمنتجات التحلل، والتي تعد بالغة الأهمية لفهم التوازن الحراري وتطور المنحنى الضوئي. علاوة على ذلك، تظل المساهمات النسبية لمواقع فيزيائية فلكية مختلفة في نواتج عملية الـ r غير مؤكدة بسبب محدودية البيانات التجريبية للنظائر الغنية بالنيوترونات.

المنهجية
يقدم المؤلفون حساباً من المبادئ الأولى لأطياف الانبعاث الناتجة عن تحلل $\beta$ لنوى عملية الـ r، من خلال ربط بيانات البنية النووية التفصيلية بشبكة تفاعلات نووية.

• الحساب الطيفي: تستخدم الدراسة إطار عمل إحصائياً متعدد المراحل (بناءً على M. Mumpower et al. 2025b) لحساب أطياف الانبعاث للإلكترونات، والنيوترينوهات المضادة، وأشعة غاما، والنيوترونات المتأخرة عن تحلل $\beta$. وبالنسبة لدالة قوة $\beta$، يستخدم المؤلفون البيانات التجريبية حيثما توفرت، وإلا فإنه يتم استخدام تقريب "الارتباط العشوائي شبه الجسيمي" (QRPA). ويتم اشتقاق أطياف الإلكترونات والنيوترينوهات بافتراض أن إجمالي طاقة التحلل تتوزع بينهما، بينما تُحسب انبعاثات أشعة غاما والنيوترونات باستخدام نهج "هاوزر-فايشباخ" الإحصائي لنمذجة إزالة الإثارة ونسب التفرع.
• شبكة التفاعلات: تُربط هذه الأطياف الفردية بشبكة التفاعلات PRISM (الروتينات المحمولة لنمذجة التخليق النووي المتكامل) (الإصدار 1.6.0). وتحاكي الشبكة التخليق النووي باستخدام نموذج "قطرة المدى المحدود" (Finite Range Droplet Model) للمدخلات النووية، وكود CoH3 لمعدلات الالتقاط الإشعاعي والانشطار.
• المسار الفيزيائي الفلكي: يتتبع المحاكاة تطور مسار محدد (المسار b من M. R. Mumpower et al. 2025a) الذي يمثل أكين (cocoon) انفجار أشعة غاما. وقد اختير هذا المسار لقدرته على إعادة إنتاج بقايا عملية الـ r الشمسية، لا سيما في منطقة الأكتينيدات. تتبع الظروف الفيزيائية تمدداً أديباتياً (كظومياً) بكثافة ابتدائية قدرها $3.2 \times 10^4$ جم/سم$^3$، ودرجة حرارة ابتدائية قدرها 2 جيجا كلفن، وكسر إلكتروني $Y_e = 0.034$.
• حساب التدفق: يتم تحليل إجمالي أطياف الانبعاث $S^{(i)}(E, t)$ إلى تدفق التفاعل $F(t)$ (التحللات في الثانية لكل وحدة كتلة) والطيف الانبعاثي لكل تحلل $S(E)$. وهذا يسمح بحساب إجمالي أعداد الجسيمات، وتدفقات الطاقة، ومتوسط الطاقات كدوال في الزمن.

النتائج الرئيسية
تكشف الدراسة أن توزيعات الانبعاث تعتمد بشدة على الزمن، وهي غير حرارية وتمتلك طاقات متوسطة كبيرة:

• الطبيعة غير الحرارية: الطاقات المتوسطة للجسيمات المنبعثة أعلى بكثير مما هو متوقع في حالة التوازن الحراري (والذي سيكون في نطاق الكيلو إلكترون فولت keV). خلال الثواني القليلة الأولى، تمتلك الإلكترونات والنيوترينوهات طاقات متوسطة تبلغ $\sim$5 ميجا إلكترون فولت، تنخفض إلى 1.5–2.5 ميجا إلكترون فولت عند 10 ثوانٍ، وتظل فوق 1 ميجا إلكترون فولت حتى يتوقف التقاط النيوترونات ($\sim$100 ثانية). في المقابل، تحافظ أشعة غاما والنيوترونات على طاقات متوسطة أقل ($\lesssim$1 ميجا إلكترون فولت) ولكنها تظل غير حرارية أيضاً.
• التطور الطيفي:
  • الإلكترونات والنيوترينوهات: أطيافها متطابقة تقريباً في الشكل، وتمتد إلى طاقات عالية ($\sim$10–20 ميجا إلكترون فولت) في البداية. ويستمر ذيل عالي الطاقة بشكل سلس لفترات زمنية طويلة، ويُعزى ذلك إلى نواتج الانشطار.
  • أشعة غاما: تكون ثنائية المنوال في البداية بقمة منخفضة الطاقة ($\sim$1 ميجا إلكترون فولت) وذيل عالي الطاقة، ثم يتطور الطيف إلى توزيع أحادي المنوال بحلول $\sim$1 ثانية. في الأوقات المبكرة، يكون الطيف سلساً بسبب النمذجة النظرية للنوى قصيرة العمر. وبحلول $\sim$1 يوم، تظهر خطوط طيفية محددة من النوى طويلة العمر (على سبيل المثال، الخط 2.6 ميجا إلكترون فولت من $^{208}$Tl والخطوط من $^{140}$La).
  • النيوترونات: يكون طيف النيوترونات كبيراً في البداية ولكنه ينخفض بسرعة بعد $\sim$100 ثانية مع انخفاض انبعاث النيوترونات المتأخر عن تحلل $\beta$ وانتهاء التقاط النيوترونات.
• توزيع الطاقة: توزيع طاقة التحلل ليس ثابتاً. تحمل النيوترينوهات الجزء الأكبر (40–50%) طوال فترة التطور. في الأوقات المبكرة ($\lesssim$1 ثانية)، تكون مساهمة أشعة غاما صغيرة، بينما تكون مساهمة الإلكترونات مهيمنة. أما مساهمة طاقة النيوترونات فتكون كبيرة فقط في الأوقات المبكرة جداً.

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

علاوة على ذلك، تقدر الورقة إمكانية الرصد المباشر لهذه الإشارات داخل مجرة درب التبانة:

• النيوترينوهات: بالنسبة لحدث مجري (كتلة مقذوفات $0.01 M_\odot$، مسافة 15 كيلو فرسخ فلكي)، يمكن أن ينتج عن أول 10 ثوانٍ من الانبعاث حوالي 2.2 حدثاً قابلاً للكشف في كاشف بحجم "هايبر كاميوكاندي" (Hyper-Kamiokande) عبر التحلل البيتا العكسي.
• أشعة غاما: رغم أنها تكون معتمة في البداية، يُتوقع أن تصبح أشعة غاما قابلة للرصد بعد أسابيع إلى أشهر من الحدث مع انخفاض عتامة المقذوفات. ويمكن رصد خطوط طيفية محددة (مثل تلك الناتجة عن $^{208}$Tl) بشكل مباشر، مما يوفر وسيلة لسبر أغوار التخليق النووي بشكل مستقل عن نمذجة المنحنى الضوئي للكيلونوفا.

تخلص الدراسة إلى أنه بينما يتطلب تحديد أوقات الرصد الدقيقة إجراء عمليات نقل إشعاعي كاملة، فإن الكشف المباشر عن النيوترينوهات وأشعة غاما من أحداث عملية الـ r المجرية يمثل مساراً واعداً لتجاوز غموض النمذجة وسبر عملية تكوين العناصر الثقيلة بشكل مباشر.