$r$-process Heating Feedback on Disk Outflows from Neutron Star Mergers

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

تُعد اندماجات النجوم النيوترونية (BNSMs) المواقع الفيزيائية الفلكية الأساسية لعملية التقاط النيوترونات السريعة (r-process)، والتي تُصنّع العناصر الثقيلة وتغذي ومضات الكيلونوفا العابرة. وبينما تنجح المحاكاة العددية الحالية لتدفقات الأقراص في مرحلة ما بعد الاندماج في نمذجة التسخين اللزوجي وإعادة اتحاد جسيمات ألفا، إلا أنها تتجاهل عادةً التغذية الراجعة للطاقة الناتجة عن تسخين عملية r-process.

• الفجوة: غالبًا ما تتعامل المحاكاة القياسية مع التسخين النووي كخطوة معالجة لاحقة أو تستخدم توصيفات مبسطة غير معتمدة على التركيب الكيميائي. ومع ذلك، فإن إطلاق الطاقة من تخليق نووي r-process حساس للغاية للتركيب النووي (تحديدًا الكسر الإلكتروني، $Y_e$) والتاريخ الحراري للمواد المقذوفة.
• النتيجة: يؤدي تجاهل هذه التغذية الراجعة إلى تنبؤات غير دقيقة بكتلة المواد المقذوفة، وسرعتها، وشكلها المورفولوجي، وهي أمور بالغة الأهمية لتفسير الملاحظات متعددة الرسائل (الموجات الثقالية ونظائرها الكهرومغناطيسية) ونمذجة منحنيات ضوء الكيلونوفا.

2. المنهجية

طوّر المؤلفون توصيفًا جديدًا ومتسقًا ذاتيًا لربط تسخين r-process مباشرة في محاكاة الهيدروديناميكا اللزوجية طويلة المدى باستخدام كود FLASH.

• إعداد المحاكاة:
  • الهندسة: هيدروديناميكا لزوجية ثنائية الأبعاد متماثلة محورياً.
  • الشروط الأولية: تم نمذجتها بناءً على حدث GW170817، وتتكون من ثقب أسود دوار ($2.65 M_\odot$، دوران $0.8$) محاط بقرص تراكمي ($0.1 M_\odot$).
  • الفيزياء: تشمل تسرب النيوترينو (مخطط من 3 أنواع)، ونموذج لزوجة قرص $\alpha$ (تم اختباره مع $\alpha = 0.03$ و $0.06$)، ومعادلة حالة هيلمهولتز.
• توصيف التسخين (الابتكار الجوهري):
  • جسيمات التتبع (Tracer Particles): أدخل المؤلفون حوالي $10^6$ من جسيمات التتبع ذات "الذاكرة" لتتبع تاريخ الكسر الإلكتروني ($Y_e$) للسائل.
  • تعريف $Y_{e, T6}$: عندما يمر جسيم تتبع لأول مرة عبر درجة حرارة تبلغ $6$جيجا كلفن، يتم تسجيل كسر الإلكتروني الخاص به كـ$Y_{e, T6}$. وتعمل هذه القيمة كوسيط للتركيب النووي الذي سيحرك عملية r-process.
  • جدول معدل التسخين: تم إنشاء جداول لمعدلات إطلاق الطاقة النووية ($\dot{\epsilon}_{nuc}$) مسبقة الحساب لـ 50 قيمة مختلفة من $Y_e$ الأولية (من $0.01$إلى$0.50$) بناءً على مسارات توسع معلمية (الإنتروبيا$s=30 k_B$/نوية، والزمن$\tau=10$ مللي ثانية).
  • طريقة السحابة في الخلية (CIC): لربط بيانات جسيمات التتبع اللاغرانجية بشبكة أويلر، استخدم المؤلفون طريقة CIC. وتسمح هذه الطريقة للجسيمات بترسيب معدلات التسخين على خلايا الشبكة بناءً على درجة الحرارة المحلية و $Y_{e, T6}$ المسجل لكل جسيم.
  • العتبات: يتم تطبيق التسخين فقط عندما تنخفض درجة حرارة السائل المحلية تحت عتبة معينة ($T_{heat} = 4$ جيجا كلفن أو $6$ جيجا كلفن)، مما يقرب بداية تخليق r-process. أما المواد المتدفقة للداخل ($v_r < 0$) فتتلقى تسخينًا سالبًا (تبريدًا) لمنع تراكم الطاقة الاصطناعي في الدورات الحملية.
• تنوع النماذج: تم تشغيل عشرة نماذج، مع تغيير اللزوجة ($\alpha$)، وعتبات التسخين ($T_{heat}$)، ودقة فئات $Y_e$، وتضمين تبريد التدفق للداخل. استُخدم نموذج مقارنة (al03K) يستخدم معدل تسخين يعتمد على الزمن ومنتظم مكانيًا (بناءً على Klion et al. 2022) لاختبار أهمية الاعتماد على التركيب الكيميائي.

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

1. الربط المتسق ذاتيًا: تطوير طريقة قوية لربط تسخين r-process ديناميكيًا بالهيدروديناميكا باستخدام جسيمات التتبع ذات الذاكرة وخوارزمية CIC، مما يغني عن الحاجة إلى شبكات تفاعل نووي كاملة أثناء المحاكاة.
2. الاعتماد على التركيب الكيميائي: إثبات أن معدل التسخين يجب أن يعتمد على التاريخ النووي المحدد ($Y_e$) لعنصر السائل، وليس فقط على الزمن أو درجة الحرارة.
3. التأثير على الديناميكا: تحديد كيفية تغيير تسخين r-process للتطور العالمي لتدفقات الأقراص، وتحديدًا فيما يتعلق بالكتلة غير المرتبطة، وتوزيع السرعات، والشكل المورفوري (المورفولوجي).

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

• زيادة كتلة المقذوفات: يؤدي تضمين تسخين r-process إلى زيادة إجمالي كتلة مقذوفات القرص غير المرتبطة بنسبة تقارب 10% مقارنة بالنماذج المرجعية (التي تتضمن فقط إعادة اتحاد جسيمات ألفا).
• تعزيز السرعة:
  • يتم تعزيز السرعة الشعاعية للمقذوفات الغنية بالنيوترونات ($Y_e < 0.25$) بعامل قدره اثنين (على سبيل المثال، من $\sim 0.03c$ إلى $\sim 0.06c$).
  • يحدث هذا نتيجة للتسخين المكثف والمبكر في المناطق ذات $Y_e$ المنخفضة (القرص الخارجي)، مما يؤدي إلى تسريع هذه المكونات بسرعة.
  • المكونات ذات $Y_e$ الأعلى تشهد تسخينًا متأخرًا وأضعف بسبب فترات الإقامة الأطول في القرص الداخلي الأكثر سخونة.
• التغيرات المورفولوجية:
  • يدفع تسخين r-process التدفق نحو شكل مورفولوجي أكثر كروية. حيث يعمل الزيادة العالمية في الضغط الحراري على تنعيم عدم التماثل واسع النطاق وكبح الحركات الحملية هامشية الارتباط بالقرب من الثقب الأسود.
  • النماذج التي تفتقر إلى التسخين تحتفظ بمزيد من المواد ذات السرعة المنخفضة والحركة الحملية.
• المقارنة مع التسخين المنتظم:
  • النموذج الذي استخدم معدل تسخين منتظم مكانيًا ومعتمدًا على الزمن (al03K) قلل بشكل كبير من تقدير التغذية الراجعة. فقد أنتج طاقة تسخين إجمالية أقل وفشل في إعادة إنتاج السرعات العالية ومعدلات إخراج الكتلة التي أظهرتها النماذج المعتمدة على التركيب الكيميائي.
  • وهذا يؤكد أن التباين المكاني والزماني للتسخين (المدفوع بالتركيب الكيميائي) أمر حاسم للديناميكا الدقيقة.
• تأثيرات اللزوجة: تؤدي اللزوجة الأعلى ($\alpha=0.06$) إلى تطور أسرع وكتل مقذوفات أعلى قليلاً، لكن التأثير النسبي لتسخين r-process يظل كبيرًا عبر مختلف مستويات اللزوجة.

5. الأهمية

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

في الختام، يوضح Ma et al. أن تسخين r-process ليس مجرد تصحيح طفيف، بل هو آلية فيزيائية مهيمنة تعمل على تسريع تدفقات الأقراذ بشكل كبير، وتزيد من كتلتها، وتُشكل هندستها الكروية، مما يغير بشكل جذري البصمات المرصودة لاندماجات النجوم النيوترونية.