Search for neutrino emission from blazar $γ$-ray flares accounting for possible neutrino time delays

ملخص تقني: البحث عن انبعاث النيوترينو من توهجات أشعة غاما في البلازارات مع مراعاة التأخيرات الزمنية المحتملة للنيوترينو

بيان المشكلة
لا يزال أصل تدفق النيوترينو عالي الطاقة المكتشف بواسطة IceCube غير محدد بشكل كبير، على الرغم من وجود تلميحات لارتباطات مع مصادر محددة مثل TXS 0506+056 وNGC 1068. لطالما اقتُرحت البلازارات (نوى مجرات نشطة ذات نفاثات نسبية موجهة نحو الأرض) كمصادر للنيوترينو. ومع ذلك، فإن عمليات البحث القياسية عن الارتباطات بين البلازارات والنيوترينو تفترض عادةً التزامن الزمني بين وصول النيوترينو والتوهجات الكهرومغناطيسية الكبرى (تحديداً أشعة غاما).

تشير الاعتبارات النظرية إلى أن هذا التزامن قد لا يتحقق. فإذا تم إنتاج النيوترينوهات عبر تفاعلات $p\gamma$، فإن البروتونات المسؤولة يجب أن تُسرّع إلى طاقات ($E_p \sim 20 E_\nu$) أعلى بكثير من تلك المطلوبة للإلكترونات التي تنتج توهجات أشعة غاما. وبالتالي، فإن مقياس زمن تسريع البروتونات ($t'_{acc}$) و/أو مقياس زمن فقدان طاقة البروتون عبر تفاعلات $p\gamma$ ($t'_{p\gamma}$) في إطار النفاث قد يكون أطول بعدة مراتب قدرية من مقياس زمن تبريد الإلكترونات. وهذا يعني أن النيوترينوهات عالية الطاقة ($E_\nu \gtrsim 100$ TeV) يمكن أن تصل بعد أشهر أو سنوات من ذروة توهج أشعة غاما الرئيسي، وهذا التأخير يعتمد على تأخير زمن الإطار النفاث $t'_{delay}$ الذي يتم تحويله إلى إطار المراقب بواسطة الإنزياح نحو الأحمر ($z$) وعامل دوپلر ($D$).

المنهجية
يقوم المؤلفون بإجراء بحث للارتباط المتبادل بين فهرس IceCat-1 العام لتنبيهات النيوترينو عالية الطاقة وعينتين متميزتين من البلازارات، مع مراع تدخّل تأخير زمني محتمل.

1. مجموعات البيانات:

   • النيوترينوهات: 348 تنبيهاً من فهرس IceCat-1 (باستثناء التنبيهات المتزامنة مع إشارات IceTop).
   • عينات المصادر:
     • H21+: 294 مصدراً من قاعدة بيانات MOJAVE (Homan et al. 2021) ذات الإنزياح نحو الأحمر وعوامل دوپلر معروفة، ومنحنيات ضوئية لـ Fermi-LAT.
     • R24+: 295 مصدراً من فهرس CGRaBS (Rodrigues et al. 2024b) مع عوامل دوپلر مشتقة من نمذجة الطيف الطاقي (SED).
   • المنحنيات الضوئية: منحنيات ضوئية لـ Fermi-LAT بتبويب شهري من المستودع العام (LCR).

2. الإطار الإحصائي:

   • الارتباط المكاني: يُعتبر البلازار مرتبطاً بالنيوترينو إذا وقع ضمن منطقة عدم يقين زاوي موسعة (صندوق خطأ النيوترينو + عدم اليقين المنهجي $\Delta = 0.78^\circ$).
   • مخطط الأوزان: لكل تنبيه، يتم حساب وزن $w_\nu$ بناءً على قوة الإشارة ($s_\nu$)، وعقوبة لعدم اليقين الزاوي الكبير ($\rho_\nu$)، وعامل زمني ($G_{\nu b}$).
   • العامل الزمني: الابتكار الجوهري هو الدالة $G_{\nu b}$، التي تربط زمن وصول النيوترينو $t_\nu$ بزمن الحد الأقصى العالمي للمنحنى الضوئي لأشعة غاما للبلازار قبل الوصول ($t^{\gamma \max}_b$). تقوم الدالة بمراعاة تأخير زمن الإطار النفاث $t'_{delay}$، والذي يتم تحويله إلى إطار المراقب كـ $t_{delay} = t'_{delay}(1+z)/D$. وتخصص الدالة وزناً عالياً إذا تطابق فرق الزمن مع التأخير المتوقع ضمن هامش معين، وتطبق عقوبة غاوسية بخلاف ذلك.
   • المسح: يقوم التحليل بمسح نطاق واسع من تأخيرات الإطار النفاث العالمية $t'_{delay}$ من شهر واحد إلى حوالي $3.5 \times 10^4$ يوماً.
   • اختبار المعنوية: يتم حساب إحصائية الاختبار العالمي (TS) للبيانات الحقيقية ومقارنتها بتوزيع قيم TS الناتجة من 49,999 مجموعة بيانات وهمية حيث يتم تعشية أوقات وصول النيوترينو وصعودها المستقيم. ويتم حساب القيم الاحتمالية (P-values) وتصحيحها لـ "تأثير البحث في كل مكان" (عوامل التجارب) عبر نطاق التأخير الممسوح.

3. اختبارات المتانة:

   • أُجريت اختبارات إضافية باستبعاد وزن تدفق أشعة غاما (الاعتماد فقط على قوة الإشارة والتوقيت) واستخدام وزن بوليني (Boolean) بناءً على ما إذا كان المصدر في حالة توهج (تدفق > المئوي 68) بدلاً من الحد الأقصى العالمي.
   • أُجري اختبار إضافي يتضمن بيانات التدفق الراديوي لعينة H21+.

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

• المعنوية قبل تصحيح التجارب: حدد البحث ارتباطاً قبل التصحيح بمعنوية تقارب $2\sigma$ عند قيم تأخير محددة:
  • لعينة H21+: الحد الأدنى لـ $p \approx 0.037$ عند $t'_{delay} \approx 7.7 \times 10^3$ يوماً.
  • لعينة R24+: الحد الأدنى لـ $p \approx 0.026$ عند $t'_{delay} \approx 1.9 \times 10^3$ يوماً.
• المعنوية بعد تصحيح التجارب: بعد تصحيح التجارب المتعددة التي أجريت عبر النطاق الواسع لـ $t'_{delay}$، لم تكن الارتباطات ذات دلالة إحصائية:
  • $p_{post-trial} \approx 0.14$ (1.06$\sigma$) لـ H21+.
  • $p_{post-trial} \approx 0.11$ (1.21$\sigma$) لـ R24+.
• الارتباطات المهيمنة: جاءت أعلى المساهمات الفردية في TS في كلا الاختبارين من ارتباط التنبيه IC120523B بالبلازار 3C 454.3 (المرتبط بتوهجه في عام 2010) وIC150926A بـ 3C 279. ويعود الاختلاف في أوقات التأخير المثلى بين العينتين بشكل كبير إلى التباينات في تقديرات عامل دوپلر لهذه المصادر المحددة بين فهرسي H21+ وR24+.
• الاختبارات البديلة: لم تسفر الاختبارات التي استخدمت مخططات أوزان بديلة (تجاهل مقدار التدفق أو استخدام حالات التوهج البولينية) واختبار التدفق الراديوي عن أي ارتباطات معنوية، حيث كانت القيم الاحتمية بعد التصحيح أعلى من الاختبار الأساسي.

الأهمية والاستنتاجات
يسجل البحث أول محاولة للبحث عن ارتباطات نيوترينو-بلازار مع مراعاة تأخير زمني عالمي في الإطار النفاث ناتج عن تسريع البروتونات البطيء أو فقدان الطاقة. وبينما وجد التحليل تلميحاً لارتباط بمعنوية $\sim 2\sigma$ عند قيم معينة، فإن هذه المعنوية تتلاشى بعد تصحيح التجارب.

يخلص المؤلفون إلى أنه لا يوجد دليل معنوي على وجود ارتباط بين تنبيهات IceCat-1 والتوهجات الكبرى لأشعة غاما في عينات H21+ أو R24+، حتى عند السماح بفترات تأخير زمنية. ويناقش المؤلفون عدة أسباب محتملة لهذه النتيات الصفرية:

1. قد تكون البلازارات مصادر غير فعالة لنيوترينوهات $\sim 100$ TeV.
2. قد تكون البلازارات فعالة فقط عند طاقات أعلى بكثير ($\gtrsim 100$ PeV)، حيث تقع ذروة طيف النيوترينو.
3. افتراض وجود تأخير زمني عالمي في الإطار النفاث $t'_{delay}$ قد يكون غير صالح بسبب التباينات بين المصدر والآخر في الظروف الفيزيائية.
4. يمكن لعدم اليقين الكبير في تقديرات عامل دوپلر ($D$) بين الفهارس والطرق المختلفة أن يحجب الارتباطات عند تحويل التأخيرات إلى إطار المراقب.
5. قد يكون تسريع البروتونات وتبديدها منفصلين عن العمليات الإلكترونية، مما يعني أن إنتاج النيوترينو ليس مرتبطاً سببياً بتوهجات الأشعة الكهرومغناطيسية المحددة المستخدمة كوكلاء.

يقترح المؤلفون أن التحسينات المستقبلية في الحساسية ستتطلب مجموعات بيانات أكبر (مثل IceCat-2، وفهارس جديدة لعامل دوپلر) وتضمين المنحنيات الضوئية الراديوية، حيث يُعرف أن التوهجات الراديوية تتأخر عن توهجات أشعة غاما.