Impact of the equation of state on core collapse supernovae I: the low-$T/|W|$ instability

ملخص تقني: تأثير معادلة الحالة على انهيار النواة المستعر الأعظم I: عدم الاستقرار منخفض $T/|W|$

بيان المشكلة
تعد المستعرات العظمى لانهيار النواة (CCSNe) سريعة الدوران مصادر واعدة للانبعاثات متعددة الرسائل، حيث يمكن للديناميكيات غير المتماثلة محورياً في النجم النيوتروني الأولي (PNS) المتشكل حديثاً أن تترك بصمات مميزة على موجات الجاذبية (GWs) والنيوترينوهات. ومن أوجه عدم اليقين الحرجة في نمذجة هذه الأحداث هو معادلة الحالة (EOS) النووية للمادة الكثيفة، والتي تحكم انكماش النجم النيوتروني الأولي، وتطوره الحراري، وكثافته القصوى. وبينما استكشفت الدراسات السابقة تأثير معادلة الحالة على ديناميكيات الانفجار وانبعاث موجات الجاذبية، فإن التأثير المحدد لمعادلات الحالة النووية المختلفة على تطور عدم الاستقرار منخفض $T/|W|$ (LTWI) — وهو عدم استقرار مدفوع بالقص في النجوم النيوترونية الأولية سريعة الدوران — وعلى توقيعاته المرتبطة متعددة الرسائل، لم يتم توصيفه بشكل منهجي بعد. ويتميز الـ LTWI عن عدم الاستقرار النمطي لنمط القضيب (bar-mode instability) بكونه يمكن أن يتطور عند نسب أقل بكثير من الطاقة الدورانية إلى طاقة الجاذبية ($T/|W| \sim 0.04-0.06$)، وهي ظروف يمكن تحقيقها في حالات انهيار النواة سريعة الدوران.

المنهجية
أجرى المؤلفون سلسلة من محاكاة الديناميكا المائية المغناطيسية (MHD) ثلاثية الأبعاد (3D) باستخدام كود Aenus-ALCAR. ركزت الدراسة على انهيار نجم "وولف-رايت" (Wolf–Rayet) منفرد وسريع الدوران (35OC، $M_{\text{ZAMS}} = 35 M_\odot$) لعزل تأثيرات معادلة الحالة.

• الأسلاف والمجال المغناطيسي: استخدمت المحاكاة ملف الدوران الأصلي لنموذج 35OC ولكن مع استبدال المجال المغناطيسي لما قبل الانفجار بمجال ثنائي القطب محدد تحليلياً وضعيفاً ($B \sim 10^6$ G) لضمان طوبولوجيا مغناطيسية منتظمة مع تجنب كبح الـ LTWI الملاحظ في النماذج السابقة ذات المجالات الأقوى.
• معادلات الحالة: تم استخدام خمس معادلات حالة نووية متميزة ذات درجة حرارة نهائية: LS220، وSLy4، وBHB$\Lambda\phi$، وBBSk3، وTSO. تغطي هذه النماذج نطاقاً من الصلابة، ومعاملات طاقة التماثل، وتتضمن معالجات مختلفة للهيبيرونات والتفاعلات النووية (نهج Skyrme مقابل نهج المجال المتوسط النسبي).
• الإعداد العددي: تم تطوير المحاكاة في تناظر محوري حتى 5 مللي ثانية قبل ارتداد النواة، ثم نُقلت إلى شبكات ثلاثية الأبعاد. امتدت الشبكة إلى $4 \times 10^{10}$ سم مع دقة عالية في المناطق الداخلية ($5 \times 10^4$ سم). تمت معالجة انتقال النيوترينو عبر مخطط طيف ذي لحظتين.
• أدوات التحليل: تم استخراج إشارات موجات الجاذبية باستخدام صيغة الرباعي الأقطاب القياسية مع تفكيك القشرة. تم توصيف عدم الاستقرار باستخدام تفكيك التوافقات الكروية للكثافة، والتحليل الفوري (Fourier)، وتفكيك النمط التجريبي التجميعي (EEMD) لعزل دوال الأنماط الجوهرية (IMFs). كما تم تتبع "صلابة" النجم النيوتروني الأولي باستخدام عدد لوف (Love number) المدمر ($\kappa_2$) المحسوب من المقاطع متوسطة التناظر الكروي.

المساهمات والنتائج الرئيسية

1. متانة الـ LTWI: تطور الـ LTWI في جميع نماذج معادلة الحالة الخمس المدروسة، مما يؤكد أن حدوثه أمر متين لهذا السلف سريع الدوران. يظهر عدم الاستقرار في شكل أنماط حلزونية واسعة النطاق ($m=1, 2$) تولد انبعاثات موجات جاذبية شبه دورية وتعدل لومات النيوترينو.
2. الاختلافات المعتمدة على معادلة الحالة: بينما يعد ظهور عدم الاستقرار أمراً متيناً، فإن خصائصه المحددة تختلف بشكل كبير مع معادلة الحالة:
   • البداية والمورفولوجيا: تراوحت أوقات البداية بين $\sim 50$ مللي ثانية إلى $\sim 100$ مللي ثانية بعد الارتداد. وتنوع الهيكل السمتي المهيمن؛ فعلى سبيل المثال، أظهر نموذج SL هيكلاً حلزونياً مستمراً بذراعين، بينما أظهرت النماذج الأخرى انتقالات بين تكوينات بذراع واحدة أو اثنتين أو ثلاث.
   • الارتباط الترددي: وُجد ارتباط خطي واضح بين تردد موجات الجاذبية المهيمن المرتبط بالـ LTWI والصلابة الفعالة (مقلوب عدد لوف المدي $\kappa_2$) للنجم النيوتروني الأولي. أنتجت معادلات الحالة الأكثر صلابة (أقل $\kappa_2$) ترددات موجات جاذبية أعلى (تصل إلى $\sim 370$ هرتز)، بينما أنتجت المعادلات الأكثر ليونة ترددات أقل.
   • تعديل النيوترينو: أحدث الـ LTWI تعديلات شبه دورية مميزة في لومات النيوترينو، لا سيما في الاتجاهات القريبة من المستوى الاستوائي. ظهرت هذه التعديلات عند ترددات ($\sim 150-180$ هرتز) متميزة عن التوصيل المفاجئ وSASI.
3. تطور ما قبل الارتداد: أثر اختيار معادلة الحالة على الطبقية والهيكل الدوراني للنواة الداخلية قبل الارتداد، مما أدى إلى اختلافات في نصف قطر النواة الداخلية، وكتلتها، وملفات الدوران التفاضلي عند لحظة ارتداد النواة، حتى مع وجود ملفات دوران أولية متطابقة.
4. تفكيك الإشارة: أثبتت الدراسة أن تقنية EEMD تعزل بفعالية مكون إشارة موجات الجاذبية المرتبط بالـ LTWI (تحديداً IMF 6)، وتفصله عن إشارات الارتداد والأنماط التذبذبية الأخرى.
5. القابلية للكشف:
   • موجات الجاذبية: إشارات موجات الجاذبية قابلة للكشف بواسطة الكواشف الحالية (LIGO، Virgo، KAGRA) للأحداث المجرية، وبالكواشف من الجيل التالي (Einstein Telescope، Cosmic Explorer) للأحداث في كامل المجموعة المحلية. يعمل التناظر الاستوائي للنظام على كبح الاستقطاب $\times$ للمراقبين الاستوائيين، مما يقلل مسافات الكشف بمعامل قدره اثنين مقارنة بالمشاهدات القطبية.
   • النيوترينوهات: بينما تختلف طاقة النيوترينو الإجمالية بشكل طفيف فقط ($\sim 10\%$) عبر معادلات الحالة، فإن معادلة الحالة تؤثر على درجة حرارة الغلاف النيوتروني وسعة ومدة التعديلات المدفوعة بـ LTWI. ومن المتوقع أن تكون توقيعات التعديل قابلة للرصد في الكواشف ذات الإحصائيات العالية مثل IceCube للأحداث المجرية.

الأهمية والادعاءات
يدعي البحث أن تردد إشارة موجات الجاذبية الناتجة عن عدم استقرار $T/|W|$ المنخفض يعمل كأداة تشخيصية محتملة لمعادلة الحالة الكثيفة للمادة في حالات انهيار النواة سريعة الدوران. وهذا يقدم قيداً تكميلياً للمعلومات التي تحملها إشارات النيوترينو. وتؤكد الدراسة أنه بينما يعد الـ LTWI سمة متينة للنوى سريعة الدوران، فإن بصماته متعددة الرسائل المحددة (الترددات، وفترات الحياة، والمورفولوجيا) حساسة للفيزياء الدقيقة للمادة الكثيفة المشفرة في معادلة الحالة.

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