TRINITY: A coupled model of winds, radiation, and photoionised gas in molecular clouds. I. Methods and validation

ملخص تقني: TRINITY: نموذج مقترن للرياح، والإشعاع، والغاز المتأين ضوئياً في السحب الجزيئية – الجزء الأول: المنهجية والتحقق

المشكلة والسياق
تتشتت السحب الجزيئية العملاقة (GMCs) خلال 1-5 مليون سنة من ظهور النجوم الضخمة، عادةً قبل وصول أول مستعرات عظمى ناتجة عن انهيار النوى. تشير المسوحات الرصدية (مثل PHANGS وLEGUS) إلى أن آليات التغذية الراجعة ما قبل المستعر الأعظم (stellar feedback)—وهي الرياح النجمية، وضغط الإشعاع، وضغط الغاز المتأين ضوئياً ($P_{\text{H II}}$)—هي المحركات الأساسية لتطهير السحب. ومع ذلك، لا يزال مدى هيمنة هذه القنوات النسبية، وتوقيتها، واعتمادها على خصائص السحابة (الكتلة، ملف الكثافة، كفاءة تكوين النجوم) أمراً غير مؤكد. غالباً ما تتعامل الأدوات شبه التحليلية الحالية مع آليات التغذية الراجعة هذه بمعزل عن بعضها البعض أو تعتمد على ملفات كثافة مبسطة، بينما تعد المحاكاة ثلاثية الأبعاد مكلفة حاسوبياً للمسوحات الإحصائية. هناك حاجة إلى نموذج أمامي (forward model) فعال ومتسق ذاتياً يربط جميع قنوات التغذية الراجعة ما قبل المستعر الأعظم بهياكل سحابية متنوعة لرسم خريطة هيمنة التغذية الراجعة عبر فضاء المعلمات.

المنهجية
يقدم المؤلفون Trinity، وهو كود ذو غلاف رقيق (thin-shell code) أحادي البعد ومتماثل كروياً، يقوم بتطوير فقاعة مدفوعة بالتغذية الراجعة وغلاف متراكم داخل سحابة جزيئية. يحل Trinity محل إطار عمل warpfield ولكنه يقدم تغييرين هيكليين حاسمين:

1. ضغط الغاز المتأين ضوئياً الصريح: يتم تضمين $P_{\text{H II}}$ صراحة في معادلة الحركة للغلاف، بدلاً من معاملته فقط كشرط حدودي أو تجاهله.
2. وصف مرحلي للقيادة (Phase-Aware Driving Prescription): يختار الكود ديناميكياً مصطلح ضغط القيادة بناءً على المرحلة التطورية لتجنب الازدواج في حساب القوى:
   • مرحلة القيادة بالدفع الطاقي (Energy-driven phase): يكون ضغط القيادة هو الحد الأقصى بين الضغط الحراري للفقاعة الساخنة ($P_b$) وضغط الغاز المتأين ضوئياً ($P_{\text{H II}}$).
   • مرحلة الانتقال (Transition phase): يكون ضغط القيادة هو الحد الأقصى بين $P_b$ المتبقي ومجموع $P_{\text{H II}}$ وضغط الاندفاع (ram pressure) للرياح المحقونة حديثاً ومخلفات المستعر الأعظم ($P_{\text{ram}}$).
   • مرحلة القيادة بالزخم (Momentum-driven phase): بمجرد تبريد الفقاعة الساخنة ($E_b \to 0$)، يكون ضغط القيادة هو مجموع $P_{\text{H II}} + P_{\text{ram}}$.

يقوم النموذج بحل نظام مقترن من المعادلات التفاضلية العادية لنصف قطر الفقاعة ($R_b$)، والسرعة ($v_b$)، والطاقة الحرارية ($E_b$). ويتضمن ما يلي:

• الشروط الابتدائية: يدعم توزيعات الكثافة المنتظمة، والأسية الجزئية (power-law)، وكروات بونور-إيبرت (Bonnor–Ebert).
• نماذج الفيزياء: يتضمن الرياح النجمية، المستعرات العظمى، ضغط الإشعاع المباشر والمُعاد معالجته بواسطة الغبار، الجاذبية الذاتية للغلاف والعنقود، والتبخر الموصلي من الغلاف إلى الفقاعة الساخنة.
• الاقترن الإشعاعي-الغلافي: تقوم عملية حساب هيكل الغلاف الهيدروستاتيكي بتحديد موقع جبهة التأين ($R_{\text{if}}$)، وأجزاء الإشعاع الممتصة، والعمق البصري للأشعة تحت الحمراء بشكل متسق ذاتياً.
• التبريد: يعالج الانتقال من حالة التوازن التأيني التصادمي (CIE) إلى الأنظمة غير التصادمية (non-CIE)، ويحسب خسائر التبريد لتحديد الانتقال من الطاقة إلى الزخم.

النتائج الرئيسية
قام المؤلفون بالتحقق من صحة الكود مقابل الحدود التحليلية (رياح Weaver، توسع Spitzer، ونموذج المكنسة/snowplow) وتشغيلات warpfield، ثم عرضوا النتائج لسحب تتراوح كتلتها بين $M_{\text{cloud}} = 10^5–10^{6.5} M_\odot$، وكثافات نوى $n_{\text{core}} = 10^3–10^4 \text{ cm}^{-3}$، وكفاءة تكوين نجوم $\epsilon = 0.01–0.30$.

1. تأثير الغاز المتأين ضوئياً: يؤدي تضمين $P_{\text{H II}}$ إلى تغيير الديناميكيات بشكل ملموس. في تجربة تحقق مرجعية ($M_{\text{cloud}}=10^6 M_\odot, \epsilon=0.01$)، أدى تضمين $P_{\text{H II}}$ إلى زيادة نصف قطر الغلاف بنسبة $\sim 17\%$ عند 10 ملايين سنة مقارنة بتشغيل warpfield المكافئ الذي تم فيه تعطيل $P_{\text{H II}}$. وتظل $P_{\text{H II}}$ مهمة ديناميكياً حتى بعد الانتقال إلى التطور المهيمن بالزخم.
2. الاعتماد على ملف الكثافة: يحدد هيكل السحابة النتيجة حتى عندما تكون التجمعات النجمية ثابتة. بالنسبة لـ M_{\text{cloud} و $n_{\text{core}}$ و $\epsilon$ ثابتة:
   • الملفات المتجانسة ($\rho \propto r^0$) والضحلة ($\rho \propto r^{-1}$) تعيد الانهيار.
   • الملف الحاد ($\rho \propto r^{-2}$) يسمح للغلاف بالتسارع عبر تدرج الكثافة والاستمرار في التوسع.
   • كروات بونور-إيبرت تعيد الانهيار في الحالات التي قد تتشتت فيها السحب المتجانسة، أو تتشتت في وقت لاحق بشكل ملحوظ.
3. الجداول الزمنية للتشتت: في شبكة من النماذج حيث يكون التشتت مضموناً، تتشتت سحب بونور-إيبرت متأخرة بنحو 55% عن السحب المتجانسة. وينشأ هذا الفرق المنهجي لأن هياكل بونور-إيبرت توزع نفس الكتلة على نصف قطر فيزيائي أكبر، مما يتطلب من الغلاف قطع مسافة أبعد.
4. الهيمنة ما قبل المستعر الأعظم: في مجال المعلمات المختبر، يتم تحديد تشتت السحابة أو إعادة انهيارها بواسطة التغذية الراجعة ما قبل المستعر الأعظم (الرياح، الإشعاع، $P_{\text{H II}}$) قبل وصول أول مستعر أعظم ناتج عن انهيار النواة (عادةً $>3$ مليون سنة).
5. هروب فوتونات لايمان-كونتينيووم (LyC): يتتبع الكود كسر هروب فوتونات لايمان-كونتينيووم ($f_{\text{LyC}}^{\text{esc}}$) كـ "ساعة" منفصلة عن التطهير الديناميكي، مما يميز الوقت الذي تغير فيه التغذية الراجعة ديناميكيات الغلاف عن الوقت الذي تبدأ فيه الفوتونات بالتسرب إلى الوسط بين النجمي المنتشر.

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

• كيفية اقتران قنوات التغذية الراجعة بالغلاف (تحديداً دور $P_{\text{H II}}$).
• توقيت انتقال تبريد الفقاعة الساخنة.
• البنية الكثافية للسحابة الأم.

ويؤكد البحث أن افتراض سحابة متجانسة نموذجية أو منحدر قوة واحدة (power-law) يؤدي إلى أخطاء منهجية في أوقات الظهور المستنتجة وعتبات التشتت. ومن خلال تقديم أداة تربط الرياح والإشعاع والتأين الضوئي بشكل متسق ذاتياً مع ملفات كثافة مرنة، يتيح Trinity دراسة مناطق H II المحللة ويوفر توصيفات فرعية للنماذج على مستوى المجرة فيما يتعلق بتسرب الفوتونات المؤينة وتشتت السحب.

القصور
يشير المؤلفون صراحةً إلى أن التماثل الكروي أحادي البعد يغفل الهياكل الزاوية (مثل تدفقات الشامبين أو انفجارات "البلستر") ولا يمكنه حل عدم الاستقرار مثل "كيلفن-هلمهولتز" أو "رايلي-تايلور" التي قد تعزز التبريد. كما أن معالجة $P_{\text{H II}}$ تستخدم نصف قطر واحداً لنظام الفقاعة-جبهة التأين بدلاً من تطوير $R_b$ و $R_{\text{if}}$ كمتغيرات ديناميكية منفصلة. بالإضافة إلى ذلك، يتجاهل النموذج حالياً الدعم المغناطيسي، وتدرجات ضغط الأشعة الكونية، وتطور الكسر المغطي (الجلود المتسربة). وقد تم تحديد هذه النقاط كأهداف للتطوير المستقبلي.