BeeTheory – Symulacja 4-składnikowa – SPARC 2025

Cztery komponenty geometryczne,
Jedno uniwersalne prawo:
20 galaktyk SPARC

Każdą galaktykę rozkładamy na cztery fizyczne komponenty – cienki dysk, gruby dysk, wybrzuszenie, pierścień gazowy – każdy z własną geometrią i skalą. Wszystkimi rządzi jedno prawo BeeTheory: $K_i = K_0/R_i$, $\ell_i = c\cdot R_i$. Jeden wolny parametr $K_0$ pasuje jednocześnie do 19 galaktyk.

0. Wyniki – Pierwszy

Dopasowanie 4-składnikowej teorii Bee – 20 galaktyk SPARC, dopasowanie K₀, c od Drogi Mlecznej

Dekompozycja każdej galaktyki na cienki dysk, gruby dysk, wybrzuszenie, jeśli jest obecne, i pierścień gazowy HI – każdy traktowany jako niezależne źródło BeeTheory z własną skalą i długością koherencji – daje 17/20 galaktyk w granicach 20% obserwowanej płaskiej prędkości rotacji $V_f$, z medianą błędu 7,4% dla 18 galaktyk rdzenia, z wyłączeniem dwóch strukturalnych wartości odstających CamB i NGC 3741.

Wynik jest bezpośrednio porównywalny z modelem 1-dyskowym, który daje 18/20, ale fizycznie bogatszy: pole ciemnej masy teraz poprawnie odzwierciedla geometrię każdego składnika barionowego.

Jedna uniwersalna stała $K_0 = 0,3759$ pozostała niezmieniona w stosunku do dopasowania 1-dyskowego. Dodanie trzech nowych źródeł geometrycznych nie wymaga ponownego dostrojenia sprzężenia fundamentalnego. Każdy składnik po prostu wnosi własne ciemne pole BeeTheory proporcjonalne do jego masy i odwrotnie do jego skali.

17 / 20W granicach20% $V_f$

7,4%Błądśredni, 18 rdzeni

$K_0 = 0,3759$Uniwersalny, bez zmian

4Składnikigeometryczne

1Wolnyparametr, K₀

2Strukturalnewartości odstające

Bezpośrednie porównanie: model 1-dyskowy vs 4-komponentowy

Kryterium	Model 1-dyskowy	Model 4-komponentowy	Werdykt
W granicach 20%	18 / 20	17 / 20	Porównywalny
Błąd mediany, rdzeń 18	6.8%	7.4%	Bardzo blisko
$K_0$	0.3759	0.3759, tak samo	Potwierdzony uniwersalny
Średni błąd galaktyk wybrzuszenia	-10.0%	-10.0%	Taki sam – model wybrzuszenia niewystarczający
Korekta bogata w gaz	Nieobecny	Uwzględniono, źródło pierścieniowe	Gaz ma teraz swój wkład
Rozkład fizyczny	Brak	Pełny, 4 składniki	Bardziej realistyczny
Darmowe parametry	1, $K_0$	1, $K_0$	Ta sama parsymonia

Kluczowe odkrycie – K₀ jest potwierdzonym uniwersalnym rozkładem

Fakt, że $K_0 = 0,3759$ jest identyczny zarówno w dopasowaniu 1-dyskowym, jak i 4-składnikowym – pomimo tego, że model 4-składnikowy zawiera trzy dodatkowe źródła – jest najsilniejszym wewnętrznym sprawdzianem spójności struktury BeeTheory.

Oznacza to, że pole ciemnej masy generowane na jednostkę masy jest naprawdę uniwersalne, niezależnie od tego, czy masa ta znajduje się w cienkim dysku młodych gwiazd, grubym dysku starych gwiazd, zwartym kulistym wybrzuszeniu czy pierścieniu gazu HI. Geometria, poprzez $R_i$, moduluje amplitudę pola; stała sprzężenia $K_0$ nie zmienia się.

1. Filozofia modelowania – jedno prawo, cztery geometrie

Głównym postulatem BeeTheory jest to, że każdy element masy $dV$ emituje pole falowe rozpadające się jako $e^{-\alpha D}/D^2$ w przestrzeni 3D. Amplituda sprzężenia i długość koherencji zależą od skali geometrycznej struktury źródła, a nie od rodzaju materii.

Uniwersalne prawo skalowania teorii pszczół – dotyczy identycznie wszystkich 4 składników $$\boxed{K_i = \frac{K_0}{R_i}, \qquad \ell_i = c\cdot R_i, \qquad \alpha_i = \frac{1}{\ell_i}}$$. $$K_0 = 0.3759\;(\text{dimensionless}), \qquad c = c_\text{disk}\text{ lub }c_\text{sph}\text{ według geometrii}$$

Współczynnik koherencji $c$ przyjmuje dwie wartości, określone na podstawie analizy dwustrefowej Drogi Mlecznej:

Dysk / Pierścień Źródła

$c_\text{disk} = 3,17$

Planarna geometria 2D
$\ell = 3,17 \times R_d$

Źródła sferyczne

$c_\text{sph} = 0,41$

Zwarta geometria 3D
$\ell = 0,41 \times r_b$

Współczynnik, stały

$c_\text{sph}/c_\text{disk} = 0,129$

Źródła sferyczne mają
$7,7 \times$ krótszą koherencję

Pochodzenie

Droga Mleczna

$\chi^2/\text{dof} = 0,24$
Kalibracja dwóch reżimów

2. Cztery składniki – wzory i skale

① Cienki dysk gwiezdny

Dysk wykładniczy – płaski 2D

Dominujący składnik gwiezdny. Zawiera młode gwiazdy, ramiona spiralne i Słońce. Modelowany jako dysk wykładniczy o promieniu skali $R_d$ bezpośrednio z fotometrii SPARC. Zawiera 75% masy gwiezdnej innej niż kulista.

Σ_thin(R) = Σ₀_thin - exp(-R/Rd)

K_thin = K₀/Rd, ℓ_thin = c_disk - Rd

Grubość dysku gwiezdnego ②

Dysk wykładniczy – 2D, bardziej rozciągnięty

Starsza, gorętsza kinematycznie populacja gwiazd. Bardziej rozciągnięty w pionie niż cienki dysk; w płaszczyźnie poziomej modelowany ze skalą $R_{d,\text{thick}} = 1,5R_d$ i 25% masy gwiezdnej.

Σ_thick(R) = Σ₀_thick - exp(-R / 1.5Rd)

K_thick = K₀/(1.5Rd), ℓ_thick = 1.5-c_disk-Rd

③ Wybrzuszenie, jeśli występuje

Sferyczny wykładniczy – kompaktowy 3D

Występuje tylko wtedy, gdy typ Hubble’a $T \leq 5$ i jest zidentyfikowany morfologicznie. Ułamek masowy $f_b(T)$ pochodzi ze standardowego rozkładu morfologicznego. Skala $r_b = 0,5R_d$. Używa $c_\text{sph} = 0.41$ – krótka spójność, intensywne pole wewnętrzne.

ρ_bulge(r) = ρ₀ - exp(-r / rb)

K_bulge = K₀/rb, ℓ_bulge = c_sph - rb

④ Dysk pierścienia gazowego HI

Profil pierścienia – 2D z centralnym otworem

Gazowy dysk HI ma centralną dziurę i rozciąga się do $R_\text{HI} \około 1,7R_d$. Jest on modelowany za pomocą profilu pierścieniowego $\Sigma \propto \exp(-R_m/R – R/R_\text{gas})$, tworząc centralny deficyt i naturalny pik. Masa gazu wynosi $M_\text{gas} = 1,33M_\text{HI}$, w tym helu.

Σ_gas(R) ∝ exp(-0.5-Rgas/R - R/Rgas)

K_gas = K₀/Rgas, ℓ_gas = c_disk - Rgas, Rgas = 1.7Rd

Poszczególne równania gęstości ciemnej

① Cienki dysk – gęstość ciemnej teorii pszczół $$\rho_\text{dark}^\text{thin}(r) = \frac{K_0}{R_d}\int_0^{8R_d}\!\Sigma_{0,\text{thin}}\,e^{-R’/R_d}\cdot\frac{(1+\alpha_d D)\,e^{-\alpha_d D}}{D^2}\cdot 2\pi R’\,dR’$$. $$\alpha_d = \frac{1}{c_\text{disk}\,R_d}, \quad D = \sqrt{r^2+R’^2}$$.

Gruby dysk – to samo jądro, inna skala $$\rho_\text{dark}^\text{thick}(r) = \frac{K_0}{1.5R_d}\int_0^{12R_d}\!\Sigma_{0,\text{thick}}\,e^{-R’/(1.5R_d)}\cdot\frac{(1+\alpha_k D)\,e^{-\alpha_k D}}{D^2}\cdot 2\pi R’\,dR’$$. $$\alpha_k = \frac{1}{c_\text{disk}\cdot 1.5R_d}$$.

③ Wybrzuszenie – źródło sferyczne 3D, krótsza koherencja $$\rho_\text{dark}^\text{bulge}(r) = \frac{K_0}{r_b}\int_0^{6r_b}\!\rho_{0,b}\,e^{-r’/r_b}\cdot\frac{(1+\alpha_b D)\,e^{-\alpha_b D}}{D^2}\cdot 4\pi r’^2\,dr’$$ $$\alpha_b = \frac{1}{c_\text{sph}\,r_b} = \frac{1}{0.41\,r_b}, \quad r_b = \max(0.5\,R_d,\;0.3\,\text{kpc}) $$

④ Pierścień gazowy – profil pierścienia z otworem centralnym $$\rho_\text{dark}^\text{gas}(r) = \frac{K_0}{R_\text{gas}}\int_0^{6R_\text{gas}}\!\Sigma_\text{gas}(R’)\cdot\frac{(1+\alpha_g D)\,e^{-\alpha_g D}}{D^2}\cdot 2\pi R’\,dR’$$. $$\Sigma_\text{gas}(R) \propto \exp\!\left(-\frac{0.5\,R_\text{gas}}{R} – \frac{R}{R_\text{gas}}\right), \quad R_\text{gas} = 1.7\,R_d$$

Całkowita gęstość ciemna – superpozycja wszystkich 4 składników $$\rho_\text{dark}(r) = \rho^\text{thin}(r) + \rho^\text{thick}(r) + \rho^\text{bulge}(r) + \rho^\text{gas}(r)$$ $$V_c(R) = \sqrt{\frac{G\!\left[M_\text{bar}(<R) + M_\text{dark}(<R)\right]}{R}}, \qquad M_\text{dark}(<R) = \int_0^R 4\pi r^2\,\rho_\text{dark}(r)\,dr$$.

3. Wszystkie parametry

$K_0$ sprzężenie uniwersalne0.3759

$c_\text{disk}$ dysk/pierścień3.17

$c_\text{sph}$ wybrzuszenie0.41

Frakcja cienkiego dysku $0.75(1-f_b)$

Ułamek grubego dysku$0.25(1-f_b)$

Skala grubego dysku$R_{d,\text{thick}} = 1.5R_d$

Skala wybrzuszenia$r_b= \max(0.5R_d,\,0.3\,\text{kpc})$

Skala pierścienia gazowego$R_\text{gas}= 1.7R_d$

Ułamek wybrzuszenia $f_b(T)$T≤0: 40%, T=2: 30%, T=3: 20%, T=4: 12%, T=5: 5%

Masa gazu$M_\text{gas} = 1.33M_\text{HI}$

$\Upsilon_\star$ gwiezdny M/L$0.5M_\odot/L_\odot$ przy 3.6 µm

Promień oceny$R_\text{eval}= 5R_d$

Jedyny wolny parametr

$K_0 = 0,3759$ to jedyny parametr dopasowany do danych SPARC, z wyłączeniem CamB. Wszystkie pozostałe wielkości – $c_\text{disk}$, $c_\text{sph}$, frakcje dyskowe i współczynniki skali – pochodzą z dwustopniowej kalibracji Drogi Mlecznej lub ze standardowych modeli populacji gwiazd. Model ma dokładnie 1 stopień swobody dla 19 galaktyk.

4. Przewidywania – wszystkie 20 galaktyk

Galaktyka	$R_d$	$f_b$	$f_\text{gas}$	$V_f$ obs	$V_\text{bar}$	$V_\text{ciemny}$	$V_\text{BT}$	Błąd	Status

$V_\text{BT}$ vs $V_f$ zaobserwowane – model 4-składnikowy

W granicach 20% 20-35% Wartości odstające Idealne 1:1 ±20%

Podział składowych ciemnej prędkości na galaktyki – skumulowany wkład $V_\text{dark}^2$

Ciemny cienki dysk Ciemny gruby dysk Ciemne wybrzuszenie Ciemny pierścień gazowy

5. Wnioski

Jakie 4 komponenty dowodzą, że 1 komponent nie może tego zrobić?

K₀ jest naprawdę uniwersalna. Stała sprzężenia nie zmienia się niezależnie od tego, czy źródło jest cienkim dyskiem młodych gwiazd, grubym dyskiem starych gwiazd, zwartym kulistym wybrzuszeniem, czy pierścieniem gazu HI. Jest to właściwość interakcji fala-masa, a nie typu składnika barionowego.

Pierścień gazowy generuje największe ciemne pole w galaktykach bogatych w gaz. W NGC 3621, gdzie $f_\text{gas} = 0,82$, pierścień gazowy przyczynia się do 68% całkowitej ciemnej prędkości – więcej niż dysk gwiezdny. BeeTheory prawidłowo przewiduje, że tam, gdzie znajdują się bariony, podąża za nimi ciemna masa, niezależnie od ich stanu fizycznego.

Pozostałe pozostałości wskazują na te same dwie przyczyny, co poprzednio. 7 galaktyk wybrzuszenia jest nadal niedoszacowanych średnio o około 10%, a dwie odstające galaktyki – CamB i NGC 3741 – wymagają niezależnego modelowania gazu za pomocą $R_\text{HI}$ z obserwacji radiowych, a nie skalowanego przybliżenia $1.7R_d$.

Dlaczego 4 składniki dają 17/20, podczas gdy 1 dysk daje 18/20?

Model 4-składnikowy wykorzystuje wartość Drogi Mlecznej $c_\text{disk} = 3,17$, podczas gdy zoptymalizowany model 1-dyskowy wykorzystuje $c = 6,40$, dopasowany na SPARC. Mniejsze $c$ oznacza krótszą długość koherencji i mniej ciemnego pola przy dużym $r$, co nieco zaniża przewidywania dla kilku galaktyk.

To napięcie między kalibracją Drogi Mlecznej a optymalnym modelem SPARC jest samo w sobie ważnym wynikiem naukowym: sugeruje, że $c$ może słabo zależeć od typu galaktyki lub że skalowanie gazu $1.7R_d$ niedoszacowuje rzeczywistego zasięgu gazu w bogatych w gaz galaktykach SPARC. Model 4-składnikowy jest fizycznie bardziej uczciwy, nawet jeśli numerycznie jest nieco mniej dokładny dla tej konkretnej metryki.

Dane: Lelli, F., McGaugh, S. S., Schombert, J. M., SPARC, AJ 152, 157 (2016). BeeTheory: Dutertre (2023), rozszerzona 2025. Frakcje wybrzuszenia: Moster et al. (2010), kalibracja morfologiczna. Stosunek HI/dysk gwiazdowy: Broeils & Rhee (1997), Lelli et al. (2014). Frakcja grubego dysku: Bland-Hawthorn & Gerhard (2016).