Exploring the High-Redshift 21-cm Signal via Self-Consistent Simulations using Artificial Neural Network Emulation

ملخص تقني: استكشاف إشارة الـ 21 سم عند الانزياح الأحمر العالي باستخدام محاكاة متسقة ذاتياً عبر محاكاة الشبكة العصبية الاصطناعية

بيان المشكلة
يعد فهم الفجر الكوني (Cosmic Dawn - CD)، وهي حقبة النجوم والمجرات الأولى ($z \sim 15-60$)، أمراً حاسماً لتتبع التطور الكوني. وبينما تتقدم المسابير الرصدية مثل تلسكوب جيمس ويب الفضائي (JWST) والمصفوفات الراديوية (HERA، LOFAR، MWA)، فإن تفسير بيانات خط الـ 21 سم يتطلب نماذج نظرية قوية. تواجه النهج الحالية معضلة ثنائية: النماذج التحليلية المبسطة تتسم بالكفاءة الحسابية ولكنها تعتمد على فرضيات مقيدة، بينما النماذج الهيدروديناميكية المعقدة توفر دقة عالية ولكنها محدودة في أحجام صغيرة ونطاقات بارامترية ضيقة بسبب التكلفة الحسابية. علاوة على ذلك، تفتقر العديد من النماذج شبه العددية إلى معالجات متسقة ذاتياً لتكون النجوم (SF) على المقاييس الصغيرة والمعايرة بالديناميكا المائية، وتاريخ اندماج هالات المادة المظلمة (DM) العشوائي، والتغذية الراجعة الإشعاعية واسعة النطاق. تعالج هذه الورقة الحاجة إلى إطار عمل شبه عددي واسع النطاق ومتسق ذاتياً لسد هذه الفجوة لتوفير تنبؤات دقيقة لإشارة الـ 21 سم العالمية وطيف القدرة.

المنهجية
يقدم المؤلفون إطار عمل مطور للمحاكاة شبه العددية، بناءً على عمل Feathers et al. (2025, F25). تتضمن المنههجية الجوهرية ما يلي:

1. حجم المحاكاة: يتم توليد حجم واسع النطاق قدره $192^3$ ميجابارسيك$^3$ (مكون من $64^3$ خلية بحجم $(3 \text{ Mpc})^3$) باستخدام 21cmvFAST. يتضمن هذا الحجم تقلبات مكانية في كثافة المادة ($\delta$) وسرعات الباريون-المادة المظلمة ($v_{bc}$)، وهي ضرورية لتنبؤات تكون النجوم (SF) الدقيقة.
2. تطور هالات المادة المظلمة: يستخدم إطار العمل صيغة "إكستند بريس-شيسر" (EPS) لبناء أشجار اندماج هالات المادة المظلمة العشوائية. يتم توليد مكتبة من $\sim 10^4$ شجرة اندماج لـ 40 فئة كتلة هالة ($10^{5.6} - 10^{9.5} M_\odot$). تُخصص هذه الأشجار للخلايا بناءً على الكثافة المحلية الزائدة، حيث تعمل كدعامة لتكون النجوم.
3. محاكاة الشبكة العصبية (NN Emulation): للتغلب على القيود الحسابية لتشغيل نموذج شبه تحليلي فرعي (SAM) معقد لكل خلية، استخدم المؤلفون محاكيات الشبكات العصبية الاصطناعية.
   • بيانات التدريب: تم تدريب الشبكات العصبية على نموذج SAM معاير لمحاكاة هيدروديناميكية حديثة: محاكاة Renaissance (O'Shea et al. 2015) ومحاكاة AEOS (Brauer et al. 2025).
   • وصف تكون النجوم (SF): يتضمن النموذج المرجعي وصفات تكون النجوم من Hazlett et al. (2025a,b)، ويشمل:
     • تكون نجوم PopIII: تمت نمذجتها باستخدام دالة التوزيع التراكمي للكتل النجمية، مما ينتج عنه متوسط أربعة نجوم بكتلة $35.6 M_\odot$ لكل هالة عند تجاوز عتبة الكتلة الحرجة ($M_{crit}$).
     • تكون نجوم PopII: تمت نمذجتها كعملية ذات مرحلتين: مرحلة "انفجارية" أولية تليها مرحلة "مستقرة" بمجرد وصول درجة حرارة الهالة الفيريال إلى $T_{vir, steady} \approx 1.76 \times 10^4$ كلفن. يتضمن ذلك فترة سكون ($t_{quies}$) تتبع الانفجارات.
     • الكتلة الحرجة ($M_{crit}$): تُحسب مع مراعاة تغذية "ليمان-ويرنر" (LW) الراجعة والحماية الذاتية لـ $H_2$ (Kulkarni et al. 2021).
   • الكفاءة: تقلل محاكاة الشبكة العصبية التكلفة الحسابية لتوليد تنبؤات كثافة معدل تكون النجوم (SFRD) بمقدار $\sim 200$ ضعف مقارنة بالنموذج الكامل لـ SAM، مما يسمح بمحاكاة الحجم الكامل في حوالي 21 ساعة من وحدة المعالجة المركزية (CPU-hours).
4. حساب إشارة الـ 21 سم: تتم معالجة كثافة معدل تكون النجوم (SFRD) المحاكية لاحقاً لحساب درجة حرارة السطوع للـ 21 سم ($\delta T_b$) وطيف القدرة ($\Delta^2(k)$). يتضمن ذلك:
   • حساب خلفية "ليمان-$\alpha$" ($J_\alpha$) لتحديد معامل اقتران "ووثويسن-فيلد" (W-F) ($x_\alpha$).
   • حل درجة حرارة السبن ($T_S$) والحرارة الحركية ($T_K$) للوسط بين المجري (IGM)، مع مراعاة تسخين "كومبتون" وتسخين الأشعة السينية (المفترض أنه يتتبع عملية تكون النجوم).
   • الحل التكراري لـ $T_S$ لتحديد $\delta T_b$ وطيف القدرة الخاص بها.

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

• معايرة مبتكرة: هذه هي أول محاكاة شبه عددية واسعة النطاق توحد بشكل متسق ذاتياً بين أشجار اندماج هالات المادة المظلمة العشوائية ووصفات تكون النجوم المعايرة صراحةً بالمحاكاة الهيدروديناميكية (Renaissance و AEOS).
• نمذجة PopII ثنائية المرحلة: يقدم إطار العمل انتقالاً من الحالة "الانفجارية" إلى الحالة "المستقرة" لتكون نجوم PopII، والمعايرة بالنتائج الهيدروديناميكية، وهو ما يختلف بشكل كبير عن نماذج الكفاءة الثابتة السابقة.
• تاريخ الاندماج العشوائي: يسمح تضمين تاريخ اندماج هالات المادة المظلمة العشوائي بمعالجة واقعية للانتقال من تكون نجوم PopIII إلى PopII، بما في ذلك فترة التأخير ($t_{delay}$) وتأثيرات $v_{bc}$ على $M_{crit}$.
• الكفاءة الحسابية: تسمح بنية الشبكة العصبية المحسنة ومجموعة التدريب المختصرة بتشغيل عمليات متعددة بفرضيات فيزيائية فلكية متغيرة (مثل كتل PopIII، و $t_{delay}$، و $M_{crit}$) بكفاءة عالية.

النتائج
يعرض المؤلفون النتائج من نموذج مرجعي وعدة تجارب بديلة (متغيرة في كتلة PopIII، و $M_{crit}$، و $t_{delay}$):

• الإشارة العالمية ($\delta T_b$):
  • يتنبأ النموذج المرجعي بحد أدنى لامتصاص $\delta T_b \approx -146$ mK عند $z \approx 20.25$.
  • تأثير PopIII مقابل PopII: تهيمن تكون نجوم PopIII على إشارة الـ 21 سم عند $z \gtrsim 34$ وعلى المقاييس الكبيرة. بينما تهيمن تكون نجوم PopII عند $z \lesssim 20$ وعلى المقاييس الأصغر ($k \gtrsim 0.2$ Mpc$^{-1}$ عند $z \sim 20-34$).
  • فترة التأخير ($t_{delay}$): تؤثر فترة التأخير بين PopIII و PopII بشكل كبير على الإشارة. يؤدي التأخير الأطول (مثلاً 100 مليون سنة ضوئية) إلى فجوة امتصاص متأخرة وأعمق وأضيق مقارنة بفترة التأخير المرجعية البالغة 10 ملايين سنة ضوئية.
  • المقارنة: يعطي النموذج المرجعي فجوات امتصاص أبكر وأعمق من العديد من النماذج السابقة (مثل Muñoz et al. 2022) بسبب الجمع بين أشجار اندماج المادة المظلمة، والحماية الذاتية لـ $H_2$، وكفاءات تكون النجوم المعايرة.
• طيف القدرة ($\Delta^2(k)$):
  • يمكن رصد أطياف القدرة بواسطة HERA بدءاً من $z \approx 25$ عند $k \sim 0.2$ Mpc$^{-1}$، وتتوسع لتصل إلى $k \lesssim 0.9$ Mpc$^{-1}$ بحلول $z=15$ تحت افتراضات متوسطة للمقدمات (1080 ساعة رصد).
  • تظهر النماذج ذات تاريخ الاندماج العشوائي قدرة معززة في المقاييس الكبيرة عند الارتفاعات الحمراء العالية مقارنة بالنماذج التي تفتقر إليها (مثل نموذج Muñoz et al. 2022)، مدفوعة بتقلبات $v_{bc}$ التي تؤثر على $M_{crit}$.
• دقة الشبكة العصبية: تعيد المحاكيات إنتاج نتائج SAM بدقة متوسطة تبلغ $\sim 8.3\%$ لـ PopII و $\sim 5.0\%$ لـ PopIII عبر جميع الارتفاعات الحمراء والنماذج.

الأهمية والادعاءات
تدعي الورقة أن الجمع بين معايرة تكون النجوم الهيدروديناميكية، ونموذج كتلة الهالة الحرجة الذي يأخذ في الاعتبار الحماية الذاتية لـ $H_2$، وتاريخ اندماج هالات المادة المظلمة العشوائي، يؤدي إلى تكون نجوم أبكر وكثافات أعلى لمعدل تكون النجوم (SFRDs) عبر الفجر الكوني مقارنة بالأعمال شبه العددية السابقة.

ويؤكد المؤلفون ما يلي:

1. القابلية للرصد: إشارة الـ 21 سم المرجعية قابلة للرصد المحتمل بواسطة HERA عند $z \lesssim 25$.
2. القيود على الفيزياء: إن عدم الرصد عند $z \gtrsim 20$ قد يشير إلى أن فترة التأخير ($t_{delay}$) بين PopIII و PopII تبلغ $\gtrsim 30$ مليون سنة ضوئية.
3. الضرورة المنهجية: تتطلب المعالجة الصحيحة للانتقال بين PopIII و PopII تضمين تاريخ اندماج هالات المادة المظلمة، حيث تفشل النماذج التحليلية المبسطة (مثل تكاملات HMF) في التقاط العشوائية وتوقيت هذا الانتقال.
4. التمايز: بينما قد يكون التمييز بين الاختلافات الطفيفة في كفاءة تكون النجوم صعباً، فإن حساسية طيف القدرة لـ HERA كافية للتمييز بين النماذج ذات قيم $t_{delay}$ المختلفة، وللتمييز بين نموذج المؤلفين (الذي يستخدم أشجار الاندماج والمعايرة الهيدروديناميكية) والنماذج التي تفتقر إلى هذه الميزات (مثل واقعة Muñoz et al. 2022).

يضع هذا العمل نفسه كأداة جديدة ومكملة للأكواد شبه العددية الموجودة (مثل 21cmFAST و Zeus21)، مقدمًا نهجًا أكثر استنادًا إلى الفيزياء للتنبؤ بإشارة الـ 21 سم للحملات الرصدية الحالية والمستقبلية.