Black Hole Spin-down in Collapsars in 3D Neutrino Transport GRMHD Simulations

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

يُعتقد أن انفجارات أشعة غاما طويلة الأمد (LGRBs) تُغذى بواسطة "الانهيارات النجمية" (collapsars)—وهي نجوم ضخمة تنهار نوياتها لتتحول إلى ثقوب سوداء (BHs)، مكونةً أقراصاً تراكمية تطلق نفاثات نسبوية عبر آلية بلاندفورد-زناجيك (BZ). وتعتمد قدرة هذه النفاثات ($L_{BZ} \propto a^2 \Phi^2$) بشكل حاسم على دوران الثقب الأسود ($a$) والتدفق المغناطيسي ($\Phi$).

أشارت محاكاة الديناميكا المائية المغناطيسية النسبية العامة (GRMHD) السابقة للأقراص غير الإشعاعية (الأديباتية) إلى أن الأقراض ذات التراكم المغناطيسي الأقصى (MADs) تعمل على تقليل دوران الثقب الأسود بسرعة إلى قيم دوران توازن منخفضة جداً ($a_{eq} \approx 0.035 - 0.07$). وهذه الدورات المنخفضة من شأنها إنتاج نفاثات ضعيفة جداً بحيث لا تستطيع اختراق النجم السلف، مما يتناقض مع الملاحظات التي ترصد انفجارات أشعة غاما قوية. علاوة على ذلك، تشير ملاحظات الموجات الثقالية (GW) لاندماج الثقوب السوداء الثنائية إلى توزيعات دوران منخفضة، ولكن تطور دوران الثقوب السوداء المولدة من الانهيارات النجمية في بيئات واقعية ذات تبريد مجهول.

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

2. المنهجية

قام المؤلفون بإجراء مجموعة من محاكاة GRMHD ثلاثية الأبعاد باستخدام الكود المعزز بالمعالجات الرسومية $\nu$H-AMR. وتشمل الميزات المنهجية الرئيسية ما يلي:

• نقل النيوترينو: تم تنفيذ مخطط M1 ثنائي العزم ومتكامل الطاقة ("الرمادي" - gray) لثلاثة أنواع من النيوترينوهات ($\nu_e, \bar{\nu}_e, \nu_x$). وهذا يسمح بنمذجة ذاتية متسقة لتبريد النيوترينو، مما يقلل من سمك القرص ويغير الخصائص الثرموديناميكية.
• نماذج النجوم السلف: تمت محاكاة ملفين مختلفين لكثافة النجوم لتغطية نطاق من معدلات التراكم الكتلي ($\dot{m}$):
  • نماذج p#: ملفات كثافة ذات قانون القوة ($\rho \propto r^{-1.5}$)، وهي نموذجية للانهيارات النجمية، وتنتج $\dot{m} \sim 0.1 - 1 M_\odot/s$.
  • نماذج c#: نوى ذات كثافة ثابتة تليها ملفات أكثر حدة، وتنتج معدلات تراكم عالية في البداية ($\dot{m} \gtrsim 10 M_\odot/s$).
• الشروط الأولية: بدأت المحاكاة بثقوب سوداء ذات دورانات أولية متفاوتة ($a = 0.0$ إلى $0.8$) وتكوينات مجال مغناطيسي مهيأة لتشكيل حالة الـ MAD.
• إطار التحليل: قام المؤلفون بحساب معامل الدوران ($s$) وفصل التدفقات إلى مكونات هيدروديناميكية (HD)، وكهرومغناطيسية (EM)، ونيوترينوية (NU) لاشتقاق نموذج شبه تحليلي لتطور الدوران. كما قاموا باستقراء هذه النتائج على مدار العمر النموذجي لمحرك انفجار غاما ($t_{eng} \sim 40$ ثانية).

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

• أول محاكاة لـ MAD مبردة بالنيوترينو: تعد هذه الدراسة الأولى التي تجمع بين GRMHD ثلاثية الأبعاد كاملة مع نقل النيوترينو لاستقصاء تطور دوران الثقب الأسود في الانهيارات النجمية.
• مراجعة دوران التوازن: تثبت الدراسة أن تبريد النيوترينو يغير دوران التوازن بشكل كبير، حيث وجدت قيمة $a_{eq} \approx 0.13$. وهذه القيمة أعلى بمرتين إلى أربع مرات من القيم التي وجدت في محاكاة MAD غير الإشعاعية ($a_{eq} \approx 0.035 - 0.07$).
• تحديد آلية تقليل الدوران كمياً: يوضح المؤلفون كمياً أنه بينما تمارس تدفقات النيوترينو نفسها عزماً مهملًا على الثقب الأسود، فإن تبريد النيوترينو أمر حاسم لأنه يعمل على ترقيق القرص. الأقراص الأرق تقلل من كفاءة عزم تقليل الدوران الذي يمارسه المجال المغناطيسي، مما يسمح للثقب الأسود بالاحتفاظ بدوران أعلى.
• تحليل استقرار النفاثات: يربط البحث بين دوران الثقب الأسود المنخفض وعدم استقرار النفاثات (عدم استقرار الالتواء - kink instability)، موضحاً لماذا تنتج النماذج ذات الدوران المنخفض نفاثات متغيرة أو فاشلة، خاصة في النجوم السلف ذات ملفات الكثافة الضحلة.

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

• دوران التوازن ($a_{eq}$):
  • في حالة الـ MAD، يتطور دوران الثقب الأسود بسرعة إلى حالة توازن حيث يكون معامل الدوران $s(a_{eq}) = 0$.
  • دوران التوازن المستخلص هو $a_{eq} \approx 0.13$.
  • هذه القيمة تتوافق مع توزيعات الدوران المستنتجة من رصد LIGO/Virgo/KAGRA لاندماج الثقوب السوداء الثنائية (التي تتركز حول $a \lesssim 0.3$).
• تعزيز قدرة النفاثات:
  • نظرًا لأن قدرة النفاث تتناسب مع $a^2$، فإن دوران التوازن الأعلى ($0.13$مقابل$0.07$) يؤدي إلى نفاثات أقوى بمقدار 4 إلى 16 مرة مقارنة بنماذج MAD غير الإشعاعية.
  • هذه النفاثات المعززة قادرة على اختراق النجم السلف بنجاح، مما يحل مشكلة "النفاث الفاشل" في نماذج الدوران المنخفض السابقة.
• مقاييس زمن تطور الدوران:
  • بالنسبة لمدد انفجارات غاما النموذجية ($t \gtrsim 30$ ثانية)، تمتلك الثقوب السوداء في كل من النماذج ذات معدلات التراكم العالية والمنخفضة وقتاً كافياً لتراكم كتلة كافية للوصول إلى $a_{eq} \approx 0.13$.
  • الأحداث الأقصر أو ذات التراكم الأقل قد تترك خلفها ثقوباً سوداء بدوران أعلى قليلاً ($a > 0.13$).
• دور النيوترينو:
  • عزوم النيوترينو المباشرة على الثقب الأسود مهملة.
  • الدور الأساسي للنيوترينوهات هو التبريد، الذي يقلل من نسبة سمك القرص إلى نصف قطره ($h/r$). الأقراص الأرق تكون أقل فعالية في نقل الزخم الزاوي بعيداً عن الثقب الأسود عبر العزوم المغناطيسية، مما يؤدي إلى دوران توازن أعلى.
• التقلب وعدم الاستقرار:
  • النماذج ذات الدورات الأولية المنخفضة ($a=0.1$) وملفات الكثافة الضحلة (p#) تظهر تقلباً كبيراً في التدفق المغناطيسي وكفاءة النفاث.
  • ينتج هذا عن عدم استقرار الالتواء العالمي في النفاثات الأضعف، مما يؤدي إلى عدم محاذاة القرص والخروج المؤقت من حالة الـ MAD. قد تفسر هذه الآلية انفجارات غاما خافتة أو فاشلة.

5. الأهمية

• حل مفارقة "الدوران المنخفض": توفق هذه الدراسة بين الحاجة إلى قدرة نفاث عالية (تتطلب دوراناً أعلى) وبين ملاحظات الموجات الثقالية للدوران المنخفض للثقوب السوداء الثنائية. فهي تقترح أن الثقوب السوداء الناتجة عن الانهيارات النجمية تستقر طبيعياً عند $a \approx 0.13$، وهي قيمة منخفضة بما يكفي لتتوافق مع بيانات الموجات الثقالية ولكنها عالية بما يكفي لتوليد انفجارات أشعة غاما ناجحة.
• الآثار على أسلاف الموجات الثقالية: إذا كانت انفجارات غاما هي أسلاف اندماج الثقوب السوداء الثنائية، فإن توزيع الدوران المتوقع ($a \sim 0.1 - 0.2$) يتوافق تماماً مع دراسات تعداد الموجات الثقالية الحالية، مما يدعم نموذج الانهيار النجمي لجزء من اندماجات الثقوب السوداء الثنائية.
• الحدود العليا لللمعان في انفجارات غاما: يوفر التراجع السريع في دوران الثقوب السوداء ذات الدوران العالي الأولي ($a \sim 1$) إلى $a \approx 0.13$ خلال الثواني الأولى من نشاط المحرك آلية فيزيائية للحد الأقصى المرصود في لمعان انفجارات غاما. حتى لو انهار نجم بقلب ذي دوران أقصى، فإن قدرة النفاث تُقيد بسبب التراجع السريع في الدوران قبل اختراق النفاث للنجم.
• النمذجة المستقبلية: يسلط هذا العمل الضوء على ضرورة تضمين نقل النيوترينو في محاكاة الانهيارات النجمية للتنبؤ بدقة بطاقة النفاثات وظروف تخليق العناصر في عملية r-process.

في الختال، تثبت هذه الورقة أن نماذج MAD المبردة بالنيوترينو تسمح للثقوب السوداء في الانهيارات النجمية بالحفاظ على دوران توازن أعلى ($a_{eq} \approx 0.13$) مما كان يُعتقد سابقاً، مما يمكن من إنتاج انفجارات أشعة غاما قوية وناجحة مع البقاء متسقة مع ملاحظات الموجات الثقالية للثقب السوداء الثنائية.