High-energy Multi-messenger Emission from Galaxy Clusters in the Local Universe

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

بيان المشكلة
لا يزال أصل الخلفيات المنتشرة للنيوترينو وأشعة غاما عالية الطاقة مشكلة لم تُحل بعد في الفيزياء الفلكية. وبينما تُعد النوى المجرية النشطة (AGN) والمجرات ذات النبض النجمي (starburst galaxies) مساهمين معروفين، فإن الدور المحدد لعناقيد المجرات كمصادر لهذه الرسائل الكونية لم يتم ضبطه بشكل كامل بعد. وتتميز عناقيد المجرات بغاز الوسط بين المجري (ICM) الساخن ($\sim 10^8$ كلفن) ومجالات مغناطيسية قوية ($\sim 1$ ميكرو تسلا)، وهي قادرة نظرياً على تسريع الأشعة الكونية (CRs) عبر صدمات الاندماج، وصدمات التراكم، والاضطرابات. يمكن لهذه البروتونات المتسارعة أن تتفاعل مع الوسط بين المجري لإنتاج أشعة غاما ونيوترينو ثانوية عبر العمليات الهادرونية. ومع ذلك، غالباً ما اعتمدت الدراسات السابقة على نماذج مبسطة ذات تماثل كروي أو افترضت حقن الأشعة الكونية حصرياً من مراكز العناقيد، مما قد يؤدي إلى تمثيل خاطئ للتوزيع المكاني وإجمالي تدفق الانبعاثات المنتشرة. علاوة على ذلك، يظل التمييز بين الانبعاث المنتشر للعنقود وانبعاث المصدر النقطي (مثل NGC 1275 أو M87) أمراً صعباً بسبب محدودية الدقة الزاوية في الأرصاد الحالية.

المنهجية
يستخدم هذا العمل إطاراً عددياً صارماً يجمع بين محاكاة الديناميكا المائية المغناطيسية (MHD) عالية الدقة وطرق مونت كارلو لنمذجة الانبعاثات متعددة الرسائل من عناقيد المجرات في الكون المحلي ($z \lesssim 0.1$).

1. محاكاة MHD (مشروع SLOW): تستخدم الدراسة محاكاة "Simulating the LOcal Web" (SLOW)، والتي تتضمن شروطاً ابتدائية مقيدة رصدياً بناءً على كتالوج CosmicFlow. وبخلاف النماذج السابقة، توفر SLOW وصفاً ثلاثي الأبعاد متسقاً ذاتياً لكثافة الغاز، ودرجة الحرارة، وتوزيعات المجال المغناطيسي داخل العناقيد (تحديداً العذراء، وبرسيوس، وكوما)، مما يلتقط المورفولوجيا غير الكروية والبنى التحتية الناجمة عن عمليات الاندماج.
2. حقن الأشعة الكونية: تم اعتماد نموذج حقن مدفوع فيزيائياً حيث يتبع حقن الأشعة الكونية التوزيع المكاني لكثافة غاز العنقود. يأخذ هذا النهج في الاعتبار تسريع الأشعة الكونية في المناطق عالية الكثافة المرتبطة بالاضطرابات والصدمات في جميع أنحاء حجم العنقود، بدلاً من قصر الحقن على قلب العنقود فقط. يتم تحديد إجمالي ميزانية طاقة الأشعة الكونية من خلال النسبة $X_{CR} = E_{CR}/E_{cluster} \simeq 0.01$، وهو ما يتوافق مع الحدود العليا المرصودة.
3. نقل الأشعة الكونية والإنتاج الثانوي: تمت محاكاة انتشار الأشعة الكونية (المفترضة أنها بروتونات) داخل الوسط بين المجري (ICM) والوسط بين المجرات (IGM) باستخدام كود مونت كارلو CRPropa. يقوم الكود بحل معادلة باركر للنقل، مع مراعاة الحمل، والانتشار المكاني، وانتشار الزخم، وفقدان الطاقة. تشمل قنوات التفاعل الرئيسية تصادمات بروتون-بروتون ($pp$) (السائدة تحت$10^{17}$ إلكترون فولت)، وإنتاج البيون الضوئي، والتتالي الكهرومغناطيسي (تشتت كومبتون العكسي مع خلفية الإشعاع الكوني الميكروي EBL والتدفق الكوني الميكروي CMB، وخسائر السينكروترون). يتم توليد أشعة غاما والنيوترينو الثانوية بشكل متسق ذاتياً من هذه التفاعلات.

المساهمات الرئيسية

• نمذجة بيئية واقعية: يسمح دمج محاكاة SLOW MHD بنمذجة انتشار الأشعة الكونية في بيئات عناقيد واقعية وغير كروية، مما يتجاوز افتراضات التماثل الكروي المبسطة للدراسات التحليلية أو شبه التحليلية السابقة.
• مخطط حقن ممتد: يمثل اعتماد نموذج حقن أشعة كونية ممتد مكانياً، يتتبع ملف كثافة الغاز، تحسيناً منهجياً كبيراً مقارنة بالأعمال السابقة التي افترضت حقناً مركزياً فقط. يؤدي هذا إلى تقليل زمن احتجاز الأشعة الكونية الفعال ومعدلات التفاعل في المركز، مما يؤدي إلى ملف انبعاث أكثر انتشاراً.
• تقدير التدفق الشامل: تحسب الدراسة كلاً من تدفقات العناقيد الفردية (العذراء، برسيوس، وكوما) والتدفق المنتشر التراكمي من كامل الجمهرة من العناقيد الضخمة ($M \gtrsim 5 \times 10^{13} M_\odot$) ضمن مسافة 500 ميجابارسيك مجري.

النتائج

• العناقيد الفردية: تدفقات أشعة غاما المنتشرة المتوقعة لعناقيد العذراء، وبرسيوس، وكوما هي أقل بكثير من الحدود العليا الحالية التي أبلغ عنها تعاون LHAASO. كما أن التدفقات المتوقعة أقل من انبعاث أشعة غاما المرصود بواسطة MAGIC للمصدر المركزي NGC 1275 في برسيوس، وهو ما يتوافق مع التوقع بأن MAGIC يرصد مصدراً نقطياً بينما يتنبأ النموذج بانبعاث منتشر.
• التدفق التراكمي: إن تدفق أشعة غاما المتكامل من جميع العناقيد ضمن 500 ميجابارسيك أقل بكثير من خلفية أشعة غاما المنتشرة التي رصدها Fermi-LAT. وتشير الدراسة إلى أنه عند الطاقات تحت $\sim 100$ جيجا إلكترون فولت، تهيمن مصادر أخرى (البلازارات، والنوى المجرية غير المتوافقة محورياً، والمجرات ذات النبض النجمي) على الخلفية.
• انبعاث النيوترينو: تدفق النيوترينو المتوقع لجميع النكهات من جمهرة العناقيد المحلية أقل بكثير من التدفق المنتشر الذي رصده IceCube، ولكنه يقارب الحدود العليا الموجودة لـ IceCube خصيصاً للعناقيد المجرية ($M > 10^{14} M_\odot$). وُجد أن تدفق النيوترينو المتوقع من عنقود العذراء يقارب الحساسية المتوقعة لكاشف KM3NeT/ARCA.
• الخصائص الطيفية: مؤشرات الطيف ($\alpha \approx 2.0 - 2.5$) المتوقعة متسقة مع التوقعات للتسريع بواسطة صدمات الاندماج والتراكم. تشير النتائج إلى أن العناقيد الأقل كتلة تساهم بكفاءة في انبعاث أشعة غاما عند $\sim 10$ جيجا إلكترون فولت، بينما تهيمن الأنظمة الأكثر ضخمة على تدفق النيوترينو عالي الطاقة ($\gtrsim 10$ تيرا إلكترون فولت) بسبب قدراتها الفائقة على احتجاز الأشعة الكونية.

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

تخلص الدراسة إلى أنه بينما من غير المرجح أن تهيمن عناقيد المجرات على إجمالي خلفية النيوترينو عالية الطاقة المنتشرة، إلا أنها قد تشكل جزءاً غير مهمل منها. يوفر هذا العمل إطاراً نظرياً يحفز عمليات البحث المخصصة للنيوترينوهات المنتشرة من عنقود العذراء باستخدام KM3NeT وIceCube-Gen2، ولأشعة غاما المنتشرة باستخدام CTA وLHAASO. ويؤكد المؤلفون أن نتائجهم تتماشى مع آفاق المراصد متعددة الرسائل القادمة، مما يوفر اختباراً في الوقت المناسب لعمليات تسريع الأشعة الكونية الهادرونية في البنى واسعة النطاق.