Relativistic jets from millisecond proto-magnetars

ملخص تقني: النفاثات النسبية من النجوم المغناطيسية الأولية ذات الميلي ثانية

بيان المشكلة
لطالما اقتُرح أن النجوم المغناطيسية الأولية ذات الميلي ثانية (millisecond proto-magnetars) — وهي نجوم نيوترونية سريعة الدوران ($P \sim 1$ ms) وذات مجال مغناطيسي قوي ($B \gtrsim 10^{15}$ G) تتشكل في حالات المستعرات العظمى نتيجة انهيار القلب، أو اندماج النجوم النيوترونية، أو الانهيار المستحث بالتراكم في الأقزام البيضاء (AIC) — تمثل المحركات المركزية لانفجارات أشعة غاما (GRBs). ومع ذلك، لا يزال هناك تحدٍ نظري حاسم: خلال الثواني القليلة الأولى بعد الولادة، تبعث هذه الأجرام تدفقات نيوترينو مكثفة تعمل على تسخين السطح ودفع رياح غنية بالباريونات. تقوم هذه "الرياح النيوترينية" بتحميل التدفق بالمواد، مما قد يؤدي إلى خفض المغناطيسية ($\sigma$) والحد من عامل لورنتز الكتلي ($\Gamma_\infty$) الذي يمكن تحقيقه. تنبأت النماذج شبه التحليلية السابقة، التي تمت معايرتها بناءً على محاكاة أحادية البعد وافتراضت هندسات "القطب المنفرد المنقسم" (split-monopole) متوسطة الزاوية، بأن $\sigma \lesssim 1$ في الأوقات المبكرة، مما يشير إلى أن النجوم المغناطيسية الأولية لا يمكنها إطلاق النفاثات فائقة النسبية ($\Gamma \gtrsim 10-100$) المطلوبة لتشغيل انفجارات أشعة غاما قصيرة المدى. ولم يكن قد تم التعامل بشكل كامل مع مسألة ما إذا كانت التأثيرات متعددة الأبعاد يمكن أن تتغلب على تحميل الباريونات هذا لإنتاج نفاثات نسبية في غضون ثوانٍ من التكون.

المنهجية
يقدم المؤلفون أول محاكاة ثلاثية الأبعاد باستخدام الديناميكا المائية المغناطيسية النسبية العامة (GRMHD) لرياح النجوم المغناطيسية الأولية المسخنة بالنيوترينو. تستخدم الدراسة نسخة معدلة من كود GRHydro مفت المصدر (جزء من Einstein Toolkit) وتعتمد على مخطط حجم نهائي (finite-volume scheme) مع إعادة بناء WENO-Z ومحلل ريمان تقريبي من نوع HLLE. الميزات العددية الرئيسية تشمل:

• نقل النيوترينو: تم استخدام تقريب العزم الواحد ($M_0$) لمعادلة بولتزمان النسبية العامة لنمذجة تسخين وتبريد النيوترينو بشكل متسق ذاتياً، بما في ذلك عمليات $\beta$ ذات التيارات المشحونة، وتلاشي الأزواج، والتشتت.
• تطور المجال المغناطيسي: يتم تطوير المجال المغناطيسي عبر الجهد الكهرومغناطيسي الرباعي باستخدام مخطط نقل القيود (constraint transport) في مقياس لورنز العلوي للحفاظ على شرط "سولينويدال" بدقة الآلة.
• الشروط الأولية: تفترض النماذج نجماً نيوترونياً بكتلة $1.4 M_\odot$ ودرجة حرارة مركزية تبلغ $20$MeV. تمتد المحاكاة عبر نطاق من فترات الدوران الأولية ($P \approx 1 - 18$ ms) وقوة مجال مغناطيسي قطبي ($B_P \approx 10^{13} - 3 \times 10^{15}$ G).
• الإعداد: تستخدم المحاكاة تناظر الانعكاس عبر المستوى الاستوائي وتستخدم تكثيف الشبكة التكيفي (Berger-Oliger) حيث تقوم الشبكة الداخلية بحل سطح النجم المغناطيسي الأولي. تم تطوير النماذج أولاً للوصول إلى حالة مستقرة هيدروديناميكية مدفوعة بالنيوترينو قبل بدء إدخال مجال ثنائي قطب قطبي.

المساهمات والنتائج الرئيسية
تكشف المحاكاة أن التأثيرات متعددة الأبعاد تغير هيكل الرياح بشكل جوهري مقارنة بالنماذج أحادية البعد أو متوسطة الزاوية، لا سيما في النماذج سريعة الدوران ($P \approx 1$ ms) وعالية المغناطيسية ($B_P \approx 3 \times 10^{15}$ G):

1. التدرج الزاوي للتدفق: التدفق ليس موحداً بل ينقسم إلى ثلاثة مناطق متميزة:

   • الحزام الاستوائي ($\theta \gtrsim 30^\circ$): تعزز القوى الطاردة المركزية فقدان الكتلة بشكل كبير بالقرب من خط الاستواء الدوراني. وهذا يخلق رياحاً كثيفة دون نسبية السرعة ($v \sim 0.1c$) مع تحميل باريوني عالٍ.
   • النفاث القطبي ($\theta \lesssim 10^\circ$): تتعرض خطوط المجال القطبية لتعزيز ضئيل من الناحية الطاردة المركزية لتحميل الكتلة. وبناءً على ذلك، تظل هذه المنطقة "نظيفة" (تحميل باريوني منخفض) ومغناطيسية للغاية.
   • منطقة الانتقال: خطوط العرض المتوسطة حيث تتسارع الرياح بأكبر قدر من الكفاءة.

2. التركيز الذاتي للنفاث: تعمل الرياح الاستوائية الكثيفة دون النسبية كوسط حاصر، مما يعمل بشكل طبيعي على تركيز التدفق القطبي الأقل تحميلاً بالباريونات. يحدث هذا التركيز الذاتي عند أنصاف أقطار صغيرة ($\lesssim 100$ كم) دون الحاجة إلى الحصر بواسطة مقذوفات المستعر الأعظم أو مقذوفات الاندماج.

3. تحقيق الظروف فائقة النسبية: في النماذج الأكثر دورانًا ومغناطيسية (مثل النموذج B3P1)، يحقق النفاث القطبي قيم مغناطيسية ذروية تبلغ $\sigma \sim 30 - 100$. وهذا كافٍ للوصول إلى عوامل لورنتز كتلية $\Gamma_\infty \sim 100$ على مقاييس أكبر. يبلغ نصف زاوية فتح النفاث حوالي $8^\circ$، وهو ما يتوافق مع التقديرات الرصدية لانفجارات أشعة غاما القصيرة والطويلة.

4. توزيع الطاقة: يتم تقسيم طاقة الرياح الإجمالية بين نفاث قطبي فائق النسبية ورياح استوائية واسعة الزاوية دون نسبية. تحمل الأخيرة معظم الطاقة الحركية والكتلة، مما قد يغذي مقذوفات المستعر الأعظم أو الكيلونوفا، بينما يغذي الأول انفجار أشعة غاما (GRB).

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

• تعدد الأبعاد أمر ضروري: لا يمنع تحميل الباريونات تشكل النفاثات النسبية عند أخذ الاعتماد على خط العرض للقوى الطاردة المركزية وتسخين النيوترينو في الاعتبار.
• سيناريو المحرك الموحد: يمكن لمحرك مركزي واحد تشغيل كل من انفجار أشعة غاما قصير المدى (عبر النفاث القطبي) وانفجار المستعر الأعظم/الكيلونوفا المرتبط به (عبر الرياح الاستوائية)، مما يفسر التقسيم المرصود بين طاقات GRB المصححة حسب التوجيه وطاقات المقذوفات الحركية.
• الاتساق الزمني: تعمل هذه الآلية ضمن الثواني الأولى من حياة النجم المغناطيسي الأولي، وهو ما يتوافق مع مرحلة الانبعاث الفوري لانفجارات أشعة غاما القصيرة.

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