Gravothermal Collapse: Robust Against Baryonic Feedback

ملخص تقني: الانهيار الحراري الجاذبي في هالات المادة المظلمة ذات التفاعل الذاتي (SIDM) تحت تأثير التغذية الراجعة الباريونية

بيان المشكلة
يتنبأ نموذج المادة المظلمًا الباردة (CDM) القياسي بوجود نتوءات كثافة مركزية شديدة الانحدار في هالات المادة المظلمة، وهو ما يتعارض غالبًا مع النوى الضحلة المرصودة في المجرات القزمة. وبينما يمكن لآليات التغذية الراجعة الباريونية (مثل التدفقات الخارجة الناجمة عن المستعرات العظمى) ضخ الطاقة لتحويل النتوءات إلى نوى، إلا أن الأرصاد الأخيرة تكشف عن تنوع واسع في توزيعات المادة المظلمة — يتراوح من النوى الضحلة إلى النتوءات شديدة الكثافة — وهو ما يكافح نموذج CDM لتفسيره. علاوة على ذلك، حددت عمليات العدسة الجاذبية القوية وأرصاد التيارات النجمية وجود مضطربات مدمجة ذات كثافة أعلى بكثير من توقعات نموذج CDM التقليدية. وفي إطار المادة المظلمة ذات التفاعل الذاتي (SIDM)، تؤدي التفاعلات الذاتية للمادة المظلمة إلى تطور حراري جاذبي، مما ينتج بشكل طبيعي كلاً من النوى منخفضة الكثافة وحالات الانهيار عالية الكثافة. ومع ذلك، اعتمدت الدراسات السابقة إلى حد كبير على محاكاة تعتمد على SIDM فقط أو جهود باريونية ثابتة، وبما أن الجهود الباريونية الواقعية التي تشكلها التغذية الراجعة هي ديناميكية وغير ثابتة (non-adiabatic)، فمن غير الواضح ما إذا كان التنبؤ الجوهري لنموذج SIDM — وهو الانهيار الحراري الجاذبي — سيظل صامدًا أمام حقن الطاقة والتقلبات المحتملة الناتجة عن التغذية الراجعة الباريونية القوية.

المنهجية
يستخدم المؤلفون محاكاة $N$-body محكومة باستخدام كود GADGET-2 مع وحدة SIDM لاختبار مدى صمود الانهيار الحراري الجاذبي تحت تأثير تغذية راجعة باريونية قوية.

• الظروف الأولية: تم إجراء محاكاة لهالتي مادة مظلمة بكتلة $M_{200} = 2 \times 10^9 M_\odot$: هالة ذات تركيز عالٍ ($c_{200} = 42$، تمثل تقلبًا بمقدار $3\sigma$)، وهالة ذات تركيز متوسط ($c_{200} = 15$، الوسيط الكوني). تتبع كلتا الهالتين ملفات تعريف نافارو-فرينك-وايت (NFW).
• معلمات SIDM: تم تثبيت مقطع التفاعل الذاتي لكل وحدة كتلة عند $\sigma/m = 50 \text{ cm}^2/\text{g}$، وهو ما يتوافق مع النماذج التي تفسر المضطربات عالية الكثافة في العدسات الجاذبية القوية ومجرة "Crater II" المنتشرة.
• نموذج التغذية الراجعة: تم تنفيذ جهد متذبذب شبه تحليلي لمحاكاة التغذية الراجعة الباريونية. تمت إضافة جهد "بلومر" مركزي بكتلة متغيرة زمنيًا $M_{ext}(t)$ إلى حساب القوة. يتضمن النموذج "تغذية راجعة قوية" (حيث تتذبذب $M_{ext}$ بين قيم موجبة وسالبة، مما يخلق تسارعًا تنافريًا) و"تغذية راجعة ضعيفة" (تتذبذب بين الصفر والقيم الموجبة). تم ضبط فترة التذبذب عند $P = 0.2$ مليار سنة، مع سعة كتلة قصوى $M_{max} = 10^7 M_\odot$ ونصف قطر مقياس $a = 0.1$ كيلو فرسخ فلكي (kpc).
• السيناريوهات: تقارن الدراسة تطور SIDM-only مقابل حالات التغذية الراجعة القوية/الضعيفة المستمرة، والتغذية الراجعة الدورية (النشطة لفترات محددة، مثل 0–4 مليار سنة أو 8–10 مليار سنة) متبوعة إما بالإزالة أو الاستبدال بجهد ثابت.

النتائج الرئيسية
تكشف المحاكاة أن صمود الانهيار الحراري الجاذبي يعتمد بشكل حاسم على المنافسة بين مقياس زمن التوازن الحراري لـ SIDM ($t_{th}$) وفترة تذبذب التغذية الراجعة ($P$).

1. الهالات عالية التركيز ($c_{200} = 42$):

   • مقياس زمن التوازن الحراري قصير ($t_{th} \sim 0.01\text{--}0.04$ مليار سنة)، وهو أقصر بكثير من فترة التغذية الراجعة ($P = 0.2$ مليار سنة).
   • بناءً على ذلك، تهيمن ديناميكيات SIDM على التغذية الراجعة الباريونية.
   • تحت تأثير التغذية الراجعة القوية، يتأخر الانهيار بشكل طفيف فقط (بمقدار ~1 مليار سنة) ولا يتوقف أبدًا. تصل الكثافة المركزية في النهاية إلى مستويات تتوافق مع محاكاة SIDM-only والقيود الرصدية للمضطربات الكثيفة (مثل JVAS B1938+666).
   • التغذية الراجعة الضعيفة لها تأثير ضئيل أو تسرع الانهيار قليلاً بسبب الانكماش أثناء مراحل التدفق الداخل.

2. الهالات متوسطة التركيز ($c_{200} = 15$):

   • مقياس زمن التوازن الحراري أطول ($t_{th} > P$ لمعظم فترة التطور)، مما يجعل الهالة عرضة للتسخين الناتج عن التغذية الراجعة.
   • التغذية الراجعة القوية: خلال مرحلة توسع النواة، يمكن للتغذية الراجعة القوية أن تولد بسرعة نوى كبيرة (تصل إلى ~0.8 كيلو فرسخ فلكي) وتقلل الكثافة المركزية إلى ~25% من قيمتها الأولية. ومع ذلك، بمجرد توقف التغذية الراجعة، يستأنف الانهيار الحراري الجاذبي مساره، وإن كان بمعدل أبطأ.
   • التغذية الراجعة الدورية: يعتمد ملف الكثافة النهائي بشكل كبير على تاريخ التغذية الراجعة.
     • إذا حدثت التغذية الراجعة مبكرًا (0–4 مليار سنة) وتمت إزالتها، ينهار الهالة في وقت لاحق مقارنة بحالة SIDM-only.
     • إذا تم استبدال الجهد المتذبذب بجهد ثابت بعد التغذية الراجعة، يتسارع الانهيار.
     • إذا حدثت التغذية الراجعة متأخرًا (8–10 مليار سنة) خلال مرحلة الانهيار، فإنها تعطل العملية مؤقتًا لكنها لا تعكسها.
   • ينتج هذا الحساسية تنوعًا واسعًا في كثافات المركز النهائية، تتراوح من النوى الضحلة إلى النتوءات الكثيفة، اعتمادًا على الجدول الزمني المحدد للتغذية الراجعة.

الأهمية والادعاءات
تثبت الورقة أن الانهيار الحراري الجاذبي هو تنبؤ قوي لنموذج SIDM، ويظل صامدًا حتى أمام التغذية الراجعة الباريونية المتذبذبة وشديدة القوة.

• بالنسبة للهالات عالية التركيز: يضمن مقياس زمن التوازن الحراري القصير استمرار الانهيار بغض النظر عن شدة التغذية الراجعة. وهذا يعزز تفسير المضطربات الكثيفة في العدسات الجاذبية القوية والأجسام المدمجة كـ هالات SIDM منهارة النواة. كما يدعم سيناريوهات SIDM لإنبات الثقوب السوداء فائقة الكتلة في الكون المبكر، حتى في البيئات ذات تكوين النجوم المتفجر.
• بالنسبة للهالات متوسطة التركيز: يمكن للتفاعل بين التغذية الراجعة والتفاعلات الذاتية أن يفسر التنوع الملحوظ في توزيعات المادة المظلمة (من النوى الضحلة إلى النتوءات الكثيفة) داخل الهالات ذات الكتلة المتشابهة. الحالة النهائية لا تتحدد فقط بكتلة الهالة أو التركيز الأولي، بل يتم تعديلها بواسطة التاريخ الدوري للتغذية الراجعة الباريونية.
• الاستنتاج: تشير النتائج إلى أن الانهيار الحراري الجاذبي يظل "علامة فارقة" (smoking-gun) لنموذج SIDM. إن التنوع المتوقع في ملفات الكثافة، المدفوع بالجمع بين التفاعلات الذاتية وتاريخ التغذية الراجعة، يوفر إطارًا قابلًا للاختبار للأرصاد المستقبلية لمنحنيات دوران المجرات وطبيعة المادة المظلمة.