Mass distribution of neutron stars in binary systems

ملخص تقني: توزيع كتلة النجوم النيوترونية في الأنظمة الثنائية

بيان المشكلة
تشير الملاحظات إلى أن توزيع كتلة النجوم النيوترونية (NSs) المعروفة في الأنظمة الثنائية يظهر نمطاً ثنائي المنوال، مع قمم حول 1.4 كتلة شمسية و $\sim$1.8 كتلة شمسية. وبينما يُعزى ذروة الـ 1.4 كتلة شمسية غالباً إلى انهيار النواة القياسي في المستعرات العظمى أو تطور الغلاف المشترك، فإن أصل ذروة الـ $\sim$1.8 كتلة شمسية لا يزال محل نقاش. لقد كانت النماذج النظرية السابقة التي حاولت تفسير نمو كتلة النجم النيوتروني عبر التنامي محدودة بالحد الإدینگتون (Eddington limit)، أو لم تأخذ في الحسبان بشكل كامل أنظمة التنامي فائقة الإدینگتون. علاوة على ذلك، أعاقت الشكوك المتعلقة بكفاءة انتقال الكتلة، وتأثير المروحة (propeller effect)، وتطور المجال المغناطيسي، التنبؤات الدقيقة للتوزيع النهائي لكتلة النجوم النيوترونية المتنامية. تهدف هذه الدراسة إلى استقصاء ما إذا كان نموذج التنامي فائق الإدینگتون يمكن أن ينجح في إعادة إنتاج التوزيع ثنائي المنوال الملحوظ لكتلة النجوم النيوترونية في الأنظمة الثنائية.

المنهجية
استخدم المؤلفون نهجاً هجيناً لتخليق الجمهرة يجمع بين كودين متميزين:

1. تطور النجوم الثنائية (BSE): استُخدم لمحاكاة تشكل وتطور أولي لـ $10^7$ من الأنظمة الثنائية الأولية. أنتج هذا الكود جمهرة من النجوم النيوترونية الثنائية ذات رفاق من التسلسل الرئيسي، مع دمج دالات الكتلة الأولية، ونسب الكتلة، واللامركزية، وسرعات الركلة الوليدية. تم ترشيح الأنظمة لتحديد تلك القادرة على انتقال كتلة مستقر، مع استبعاد الأنظمة المعرضة لتطور غلاف مشترك غير مستقر أو معدلات انتقال كتلة قصوى ($>10^4 \dot{M}_{Edd}$).
2. نماذج التجارب في الفيزياء الفلكية النجمية (MESA): خضعت النجوم النيوترونية الثنائية المختارة من مخرجات BSE لحسابات تطورية تفصيلية باستخدام MESA (الإصدار r24.08.1). ركزت هذه المرحلة على عملية التنامي تحت ظروف تدفق "روتش" (Roche lobe overflow) المستقر.

النماذج الفيزيائية والافتراضات الرئيسية:

• التنامي فائق الإدینگتون: على عكس النماذج السابقة المحدودة بحد الإدینگتون، تدمج هذه الدراسة فيزياء التنامي فائق الإدینگتون كما وصفها ليبونوف (1987). يميز النموذج بين الأنظمة دون الحرجة وفائقة الحرجة بناءً على العلاقة بين نصف قطر ألففين ($R_a$) ونصف قطر التكور ($R_s$).
• كفاءة التنامي: يتم تحديد كسر الكتلة المنتقلة التي يتنامى عليها النجم النيوتروني ($f_{mt,a}$) بناءً على حالة التنامي. في النظام فائق الحرجة، يتم نفخ المادة الزائدة بواسطة الضغط الإشعاعي بالقرب من النجم النيوتروني، مما يؤدي إلى حمل زخم زاوي مداري نوعي بعيداً. يُفترض وجود كفاءة تنامي قدرها 0.5 للكتلة المتبقية.
• تطور المجال المغناطيسي: يتلاشى المجال المغناطيسي للنجم النيوتروني كدالة للكتلة المتنامية ($\Delta M$) وفق قانون أسي، حيث ينخفض من $5 \times 10^{12}$ غاوس إلى حد أدنى قدره $5 \times 10^8$ غاوس.
• تطور الدوران: يتتبع النموذج فترة دوران النجم النيوتروني عبر ثلاث مراحل: الراديوي (لا يوجد تنامي)، المروحة (طرد المادة)، والتنامي (زيادة السرعة الدورانية). يعتمد الانتقال بين هذه المراحل على المواقع النسبية لأسطوانة الضوء، ونصف قطر التدوير، ونصف قطر المغناطيسية.
• فقدان الزخم الزاوي: تأخذ المحاكاة في الاعتبار إشعاع الموجات الثقالية، والكبح المغناطيسي، وفقدان الزخم الزاوي الناتج عن المادة المنفوخة وتكون القرص المحيط بالثنائي.

النتائج الرئيسية

1. آلية نمو الكتلة: تُظهر الدراسة أن نمو كتلة النجم النيوتروني يعتمد بشدة على الفترة المدارية للنظام الثنائي وكتلة النجم المانح. في نموذج فائق الإدิงتون، يمكن للنجوم النيوترونية أن تظل في حالة توازن تنامي لفترات ممتدة، مما يسمح بتراكم كتلة كبير (نمو يصل إلى $\sim$1.0 $M_\odot$ من كتلة ولادة قدرها 1.4 $M_\odot$). وهذا يتناقض مع نموذج لي وآخرون (2021) (المحدود بالإدิงتون)، حيث تدخل النجوم النيوترونية غالباً في مرحلة "مروحة" طويلة الأمد تحد من نمو الكتلة.
2. إعادة إنتاج ثنائية المنوال: نجحت نتائج تخليق الجمهرة في إعادة إنتاج التوزيع ثنائي المنوال الملحوظ:
   • ذروة الـ $\sim$1.8 $M_\odot$: تنشأ هذه الذروة أساساً من الأنظمة الثنائية التي تكون فيها كتل النجوم المانحة $M_d < \sim$1.6 $M_\odot$ والفترات المدارية $P_{orb} < 20$ يوماً. في هذه الأنظمة، تسمح معدلات انتقال الكتلة بتنامٍ فعال وشبه متوازن (دون الإدิงتون) على فترات زمنية طويلة.
   • ذروة الـ 1.4 $M_\odot$: تُعزى هذه الذروة إلى أنظمة النجوم النيوترونية التي تمر بتطور الغلاف المشترك (غالباً ما يشمل مانحين أكثر كتلة أو تكوينات مدارية محددة)، حيث يكون نمو الكتلة في حده الأدنى.
   • الذيل عالي الكتلة: الأنظمة ذات المانحين متوسطي الكتلة ($1.6 \le M_d \le 3.0$ $M_\odot$) والفترات المدارية القصيرة يمكنها التنامي كمية كافية من المادة لتتجاوز 2.1 $M_\odot$.
3. توزيع فترة الدوران: يتنبأ النموذج بفترات دوران للنجوم النيوترونية تتركز بين 1 مللي ثانية و15 مللي ثانية للثنائيات ذات الكتلة المنخفضة في الأشعة السينية (LMXBs). النتائج أكبر قليلاً من بعض البيانات الرصدية، وهو ما يعزوه المؤلفون إلى افتراض حد أدنى للمجال المغناطيسي قدره $5 \times 10^8$ غاوس، مما يمنع قرص التنامي من الامتداد بالقرب الشديد من سطح النجم النيوتروني.

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

القصور والشكوك
يقر المؤلفون بوجود عدة شكوك في نموذجهم:

• كفاءة التنامي: الافتراض بأن الكفاءة هي 0.5 ومعالجة فقدان الكتلة بالقرب من النجم النيوتروني هي معالجات مبسطة.
• تلاشي المجال المغناطيسي: معدل التلاشي وقوة المجال الدنيا هي بارامترات تؤثر بشكل كبير على تطور الدوران وكفاءة التنامي.
• معادلة الحالة (EoS): تم تثبيت عزم القصور الذاتي عند $10^{45}$ جم سم$^2$، وهو ما يقرب من الحد الأدنى. قد تؤدي الاختلافات في معادلة الحالة للنجم النيوتروني إلى تغيير عزم القصور الذاتي، مما يؤثر على أزمنة زيادة السرعة الدورانية والانتقال إلى نظام المروحة.
• القنوات المستبعدة: لا تحدد الدراسة مساهمة النجوم النيوترونية المتكونة عبر الانهيار الناجم عن التنامي (AIC) أو تلك التي ولدت بتوزيعات كتلة ثنائية المنوال، حيث تركز بدلاً من ذلك على النجوم النيوترونية المتكونة عبر الانهيار النووي القياسي والتطور الثنائي اللاحق.