Galaxy assembly bias in cosmological hydrodynamical simulations: a comparison between SIMBA and IllustrisTNG

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

يُعد التوزيع المكاني للمجرات ملاحظة أساسية في علم الكونيات، حيث يتتبع البنية واسعة النطاق (LSS) للمادة المظلمة الكامنة. تفترض النماذج القياسية، مثل توزيع إشغال الهالات (HOD)، أن عدد المجرات في هالة المادة المظلمة يعتمد فقط على كتلة الهالة. ومع ذلك، فإن انحياز تجميع المجرات (GAB) يتحدى هذا الافتراض، مفترضاً أن تكتل المجرات يعتمد على خصائص ثانوية للهالة (مثل زمن التكوين، والتركيز) وعوامل بيئية (مثل الكثافة المحلية، والمجالات المدية) تتجاوز مجرد الكتلة.

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

1. قياس قوة وتطور GAB عبر الزمن الأحمر في محاكاتين رئيسيتين: SIMBA و IllustrisTNG (TNG).
2. تحديد مدى كون GAB مدفوعاً بالكثافة الزائدة للمادة الملساء ($\delta_5$) مقابل تباين المد والجزر ($\alpha_5$).
3. استقصاء دور الهالات منخفضة الكتلة مقابل عالية الكتلة في توليد إشارة GAB.
4. تقديم بارامتر (معلمة) لدمج GAB المعتمد على الكثافة في نماذج HOD.

2. المنهجية

يستخدم المؤلفون أكبر النسخ المتاحة من محاكات هيدروديناميكية كونية:

• SIMBA (m100n1024): يستخدم كود GIZMO مع هيدروديناميكا الكتلة النهائية غير المتشكلة (MFM). ويتميز بتنفيذات محددة للرياح المجرية وتغذية AGN الراجعة (النفثات ثنائية القطب).
• IllustrisTNG (TNG300-1): يستخدم كود AREPO ذو الشبكة المتحركة. ويتميز بالتغذية الراجعة النجمية الحركية وتغذية AGN متعددة الأنماط (حرارية، حركية، إشعاعية).

الخطوات المنهجية الرئيسية:

• اختيار المجرات: يتم اختيار المجرات حسب الكتلة النجمية لمطابقة كثافة العدد بين المحاكاتين ($N \approx 0.043 \, h^3 \text{Mpc}^{-3}$ عند $z=0$).
• تقنية الخلط (Shuffling): لعزل GAB، يستخدم المؤلفون طريقة "الخلط". حيث يتم إعادة تعيين المجرات عشوائياً إلى هالات من نفس فئة الكتلة. وهذا يزيل أي اعتماد لتعشيش المجرات على خصائص الهالة بخلاف الكتلة.
  • إجمالي GAB: يُقاس كنسبة دالة الارتباط المسقطة $w_p(r_p)$ للعينة الأصلية إلى العينة المخلطة: $[w_p(r_p)/w_p(r_p)_{\text{shuffled}} - 1]$.
• الخلط المزدوج: لعزل الدوافع البيئية، يتم خلط المجرات ضمن فئات من كتلة الهالة و متغير بيئي محدد:
  • الكثافة الزائدة الملساء ($\delta_5$): كثافة المادة الملساء على مدى $5 \, h^{-1}\text{Mpc}$.
  • تباين المد والجزر ($\alpha_5$): مقياس لشكل مجال المد والجزر، والذي تم حسابه ليكون غير مرتبط بـ $\delta_5$.
• تحليل HOD: يتم قياس توزيع إشغال الهالة للهالات في الربعين العلوي والسفلي من $\delta_5$ و $\alpha_5$ لمعرفة كيف تعدل البيئة من تعشيش المجرات المركزية والتابعة.
• الخلط المعتمد على الكتلة: يتم تقييد الخلط للهالات منخفضة الكتلة ($M_h < 10^{11.5} M_\odot$) والهالات عالية الكتلة بشكل منفصل لتحديد نظام الكتلة المسؤول عن GAB.

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

• المقارنة بين المحاكات: أول مقارنة مباشرة لتطور GAB بين SIMBA و TNG، مما يسلط الض الضوء على كيفية تغيير الفيزياء البارونية المختلفة (تحديداً تغذية AGN الراجعة) لاتجاهات GAB.
• التفكيك البيئي: التحديد الكمي للمساهمات النسبية للكثافة الزائدة المحلية ($\delta_5$) وتباين المد والجزر ($\alpha_5$) في إشارة GAB الإجمالية.
• هيمنة الكتلة المنخفضة: إثبات أن GAB مدفوع بشكل حصري تقريباً بواسطة الهالات منخفضة الكتلة، وخاصة تلك القريبة من كتلة "بدء التشغيل" لتكوين المجرات.
• بارامترية HOD: تطوير صيغة ملائمة تعدل معاملات HOD القياسية (تحديداً ميول "بدء التشغيل" للمجرات المركزية والتابعة) كدالة خطية لرتبة البيئة، مما يسمح بإنشاء كتالوجات وهمية تشمل GAB.

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

أ. تطور GAB عبر الزمن الأحمر

• الزمن الأحمر العالي ($z > 2$): تُظهر كلتا المحاكات GAB ضئيلاً عند $z=5$، يرتفع إلى حوالي $5\%$ عند $z=2$.
• التباعد عند الزمن الأحمر المنخفض ($z < 2$):
  • SIMBA: ينخفض GAB بعد $z=2$، ليصل إلى 0% تقريباً عند $z=0$.
  • TNG: يستمر Gable في الارتفاع، ليصل إلى $\sim 10\%$ عند $z=0$.
• التفسير: يتزامن التباعد مع بدء تغذية AGN الراجعة الكبيرة. تعمل تغذية النفثات ثنائية القطب في SIMBA بكفاءة على كبح تكوين النجوم في الهالات ذات مقياس المجموعات وتفريغ الغاز، مما يقلل من الاعتماد البيئي لتكوين المجرات عند الزمن الأحمر المنخفض. أما في TNG، فإن اتجاه النفثات العشوائي وتوصيل التغذية الراجعة المختلف يؤديان إلى استمرار أقوى للاعتماد البيئي.

ب. الدوافع البيئية ($\delta_5$ مقابل $\alpha_5$)

• الكثافة الزائدة ($\delta_5$): تمثل الجزء الأكبر من إشارة GAB في كلتا المحاكات.
• تباين المد والجزر ($\alpha_5$):
  • SIMBA: مساهمته ضئيلة.
  • TNG: يساهم بجزء صغير ولكنه غير مهمل، وإن كان لا يزال ثانوياً مقارنة بـ $\delta_5$.
• الاستنتاج: المحرك الأساسي لـ GAB هو الكثافة الزائدة للمادة المحلية، وليس شكل مجال المد والجزر.

ج. HOD ونظام الكتلة

• تباين HOD: في TNG، تزيد البيئات ذات $\delta_5$ العالية بشكل كبير من احتمالية استضافة كل من المجرات المركزية والتابعة في الهالات منخفضة الكتلة. في SIMBA، التأثير يتركز بشكل أساسي على المجرات المركزية.
• الهالات منخفضة الكتلة: استعادة كامل إشارة Gable تقريباً عند تقييد الخلط للهالات منخفضة الكتلة ($M_h < 10^{11.5} M_\odot$). تساهم الهالات عالية الكتلة بشكل مهمل.
• الآلية: الهالات منخفضة الكتلة في المناطق عالية الكثافة هي أكثر عرضة لاستضافة المجرات. أما الهالات عالية الكتلة فهي دائماً تقريباً في مناطق كثيفة، مما يترك تشتتاً ضئيلاً في بيئتها، وبالتالي لا توجد إشارة GAB.
• عتبة الكتلة النجمية: رفع عتبة الكتلة النجمية في TNG (إلى $10^{10} M_\odot$) يقلل قليلاً من سعة GAB الإجمالية ولكنه يحافظ على الاتجاه بأن $\delta_5$ يقود الإشارة عبر الهالات منخفضة الكتلة.

د. بارامترية HOD

يقترح المؤلفون نموذج HOD معدلاً حيث تعتمد الميول ($\mu$ و $\nu$) التي تحكم "بدء تشغيل" الجمهرات المركزية والتابعة خطياً على رتبة الكثافة المحلية ($r_\delta$):
$$\mu = \mu_{\text{tot}} - C_\delta (r_\delta - 0.5)$$
هذا يسمح للكتالوجات الوهمية بإعادة إنتاج التكتل المعتمد على الكثافة كما هو مشهد في المحاكات الهيدروديناميكية.

5. الأهمية

• علم الكونيات الدقيق: يمكن لـ GAB أن يغير سعة تكتل المجرات بنسبة تصل إلى 10%. تجاهل هذا التأثير في المسوح المستقبلية (مثل DESI، Euclid) قد يؤدي إلى انحياز في تقدير المعلمات الكونية.
• الحساسية للفيزياء البارونية: يوضح الفرق الصارخ في تطور GAB بين SIMBA و TNG أن GAB حساس للغاية للتنفيذ المحدد للفيزياء البارونية (خاصة تغذية AGN الراجعة). هذا يشير إلى أن GAB يمكن أن يعمل كمسبار لتقييد نماذج التغذية الراجعة.
• تحسين النمذجة: توفر البارامترية المقدمة أداة عملية لتوليد الكتالوجات الوهمية القائمة على HOD، متجاوزة افتراض أن كتلة الهالة هي المحدد الوحيد لتعشيش المجرات.
• التوتر الرصدي: تتناقض نتيجة الدور الضئيل لتباين المد والجزر في المحاكات مع بعض الادعاءات الرصدية الأخيرة، مما يسلط الضوء على الحاجة لمزيد من الاستقصاء في الطبيعة التجريبية لانحياز التجميع.

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