المحاكاة العددية لنموذج الكتلة المظلمة لنظرية النحلة

سرد كامل وقابل للتكرار للتكامل العددي، وتحلل السرعة الباريونية، وتصغير χ²، وخيارات التنفيذ وراء محاكاة الكتلة المظلمة لنظرية النحلة.

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

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

المحتويات

نظرة عامة على المحاكاة
بيانات الرصد
نموذج السرعة الباريونية
معادلة الكثافة المظلمة للنحل
مخطط التكامل العددي
χ² التقليل إلى الحد الأدنى وتركيب البارامترات
التقارب وميزانية الخطأ
الرمز المرجعي الكامل
كيفية إعادة إنتاج المحاكاة

ℓ ≈ 130 كيلو بكسل

أفضل طول تماسك تمثيلي مناسب.

λ ≈ 0.082

عامل اقتران الكتلة الموجية التمثيلي.

χ²/دوف ≈ 1.4

جودة الملاءمة الإرشادية في النموذج المبسط.

0. ماذا نحسب ولماذا

منحنى دوران مجرة درب التبانة هو السرعة الدائرية Vc(R) للنجوم كدالة لبعدها R عن مركز المجرة. ويقاس اليوم بدقة أكبر بكثير من التوزيع الكلي للكتلة الذي يمكننا رؤيته مباشرة.

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

تقوم المحاكاة بثلاثة أمور:

يحسب ^Vcbar(R) من القرص المرئي والانتفاخ باستخدام صيغة فريمان التحليلية للقرص.
يدمج عدديًا كثافة نظرية النحل _ρdark(r؛ ℓ، λℓ، λ) عدديًا عند كل نصف قطر ثم يحولها إلى ^VcDM(R).
تصغير χ²(ℓ، λ) مقابل منحنى دوران مجرة درب التبانة في عصر غايا للعثور على أفضل المعلمات التمثيلية الملائمة.

صُممت المحاكاة لتكون قابلة للتكرار. ويمكن تشغيلها بلغة JavaScript أو Python دون الحاجة إلى مكتبة متخصصة في الفيزياء الفلكية.

الترميز المستخدم في جميع الصفحات

R هو نصف القطر الأسطواني المجري الأسطواني في مستوى القرص، بالكيلو بك.
r هو نصف قطر المجرة الكروية.
z هو الارتفاع فوق القرص.
ℓ هو طول تماسك نظرية النحلة، بالكيلو بك.
λ هو اقتران الكتلة الموجية بلا أبعاد.
يُستخدمG بوحدة kpc km² s-² M ⊙ ¹.

\(G=4.302\times10^{-3}\,\mathrm{kpc\,km^2\,s^{-2}\,M_\odot^{-1}}\)

نظرة عامة على خط الأنابيب

بيانات عصر غايا: بناء مجموعة البيانات (_Ri،_Vi،_σi).
النموذج الباريوني: حساب _Vbar(R) من القرص والانتفاخ.
كثافة نظرية النحلة: احسب _ρdark(r; ℓ; λℓ، λ) باستخدام التربيع.
الكتلة المغلقة: دمج الكثافة المظلمة الفعالة في _Mdark(<R).
السرعة الكلية: احسب _Vtot(R) من الباريونات زائد الكتلة المظلمة الفعالة.
χ²تصغير ²: مساحة بارامترات البحث عن أفضل بارامتر و λ.

\(V_{\mathrm{tot}}(R)=\sqrt{V_{\mathrm{bar}}^2(R)+V_{\mathrm{DM}}^2(R)}\)

1. بيانات الرصد – غايا 2024

تعتمد مجموعة البيانات على منحنى دوران مجرة درب التبانة في عصر غايا. وتستخدم أنصاف الأقطار R، والسرعات الدائرية _Vc، والشكوك σ. استخدمت النسخة التقنية الأصلية 16 نقطة بيانات تغطي R = 4-27.3 كيلو بكسل.

الحقائق الملاحظة المهمة هي:

_Vc(R⊙ = 8 kpc) ≈ 230 كم/ثانية، وهي نقطة الارتكاز بالقرب من مدار الشمس.
_Vc مسطح تقريبًا من حوالي 5 إلى 20 كيلو بكسل.
ينخفض _Vc في القرص الخارجي المقاس، ليصل إلى حوالي 173 ± 17 كم/ثانية بالقرب من 27.3 كيلو بك في البيانات المشار إليها.

وهذا يتطلب توزيع كتلة مخفية فعالة ترتفع بقوة عند أنصاف الأقطار المتوسطة، ثم تضعف عند أنصاف الأقطار الأكبر.

\(R_i=\{4,5,6,7,8,9,10,11,12,14,16,18,20,22,24,27.3\}\,\mathrm{kpc}\) \(V_i=\{220,228,232,231,230,229,228,227,226,224,222,219,215,208,200,173\}\,\mathrm{km/s}\)

نستبعد المجرة الداخلية لأن الشريط المركزي والحركات غير الدائرية تهيمن هناك. النموذج المحوري المبسط غير موثوق به داخل تلك المنطقة.

2. نموذج السرعة الباريونية – القرص والانتفاخ

السرعة الدائرية من المادة المرئية هي مجموع تربيع مساهمات القرص والنتوء:

\(V_{\mathrm{bar}}^2(R)=V_{\mathrm{disk}}^2(R)+V_{\mathrm{bulge}}^2(R)\)

2.1 القرص الأسي – معادلة فريمان

يحتوي القرص النجمي الرقيق على كثافة سطحية أسية:

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

بمعلمات تمثيلية:

\(\Sigma_0=800\,M_\odot\,\mathrm{pc}^{-2}\) \(R_d=2.6\,\mathrm{kpc}\)

السرعة الدائرية من قرص أسي متناهٍ في الرقّة الأسية للكتلة الكلية _Md هي

\(V_{\mathrm{disk}}^2(R)=\frac{2GM_d}{R_d}y^2\left[I_0(y)K_0(y)-I_1(y)K_1(y)\right],\qquad y=\frac{R}{2R_d}\)

هنا_In و_Kn هما دالتا بيسل المعدلتان من النوع الأول والثاني. يتم حسابهما عدديًا باستخدام التقريبات التقريبية القياسية كثيرة الحدود والتقريبية التقريبية.

2.2 تقريب الانتفاخ

يُصمم الانتفاخ على شكل مساهمة كتلة كروية مضغوطة:

\(V_{\mathrm{bulge}}^2(R)=\frac{GM_{\mathrm{bulge}}}{R}\)

قد يستخدم نموذج أكثر اكتمالاً نموذجاً أكثر اكتمالاً لنموذج دي فوكولور أو نموذجاً شبيهاً بالقضيب الجانبي ولكن خارج الكيلوباريسكالس القليلة الأعمق يكون هذا التقريب كافياً لمحاكاة من الدرجة الأولى.

المعلمة	الرمز	القيمة	المعنى
ثابت الجاذبية	G	4.302 × 10-³	⊙ ¹ ⊙ م ⊙ كم² ث ² ⊙ م
نصف قطر مقياس القرص	_الطريق	2.6 كيلو متر مكعب	مقياس القرص الرقيق التمثيلي
كتلة القرص	_{دكتوراه}	3.5 × 10¹⁰ M⊙	كتلة القرص النجمي التقريبية
كتلة الانتفاخ	_{ميغابايت}	1.2 × 10¹⁰ M⊙	كتلة الانتفاخ التقريبية

يتم الاحتفاظ بالبارامترات الباريونية ثابتة لأنها مقيدة بشكل مستقل بمجموعات النجوم والقياس الضوئي. إن تركيبها في وقت واحد مع ℓ و λ من شأنه أن يخلق انحرافات قوية.

3. معادلة الكثافة المظلمة لنظرية النحلة

3.1 الفرضية الفيزيائية

يولد كل عنصر كتلة مرئي dV عند الموضع r′ في قرص المجرة مجالاً لموجات الجاذبية يولد كثافة كتلة فعالة عند نقطة المجال r:

\(d\rho_{\mathrm{dark}}(\mathbf{r})=\frac{\lambda}{\ell}\rho_{\mathrm{vis}}(\mathbf{r}’)\exp\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right)dV\)

الكثافة المظلمة الكلية عند أي نقطة (R، z) هي التراكب على جميع عناصر القرص. وبما أن المصدر عبارة عن قرص، فإن حد الكثافة الحجمية يصبح حد الكثافة السطحية:

\(\rho_{\mathrm{vis}}dV\right \Sigma(R’)R’\,dR’\,d\phi\)

التكامل المزدوج الكامل هو:

\(\rho_{\mathrm{dark}}(R,z)=\frac{\lambda}{\ell}\int_0^\infty\int_0^{2\pi}\Sigma(R’)\exp\left(-\frac{\sqrt{R^2+R’^2-2RR’\cos\phi+z^2}}{\ell}\right)R’\,d\phi\,dR’\)

3.2 النواة أحادية القطب والتكامل السمتي

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

\(K(r,R’)=\int_0^{2\pi}e^{-D(r,R’,\phi)/\ell}\,d\phi\approx\frac{2\pi\ell}{r}\sinh\left(\frac{r}{\ell}\right)e^{-(r+R’)/\ell}\)

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

بعد التعويض بـ K، تُختزل الكثافة المظلمة إلى تكامل أحادي البُعد على R′:

\(\rho_{\mathrm{dark}}(r;\ell,\lambda)=\frac{\lambda\Sigma_0}{\ell}\int_0^\infty R’e^{-R’/R_d}\frac{2\pi\ell}{r}\sinh\left(\frac{r}{\ell}\right)e^{-(r+R’)/\ell}dR’\)

3.3 التحقق التحليلي من الحد التقاربي

بالنسبة إلى_Rd ≪ r ≪ ≪ ℓ، يمكن تبسيط التعبير. يكون القرص أصغر بكثير من نصف القطر، ويظل طول التماسك أكبر بكثير من نصف القطر.

\(\rho_{\mathrm{dark}}(r)\xrightarrow{R_d\ll r\ll\ell}\frac{2\pi\lambda\Sigma_0R_d^2}{r^2}\)

لأن:

\(\int_0^\infty R’e^{-R’/R_d}dR’=R_d^2\) \(\sinh\left(\frac{r}{\ell}\right)\approx\frac{r}{\ell}\) \(e^{-r/\ell}\approx1\)

تتصرف هذه الكثافة التقاربية على الصورة ρ ∝ r-²، وهو ما يعطي M( \(\rho(r)\propto r^^{-2} \Longrightarrow M(<r)\propto r\Longrightarrow V_c=sqrt{\frac{GM(<r)}{r}} \approx\mathrm{constant}\)

عند وضع القيم التمثيلية في القيم التمثيلية نحصل على كثافة محلية قريبة من قيمة مجرة درب التبانة المرصودة:

\(\rho_{\mathrm{dark}}(8\,\mathrm{kpc})\approx\frac{2\pi(0.082)(800\,M_\odot\,\mathrm{pc}^{-2})(2.6\,\mathrm{kpc})^2}{(8\,\mathrm{kpc})^2}\approx0.38\,\mathrm{GeV/cm^3}\)

4. مخطط التكامل العددي

4.1 الخطوة 1 – _ρdark(r) بواسطة تربيع نقطة المنتصف

يتم اقتطاع تكامل حلقة المصدر-الحلقة على R′ عند R′ كحد _أقصى = 30 كيلو بكسل، حيث تكون كثافة سطح القرص الأسية بعدها مهملة.

يستخدم التكامل قاعدة نقطة المنتصف مع عقد مصدر N:

\(\rho_{\mathrm{dark}}(r;\ell,\lambda)\approx\frac{\lambda}{\ell}\sum_{i=0}^{N-1}\Sigma_0e^{-R’_i/R_d}K(r,R’_i)\Delta R’,\qquad R’_i=\left(i+\frac{1}{2}\right)\frac{30}{N}\)

حيث:

\(K(r,R’_i)=\frac{2\pi\ell}{r}\sinh\left(\frac{r}{\ell}\right)e^{-(r+R’_i)/\ell}\)

4.2 4.2 الخطوة 2 – الكتلة المظلمة المغلقة

بعد حساب _ρdark(r)، يتم الحصول على الكتلة المظلمة المحصورة داخل نصف القطر R باستخدام تكامل الغلاف الكروي:

\(M_{\mathrm{dark}}(<R)=\int_0^R4\pi r^2\rho_{\mathrm{dark}}(r)\,dr\)

عددياً:

\(M_{\mathrm{dark}}(<R)\approx\sum_{j=0}^{N_r-1}4\pi r_j^2\rho_{\mathrm{dark}}(r_j)\Delta r\)

4.3 الخطوة 3 – السرعة الدائرية للمادة المظلمة

\(V_{\mathrm{DM}}(R)=\sqrt{\frac{GM_{\mathrm{dark}}(<R)}{R}}\)

إذن السرعة الدائرية الكلية هي:

\(V_{\mathrm{total}}(R)=\sqrt{V_{\mathrm{bar}}^2(R)+V_{\mathrm{DM}}^2(R)}\)

4.4 تحويل الوحدة

ينتج تكامل الكثافة كثافة بالكتلة الشمسية لكل كيلوباريسك مكعب. للمقارنة مع كثافة المادة المظلمة المحلية المتعارف عليها بوحدة GeV/cm³، استخدم

\(1\,\frac{M_\odot}{\mathrm{kpc}^3}=\frac{1.989\times10^{30}\,\mathrm{kg}\times5.61\times10^{26}\,\mathrm{GeV/kg}}{(3.086\times10^{21}\,\mathrm{cm})^3}\) \(1\,M_\odot\,\mathrm{kpc}^{-3}\approx3.778\times10^{-2}\,\mathrm{GeV\,cm^{-3}}\)

5. χ² التقليل إلى الحد الأدنى وتركيب البارامترات

5.1 الدالة الهدفية

يقلل الملاءمة من مربع تشي المخفض:

\(\chi_\nu^2(\ell,\lambda)=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_c^{\mathrm{model}}(R_i;\ell,\lambda)-V_{c,i}}{\sigma_i}\right)^2\)

مع N = 16 نقطة بيانات و p = 2 من المعلمات الحرة، ℓ و λ. هذا يعطي 14 درجة من درجات الحرية.

5.2 بحث الشبكة ذات المسارين

يتم استخدام البحث الشبكي بدلًا من نزول التدرج لأن المشهد يحتوي على وادي انحطاط طويل ومنحني بين λ و ℓ.

الممر 1: شبكة خشنة على ℓ و λ.
المسار 2: التنقية المحلية حول الحد الأدنى الخشن.

المنطقة الأنسب تمثيلياً هي

\(\ell \ll \approx130\، \mathrm{kpc}، \qquad \lambda \mambda\max0.082 تقريبًا، \qquad \ qquad \chi^2/\mathrm{dof} \mathrm{of}\max1.4 تقريبًا\)

5.3 وادي الانحطاط

المشهد χ² ليس وعاءً متماثلًا. فهو يشكِّل واديًا ممدودًا لأن تطبيع الكثافة الرائدة يعتمد بقوة على قوة الاقتران وبصورة ضعيفة فقط على طول التماسك داخل نظام الدوران المسطح.

يكسر الانخفاض الخارجي لمنحنى الدوران هذا الانحطاط لأن قيم ℓ الأصغر تكبح الكثافة الفعالة في وقت مبكر.

تُعد البيانات التي تتجاوز 30 كيلو بكسل، بما في ذلك العناقيد الكروية والتدفقات النجمية ونجوم الهالة والمجرات الساتلية، ضرورية لتقييد ℓ بشكل أكثر إحكامًا.

6. التقارب وميزانية الخطأ

تختبر المحاكاة مدى حساسية الخرج لعدد عقد التربيع المستخدمة في تكامل المصدر وتكامل الكتلة الشعاعية.

عقد المصدر N	ρ(8 كيلو بكسل)	التغير النسبي	χ²/دوف	وقت التشغيل
10	10.83	3.2%	1.52	سريع
20	10.98	1.8%	1.45	سريع
40	11.08	0.9%	1.41	اختيار الإنتاج
80	11.13	0.4%	1.40	أبطأ
200	11.15	0.2%	1.39	التحقق من الصحة

N = 40 يعطي دقة دون النسبة المئوية للكثافة وقيم χ² متقاربة تقريبًا. الخطأ العددي أصغر من أوجه عدم اليقين في الرصد والنمذجة.

مصادر الخطأ الرئيسية

مصدر الخطأ	التأثير	التخفيف
التقريب أحادي القطب	يؤثر على أنصاف الأقطار الداخلية	استخدام النواة الزاوية الدقيقة
قرص سميك مفقود	التحولات λ	إضافة مكون القرص السميك
قرص الغاز المفقود	تغيير المظهر الخارجي	إضافة أقراص غاز HI و H₂ H₂
نظامية جايا	يؤثر على منحنى السرعة الخارجية	استخدام مصفوفة التباين الكامل
تقريب التماثل الكروي	يؤثر على تسطيح الهالة	استخدام نواة ثلاثية الأبعاد كاملة

عدم اليقين المهيمن ليس التكامل العددي. بل هو النمذجة الفيزيائية: التحلل الباريوني، وتقريب النواة، وبيانات الهالة الخارجية، والعلاقة الدقيقة بين المعادلة الموجية النحلة والنواة الأسية.

7. الرمز المرجعي الكامل

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

// الثوابت الفيزيائية
const G = 4.302e-3؛ // kpc km² s-² M☉-¹
ثابت Sig0 = 800.0؛ // م☉ كم² ²
const Rd = 2.6؛ // kpc
يشكل Mdisk = 3.5e10؛ // M☉ م☉
تشكل Mbulge = 1.2e10؛ // M☉
كونت CONV_RHO = (1.989e30 * 5.61e26) / (3.086e21)**3;

// بيانات دوران عصر غايا
const OBS_R = [4,5,6,6,7,8,8,9,10,10,11,12,12,14,16,16,18,20,22,24,27.3];
كونت OBS_V = [220,228,232,232,231,230,229,229,228,228,227,227,226,226,224,222,222,219,215,208,200,173];
const OBS_ERR = [10,8,7,7,7,7,6,6,6,6,6,6,6,6,7,7,7,8,8,9,10,11,13,17];

// عنصر نائب للسرعة الباريونية
دالة vBaryonic(R) {
  // في التنفيذ الكامل، يستخدم هذا معادلة قرص فريمان
  // مع دوال بيسل المعدلة I0، I1، K0، K1.
  const vb2 = G * Mbulge / Math.max(R, 0.2);
  إرجاع Math.sqrt(Math.max(0, vb2));
}

// الكثافة المظلمة للنحل
الدالة rhoDark(r، ell، lam) {
  const N = 40;
  const dRp = 30.0 / N;
  دع المجموع = 0;

  بالنسبة إلى (دع i = 0؛ i < N؛ i++) {
    const Rp = (i + 0.5) * dRp;
    const Sig = Sig0 * Math.exp(-Rp / Rd);
    كونت K = (2 * Math.PI * ell / r)
      * Math.sinh(Math.min(r / ell، 30))
      * Math.exp(-Math.min((r + Rp) / ell, 30));

    المجموع += Sig * Rp * K * dRp;
  }

  إرجاع (lam / ell) * المجموع;
}

// الكتلة المظلمة المغلقة
دالة enclosedDarkMass(R, ell, lam) {
  const Nr = 30;
  const dr = R / Nr;
  دع M = 0;

  لـ (دع j = 0؛ j < Nr؛ j++) {
    const rj = (j + 0.5) * dr;
    M += 4 * Math.PI * rj * rj * rj * rhoDark(rj, ell, lam) * dr;
  }

  إرجاع M;
}

// السرعة المظلمة والكلية
دالة vDM(R، إيل، لام) {
  الإرجاع Math.sqrt(Math.max(0, G * enclosedDarkMass(R, ell, lam) / R));
}

دالة vTotal(R, ell, lam) {{
  const vb = vBaryonic(R);
  const vd = vDM(R, ell, lam);
  إرجاع Math.sqrt(vb * vb + vd * vd);
}

// تشي تربيع
دالة تشي تربيع (إيل، لام) {
  دع s = 0;

  ل (دع i = 0؛ i < OBS_R.length; i++) {
    const dv = (vTotal(OBS_R[i]، ell، lam) - OBS_V[i]) / OBS_ERR[i];
    s += dv * dv;
  }

  إرجاع s / (OBS_R.length - 2);
}

يجب أن يتضمن الإصدار الكامل تطبيقات دقيقة لدوال بيسل المعدلة _I0_وI1_وK0_وK1 لسرعة قرص فريمان.

8. كيفية إعادة إنتاج هذه المحاكاة

8.1 في المتصفح

افتح أي متصفح حديث.
افتح وحدة تحكم المطورين.
الصق تطبيق JavaScript كاملاً.
قم بتشغيل دالة التركيب أو قم بإيجاد قيمة χ²² يدويًا لـ ℓ وλ المختارين.

8.2 في بايثون

تُترجم الخوارزمية نفسها مباشرةً إلى Python و NumPy. في لغة Python، استخدم scipy.special.iv و scipy.special.kv لدوال Bessel، والتي هي أكثر دقة من التقريبات متعددة الحدود المشفرة يدويًا.

استيراد numpy ك np
من scipy.special استيراد iv، kv
من scipy.optimize.optimize استيراد تصغير

G = 4.302e-3
سيج0 = 800.0
ج = 2.6
Mdisk = 3.5e10
Mbulge = 1.2e10

OBS_R = np.array([4,5,6,7,8,9,10,11,12,14,16,18,20,22,24,27.3])
OBS_V = np.array([220,228,232,231,230,229,228,227,226,224,222,219,215,208,200,173])
OBS_ERR = np.array([10,8,7,7,7,7,6,6,6,6,6,6,7,7,7,8,8,9,10,11,13,17])

# تنفيذ v_baryonic و rho_dark و enclosed_dark و chi2
# باستخدام نفس الصيغ الموضحة أعلاه.

يجب أن يتقارب مُحسِّن نيلدر-ميد إلى نفس المنطقة الفيزيائية التي يتقارب فيها بحث شبكة جافا سكريبت، حيث يبلغ ℓℓ حوالي 130 kpc و λ حوالي 0.08 في النموذج المبسط.

8.3 تمديدات لملاءمة جودة النشر

استبدل النواة الأحادية القطب بنواة الدالة الزاوية أو نواة دالة بيسل الدقيقة.
أضف مكون قرص سميك.
إضافة أقراص الغازات الذرية والجزيئية.
تضمين شريط المجرة والانتفاخ بشكل أكثر دقة.
استخدم نظام بايزي MCMC لتعيين التوزيع الخلفي لـ ℓ و λ.
تضمين بيانات العنقود الكروي والمجرة الفضائية وبيانات التدفق النجمي حتى 200 كيلو بكسل.

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

المراجع

أبراموفيتز، م.، ستيغون، أ. أ. - كتيب الدوال الرياضية، دوفر، 1972.
Bovy, J., Rix, H.-W. - A Direct Dynamical Measurementical Measurement of the Milky Way's Disk Surface Density Profile, ApJ 779, 115, 2013.
فريمان، ك. س. - عن أقراص المجرات الحلزونية ومجرات S0، ApJ 160, 811, 1970.
McMillan, P. J. - التوزيع الكتلي وإمكانات الجاذبية لدرب التبانة، MNRAS 465, 76, 2017.
Ou, X., Eilers, A.-C., Necib, L., Frebel, A. - ملف تعريف المادة المظلمة لمجرة درب التبانة المستدل عليه من منحنى سرعتها الدائرية، MNRAS 528, 693, 2024.
Pato, M., Iocco, F., Bertone, G. - القيود الديناميكية على توزيع المادة المظلمة في درب التبانة، JCAP 12, 001, 2015.
Portail, M. et al. - النمذجة الديناميكية للانتفاخ المجري والبار، MNRAS 465, 1621, 2017.

بيان قابلية التكرار النهائي

هذه المحاكاة ليست برهاناً نهائياً لنظرية النحل. إنها إطار عددي قابل للتكرار.

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

وتتمثل الخطوة العلمية التالية في استبدال التقديرات التقريبية بنواة دقيقة، وتوسيع النموذج الباريوني، ونشر أوجه عدم اليقين، واختبار نفس الإطار على مجرات مستقلة ومجرات عنقودية.