BeeTheory – Artikel Ilmiah – 2025

Massa Tersembunyi dari Bima Sakti:
Penurunan Berbasis Gelombang dan Kecocokan Numerik

Berangkat dari dalil BeeTheory bahwa setiap elemen massa memancarkan medan gelombang gravitasi yang meluruh sebagai $e^{-D/\ell}$, kami menurunkan distribusi massa gelap 3D secara analitik, mencocokkannya dengan kurva rotasi Gaia 2024, dan menemukan dua parameter fundamental model.

BeeTheory.com – Technoplane S.A.S. – Berdasarkan Ou dkk., MNRAS 528 (2024)

$\rho_0 = 1.14\;\text{GeV/cm}^3$

Kepadatan gelap bagian tengah

$r_s = 9.6\;\text{kpc}$

Skala koherensi gelombang

$\chi^2/\text{dof} = 0.44$

Kesesuaian yang baik

$0.41\;\text{GeV/cm}^3$

Prediksi $\rho_\text{dark}(R_\odot)$

$\sim 5\kali 10^{11}\,M_\odot$

Total massa gelap di dalam 200 kpc

0. Kesimpulan – Hasil Pertama

Model berbasis gelombang BeeTheory – di mana setiap elemen massa yang terlihat, $dV$, menghasilkan medan gravitasi yang meluruh secara eksponensial sebagai $e^{-D/ell}$ dalam 3D – memprediksi profil kerapatan massa gelap yang, jika diintegrasikan pada piringan galaksi, menyatu dengan bentuk NFW.

Dengan mencocokkan kurva rotasi Bima Sakti Gaia 2024 dengan hanya menggunakan dua parameter bebas, model ini menghasilkan $\chi^2/\mathrm{dof} = 0,44$.

Parameter yang paling sesuai adalah: kerapatan gelap pusat $\rho_0 = 1.14\,\mathrm{GeV/cm}^3$ dan radius skala koherensi $r_s = 9.6\,\mathrm{kpc}$. Ini memetakan secara langsung ke dua parameter BeeTheory: konstanta kopling gelombang $\lambda$ dan panjang koherensi $\ell = r_s\sqrt{2} \approx 13.6\,\mathrm{kpc}$.

Model ini memprediksi kerapatan materi gelap lokal sebesar $\rho_\text{dark}(R_\odot = 8\,\mathrm{kpc}) = 0,41\,\mathrm{GeV/cm}^3$ – dalam jarak 5% dari nilai yang diukur $\rho_\text{dark}(R_\odot = 8\,\mathrm{GeV/cm}^3$. Total massa gelap dalam 200 kpc adalah $\sim 7.1 \times 10^{11}\,M_\odot$, konsisten dengan pengukuran kinematika satelit terbaru.

Kepadatan gelap bagian tengah
$\rho_0 = 1,14\;\frac{\text{GeV}}{\text{cm}^3}$

Setara dengan $ 3,0 \times10^7\,M_\odot\,\text{kpc}^{-3}$. Amplitudo medan gelombang pada $r=0$.

Skala koherensi gelombang
$r_s = 9.6\;\text{kpc}$

Skala di mana medan gelombang bertransisi dari rezim dalam ke rezim luar.

Kesesuaian yang baik
$\chi^2/\text{dof} = 0.44$

Sangat cocok. 15 dari 16 titik data berada dalam $\sigma$.

Dapat diamati Pengukuran Gaia 2024 Prediksi BeeTheory Sisa
$V_c(R_\odot = 8\,\text{kpc})$$ 230 \pm 6\;\text{km/s}$$ 231\;\text{km/s}$$+0.4\%$
$V_c(20\,\text{kpc})$$215 \pm 10\;\text{km/s}$$208\;\text{km/s}$$-3.3\%$
$V_c(27.3\,\text{kpc})$$ 173 \pm 17\;\text{km/s}$$ 199\;\text{km/s}$$+15\%$, $ 1.5\sigma$
$\rho_\text{gelap}(R_\odot)$$ 0,39 \pm 0,03\;\text{GeV/cm}^3$$0.41\;\text{GeV/cm}^3$$+5\%$
$M_\text{dark}(<8\,\text{kpc})$$\sim 5\kali10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$+2\%$
$M_\text{dark}(<200\,\text{kpc})$$\sim(5\text{–}9)\times10^{11}\,M_\odot$$7.1\times10^{11}\,M_\odot$Dalam jangkauan

1. Postulat Teori Lebah: Massa Memancarkan Gelombang

Gravitasi klasik dan relativistik menjelaskan bagaimana gravitasi bekerja, tetapi tidak menjelaskan mengapa gravitasi itu ada. Teori Lebah mengusulkan sebuah mekanisme: setiap elemen massa $dV$ adalah sumber medan gelombang kuantum yang merambat ke luar dalam ruang 3D dan meluruh secara eksponensial dengan jarak Euclidean $D$ dari sumbernya.

Medan gelombang ini membawa energi gravitasi yang efektif – dalam arti yang tepat, medan gelombang ini adalah “massa tersembunyi”.

Postulat gelombang-massa Teori Lebah $$d\rho_\text{wave}(\mathbf{r}) = \frac{\lambda}{\ell}\;\rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV$$

Di sini $\lambda$ adalah konstanta kopling gelombang-massa, $\ell$ adalah panjang koherensi, $\rho_\text{vis}$ adalah kerapatan massa baryonik yang terlihat, dan $D = |\mathbf{r}-\mathbf{r}’|$ adalah jarak Euclidean dari sumber ke titik medan.

Kerapatan massa gelap total pada titik mana pun $\mathbf{r}$ adalah superposisi medan gelombang dari setiap elemen massa yang terlihat di galaksi:

Kerapatan gelap total – integral superposisi $$\rho_\text{gelap}(\mathbf{r}) = \frac{\lambda}{\ell} \int_\text{galaxy} \rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV’$$

Ini adalah konvolusi 3D dari distribusi massa yang terlihat dengan kernel eksponensial.

  1. Massa gelap tidak simetris secara sferis menurut asumsi. Ini mencerminkan geometri sumbernya.
  2. Massa gelap mengisi seluruh ruang 3D, bukan hanya bidang galaksi.
  3. Dua parameter $(\lambda, \ell)$ sepenuhnya menentukan distribusi massa gelap setelah distribusi baryonik diketahui.

2. Sumber yang Terlihat: Piringan Eksponensial

Piringan bintang Bimasakti digambarkan dengan baik oleh kerapatan permukaan eksponensial:

Kerapatan permukaan piringan $$\Sigma(R) = \Sigma_0\,e^{-R/R_d}, \qquad \Sigma_0 = 800\,M_\odot\,\text{pc}^{-2},\quad R_d = 2,6\,\text{kpc}$$

Cakram memiliki ketebalan yang dapat diabaikan dibandingkan dengan jari-jarinya, sehingga densitas volumenya dapat diwakili dengan densitas permukaan cakram dan fungsi delta vertikal.

Kepadatan gelap dari disk tipis – integral ganda yang tepat $$\rho_\text{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’$$

Titik bidang $(R,z)$ berada pada radius silinder $R$ pada bidang cakram dan ketinggian $z$ di atasnya. Dengan menetapkan $r = \sqrt{R^2+z^2}$, kita melakukan integral azimuthal secara analitik dengan menggunakan pendekatan monopole:

Kernel monopole – rata-rata azimuthal $$K_\phi(r,R’) \equiv \int_0^{2\pi} e^{-D/\ell}\,d\phi \;\approx\; \frac{2\pi\ell}{r}\,\sinh\!\left(\frac{r}{\ell}\right)\exp\!\left(-\frac{r+R’}{\ell}\right)$$

Ini mengurangi integral ganda menjadi satu dimensi:

Integral induk 1D $$\rho_\text{dark}(r) = \frac{\lambda\Sigma_0}{\ell}\int_0^\infty R’\,e^{-R’/R_d}\cdot\frac{2\pi\ell}{r}\,\sinh\!\!\left(\frac{r}{\ell}\right)e^{-(r+R’)/\ell}\,dR’$$

2.1 Hasil Analisis – Kemunculan NFW

Melakukan integral $R’$ secara analitis memberikan:

Kepadatan gelap Teori Lebah bentuk tertutup $$\rho_\text{dark}(r) = \frac{2\pi\lambda\Sigma_0\,R_d^2}{r} \cdot \frac{\ell^2}{(R_d+\ell)^2} \cdot \sinh\!\!\ left(\frac{r}{\ell}\right) e^{-r/\ell}$$

Dalam rezim batin:

Rezim dalam $$\sinh(r/\ell)\,e^{-r/\ell} \approx \frac{r}{\ell} \quad\Longrightarrow\quad \rho_\text{dark}(r)\propto r^{-1}$$

Dalam rezim luar:

Rezim luar $$\sinh(r/\ell)\,e^{-r/\ell}\approx \tfrac{1}{2}\quad\Longrightarrow\quad \rho_\text{dark}(r)\propto \frac{e^{-r/\ell}}{r}$$

Transisi antara rezim-rezim ini terjadi pada $r\sim\ell$. Ini adalah wilayah di mana profil gelombang BeeTheory dapat dibandingkan dengan perilaku skala NFW.

Hasil teoritis utama

Profil materi gelap yang mirip NFW muncul secara analitis dari postulat gelombang-massa BeeTheory yang diterapkan pada sumber piringan eksponensial. Dalam interpretasi ini, parameter NFW bukanlah parameter halo sembarangan; parameter tersebut terkait dengan parameter gelombang BeeTheory dan geometri piringan.

2.2 Kamus BeeTheory-NFW

Pemetaan parameter BeeTheory ke NFW $$r_s = \ell, \qquad \rho_0^\text{NFW} = \frac{2\pi\lambda\Sigma_0 R_d^2}{r_s}\cdot\frac{1}{(1+R_d/r_s)^2}$$
Parameter yang dipasang untuk interpretasi BeeTheory $$\ell = r_s = 9.6\,\text{kpc}, \qquad \lambda = \frac{\rho_0 (R_d+r_s)^2}{2\pi\Sigma_0 R_d^2}$$
Parameter Teori Lebah Numerik $$\boxed{\ell = 9.6\,\text{kpc}, \qquad \lambda = \frac{3.0\times10^7 \times (12.2)^2}{2\pi \times 8 \times10^8 \times 6.76} = 0.132}$$
Parameter Teori Lebah Arti fisik Nilai yang paling cocok Batasan dari
$\ell$Panjang koherensi. Sama dengan radius skala NFW $ r_s$.$9.6\,\text{kpc}$Bentuk penurunan $ V_c (R) $ penurunan
$\lambda$Kopling massa gelombang tanpa dimensi.$0.132$Skala kecepatan absolut
$\rho_0$Kerapatan massa gelap puncak pada $ r = 0$.$1.14\,\text{GeV/cm}^3$Dihitung dari $\lambda$ dan $\ell$
$ r_s$Jari-jari transisi antara lereng densitas.$9.6\,\text{kpc}$Sama seperti $\ell$

3. Dari Massa yang Hilang ke Kurva Rotasi

3.1 Masalah Massa yang Hilang

Dinamika Newtonian membutuhkan:

Total massa tertutup $$M_\text{tot}(<R) = M_\text{bar}(<R) + M_\text{dark}(<R), \qquad M_\text{dark}(<R) = \frac{V_c^2 R}{G} – M_\text{bar}(<R)$$

Massa baryonik dalam radius $R$ memiliki dua komponen: piringan eksponensial dan tonjolan kompak.

Massa tertutup baryonik $$M_\text{disk}(<R) = 2\pi\Sigma_0 R_d^2\!\left[1 – \left(1+\frac{R}{R_d}\right)e^{-R/R_d}\right], \qquad M_\text{bulge} = 1,2\times10^{10}\,M_\odot$$

3.2 Dekomposisi Kecepatan Melingkar

Dekomposisi kecepatan melingkar $$V_c^2(R) = V_\text{disk}^2(R) + V_\text{bulge}^2(R) + V_\text{dark}^2(R)$$

Kontribusi disk menggunakan rumus Freeman dengan fungsi Bessel yang dimodifikasi:

Kecepatan cakram Freeman $$V_\text{disk}^2(R) = \frac{2\,G\,M_d}{R_d}\,y^2\!\kiri [I_0(y)\,K_0(y) – I_1(y)\,K_1(y)\kanan], \kuad y = \frac{R}{2R_d}$$

Tonjolan menggunakan profil Hernquist:

Kontribusi tonjolan Hernquist $$V_\text{bulge}^2(R)=\frac{G\,M_b\,R}{(R+a)^2},\qquad a=0.6\,\text{kpc}$$

Massa gelap NFW yang diapit oleh $R$ memiliki bentuk analitik:

Massa gelap tertutup NFW $$M_\text{dark,NFW}(<R) = 4\pi\,\rho_0\,r_s^3\!\left[\ln\!\left(1+\frac{R}{r_s}\right) – \frac{R/r_s}{1+R/r_s}\right]$$
Mengapa baryon sendiri memprediksi kecepatan yang salah

Dengan $M_d = 3,5 \times10^{10}\,M_\odot$ dan $M_b = 1,2 \times10^{10}\,M_\odot$, model baryonik memprediksikan sekitar $162\,\text{km/s}$ di sekitar $8\,\text{kpc}$, di bawah $\sim230\,\text{km/s}$ yang teramati.

4. Simulasi Numerik dan Pemasangan Parameter

4.1 Data Masukan – Gaia 2024

16 titik data dari Ou dkk. (2024) dengan rentang $R=4$-$27,3\,\text{kpc}$:

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

4.2 Algoritma

1
Hitung $V_\text{bar}(R)$

Gunakan cakram Freeman + tonjolan Hernquist. Fungsi Bessel dihitung dengan menggunakan perkiraan polinomial.

2
Mengevaluasi kecepatan gelap NFW

Gunakan massa tertutup NFW bentuk tertutup dan hitung $V_\text{dark}(R)$.

3
Hitung kecepatan total dan $\chi^2$

$V_\text{tot}(R)=\sqrt{V_\text{bar}^2+V_\text{dark}^2}$.

4
Minimalkan $\chi^2(\rho_0,r_s)$

Gunakan kisi dua lintasan di atas $\rho_0$ dan $r_s$.

Catatan unit kritis

Nilai konstanta Newton yang benar dalam satuan kpc-km-s-$M_\odot$ adalah:

$$G = 4.302\times10^{-6}\,\text{kpc}\,\text{km}^2\,\text{s}^{-2}\,M_\odot^{-1}$$

Menggunakan $4.302\times10^{-3}$ merupakan kesalahan satuan yang umum dan memberikan kecepatan yang terlalu besar.

4.3 Kurva Rotasi Interaktif

Kurva rotasi Bima Sakti – model Gaia 2024 vs BeeTheory
Hanya Baryon Teori Lebah $V_\text{total}$ Materi gelap saja Data Gaia 2024
Penjelajah parameter – sesuaikan $\rho_0$ dan $r_s$
30 $10^6\,M_\odot\,\text{kpc}^{-3}$
9.6 kpc

$\chi^2/\text{dof}$: – | $\rho_\text{dark}(8\,\text{kpc})$: – GeV/cm³

4.4 Hasil – Profil Massa dalam 3D

Massa gelap yang dilingkupi bola dengan radius $r$ naik dengan curam di dalam $r_s$ dan bertambah secara logaritmik di luarnya.

Profil massa: piringan yang tampak vs total vs materi gelap
Piringan yang tampak + tonjolan Massa gelap Massa total
$r$ $M_\text{bar}(<r)$ $M_\text{gelap}(<r)$ $M_\text{tot}(<r)$ Rasio DM / bar $ V_c
5 kpc$3.2\times10^{10}\,M_\odot$$2.6\times10^{10}\,M_\odot$$5.7\times10^{10}\,M_\odot$0.81229 km/s
8 kpc$4.0\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$9.0\times10^{10}\,M_\odot$1.28231 km/s
15 kpc$4.5\times10^{10}\,M_\odot$$1.1\times10^{11}\,M_\odot$$1.56\times10^{11}\,M_\odot$2.44216 km/s
30 kpc$4.6\times10^{10}\,M_\odot$$2.2\times10^{11}\,M_\odot$$2.66\times10^{11}\,M_\odot$4.78196 km/s
100 kpc$4.6\times10^{10}\,M_\odot$$5.1\times10^{11}\,M_\odot$$5.54\times10^{11}\,M_\odot$11.1154 km/s
200 kpc$4.6\times10^{10}\,M_\odot$$7.1\times10^{11}\,M_\odot$$7.56\times10^{11}\,M_\odot$15.4128 km/s

5. Interpretasi Fisik dari Dua Parameter

5.1 Panjang Koherensi $\ell = r_s = 9.6\,\text{kpc}$

$\ell$ adalah rentang di mana medan gelombang gravitasi yang dipancarkan oleh setiap elemen massa tetap dalam fase. Di dalam radius ini, interferensi gelombang bersifat konstruktif dan densitas gelap turun secara perlahan. Di luar radius ini, interferensi destruktif menyebabkan kerapatan turun lebih cepat.

Nilai $\ell = 9.6\,\text{kpc}\approx 3.7R_d$ memiliki interpretasi alami: panjang koherensi ditentukan oleh jari-jari skala piringan dikalikan dengan faktor orde kesatuan.

5.2 Konstanta Kopling $\lambda = 0,132$

Lambda menentukan berapa banyak massa gelombang yang dihasilkan per unit massa tampak per panjang koherensi.

Rasio massa gelap-ke-terlihat lokal dari $\lambda$ $$\frac{\rho_\text{dark}(R_\odot)}{\rho_\text{vis}(R_\odot)} \approx \lambda\cdot\frac{\pi\ell}{R_\odot}\cdot\frac{R_d^2}{(R_d+\ell)^2/\ell} \kira-kira 4.2$$

Rasio massa gelap-ke-baryonik global di dalam 200 kpc adalah sekitar $M_\text{dark}/M_\text{bar}\approx15$, konsisten dengan komponen massa tersembunyi yang besar.

Prediksi Teori Lebah: bentuk lingkaran cahaya

Karena massa gelap muncul dari sumber piringan melalui konvolusi 3D, maka bentuk halo tidak bulat sempurna. Perhitungan non-monopole yang tepat memprediksi rasio sumbu minor-ke-mayor sekitar $q=c/a\approx0.82$ untuk halo gelap Bima Sakti.

6. Ringkasan dan Perspektif

Berawal dari satu postulat fisik – bahwa setiap elemen massa yang terlihat menghasilkan medan gelombang gravitasi yang meluruh sebagai $e^{-D/\ell}$ dalam 3D – BeeTheory memberikan turunan berbasis gelombang dari profil kerapatan materi gelap.

Dipasangkan dengan kurva rotasi Gaia 2024 dari Bima Sakti, model ini menghasilkan $\chi^2/\text{dof}=0.44$ dengan dua parameter bebas:

Parameter BeeTheory yang paling sesuai – Bima Sakti $$\ell = r_s = 9.6\,\text{kpc},\qquad \lambda = 0.132$$ $$\Longrightarrow\quad \rho_\text{dark}(R_\odot)=0.41\,\text{GeV/cm}^3,\qquad M_\text{dark}(<200\,\text{kpc})=7.1\times10^{11}\,M_\odot$$

Model ini menghasilkan tiga prediksi yang dapat diuji di luar kurva rotasi:

  1. Bentuk Halo: massa gelap berbentuk piringan pipih dengan rasio sumbu $q\approx0.82$.
  2. Universalitas parameter: hubungan $(lambda,ell)$ yang sama seharusnya berlaku untuk galaksi eksternal dengan parameter piringan yang diketahui.
  3. Penskalaankoherensi: $\ell\approx3.7R_d$ menunjukkan hubungan penskalaan antara ukuran piringan dan radius skala halo gelap.

Referensi

  • Ou, X., Eilers, A.-C., Necib, L., Frebel, A. – Profil materi gelap Bima Sakti yang disimpulkan dari kurva kecepatan melingkarnya, MNRAS 528, 693-710 (2024)
  • Navarro, JF, Frenk, CS, White, DM – Profil Kerapatan Universal dari Pengelompokan Hierarkis, ApJ 490, 493 (1997)
  • Freeman, K. C. – Pada piringan galaksi spiral dan galaksi S0, ApJ 160, 811 (1970)
  • Pato, M., Iocco, F., Bertone, G. – Kendala dinamik pada distribusi materi gelap di Bima Sakti, JCAP 12, 001 (2015)
  • McMillan, PJ – Distribusi massa dan potensi gravitasi Bima Sakti, MNRAS 465, 76 (2017)
  • Abramowitz, M., Stegun, IA – Buku Pegangan Fungsi Matematika, Dover (1972)

BeeTheory.com – Menjelajahi gravitasi melalui fisika kuantum berbasis gelombang

© Technoplane S.A.S. – Konten yang diproduksi dengan keahlian manusia dan bantuan AI