نظرية النحلة – الأسس – المذكرة الفنية الرابعة عشرة
الخطوة 1 – درب التبانة:
تطبيق نواة يوكاوا نظرية النحلة
يتم تطبيق المنهجية الواردة في الملاحظة الثانية عشرة على مجرة درب التبانة باستخدام نواة موجة يوكاوا الصريحة \\mathcal{K}(D) = K_0\،(1+\ألفا D)\،E^{-\ألفا D}/D^2$، مع إجراء التكامل بشكل منفصل لكل مكون باريوني وفقاً لهندسته. تُقارن النتيجة نقطة بنقطة بمنحنى دوران غايا 2024 و”الكتلة المفقودة” للنموذج القياسي. تحدد هذه المذكرة خط الأساس لإطار العمل كما هو مطبق مع الصيغة الهندسية الجديدة بالكامل.
1. النتيجة أولاً
نتيجة خط الأساس – درب التبانة مع نواة يوكاوا الصريحة
منحنى الدوران. يستنسخ النموذج سرعة غايا 2024 في حدود 2 كم/ثانية عند R = 4 كم/ثانية عند R = 4 كم/ثانية، لكنه يبالغ في التنبؤ بمقدار 33 كم/ثانية عند نصف القطر الشمسي (R = 8 كم/ثانية عند R = 27.3 كم/ثانية عند R = 27.3 كم/ثانية.
الكتلة المفقودة. تتطابق كتلة المجال الموجيلنظرية BeeTheory مع “الكتلة المفقودة” للنموذج القياسي في حدود 5% عند R = 4$ kpc، لكنها تتجاوزها بمعامل 2.2 عند R = 27.3$ kpc. ينتج النموذج كتلة مجال مظلم أكثر من اللازم عند أنصاف أقطار كبيرة.
الكثافة المحلية. \rho_text_{wave}(R_odot) = 0.72$ GeV/cm³، مقارنةً بالنطاق المرصود 0.39$- 0.45$ GeV/cm³. زيادة في التنبؤ بمعامل 1.7 تقريبًا.
ما يعنيه هذا
تُنتج صيغة يوكاوا الصريحة منحنى دوران مسطحًا للغاية عند أنصاف الأقطار الكبيرة. طول اضمحلال الحقل الموجي $\ell$ طويل جدًّا، مما يسمح للحقل الموجي بمواصلة المساهمة بكتلة خارج القرص المرئي. هذا هو خط الأساس الهيكلي قبل دمج تنقيح الكثافة السطحية المحدد في الملاحظة الحادية عشرة.
2. ما شرعنا في حسابه
إن مجرة درب التبانة هي حالة الاختبار الطبيعية لأنها المجرة التي تمت معايرة الاقتران العالمي $لامبدا$ عليها في الأصل، ولأنه توجد عمليتا رصد مستقلتان للمقارنة:
(أ) منحنى الدوران $V_c(R)$ من Gaia 2024 (Ou et al., MNRAS 528)، والذي يقيس السرعة الدائرية عند عشرة أنصاف أقطار من $R = 2$ kpc إلى $R = 27.3$ kpc مع شكوك إحصائية تتراوح بين 7$-$17$ كم/ثانية. هذه هي السرعة التي يجب أن تستنسخها BeeTheory، من خلال الجمع بين المساهمة الباريونية $V_\\text{Bar}(R)$ مع مساهمة المجال الموجي $V_\text{wave}(R)$.
(ب) “الكتلة المفقودة” $M_text{missing}(. يستدعي التفسير القياسي المادة المظلمة الجسيمية لتوفير هذه الكتلة. تتنبأ BeeTheory بدلاً من ذلك بأن هذه الكتلة هي المجال الموجي $M_\النص{الموجي}(
(ج) كثافة المادة المظلمة المحلية في الشمس، مقيسة عند \rho \rho \approx 0.39$- 0.45$ GeV/cm³ من الدراسات الكينماتيكية للمحيط الشمسي. تتنبأ BeeTheory بقيمة \rho_\\text{wave}(R_\\odot)$ من نفس حساب المجال الموجي.
إن الاتفاق (أو عدم الاتفاق) على هذه الملاحظات الثلاثة يختبر ثلاثة جوانب مختلفة من النموذج: شكل منحنى الدوران، وملف الكتلة المحصورة فيه، وتطبيع الكثافة المحلية.
3. النواة الموجية – الشكل الصريح
يولِّد كل عنصر كتلة باريوني مجالًا موجيًّا ثنائي الكتلة، مع كثافة عند نقطة مجال تفصلها مسافة $D$ تُعطى بالعلاقة
نواة موجة شكل يوكاوا
\$$\mathcal{K}(D) \\;=\\؛ K_0 \cdot \frac{(1 + \ \ \ alpha D)\,e^{-\alpha D}}{D^2}، \qquad \alpha = \frac{1}{\ell}$$$$$$
يضمن التخميد الأسي $$e^^{- \ألفا D}$ أن يكون للحقل الموجي كتلة كلية محدودة – وبدون ذلك فإن تكامل الحقل الموجي سيتباعد عند ما لا نهاية. يأتي العامل المسبق $(1 + ألفا D) $ من الدالة الموجية المنتظمة في الملاحظة I؛ فهو يجعل البنية المكانية للنواة شبه نيوتونية عند $D ll ell$، ويجعلها مكبوتة أسيًّا عند $D gg ell$.
الطول المميز $\ell$$ يعتمد على المكون الذي يولد المجال:
| المكوّن | طول التماسك $ \ell$ (كيلو متر مكعب) | مقياس هندسي |
|---|---|---|
| انتفاخ (3D هيرنكويست) | $\ell_b = c_\\{sph}\\,r_b = 0.41 \times 0.61$ | $0.25$ |
| قرص رقيق | $\ell_\نص_نص_{رقيق} = ج_نص_نص_قرص_،R_d = 3.17 \times 2.6$ | $8.24$ |
| قرص سميك | $ \ell_نص{ سميك} = ج_نص{قرص} \،(1.5\،R_d)$ | $12.36$ |
| حلقة الغاز | $\ell_نص{غاز} = ج_نص{قرص}\،(1.7\،R_d)$ | $14.01$ |
| أذرع حلزونية | $ \ell_نص{ذراع} = ج_نص{ذراع} \,R_d = 2.0 \times 2.6$ | $5.20$ |
4. هندسة التكامل، كل مكون على حدة
لكل مكوِّن، نقوم بتكامل كثافة المصدر الملتف بالنواة باستخدام عنصر الحجم المناسب للهندسة. والنتيجة هي كثافة المجال الموجي \\rho_\\\text{wave}^{(i)}(r)}(r)$ عند نقطة المجال.
4.1 انتفاخ – تكامل الغلاف الكروي
الانتفاخ عبارة عن توزيع كروي ثلاثي الأبعاد. تحتوي كل قذيفة رقيقة نصف قطرها \r’$ على كتلة \rho_b(r’)\r’b(r’)\r،4\pi r’^2\r’$ وتساهم في المجال على مسافة شعاعية \r$ من المركز. في التقريب الأحادي القطب، يكون الفاصل الفعال بين نقطة المجال ونقطة عامة على الغلاف $$D = \sqrt{r^2 + r’^2}$:
$$ \rho_rho_نص \{موجة}^{(ب)} (ص) \؛ = \\\؛ \int_0^^{r_b} \rho_b(r’) \cdot K_0\\,\frac{(1 + \ألفا_b D)\، ه^^{-\ألفا_b D}}{D^2} \cdot 4\pi r’^2 \، د’، \quad D = \sqrt{r^2 + r’^2}$$$$
مع $ \rho_b(r’) = M_b r_b / [2\pi r'(r’+r_b)^3]$ (Hernquist) و $ \alpha_b = 1\0.25 = 4.0$ kpc$ ^{-1}$. يتم قطع التكامل عند 6 \ r_b$، حيث تكون الكثافة بعد ذلك مهملة عدديًا.
4.2 الأقراص الرقيقة والسميكة – تكامل الحلقات متحدة المركز
كل قرص هو توزيع محوري رقيق غير متماثل. يتحلل القرص إلى حلقات متحدة المركز عند نصف قطر المجرة $R’$ وعرضها $dR’$، كل منها يحمل كتلة $\Sigma(R’)\،2\pi R’\،dR’$. تتطلب المساهمة في نقطة المجال عند نصف القطر $r$ (في مستوى القرص) الفصل الفعال $D = \sqrt{r^2 + R’^2}$ في نفس التقريب أحادي القطب:
$$$\rho_rho_\نص{موجة}^{(\نص{رقيقة})} (ص) \؛ =\\؛ \int_0^^{8R_d} \سيغما_نص_نص_{ر}(R’) \cdot K_0\\، \frac{(1 + \ألفا_نص_{ر} د)\، ه^{- \ألفا_نص_{ر} D}}}{D^2} \cdot 2\pi R’\، dR’، \R’، \quad D = \sqrt{r^2 + R’^2}$$$$
مع $\\سيغما_نص_نص_{سميك}(R’) = (M_\\نص_{سميك}/2\pi R_d^2)، و$\\alpha_\نص_نص_{سميك} = 1/8.24 = 0.121$ kpc$^{-1}$$. يتطابق القرص السميك مع مقياسه الخاص: $\سيجما_نص_\نص_{سميك}$، و$R_نص_\نص_{سميك} = 3.9$ كيلو بكسل، و$\ألفا_نص_نص_{سميك} = 0.081$ كيلو بكسل ^{-1}$$.
4.3 حلقة الغاز – تكامل الحلقة مع الاستنفاد المركزي
يحتوي توزيع الغاز على ثقب مركزي (HI مهمل عند R \sim 2 كيلومتر مكعب) ويمتد أبعد من القرص النجمي. ويلتقط المظهر الجانبي $\سيغما_نص_{غاز}(R’) = \سيغما_0\\، \إكسب (-R_\text_{hole}/R’ – R’/R_g)$ كلتا الميزتين:
$$\rho_\text{wave}^{(\text{gas})}(r) \;=\; \int_{R_\text{hole}}^{8R_g} \سيغما_نص_نص_{غاز}(R’) \cdot K_0\\، \frac{(1 + \ألفا_نص_{غاز} D)\، e^^{-\ألفا_نص_{غ} D}}}{D^2} \cdot 2\pi R’ \، dR’$$$
مع $R_g = 4.42$ kpc، و$R_r_g = 4.42$ kpc، و$R_r_g = 2.21$ kpc، و$$ألفا_r_g = 0.071$ kpc ^{-1}$$. يعكس طول التماسك الأطول توزيع الغاز الأكثر امتدادًا.
4.4 الأذرع الحلزونية – التكامل الحلزوني مع سعة أقل ونواة أضيق
The spiral arms carry $10\%$ of the thin-disk surface density and have their own coherence length $\ell_\text{arm} = 5.2$ kpc, narrower than that of the disk to reflect the azimuthal concentration of the arms:
$$\rho_rho_نص_{موجة}^{(\نص_{ذراع})} (ص)؛ = \\\؛ \int_0^^{8R_d} 0.10\، \سيغما_نص_{ذراع}(R’) \cdot K_0\، \frac{(1 + \ألفا_نص_{ذراع} D)\، \e^{-\ألفا_نص_{ذراع} D}}}{D^2} \cdot 2\pi R’ \، dR’$$$
4.5 إجمالي كثافة المجال الموجي
$$$ \rho_rho_\{نص}(r) \؛ = \\؛ \lambda \,\sum_{i} \rho_\\text_{wave}^{(i)}(r)، \quad \lambda = 0.189$$$$$
5. من كثافة الموجة إلى منحنى الدوران
بمجرد معرفة \\rho_\\text{wave}(r)$ عند كل نصف قطر يتم الحصول على إجمالي كتلة المجال الموجي المحصورة بالتكامل الشعاعي:
$$$$M_نص_نص_{موجة}(R) \؛ = \؛ \int_0^^{R} 4\pi r^2 \rho_نص_{موجة}(r) \، dr$$$$$
ثم تجمع السرعة الدائرية المتوقعة بين مساهمات المجال الباريوني ومساهمات المجال الموجي في التربيع، من خلال العلاقة النيوتونية:
$$$$V_c^2(R) \\؛ =\\؛ V_c^2(R) \؛ +\؛ \frac{G\\، M_\_text{wave}(R)}{R}$$$$
6. منحنى الدوران – نتائج كل نقطة على حدة
| R$ (كيلو متر مكعب) | $V_نص_نص_{بار}$ (كم/ثانية) | $M_\text{wave}/10^{10}$ | $V_نص/{موجة}$ (كم/ثانية) | $V_نص{tot}$ (كم/ثانية) | $_V_text{obs}$ (جايا 2024) | $\دلتا دولار |
|---|---|---|---|---|---|---|
| 2.0 | 157.8 | 0.67 | 120.0 | 198.3 | 250 ± 12 | -51.7 |
| 4.0 | 164.1 | 2.54 | 165.4 | 233.0 | 235 ± 10 | -2.0 |
| 6.0 | 163.5 | 5.17 | 192.5 | 252.6 | 230 ± 8 | +22.6 |
| 8.0 | 157.3 | 8.13 | 209.1 | 261.7 | 229 ± 7 | +32.7 |
| 10.0 | 148.7 | 11.18 | 219.3 | 265.0 | 224 ± 8 | +41.0 |
| 12.0 | 139.6 | 14.15 | 225.2 | 265.0 | 217 ± 9 | +48.0 |
| 15.0 | 126.7 | 18.29 | 229.0 | 261.8 | 208 ± 10 | +53.8 |
| 20.0 | 109.3 | 24.10 | 227.7 | 252.6 | 195 ± 12 | +57.6 |
| 25.0 | 96.7 | 28.54 | 221.6 | 241.8 | 180 ± 15 | +61.8 |
| 27.3 | 92.0 | 30.18 | 218.1 | 236.7 | 173 ± 17 | +63.7 |
يتطابق النموذج مع الملاحظة بشكل ممتاز عند R = 4 كم/ثانية (Δ = -2 كم/ثانية) ولكن التنبؤات الزائدة تتزايد مع زيادة نصف القطر. عند نصف القطر الشمسي، يبلغ التوقع الزائد +33 كم/ثانية (4.7 σ أعلى من عدم اليقين في غايا). عند الحد الخارجي $R = 27.3$ kpc، يصل التوقع الزائد إلى +64 كم/ثانية (3.8σ). المنحنى المتوقع مسطح للغاية – يستمر المجال الموجي في المساهمة بكتلة خارج القرص المرئي، لأن القطع الأسي عند $D \sim \ell$ يسمح بذلك.
7. الكتلة الموجية مقابل “الكتلة المفقودة” في النموذج القياسي
بالنسبة لكل نصف قطر نقارن بين ثلاث كميات: الكتلة الباريونية المحصورة (المادة المرئية فقط)، والكتلة الديناميكية التي تتطلبها السرعة المرصودة (قانون نيوتن المطبق على $V_\text{obs}$)، وكتلة المجال الموجي للنموذج البياني. والفرق بين الثانية والأولى هو ما يسميه النموذج القياسي “الكتلة المفقودة”:
$$$$M_نص{م_نص{مفقود}(<R) \؛ =\\؛ \frac{R\، V_نص{مطلوب}^2(R)}{G} \\؛ -\\؛ M_\\نص{بار}(<R)$$$
تخبرنا النسبة $$M_نص{موجة}/نص{موجة}$ بمدى إحلال مجال موجة نظرية النحل محل المادة المظلمة الجسيمية على أساس نصف قطر بنصف قطر:
| R$ (كيلو متر مكعب) | $M_\\نص{بار}(| $M_\\نص{دين}( | $M_\\نص_{مفقود}$$ | $M_نص_نص_{موجة}$ (BT) | نسبة $$M_\\نص{موجة}$$/م_نص{مخطئ}$$ | |
|---|---|---|---|---|---|
| 2.0 | 1.16e+10 | 2.91e+10 | 1.75e+10 | 6.69e+09 | 0.38 |
| 4.0 | 2.51e+10 | 5.13e+10 | 2.63e+10 | 2.54e+10 | 0.97 |
| 6.0 | 3.73e+10 | 7.38e+10 | 3.65e+10 | 5.17e+10 | 1.42 |
| 8.0 | 4.60e+10 | 9.75e+10 | 5.15e+10 | 8.13e+10 | 1.58 |
| 10.0 | 5.14e+10 | 1.17e+11 | 6.52e+10 | 1.12e+11 | 1.71 |
| 12.0 | 5.44e+10 | 1.31e+11 | 7.70e+10 | 1.41e+11 | 1.84 |
| 15.0 | 5.60e+10 | 1.51e+11 | 9.49e+10 | 1.83e+11 | 1.93 |
| 20.0 | 5.56e+10 | 1.77e+11 | 1.21e+11 | 2.41e+11 | 1.99 |
| 25.0 | 5.43e+10 | 1.88e+11 | 1.34e+11 | 2.85e+11 | 2.13 |
| 27.3 | 5.37e+10 | 1.90e+11 | 1.36e+11 | 3.02e+11 | 2.22 |
القراءة الكمية
عند 4$R = 4$ kpc يتطابق الحقل الموجي بشكل أساسي مع الكتلة المفقودة (النسبة 0.97$). بين $R = 6$ و $R = 8$ kpc يتجاوز النموذج بالفعل الكتلة المفقودة بنسبة 40-60%. أما ما بعد R = 15$ kpc، فإن كتلة المجال الموجي تبلغ تقريبًا ضعف ما يستدعيه النموذج القياسي من المادة المظلمة. وينتج النموذج كتلة إضافية عند أنصاف أقطار كبيرة – وهو بالضبط عَرَضُ طول التماسك $$\ll$ الطويل جداً بالنسبة للقرص النجمي المرئي.
8. مساهمات المكوِّنات في المجال الموجي عند نصف القطر الشمسي
يُظهر تقييم مساهمة كل مكوِّن في $\rho_\\text{wave} (R_\odot = 8\، \text{kpc})$ أي مصدر باريوني يهيمن على الحقل الموجي هناك:
| المكوّن | $ \rho_rho_نص_{موجة} (R_\odot) $ ($ M_\odot$/kpc³) | جزء من الإجمالي |
|---|---|---|
| قرص نجمي رقيق | 6.05 دولار في 10^7 دولار في 10^7 دولار | 60.6% |
| قرص نجمي سميك | 1.91 \10^7$ دولار أمريكي | 19.1% |
| حلقة الغاز | 1.62 دولارًا أمريكيًا في 10^7 دولار أمريكي | 16.2% |
| أذرع حلزونية | 4.15 دولار أمريكي في 10^6 دولار أمريكي | 4.1% |
| انتفاخ | 1.55 دولار في 10^^{-5}$ دولار | 0 دولار أمريكي/دولار أمريكي |
يهيمن القرص النجمي الرقيق على الحقل الموجي عند موقع الشمس (60%)، ويساهم القرص السميك وحلقة الغاز بالتساوي تقريبًا (16-19%). يساهم الانتفاخ بشكل لا يُذكر لأن $\ell_b = 0.25$ kpc أصغر بكثير من $R_\\odot = 8$ kpc – فالقمع الأسي في النواة يقتل مساهمة الانتفاخ عند هذه المسافة.
وبتحويل الكثافة الكلية إلى وحدات فيزياء الجسيمات نحصل على \\rho_\text\{wave} (R_\\odot) = 0.717$ GeV/cm³، لمقارنتها بالقياس الحركي الذي يتراوح بين 0.39$- 0.45$ GeV/cm³ من Read 2014 والتحليلات اللاحقة. يتخطى التنبؤ الكثافة المحلية المرصودة بمعامل 1.6$-$1.8$-$1.8$ – بما يتوافق مع التنبؤ الزائد لمنحنى الدوران عند نفس نصف القطر.
9. ما يحدده خط الأساس هذا
تعمل الآلية من حيث المبدأ
حول $R = 4 kpc$ – أي الجسم المركزي للقرص – يساوي الحقل الموجي المتكامل الكتلة المفقودة للنموذج القياسي في حدود 5%، ويتم إعادة إنتاج منحنى الدوران في حدود 2 كم/ثانية. تُنتج النواة الموجية، المطبقة على الباريونات المرئية، كتلة جاذبية مماثلة كمياً للمادة المظلمة الجسيمية عند نصف القطر هذا. ليست هناك حاجة إلى جسيم جديد؛ فالمجال الموجي للمادة المرئية يمثل الجاذبية المفقودة.
لكن المنحنى مسطح للغاية عند أنصاف الأقطار الكبيرة
خارج القرص المركزي، يبالغ النموذج في توقع سرعة الدوران بمقدار ينمو بشكل رتيب مع نصف القطر. يستمر الحقل الموجي بتراكم الكتلة خارج القرص المرئي لأن طول تماسك النواة $\ell_\text{thin} = 8.24$ كيلو بكسل يماثل حجم القرص نفسه، مما يسمح بمساهمات كبيرة عند $D = 15$-$$25$ كيلو بكسل. في المقابل، ينخفض منحنى دوران غايا قليلاً بعد R \sim 10 كيلو بكسل – وهي ميزة لا تعيد الصيغة الحالية إنتاجها.
خط أساس وليس إجابة نهائية
تحدد هذه العملية الحسابية خط الأساس للنموذج مع اعتماد $\ell_i$ خطيًّا على $R_d$ وحده. حدد التشخيص في الملاحظة XI أن $\Sigma_d$ – كثافة السطح المركزية – يجب أن تدخل في تحديد $\ell_i$ لتصحيح المنحنى عند أنصاف الأقطار الكبيرة. فكلما كانت كثافة القرص أكثر كثافة، يجب أن تكون استجابة الموجة أكثر محلية. دمج هذا التنقيح هو موضوع الملاحظات اللاحقة. خط أساس مجرة درب التبانة المذكور هنا هو ما يجب أن تحسنه تلك الملاحظات.
10. ملخص
1. يتم حساب منحنى دوران مجرة درب التبانة عن طريق تكامل كل مكون باريوني مقابل نواة موجة يوكاوا $mathcal{K}(D) = K_0(1 + ألفا D)، ه^^{-alpha D}/D^2$، مع هندسة مناسبة: قذائف كروية للانتفاخ، وحلقات متحدة المركز للأقراص والغاز والأذرع الحلزونية.
2. عند $R = 4 كيلو متر مكعب، تتوافق كتلة المجال الموجي لنظرية BeeTheory مع “الكتلة المفقودة” للنموذج القياسي في حدود 5٪ (النسبة 0.97)، والسرعة المتوقعة تتطابق مع سرعة غايا 2024 في حدود 2 كم/ثانية.
3. عند نصف القطر الشمسي (R = 8 كم/ثانية = 8 كم/ثانية)، يتنبأ النموذج بسرعة دوران زائدة بمقدار 33 كم/ثانية +33 كم/ثانية وكثافة المادة المظلمة المحلية بمعامل 1.6 – وكلاهما متوافق مع الآخر.
4. وفيما يتجاوز R = 15 كيلو متر مكعب، فإن كتلة المجال الموجي المتوقعة تتجاوز الكتلة المفقودة للنموذج القياسي بمعامل 2 أو أكثر. لا ينخفض منحنى الدوران المتوقع كما تتطلب بيانات غايا.
5. يهيمن القرص النجمي الرقيق على المجال الموجي عند موقع الشمس (60٪ من $\rho_\text{wave}$). ويساهم الانتفاخ بشكل ضئيل. يتوافق التحلل مع هندسة التكامل الموصوفة.
6. إن التنبؤ الزائد عند $R$ كبير $R هو التوقيع الهيكلي لـ $\ell_i$ كونه طويلًا جدًا. لاحظ XI المحدد أن $\Sigma_d$ يجب أن يدخل في صيغة طول التماسك. إن تنقيح $\ell_i$ عبر $\Sigma_d$ هو الخطوة التالية.
المراجع. Ou, X. et al. – ملف تعريف المادة المظلمة لمجرة درب التبانة المستدل عليه من منحنى السرعة الدائرية، MNRAS 528, 693 (2024). منحنى دوران غايا 2024. – Bland-Hawthorn, J., Gerhard, O. – المجرة في السياق، ARA&A 54, 529 (2016). MW التحلل الهيكلي. – Read, J. I. – كثافة المادة المظلمة المحلية، J. Phys. G. 41, 063101 (2014). قياسات كثافة المادة المظلمة المحلية. – فريمان، ك. س. – على أقراص المجرات الحلزونية ومجرات S0، ApJ 160, 811 (1970). – هيرنكويست، ل. – نموذج تحليلي للمجرات الكروية والانتفاخات، ApJ 356، 359 (1990). – Dutertre، X. – نظرية النحل™: النمذجة المستندة إلى الموجة للجاذبية، الإصدار 2، BeeTheory.com (2023).
موقع BeeTheory.com – الجاذبية الكمية القائمة على الموجات – تطبيق الخطوة 1 – © Technoplane S.A.S 2026