Cyclotron lines in subcritical X-ray pulsars: Monte Carlo simulations reveal the origin of the observed variability

ملخص تقني: خطوط السيكلوترون في النباضات السينية تحت الحرجة

بيان المشكلة
تُظهر النباضات السينية (XRPs) ميزات تشتت سيكلوتروني رنيني (CRSFs) في أطيافها، والتي تعمل كأدوات تشخيص مباشرة لقوة المجال المغناطيسي للنجم النيوتروني (NS). وتكشف الملاحظات عن تغير ملحوظ في طاقة مركز خط السيكلوترون ($E_{CRSF}$) وعرض الخط ($\sigma_{CRSF}$) كدالة للامعان الأشعة السينية ($L$) وطور النبضة. في نظام اللمعان تحت الحرج ($L \lesssim 10^{37}$ erg s$^{-1}$)، يُلاحظ غالباً وجود ارتباط إيجابي بين $E_{CRSF}$ و $L$، وهو ما يتناقض مع الارتباط السلبي المشاهد في المصادر فوق الحرجة. وبينما تقترح النماذج النوعية أن هذا ينشأ من إزاحة دوپلر للفوتونات التي تتشتت بعيداً عن البلازما المتدفقة للداخل في قناة التراكم، إلا أن معالجة دقيقة وكمية لنقل الإشعاع لهذا التغير كانت تفتقر إليها الدراسات. وتحديداً، يتطلب التفاعل بين التشتت الرنيني، وتأثيرات ضغط الإشعاع على سرعة التدفق، وهندسة الرؤية، محاكاة تفصيلية لتفسير الاتجاهات المرصودة.

المنهجية
قام المؤلفون بتطوير كود مونت كارلو (MC) نسبي لإجراء حسابات تفصيلية لنقل الإشعاع داخل قناة التراكم فوق البقعة الساخنة للنجم النيوتروني. يعتمد النموذج الإطار الفيزيائي المقترح من قبل Mushtukov et al. (2015b) للمصادر تحت الحرجة:

• الهندسة: قناة تراكم أسطوانية شبه نهائية محاذية للمحور المغناطيسي ($B = B\hat{z}$).
• ملفات البلازما: تُشتق تعبيرات تحليلية لكثافة الإلكترونات ($n_e$) والسرعة الكلية ($\beta$) من حفظ الكتلة والزخم، مع دمج تباطؤ البلازما المتدفقة للداخل نتيجة ضغط الإشعاع. يعتمد ملف السرعة على اللمعان بالنسبة لـ لمعان حرج $L^*$.
• نقل الإشعاع: يتم حقن فوتونات البذرة بشكل متماثل المناحي (isotropically) من البقعة الساخنة بطيف جسم أسود. تنتشر حزم الفوتونات عبر البلازما، وتخضع لتشتت كومبتون.
• فيزياء التشتت: يستخدم الكود مقاطع عرضية تقريبية للتشتت المغناطيسي الرنيني متوسط الاستقطاب (Loudas et al. 2021) صالحة لـ $b \ll 1$ و $E \ll m_ec^2$. يتم التعامل مع التشتت في إطار rest frame للإلكترون، مع تطبيق تحويلات لورنتز لمراعاة الحركة الكلية للبلازما.
• الشروط الحدية: تم اختبار مخططين: "الانعكاس" (انعكاس هندسي عند القاعدة، مع حفظ الطاقة) و"التحول الحراري" (الامتصاص وإعادة الانبعاث). تستخدم النتائج الأولية مخطط الانعكاس.
• المعلمات: أُجريت المحاكاة لنجم نيوتروني نموذجي ($M=1.4 M_\odot$, $R=10^6$ cm, $B \approx 10^{12}$ G) وطُبقت على الحالة المحددة للنباض السيني GX 304−1. تم تغيير لمعان التراكم عبر النطاق $L/L^* \in (0, 1)$.

المساهمات والنتائج الرئيسية
تقدم الدراسة أول تحقق كمي باستخدام مونت كارلو لنموذج Mushtukov et al. (2015b)، مما أسفر عن النتائج التالية:

1. مورفولوجيا الطيف: تظهر الأطياف الناشئة ميزة CRSF بارزة وغير متماثلة مصحوبة بجناح أزرق عريض. يظهر خط CRSF إزاحة نحو الأحمر بشكل منهجي بالنسبة لطاقة السيكلوترون الكلاسيكية المحلية ($E_c$). وتنتج هذه الإزاحة نحو الأحمر عن تأثير دوپلر بين إطار المختبر وإطار الإلكترون المتدفق للداخل.
2. الاعتماد على اللمعان: يُظهر كل من $E_{CRSF}$ و $\sigma_{CRSF}$ ارتباطاً إيجابياً مع اللمعان لجميع زوايا الرؤية. مع زيادة $L$، يعمل ضغط الإشعاع على إبطاء التدفق المتدفق للداخل بشكل أكثر فعالية، مما يقلل من السرعة الكلية ($\beta$). هذا الانخفاض في السرعة يقلل من إزاحة دوپلر نحو الأحمر، مما يحرك $E_{CRSF}$ لتقترب من $E_c$. وفي الوقت نفسه، يؤدي زيادة العمق البصري عند اللمعانات الأعلى إلى تعميق الخط وتوسيع الجناح الأزرق.
3. الاعتماد على زاوية الرؤية: تتغير خصائص CRSF بشكل كبير مع زاوية الرؤية $\theta$ (بالنسبة للمحور المغناطيسي).
   • تنخفض $E_{CRSF}$ مع زيادة $\cos \theta$ (أي عندما تصبح الرؤية أكثر "مواجهة للوجه"). ويرجع ذلك إلى اعتماد إزاحة دوپلر على الزاوية بين زخم الفوتون وسرعة التدفق.
   • يزداد $\sigma_{CRSF}$ مع $\cos \theta$. حيث ترصد الرؤية المواجهة للوجه نطاقاً أوسع من تدرجات السرعة على طول خط البصر، مما يؤدي إلى زيادة توسيع دوپلر.
   • بناءً على ذلك، يتنبأ النموذج بوجود ارتباط عكسي بين $E_{CRSF}$ و $\sigma_{CRSF}$ خلال دورة النبض (التباين المستخلص عبر الطور)، بينما تُظهر الأطياف متوسطة الطور ارتباطاً إيجابياً مع اللمعان.
4. التطبيق على GX 304−1: عند تطبيقه على البيانات الرصدية لـ GX 304−1، نجح النموذج في إعادة إنتاج الارتباطات الإيجابية لـ $E_{CRSF}$ و $\sigma_{CRSF}$ مع اللمعان عبر ما يقرب من رتبة مقدار واحدة. ويتحقق أفضل مطابقة لزوايا الرؤية حيث تُرى قناة التراكم بشكل أساسي من الجانب ($\cos \theta \approx 0.0-0.1$)، وهو ما يتوافق مع تكوينات هندسية محددة (مثل دوران محاذٍ تقريباً عند ميل مرتفع أو دوران متعامد تقريباً عند ميل منخفض).

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

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

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