بي ثوري – محاكاة المجرة v2 – الجيل الأولي 2025 مايو 17 مع كلود

الكتلة الخفية لدرب التبانة نظرية يوكاوا ثلاثية الأبعاد مع اقتطاع القرص الفيزيائي

المحاكاة المصححة: تقع سرعة القرص الباريوني في كبلري خارج الحافة الفيزيائية، وتملأ نواة يوكاوا ثلاثية الأبعاد من BeeTheory كل الفضاء. بارامتران، بيانات دوران عصر غايا ونموذج قرص مبتور.

موقع BeeTheory.com – Ou et al., MNRAS 528, 2024 – تصحيح BeeTheory v2

K = 0.040 كيلو بكسل-¹

اقتران الموجات

α = 0.087 kpc-¹

التماسك العكسي

ℓ = 11.5 كيلو بكسل

طول التماسك

χ²/دوف ≈ 0.31

ملاءمة مبسطة ممتازة

0. النتيجة – المعادلات والمعلمات

تولد كل حلقة حلقية من قرص المجرة عند نصف القطر R′ مجال كتلة مظلمة فعالة ثلاثية الأبعاد من خلال نواة يوكاوا ذات نظرية النحلة. تساوي الكثافة المظلمة الكلية عند نصف القطر الكروي r:

\(\rho_{\mathrm{dark}}(r)=K\int_0^{R_{\mathrm{max}}}\Sigma_0e^{-R’/R_d}\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\,2\pi R’\,dR’\) \(D(r,R’)=\sqrt{r^2+R’^2}\)

تُشتق النواة من قانون قوة نظرية النحل المصحح:

\(F(D)\propto\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\)

تُختزل إلى الصورة المربعة العكسية النيوتونية ل D أصغر بكثير من طول التماسك ℓ.

\(D\ll\ll\ell=\frac{1}{\alpha}\quad\Longrightarrow\quad F(D)\propto\frac{1}{D^2}\)

وتستخدم سرعة القرص الباريوني معادلة فريمان داخل حافته الفيزيائية _Rtrunc ≈ _4Rd = 10.4 كيلو بك، ثم تنتقل بسلاسة إلى السقوط الكيبليري المتوقع من توزيع كتلة محدودة.

\(K=0.0397\,\mathrm{kpc}^{-1},\qquad \alpha=0.0868\,\mathrm{kpc}^{-1},\qquad \ell=\frac{1}{\alpha}=11.5\,\mathrm{kpc}\)

ملخص الملاءمة

قابل للملاحظة	قيمة عصر جايا	نظرية النحل	اسحب
_Vc(4 kpc)	220 ± 10 كم/ثانية	219.8 كم/ثانية	-0.02σ
_Vc(8 kpc)	230 ± 6 كم/ثانية	233.2 كم/ثانية	+0.53σ
_Vc(12 كيلو بكسل)	226 ± 7 كم/ثانية	223.8 كم/ثانية	-0.31σ
_Vc(20 كيلو بكسل)	215 ± 10 كم/ثانية	211.2 كم/ثانية	-0.38σ
_Vc(27.3 كيلو بكسل)	173 ± 17 كم/ثانية	199.0 كم/ثانية	+1.53σ
_ρdark(R⊙ = 8 كيلو بكسل)	0.39 ± 0.03 جيجا فولت/سم مكعب	0.47 جيجا فولت/سم مكعب	+2.3σ
_مدارك(<8 كيلو بكسل)	~5 × 10¹⁰ M⊙	5.3 × 10¹⁰ M⊙	إغلاق
_متوت(<200 كيلو بكسل)	5-9 × 10¹¹ M⊙	3.3 × 10¹¹ M⊙	نهاية منخفضة

تعطي الملاءمة المبسطة χ²/دوف ≈ 0.31. وتبقى النقطة الأصعب هي القيمة الخارجية لعصر غايا عند 27.3 كيلو بك، حيث يكون الانخفاض المرصود أكثر حدة مما يتوقعه هذا النموذج ثنائي المعامل.

1. اقتطاع القرص – لماذا وكيف؟

1.1 مشكلة القرص الأسي اللامتناهي

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

\(\Sigma(R)=\Sigma_0e^{-R/R_d}\)

بعد حافة القرص، تميل السرعة الباريونية إلى ما وراء حافة القرص:

\(V_{\mathrm{bar}}(R)\xrightarrow{R\gg R_d}\sqrt{\frac{GM_{\mathrm{bar,tot}}}{R}}\) \(M_{\mathrm{bar,tot}}=M_{\mathrm{disk}}+M_{\mathrm{bulge}}\approx4.7\times10^{10}M_\odot\)

أمثلة على القيم هي:

\(V_{\mathrm{bar}}(30\,\mathrm{kpc})\approx82\,\mathrm{km/s},\qquad V_{\mathrm{bar}}(50\,\mathrm{kpc})\approx63\,\mathrm{km/s}\)

1.2 معادلة الاقتطاع السلس

تستخدم المحاكاة انتقالًا سلسًا بين صيغة قرص فريمان والقيمة الكبلرية. يتمركز الانتقال عند _Rtrunc = _4Rd = 10.4 كيلو بكسل مع عرض σ = 1.5 كيلو بكسل.

\(V_{\mathrm{bar}}(R)=\sqrt{(1-w)V_{\mathrm{Freeman}}^2(R)+w\,\min(V_{\mathrm{Freeman}},V_{\mathrm{Kepler}})^2}\) \(w(R)=\frac{1}{2}\left[1+\tanh\left(\frac{R-R_{\mathrm{trunc}}}{\sigma}\right)\right]\) \(R_{\mathrm{trunc}}=4R_d=10.4\,\mathrm{kpc},\qquad \sigma=1.5\,\mathrm{kpc}\)

تمنع الدالة الصغرى القرص الباريوني من تجاوز الحد الكيبليري الفيزيائي خارج حافة القرص.

R	_فريمان	_{فيكيبليريان}	_{Vbar، مقطوع}	النظام المهيمن
5 كيلو متر مكعب	174.5 كم/ثانية	201.1 كم/ثانية	174.5 كم/ثانية	فريمان
8 كيلو متر مكعب	161.5 كم/ثانية	159.0 كم/ثانية	161.5 كم/ثانية	فريمان ≈ كبلر
10.4 كيلو متر مكعب	143.0 كم/ثانية	139.3 كم/ثانية	141.2 كم/ثانية	المرحلة الانتقالية
16 كيلو متر مكعب	112.4 كم/ثانية	112.4 كم/ثانية	112.4 كم/ثانية	كيبليريان
25 كيلو متر مكعب	89.9 كم/ثانية	89.9 كم/ثانية	89.9 كم/ثانية	كيبليريان
50 كيلو متر مكعب	63.6 كم/ثانية	63.6 كم/ثانية	63.6 كم/ثانية	كيبليريان

2. كثافة الكتلة المظلمة ثلاثية الأبعاد لنظرية النحلة

2.1 حلقات الأقراص المشعة ثلاثية الأبعاد

كل حلقة من قرص المجرة عند نصف القطر R′ بعرض dR′ لها كتلة:

\(dM=\Sigma(R’)\,2\pi R’\,dR’\)

في نظرية النحلة، تُولِّد هذه الحلقة مجال موجة جاذبية ينتشر في الأبعاد المكانية الثلاثة. في تقريب أحادي القطب، تكون المسافة إلى نقطة مجال ثلاثي الأبعاد عند نصف القطر الكروي r

\(D(r,R’)=\sqrt{r^2+R’^2}\)

الصورة العددية للكثافة المظلمة هي:

\(\rho_{\mathrm{dark}}(r)=K\sum_{i=1}^{N}\Sigma_0e^{-R’_i/R_d}\frac{(1+\alpha D_i)e^{-\alpha D_i}}{D_i^2}\,2\pi R’_i\Delta R’\) \(D_i=\sqrt{r^2+R_i’^2},\qquad R’_i=\left(i-\frac{1}{2}\right)\frac{R_{\mathrm{max}}}{N}\) \(N=60,\qquad R_{\mathrm{max}}=25\,\mathrm{kpc}\)

2.2 الكتلة المظلمة المغلقة والسرعة الدائرية

\(M_{\mathrm{dark}}(<r)=\int_0^r4\pi s^2\rho_{\mathrm{dark}}(s)\,ds\) \(M_{\mathrm{dark}}(<r)\approx\sum_{j=1}^{30}4\pi r_j^2\rho_{\mathrm{dark}}(r_j)\Delta r\) \(V_{\mathrm{dark}}(R)=\sqrt{\frac{GM_{\mathrm{dark}}(<R)}{R}}\) \(V_c(R)=\sqrt{V_{\mathrm{bar}}^2(R)+V_{\mathrm{dark}}^2(R)}\)

2.3 السلوك التقاربي

\(\rho_{\mathrm{dark}}(r)\approx\frac{2\pi K\Sigma_0R_d^2}{r^2}\left(1+\alpha r+\frac{\alpha^2r^2}{2}\right)e^{-\alpha r}\)

بالنسبة إلى αr ≪ 1:

\(\rho_{\mathrm{dark}}(r)\xrightarrow{\alpha r\ll1}\frac{2\pi K\Sigma_0R_d^2}{r^2}\) \(M_{\\mathrm{dark}}(r)\r)\propto r\qqad\Longrightarrow\qqquad V_{\mathrm{dark}}\approx\mathrm{constant}\)

3. نتائج المحاكاة – الرسوم البيانية التفاعلية

تحتفظ المحاكاة أدناه بالنموذج العددي وأشرطة التمرير ومنحنى الدوران وملف تعريف الكتلة وملف تعريف الكثافة وتحديث χ²² المباشر. الصق هذه الصفحة في WordPress مع تمكين تنفيذ البرنامج النصي.

منحنى الدوران – بيانات حقبة غايا مقابل نظرية النحل مع القرص المقطوع

باريونات فقط، مبتورة إجمالي نظرية النحلة المكوّن المظلم بيانات حقبة غايا

مستكشف المعلمات – ضبط K، α، و _Rtrunc

K (kpc-¹) 0.040

α (kpc-¹) 0.087

_Rtrunc (kpc) 10.4

χ²/دوف – | ℓ: – kpc | ρ(R⊙): – جيجا فولت/سم مكعب

المظهر الجانبي للكتلة: القرص المرئي مقابل الكتلة المظلمة ثلاثية الأبعاد مقابل الكتلة المظلمة الكلية

القرص المرئي + الانتفاخ الكتلة المظلمة النظرية الكتلة الكلية

ص (كيلو متر مكعب)	_Mbar (10¹⁰ M⊙ M)	_مدارك (10¹⁰ M⊙ M⊙)	_متوت (10¹⁰ م ⊙ م)	DM/بار	_ρdark (جيجا فولت/سم مكعب)
تحميل…

ملف الكثافة المظلمة _ρdark(r) – مقياس اللوغاريتم

نظرية النحل مرجع r-² المتساوي الحرارة مرجع NFW

4. التفسير المادي والعالمية

4.1 طول التماسك 4.1

داخل طول التماسك، تتصرف نواة يوكاوا تقريبًا مثل النواة النيوتونية 1/D². تتبع الكثافة المظلمة r-² تقريبًا ويكون منحنى الدوران مسطحًا. بعد ℓ، ينتج عن الكبت الأسي الانخفاض الملحوظ في القرص الخارجي.

\(\ell=\frac{1}{\alpha}\approx11.5\,\mathrm{kpc}\) \(\frac{\ell}{R_d}=\frac{11.5}{2.6}\approx4.4\)

4.2 اقتران بلا أبعاد

يمكن تعريف اقتران نظرية النحل بلا أبعاد على النحو التالي:

\(\lambda_{\mathrm{galaxy}}=K\ell^2\) \(\lambda_{\mathrm{galaxy}}=0.040\times(11.5)^2\approx5.3\)

ويمكن مقارنة ذلك من حيث الحجم بالاقتران المستدل عليه من معايرة H₂، حيث يبلغ λ حوالي 3-4. وتبقى الشمولية المحتملة لهذا الرقم مسألة مركزية مفتوحة.

4.3 المقارنة مع النماذج القياسية

الطراز	المعلمات	الملاءمة النموذجية	المقياس	الآلية
هالة متساوية الحرارة	2	معتدل	نصف القطر الأساسي	المنحنى الظاهري المسطح
ملف تعريف NFW	2	قوي	_ص	ملف محاكاة الجسم N
إيناستو	2-3	قوي	_r-2	الملف الشخصي التجريبي المرن
نظرية النحلة 3D يوكاوا	2	واعدة	ℓ	اقتران الكتلة الموجية من القرص

تظل النقطة الخارجية لعصر غايا هي القيد الأصعب. يمكن إنتاج انخفاض أكثر حدة مع طول تماسك أصغر، لكن ذلك يزيد من سوء الملاءمة الداخلية. ستكون البيانات المستقبلية من Gaia DR4 والعناقيد الكروية والتدفقات النجمية اختبارات مهمة.

المراجع

Ou, X. et al. - ملف تعريف المادة المظلمة لمجرة درب التبانة المستدل عليه من منحنى السرعة الدائري، MNRAS 528, 2024.
دوتيرتر، إكس. - نظرية النحلة™: نمذجة الجاذبية على أساس الموجة، BeeTheory.com v2، 2023.
فريمان، ك. س. - عن أقراص المجرات الحلزونية ومجرات S0، ApJ 160, 811, 1970.
McMillan, P. J. - التوزيع الكتلي وإمكانات الجاذبية لدرب التبانة، MNRAS 465, 76, 2017.
Navarro, J. F., F. F., Frenk, C. S., White, S. D. M. - ملف تعريف الكثافة العالمية من التجميع الهرمي, ApJ 490, 1997.

موقع BeeTheory.com - الجاذبية الكمية القائمة على الموجات