The Radius of the Neutron Star PSR J0614-3329 from NICER Data
تقدم هذه الورقة تحليلاً قائماً على مرصد NICER للنجم النيوتروني PSR J0614-3329، والذي يستنتج نصف قطر محيطي استوائي يتراوح بين 9.88 و12.77 كم باستخدام نموذج البقعة الساخنة الثلاثية، مع الإشارة إلى التحديات في الملاءمة المشتركة لبيانات NICER وXMM-Newton ومناقشة القيود الناتجة على خصائص المادة الكثيفة.
المؤلفون الأصليون:M. C. Miller, A. J. Dittmann, I. M. Holt, F. K. Lamb, C. Chirenti, Z. Arzoumanian, J. Berteaud, S. Bogdanov, K. C. Gendreau, A. K. Harding, W. C. G. Ho, C. Kalapotharakos, S. M. Morsink, P. S. Ray, R.M. C. Miller, A. J. Dittmann, I. M. Holt, F. K. Lamb, C. Chirenti, Z. Arzoumanian, J. Berteaud, S. Bogdanov, K. C. Gendreau, A. K. Harding, W. C. G. Ho, C. Kalapotharakos, S. M. Morsink, P. S. Ray, R. A. Remillard, Z. Wadiasingh, M. T. Wolff
المؤلفون الأصليون: M. C. Miller, A. J. Dittmann, I. M. Holt, F. K. Lamb, C. Chirenti, Z. Arzoumanian, J. Berteaud, S. Bogdanov, K. C. Gendreau, A. K. Harding, W. C. G. Ho, C. Kalapotharakos, S. M. Morsink, P. S. Ray, R. A. Remillard, Z. Wadiasingh, M. T. Wolff
في أعماق الكون، مختبئةً داخل بقايا النجوم المنفجرة، تكمن أجرام شديدة الكثافة لدرجة أن ملعقة صغيرة واحدة من مادتها ستزن بقدر جبل على الأرض. هذه هي النجوم النيوترونية، وهي النوى المنهارة للنجوم الضخمة التي نفد وقودها. في داخلها، تُسحق المادة إلى حالة لا يمكن إعادة إنتاجها في أي مختبر على الأرض؛ حيث تتفكك الذرات، وتُجبر البروتونات والإلكترونات على الاندماج معاً لتشكل بحراً من النيوترونات. لطالما أراد الفيزيائيون معرفة كيف تتصرف هذه المادة الغريبة بدقة تحت هذا الضغط الهائل. هل تظل سائلاً بسيطاً، أم تتحول إلى شيء أكثر غرابة، مثل بلورة فائقة الكثافة أو حساء من الجسيمات الغريبة؟ تكمن الإجابة في حجم هذه النجوم؛ فإذا تمكنا من قياس مدى كبر حجم نجم نيوتروني لوزن معين، يمكننا استنتاج القواعد التي تحكم المادة بداخله.
لسنوات، كافح علماء الفلك للحصول على قياس دقيق لأحجام هذه النجوم. يكمن التحدي في أن النجوم النيوترونية صغيرة، وبعيدة، وغالباً ما تكون محجوبة خلف سحب من الغبار أو مستترة بسبب جاذبيتها الشديدة. ومع ذلك، فإن جيلاً جديداً من التلسكوبات الفضائية قد غير قواعد اللعبة. إحدى هذه الأدوات، وهي "مستكشف التركيب الداخلي للنجم النيوتروني" (NICER)، تدور حول الأرض وتراقب هذه النجوم بدقة مذهلة. فهي لا تكتفي بالتقاط الصور فحسب، بل تعدّ جسيمات الأشعة السينية الفردية القادمة من سطح النجم، وتضبط توقيت وصولها بجزء من مليار من الثانية. ومن خلال مراقبة كيفية تغير الضوء المنبعث من البقع الساخنة على سطح النجم أثناء دورانه، يمكن للعلماء استخدام تشوه الزمكان الناتج عن جاذبية النجم لحساب حجمه.
في دراسة حديثة، وجه فريق من الباحثين اهتمامهم إلى نجم نيوتروني دوار محدد يُعرف باسم PSR J0614−3329. يقع هذا النجم على بعد حوالي 600 سنة ضوئية ويدور 318 مرة في كل ثانية. وقد حددت عمليات الرصد الراديوية السابقة بالفعل كتلته لتكون حوالي 1.4 إلى 1.5 ضعف كتلة شمسنا، مما يجعله مرشحاً مثالياً لهذا النوع من الدراسات. جمع الباحثون بيانات لسنوات من جهاز NICER، مع تصفية التداخل الناتج عن الغلاف الجوي للأرض وضوضاء الفضاء الخلفية لعزل الإشارة الخافتة القادمة من النجم. ثم بنوا نماذج حاسوبية مفصلة لمحاكاة شكل منحنى ضوء الأشعة السينية إذا كان للنجم أحجام مختلفة وأنماط مختلفة من البقع الساخدة على سطحه.
اختبر الفريق عدة سيناريوهات مختلفة؛ حيث نظروا في نماذج يكون فيها للنجم بقعتان ساخنتان دائريتان، وبقعتان بيضاويتان، وثلاث بقع دائرية. وبعد إجراء ملايين عمليات المحاكاة لمعرفة أي نموذج يطابق البيانات الفعلية بشكل أفضل، وجدوا أن النموذج الذي يحتوي على ثلاث بقع دائرية ساخنة هو الأكثر دقة. يشير هذا النموذج إلى أن نصف قطر النجم يتراوح بين 9.88 و12.77 كيلومتراً. ولوضع هذا في المنظور، إذا أخذت كرة من المادة بكتلة شمسنا وضغطتها حتى تصبح في حجم مدينة صغيرة فقط، فستصبح نجماً نيوترونياً. كما حاول الباحثون دمج بيانات NICER مع ملاحظات قديمة من تلسكوب آخر وهو XMM-Newton، لكنهم وجدوا أن الطريقة القياسية للتعامل مع الضوضاء الخلفية في البيانات القديمة جعلت النتائج المدمجة غير موثوقة، لذا قرروا الاكتفاء ببيانات NICER عالية الدقة وحدها لاستنتاجهم النهائي.
تعد نتائج هذه الدراسة مهمة لأنها تضيف نقطة بيانات دقيقة جديدة إلى الخريطة المتنامية لأحجام النجوم النيوترونية. فعندما دمج الباحثون نتائجهم الخاصة بـ PSR J0614−33−32 مع بيانات من نجوم نيوترونية أخرى، وجدوا أن المادة داخل هذه النجوم ليست صلبة كما توقعت بعض النظريات. تشير القياسات الجديدة إلى أن الوزن الأقصى الذي يمكن أن يتحمله النجم النيوتروني قبل الانهيار إلى ثقب أسود هو أقل قليلاً مما كان يُعتقد سابقاً، وأن النجم النيوتروني النموذجي الذي يمتلك كتلة شمسنا هو أصغر قليلاً من التقديرات السابقة. يساعد هذا في تضييق نطاق السلوكيات المحتملة للمادة عند كثافات تفوق بكثير ما يمكننا إنتاجه على الأرض. وبينما لا تكشف الدراسة عن التركيب الدقيق لنواة النجم، إلا أنها تستبعد العديد من الاحتمالات الأكثر تطرفاً وتقربنا من فهم القوانين الأساسية التي تحكم أكثر المواد كثافة في الكون. ويقف هذا العمل كشهادة على كيف يمكن لعدّ جزيئات الضوء الفردية من نجم دوار بعناية أن يكشف أسرار أكثر البيئات تطرفاً في الوجود.
ملخص تقني: نصف قطر النجم النيوتروني PSR J0614−3329 من بيانات NICER
المشكلة والدوافع إن التركيب الداخلي للنجوم النيوترونية، وتحديداً معادلة الحالة (EOS) للمادة الباردة والمتوازنة كيميائياً عند كثافات تتجاكل كثافة النواة (n0≈0.16 fm−3)، لا يمكن محاكاتها في المختبرات الأرضية أو اشتقاقها من المبادئ الأولية. وبناءً على ذلك، تُعد الملاحظات الفلكية التي تضع قيوداً على كتلة وأبعاد هذه الأجرام أمراً ضرورياً. وبينما قدم مرصد (NICER) سابقاً تقديرات لنصف قطر عدة نباضات (مثل PSR J0030+0451 و PSR J0740+6620)، تقدم هذه الورقة تحليلاً جديداً لبيانات NICER للنباض PSR J0614−3323 الذي يبلغ تردده 318 هرتز. وتعتبر هذه الكتلة هدفاً حاسماً لتنقيح قيود معادلة الحالة، حيث أن كتلتها تبلغ حوالي 1.4–1.5M⊙ وفقاً للتوقيت الراديوي.
المنهجية قام المؤلفون بتحليل بيانات أداة التوقيت بالأشعة السينية (XTI) التابعة لـ NICER والتي جُمعت بين مارس 2018 ومايو 2023، مما نتج عنه تعرض فعال قدره 1.09 ميجا ثانية. كما تضمن التحليل بيانات أرشيفية من مرصد XMM-Newton (جهاز EPIC) من عام 2010، رغم أن المطابقة المشتركة كشفت عن تباينات منهجية في نمذجة الخلفية لبيانات XMM-Newton.
يتبع نهج النمذجة أطر الاستدلال البايزي (Bayesian inference) المستخدمة في تحليلات NICER السابقة:
نموذج الانبعاث السطحي: ينمذج المؤلفون سطح النجم النيوتروني بوجود مناطق ساخنة ("بقع ساخنة") تتعرض للقصف بجسيمات نسبوية من الغلاف المغناطيسي. تُعامل هذه البقع كمناطق دائرية أو بيضاوية ذات درجة حرارة موحدة. ويتم نمذجة طيف الانبعاث المحلي باستخدام غلاف جوي للهيدروجين غير مغناطيسي ومتأين بالكامل (Ho & Lai 2001).
تتبع الأشعة (Ray Tracing): تم استخدام "تقريب شوابزشيلد المفلطح" (oblate Schwarzschild approximation) لأخذ التسطح الناتج عن دوران النجم في الاعتبار، مع تقريب الزمكان الخارجي كزمكان شوابزشيلد.
الاستدلال الإحصائي: يستخدم التحليل عينة pocoMC (و emcee للتدقيق) لاستكشاف توزيعات الاحتمال البعدي لمعلمات النموذج. تشمل المعلمات الرئيسية الكتلة الثقالية (M)، ونصف القطر المحيطي الاستوائي (Re)، وهندسة البقعة (الموقع، الحجم، درجة الحرارة)، وميل المراقب، والامتصاص بين النجمي (NH).
مقارنة النماذج: تم اختبار ثلاث فئات من النماذج: بقعتان دائريتان، بقعتان بيضاويتان، وثلاث بقع دائرية.
النتائج الرئيسية
النموذج المفضل: تدعم الأدلة البايزية بقوة النموذج الذي يحتوي على ثلاث بقع ساخنة دائرية ذات درجة حرارة موحدة (ΔlnZ>10 مقارنة بالتكوينات الأخرى).
قيود نصف القطر والكتلة: يعطي أفضل نموذج ملائم للبقع الثلاث نطاقاً موثوقاً متماثلاً بنسبة 68% لنصف القطر المحيطي الاستوائي يتراوح بين 9.88–12.77 كم. أما توزيع الكتلة البعدي فهو مهيمن عليه من قبل "الفرضية المسبقة" (prior) المستمدة من التوقيت الراديوي، ويظل متسقاً مع نطاق ∼1.4–1.5M⊙، بينما يتم تحديد توزيع نصف القطر بقوة من خلال بيانات الموجة السينية.
الهندسة: تتميز الهندسة المفضلة بمناطق انبعاث غير متقابلة، مما يشير إلى تكوين مجال مغناطيسي أكثر تعقيداً من مجرد ثنائي قطب مركزي بسيط، وقد يتضمن مكونات متعددة الأقطاب.
جودة المطابقة: يعطي أقصى حد للمطابقة (maximum-likelihood fit) قيمة χ2/dof=5154.89/5252 (احتمالية 83%) وقيمة χ2/dof بولومترية قدرها 26.42/27 (احتمالية 50%)، مما يشير إلى مطابقة كافية دون وجود بقايا أو قيم شاذة كبيرة.
تباين XMM-Newton: تؤدي المطابقات المشتركة لبيانات NICER و XMM-Newton إلى نقص منهجي في توقع عدد العدات (counts) لـ XMM-Newton عند تطبيق افتراضات الخلفية القياسية (المستمدة من الحقول المحيطة)، مما يشير إلى وجود مشكلات في نمذجة الخلفية لهذا الهدف تحديداً في بيانات XMM.
الآثار المترتبة على معادلة الحالة (EOS) يدمج المؤلفون قياس نصف قطر PSR J0614−3329 مع القيود الحالية، بما في ذلك:
قياسات الكتلة لـ PSR J1614−2230 و PSR J0348+0432.
قابلية التشوه المد وجذبي من GW170817 و GW190425.
قياسات الكتلة/نصف القطر المستندة إلى NICER لكل من PSR J0030+0451، و PSR J0740+6620، و PSR J0437−4715.
قيد على طاقة التماثل النووي عند كثافة التشبع (S=32±2 MeV).
يتم نمذجة الـ EOS تحت نصف كثافة التشبع باستخدام EOS QHC21 الواقعي مجهرياً (Kojo et al. 2022)، مع التوسع إلى كثافات أعلى باستخدام العمليات الغاوسية (Gaussian processes).
الأهمية والادعاءات يدعي البحث أن إدراج قياس نصف قطر PSR J0614−3329 (الذي يتمركز عند نصف قطر أصغر من قياسات النباضات السابقة ذات الـ ∼1.4M⊙) يؤدي إلى انخفاض طفيف في الكتلة القصوى المستقرة المقدرة لنجم غير دوار (MTOV) وفي نصف القطر المرجعي لنجم بكتلة 1.4M⊙.
الكتلة القصوى: انخفض وسيط MTOV من 2.20 M⊙ (في العمل السابق) إلى 2.15 M⊙ (في هذا العمل).
نصف القطر المرجعي: انخفض وسيط نصف القطر لنجم بكتلة 1.4M⊙ من 12.56 كم إلى 12.44 كم.
يخلص المؤلفون إلى أن هذه النتائج متسقة بشكل عام مع تحليلات NICER السابقة وتستمر في تنقيح فهمنا لمعادلة الحالة للمادة الباردة والمتوازنة كيميائياً عند الكثافات العالية. كما أشاروا صراحة إلى أنه بينما يمكن لتقدير الكتلة العالية لـ PSR J0952−0607 (2.32±0.11M⊙) أن يضع قيوداً إضافية على الـ EOS، فإن عدم اليقين المنهجي في نمذجة النباضات "العنكبوتية" (spider pulsars) يستوجب الحذر قبل إدراجها الفوري. وتؤكد الدراسة على فائدة بيانات NICER المعتمدة على الطور (phase-resolved) في تقييد بنية النجوم النيوترونية دون الاعتماد على فيزياء التنامي.