Formation of heavy double neutron stars II: the role of heavy first-born neutron stars and low metallicity
تُثبت هذه الورقة أن النجوم النيوترونية المزدوجة الثقيلة، مثل GW190425، هي نواتج نادرة (حوالي 0.5% من الجمهرة) تتشكل أساساً عبر القنوات القياسية بدلاً من آليات متميزة، وهي تمثل مجرد الذيل عالي الكتلة لمجتمع النجوم النيوترونية المزدوجة القياسي دون اعتماد كبير على المعدنية.
في أعماق الهمهمة الهادئة للكون، يمكن لنجمين ميتين أن يلتفا حول بعضهما البعض في دوامة، محبوسين في عناق أخير وقاتل. هذه هي النجوم النيوترونية المزدوجة، وهي النوى المنهارة لنجوم ضخمة استنفدت وقودها وانفجرت. وعندما تصطدم، ترسل تموجات عبر الزمكان تسمى موجات الجاذبية، وهي ظاهرة أحدثت ثورة في كيفية استماعنا إلى الكون. لسنوات، قام علماء الفلك بفهرسة هذه الاصطدامات، ووجدوا أن معظمها يتضمن أزواجًا بوزن إجمالي يبلغ حوالي ضعفين ونصف ضعف كتلة شمسنا. لكن حدثًا واحدًا، تم رصده في عام 2019، كسر هذا النمط. هذا الاصطدام، الذي سمي GW190425، تضمن زوجًا من النجوم النيوترونية ثقيلًا لدرجة أن كتلتهما المجتمعة بلغت حوالي 3.4 ضعف كتلة الشمس. كان هذا الاكتشاف لغزًا لأن مثل هذا الزوج الثقيل لم يُرَ قط في مجرتنا، درب التبانة، حيث تتبع علماء الفلك بعناقة هذه الأنظمة عشرات المرات باستخدام التلسكوبات الراديوية. إن وجود هذا الثنائي الثقيل يتحدى فهمنا لكيفية ولادة وتطور هذه الأزواج النجمية، مما يجبر العلماء على التساؤل عما إذا كانت هناك آلية خفية تخلقها، أو ما إذا كانت نظرياتنا الحالية تفتقد قطعة حاسمة من القصة.
ولحل هذا الغموض، لجأ فريق من الباحثين إلى أقوى الأدوات المتاحة لمحاكاة حياة النجوم: نماذج حاسوبية مفصلة تتبع فيزياء التطور النجمي من البداية إلى النهاية. لقد ركزوا على اللحظة المحددة التي يبدأ فيها نظام ثنائي، يحتوي بالفعل على نجم نيوتروني واحد، رقصته الأخيرة. في هذه الأنظمة، يبدأ نجم ثانٍ، مجرد من طبقاته الخارجية ليكشف عن قلب هليوم ساخن وكثيف، في نقل المادة إلى رفيقه النجم النيوتروني. أراد الباحثون معرفة ما إذا كانت هذه العملية يمكن أن تنتج بشكل طبيعي الأزواج الثقيلة التي تظهر في بيانات موجات الجاذبية، أم أنها تتطلب سيناريوهات غريبة وغير مرجحة. لقد بنوا شبكة واسعة من المحاكاة، واختبروا كيف تتغير النتيجة عندما يكون النجم النيوتروني الأول أثقل أو أخف، وعندما يتشكل النظام في بيئات ذات تركيبات كيميائية مختلفة، تُعرف بالمعادنية. ومن خلال تشغيل هذه النماذج، استطاعوا مراقبة مصير آلاف الأزواج المحتملة ورؤية أي منها سيصمد ليصبح عمليات الاندماج الثقيلة التي ترصدها المراصد.
كشفت نتائج هذه المحاكاة المكثفة أن تكوين مثل هذه الأزواج الثقيلة هو حدث نادر، يحدث في حوالي نصف بالمائة فقط من جميع أنظمة النجوم النيوترونية المزدوجة. وخلافًا لبعض النظريات التي اقترحت أن هذه الأنظمة الثقيلة يجب أن تأتي من قناة اندماج سريعة خاصة أو تتطلب نوعًا محددًا جدًا من النجوم الضخمة، وجد الباحثون أنها تنشأ من نفس العملية القياسية التي تخلق الأزواج الأخف وزنًا التي نراها في مجرتنا. في هذا المسار القياسي، يتمدد نجم الهليوم ويصب المادة على النجم النيوتروني في تدفق مستقر، لينهار في النهاية إلى نجم نيوتروني ثانٍ. أظهرت عمليات المحاكاة أنه حتى مع أكثر النجوم النيوترونية المولودة أولاً ضخامة أو أكثر نجوم الهليوم ضخامة، فإن قناة "الاندماج السريع"، حيث يحدث انهيار ثانٍ غير مستقر بسرعة، لا تساهم في تكوين هذه الأزواج الثقيلة عند معدنية شمسنا. بدلاً من ذلك، فإن الأنظمة الثقيلة هي ببساطة ذيل الكتلة العالية للمجتمع الطبيعي، أي الاستثناءات الإحصائية لعملية شائعة.
كما بحثت الدراسة فيما إذا كانت البيئة الكيميائية لميلاد النجم، وتحديدًا كمية العناصر الثقيلة، تلعب دورًا. في الكون، تولد النجوم في مناطق تحتوي على عناصر ثقيلة أقل تفقد كتلة أقل بفعل الرياح النجمية، مما قد يترك خلفه نوى أثقل. اختبر الباحثون ذلك من خلال محاكاة أنظمة في بيئات تحتوي على عشر مرات أقل ومائة مرة أقل من العناصر الثقيلة مقارً بشمسنا. وبينما وجدوا أن انخفاض المعدنية يسمح بوجود نجوم نيوترونية ثانية أكثر ضخامة قليلاً، إلا أن الجزء الإجمالي للنجوم النيوترونية المزدوجة الثقيلة لم يتغير بشكل كبير. كان الفرق صغيرًا للغاية لدرجة أنه يشير إلى أن تكوين هذه الأنظمة الثقيلة لا يعتمد بقوة على التركيب الكيميائي للمجرة. وتتعارض هذه النتيجة مع بعض الأفكاء السابقة التي تقول إن الأزواج الثقيلة هي فئة فرعية متميزة تتشكل فقط في الكون المبكر فقير المعادن.
في النهاية، تخلص الورقة البحثية إلى أن النجم النيوتروني المزدوج الثقيل GW190425 لا يتطلب تفسيرًا جديدًا أو غريبًا. فمن المرجح أنه مجرد تباين نادر وطبيعي لعملية التكوين القياسية. إن حقيقة أننا لم نرَ مثل هذه الأزواج الثقيلة في مجرتنا هي على الأرجح بسبب الندرة الشديدة لهذا الحدث وصعوبة رصدهم باستخدام التلسكوبات الراديوية، وليس بسبب اختلاف جوهري في كيفية صنعهم. ويقترح الباحثون أنه مع زيادة حساسية كواشف موجات الجاذبية، قد نجد أن هذه الأزواج الثقيلة أكثر شيوعًا في الكون البعيد منها في جوارنا المحلي، ببساطة لأن رصدها أصعب عندما تكون شابة وسريعة الدوران. في الوقت الحالي، يتم حل لغز GW190425 ليس من خلال العثات على قانون جديد للفيزياء، بل من خلال الاعتراف بأن الطبيعة تنتج أحيانًا بطلاً ثقيلاً من نفس الوصفة التي تصنع بقية الفريق.
بيان المشكلة يمثل اكتشاف مصدر موجات الجاذبية GW190425، وهو نظام نجمي نيتروني مزدوج (DNS) بكتلة إجمالية تبلغ 3.4−0.1+0.3M⊙، تحدياً كبيراً لنظريات التكوين الحالية. لم يتم رصد مثل هذه الأنظمة الثقيلة من النجوم النترونية المزدوجة (Mtot≥3M⊙) في قرص مجرة درب التبانة، حيث يهيمن على الجمهرة الأنظمة التي تقع في نطاق 2.3–2.9M⊙. وبينما تم اقتراح عدة سيناريوهات تكوين لتفسير GW190425 — تتراوح بين قنوات "الاندماج السريع" (fast-merger) التي تتضمن انتقال كتلة غير مستقر، ومسارات تتطلب نجوماً نترونية أولى (NS1) ضخمة أو أسلاف نجوم هيليوم (He) ضخمة — إلا أن هذه الفرضيات غالباً ما تم اختبارها باستخدام افتراضات فيزيائية مستقلة. هناك نقص في وجود إطار عمل متسق ذاتياً يعمل في آن واحد على معايرة الجمهرة المرصودة للنجوم النترونية المزدوجة في المجرة واختبار مدى جدوى هذه القنوات البديلة لتكوين نجوم نترونية مزدوجة ثقيلة عبر مختلف مستويات المعدنية.
المنهجية توسع هذه الدراسة العمل الذي قام به Nair & Stevenson (2025, الورقة الأولى) من خلال توظيف إطار تطور ثنائي مفصل باستخدام كود التطور النجمي أحادي البعد MESA (الإصدار 23.05.01). يقوم المؤلفون بنمذجة تطور ما بعد الغلاف المشترك (post-common envelope) لثنائيات نجم الهيليوم + النجم النتروني (He star + NS)، مع اعتبار نظام (He star + NS) كنقطة البداية.
تشمل الميزات المنهجية الرئيسية ما يلي:
توسيع فضاء المعلمات: على عكس الورقة الأولى، التي استخدمت كتلة ثابتة للنجم النتروني الأول (NS1) قدرها 1.4M⊙، يستكشف هذا العمل نطاقاً واسعاً من كتل النجوم النترونية المولودة أولاً (1.1–1.9M⊙) وكتل نجوم الهيليوم الأولية (2.5–10M⊙).
تغيرات المعدنية: تُجرى المحاكاة عند معدنية شمسية (Z=Z⊙=0.02) وفي بيئتين منخفضتي المعدنية (0.1Z⊙ و 0.ط0.01Z⊙) لتقييم تأثير الرياح النجمية على كتل البقايا.
انتقال الكتلة والرياح: تفترض النماذج تراكم كتلة محدود بحد إدنغتون (M˙Edd=3×10−8M⊙yr−1) على النجم النتروني. وتتم نمذجة فقدان كتلة الرياح لنجوم الهيليوم باستخدام صيغة Nugis & Lamers (2000) عند المعدنية الشمسية، وصيغة معدلة من Belczynski et al. (2010a) مع اعتماد Vink & de Koter (2005) على المعدنية (β=0.86) للنماذج منخفضة المعدنية.
تخليق الجمهرة: لتقي lòng المقارنة مع الأرصاد، يطبق المؤلفون مخططات ترجيح بناءً على التوزيعات المرصودة لكتل نجوم الهيليوم الأولية (قانون القوة)، والفترات المدارية (مسطحة في اللوغاريتم)، وكتل NS1 (نموذج التوزيع الطبيعي المزدوج للنباضات المعاد تدويرها ونموذج قانون القوة للبدء).
ديناميكيات ما بعد المستعر الأعظم (Post-SN): يتم تحديد تكوين النجم النتروني الثاني (NS2) وبقاء الثنائي باستخدام صيغة معدلة من Mandel & Müller (2020) لكتلة البقايا والركلات الولادية (natal kicks)، والمعايرة لإعادة إنتاج جمهرة النجوم النترونية المزدوجة في المجرة.
المساهمات الرئيسية والنتائج تختبر الورقة بشكل منهجي أربعة قنوات تكوين مقترحة للنجوم النترونية المزدوجة الثقيلة (Mtot≥3M⊙):
التكوين القياسي: تجاوز روتش المستقر من النوع (Case BB Roche-lobe overflow) مع MNS1<1.7M⊙.
مسار نجم الهيليوم الضخم:MHe≥8.8M⊙ مع تباين في كتل NS1، وتجنب المزيد من انتقال الكتلة.
مسار NS1 الضخم:MNS1≥1.7M⊙ يخضع لعملية انتقال كتلة مستقرة من النوع (Case BB/BC RLO).
الثنائيات منخفضة المعدنية: التكوين عند Z<Z⊙.
تشير النتائج إلى:
ندرة النجوم النترونية المزدوجة الثقيلة: تعتبر النجوم النترونية المزدوجة الثقيلة نادرة، حيث تشكل حوالي 0.5% فقط من إجمالي جمهرة النجوم النترونية المزدوجة عند المعدلية الشمسية.
هيمنة القناة القياسية: "قناة التكوين القياسي" هي المسار المهيمن للنجوم النترونية المزدوجة الثقيلة، حيث تساهم بنسبة ∼0.5% من الجمهرة. أما مساري "نجم الهيليوم الضخم" و"NS1 الضخم" فيساهمان بنسب ضئيلة للغاية (∼6×10−5% و ∼2×10−4% على التوالي) لأن نجوم الهيليوم الضخمة ونجوم NS1 الضخمة نادرة جوهرياً في التوزيعات المعتمدة.
عدم وجود قناة الاندماج السريع عند المعدلية الشمسية: وجدت الدراسة عدم وجود أي مساهمة من قناة "الاندماج السريع" (انتقال كتلة غير مستقر يؤدي إلى مرحلة غلاف مشترك ثانية) عند المعدلية الشمسية.
الاعتماد على المعدنية: يظهر كسر النجوم النترونية المزدوجة الثقيلة اعتماداً ضعيفاً جداً على المعدنية. وبينما تنتج النماذج منخفضة المعدنية (0.1Z⊙) نجوم NS2 أكثر كتلة قليلاً بسبب انخفاض فقدان الكتلة عبر الرياح، إلا أن الكسر الإجمالي للنجوم النترونية المزدوجة الثقيلة يظل مشابهاً للمعدلية الشمسية (∼4.8% مقابل ∼5%). ومع ذلك، فإن الآلية تختلف: عند المعدنية المنخفضة، تتكون النجوم النترونية المزدوجة الثقيلة بشكل أساسي من خلال انتقال كتلة غير مستقر، بينما تتكون عند المعدلية الشمسية عبر انتقال كتلة مستقر.
اتساق الجمهرة: تتوافق جمهرة النجوم النترونية المزدوجة المحاكات، عند ترجيحها بفضاء المعلمات الموسع، بشكل عام مع توزيع الكتلة الكلية والارتباط بين الفترة المدارية واللامركزية (eccentricity) المرصود في المجرة، رغم أنها تنتج بشكل طفيف فائضاً في الأنظمة ذات Mtot≈2.55M⊙.
الأهمية والاستنتاجات يخلص المؤلفون إلى أن النجوم النترونية المزدوجة الثقيلة لا تشكل جمهرة فرعية منفصلة تتطلب فيزياء تكوين غريبة. بدلاً من ذلك، هي تمثل الذيل عالي الكتلة لجمهرة النجوم النترونية المزدوجة القياسية المتكونة عبر التطور الثنائي المعزول. تشير الدراسة إلى أن قناة التكوين القياسية (Tauris et al. 2017) كافية لتفسير وجود أنظمة تشبه GW190425، بشرما إذا تضمن فضاء المعلمات النطاق الكامل لكتل NS1 والمعدنيات.
تشير الورقة إلى أنه بينما تعيد نماذجنا إنتاج جمهرة النجوم النترونية المزدوجة المرصودة في المجرة، فإنها لا تفسر تماماً معدل الاندماج العالي لـ GW190425 مقارلة بقلة وجود مثل هذه الأنظمة في المجال المجري. يقترح المؤلفون أنه إذا استمرت ملاحظات الموجات الثقالية المستقبلية في الكشف عن نسبة كبيرة من النجوم النترونية المزدوجة الثقيلة، فقد يعني ذلك أن جزءاً كبيراً من النجوم النترونية المزدوجة يولد ثقيلاً ولكنه يظل غير مكتشف في المسوحات الراديوية بسبب تأثيرات الاختيار (مثل نقص إعادة التدوير، أو الأعمار القصيرة، أو تأثيرات التسارع). ويؤكد هذا العمل على أهمية الأطر المتسقة ذاتياً للتمييز بين قنوات التكوين المختلفة، ويسلط الضوء على أن المعدنية تلعب دوراً ثانوياً في كسر النجوم النترونية المزدوجة الثقيلة، وإن كانت تؤثر على استقرار عملية انتقال الكتلة المؤدية لتكوينها.