نظرية النحل – الأسس – المذكرة الفنية الثانية عشرة

إضفاء الطابع الرسمي:
حساب نظرية النحل على نطاق المجرة

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

1. الحساب في مخطط واحد

سلسلة أحادية الاتجاه

Observed photometry $\;\longrightarrow\;$ Baryonic decomposition $(\rho_\text{bar})$
\\Big\downarrow$
التفاف المجال الموجي \\؛ \ لونغ رايتارو \؛ \ كثافة الموجة \(\rho_\\text{موجة})$$
$\Big\downarrow\$
تكامل الكتلة $\;\\\longrightarrow\\\$ كتلة الموجة المغلقة $(M_\\\text\{موجة})$$
$\Big\downarrow\$
العلاقة النيوتونية $\;\\\longrightarrow\\$ منحنى الدوران المتوقع $(V_c)$

لا توجد خطوة مقلوبة. لا يُستخدَم منحنى الدوران $(V_c(R)$ كمدخل.

2. مدخلات الرصد

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

الرمز	الكمية	المصدر
$T$	نوع هابل المورفولوجي	فهرس (دي فوكولور وآخرون 1991، SPARC)
$ R_d$	طول مقياس القرص النجمي (kpc)	القياس الضوئي لسبيتزر 3.6 ميكرومتر (SPARC)
$\سيجما_دولار	سطوع سطح القرص المركزي ($L_\odot/\\text{pc}^2$)	القياس الضوئي لسبيتزر 3.6 ميكرومتر (SPARC)
$_M_نص/{HI}$$	إجمالي كتلة الهيدروجين الذري ($M_odot$)	أرصاد راديو 21 سم (SPARC)
$\أوبسيلون_ستار_دولار	نسبة الكتلة النجمية إلى الضوء عند 3.6 ميكرومتر	عالمي ثابت: 0.5 دولار أمريكي/مليون دولار أمريكي/دولار أمريكي (McGaugh 2014)

بالنسبة لمجرة درب التبانة، يتم استبدال $R_d$ و$Sigma_d$ و$M_text{HI}$ بالقيم المماثلة المحددة من المسوحات النجمية الداخلية (Bovy & Rix 2013) وخرائط 21 سم. يتم استخدام نفس متجه المدخلات الخماسية الكمية.

3. التفكيك الباريوني – خمسة مكونات هندسية

من المدخلات الرصدية الخمسة، يتم تقسيم الكتلة الباريونية إلى خمسة مكونات هندسية متميزة. ويحمل كل مكوِّن ملف الكثافة الخاص به ومقياسه المميز.

3.1 إجمالي الكتل النجمية والغازية

$$$$$M_\\نجم \\، =\\؛ 2\pi\، R_d^2\، \سيغما\، \أبسلون\\ نجم$$$$

$$$$M_\\نص{غاز} \\؛ =\\؛ 1.33\، م \ نص \{HI} \\qquad \\نص \{(تصحيح؛ أرنيت 1996)}$$$$$

3.2 كتل المكونات والموازين

المكوّن	القداس	المقياس	التفعيل
انتفاخ	$ M_b = 0.20\\، M_\ النجم$	\r_b = \max(0.5\،R_d,\,0.3\text{ kpc})$	إذا كان $T \leq 4$
قرص رقيق	\M_نص_نص_نص_رئيسي = 0.75\، (M_\ النجم – M_b)$	$ R_d$	دائماً
قرص سميك	\M_نص_نص_{سميك} = 0.25\، (M_\ النجم – M_b)$	1.5 دولار أمريكي، R_d$	دائماً
حلقة الغاز	$ M_\text{غاز} = 1.33\\، M\\text{HI}$	$R_g = 1.7\،R_d$ (Broeils & Rhee 1997)	دائماً
أذرع حلزونية	\M_نص{ذراع}$ = 0.10\، M_نص{ذراع}$ (فعال)	$ R_d$ (يتبع القرص الرقيق)	دائماً

3.3 ملامح الكثافة

انتفاخ (3D هيرنكويست)

$$$\rho_b(r) \;؛ =\\؛ \frac{M_b\\,r_b}{2\pi\,r\،(r + r_b)^3}$$$$$

الأقراص النجمية الرقيقة والسميكة (2D الأسية)

$$\Sigma_\text{thin}(R) \;=\; \frac{M_\text{thin}}{2\pi\,R_d^2}\,e^{-R/R_d}$$

$$\Sigma_\text{thick}(R) \;=\; \frac{M_\text{thick}}{2\pi\,(1.5\,R_d)^2}\,e^{-R/(1.5R_d)}$$

حلقة الغاز (2D أسي ثنائي الأبعاد مع ثقب مركزي)

\$$$\سيغما_نص_{غاز}(R) \؛ =\؛ \\frac{M_نص_نص_{غ}}{2\pi\\، R_g^2}\، \\\إكسبر\! \ليسار (- \frac{R_نص_نص_{ثقب}}{R} – \frac{R}{R_R_g}\ يمين)، \quad R_\\نص_{ثقب} = 0.5\، R_g$$$$

فائض الذراع الحلزوني (2D، يتبع القرص الرقيق)

$$$\سيغما_نص_{ذراع}(R) \؛ = \؛ 0.10\؛ \سيغما_نص_{ذراع}(R)$$

4. نواة الموجة

يولد كل عنصر كتلة باريوني مجالاً موجيًا نحلًا. الحقل عند نقطة $vec{r}$ الناتج عن عنصر مصدر عند $vec{r}،’$ يفصل بينهما $D = |vec{r} – vec{r}،’|$ تحكمه نواة من نوع يوكاوا مشتقة من الدالة الموجية المنتظمة في الملاحظة الأولى:

نواة موجة النظرية

\$$\mathcal{K_K}_i(D) \\؛ =\\؛ K_0\؛ \frac{(1 + \alpha_i\، D)\، e^^{-\alpha_i\، D}}{D^2}، \qquad \alpha_i \؛ =\؛ \frac{1}{\ll_i}$$$$

هنا $K_0$ هو السعة الشاملة للكتلة الموجية (عدد واحد بلا أبعاد) و $\ell_i$ هو طول تماسك المكون $i$. ترمز النواة إلى سلوك شبه نيوتوني 1$/دي^2$ عند المسافات القصيرة، ويتم تعديلها بقطع أسي عند مقاييس تتجاوز $\ell_i$. تضمن الصيغة $(1 + \ألفا D) \، ه^^{- \ألفا D}$ الاستمرارية والكتلة الكلية المحدودة عند اللانهاية.

4.1 أطوال ترابط المكونات

يتم تعيين طول التماسك لكل مكوّن بمقياسه الهندسي الطبيعي، مضروبًا في ثابت بلا أبعاد خاص بأبعاده:

المكوّن	طول التماسك	الثابت الهندسي
انتفاخ (كرة ثلاثية الأبعاد)	$ \ell_b = c_\\text{sph} \\,r_b$	$_ج_نص{sph}$ دولار أمريكي
قرص رقيق (2D)	$\ell_\نص_نص_{رقيقة} = ج_نص_{قرص}\، ر_د$	$_ج_نص/{قرص}$$
قرص سميك (2D)	$ \ell_نص{ سميك} = ج_نص{قرص} \،(1.5\،R_d)$	$_ج_نص/{قرص}$$
حلقة الغاز (2D)	$\ell_نص{غاز} = ج \ نص{قرص}\، R_g$	$_ج_نص/{قرص}$$
أذرع حلزونية (ثنائية الأبعاد، مركزة سمتيًا)	$\ell_نص_نص_{ذراع} = ج_نص_{ذراع}\، R_d$	$_ج_نص_{ذراع}$

إن الثوابت الهندسية الثلاثة $(c_\{sph}، \، c\{disk}، \، c\{text{arm})$ هي ثوابت عامة – فهي لا تختلف من مجرة إلى أخرى. وبالإضافة إلى سعة الكتلة الموجية العالمية $K_0$ واقتران المجال الموجي $\lambda$، فإنها تشكل المجموعة الكاملة من المعلمات على المستوى النظري.

5. التفاف المجال الموجي – معادلات تكاملية لكل مكون

كثافة المجال الموجي عند موضع $\vec{r}$ هي التقاء توزيع المصدر الباريوني مع نواة الموجة. بالنسبة إلى نظام متماثل مجري (متماثل محوريًا، تقريب أحادي القطب)، يساهم كل مكون باريوني بشكل إضافي:

الكثافة الكلية لمجال الموجة عند نصف القطر $$r$

\$$ \rho_\\text{wave}(r) \؛ =\\؛ \lambda \\؛ \sum_{i \\in \{{{نص،{نص، سميك، غاز، غاز، ذراع، انتفاخ}\}} \\rho_\\نص \{موجة}^{(i)}(r)$$

فيما يلي التكاملات الخمسة مكتوبة أدناه، تكامل واحد لكل مكوِّن. يحوِّل كل تكامل توزيع الكتلة الباريونية إلى توزيع كتلة المجال الموجي عند نفس النقطة المكانية.

5.1 الانتفاخ – تكامل الغلاف ثلاثي الأبعاد

\$$\rho_rho_نص{موجة}^{(ب)}(r)\؛ =\؛ \int_0^^{r_r_نص{ماكس}} \rho_b(r’)؛ \mathcal{K}_b\\!\r_b(\sqrt{r^2 + r’^2}\right)؛ \4\pi r’^2\، dr’$$$

يتم التكامل على قذائف كروية متحدة المركز نصف قطرها $r’$$. وترى نقطة المجال عند نصف القطر $r$ من المركز كل غلاف عند فاصل فعال $D = \sqrt{r^2 + r’^2}$ في التقريب أحادي القطب. ويمتد التكامل إلى \r_text{max} = 6\، r_b$، حيث تكون كثافة الانتفاخ بعد ذلك مهملة عدديًا.

5.2 القرص الرقيق – تكامل الحلقة ثنائية الأبعاد

$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{R_\text{max}} \سيغما_نص_{رقيقة}(R’)؛ \\mathcal{K}_نص_{رقيقة}\! \\ليسار (\sqrt{r^2 + R’^2}\يمين)؛ 2\pi R’\،dR’$$$

يتحلل القرص إلى حلقات متحدة المركز نصف قطرها $R’$ وعرضها متناهي الصغر $dR’$، تحمل كل منها كتلة سطحية $\سيغما_نص{ر}(R’)\،2\pi R’\،dR’$. ينطبق التقريب الأحادي القطب نفسه: يتلقى المجال الموجي عند نصف القطر $ r$ من المركز مساهمات من كل حلقة عند الفصل الفعال $D = \sqrt{r^2 + R’^2}$. نطاق التكامل هو $R_R_text{max} = 8 \، R_d$.

5.3 القرص السميك – تكامل الحلقة ثنائية الأبعاد

$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}} \سيغما_نص_{سميك}(R’)؛ \mathcal{K}_نص_{سميك}\! \ يسار (\sqrt{r^2 + R’^2}\ يمين)؛ 2\pi R’\,dR’$$$

مطابق لتكامل القرص الرقيق، مع $\سيجما_\نص_\{ثخين}(R’)$ ككثافة المصدر ومعامل نواة $\ألفا_\\نص_\{ثخين} = 1/(c_\\text_{disk}\، \cdot 1.5\، R_d)$. ويؤدي المدى الشعاعي الأوسع للقرص السميك إلى نطاق أوسع قليلاً من التماسك الموجي.

5.4 الحلقة الغازية – تكامل الحلقة ثنائية الأبعاد مع استنفاد مركزي

$$\rho_\text{wave}^{(\text{gas})}(r) \;=\; \int_{R_\text{hole}}^{R_\text{max}} \Sigma_\text{gas}(R’)\;\mathcal{K}_\text{gas}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$

يتضمن توزيع الغاز ثقبًا مركزيًا يلتقطه القطع الشعاعي عند R_\text{hole} = 0.5\،R_g$ في الحد الأدنى للتكامل. خارج هذا القطع، يمتد الغاز إلى أبعد من القرص النجمي؛ وينعكس هذا في المقياس المميز الأكبر R_g = 1.7\، R_d$، الذي يغذي طول التماسك $\ell_\\text_{gas} = c_\\text{disk}\، R_g$.

5.5 فائض الذراع الحلزوني – تكامل الحلقات ثنائية الأبعاد مع سعة مخفضة

$$\rho_\text{wave}^{(\text{arm})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{arm}(R’)\;\mathcal{K}_\text{arm}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$

يتم التعامل مع الأذرع الحلزونية على أنها تعزيز متوسط محوري لكثافة سطح القرص الرقيق عند مستوى 10\%$$، مع طول التماسك الخاص بها $\ell_\\text\{arm} = c_\\text\{arm}\،R_d$. وبالتالي فإن النواة أضيق من نواة القرص الرقيق، مما يعكس التركيز السمتي للبنية الحلزونية.

6. كتلة الموجة المغلقة ومنحنى الدوران المتوقع

بمجرد معرفة الكثافة الكلية للمجال الموجي $\rho_\\text{wave}(r)$، يتم الحصول على كتلة المجال الموجي المحصورة داخل كرة نصف قطرها $R$ بالتكامل الشعاعي:

كتلة المجال الموجي المغلق

$$$$M_نص_نص_{موجة}(R) \؛ = \؛ \int_0^^{R} 4\pi\، r^2\، \rho_نص_{موجة}(r)، \rho_نص_{موجة}(r)، \rho_نص_موجة}(r)، \r$$$$$

ثم تأتي السرعة الدائرية المتوقعة عند نصف القطر $R$ من العلاقة النيوتونية، التي تجمع بين مساهمات المجال الباريوني والمجال الموجي في التربيع:

السرعة الدائرية المتوقعة

$$$$V_c^2(R) \\؛ =\\؛ V_c^2(R) \؛ +\\؛ \frac{G\\، M_\\text{موجة}(R)}{R}$$$$$

إن السرعة الباريونية $V_\\text{Bar}(R)$ هي نفسها المجموع التربيعي للمساهمات من المكونات الأربعة الشبيهة بالقرص (صيغة فريمان 1970 لكل ملف أسي) والانتفاخ (صيغة الكتلة المغلقة لهيرنكويست):

$$$$$V_نص{بار}^2(R) \؛ =\؛ V_نص{بولج}^2 + V_نص{رقيق}^2 + V_نص{ثخين}^2 + V_نص{غاز}^2$$$$$$

حيث يكون كل $V_i(R)$ هو السرعة الدائرية النيوتونية القياسية للتوزيع الكتلي المناظر.

7. المعلمات على المستوى النظري

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

الرمز	المعنى	الدور
$K_0$	سعة الكتلة الموجية	يضبط المقياس غير البعدي لنواة الموجة
$_ج_نص{sph}$ دولار أمريكي	ثابت هندسي ثلاثي الأبعاد	النسبة $\ell\r\r_r\r_text{scale}$ للمصادر الكروية (الانتفاخ)
$_ج_نص/{قرص}$$	ثابت هندسي ثنائي الأبعاد 2D	النسبة $\ell/R\\\text{scale}$ لمصادر الأقراص والحلقات
$_ج_نص_{ذراع}$	ثابت هندسي حلزوني حلزوني	النسبة $ \ell/Rell_d$ لفائض الذراع المركزة سمتيًا
دولار أمريكي	اقتران المجال الموجي العالمي	قياس كثافة المجال الموجي الكلي

عالمية المعلمات

جميع البارامترات الخمسة عالمية. تنطبق نفس القيم العددية على مجرة درب التبانة، وعلى المجرات القزمة غير المنتظمة، وعلى اللوالب الضخمة. وتدخل المعلومات الخاصة بالمجرة فقط من خلال المدخلات الرصدية الخمسة $(T، \، R_d، R_d، \، Sigma_d، \، M_text{HI}، \، Upsilon_\star)$. لا يحتوي النموذج على أي بارامتر قابل للضبط لكل مجرة.

8. الطبيعة الأحادية الاتجاه للحساب

سلسلة مفتوحة – لا توجد تغذية مرتدة

يتدفق الحساب بأكمله من المدخلات إلى المخرجات، في اتجاه واحد. تحدد الملاحظات الضوئية وملاحظات 21 سم التحلل الباريوني. يحدد التحلل الباريوني كثافة المجال الموجي. تُحدِّد كثافة المجال الموجي الكتلة الموجية المحصورة. تُحدِّد كتلة الموجة المغلقة منحنى الدوران المتوقَّع. لا يؤثِّر منحنى الدوران في أي مرحلة على أي خطوة سابقة من الحساب.

تترتب على هذه الأحادية الاتجاه ثلاث نتائج مهمة.

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

(ب) لا يحتوي النموذج على آلية لتعديل كل مجرة على حدة. يجب أن يأتي كل تعديل في تنبؤ منحنى الدوران من تعديل متجه المدخلات $(T,\\، R_d، R_d، \، Sigma_d، M_text{HI}، \، Upsilon_\star)$ أو من تغيير في المعلمات على مستوى النظرية العامة $(K_0، \، c_text{sph}، \، c_text{disk}، \، c_text{disk}، \، c_text{arm}، \، \،lambda)$.

(ج) معايرة $lambda$ على مجرة مرجعية ليست هي نفسها التي تلائم منحنى دوران تلك المجرة. فالمعايرة تحدِّد رقمًا عالميًّا واحدًا؛ ومنحنى الدوران عند جميع أنصاف الأقطار الأخرى للمجرة المرجعية، ومنحنيات دوران جميع المجرات الأخرى، هي بعد ذلك تنبؤات صارمة للإطار المعاير.

9. دور الكثافة السطحية المركزية (مراجعة الملاحظة الحادية عشرة)

حدد التشخيص الوارد في الملاحظة الحادية عشرة أن خطأ التنبؤ المتبقي يرتبط بقوة بكثافة السطح الباريوني المركزية $\Sigma_d$، بشكل مستقل عن طول مقياس القرص $R_d$. الصيغة المعروضة أعلاه هي نسخة النموذج قبل دمج هذه النتيجة – فهي تستخدم فقط $R_d$ في تعبيرات طول التماسك $\ell_i = c_i\،R_d$.

حيث سيدخل التنقيح

في النموذج المحسّن، ستعتمد أطوال التماسك $\ll_i$ على كل من $R_d$ و $\Sigma_d$، مع استبدال العلاقة الخطية الصارمة $\ell_i = c_i\\، R_d$ بدالة $\ell_i = c_i\، R_d\، \phi(\Sigma_d\\Sigma_\{ref})$ التي تمتص البقايا المحددة في الملاحظة الحادية عشرة. سيتم تحديد الشكل الوظيفي لـ \phi$ ومعلماته في الملاحظات اللاحقة، أولاً على مجموعة المعايرة المكونة من 22 مجرة، ثم يتم التحقق من صحتها عن طريق التنبؤ الأعمى على عينة SPARC المتبقية.

يتم الحفاظ على البنية أحادية الاتجاه للحساب من خلال هذا التنقيح: $\Sigma_d$ هو مدخل رصد، وتغذي أطوال التماسك المعدلة نفس تكاملات الالتفاف، ويظهر منحنى الدوران كما كان من قبل. تمت إضافة رابط تشغيلي واحد فقط – وهو اعتماد $\ell_i$ على مرصود ثانٍ.

10. ملخص المنهجية

1. المدخلات. خمسة مواد قابلة للرصد لكل مجرة: نوع هابل $T$، ومقياس القرص $R_d$، والسطوع السطحي $\Sigma_d$، وكتلة HI $M_\text{HI}$، ونسبة الكتلة النجمية الشاملة إلى الضوء $Upsilon_\star$.

2. التحلل الباريوني. خمسة مكونات: الانتفاخ (إذا كان $T \leq 4$)، القرص الرقيق، القرص السميك، الحلقة الغازية، فائض الذراع الحلزوني. يحمل كل منها ملف تعريف كثافة تحليلي.

3. نواة الموجة. صيغة يوكاوا العالمية من نوع Yukawa $\mathcal{K}K_i(D) = K_0\،(1 + \alpha_i D)\، هـ^{-\alpha_i D}/D^2$ مع طول التماسك $\ell_i = c_i\، R_\\text{scale}$ الذي يحدده المدى الهندسي لكل مكون.

4. الالتفاف. يولِّد كل مكوِّن كثافة مجال الموجة من خلال تكامل أحادي البُعد على الحلقات (مكوِّنات ثنائية الأبعاد) أو الأصداف (انتفاخ ثلاثي الأبعاد). وكثافة المجال الموجي الكلية هي مجموع المكونات الخمسة، مقيسة بالاقتران الكلي $\lambda$.

5. الناتج. يتم تكامل الكتلة الموجية المغلقة $M_\\نص{موجة}(R)$ ودمجها مع السرعة الباريونية $V_\نص{موجة}(R)$ للحصول على منحنى الدوران المتوقع $V_c(R)$.

6. معلمات على مستوى النظرية. $(K_0,\0,\، ج_نص{sph},\، ج_نص{قرص},\، ج_نص{ذراع}، \، \، \، \لامبدا)$ – شامل، لا يوجد ضبط لكل مجرة. سيضيف التنقيح قيد الدراسة اعتمادًا على $\Sigma_d$.

7. الاتجاه. المدخلات → الباريونات → المجال الموجي → منحنى الدوران. بدون تغذية راجعة. منحنى الدوران هو تنبؤ وليس مطابقة.

المراجع. Lelli, F., McGaugh, S. S., Schombert, J. M. – SPARC: نماذج الكتلة ل 175 مجرات قرصية مع قياس ضوئي لسبيتزر ومنحنيات دوران دقيقة، AJ 152, 157 (2016). – فريمان، ك. س. – على أقراص المجرات الحلزونية ومجرات S0، ApJ 160, 811 (1970). – Hernquist, L. – نموذج تحليلي للمجرات الكروية والانتفاخات، ApJ 356، 359 (1990). – Broeils, A. H., Rhee, M.-H. – ملاحظات قصيرة 21-cm WSRT للمجرات الحلزونية وغير المنتظمة، A&A 324، 877 (1997). – McGaugh، S. S. – القانون الثالث للدوران المجرّي، المجرات 2، 601 (2014). – Bovy, J., Rix, H.-W. – قياس ديناميكي مباشر لملف كثافة سطح قرص مجرة درب التبانة، ApJ 779, 115 (2013). – Arnett, D. – Supernovae and Nucleosynthesis, Princeton (1996). – Dutertre, X. – نظرية النحل™: النمذجة المستندة إلى الموجة للجاذبية، الإصدار 2، BeeTheory.com (2023).

إضفاء الطابع الرسمي:حساب نظرية النحل على نطاق المجرة