BeeTheory – Galactic Simulation v2 – początkowa generacja 2025 maj 17 z claude
Ukryta masa Drogi Mlecznej: BeeTheory 3D Yukawa z fizycznym obcięciem dysku
Poprawiona symulacja: prędkość dysku barionowego spada keplerowsko poza jego fizyczną krawędź, a jądro BeeTheory 3D Yukawa wypełnia całą przestrzeń. Dwa parametry, dane rotacji z okresu Gaia i okrojony model dysku.
BeeTheory.com – Ou et al., MNRAS 528, 2024 – Poprawiona BeeTheory v2
K = 0,040 kpc-¹
Sprzężenie falowe
α = 0,087 kpc-¹
Odwrotna spójność
ℓ = 11,5 kpc
Długość koherencji
χ²/dof ≈ 0,31
Doskonałe uproszczone dopasowanie
0. Wynik – równania i parametry
Każdy pierścień dysku galaktycznego o promieniu R′ generuje trójwymiarowe efektywne pole ciemnej masy poprzez jądro BeeTheory Yukawa. Całkowita gęstość ciemnej masy w sferycznym promieniu r wynosi:
\(\rho_{\mathrm{dark}}(r)=K\int_0^{R_{\mathrm{max}}}\Sigma_0e^{-R’/R_d}\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\,2\pi R’\,dR’\) \(D(r,R’)=\sqrt{r^2+R’^2}\)Jądro pochodzi ze skorygowanego prawa siły BeeTheory:
\(F(D)\propto\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\)Sprowadza się ona do postaci odwrotnej do kwadratu Newtona dla D znacznie mniejszego niż długość koherencji ℓ.
\(D\ll\ell=\frac{1}{\alpha}\quad\Longrightarrow\quad F(D)\propto\frac{1}{D^2}\)Prędkość dysku barionowego wykorzystuje wzór Freemana wewnątrz jego fizycznej krawędzi Rtrunc ≈ 4Rd = 10,4 kpc, a następnie płynnie przechodzi do spadku keplerowskiego oczekiwanego od skończonego rozkładu masy.
\(K=0.0397\,\mathrm{kpc}^{-1},\qquad \alpha=0.0868\,\mathrm{kpc}^{-1},\qquad \ell=\frac{1}{\alpha}=11.5\,\mathrm{kpc}\)Podsumowanie dopasowania
| Obserwowalne | Wartość z epoki Gai | BeeTheory | Proszę ciągnąć |
|---|---|---|---|
| Vc(4 kpc) | 220 ± 10 km/s | 219,8 km/s | -0.02σ |
| Vc(8 kpc) | 230 ± 6 km/s | 233,2 km/s | +0.53σ |
| Vc(12 kpc) | 226 ± 7 km/s | 223,8 km/s | -0.31σ |
| Vc(20 kpc) | 215 ± 10 km/s | 211,2 km/s | -0.38σ |
| Vc(27.3 kpc) | 173 ± 17 km/s | 199,0 km/s | +1.53σ |
| ρdark(R⊙ = 8 kpc) | 0,39 ± 0,03 GeV/cm³ | 0,47 GeV/cm³ | +2.3σ |
| Mdark(<8 kpc) | ~5 × 10¹⁰ M⊙ | 5.3 × 10¹⁰ M⊙ | blisko |
| Mtot(<200 kpc) | 5-9 × 10¹¹ M⊙ | 3.3 × 10¹¹ M⊙ | low end |
Uproszczone dopasowanie daje χ²/dof ≈ 0,31. Najtrudniejszym punktem pozostaje najbardziej zewnętrzna wartość z epoki Gaia przy 27,3 kpc, gdzie obserwowany spadek jest ostrzejszy niż przewiduje ten dwuparametrowy model.
1. Obcinanie dysku – dlaczego i w jaki sposób
1.1 Problem z nieskończonym dyskiem wykładniczym
Wzór na dysk Freemana zakłada wykładniczą gęstość powierzchniową rozciągającą się do nieskończoności. Matematycznie nigdy nie osiąga ona zera, ale fizycznie dysk gwiezdny Drogi Mlecznej ma skończony zasięg. Poza efektywną krawędzią gwiezdną, zamknięta masa barionowa jest zasadniczo stała, a wkład prędkości musi spadać w przybliżeniu jako keplerowskie pole punkt-masa.
\(\Sigma(R)=\Sigma_0e^{-R/R_d}\)Za krawędzią dysku prędkość barionowa dąży w kierunku:
\(V_{\mathrm{bar}}(R)\xrightarrow{R\gg R_d}\sqrt{\frac{GM_{\mathrm{bar,tot}}}{R}}\) \(M_{\mathrm{bar,tot}}=M_{\mathrm{disk}}+M_{\mathrm{bulge}}\approx4.7\times10^{10}M_\odot\)Przykładowe wartości to:
\(V_{\mathrm{bar}}(30\,\mathrm{kpc})\approx82\,\mathrm{km/s},\qquad V_{\mathrm{bar}}(50\,\mathrm{kpc})\approx63\,\mathrm{km/s}\)1.2 Formuła gładkiego obcinania
Symulacja wykorzystuje płynne przejście między formułą dysku Freemana a wartością keplerowską. Przejście jest wyśrodkowane na Rtrunc = 4Rd = 10.4 kpc z szerokością σ = 1.5 kpc.
\(V_{\mathrm{bar}}(R)=\sqrt{(1-w)V_{\mathrm{Freeman}}^2(R)+w\,\min(V_{\mathrm{Freeman}},V_{\mathrm{Kepler}})^2}\) \(w(R)=\frac{1}{2}\left[1+\tanh\left(\frac{R-R_{\mathrm{trunc}}}{\sigma}\right)\right]\) \(R_{\mathrm{trunc}}=4R_d=10.4\,\mathrm{kpc},\qquad \sigma=1.5\,\mathrm{kpc}\)Funkcja minimum zapobiega przekroczeniu przez dysk barionowy fizycznej granicy keplerowskiej poza krawędzią dysku.
| R | VFreeman | VKeplerian | Vbar, obcięty | System dominujący |
|---|---|---|---|---|
| 5 kpc | 174,5 km/s | 201,1 km/s | 174,5 km/s | Freeman |
| 8 kpc | 161,5 km/s | 159,0 km/s | 161,5 km/s | Freeman ≈ Kepler |
| 10.4 kpc | 143,0 km/s | 139,3 km/s | 141,2 km/s | Przejście |
| 16 kpc | 112,4 km/s | 112,4 km/s | 112,4 km/s | Keplerian |
| 25 kpc | 89,9 km/s | 89,9 km/s | 89,9 km/s | Keplerian |
| 50 kpc | 63,6 km/s | 63,6 km/s | 63,6 km/s | Keplerian |
2. Teoria Bee 3D gęstości ciemnej masy
2.1 Pierścienie dyskowe promieniujące w 3D
Każdy pierścień dysku galaktycznego o promieniu R′ i szerokości dR′ ma masę:
\(dM=\Sigma(R’)\,2\pi R’\,dR’\)W teorii BeeTheory pierścień ten generuje pole fal grawitacyjnych, które rozchodzi się we wszystkich trzech wymiarach przestrzennych. W przybliżeniu monopolowym odległość do punktu pola 3D o promieniu sferycznym r wynosi:
\(D(r,R’)=\sqrt{r^2+R’^2}\)Postać liczbowa gęstości ciemnej jest następująca:
\(\rho_{\mathrm{dark}}(r)=K\sum_{i=1}^{N}\Sigma_0e^{-R’_i/R_d}\frac{(1+\alpha D_i)e^{-\alpha D_i}}{D_i^2}\,2\pi R’_i\Delta R’\) \(D_i=\sqrt{r^2+R_i’^2},\qquad R’_i=\left(i-\frac{1}{2}\right)\frac{R_{\mathrm{max}}}{N}\) \(N=60,\qquad R_{\mathrm{max}}=25\,\mathrm{kpc}\)2.2 Zamknięta ciemna masa i prędkość kołowa
\(M_{\mathrm{dark}}(<r)=\int_0^r4\pi s^2\rho_{\mathrm{dark}}(s)\,ds\) \(M_{\mathrm{dark}}(<r)\approx\sum_{j=1}^{30}4\pi r_j^2\rho_{\mathrm{dark}}(r_j)\Delta r\) \(V_{\mathrm{dark}}(R)=\sqrt{\frac{GM_{\mathrm{dark}}(<R)}{R}}\) \(V_c(R)=\sqrt{V_{\mathrm{bar}}^2(R)+V_{\mathrm{dark}}^2(R)}\)2.3 Zachowanie asymptotyczne
\(\rho_{\mathrm{dark}}(r)\approx\frac{2\pi K\Sigma_0R_d^2}{r^2}\left(1+\alpha r+\frac{\alpha^2r^2}{2}\right)e^{-\alpha r}\)Dla αr ≪ 1:
\(\rho_{\mathrm{dark}}(r)\xrightarrow{\alpha r\ll1}\frac{2\pi K\Sigma_0R_d^2}{r^2}\) \(M_{\mathrm{dark}}(<r)\propto r\qquad\Longrightarrow\qquad V_{\mathrm{dark}}\approx\mathrm{stała}\)3. Wyniki symulacji – wykresy interaktywne
Poniższa symulacja zachowuje model numeryczny, suwaki, krzywą rotacji, profil masy, profil gęstości i aktualizację χ² na żywo. Proszę wkleić tę stronę do WordPressa z włączonym wykonywaniem skryptów.
χ²/dof: – | ℓ: – kpc | ρ(R⊙): – GeV/cm³
| r (kpc) | Mbar (10¹⁰ M⊙) | Mdark (10¹⁰ M⊙) | Mtot (10¹⁰ M⊙) | DM/bar | ρdark (GeV/cm³) |
|---|---|---|---|---|---|
| Ładowanie… | |||||
4. Interpretacja fizyczna i uniwersalność
4.1 Długość koherencji
Wewnątrz długości koherencji jądro Yukawy zachowuje się prawie jak newtonowskie jądro 1/D². Gęstość ciemności podąża w przybliżeniu za r-², a krzywa rotacji jest płaska. Powyżej ℓ, wykładnicze tłumienie powoduje spadek obserwowany w zewnętrznym dysku.
\(\ell=\frac{1}{\alpha}\approx11.5\,\mathrm{kpc}\) \(\frac{\ell}{R_d}=\frac{11.5}{2.6}\approx4.4\)4.2 Sprzężenie bezwymiarowe
Bezwymiarowe sprzężenie BeeTheory można zdefiniować jako:
\(\lambda_{\mathrm{galaxy}}=K\ell^2\) \(\lambda_{\mathrm{galaxy}}=0.040\times(11.5)^2\approx5.3\)Jest to rząd wielkości porównywalny ze sprzężeniem wynikającym z kalibracji H₂, gdzie λ wynosi około 3-4. Ewentualna uniwersalność tej liczby w skali pozostaje głównym otwartym pytaniem.
4.3 Porównanie z modelami standardowymi
| Model | Parametry | Typowe dopasowanie | Skala | Mechanizm |
|---|---|---|---|---|
| Aureola izotermiczna | 2 | Umiarkowany | promień rdzenia | Płaska krzywa fenomenologiczna |
| Profil NFW | 2 | Silny | rs | Profil symulacji N-ciała |
| Einasto | 2-3 | Silny | r-2 | Elastyczny profil empiryczny |
| BeeTheory 3D Yukawa | 2 | Obiecujący | ℓ | Sprzężenie fala-masa z dysku |
Najbardziej zewnętrzny punkt z ery Gaia pozostaje najtrudniejszym ograniczeniem. Ostrzejszy spadek można uzyskać przy mniejszej długości koherencji, ale pogarsza to dopasowanie wewnętrzne. Przyszłe dane z Gaia DR4, gromad kulistych i strumieni gwiazd będą ważnymi testami.
Referencje
- Ou, X. et al. - The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 2024.
- Dutertre, X. - Bee Theory™: Wave-Based Modeling of Gravity, BeeTheory.com v2, 2023.
- Freeman, K. C. - On the disks of spiral and S0 galaxies, ApJ 160, 811, 1970.
- McMillan, P. J. - Rozkład masy i potencjał grawitacyjny Drogi Mlecznej, MNRAS 465, 76, 2017.
- Navarro, J. F., Frenk, C. S., White, S. D. M. - A Universal Density Profile from Hierarchical Clustering, ApJ 490, 1997.
BeeTheory.com - Kwantowa grawitacja oparta na falach
© Technoplane S.A.S. - 2025 r.