BeeTheory – Fondasi – Catatan Teknis XII
Formalisasi:
Perhitungan Teori Lebah Berskala Galaksi
Catatan ini memformalkan kerangka kerja BeeTheory yang diterapkan pada galaksi piringan. Catatan ini merinci input pengamatan, dekomposisi geometris dari distribusi baryonik, persamaan integral yang mendefinisikan medan gelombang untuk setiap komponen, dan rantai operasi yang menghasilkan kurva rotasi yang diprediksi. Prosedurnya sangat searah: struktur baryonik yang diamati menentukan medan gelombang, yang menentukan kurva rotasi – tidak pernah sebaliknya.
1. Komputasi dalam satu diagram
Rantai searah
Fotometri yang diamati $\;\longrightarrow\;$ Dekomposisi baryonik $(\rho_\text{bar})$
$\big\panah bawah\$
Konvolusi medan gelombang $\;\longrightarrow\;$ Kerapatan gelombang $(\rho_\text{gelombang})$
$\big\downarrow$
Integrasi massa $\;\longrightarrow\;$ Massa gelombang tertutup $(M_\text{gelombang})$
$\big\downarrow$
Hubungan Newton $\;\longrightarrow\;$ Kurva rotasi yang diprediksi $(V_c)$
Tidak ada langkah yang terbalik. Kurva rotasi $V_c(R)$ tidak pernah digunakan sebagai input.
2. Masukan observasi
Untuk setiap galaksi, komputasi membutuhkan lima variabel yang telah dipublikasikan. Ini adalah satu-satunya besaran spesifik galaksi; yang lainnya dihitung dari besaran tersebut. Tidak ada pencocokan terhadap kurva rotasi yang dilakukan pada tahap ini.
| Simbol | Kuantitas | Sumber |
|---|---|---|
| $T$ | Tipe morfologi Hubble | Katalog (de Vaucouleurs et al. 1991, SPARC) |
| $ R_d $ | Panjang skala cakram bintang (kpc) | Fotometri Spitzer 3,6 µm (SPARC) |
| $\Sigma_d$ | Kecerahan permukaan disk tengah ($L_\odot/\text{pc}^2$) | Fotometri Spitzer 3,6 µm (SPARC) |
| $M_\text{HI}$ | Massa hidrogen atom total ($M_odot$) | Pengamatan radio 21 cm (SPARC) |
| $\Upsilon_\star$ | Rasio massa terhadap cahaya bintang pada 3,6 µm | Tetap universal: $ 0,5 \,M_\odot/L_\odot$ (McGaugh 2014) |
Untuk Bima Sakti, $R_d$, $Sigma_d$, dan $M_text{HI}$ digantikan dengan nilai analog yang ditentukan dari survei bintang internal (Bovy & Rix 2013) dan peta 21 cm. Vektor masukan lima kuantitas yang sama digunakan.
3. Dekomposisi baryonik – lima komponen geometris
Dari lima input pengamatan, massa baryonik dipartisi menjadi lima komponen geometris yang berbeda. Setiap komponen memiliki profil kerapatan dan skala karakteristiknya masing-masing.
3.1 Total massa bintang dan gas
$$M_\star \;=\; 2\pi\,R_d^2\,\Sigma_d\,\Upsilon_\star$$
$$M_\text{gas} \;=\; 1.33\,M_\text{HI} \qquad \text{(Koreksi He; Arnett 1996)}$$
3.2 Massa dan timbangan komponen
| Komponen | Massa | Skala | Aktivasi |
|---|---|---|---|
| Tonjolan | $M_b = 0.20\,M_\star$ | $r_b = \max(0.5\,R_d,\,0.3\text{ kpc})$ | Jika $T \leq 4$ |
| Disk tipis | $M_\text{tipis} = 0.75\,(M_\star – M_b)$ | $ R_d $ | Selalu |
| Disk tebal | $M_\text{tebal} = 0.25\,(M_\star – M_b)$ | $ 1,5 \, R_d $ | Selalu |
| Cincin gas | $M_\text{gas} = 1.33\,M_\text{HI}$ | R_g = 1,7 \,R_d $ (Broeils & Rhee 1997) | Selalu |
| Lengan spiral | $M_\text{lengan} = 0.10\,M_\text{tipis}$ (efektif) | $R_d$ (mengikuti disk tipis) | Selalu |
3.3 Profil kepadatan
Tonjolan (Hernquist 3D)
$$\rho_b(r) \;=\; \frac{M_b\,r_b}{2\pi\,r\,(r + r_b)^3}$$
Piringan bintang tipis dan tebal (eksponensial 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)}$$
Cincin gas (eksponensial 2D dengan lubang tengah)
$$\Sigma_\text{gas}(R) \;=\; \frac{M_\text{gas}}{2\pi\,R_g^2}\,\exp\!\left(-\frac{R_\text{lubang}}{R} – \frac{R}{R_g}\right), \quad R_\text{lubang} = 0,5\,R_g$$
Kelebihan lengan spiral (2D, mengikuti cakram tipis)
$$\Sigma_\text{lengan}(R) \;=\; 0.10\;\Sigma_\text{tipis}(R)$$
4. Kernel gelombang
Setiap elemen massa baryonik menghasilkan medan gelombang Teori Lebah. Medan pada titik $vec{r}$ yang dihasilkan oleh elemen sumber pada $vec{r},’$ yang dipisahkan oleh $D = |vec{r} – vec{r},’|$ diatur oleh kernel tipe Yukawa yang berasal dari fungsi gelombang teregulasi pada Catatan I:
Kernel gelombang BeeTheory
$$\mathcal{K}_i(D) \;=\; K_0\,\frac{(1 + \alpha_i\,D)\,e^{-\alpha_i\,D}}{D^2}, \qquad \alpha_i\;=\; \frac{1}{\ell_i}$$
Di sini $K_0$ adalah amplitudo massa gelombang universal (angka tanpa dimensi tunggal) dan $\ell_i$ adalah panjang koherensi komponen $i$. Kernel mengkodekan perilaku quasi-Newtonian $1/D^2$ pada pemisahan pendek, dimodulasi oleh cutoff eksponensial pada skala di luar $\ell_i$. Bentuk $(1 + \alpha D)\,e^{-\alpha D}$ memastikan kontinuitas dan total massa tertutup yang terbatas pada tak terhingga.
4.1 Panjang koherensi komponen
Panjang koherensi setiap komponen diatur oleh skala geometris alami, dikalikan dengan konstanta tak berdimensi yang spesifik untuk dimensinya:
| Komponen | Panjang koherensi | Konstanta geometris |
|---|---|---|
| Tonjolan (bola 3D) | $\ell_b = c_\text{sph}\,r_b$ | $c_\text{sph}$ |
| Disk tipis (2D) | $\ell_\text{tipis} = c_\text{disk}\,R_d$ | $c_\text{disk}$ |
| Disk tebal (2D) | $\ell_\text{tebal} = c_\text{disk}\,(1.5\,R_d)$ | $c_\text{disk}$ |
| Cincin gas (2D) | $\ell_\text{gas} = c_\text{disk}\,R_g$ | $c_\text{disk}$ |
| Lengan spiral (2D, terkonsentrasi secara azimuthal) | $\ell_\text{lengan} = c_\text{lengan}\,R_d$ | $c_\text{arm}$ |
Tiga konstanta geometris $(c_\text{sph},\,c_\text{disk},\,c_\text{arm})$ bersifat universal – tidak berbeda dari satu galaksi ke galaksi lainnya. Bersama-sama dengan amplitudo massa-gelombang global $K_0$ dan kopling medan-gelombang $\lambda$, mereka membentuk satu set lengkap parameter tingkat teori.
5. Konvolusi medan gelombang – persamaan integral per komponen
Kerapatan medan gelombang pada posisi $\vec{r}$ adalah konvolusi dari distribusi sumber baryonik dengan kernel gelombang. Untuk sistem simetris galaksi (simetris aksial, pendekatan monopolar), setiap komponen baryonik berkontribusi secara aditif:
Kepadatan medan gelombang total pada radius $r$
$$\rho_\text{gelombang}(r) \;=\; \lambda \;\jumlah_{i \in \{\text{tipis, tebal, gas, lengan, tonjolan}\}} \rho_\text{wave}^{(i)}(r)$$
Kelima integral tersebut dituliskan di bawah ini, satu per komponen. Setiap integral mengubah distribusi massa baryonik menjadi distribusi massa medan gelombang pada titik spasial yang sama.
5.1 Tonjolan – Integrasi cangkang 3D
$$\rho_\text{wave}^{(b)}(r) \;=\; \int_0^{r_\text{max}} \rho_b(r’)\;\mathcal{K}_b\!\left(\sqrt{r^2 + r’^2}\right)\;4\pi r’^2\,dr’$$
Integrasi dilakukan pada cangkang bola konsentris dengan jari-jari $r’$. Titik medan pada radius $r$ dari pusat melihat setiap cangkang pada pemisahan efektif $D = \sqrt{r^2 + r’^2}$ dalam pendekatan monopolar. Integrasi meluas hingga $r_\text{max} = 6\,r_b$, di luar itu kerapatan tonjolan secara numerik dapat diabaikan.
5.2 Disk tipis – Integrasi cincin 2D
$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{tipis}(R’)\;\mathcal{K}_\text{tipis}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
Piringan diuraikan menjadi cincin konsentris dengan jari-jari $R’$ dan lebar tak terhingga $dR’$, masing-masing membawa massa permukaan $\Sigma_\text{tipis}(R’)\,2\pi R’\,dR’$. Pendekatan monopolar yang sama berlaku: medan gelombang pada radius $r$ dari pusat menerima kontribusi dari setiap cincin pada pemisahan efektif $D = \sqrt{r^2 + R’^2}$. Rentang integrasi adalah $R_\text{max} = 8\,R_d$.
5.3 Disk tebal – Integrasi cincin 2D
$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{tebal}(R’)\;\mathcal{K}_\text{tebal}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
Identik dengan integrasi cakram tipis, dengan $\Sigma_\text{tebal}(R’)$ sebagai kerapatan sumber dan parameter kernel $\alpha_\text{tebal} = 1/(c_\text{disk}\, \cdot 1.5\, R_d)$. Jangkauan radial yang lebih luas dari cakram tebal menghasilkan rentang koherensi gelombang yang sedikit lebih luas.
5.4 Cincin gas – Integrasi cincin 2D dengan penipisan pusat
$$\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’$$
Distribusi gas memiliki lubang pusat, yang ditangkap oleh cutoff radial pada $R_\text{hole} = 0.5\,R_g$ di batas bawah integrasi. Di luar batas ini, gas meluas lebih jauh dari piringan bintang; ini tercermin dalam skala karakteristik yang lebih besar $R_g = 1.7\,R_d\,R_g\, yang menjadi masukan bagi panjang koherensi $\ell_\text{gas} = c_\text{disk}\,R_g\.
5.5 Kelebihan lengan spiral – Integrasi cincin 2D dengan amplitudo yang dikurangi
$$\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’$$
Lengan spiral diperlakukan sebagai peningkatan rata-rata secara aksial dari densitas permukaan cakram tipis pada tingkat $10\%$, dengan panjang koherensinya sendiri $\ell_\text{lengan} = c_\text{lengan}\,R_d$. Oleh karena itu, kernel ini lebih sempit daripada kernel cakram tipis, yang mencerminkan konsentrasi azimuthal dari struktur spiral.
6. Massa gelombang tertutup dan kurva rotasi yang diprediksi
Setelah kerapatan medan gelombang total $\rho_\text{wave}(r)$ diketahui, massa medan gelombang yang tertutup di dalam bola berjari-jari $R$ diperoleh dengan integrasi radial:
Massa medan gelombang tertutup
$$M_\text{gelombang}(R) \;=\; \int_0^{R} 4\pi\,r^2\,\rho_\text{gelombang}(r)\,dr$$
Kecepatan melingkar yang diprediksi pada radius $R$ kemudian mengikuti hubungan Newtonian, menggabungkan kontribusi baryonik dan medan gelombang dalam kuadratur:
Prediksi kecepatan melingkar
$$V_c^2(R) \;=\; V_\text{bar}^2(R) \;+\; \frac{G\,M_\text{gelombang}(R)}{R}$$
Kecepatan baryonik $V_\text{bar}(R)$ adalah jumlah kuadrat kontribusi dari empat komponen seperti piringan (rumus Freeman 1970 untuk setiap profil eksponensial) dan tonjolan (rumus massa terlampir Hernquist):
$$V_\text{bar}^2(R) \;=\; V_\text{bulge}^2 + V_\text{thin}^2 + V_\text{thick}^2 + V_\text{gas}^2$$
di mana setiap $V_i(R)$ adalah kecepatan melingkar Newtonian standar dari distribusi massa yang sesuai.
7. Parameter tingkat teori
Kerangka kerja BeeTheory yang lengkap, seperti yang diterapkan pada galaksi, berisi lima parameter pada tingkat teori. Parameter-parameter ini bersifat universal: tidak berbeda antara satu galaksi dengan galaksi lainnya.
| Simbol | Arti | Peran |
|---|---|---|
| $K_0$ | Amplitudo massa gelombang | Mengatur skala tanpa dimensi dari kernel gelombang |
| $c_\text{sph}$ | Konstanta geometris 3D | Rasio $\ell/r_\text{skala}$ untuk sumber bola (tonjolan) |
| $c_\text{disk}$ | Konstanta geometris 2D | Rasio $\ell/R_\text{scale}$ untuk sumber disk dan ring |
| $c_\text{arm}$ | Konstanta geometris spiral | Rasio $\ell/R_d$ untuk kelebihan lengan yang terkonsentrasi secara azimuthal |
| $\lambda$ | Kopling medan gelombang global | Mengukur kepadatan medan gelombang total |
Universalitas parameter
Kelima parameter tersebut bersifat global. Nilai numerik yang sama berlaku untuk Bimasakti, bintang katai tak beraturan, dan spiral masif. Informasi spesifik galaksi hanya masuk melalui lima masukan observasi $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\bintang)$. Model ini tidak mengandung parameter yang dapat disetel per galaksi.
8. Sifat komputasi yang searah
Rantai terbuka – tidak ada umpan balik
Seluruh komputasi mengalir dari input ke output, dalam satu arah. Pengamatan fotometrik dan 21 cm menentukan dekomposisi baryonik. Penguraian baryonik menentukan kerapatan medan gelombang. Kerapatan medan gelombang menentukan massa gelombang tertutup. Massa gelombang tertutup menentukan kurva rotasi yang diprediksi. Kurva rotasi tidak mempengaruhi langkah sebelumnya dalam perhitungan.
Kesatuan arah ini memiliki tiga konsekuensi penting.
(a) Setelah lima parameter tingkat teori ditetapkan, kurva rotasi adalah prediksi yang ketat, bukan pencocokan. Perbandingan dengan kurva rotasi yang diamati adalah pengujian, bukan kalibrasi.
(b) Model ini tidak memiliki mekanisme untuk penyesuaian galaksi per galaksi. Setiap modifikasi prediksi kurva rotasi harus berasal dari modifikasi vektor masukan $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$ atau dari perubahan pada parameter tingkat teori universal $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$.
(c) Mengkalibrasi $lambda$ pada galaksi referensi tidak sama dengan menyesuaikannya dengan kurva rotasi galaksi tersebut. Kalibrasi menentukan satu angka global; kurva rotasi pada semua jari-jari galaksi referensi, dan kurva rotasi semua galaksi lainnya, kemudian merupakan prediksi yang tepat dari kerangka yang telah dikalibrasi.
9. Peran densitas permukaan pusat (Revisi Catatan XI)
Diagnostik Catatan XI mengidentifikasi bahwa kesalahan prediksi residual berkorelasi kuat dengan kerapatan permukaan baryonik pusat $\Sigma_d$, terlepas dari panjang skala piringan $R_d$. Formalisasi yang disajikan di atas adalah versi model sebelum temuan ini digabungkan – model ini hanya menggunakan $R_d$ dalam ekspresi panjang koherensi $\ell_i = c_i\,R_d$.
Di mana penyempurnaan akan masuk
Dalam model yang telah disempurnakan, panjang koherensi $\ell_i$ akan bergantung pada $\R_d$ dan $\Sigma_d$, menggantikan hubungan linier yang ketat $\ell_i = c_i\,R_d$ dengan sebuah fungsi $\ell_i = c_i\,R_d\,\phi(\Sigma_d/\Sigma_\text{ref})$ yang menyerap sisa yang teridentifikasi dalam Catatan XI. Bentuk fungsional dari $phi$ dan parameternya akan ditentukan dalam catatan berikutnya, pertama pada set kalibrasi 22 galaksi, kemudian divalidasi dengan prediksi buta pada sampel SPARC yang tersisa.
Struktur searah dari komputasi dipertahankan oleh penyempurnaan ini: $\Sigma_d$ adalah input observasi, panjang koherensi yang dimodifikasi dimasukkan ke dalam integral konvolusi yang sama, dan kurva rotasi muncul seperti sebelumnya. Hanya satu hubungan operasional yang ditambahkan – ketergantungan $\ell_i$ pada pengamatan kedua.
10. Ringkasan metodologi
1. Masukan. Lima pengamatan per galaksi: tipe Hubble $T$, skala cakram $R_d$, kecerlangan permukaan $\Sigma_d$, massa HI $M_\text{HI}$, dan rasio massa terhadap cahaya bintang universal $\Upsilon_\bintang$.
2. Dekomposisi baryonik. Lima komponen: tonjolan (jika $T \leq 4$), piringan tipis, piringan tebal, cincin gas, kelebihan lengan spiral. Masing-masing membawa profil densitas analitik.
3. Kernel gelombang. Bentuk tipe Yukawa universal $\mathcal{K}_i(D) = K_0\,(1 + \alpha_i D)\, e^{-\alpha_i D}/D^2$ dengan panjang koherensi $\ell_i = c_i\, R_\text{skala}$ yang ditentukan oleh luas geometris dari setiap komponen.
4. Konvolusi. Setiap komponen menghasilkan kerapatan medan gelombang melalui integral satu dimensi pada cincin (komponen 2D) atau cangkang (tonjolan 3D). Kerapatan medan gelombang total adalah jumlah dari lima komponen, yang diskalakan dengan kopling global $\lambda$.
5. Keluaran. Massa gelombang terlampir $M_\text{wave}(R)$ diintegrasikan dan digabungkan dengan kecepatan baryonik $V_\text{bar}(R)$ untuk menghasilkan kurva rotasi yang diprediksi $V_c(R)$.
6. Parameter tingkat teori. $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$ – universal, tidak ada penyetelan per galaksi. Penyempurnaan yang sedang diteliti akan menambah ketergantungan pada $\Sigma_d$.
7. Arah. Input → baryon → medan gelombang → kurva rotasi. Tidak ada umpan balik. Kurva rotasi adalah prediksi, bukan kecocokan.
Referensi. Lelli, F., McGaugh, S. S., Schombert, J. M. – SPARC: Model Massa untuk 175 Galaksi Cakram dengan Fotometri Spitzer dan Kurva Rotasi yang Akurat, AJ 152, 157 (2016). – Freeman, K. C. – Pada piringan galaksi spiral dan galaksi S0, ApJ 160, 811 (1970). – Hernquist, L. – Model analitik untuk galaksi bola dan tonjolan, ApJ 356, 359 (1990). – Broeils, AH, Rhee, M.-H. – Pengamatan WSRT pendek 21 cm pada galaksi spiral dan tak beraturan, A&A 324, 877 (1997). – McGaugh, S. S. – Hukum ketiga rotasi galaksi, Galaxies 2, 601 (2014). – Bovy, J., Rix, H.-W. – Pengukuran dinamik langsung dari profil kerapatan permukaan piringan Bima Sakti, ApJ 779, 115 (2013). – Arnett, D. – Supernova dan Nukleosintesis, Princeton (1996). – Dutertre, X. – Teori Lebah ™: Pemodelan Gravitasi Berbasis Gelombang, v2, BeeTheory.com (2023).
BeeTheory.com – Gravitasi kuantum berbasis gelombang – Metodologi galaksi – © Technoplane S.A.S. 2026