BeeTheory – Foundations – Uwaga techniczna XIV
Krok 1 – Droga Mleczna:
Zastosowanie jądra Yukawy BeeTheory
Metodologia przedstawiona w Nocie XII została zastosowana do Drogi Mlecznej przy użyciu jawnego jądra fali Yukawy $\mathcal{K}(D) = K_0\,(1+\alpha D)\,e^{-\alpha D}/D^2$, z całkowaniem przeprowadzonym oddzielnie dla każdego składnika barionowego zgodnie z jego geometrią. Wynik jest porównywany punkt po punkcie z krzywą rotacji Gaia 2024 i „brakującą masą” modelu standardowego. Niniejsza notatka ustanawia podstawy ram zastosowanych z nowym, w pełni geometrycznym sformułowaniem.
1. Wynik pierwszy
Wynik bazowy – Droga Mleczna z jawnym jądrem Yukawy
Krzywa rotacji. Model odtwarza prędkość Gaia 2024 z dokładnością do 2 km/s przy $R = 4$ kpc, ale zawyża przewidywania o $+33$ km/s przy promieniu Słońca ($R = 8$ kpc) i o $+64$ km/s przy $R = 27.3$ kpc. $\chi^2/\text{dof} = 1.27$.
Brakująca masa. Masa pola falowegoBeeTheory odpowiada „brakującej masie” modelu standardowego z dokładnością do 5% przy $R = 4$ kpc, ale przekracza ją o czynnik 2,2 przy $R = 27,3$ kpc. Model wytwarza zbyt dużą masę ciemnego pola przy dużych promieniach.
Lokalna gęstość. $\rho_\text{wave}(R_\odot) = 0.72$ GeV/cm³, w porównaniu do obserwowanego zakresu $0.39$-0.45$ GeV/cm³. Przewidywania są zawyżone o około 1,7 raza.
Co to oznacza
Wyraźne sformułowanie Yukawy tworzy krzywą rotacji, która jest zbyt płaska przy dużych promieniach. Długość rozpadu pola falowego $\ell$ jest zbyt długa, pozwalając polu falowemu na utrzymywanie masy poza widocznym dyskiem. Jest to strukturalny punkt odniesienia przed wprowadzeniem udoskonalenia gęstości powierzchniowej, o którym mowa w Uwadze XI.
2. Co zamierzamy obliczyć
Droga Mleczna jest naturalnym przypadkiem testowym, ponieważ jest galaktyką, na której pierwotnie skalibrowano globalne sprzężenie $lambda$, a także dlatego, że istnieją dwie niezależne obserwacje, z którymi można je porównać:
(a) Krzywa rotacji $V_c(R)$ z Gaia 2024 (Ou et al., MNRAS 528), która mierzy prędkość kołową na dziesięciu promieniach od $R = 2$ kpc do $R = 27.3$ kpc z niepewnością statystyczną $7$-17$ km/s. Jest to prędkość, którą BeeTheory musi odtworzyć, łącząc wkład barionowy $V_\text{bar}(R)$ z wkładem pola falowego $V_\text{wave}(R)$. Jest to prędkość, którą BeeTheory musi odtworzyć, łącząc wkład barionowy $V_\text{bar}(R)$ z wkładem pola falowego $V_\text{wave}(R)$.
(b) „Brakująca masa” $M_text{missing}(. Standardowa interpretacja odwołuje się do cząstek ciemnej materii, aby zapewnić tę masę. BeeTheory przewiduje natomiast, że masa ta jest polem falowym $M_\text{wave}(
(c) Lokalna gęstość ciemnej materii na Słońcu, zmierzona przy $\rho \ około 0,39$-0,45$ GeV/cm³ na podstawie badań kinematycznych okolic Słońca. BeeTheory przewiduje wartość $\rho_\text{wave}(R_\odot)$ na podstawie tych samych obliczeń pola falowego.
Zgodność (lub niezgodność) tych trzech obserwabli testuje trzy różne aspekty modelu: jego kształt krzywej rotacji, profil masy zamkniętej i lokalną normalizację gęstości.
3. Jądro fali – postać jawna
Każdy element masy barionowej generuje pole falowe BeeTheory, którego natężenie w punkcie pola oddzielonym odległością $D$ jest określone przez:
Jądro fali kształtu Yukawy
$$\mathcal{K}(D) \;=\; K_0 \cdot \frac{(1 + \alpha D)\,e^{-\alpha D}}{D^2}, \qquad \alpha = \frac{1}{\ell}$$
Wykładnicze tłumienie $e^{-\alpha D}$ zapewnia, że pole falowe ma skończoną masę całkowitą – bez niego całka pola falowego rozchodziłaby się w nieskończoności. Współczynnik wstępny $(1 + alfa D)$ pochodzi z uregulowanej funkcji falowej z uwagi I; wraz z wykładnikiem sprawia on, że struktura przestrzenna jądra jest quasi-newtonowska przy $D ll ell$ i wykładniczo tłumiona przy $D gg ell$.
Długość charakterystyczna $\ell$ zależy od komponentu generującego pole:
| Komponent | Długość koherencji $\ell$ (kpc) | Skala geometryczna |
|---|---|---|
| Wybrzuszenie (3D Hernquist) | $\ell_b = c_\text{sph}\,r_b = 0,41 \times 0,61$ | $0.25$ |
| Cienki dysk | $\ell_\text{thin} = c_\text{disk}\,R_d = 3,17 \times 2,6$ | $8.24$ |
| Gruby dysk | $\ell_\text{thick} = c_\text{disk}\,(1.5\,R_d)$ | $12.36$ |
| Pierścień gazowy | $\ell_\text{gas} = c_\text{disk}\,(1.7\,R_d)$ | $14.01$ |
| Ramiona spiralne | $\ell_\text{arm} = c_\text{arm}\,R_d = 2,0 \times 2,6$ | $5.20$ |
4. Geometria integracji, komponent po komponencie
Dla każdego komponentu całkujemy gęstość źródła sprzężoną z jądrem, używając odpowiedniego elementu objętości dla geometrii. Wynikiem jest gęstość pola falowego $\rho_\text{wave}^{(i)}(r)$ w punkcie pola.
4.1 Integracja wybrzuszenia z powłoką sferyczną
Wybrzuszenie jest trójwymiarowym rozkładem sferycznym. Każda cienka powłoka o promieniu $r’$ zawiera masę $\rho_b(r’)\,4\pi r’^2\,dr’$ i wnosi wkład do pola w odległości promieniowej $r$ od środka. W przybliżeniu monopolarnym efektywna separacja między punktem pola a ogólnym punktem na powłoce wynosi $D = \sqrt{r^2 + r’^2}$:
$$\rho_\text{wave}^{(b)}(r) \;=\; \int_0^{6r_b} \rho_b(r’) \cdot K_0\,\frac{(1 + \alpha_b D)\,e^{-\alpha_b D}}{D^2} \cdot 4\pi r’^2 \, dr’, \quad D = \sqrt{r^2 + r’^2} $$
Z $\rho_b(r’) = M_b r_b / [2\pi r'(r’+r_b)^3]$ (Hernquist) i $\alpha_b = 1/0.25 = 4.0$ kpc$^{-1}$. Całkowanie jest przerywane w punkcie $6\,r_b$, poza którym gęstość jest numerycznie zaniedbywalna.
4.2 Cienkie i grube dyski – integracja pierścienia koncentrycznego
Każdy dysk jest cienkim rozkładem osiowo-symetrycznym. Dysk jest podzielony na koncentryczne pierścienie o promieniu galaktocentrycznym $R’$ i szerokości $dR’$, z których każdy ma masę $\Sigma(R’)\,2\pi R’\,dR’$. Wkład do punktu pola na promieniu $r$ (w płaszczyźnie dysku) wymaga efektywnej separacji $D = \sqrt{r^2 + R’^2}$ w tym samym przybliżeniu monopolarnym:
$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{8R_d} \Sigma_\text{thin}(R’) \cdot K_0\,\frac{(1 + \alpha_\text{thin} D)\,e^{-\alpha_\text{thin} D}}{D^2} \cdot 2\pi R’ \, dR’, \quad D = \sqrt{r^2 + R’^2}$$
z $\Sigma_\text{thin}(R’) = (M_\text{thin}/2\pi R_d^2)\,e^{-R’/R_d}$ i $\alpha_\text{thin} = 1/8,24 = 0,121$ kpc$^{-1}$. Gruby dysk ma identyczną skalę: $\Sigma_\text{thick}$, $R_\text{thick} = 3.9$ kpc, $\alpha_\text{thick} = 0.081$ kpc$^{-1}$.
4.3 Pierścień gazowy – integracja pierścienia z centralnym zubożeniem
Rozkład gazu ma centralną dziurę (zaniedbywalne HI przy $R \lesssim 2$ kpc) i rozciąga się dalej niż dysk gwiezdny. Profil $\Sigma_\text{gas}(R’) = \Sigma_0\,\exp(-R_\text{hole}/R’ – R’/R_g)$ oddaje obie cechy:
$$\rho_\text{wave}^{(\text{gas})}(r) \;=\; \int_{R_\text{hole}}^{8R_g} \Sigma_\text{gas}(R’) \cdot K_0\,\frac{(1 + \alpha_\text{gas} D)\,e^{-\alpha_\text{gas} D}}{D^2} \cdot 2\pi R’ \, dR’$$
z $R_g = 4.42$ kpc, $R_\text{hole} = 2.21$ kpc, $\alpha_\text{gas} = 0.071$ kpc$^{-1}$. Większa długość koherencji odzwierciedla bardziej rozległy rozkład gazu.
4.4 Ramiona spiralne – integracja pierścieniowa ze zmniejszoną amplitudą i węższym jądrem
Ramiona spiralne niosą $10\%$ gęstości powierzchniowej cienkiego dysku i mają własną długość koherencji $\ell_\text{arm} = 5,2$ kpc, węższą niż w dysku, aby odzwierciedlić azymutalną koncentrację ramion:
$$\rho_\text{wave}^{(\text{arm})}(r) \;=\; \int_0^{8R_d} 0.10\,\Sigma_\text{thin}(R’) \cdot K_0\,\frac{(1 + \alpha_\text{arm} D)\,e^{-\alpha_\text{arm} D}}{D^2} \cdot 2\pi R’ \, dR’$$
4.5 Całkowita gęstość pola falowego
$$\rho_\text{wave}(r) \;=\; \lambda \,\sum_{i} \rho_\text{wave}^{(i)}(r), \quad \lambda = 0,189$$.
5. Od gęstości fali do krzywej rotacji
Gdy $\rho_\text{wave}(r)$ jest znane dla każdego promienia, całkowita masa zamkniętego pola falowego jest uzyskiwana przez całkowanie promieniowe:
$$M_\text{wave}(R) \;=\; \int_0^{R} 4\pi r^2 \rho_\text{wave}(r) \, dr$$
Przewidywana prędkość kołowa łączy następnie wkłady pola barionowego i falowego w kwadraturze, zgodnie z relacją Newtona:
$$V_c^2(R) \;=\; V_\text{bar}^2(R) \;+\; \frac{G\,M_\text{wave}(R)}{R}$$
6. Krzywa rotacji – wyniki punkt po punkcie
| $R$ (kpc) | $V_\text{bar}$ (km/s) | $M_\text{wave}/10^{10}$ | $V_\text{wave}$ (km/s) | $V_\text{tot}$ (km/s) | $V_\text{obs}$ (Gaia 2024) | $\Delta$ |
|---|---|---|---|---|---|---|
| 2.0 | 157.8 | 0.67 | 120.0 | 198.3 | 250 ± 12 | -51.7 |
| 4.0 | 164.1 | 2.54 | 165.4 | 233.0 | 235 ± 10 | -2.0 |
| 6.0 | 163.5 | 5.17 | 192.5 | 252.6 | 230 ± 8 | +22.6 |
| 8.0 | 157.3 | 8.13 | 209.1 | 261.7 | 229 ± 7 | +32.7 |
| 10.0 | 148.7 | 11.18 | 219.3 | 265.0 | 224 ± 8 | +41.0 |
| 12.0 | 139.6 | 14.15 | 225.2 | 265.0 | 217 ± 9 | +48.0 |
| 15.0 | 126.7 | 18.29 | 229.0 | 261.8 | 208 ± 10 | +53.8 |
| 20.0 | 109.3 | 24.10 | 227.7 | 252.6 | 195 ± 12 | +57.6 |
| 25.0 | 96.7 | 28.54 | 221.6 | 241.8 | 180 ± 15 | +61.8 |
| 27.3 | 92.0 | 30.18 | 218.1 | 236.7 | 173 ± 17 | +63.7 |
Model doskonale pasuje do obserwacji przy $R = 4$ kpc (Δ = -2 km/s), ale wraz ze wzrostem promienia coraz bardziej zawyża przewidywania. Przy promieniu Słońca nadprzewidywania wynoszą +33 km/s (4,7σ powyżej niepewności Gaia). Na zewnętrznej granicy $R = 27,3$ kpc, zawyżone przewidywania sięgają +64 km/s (3,8σ). Przewidywana krzywa jest zbyt płaska – pole falowe nadal wnosi masę poza widoczny dysk, ponieważ pozwala na to wykładnicze odcięcie przy $D \sim \ell$.
7. Masa falowa a „brakująca masa” modelu standardowego
Dla każdego promienia porównujemy trzy wielkości: zamkniętą masę barionową (tylko materia widzialna), masę dynamiczną wymaganą przez obserwowaną prędkość (prawo Newtona zastosowane do $V_\text{obs}$) oraz masę pola falowego BeeTheory. Różnica między drugą a pierwszą masą jest tym, co model standardowy nazywa „brakującą masą”:
$$M_\text{missing}(<R) \;=\; \frac{R\,V_\text{obs}^2(R)}{G} \;-\; M_\text{bar}(<R)$$.
Stosunek $M_text{wave}/M_text{missing}$ mówi nam, jak dobrze pole falowe BeeTheory zastępuje cząsteczkową ciemną materię na zasadzie promień po promieniu:
| $R$ (kpc) | $M_\text{bar}(| $M_\text{dyn}( | $M_\text{missing}$. | $M_\text{wave}$ (BT) | Współczynnik $M_\text{wave}$/M_\text{miss}$ | |
|---|---|---|---|---|---|
| 2.0 | 1.16e+10 | 2.91e+10 | 1.75e+10 | 6.69e+09 | 0.38 |
| 4.0 | 2.51e+10 | 5.13e+10 | 2.63e+10 | 2.54e+10 | 0.97 |
| 6.0 | 3.73e+10 | 7.38e+10 | 3.65e+10 | 5.17e+10 | 1.42 |
| 8.0 | 4.60e+10 | 9.75e+10 | 5.15e+10 | 8.13e+10 | 1.58 |
| 10.0 | 5.14e+10 | 1.17e+11 | 6.52e+10 | 1.12e+11 | 1.71 |
| 12.0 | 5.44e+10 | 1.31e+11 | 7.70e+10 | 1.41e+11 | 1.84 |
| 15.0 | 5.60e+10 | 1.51e+11 | 9.49e+10 | 1.83e+11 | 1.93 |
| 20.0 | 5.56e+10 | 1.77e+11 | 1.21e+11 | 2.41e+11 | 1.99 |
| 25.0 | 5.43e+10 | 1.88e+11 | 1.34e+11 | 2.85e+11 | 2.13 |
| 27.3 | 5.37e+10 | 1.90e+11 | 1.36e+11 | 3.02e+11 | 2.22 |
Odczyt ilościowy
Przy $R = 4$ kpc pole falowe zasadniczo odpowiada brakującej masie (stosunek 0,97$). Pomiędzy $R = 6$ i $R = 8$ kpc model już przekracza brakującą masę o 40-60%. Powyżej $R = 15$ kpc, masa pola falowego jest w przybliżeniu dwukrotnie większa niż ta, którą model standardowy określa jako ciemną materię. Model wytwarza dodatkową masę przy dużych promieniach – jest to dokładnie objaw długości koherencji $\ell$, która jest zbyt długa dla widocznego dysku gwiezdnego.
8. Składowe pola falowego w promieniu Słońca
Oszacowanie wkładu każdego składnika do $\rho_\text{wave}(R_\odot = 8\,\text{kpc})$ pokazuje, które źródło barionowe dominuje w tym polu falowym:
| Komponent | $\rho_\text{wave}(R_\odot)$ ($M_\odot$/kpc³) | Ułamek całości |
|---|---|---|
| Cienki dysk gwiezdny | $6,05 razy 10^7$. | 60.6% |
| Gruby dysk gwiezdny | $1,91 razy 10^7$. | 19.1% |
| Pierścień gazowy | $1,62 razy 10^7$. | 16.2% |
| Ramiona spiralne | $4,15 razy 10^6$. | 4.1% |
| Wybrzuszenie | $1,55 razy 10^{-5}$. | $\sim$0% |
Cienki dysk gwiezdny dominuje w polu falowym w pozycji Słońca (60%), a gruby dysk i pierścień gazowy mają mniej więcej równy wkład (16-19%). Wybrzuszenie ma pomijalny wkład, ponieważ $\ell_b = 0,25$ kpc jest znacznie mniejsze niż $R_\odot = 8$ kpc – wykładnicze tłumienie w jądrze zabija wkład wybrzuszenia w tej odległości.
Przeliczenie całkowitej gęstości na jednostki fizyki cząstek daje $\rho_\text{wave}(R_\odot) = 0,717$ GeV/cm³, co można porównać z pomiarem kinematycznym 0,39$-0,45$ GeV/cm³ z Read 2014 i późniejszymi analizami. Przewidywania zawyżają obserwowaną lokalną gęstość o czynnik 1,6-1,8$ – zgodnie z przewidywaniami krzywej rotacji dla tego samego promienia.
9. Co określa ten punkt odniesienia
Mechanizm działa w zasadzie
W okolicy $R = 4$ kpc – centralnej części dysku – zintegrowane pole falowe jest równe brakującej masie modelu standardowego z dokładnością do 5%, a krzywa rotacji jest odtwarzana z dokładnością do 2 km/s. Jądro falowe, zastosowane do widocznych barionów, wytwarza masę grawitacyjną ilościowo porównywalną z cząsteczkową ciemną materią w tym promieniu. Żadna nowa cząstka nie jest potrzebna; pole falowe widzialnej materii odpowiada za brakującą grawitację.
Krzywa jest jednak zbyt płaska przy dużych promieniach
Poza centralnym dyskiem model zawyża prędkość rotacji o wartość, która rośnie monotonicznie wraz z promieniem. Pole falowe nadal gromadzi masę poza widocznym dyskiem, ponieważ długość koherencji jądra $\ell_\text{thin} = 8,24$ kpc jest porównywalna z rozmiarem samego dysku, co pozwala na znaczący wkład przy $D = 15$-25$ kpc. Z kolei krzywa rotacji Gaia zmniejsza się nieznacznie powyżej $R \sim 10$ kpc – cecha, której obecne sformułowanie nie odtwarza.
Linia bazowa, a nie ostateczna odpowiedź
Obliczenia te stanowią punkt odniesienia dla modelu, w którym $\ell_i$ zależy liniowo od samego $R_d$. Diagnoza z Uwagi XI wykazała, że $\Sigma_d$ – centralna gęstość powierzchniowa – musi zostać uwzględniona przy wyznaczaniu $\ell_i$, aby skorygować krzywą przy dużych promieniach. Im gęstszy dysk, tym bardziej lokalna powinna być odpowiedź falowa. Uwzględnienie tego udoskonalenia jest przedmiotem kolejnych uwag. Podana tutaj linia bazowa Drogi Mlecznej jest tym, co te notatki muszą poprawić.
10. Podsumowanie
1. Krzywa rotacji Drogi Mlecznej jest obliczana poprzez całkowanie każdego składnika barionowego względem jądra fali Yukawy $mathcal{K}(D) = K_0(1 + alfa D),e^{-alfa D}/D^2$, z odpowiednią geometrią: sferyczne powłoki dla wybrzuszenia, koncentryczne pierścienie dla dysków, gazu i ramion spiralnych.
2. Przy $R = 4$ kpc, masa pola falowego BeeTheory zgadza się z „brakującą masą” modelu standardowego z dokładnością do 5% (stosunek 0,97), a przewidywana prędkość odpowiada Gaia 2024 z dokładnością do 2 km/s.
3. Przy promieniu Słońca ($R = 8$ kpc), model przeszacowuje prędkość rotacji o $+33$ km/s, a lokalną gęstość ciemnej materii o czynnik 1,6 – obie wartości są ze sobą zgodne.
4. Poza $R = 15$ kpc, przewidywana masa pola falowego przekracza brakującą masę modelu standardowego o czynnik 2 lub więcej. Przewidywana krzywa rotacji nie zmniejsza się tak, jak wymagają tego dane Gaia.
5. Cienki dysk gwiezdny dominuje w polu falowym w pozycji Słońca (60% $\rho_\text{wave}$). Wybrzuszenie ma pomijalny wkład. Rozkład jest zgodny z opisaną geometrią integracji.
6. Nadmierne przewidywanie przy dużym $R$ jest strukturalną sygnaturą $\ell_i$, która jest zbyt długa. Proszę zauważyć, że $\Sigma_d$ musi wejść do wzoru na długość koherencji. Kolejnym krokiem jest udoskonalenie $\ell_i$ poprzez $\Sigma_d$.
Referencje. Ou, X. et al. – The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 693 (2024). Krzywa rotacji Gaia 2024. – Bland-Hawthorn, J., Gerhard, O. – The Galaxy in Context, ARA&A 54, 529 (2016). Rozkład strukturalny MW. – Read, J. I. – The Local Dark Matter Density, J. Phys. G 41, 063101 (2014). Lokalne pomiary gęstości ciemnej materii. – Freeman, K. C. – On the disks of spiral and S0 galaxies, ApJ 160, 811 (1970). – Hernquist, L. – An analytical model for spherical galaxies and bulges, ApJ 356, 359 (1990). – Dutertre, X. – Bee Theory™: Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).
BeeTheory.com – Kwantowa grawitacja oparta na falach – Aplikacja Krok 1 – © Technoplane S.A.S. 2026