نظرية النحل – الأسس – المذكرة الفنية الحادية عشرة
تحديد المعلمة المفقودة:
الخطوة 1 – تحليل الارتباط المنهجي
قبل تعديل النموذج، تُشخِّص هذه الملاحظة أي متغير قابل للرصد يتنبأ بالخطأ المتبقي على أفضل وجه. بالعمل على مجموعة المعايرة المكونة من 22 مجرة الواردة في الملاحظة الثامنة، نختبر ارتباط خطأ التنبؤ بكل متغير ذي معنى فيزيائيًا، ثم بكل تركيبة ثنائية المتغيرات، لنحدد بدقة ما أغفله النموذج الحالي.
1. النتيجة أولاً
البارامتر المفقود هو كثافة السطح المركزية
كثافة السطح الباريوني المركزية $\Sigma_d$ لها أقوى ارتباط غير تافه مع خطأ التنبؤ: $r = +0.62$، وR^2 $ = 0.39$ بمفردها.
إن الجمع بين $\سيجما_د$ مع حجم القرص $R_d$ في نموذج ثنائي المتغيرات يفسر $R^2 = 0.43$ من التباين المتبقي، مقارنةً بـ $R^2 = 0.07$ مع $R_d$ وحده. تنخفض النسبة المتبقية من RMS من 19.5\%$ إلى 14.9\%$.
بعد امتصاص كل من $R_d$ و$R_d$ و$Sigma_d$، لا توجد معلومات إضافية يمكن ملاحظتها فيزيائيًا تحمل معلومات عن المتبقي.
2. الطريقة
من خلال العمل على مجموعة المعايرة المكونة من 22 مجرة (الملاحظة الثامنة)، لكل مجرة لدينا خطأ التنبؤ $$ نص{err} = (V_text{tot} – V_f)/V_f$ وقائمة بالمتغيرات الفيزيائية القابلة للقياس. نحسب ارتباطات بيرسون وسبيرمان بين الخطأ وكل متغير مرشح، ثم نختبر الانحدارات ثنائية المتغيرات على الشكل
\$$$ \$ نصّ \${ر} (\%) \%) \;=\\; أ \$ \$ \nنص \\\د \\+\; ب \$ \nنص \ س \\;+\; ج \$$$$$
حيث $X$ هو كل متغير مرشح. أفضل متغير $X$ هو المتغير الذي يزيد من التباين المفسر $R^2$ على المجرات ال 22. يتم استبعاد المتغيرات ذاتية الإحالة – تلك المستمدة من مخرجات النموذج، مثل $V_\\نص{موجة}$ أو $V_\نص{tot}$ – من البحث، لأن ارتباطها بالخطأ أمر لا معنى له.
3. الارتباطات أحادية المتغير
المتغيرات ال 24 المرشحة التي تم اختبارها، مرتبة حسب ارتباط بيرسون المطلق بالخطأ. الصفوف المظللة باللون الذهبي هي متغيرات مشتقة من النموذج نفسه (حشو)؛ أما الصفوف المظللة باللون الأحمر فهي متغيرات حقيقية ذات مرصودات فيزيائية حقيقية مع |r|> 0.5$.
| متغير | الوصف | بيرسون $ r$ | قيمة $ p$ | الأهمية |
|---|---|---|---|---|
| Vw_over_Vf | نسبة Vw / Vf | +0.974 | 0.0000 | ★★★ |
| في_ديناميكال | V_dyn = √ (GM_bar/Rd) | +0.632 | 0.0021 | ★★★ |
| لوغاريتم_سيجما_د | لوغاريتم ₁₀(Σ_d) | +0.622 | 0.0026 | ★★★ |
| م_غاز | كتلة الغاز (M_sun) | +0.609 | 0.0034 | ★★★ |
| M_HI | كتلة HI (M_sun) | +0.609 | 0.0034 | ★★★ |
| T | نوع هابل | -0.585 | 0.0053 | ★★ |
| فبار | باريوني Vبار (كم/ثانية) | +0.582 | 0.0057 | ★★ |
| M_bar_أكثر_من_Rd2 | M_bar / Rd² | +0.559 | 0.0084 | ★★ |
| فتوت | Vtot المتوقع (كم/ثانية) | +0.555 | 0.0090 | ★★ |
| فولكس فاجن | الموجة Vw (كم/ثانية) | +0.550 | 0.0098 | ★★ |
| Vbar_over_v_f | نسبة Vbar / Vf | +0.519 | 0.0158 | ★★ |
| لوغاريتم_م_غاز | لوغاريتم ₁₀(M_gas) | +0.506 | 0.0193 | ★★ |
| لوغاريتم_م_بار | لوغاريتم ₁₀(M_bar) | +0.505 | 0.0196 | ★★ |
| M_bar | الكتلة الباريونية (M_sun) | +0.498 | 0.0214 | ★★ |
| لوج_م_ستار | لوغاريتم ₁₁₀(M_star) | +0.449 | 0.0414 | ★★ |
| سيجما_د | كثافة السطح (لتر/متر مربع) | +0.426 | 0.0544 | ★★ |
| M_star_over_over_Rd2 | M_star / Rd² | +0.426 | 0.0544 | ★★ |
| م_ستار | الكتلة النجمية (M_sun) | +0.389 | 0.0815 | ★ |
قراءة الجدول
أعلى ارتباط منفرد هو $V_\\نص{موجة}/V_f = +0.974$. هذا أمر بديهي: فبحسب البناء، يتدرج الخطأ مباشرةً مع $V_\النص{الموجة}$، لذا يعكس هذا المتغير ببساطة بنية صيغة التنبؤ، وليس محركًا فيزيائيًا خارجيًا.
ومن بين المراصد الفيزيائية الحقيقية، فإن أعلى الارتباطات هي $\log(\Sigma_d) = +0.622$، و$V_\text{ديناميكي} = +0.632$، و$M_\text{غاز} = +0.609$، وT من نوع هابل = -0.585$. وترتبط هذه الإشارات الأربع فيزيائيًا: تميل الأقراص الكثيفة إلى أن تكون أكثر كتلة، ومن النوع السابق، وذات سرعة ديناميكية باريونية أعلى. والسؤال هو ما هو المحرك الأساسي.
4. تصفية المتغيرات الزائدة عن الحاجة
يرتبط العديد من المتغيرات ذات الارتباط الأعلى في حد ذاتها ارتباطًا قويًا بـ $ R_d$، وهو المتغير المعروف بالفعل أنه يقود الخطأ. السؤال هو أيهما يحمل معلومات مستقلة.
| متغير | الارتباط مع $ R_d$ | الحالة |
|---|---|---|
| $\لوج (M_\النجم) $ | $r = +0.88$ | زائدة عن الحاجة مع $ R_d$ |
| $\المدونة (M_text{bar})$$ | $r = +0.87$ | زائدة عن الحاجة مع $ R_d$ |
| $\المدونة (M_text{غاز})$$ | $r = +0.86$ | زائدة عن الحاجة مع $ R_d$ |
| نوع هابل $T$ | $r = -0.66$ | زائدة عن الحاجة جزئياً |
| $_V_نص{ديناميكي}$$ | $r = +0.50$ | مستقل جزئياً |
| $_M_نص_نص_{بار}/R_d^2$ | $r = -0.19$ | مستقل |
| $\لوغ (\سيغما_د) $ | $r = +0.10$ | مستقل |
وترتبط الكتل بـ $R_d$ بشكل مثالي تقريبًا: فالقرص الأكبر يحتوي ببساطة على المزيد من المواد الباريونية. وبالتالي فإن هذه المتغيرات تحمل في الأساس نفس المعلومات التي يحملها $R_d$ نفسه. وعلى النقيض من ذلك، فإن \”سيجما_د$ (كثافة السطح المركزية) و$M_\\text_{bar}/R_d_^2$ (متوسط كثافة السطح الباريوني) متعامدين تقريبًا مع \R_d$ في هذه العينة: فهما يلتقطان الخاصية الهيكلية لـ “مدى انضغاط المادة”، بشكل مستقل عن “مدى امتداد القرص”.
5. الخطأ مقابل كثافة السطح – التصور
رسم الخطأ مقابل $ \log_{10}(\سيجما_د) $ وحده، ملونًا حسب نوع هابل:
إن الاتجاه واضح ورتيب: المجرات ذات الكثافة السطحية المركزية الأعلى يتم التنبؤ بها بشكل منهجي أكثر من اللازم بواسطة BeeTheory، في حين أن الأقراص المنتشرة منخفضة الكثافة يتم التنبؤ بها أقل من المتوقع. يتطابق الميل المناسب البالغ 33$ نقطة مئوية لكل عقد من $\سيجما_d$ مع البيانات بشكل قوي عبر النطاق الكامل من 15 إلى 605 $L_\odot/\\text{pc}^2$.
6. نماذج ثنائية المتغيرات – مقارنة
إضافة $ R_d$ إلى كل متغير مرشح يعطي ترتيبًا أوضح. ويوضح الجدول أدناه التباين المفسر $R^2$ عند إقران $R_d$ مع كل متغير ثانٍ (باستثناء التركيبات الحشو):
| نموذج ثنائي المتغيرات | $R^2$ | متبقي RMS | الملاحظات |
|---|---|---|---|
| $\نص{ر} = a R_d + c$ (خط أساس أحادي المتغير) | 0.074 | $19.5\%$ | مرجع، لا يوجد متغير ثان |
| $$\نص{ر} = أ_ص_د + ب_نص_نص_غاز + c$ | 0.101 | $19.3\%$ | تحسن طفيف |
| $\النص{ر} = أ \ر_د + ب \لوغ م \ستار + ج$ | 0.272 | $17.3\%$ | – |
| $$\نص{ر} = أ_ص_د + ب \نص_ف_نص_بار + c$ | 0.345 | $16.4\%$ | – |
| $\نص{ر} = أ_ص_د + ب \لوغ م\نص{غاز} + c$ | 0.359 | $16.3\%$ | – |
| $\text{err} = a R_d + b T + c$$ | 0.367 | $16.2\%$ | – |
| $\نص{ر} = أ \ر_د + ب \لوغ م\نص{بار} + c$ | 0.373 | $16.1\%$ | – |
| $\نص{ر} = أ ر_د + ب\،ف\نص{ديناميكي} + ج$ | 0.402 | $15.7\%$ | قوي |
| $\نص{ر} = أ \ر_د + ب \لوغ \سيغما_د + ج$ | 0.430 | $15.3\%$ | مستقل عن $ R_d$ |
| $\\نص{ر} = a R_d + b (M_\\نص{ب}/R_d^2) + c$$ | 0.459 | $14.9\%$ | أفضل نموذج غير توتيولوجي |
أفضل نموذج ثنائي المتغيرات
\$$$ \\نص{ر} (\%) \؛ =\\؛ أ\، R_d \؛ +\\؛ ب\، \فراك{م_نص{م_نص{بار}}{R_d^2} \\؛ +\\؛ ج، \qquad R ^2 = 0.46$$$$
المتغير $M_\M\text_{bar}/R_d_^2$ هو متوسط كثافة السطح الباريوني للقرص، $\الزاوية \Sigma_\text_{bar} \rangle = M_\\text_{bar}/(\pi R_d_^2)$. وهو يحمل معلومات عن مدى انضغاط المادة المرئية، بشكل مستقل عن حجم القرص. هذا هو المتغير الذي تفشل BeeTheory حاليًا في حسابه.
7. التحقق من الإغلاق – ما يتبقى بعد احتساب $ R_d$ و $ \ Sigma_d$
إذا كان $ R_d$ و $ \log \Sigma_d$ معًا يلتقطان الخلل الهيكلي، فيجب أن تكون بقايا المطابقة ثنائية المتغيرات غير مرتبطة بكل مرصود مادي. اختبار هذا هو فحص الإغلاق الرسمي:
| متغير | الارتباط مع المتبقي | الحالة |
|---|---|---|
| $ R_d$ | $+0.00$ | عن طريق البناء |
| $\لوج \سيغما_د$ | $+0.00$ | عن طريق البناء |
| $\لوج م/ستار$ | $-0.05$ | امتصاص |
| $\المدونة M_text{Bar}$$ | $+0.07$ | امتصاص |
| $\المدخل M_text{غاز}$$ | $+0.14$ | امتصاص |
| نوع هابل $T$ | $-0.04$ | امتصاص |
| $_V_نص{ديناميكي}$$ | $+0.08$ | امتصاص |
| $_V_نص/{شريط}$ | $+0.05$ | امتصاص |
| $F_نص{غاز}$ | $+0.28$ | هامشي؛ أقل من الأهمية |
بعد احتساب $ R_d$ و $ \log \Sigma_d$، لا يحتفظ أي متغير مادي قابل للملاحظة بعلاقة كبيرة مع الخطأ المتبقي. تم التقاط المعلومات الهيكلية في الخطأ بالكامل بواسطة هذين المتغيرين. يتوافق التشتت المتبقي البالغ 15% من RMS$$ مع عدم اليقين الرصدي في معلمات مدخلات SPARC ومع التباين الجوهري من مجرة إلى مجرة الذي لا يتم التقاطه بواسطة أي من هذه الواصفات الإجمالية.
8. التفسير المادي
يستخدم نموذج BeeTheory الحالي نموذج BeeTheory طول مقياس القرص $R_d$ في موضعين: كمقياس مكاني للتوزيع الباريوني (المظهر الأسي $Sigma propto e^^{-R/R_R_d}$) وكطول تماسك نواة الموجة ($ell = c_text{disk},R_d$). تكون سعة المظهر الجانبي الباريوني $\Sigma_0$ ضمنيًا، ويتم قياسها لإعطاء الكتلة النجمية الصحيحة بمجرد تكاملها.
ما تمثله كثافة السطح فيزيائياً
إن متوسط كثافة السطح الباريوني $$langle Sigma_text_{bar} rangle = M_text{bar}/(pi R_d^2)$ هو الكتلة لكل وحدة مساحة من القرص. هناك مجرتان لهما نفس R_d$$ ولكنهما مختلفتان \R_d$، ولهما نفس المدى الهندسي ولكنهما تحتويان على كميات مختلفة من المادة المكدسة. ويتعامل النموذج الحالي مع المدى الهندسي فقط ($R_d$) على أنه ذو صلة بطول التماسك الموجي، متجاهلاً مدى تركيز المادة. هذا هو بالضبط البارامتر الذي يحدده التحليل المتبقي على أنه مفقود.
اتجاه التأثير
الارتباط موجب: يزداد الخطأ مع كثافة السطح. هذا يعني أنه بالنسبة لـ $R_d$ ثابت، فإن الأقراص الأكثر كثافة يتنبأ بها النموذج بشكل مبالغ فيه – مجال الموجة قوي جدًا بالنسبة لمنحنى الدوران. وعلى العكس، بالنسبة إلى $R_d$ معين، فإن النموذج لا يتوقع الأقراص منخفضة الكثافة المنتشرة. هناك تفسير فيزيائي معقول: يجب ألا يعتمد طول تماسك الموجة على المدى الهندسي للمصدر فحسب، بل على تركيزه أيضًا، حيث تنتج المادة الأكثر كثافة استجابة موجية أكثر محلية. ومن الطبيعي أن يؤدي ذلك إلى كبح سعة المجال الموجي في الأقراص ذات الكثافة العالية \”السيغما\” وتعزيزها في الأقراص ذات الكثافة المنخفضة.
9. ملخص الخطوة 1
1. في مجموعة المعايرة المكونة من 22 مجرة، يرتبط خطأ التنبؤ بأكبر قدر من الارتباط بكثافة السطح المركزية \Sigma_d$ ($r = +0.62$) من بين المراصد الفيزيائية الحقيقية.
2. أما المتغيرات الأخرى التي تبدو في البداية مترابطة بقوة (الكتلة النجمية وكتلة الغاز والكتلة الباريونية) فقد تبين أنها زائدة عن الحاجة إلى حد كبير مع $R_d$ (الارتباط $\gq 0.86$ مع $R_d$) وبالتالي فهي لا تحمل معلومات جديدة تذكر.
3. أفضل نموذج غير ثنائي المتغيرات هو \\نص{ر} = a\، R_d + b\، (M_\\text{bar}/R_d^2) + c$، مع R^2 = 0.46$ و RMS المتبقي 14.9\%$. المتغير الثاني هو متوسط الكثافة السطحية الباريونية للقرص.
4. بعد احتساب $ R_d$ و $ \ Sigma_d$، لا يحتفظ أي عنصر آخر يمكن ملاحظته بعلاقة هامة مع المتبقي. يتم إغلاق التشخيص.
5. وقد تم تحديد البارامتر المفقود: نموذج BeeTheory الحالي يأخذ في الحسبان المدى الهندسي للتوزيع الباريوني ($R_d$) ولكن ليس لكثافة سطحه ($\Sigma_d$). وتتمثل الخطوة التالية في دمج $\سيغما_د$ كمدخل ثانٍ لطول التماسك الموجي، ثم إعادة ضبط النموذج على مجموعة المجرات الـ 22.
المراجع. Lelli, F., McGaugh, S. S., Schombert, J. M. – SPARC: نماذج الكتلة لـ 175 مجرات قرصية مع قياس ضوئي لسبيتزر ومنحنيات دوران دقيقة، AJ 152, 157 (2016). – Pearson, K. – المساهمات الرياضية في نظرية التطور III, Phil. Trans. R. Soc. A 187, 253 (1896). معامل الارتباط. – Dutertre, X. – نظرية النحل™: النمذجة القائمة على الموجة للجاذبية، الإصدار 2، BeeTheory.com (2023).
موقع BeeTheory.com – الجاذبية الكمية القائمة على الموجات – الخطوة التشخيصية 1 – © Technoplane S.A.S 2026