Modelling the expulsion of baryons from haloes from first principles: the role of feedback and of the cosmological constant

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

1. بيان المشكلة

يعد توزيع الباريونات داخل هالات المادة المظلمة (DM) قضية مركزية في تكوين المجرات. وبينما تحتفظ عناقيد المجرات الضخمة بكسر كتلة باريونية ($f_b$) يتوافق مع المتوسط الكوني ($f_{b, \text{cosmic}} = \Omega_b/\Omega_m$)، فإن الهالات ذات الكتلة الأقل ($M_{200} \lesssim 10^{12.8} M_\odot$) تعاني من نقص ملحوظ في الباريونات. ويُعزى هذا "البروز في الباريونات المفقودة" إلى التغذية الراجعة الفيزيائية الفلكية (الرياح النجمية، المستعرات العظمى، ونفاثات النوى المجرية النشطة - AGN) التي تطرد الغاز إلى ما وراء نصف قطر فيريال (virial radius).

حددت الدراسات العددية السابقة (مثل Ayromlou et al. 2023) علاقة تجريبية بين نصف قطر الإغلاق ($R_{\text{closure}}$) — وهو المسافة من مركز الهالة حيث يعود كسر الباريونات المحتواة إلى قيمته الكونية — وبين كتلة الهالة وكسر الباريونات. ومع ذلك، كانت هذه العلاقات مجرد ملاءمات تجريبية تفتقر إلى اشتقاق فيزيائي من المبادئ الأولى. علاوة على ذلك، لم يتم التحقيق بشكل منهجي في اعتماد هذه الظاهرة على المعلمات الكونية، وتحديداً الثابت الكوني ($\Lambda$)، حيث أن معظم المحاكاة تثبت الكونية لتطابق الكون المرصود.

2. المنهجية

تطوير النموذج التحليلي

اشتق المؤلفون نموذجاً تحليلياً من المبادئ الأولى لكثافة الغاز ونصف قطر الإغلاق من خلال موازنة ثلاث قوى تؤثر على الباريونات:

1. الانهيار الثقالي: الجاذبية الذاتية للهالة.
2. الضغط الهيدروستاتيكي: تدرج ضغط الغاز.
3. التوسع الكوني/الطاقة المظلمة: القوة الطاردة الناتجة عن $\Lambda$.

الافتراضات والاشتقاقات الرئيسية:

• توازن القوى: يفترض النموذج تناظراً كروياً ويحل معادلة الحركة لقطعة من الغاز، مع دمج حد الثابت الكوني ($\Lambda c^2 r / 3$).
• معادلة الحالة: يتم التعامل مع الغاز كسائل متعدد التراتبية (polytropic fluid) مثالي ($P \propto \rho^{n+1}$).
• ملف الكثافة: في الحالة التي تكون فيها كثافة الغاز $\rho \ll \rho_{\text{crit}}$ (بعيداً عن الهالة)، يشتق المؤلفون ملف كثافة بقانون القوة: $\rho(r) \propto r^{-\eta}$.
• تعريف نصف قطر الإغلاق: يُعرف $R_{\text{closure}}$ بأنه نصف القطر حيث يكون $M_b(<R)/M(<R) = \Omega_b/\Omega_m$.
• الصيغة النهائية: من خلال دمج ملف الكثافة من نصف قطر فيريال ($R_{200}$) إلى نصف قطر الإغلاق ($R_{\text{closure}}$)، اشتقوا تعبيراً (المعادلة 22) يربط $R_{\text{closure}}$ بـ $R_{200}$، وكسر الباريونات في الهالة ($f_{b, \text{halo}}$)، وكسر الغاز ($f_{\text{gas}}$)، ومؤشر التراتبية ($n$ أو الميل $\eta$)، والمعلمات الكونية ($\Omega_m, \Omega_\Lambda, z$).

التحقق العددي

لاختبار النموذج، استخدم المؤلفون مجموعة من محاكاة EAGLE الهيدروديناميكية الكونية (Barnes et al. 2018) مع تغيير الثوابت الكونية:

• الكونيات (Cosmologies): $\Lambda = 0$ (أينشتاين-دي سيتر)، $\Lambda_0$ (القيمة المرصودة)، $10\Lambda_0$، $30\Lambda_0$، و $100\Lambda_0$.
• التغذية الراجعة: استخدمت جميع عمليات المحاكاة فيزياء فرعية متطابقة لتكوين النجوم والتغذية الراجعة (النجمية وAGN)، تمت معايرتها لتطابق الكون المرصود، وذلك لعزل تأثير تغيير $\Lambda$.
• التحليل: قاموا بحساب نصف قطر الإغلاق للهالات ذات الكتلة $M_{200} > 10^{11} M_\odot$ عبر انزياحات حمراء $0 < z < 3$. كما أجروا عملية تقليل $\chi^2$ لملاءمة معلمات النموذج التحليلي ($\eta$ و $f_{\text{gas}}$ مع بيانات المحاكاة).

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

1. الاشتقاق من المبادئ الأولى: تقدم الورقة أول اشتقاق تحليلي لنصف قطر الإغلاق بناءً على توازن الجاذبية والضغط والطاقة المظلمة، متجاوزةً مجرد الملاءمات التجريبية.
2. الاعتماد الكوني: تحدد الورقة بشكل منهجي كيف يؤثر الثابت الكوني ($\Lambda$) على إخلاء الباريونات، وهو عامل تم تجاهله سابقاً في دراسات نصف قطر الإغلاق.
3. التفسير الفيزيائي للملاءمات التجريبية: يوضح المؤلفون أن العلاقة الخطية التجريبية التي اقترحها Ayromlou et al. (2023) هي حالة حدية لنموذجهم الأكثر عمومية والمستند إلى المبادئ الأولى.
4. تحديد دور الطاقة المظلمة: تعزل الدراسة المساهمة المحددة للطاقة المظلمة في فقدان الباريونات، مفرقةً بينها وبين التغذية الراجعة الفيزيائية الفلكية.

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

• دقة النموذج: يظهر النموذج التحليلي توافقاً ممتازاً مع بيانات محاكاة EAGLE للهالات ذات $M_{200} > 10^{11} M_\odot$ في جميع الأكوان المختبرة ($0 \le \Lambda \le 30\Lambda_0$) والانزياحات الحمراء ($0 < z < 3$).
• تطور المعلمات:
  • يظهر ميل كثافة الغاز ($\eta$) وكسر الغاز ($f_{\text{gas}}$) تطوراً مع الانزياح الأحمر.
  • تأثير $\Lambda$: يؤدي زيادة $\Lambda$ إلى انخفاض في أفضل قيمة لـ $\eta$، مما يشير إلى انخفاض أكثر تدريجياً في كثافة الغاز مع نصف القطر (الغاز أكثر انتشاراً). كما يؤدي أيضاً إلى انخفاض في $f_{\text{gas}}$، مما يعني أن الهالات تحتفظ بغاز أقل في مناطقها الداخلية.
• التقريب الخطي: تعطي التقريب الخطي للنموذج (الصالح عند ثابت الانزياح الأحمر) معلمات ($\eta_{\text{approx}} \approx 1.8$، $f_{\text{gas, approx}} \approx 0.002$) ثابتة بشكل ملحوظ عبر الانزياح الأحمر والكونيات. وهذا يشير إلى أن المحرك الأساسي لتغيرات نصف قطر الإغلاق هو التغير في المعلمات الكونية ($\Omega_m, \Omega_\Lambda$) وليس التغيرات في الخصائص الديناميكية الحرارية الذاتية للغاز.
• تأثير الطاقة المظلمة:
  • في الكون المرجعي ($\Lambda_0$)، تساهم الطاقة المظلمة في إخلاء الباريونات بنسبة تقارب 30% بالنسبة لتأثيرات التغذية الراجعة الفيزيائية الفلكية وضغط الغاز (بالمقارنة بين $\Lambda_0$ وكون أينشتاين-دي سيتر EdS).
  • في الأكوان ذات $\Lambda$ العالية ($\Lambda \ge 10\Lambda_0$)، تصبح الطاقة المظلمة هي العامل المهيمن في إخلاء الباريونات. يمكن أن يمتد نصف قطر الإغلاق إلى عشرات أنصاف أقطار فيريال، ويمنع التوسع المتسارع الباريونات من إعادة التراكم.
• المقارنة مع Ayromlou et al. (2023): بينما تعمل صيغة A23 بشكل جيد لقيم $\Lambda = \Lambda_0$، إلا أنها تنحرف بشكل كبير لقيم $\Lambda$ الأعلى. نموذج المؤلفين يعمم صيغة A23، ويوفر تفسيراً مدفوعاً بالفيزياء لهذا الارتباط التجريبي.

5. الأهمية والآثار المترتبة

• المتانة النظرية: يثبت النموذج أن العلاقة بين نصف قطر الإغلاق وكتلة الهالة متينة ضد التغيرات في الثابت الكوني، بشرما بقيت فيزياء التغذية الراجعة الفيزيائية الفلكية الأساسية متسقة.
• الاستدلال الأنثروبي (المرتبط بوجود الحياة): للنتائج آثار على الحجج الأنثروبية. ففي الأكوان ذات $\Lambda$ العالية، من شأن التجريد الفعال للباريونات من الهالات أن يثبط بشدة تكوين النجوم وتكوين العناصر الثقيلة، مما قد يجعل مثل هذه الأكوان غير صالحة للحياة. وتحدد الدراسة "نافذة الصلاحية للسكن" فيما يتعلق بالاحتفاظ بالباريونات.
• التوجيه الرصدي: يوفر النموذج صيغاً تنبؤية لنصف قطر الإغلاق كدالة لكتلة الهالة، والانزياح الأحمر، والكونيات. وهذا يقدم إطاراً لتفسير الملاحظات المستقبلية من أدوات مثل JWST و ELT، والتي تهدف إلى رسم خريطة لتوزيع الباريونات في الوسط المحيط بالمجرات والوسط بين المجرات.
• القيود: تشير الدراسة إلى أن عدم اليقين في ملاءمة المعلمات يزدะد في محاكاة $\Lambda$ العالية بسبب ندرة الهالات الضخمة في صناديق المحاكاة الأصغر. ويتطلب العمل المستقبلي أحجاماً أكبر وأعداداً أكبر من الجسيمات لتنقيح هذه القيود.

باختأصار، نجحت هذه الورقة في سد الفجوة بين المحاكاة العددية والنظرية التحليلية، مظهرةً أن الطاقة المظلمة تلعب دوراً غير مهمل، وفي حالات قصوى، مهيمناً في مشكلة "الباريونات المفقودة" من خلال تسريع إخلاء الغاز من الهالات المجرية.