Probing evolution of Long GRB properties through their cosmic formation history aided by Machine Learning predicted redshifts

ملخص تقني: سبر تطور خصائص انفجارات غاما المضيئة الطويلة من خلال تاريخ تكوينها الكوني بمساعدة الإزاحة الحمراء المتوقعة عبر تعلم الآلة

بيان المشكلة
تُعد انفجارات غاما المضيئة (GRBs) مسابير حاسمة لتاريخ تشكل النجوم الكوني (SFRH) الممتد حتى عصر إعادة التأين. ومع ذلك، فإن هناك عائقاً كبيماً في الدراسات الحالية يتمثل في أن حوالي 80% من انفجارات غاما تفتقر إلى إزاحة حمراء ($z$) مقاسة، حيث يتطلب التحديد الطيفي أو الضوئي متابعة سريعة لانبعاثات الوميض اللاحق (afterglow)، وهو أمر لا يمكن تحقيقه غالباً. وتواجه الطرق الحالية لتقدير الإزاحة الحمراء باستخدام الارتباطات الداخلية لانفجارات غاما جدلية "الدورانية" (circularity arguments) (بافتراض كونية معينة لاستخلاص الارتباط، ثم استخدام هذا الارتباط لإيجاد $z$) والشكوك المنهجية. وبناءً على ذلك، فإن مجموعات البيانات المستخدمة لاستخ derivation كثافة معدل انفجارات غاما المضيئة الطويلة (LGRB-RD) هي مجموعات صغيرة ومن المحتمل أن تكون منحازة، مما يعيق الفهم الواضح لتطور انفجارات غاما وعلاقتها بمعدل تشكل النجوم (SFRD).

المنهجية
تستخدم هذه الدراسة عينة من انفجارات غاما المضيئة الطويلة (LGRBs) ذات إزاحات حمراء متوقعة بواسطة نماذج تعلم الآلة (ML)، والتي تم تدريبها على معاملات الوميض الفوري والوميض اللاحق للأشعة السينية (تحديداً ارتباطات هضبة الأشعة السينية) دون الاعتماد على افتراضات كونية. يتبع التحليل الخطوات التالية:

1. اختيار البيانات: أنشأ المؤلفون مجموعتين من البيانات:

   • العينة المتوقعة (PS): 228 انفجار غاما مضيء طويل يحتوي على هضاب أشعة سينية ولكن دون إزاحة حمراء مرصودة، حيث تم توقع $z$ باستخدام نموذج تعلم الآلة من ناريندرا وآخرون (2025) (N25). تم استبعاد القيم الشاذة ذات الإزاحة الحمراء أو خطأ اللمعان العالي.
   • العينة المجمعة (CS): 446 انفجار غاما مضيء طويل تضم العينة المتوقعة (PS) بالإضافة إلى 218 انفجار غاما مضيء طويل ذات إزاحة حمراء طيفية معروفة وهضاب أشعة سينية ("العينة المرصودة" أو OS).
   • الاستقطاعات: تم تصفية كلتا العنتين لإزالة انفجارات غاما التي تقل عن عتبة لمعان الأداة (تم تحديدها عبر اختبارات كولموغروف-سميرنوف) وانفجار غاما 060614 (وهو انفجار طويل من أصل اندماجي).

2. حساب كثافة المعدل:

   • طبق المؤلفون طريقة إيفرون-بيتروسيان (EP) لتصحيح تأثيرات الاختيار واستخلاص توزيع كثافة المعدل التراكمي ($\sigma$).
   • لمراعاة التطور المحتمل في خصائص انفجارات غاما (مثل اللمعان)، تم تطبيق عامل تصحيح $L' = L/(1+z)^k$. تم تحديد المعامل $k$ من خلال إيجاد القيمة التي يكون عندها إحصاء كيندال-$\tau$ صفراً، مما يضمن استقلال اللمعان والإزاحة الحمراء.
   • يتم حساب كثافة معدل انفجارات غاما المضيئة الطويلة ($\rho$) عن طريق اشتقاق التوزيع التراكمي والقياس بعوامل تأخذ في الاعتبار حساسية الأداة، ومجال الرؤية، ووقت التشغيل، وإجمالي عدد انفجارات غاما المكتشفة.

3. المقارنة النظرية:

   • تم بناء كثافة معدل نظري ($\rho_{theor}$) عن طريق قياس معدل تشكل النجوم (SFRD) (ماداو وديكنسون 2014) بعوامل تحويل (كسر الكولابسار، أنواع المستعرات العظمية، التوجيه الشعاعي) وحد $[\text{term}]$ تطوري $(1+z)^\delta$.
   • تختبر الدراسة ثلاثة سيناريوهات تطورية: لا تطور ($\delta=0$)، تطور التوجيه الشعاعي ($\delta=1.6$)، وحد تطوري عام حيث تتراوح $\delta$ بين -4 و 4.

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

• أهمية التطور: يكشف التحليل عن دليل كبير على وجود تطور في خصائص انفجارات غاما. قيم $k$ الأفضل ملاءمة هي $2.2^{+0.8}_{-0.7}$ للعينة المتوقعة (PS) و $2.5^{+0.4}_{-0.4}$ للعينة المجمعة (CS)، وهو ما يقابل مستويات دلالة تبلغ $\sim3\sigma$ و $\sim6\sigma$ على التوالي، مقارنة بحالة عدم التطور ($k=0$).
• اتجاهات كثافة المعدل:
  • بلغت كثافات المعدل في كلتا العنتين ذروتها في نطاق الإزاحة الحمراء $z \sim 1-2$.
  • تصل العينة المتوقعة (PS) إلى ذروتها عند إزاحة حمراء أعلى من العينة المجمعة (CS)، ويُعزى ذلك إلى انحياز نموذج تعلم الآلة نحو متوسط توزيع التدريب، والذي تم تصحيحه لاحقاً.
  • لوحظ فائض في الإزاحة الحمراء المنخفضة ($z < 1$) في العينة المجمعة (CS)، ويعزوه المؤلفون إلى مساهمات محتملة من عمليات الاندماج، أو انفجارات غاما منخفضة اللمعان (LLGRBs)، أو عيوب في الملاءمة (fitting artifacts).
• نماذج التطور:
  • لم ينجح أي نموذج تطوري منفرد (قانون القوة في $1+z$) في تفسير كثافة المعدل عبر كامل نطاق الإزحة الحمراء لأي من العنتين.
  • بينما يمكن لقيم محددة من $\delta$ أن تطابق اتجاهات الإزحة الحمراء المنخفضة أو العالية، إلا أن أحداً منها لم يستطع استيعاب التوزيع بالكامل. بالنسبة للعينة المجمعة (CS)، فإن $\delta = -0.3$ يناسب اتجاه الإزحة الحمراء المنخفضة، بينما $\delta = -0.9$ يناسب الذروة واتجاه الإزاحة الحمراء العالية.
  • تقع ذروة كثافة معدل العينة المجمعة (CS) عند $z \sim 1.6$ (زمن العودة للخلف $\sim 9.5$ مليار سنة)، بينما تبلغ ذروة SFRD الخاصة بـ MD14 عند $z \sim 1.9$ (زمن العودة للخلف $\sim 10.1$ مليار سنة). يشير المؤلفون إلى أن هذا التفاوت يقع ضمن هوامش الخطأ، لكنهم يقترحون تأثيرات محتملة من التفاعلات الثنائية التي تمد أعمار الأسلاف.

الأهمية والادعاءات
يزعم البحث أن دمج الإزاحات الحمراء المتوقعة بواسطة تعلم الآلة يمثل "خطوة حاسمة للأمام" في توسيع مجموعة البيانات المتاحة لدراسة تطور انفجارات غاما. ومن خلال تعزيز حجم العينة، أثبت المؤلفون ما يلي:

1. قضاษيا تجانس العينات: إن عدم قدرة مصطلحات تطور قانون القوة البسيطة على تفسير كثافة المعدل عبر جميع الإزحات الحمراء يشير إما إلى وجود مصطلح تطوري أكثر تعقيداً أو عدم تجانس في العينات (على سبيل المثال، مزيج من أصول الكولابسار والاندماج ضمن تصنيف الانفجارات الطويلة).
2. تتبع SFRD: بينما تتبع انفجارات غاما المضيئة الطويلة بشكل عام الـ SFRD حتى $z \sim 6$، فإن التطور المحدد لخصائصها (التوجيه، المعدنية، نوع السلف) يُدخل انحرافات لا يمكن لعمليات القياس البسيطة حلها بالكامل.
3. الاتجاهات المستقبلية: يؤكد المؤلفون أن تحليل تعلم الآلة الحالي هو الأكثر دقة حتى الآن، ولكن التحسينات المستقبلية تتطلب دمج عينات الهضبة الضوئية والأشعة السينية للتدريب، واستخدام بيانات من المهمات القادمة (SVOM، Rubin، THESEUS، Hi-Z GUNDAM) لزيادة عينة الإزاحات الحمراء المعروفة. كما يحذرون من أن دمج عينات متعددة الأطوال الموجية يتطلب حذراً بسبب الاختلافات المحتملة في مواقع الانبعاث وتأثيرات الحجب الغباري.

تخلص الدراسة إلى أنه بينما تُعد كثافات المعدلات المستمدة موثوقة في نطاق الإزاحة الحمراء المتوسطة، فإن عدم التطابق عند الإزحات الحمراء المنخفضة والعالية يسلط الضوء على الحاجة لمزيد من الاستقصاء حول تنوع الأسلاف والأصول الفيزيائية لميزات كثافة المعدل المرصودة.