Temporal evolution of the circumstellar disk orientation in the transient X-ray pulsar GRO J1008-57

ملخص تقني: التطور الزمني لتوجه القرص المحيط بالنجم في النبض الإشعاعي السيني العابر GRO J1008−57

بيان المشكلة
تُظهر الأنظمة الثنائية لنجوم Be/X-ray (BeXRBs) نوعين متميزين من الانفجارات الإشعاعية السينية: انفجارات النوع الأول (العادية) الدورية، والتي تحدث عادةً بالقرب من الحضيض المداري، وانفجارات النوع الثاني (العملاقة) العشوائية. وبينما تُعزى انفجارات النوع الأول عموماً إلى الاقتراب الأقصى للنجم النيوتروني من قرص التجريد المحيط بنجم Be، فإن الأطوار المدارية لهذه الانفجارات لا تظل مستقرة دائماً. ففي بعض الأنظمة، مثل EXO 2030+375، يُظهر الطور المداري لانفجارات النوع الأول تغيرات معقدة وطويلة الأمد، بما في ذلك التناوبات الدورية بين تأخر الطور وتقدمه، والتي غالباً ما تُحفز بواسطة انفجارات النوع الثاني. ويُعتقد أن هذه التغيرات تعكس التطور الهيكلي والفيزيائي لقرص نجم Be. ومع ذلك، لا تزال الآليات التي تقود هذه التحولات في الطور — سواء كانت تتضمن بدارية القرص (precession)، أو تذبذبات اللامركزية (Kozai-Lidov)، أو اضطرابات الكثافة — محل نقاش. ويعد GRO J1008−57، وهو نظام BeXRB ذو دور مداري طويل (249.48 يوماً) وتاريخ من انفجارات النوع الثاني المتعددة، حالة دراسية بالغة الأهمية للتحقق مما إذا كانت آليات تطور الطور الملاحظة في الأنظمة ذات الدورات الأقصر تنطبق على الأنظمة ذات المقاييس الزمنية الديناميكية المختلفة.

المنهجية
حلل المؤلفون بيانات الرصد الإشعاعي السيني طويلة الأمد لـ GRO J1008−57 التي تم الحصول عليها من أدوات Swift/BAT (15–150 keV) و MAXI/GSC (2–30 keV)، والتي تغطي الفترة من 2005 إلى 2026.

1. تحليل المنحنى الضوئي: استُخدمت المنحنيات الضوئية اليومية لتحديد 31 انفجاراً من النوع الأول (Swift) و24 (MAXI)، بالإضافة إلى أربعة انفجارات رئيسية من النوع الثاني حدثت في 2012، 2015، 2017، و2020.
2. نمذجة التوقيت: تم نمذجة الملفات الزمنية لانفجارات النوع الأول باستخدام دالة "الارتفاع السريع والاضمحلال الأسي" (FRED) المبسطة لتحديد وقت البداية ($t_{start}$) ووقت الذروة ($t_{peak}$) بدقة.
3. حساب الطور: حُسبت الأطوار المدارية بالنسبة لأوقات مرور الحضيض المستمدة من تقويم Kühnel et al. (2013).
4. النمذجة التطورية: تمت مطابقة التطور طويل الأمد للأطوار المدارية للانفجارات مقابل ثلاثة نماذج: الانجراف الخطي، الانجراف التربيعي، والنموذج الشبيه بالخطوة (step-like model). واستُخدم معيار "أكايكي للمعلومات" (AIC) إحصائياً لاختيار النموذج الأنسب.
5. الطاقة: تم تقدير إجمالي الطاقة المنبعثة (الفيض) لكل انفجار من النوع الأول عن طريق تكامل منحنيات النموذج الملائمة عبر نافذة قدرها $\pm0.15$ من الدورات المدارية حول الذروة.

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

1. تطور الطور الشبيه بالخطوة: يكشف التحليل أن الأطوار المدارية لانفجارات النوع الأول في GRO J1008−57 لا تتبع انجرافاً خطياً أو تربيعياً مستمراً. بدلاً من ذلك، تظهر تطوراً شبيهاً بالخطوة. تظل الأطوار مستقرة عبر العديد من الدورات المدارية، لكنها تمر بقفزات مفاجئة وصغيرة السعة (تتراوح بين $\sim -003$ إلى $-0.016$ في الطور) تتزامن زمنياً مع كل انفجار من انفجارات النوع الثاني الأربعة المرصودة.
2. اختيار النموذج: قدم النموذج الشبيه بالخطوة ملائمة إحصائية متفوقة مقارنة بالنماذج الخطية والتربيعية لكل من بيانات Swift و MAXI، حيث أشارت قيم $\Delta$AIC إلى أهمية تتجاوز $5\sigma$.
3. دورات الطاقة: يُظهر إجمالي الطاقة المنبعثة في انفجارات النوع الأول نمطاً منتظماً: زيادة تدريجية تؤدي إلى انفجار النوع الثاني، يتبعها انخفاض سريع ثم تعافٍ تدريجي لاحق. وهذا يشير إلى دورات من تراكم كتلة القرص، واستنزافه، ثم إعادة بنائه مدفوعة بانفجارات النوع الثاني.
4. الاستقرار بين القفزات: على عكس نظام EXO 2030+375، الذي يُظهر تحولات طور كبيرة السعة وفترات تعافٍ طويلة الأمد (مئات الأيام)، لا يوجد دليل في GRO J1008−57 على تعافٍ طويل الأمد في الطور أو انجراف مستمر بين أحداث النوع الثاني. يعود النظام إلى حالة طور جديدة ومستقرة فور انتهاء انفجار النوع الثاني مباشرة.

الأهمية والآلية المقترحة
يجادل البحث بأن التطور الشبيه بالخطوة في GRO J1008−57 يتحدى النماذج المقترحة للأنظمة الأخرى من نوع BeXRB، مثل آلية تذبذب كوزاي-ليدوف (KL)، أو تذبذبات القرص أحادية الذراع العالمية، أو بدارية القرص طويلة الأمد.

• رفض تذبذبات KL: يبلغ الجدول الزمني المقدر لتذبذب KL في GRO J1008−57 حوالي 79 عاماً، وهو ما يتجاوز بكثير النطاق الزمني للرصد. علاوة على ذلك، تتنبأ آليات KL عادةً بتحولات طور كبيرة ونمو في اللامركزية، وهو ما لم يُلاحظ هنا.
• رفض البدارية/التذبذبات: إن غياب فترات التعافي طويلة الأمد (المشاهدة في EXO 2030+375) يستبعد كون البدارية طويلة الأمد أو اضطرابات الكثافة العالمية هي المحرك الأساسي.

التفسير المقترح:
يقترح المؤلفون أن السلوك الفريد لـ GRO J1008−57 مدفوع بـ دوره المداري الطويل (249.48 يوماً)، مما يؤدي إلى تفاعلات نادرة بين النجم النيوتروني والقرص المحيط بالنجم.

• دورة إعادة بناء القرص: تُفسر انفجارات النوع الثاني كأحداث تستنزف قرص نجم Be. وبسبب طول الفترة المدارية، يمتلك النظام وقتاً كافياً بين التفاعلات لإعادة بناء كثافة القرص واستعادة هيكله الهندسي ليشبه حالته ما قبل الانفجار.
• القفزات الصغيرة: نتيجة لذلك، بعد كل حدث من النوع الثاني، يستقر القرص في تكوين مختلف قليلاً ولكنه مستقر، مما يسبب تحولاً طفيفاً ومفاجئاً في الطور المداري للانفجارات اللاحقة من النوع الأول. ثم يظل النظام في هذه الحالة المستقرة الجديدة حتى حدث الاستنزاف التالي.
• التباين مع EXO 2030+375: في المقابل، يتميز EXO 2030+375 بدورة مدارية قصيرة ($\sim46$ يوماً)، مما يؤدي إلى تفاعلات متكررة تمنع القرص من الاسترخاء الكامل، مما يسمح للاضطرابات بالتراكم وقيادة تحولات طور مستمرة وكبيرة.

الخلاصة
توسع هذه الدراسة فهم آليات انفجار BeXRB من خلال إثبات أن التطور الزمني لأطوار انفجار النوع الأول ليس عالمياً بل يعتمد على الديناميكيات المدارية للنظام. يوفر التطور الشبيه بالخطوة في GRO J1008−57 دليلاً على دورة استنزاف وإعادة بناء القرص المدفوعة بانفجارات النوع الثاني، حيث تسمح الفترة المدارية الطويلة للقرص بإعادة إنشاء هندسة مستقرة بين الأحداث، مما يؤدي إلى قفزات طور منفصلة بدلاً من الانجراف المستمر.