Evolution of proto-neutron stars to pulsars, magnetars and central compact objects

ملخص تقني: تطور النجوم النيوترونية الأولية إلى نباضات، ومغناطارات، وأجسام مدمجة مركزية

بيان المشكلة
لا يزال أصل المجالات المغناطيسية فائقة القوة المرصودة في المغناطارات ($B \sim 10^{14}–10^{15}$ G) موضوعاً للنقاش. وبينما يمكن لفرضية "المجال الأحفوري" (حفظ التدفق من النجوم السلف) أن تفسر مجالات النجوم النيوترونية التقليدية والأجسام المدمجة المركزية (CCOs)، إلا أنها تواجه صعوبة في تفسير مجالات المغناطارات دون استدعاء ظروف سلف نادرة ومتطرفة. وعلى العكس من ذلك، فإن فرضية الدينامو، التي تتطلب دورانًا سريعًا وحملًا حراريًا قويًا خلال مرحلة النجم النيوتروني الأولي (PNS)، مقبولة على نطاق واسع لتفسير المغناطارات، لكنها تواجه تحديًا: وهو أن معدل الولادة المستنتج للمغناطارات يقارب معدل المستعرات العظمى لانهيار القلب، مما يشير إلى أن تكوين المغناطارات ليس نادرًا كما توحي المتطلبات المعلمية المتطرفة. علاوة على ذلك، تشير حسابات حفظ التدفق القياسية إلى أنه حتى مع انهيار قلب حديدي، فإن المجالات الناتية ($B \sim 10^{10}$ G) غير كافية لتفسير مجالات النباضات العادية ($B \sim 10^{12}$ G)، مما يستلزم وجود آلية دينامو حتى بالنسبة للنباضات العادية. المشكلة المركزية هي التوفيق بين معدلات الولادة العالية لمختلف أنواع النجوم النيوترونية باستخدام إطار تطوري موحد.

المنهجية
يستخدم المؤلفون نموذج دينامو $\alpha-\Omega$ مدفوعًا بالقص ومختزلًا أحادي البعد، مقتبسًا من Wickramasinghe وآخرون (2014) الذين طبقوه في الأصل على اندماج الأقزام البيضاء. يتتبع النموذج تطور المجالات المغناطيسية الطولية ($B_p$) والالتفافية ($B_\phi$) خلال مرحلة النجم النيوتروني الأولي (PNS). وتشمل السمات الرئيسية للمنهجية ما يلي:

• المعادلات الحاكمة: يحل النموذج معادلات تفاضلية مقترنة لـ $B_p$ و $B_\phi$. يتم توليد المجال الالتفافي عبر تأثير $\Omega$ (الدوران التفاضلي المؤثر على المجال الطولي)، بينما يتم تجديد المجال الطولي عبر تأثير $\alpha$ (الحمل الحراري المضطرب المؤثر على المجال الالتفافي).
• تطور النجم النيوتروني الأولي (PNS): على عكس النماذج الساكنة، يدمج هذا الإطار الانكماش الديناميكي للـ PNS من نصف قطر أولي $R_0 \sim 40$ كم إلى نصف قطر نهائي $R_\infty \approx 12$ كم. ويتم تحديد تطور نصف القطر بناءً على Camelio وآخرون (2017) لمختلف الكتل ومعادلات الحالة (EoS).
• الكتلة والزخم الزاوي: يأخذ النموذج في الاعتبار نقصان الكتلة الثقالية بسبب التبريد النيوتريني وحفظ الزخم الزاوي، والذي يؤدي إلى تسريع دوران الـ PNS.
• اللزوجة وإخماد القص: يتم نمذجة إخماد الدوران التفاضلي (معدل القص $q$) عبر مزيج من لزوجة عدم استقرار المغناطوروتاسيوني (MRI) واللزوجة الحملية. يتلاشى معدل القص بمرور الوقت، مما يؤدي في النهاست إلى إنهاء الدينامو.
• فقدان التدفق والاستقرار: يتضمن النموذج حد إعادة الاتصال (reconnection) للمحاسبة على فقدان التدفق المغناطيسي (المحدد مبدئيًا بنسبة 20%)، كما يتضمن معايير الاستقرار (Braithwaite 2009) حيث تحدد عدم استقرارات المجال أزمنة اضمحلال $B_p$ و $B_\phi$.
• مسح فضاء المعلمات: يقوم المؤلفون بدمج المعادلات عددياً باستخدام طريقة Runge-Kutta، مع تغيير المعلمات الأولية بشكل منهجي بما في ذلك معدل القص ($q_0$)، والكتلة الأولية ($M_0$)، ونصف القطر الأولي ($R_0$)، والمجالات البذرية الأولية، والفترة الدورانية النهائية ($P_\infty$).

النتائج الرئيسية
تكشف المحاكاة أن التكوين النهائي للمجال المغناطيسي حساس للغاية لمعدل الدوران والقص الأولي، مما يسمح لآلية دينامو واحدة بإنتاج الطيف الكامل لأنواع النجوم النيوترونية المرصودة:

1. عتبة معدل القص: تعتمد قيمة التشبع للمجال الالتفافي ($B_{\phi, \infty}$) بقوة على معدل القص الأولي $q_0$. توجد عتبة حرجة عند $q_0 \simeq 0.23$.

   • أقل من العتبة ($q_0 < 0.23$): لا يتم تحفيز تأثير $\alpha$ بشكل كافٍ. يظل المجال الطولي قريبًا من القيمة التي يحددها حفظ التدفق ($B_p \approx 2.5 \times 10^{10}$ G)، ولا ينمو المجال الالتفافي بشكل ملحوظ. هذه المنطقة تقابل الأجسام المدمجة المركزية (CCOs).
   • فوق العتبة: يصبح الدينامو نشطًا. تؤدي معدلات القص الأعلى إلى مجالات التفافية أقوى.

2. تكوين المجموعات المختلفة:

   • المغناطارات: بالنسبة للفترات الأولية الأقصر (الدوران السريع، مثل $P_\infty \sim 2.5$ ms)، ينتج النموذج قوة تشبه المغناطارات بـ $B_p \simeq 10^{15}$ G و $B_\phi \simeq 10^{16}$ G.
   • النباضات العادية: بالنسبة للدوارات الأولية الأبطأ، ينتج النموذج مجالات نباضات نموذجية ($B_p \simeq 10^{12}$ G، $B_\phi \simeq 10^{14}$ G).
   • المغناطارات منخفضة الب-النباضات عالية الب: تنتج معدلات الدوران المتوسطة مجالات في نطاق المغناطارات منخفضة B أو النباضات عالية المجال.

3. دور حفظ التدفق: يتم تحديد قوة المجال الطولي الأولي إلى حد كبير من خلال حفظ التدفق أثناء الانهيار. ومع ذلك، فإن عملية الدينامو (تحديدًا تأثير $\alpha$) مطلوبة لتضخيم هذه المجالات إلى مستويات النباضات والمغناطارات المرصودة. يهيمن تأثير $\Omega$ بشكل عام على المجال الالتفافي.

4. الأجسام المدمجة المركزية (CCOs): يشير النموذج إلى أن الـ CCOs لا تمر بعملية دينامو قوية بسبب دورانها الأولي البطيء ($P_0 \sim 100$ ms). يتم تحديد مجالاتها الطولية بواسطة حفظ التدفق، ولكنها قد تكتسب أيضًا مجالات التفافية مضخمة فقط بواسطة تأثير $\Omega$، مما يفسر مجالات ثنائية القطب الصغيرة ولكنها تمتلك مكونات التفافية داخلية كبيرة محتملة.

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

• كلا من النباضات العادية والمغناطارات يمران على الأرجح بمرحلة دينامو؛ والاختلاف يكمن في كفاءة هذه العملية المدفوعة بالدوران الأولي.
• فرضية "المجال الأحفوري" غير كافية لتفسير مجالات النباضات العادية، مما يستلزم وجود دينامو حتى للنجوم النيوترونية غير المغناطارية.
• تمثل الأجسام المدمجة المركزية (CCOs) الطرف ذو الدوران المنخفض من هذا الطيف حيث يتم إخماد الدينامو، تاركةً المجالات تهيمن عليها عملية حفظ التدفق.

يؤكد المؤلفون أن نموذجهم أحادي البعد، رغم كونه مبسطًا، يلتقط التبعيات الجوهرية للدوران والقص والطفو، ويتوافق نوعيًا مع المحاكاة الأكثر تعقيدًا ثلاثية الأبعاد (مثل Raynaud وآخرون 2020). ويخلصون إلى أن دينامو $\alpha-\Omega$ يوفر إطارًا موحدًا لفهم تنوع المجالات المغناطيسية للنجوم النيوترونية الشابة دون الحاجة إلى ظروف سلف نادرة ومتطرفة للمغناطارات.