Study of Flat Spectrum Radio Quasars and BL Lacertae Objects as Sources of Diffusive Ultra High-Energy Cosmic Rays

ملخص تقني: دراسة الكوازارات ذات الطيف الراديوي المسطح وأجسام بلا لار (BL Lac) كمصادر للأشعة الكونية فائقة العالية الطاقة المنتشرة

بيان المشكلة
لا يزال أصل الأشعة الكونية فائقة العالية الطاقة (UHECRs, $E \gtrsim 10^{17}$ eV) يمثل قضية غير محلولة بشكل كبير في الفيزياء الفلكية الحديثة. وبينما يُقترح الأصل خارج المجري بناءً على السمات الطيفية، إلا أن تحديد المصادر بدقة يظل أمراً صعباً بسبب الانحرافات الناتجة عن المجالات المغناطيسية المجرية وخارج المجرية. تكشف البيانات الحديثة عالية الدقة من مرصد "بير Pierre Auger" وتلسكوب "Array TA" عن سمات طيفية محددة (تصلب عند "الكاحل/the ankle"، وتدرج في الانحدار، وكبح GZK) واتجاهات في التركيب (انتقال من النوى الخفيفة إلى النوى الأثقل مع زيادة الطاقة). وبينما تُعد المجرات الراديوية القريبة مرشحة قوية، فإن المساهمة المحتملة للنوى المجرة النشطة (AGN)، وتحديداً الكوازارات ذات الطيف الراديوي المسطح (FSRQs) وأجسام بلا لار (BL Lac)، تتطلب اختباراً صارماً مقابل القيود الرصدية المتعلقة بالتدفق، والتركيب، وتباين الخواص.

المنهجية
يتقصى المؤلفون فرضيات FSRQs وBL Lacs كمصادر محتملة للأشعة الكونية فائقة العالية الطاقة المنتشرة باستخدام نموذج انتشار شامل ضمن نموذج $\Lambda$CDM الكوني القياسي.

1. نماذج تطور المصادر:

   • يستخدم الدراسة دالات تطور الكثافة المعتمدة على اللمعان (LDDE) الواقعية والمستمدة من دالات لمعان أشعة غاما (LF) المرصودة من بيانات Fermi-LAT.
   • FSRQs: تم نمذجتها باستخدام دالة LDDE حيث يعتمد ذروة الإزاحة الحمراء لكثافة المصدر على اللمعان ($z_c \propto L^\alpha$).
   • BL Lacs: تم نمذجتها باستخدام توزيع أصلي لمؤشرات فوتونات (غاوسي) مرتبط باللمعان، مدمجاً مع نفس صيغة LDDE.
   • هذه النماذج تستبدل أشكال التطور الظاهراتية البسيطة (مثل القوى في $1+z$) بتطورات مدفوعة فيزيائياً بناءً على بيانات أشعة غاما المرصودة.

2. فيزياء الانتشار:

   • يتم نشر الأشعة الكونية فائقة العالية الطاقة في الفضاء خارج المجري مع مراعاة فقدان الطاقة الناتج عن التفاعلات مع الخلفية الميكروية الكونية (CMB) والخلفية الإشعاعية خارج المجرية (EBL).
   • العمليات المتضمنة: إنتاج أزواج بيته-هيتلر (Bethe-Heitler)، وإنتاج البيونات الفوتونية (للبروتونات)، والتفكك الضوئي (للنوى). تستخدم الدراسة تقريب فقدان الطاقة المستمر (CEL) للتفكك الضوئي، مع تتبع احتمالية بقاء النوى الأولية (p, He, N, Si, Fe) دون حل الشلالات الثانوية صراحةً.
   • الانتشار المغناطيسي: يحدث الانتشار في مجال مغناطيسي بين مجري منتظم ومضطرب (IGMF) بقوة جذر متوسط مربعات $B_{rms} = 1$ nG وطول ترابط $l_c = 0.5$ Mpc. يتم حساب معامل الانتشار المعتمد على الطاقة $D(E)$ بناءً على اضطراب كولموغوروف (Kolmogorov)، مع الفصل بين أنظمة الانتشار الرنيني وغير الرنيني.

3. حساب التدفق والملاءمة:

   • يتم حساب التدفق المنتشر عن طريق تكامل انبعاث المصدر عبر الإزاحة الحمراء واللمعان، مع تطبيق نظرية الانتشار لجمع المساهمات من توزيع منتظم لمصادر اصطناعية.
   • يتضمن النموذج تركيباً نووياً مختلطاً (p, He, N, Si, Fe) حيث تُعامل الكسور كمعاملات حرة.
   • يرتبط تطبيع انبعاث المصدر بلمعان أشعة غاما عبر عامل تحميل الباريون $\kappa = L_{CR}/L_\gamma$.
   • يتم ملاءمة المعلمات (مسافة فصل المصدر $d_s$، كفاءة التسريع $\eta$، والكسور النووية $\epsilon_Z$) مع البيانات الرصدية من Auger وTA (التدفق ومتوسط عمق الشلال $\langle X_{max} \rangle$).
   • يتم إجراء مقارنة إحصائية للنماذج باستخدام $\chi^2$ المختزل، ومعيار أكايكي (AIC)، ومعيار بايز (BIC).
   • يتم إنشاء خرائط كاملة للسماء للتدفق المتكامل ($10-100$ EeV) وتباين الخواص ثنائي القطب باستخدام إطار عمل HEALPix.

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

1. التدفق والتركيب:

   • يمكن لكل من سيناريوهات FSRQ وBL Lac إعادة إنتاج طيف طاقة الأشعة الكونية فائقة العالية وتدفق $\langle X_{max} \rangle$ مع تركيبات مختلطة.
   • أفضل كفاءة تسريع ملائمة هي $\eta \approx 0.1$، وهي متسقة مع المحاكاة العددية لتسريع النفاثات النسبية.
   • عامل تحميل الباريون المطلوب هو $\kappa \approx 5.65$، مما يشير إلى أن قدرة الأشعة الكونية هي قدرة مجدية طاقياً (فقط بضع مرات قدر لمعان أشعة غاما).
   • سيناريو FSRQ: يتطلب مسافة فصل للمصدر $d_s \approx 8 \pm 1$ Mpc لضبط التدفق. وهذا يعني أن كثافة المصدر أعلى بكثير من التقديرات الرصدية لدالات لمعان Fermi-LAT.
   • سيناريو BL Lac: يتطلب $d_s \approx 160 \pm 21$ Mpc، وهو ما يتوافق مع التقديرات الرصدية لكثافة BL Lac.

2. المقارنة الإحصائية:

   • قيم $\chi^2$ المختزل للبيانات المجمعة لـ Auger وTA هي 2.37 (FSRQ) و2.55 (BL Lac).
   • تفضل معايير اختيار النموذج (AIC وBIC) نموذج BL Lac ($\Delta$AIC = 4.17, $\Delta$BIC = 4.83)، مما يوفر دليلاً متوسطاً على أن BL Lac يقدم وصفاً أفضل للبيانات من FSRQ.

3. قيود تباين الخواص:

   • BL Lacs: يتوافق تباين الخواص ثنائي القطب المتوقع باستخدام معلمات ضبط التدفق ($z_{max}=0.5$) جيداً مع ملاحظات Auger.
   • FSRQs: سعة ثنائي القطب المتوقعة أعلى بكثير مما هو مرصود عند استخدام نفس المعلمات. ولمطابقة التباين المرصود، يجب تقييد نطاق الإزاحة الحمراء اصطناعياً إلى $z_{max} \approx 0.01$، وهو ما يتعارض مع متطلبات ضبط التدفق والتطور الكوني القوي لـ FSRQs.

4. التفكك الضوئي:

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

الأهمية والاستنتاجات
تخلص الورقة إلى أنه بينما يمكن لكل من FSRQs وBL Lacs إعادة إنتاج تدفق الأشعة الكونية فائقة العالية، فإن سيناريو BL Lac هو المفضل بقوة عند النظر في مجموعة القيود الرصدية الكاملة.

• تُستبعد FSRQs كمصادر مهيمنة للأشعة الكونية فائقة العالية المنتشرة لأن السيناريو يتطلب كثافة مصادر عالية بشكل غير واقعي (مسافة $d_s$ صغيرة) ويتنبأ بتباين خواص ثنائي القطب يتجاوز الحدود المرصودة ما لم يتم قطع تعداد المصادر اصطناعياً عند إزاحات حمراء منخفضة.
• تظل BL Lacs متسقة مع الملاحظات، حيث توفر وصفاً ذاتي الاتساق للتدفق، والتركيب، وتباين الخواص ثنائي القطب مع كثافة مصدر وتطور إزاحة حمراء يتوافق مع بيانات Fermi-LAT.

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