BeeTheory – Symulacja dwóch reżimów – 2025 r.
BeeTheory Galaktyczna ciemna masa: Wybrzuszenie + Dysk, dwa reżimy, cztery parametry
Krzywa rotacji Gaia 2024 ma dwa różne reżimy: zdominowany przez wybrzuszenie poniżej 5,5 kpc i zdominowany przez dysk powyżej. BeeTheory wychwytuje oba z oddzielną długością koherencji dla każdego komponentu, dając χ²/dof = 0,24.
BeeTheory.com – Ou et al. MNRAS 528, 2024 – McMillan MNRAS 465, 2017
χ²/dof = 0,24
Najlepsze dopasowanie do tej pory
ℓbulge = 0,6 kpc
Kompaktowe źródło krótkiego zasięgu
ℓdisk = 11,1 kpc
Rozszerzone źródło dalekiego zasięgu
ρ(R⊙) = 0,37 GeV/cm³
vs obserwowane 0,39 GeV/cm³
0. Wyniki – najpierw parametry i równania
Całkowita gęstość ciemnej masy w sferycznym promieniu r od Centrum Galaktyki jest sumą dwóch niezależnych pól BeeTheory: jednego ze zwartego wybrzuszenia 3D i jednego z rozszerzonego dysku 2D. Każdy składnik ma swoją własną długość koherencji.
\(\rho_{\mathrm{dark}}(r)=\rho_{\mathrm{dark,b}}(r)+\rho_{\mathrm{dark,d}}(r)\) \(\rho_{\mathrm{dark,b}}(r)=K_b\int_0^{R_b}\rho_{0,b}e^{-r’/r_b}\frac{(1+\alpha_bD)e^{-\alpha_bD}}{D^2}\,4\pi r’^2\,dr’\) \(\rho_{\mathrm{dark,d}}(r)=K_d\int_0^{R_d^{\max}}\Sigma_0e^{-R’/R_d}\frac{(1+\alpha_dD)e^{-\alpha_dD}}{D^2}\,2\pi R’\,dR’\) \(D=\sqrt{r^2+r’^2}\quad\mathrm{lub}\quad D=\sqrt{r^2+R’^2}\quad\mathrm{(aproksymacja monopolu)}\)Cztery dopasowane parametry są niezależne: długość koherencji wybrzuszenia rządzi wewnętrzną krzywą rotacji, a długość koherencji dysku rządzi krzywą zewnętrzną.
\(K_b=1.055\,\mathrm{kpc}^{-1},\qquad \alpha_b=1.634\,\mathrm{kpc}^{-1},\qquad \ell_b=\frac{1}{\alpha_b}=0.61\,\mathrm{kpc}\) \(K_d=0.02365\,\mathrm{kpc}^{-1},\qquad \alpha_d=0.0902\,\mathrm{kpc}^{-1},\qquad \ell_d=\frac{1}{\alpha_d}=11.1\,\mathrm{kpc}\)Bulge – Reżim 1
R < 5,5 kpc
Kompaktowe źródło sferyczne. Krótka koherencja oznacza, że pole falowe jest intensywne w pobliżu środka i stromo opada. Kontroluje rosnącą część krzywej rotacji, od około 220 do 232 km/s.
\(K_b=1.055\,\mathrm{kpc}^{-1}\) \(\alpha_b=1.634\,\mathrm{kpc}^{-1}\) \(\ell_b=0.61\,\mathrm{kpc}\) \(r_b=1.5\,\mathrm{kpc}\) \(M_b=1.24\times10^{10}M_\odot\) \(\lambda_b=K_b\ell_b^2=0.39\)Dysk – reżim 2
R > 5,5 kpc
Rozszerzony dysk wykładniczy. Długa spójność pozwala polu falowemu wypełnić halo w skali galaktycznej, podtrzymując płaską krzywą rotacji, a następnie powodując spadek Gaia 2024.
\(K_d=0.02365\,\mathrm{kpc}^{-1}\) \(\alpha_d=0.0902\,\mathrm{kpc}^{-1}\) \(\ell_d=11.1\,\mathrm{kpc}\) \(R_d=3.5\,\mathrm{kpc}\) \(M_d=5.47\times10^{10}M_\odot\) \(\lambda_d=K_d\ell_d^2=2.90\)Podsumowanie dopasowania
| Obserwowalne | Gaia 2024 | BeeTheory | Proszę ciągnąć |
|---|---|---|---|
| Vc(4 kpc), reżim wewnętrzny | 220 ± 10 km/s | 220,9 km/s | +0.09σ |
| Vc(6 kpc), przegięcie | 232 ± 7 km/s | 229,6 km/s | -0.35σ |
| Vc(8 kpc), okrąg słoneczny | 230 ± 6 km/s | 231,2 km/s | +0.20σ |
| Vc(16 kpc), zewnętrzny płaskowyż | 222 ± 8 km/s | 218,9 km/s | -0.38σ |
| Vc(27.3 kpc), najbardziej zewnętrzny | 173 ± 17 km/s | 195,3 km/s | +1.31σ |
| ρdark(R⊙) | 0,39 ± 0,03 GeV/cm³ | 0,372 GeV/cm³ | -0.6σ |
| Mdark(<8 kpc) | ~5 × 10¹⁰ M⊙ | 4.83 × 10¹⁰ M⊙ | blisko |
| Mtot(<200 kpc) | 5-9 × 10¹¹ M⊙ | 3.1 × 10¹¹ M⊙ | low end |
1. Odczytywanie krzywej rotacji Gai – dwa reżimy fizyczne
Krzywa rotacji Gaia DR3 ma wyraźny punkt przegięcia w pobliżu R ≈ 5.5 kpc.
- Reżim 1, R = 4-5,5 kpc: Vc wzrasta z około 220 do 232 km/s. Gradient prędkości dV/dR > 0 wskazuje na zwartą masę centralną , której ciemne pole szybko rośnie wraz z promieniem. Jest to sygnatura wybrzuszenia.
- Reżim 2, R = 5,5-27 kpc: Vc jest płaskie w pobliżu 230 km/s, a następnie powoli spada. Gradient jest początkowo prawie płaski i staje się bardziej ujemny w kierunku najbardziej oddalonego punktu Gaia. Jest to sygnatura dysku-halo.
Fizyczny powód dwóch różnych długości koherencji
Wybrzuszenie jest zwarte i skoncentrowane. Długość koherencji pola falowego jest porównywalna z fizyczną skalą samego źródła.
\(\ell_b=0.61\,\mathrm{kpc}\approx0.4r_b\)Dysk jest rozszerzony. Jego pole falowe ma znacznie większą długość koherencji, co pozwala mu utrzymać zewnętrzną krzywą rotacji w odległościach galaktycznych.
\(\ell_d=11.1\,\mathrm{kpc}\approx3.2R_d\)2. Uproszczony model barionowy – dwa składniki
Wszystkie bariony galaktyczne są asymilowane do dwóch rodzin geometrycznych: zwartego wybrzuszenia sferycznego i rozszerzonego dysku wykładniczego.
Składnik wybrzuszenia – sferyczny wykładniczy
\(\rho_b(r)=\rho_{0,b}e^{-r/r_b}\) \(r_b=1.5\,\mathrm{kpc}\) \(M_b=M_{\mathrm{bulge}}+M_{\mathrm{bar,core}}=9.23\times10^9+3.1\times10^9=1.24\times10^{10}M_\odot\)Łączna masa wybrzuszenia wynosi:
\(M_b(<r)=4\pi\rho_{0,b}r_b^3\left[2-\left(2+\frac{2r}{r_b}+\frac{r^2}{r_b^2}\right)e^{-r/r_b}\right]\)Komponent dysku – dysk wykładniczy
\(\Sigma_d(R)=\Sigma_0e^{-R/R_d}\) \(R_d=3.5\,\mathrm{kpc}\) \(\Sigma_0=\frac{M_d}{2\pi R_d^2}\) \(M_d=M_{\mathrm{thin}}+M_{\mathrm{thick,outer}}+M_{\mathrm{HI}}+M_{\mathrm{H_2}}=5.47\times10^{10}M_\odot\)Łączna masa dysku wynosi:
\(M_d(<r)=2\pi\Sigma_0R_d^2\left[1-\left(1+\frac{r}{R_d}\right)e^{-r/R_d}\right]\)Uwzględnienie całej rozszerzonej masy w pojedynczym wykładniku daje efektywny promień skali bliski 3,5 kpc. Jest to ważony masą efektywny promień skali cienkiego dysku, grubego dysku, składników HI i H₂.
Całkowita masa barionowa jest zachowana:
\(M_{\mathrm{bar,total}}=M_b+M_d=1.24\times10^{10}+5.47\times10^{10}=6.71\times10^{10}M_\odot\)3. Równania ciemnej masy BeeTheory na składnik
3.1 Ciemne pole wybrzuszenia
\(\rho_{\mathrm{dark,b}}(r)=K_b\int_0^{R_b}\rho_{0,b}e^{-r’/r_b}\frac{(1+\alpha_bD)e^{-\alpha_bD}}{D^2}\,4\pi r’^2\,dr’\) \(D=\sqrt{r^2+r’^2}\) \(K_b=1.055\,\mathrm{kpc}^{-1},\quad \alpha_b=1.634\,\mathrm{kpc}^{-1},\quad \ell_b=0.61\,\mathrm{kpc},\quad R_b=6\,\mathrm{kpc}\)3.2 Ciemne pole dysku
\(\rho_{\mathrm{dark,d}}(r)=K_d\int_0^{R_d^{\max}}\Sigma_0e^{-R’/R_d}\frac{(1+\alpha_dD)e^{-\alpha_dD}}{D^2}\,2\pi R’\,dR’\) \(D=\sqrt{r^2+R’^2}\) \(K_d=0.02365\,\mathrm{kpc}^{-1},\quad \alpha_d=0.0902\,\mathrm{kpc}^{-1},\quad \ell_d=11.1\,\mathrm{kpc},\quad R_d^{\max}=25\,\mathrm{kpc}\)3.3 Masa całkowita i zamknięta
\(\rho_{\mathrm{dark}}(r)=\rho_{\mathrm{dark,b}}(r)+\rho_{\mathrm{dark,d}}(r)\) \(M_{\mathrm{dark}}(<r)=\int_0^r4\pi s^2\rho_{\mathrm{dark}}(s)\,ds\) \(V_c(R)=\sqrt{\frac{G\left[M_{\mathrm{bar}}(<R)+M_{\mathrm{dark}}(<R)\right]}{R}}\)3.4 Podsumowanie parametrów
| Parametr | Symbol | Wartość | Jednostki | Znaczenie fizyczne |
|---|---|---|---|---|
| Sprzęgło wybrzuszeniowe | Kb | 1.055 | kpc-¹ | Amplituda masy falowej ze zwartego wybrzuszenia. |
| Spójność wybrzuszenia | αb = 1/ℓb | 1.634 | kpc-¹ | Kontroluje wzrost prędkości wewnętrznej. |
| Sprzęgło tarczowe | Kd | 0.02365 | kpc-¹ | Amplituda masy falowej z rozszerzonego dysku. |
| Spójność dysku | αd = 1/ℓd | 0.0902 | kpc-¹ | Kontroluje zewnętrzny płaskowyż i spadek. |
| Skala wypukłości | rb | 1.5 | kpc | Fizyczny promień skali kompaktowego komponentu. |
| Skala dysku | Rd | 3.5 | kpc | Efektywny promień skali dysku ważony masą. |
| Sprzęgło wybrzuszeniowe | λb = Kbℓb² | 0.39 | – | Kompaktowe źródła są mniej wydajne przy dużym promieniu. |
| Sprzęgło tarczowe | λd = Kdℓd² | 2.90 | – | Zgodne z poprzednimi dopasowaniami dysku BeeTheory. |
4. Wyniki symulacji
Poniższa symulacja zachowuje model dwuskładnikowy, niezależne suwaki wybrzuszenia i dysku, krzywą rotacji, profil masy, żywe χ², lokalną gęstość i tabelę masy.
Wybrzuszenie – reżim wewnętrzny
Dysk – reżim zewnętrzny
χ²/dof: – | ℓb: – kpc | ℓd: – kpc | ρ(R⊙): –
| r (kpc) | Mbar | Mdark, wybrzuszenie | Mdark, dysk | Mdark, całkowity | Mtotal | DM/bar | Vc |
|---|---|---|---|---|---|---|---|
| Ładowanie… | |||||||
5. Znaczenie fizyczne - co ujawniają cztery parametry
5.1 Skala długości koherencji wraz z rozmiarem źródła
Najbardziej uderzającym wynikiem dopasowania dwustopniowego jest to, że długość koherencji jest różna dla wybrzuszenia i dysku.
\(\frac{\ell_b}{r_b}=\frac{0.61}{1.5}=0.41\) \(\frac{\ell_d}{R_d}=\frac{11.1}{3.5}=3.17\)Długość koherencji dysku jest około 18 razy większa niż długość koherencji wybrzuszenia. Sugeruje to, że ℓ jest powiązana z geometrią i rozszerzeniem źródła, a nie tylko z całkowitą masą.
Możliwe prawo skalowania do przetestowania na innych galaktykach to:
\(\ell\propto R_{\mathrm{source}}^\gamma\)Zaobserwowany współczynnik wskazuje, że skalowanie może być bardziej strome niż prosta zależność kwadratowo-mnożnikowa lub liniowa.
5.2 Stałe sprzężenia i uniwersalność
\(\lambda_b=K_b\ell_b^2=1.055\times0.37=0.39\) \(\lambda_d=K_d\ell_d^2=0.02365\times123=2.91\)Bezwymiarowe sprzężenie dysku λd ≈ 3 jest zgodne z poprzednimi dopasowaniami BeeTheory. Sprzężenie wybrzuszenia λb ≈ 0,4 jest mniejsze, ponieważ zwarte źródła koncentrują swoją energię falową w pobliżu własnej powierzchni, zamiast rozprzestrzeniać ją na duże odległości galaktyczne.
Podsumowanie: co pokazuje dopasowanie dwóch reżimów
- Krzywa rotacji Gaia zawiera fizyczne informacje o dwóch różnych strukturach masy, a nie tylko o gładkim, jednoskładnikowym halo.
- Przegięcie w pobliżu 5,5 kpc oddziela wewnętrzną galaktykę zdominowaną przez wybrzuszenie od zewnętrznego halo zdominowanego przez dysk.
- BeeTheory rejestruje oba reżimy jednocześnie z czterema parametrami i osiąga χ²/dof = 0,24.
- Długości koherencji mają znaczenie fizyczne: sub-kpc dla zwartego wybrzuszenia i skala galaktyczna dla rozszerzonego dysku.
Referencje
- Ou, X., Eilers, A.-C., Necib, L., Frebel, A. - MNRAS 528, 693, 2024.
- McMillan, P. J. - MNRAS 465, 76, 2017 - referencyjny model masy galaktycznej.
- Dutertre, X. - Bee Theory™ v2, BeeTheory.com, 2023.
- Freeman, K. C. - ApJ 160, 811, 1970.
- Bland-Hawthorn, J., Gerhard, O. - ARA&A 54, 529, 2016.
BeeTheory.com - Kwantowa grawitacja oparta na falach - 2025 - © Technoplane S.A.S.