A Three-Dimensional Exploration of Magnetic Fields, Rotation, and Shock Revival in a $39 M_\odot$ Core-Collapse Supernova Progenitor

ملخص تقني: استكشاف ثلاثي الأبعاد للمجالات المغناطيسية، والدوران، وإحياء الصدمة في نجم سلف لمستعر أعظم انهياري لقلب بنجم تبلغ كتلته 39 كتلة شمسية

المشكلة والدافع
لا يزال إحياء الصدمة المتوقفة في المستعرات الأعظم ذات الانهيار اللبي (CCSNe) مشكلة معقدة تتضمن التفاعل بين التسخين النيوتريني، والدوران، والمجالات المغناطيسية، والديناميكا المائية. وبينما تعتبر آلية النيوترينو قوية في العديد من النجوم السلف، لا سيما تلك ذات الكثافة العالية، فإن الآلية المغناطيسية الدورانية (MR) تقدم مساراً بديلاً لإحياء الصدمة وتغذية الظواهر الطاقية (مثل الهايبرنوفا، والمستعرات الأعظم من النوع Ic-BL). تتطلب هذه الآلية دوراناً سريعاً ومجالات مغناطيسية قوية لتحويل الطاقة الدورانية الحرة إلى طاقة مغناطيسية عبر الالتفاف وعدم الاستقرار، مما يدفع تدفقات ثنائية القطب. ومع ذلك، فإن الدراسات السابقة غالباً ما اعتمدت على التماثل المحوري (ثنائي الأبعاد)، وهو ما قد يؤدي اصطناعياً إلى تثبيت النفاثات وتثبيط عدم الاستقرار غير المحوري مثل عدم استقرار الالتواء (kink instability). علاوة على ذلك، فإن التأثير المحدد للمجالات المغناطيسية على أزمنة إحياء الصدمة، ومورفولوجيا التدفق، واستقرار النجم النيوتروني الأولي ضد تشكل الثقب الأسود في البيئات ثلاثية الأبعاد عالية الكثافة وسريعة الدوران يتطلب مزيداً من الاستقصاء.

المنهجية
أجرى المؤلفون محاكاة هيدروديناميكية (HD) ومغناطيسية هيدروديناميكية (MHD) ثلاثية الأبعاد لنجم سلف بكتلة 39 كتلة شمسية عند مرحلة السكون (ZAMS) مع كثافة عالية ($\xi_{2.5} = 0.75$) ودوران سريع ($\Omega_c \approx 0.58$ راد ث$^{-1}$). تم تطوير النجم باستخدام MESA مع تطور متجانس شبه كيميائي. تمت محاكاة ثلاثة نماذج متميزة باستخدام كود FLASH مع حلّال Spark، مع استخدام تعزيز الشبكة التكيفي (AMR) ومخطط نقل نيوترينو متعدد الأبعاد ومعتمد على الطاقة مع إغلاق M1:

1. B1R1: نموذج دوار ومغناطيسي. تم تهيئة مجال قطبي بذرة ($B_0 \approx 5 \times 10^{11}$ غاوس)، محاذٍ لمحور الدوران.
2. B0R1: نموذج تحكم دوار وغير مغناطيسي (HD فقط).
3. B0R0: نموذج مرجعي غير دوار وغير مغناطيسي (HD فقط).

تطورت عمليات المحاكاة من مرحلة ما قبل الانهيار مروراً بارتداد القلب وصولاً إلى مرحلة ما بعد الارتداد المبكرة (حتى حوالي 1 ثانية)، باستخدام شبكة كارتيزية بدقة قصوى تصل إلى $\sim 0.61$ كم بالقرب من اللب. تم تقريب تأثيرات النسبية العامة باستخدام جهد فعال (الحالة A من Marek et al. 2006).

النتائج الرئيسية
حققت جميع النماذج الثلاثة توسعاً انفجارياً للموجة الصدمية ضمن الفترة المحاكات، لكنها أظهرت مورفولوجيا وأزمنة مختلفة تماماً:

• إحياء الصدمة والمورفولوجيا:

  • B1R1 (مغناطيسي/دوار): أظهر أسرع إحياء للصدمة ($t_{rev} \approx 150$ مللي ثانية بعد الارتداد) وطور أوضح تدفق ثنائي القطب. كان التدفق في البداية يشبه النفاث، لكنه أصبح عريضاً ومشوهاً بفعل عدم الاستقرار غير المحوري.
  • B0R0 (غير دوار): أحيا الصدمة لاحقاً ($t_{rev} \approx 250$ مللي ثانية) عبر اضطراب مدفوع بالنيوترينو. توسعت الصدمة في كتل كبيرة غير منتظمة دون اتجاه قطبي مفضل. والأهم من ذلك، انهار هذا النموذج إلى ثقب أسود بعد حوالي ثانية واحدة من الارتداد بسبب غياب الدعم الدوراني.
  • B0R0 (دوار/غير مغناطيسي): شهد أبطأ إحياء للصدمة ($t_{rev} \approx 370$ مللي ثانية). أدى الدوران إلى تثبيط الحمل الحراري وتقليل لومينانس النيوترينو، مما أدى إلى تطور أكثر اندماجاً مقارنة بـ B0R0. ولم يطور أي تدفق ثنائي القطب واضح.

• نقل الزخم الزاوي والطاقة:

  • في B1R1، أعادت إجهادات ماكسويل توزيع الزخم الزاوي بكفاءة نحو الخارج من النجم النيوتروني الأولي الداخلي، مما قلل من دوران اللب الداخلي واستخلص الطاقة الدورانية الحرة. هذه العملية استغلت خزان الدوران التفاضلي لتغذية التدفق القطبي.
  • كشفت ميزانية الزخم الزاوي أنه بينما وفرت عملية الحمل الزخم الزاوي إلى اللب الداخلي، وفرت عزوم الدوران المغناطيسية المصرف الرئيسي، مما عادل تسارع الدوران. عمل انبعاث النيوترينو كمصرف ثانوي أصغر.
  • نمت الطاقة التشخيصية ($E_{diag}$) في B1R1 في وقت مبكر وبسرعة أكبر، وهو ما يتوافق مع استخلاص الطاقة الدورانية. أشارت ميزانية الطاقة المغناطيسية إلى أن جزءاً كبيراً من الطاقة المستخلصة مر عبر حد تبديد فعال ($\dot{Q}_{diss}$)، المرتبط على الأرجح بإعادة الاتصال غير المحلولة وانتقال الشبكة الفرعية، بدلاً من الاحتفاظ بها كطاقة مغناطيسية محلولة.

• عدم الاستقرار والديناميكا ثلاثية الأبعاد:

  • تعرض التدفق المغناطيسي في B1R1 لعدم استقرار الالتواء غير المحوري (النمط $m=1$). أظهرت التشخيصات نمواً أسياً مبكراً لإزاحة مركز ثقل الطاقة المغناطيسية (زمن التضاعف $\sim 1.5$ مللي ثانية)، وهو ما يتوافق مع أزمنة الالتواء النظرية.
  • على عكس المحاكاة المحورية التي قد تتنبأ بنفاث مستقر وضيق، أظهرت المحاكاة ثلاثية الأبعاد ميل قناة التدفق (بعدة درجات) واتساعها بمرور الوقت. نجا التدفق وانتشر إلى أنصاف أقطار كبيرة، لكنه فقد توازنه في التركيز الضيق، مما يشير إلى أن التأثيرات ثلاثية الأبعاد حاسمة لنمذجة استقرار النفاث بدقة.

• التخليق النووي:

  • أظهرت المعالجة اللاحقة للمواد المقذوفة اختلافات تركيبية متميزة. أنتج B0R0 مواد غنية بالبروتونات ($Y_e > 0.5$) نموذجية للانفجارات المدفوعة بالنيوترينو، حيث صنع عناصر تصل إلى مجموعة الحديد ($A \sim 56$).
  • أنتج B1R1 مكوناً هاماً غنياً بالنيوترونات ($Y_e$ تصل إلى $\sim 0.33$) بسبب التطور السريع للتدفق، مما سمح بعملية r-ضعيفة تصل إلى $A \sim 90$. ومع ذلك، منعت البذرة المغناطيسية الأولية وتوسع الالتواء إنتاج المواد شديدة الغنى بالنيوترونات ($Y_e < 0.32$) التي لوحظت في النماذج ذات المجالات المغناطيسية الأكثر تطرفاً.

الأهمية والادعاءات
يزعم البحث أنه في حالة نجم سلف متطرف ذي كثافة عالية، يمكن للدوران السريع المقترن بالمجالات المغناطيسية القوية فتح مسار مبكر لإحياء الصدمة مدعوم مغناطيسياً، وهو مسار ثلاثي الأبعاد بطبيعته. وتتمثل الاستنتاجات الرئيسية في:

1. المجالات المغناطيسية كمحرك: الإجهادات المغناطيسية لا تشكل التدفق فحسب؛ بل تستخلص بنشاط الطاقة الدورانية الحرة، وتؤدي إلى إبطاء دوران اللب الداخلي، وتسهل إحياء الصمة في وقت أبكر من الحالات الدوارة غير المغناطيسية.
2. ضرورة البعد الثالث: التطور الطبيعي للتدفق المدفوع مغناطيسياً هو غير محوري. يعمل عدم استقرار الالتواء على توسيع وتغيير ميل التدفق، مما يمنع تشكل نفاث ضيق ومستمر. قد تبالغ المحاكاة المحورية في تقدير كل من استقرار وتمركز هذه التدفقات.
3. استقرار البقايا: يوفر الدوران دعماً كبيراً ضد التكون الفوري للثقب الأسود في النجوم السلف ذات الكثافة العالية، رغم أن الاستقرار طويل الأمد للبقايا يظل غير مؤكد خارج الفترة المحاكات.
4. مسار الطاقة: تتضمن عملية انتقال الطاقة من الدوران التفاضلي إلى التدفق تبديداً كبيراً غير محلول، مما يسلط الضوء على أهمية فيزياء الشبكة الفرعية في محاكات MHD لـ CCSN.

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