BeeTheory – Analiza statystyczna – 2025 r.

Odczytywanie dopasowania

χ²/dof = 0,475

Liczba z przedziału od 0 do ∞, która mówi Państwu, czy Państwa model pasuje do danych dobrze, zbyt dobrze lub wcale. Oto, co to oznacza w symulacji galaktycznej BeeTheory – i czego potrzebowalibyśmy, aby uzyskać naprawdę solidny wynik.

Co to jest χ²/dof?

Gdy model przewiduje wartości _Vmodel(_Ri), a dane dostarczają obserwowane wartości _Vobs(_Ri) z niepewnościami pomiarowymi _σi, zredukowany chi-kwadrat wynosi:

\(\frac{\chi^2}{\mathrm{dof}}=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_{\mathrm{model}}(R_i)-V_{\mathrm{obs}}(R_i)}{\sigma_i}\right)^2\)

gdzie N to liczba punktów danych, a p to liczba wolnych parametrów.

Tutaj: N = 16 punktów Gaia 2024, p = 2 dopasowane parametry, K i α, więc dof = 14.

Każdy człon sumy to ciąg: reszta wyrażona w jednostkach niepewności pomiaru.

Wartość 0,5 oznacza, że model różni się o pół sigma – doskonale. Przyciąganie 2,0 oznacza rozbieżność dwóch sigma – warte zbadania.

\(\frac{\chi^2}{\mathrm{dof}}\approx1\quad\Rightarrow\quad\text{model pasuje do danych na poziomie ich niepewności}\) \(\frac{\chi^2}{\mathrm{dof}}\ll1\quad\Rightarrow\quad\text\{model pasuje zbyt dobrze: niepewności mogą być przeszacowane lub nadmiernie dopasowane}\) \(\frac{\chi^2}{\mathrm{dof}}\gg1\quad\Rightarrow\quad\text{model jest błędny lub niepewności są niedoszacowane}\)

Nasz numer: 0,475

0.475

χ² / dof – N = 16 – p = 2 – dof = 14

Statystycznie doskonały

Model nie jest błędny

Niepewności mogą być nieznacznie duże

χ²/dof skala

0.475 ← us

1.0 idealny

00.511.522.53+

< 0.3
nadmierne dopasowanie

0.3-0.7
doskonały

0.7-1.3
dobry

1.3-2.0
akceptowalny

> 2.0
słaby

Wartość 0,475 znajduje się w strefie doskonałej. Model nie jest błędny – średnia wartość resztowa wynosi tylko 0,69σ dla wszystkich 16 punktów danych.

Fizycznie oznacza to, że BeeTheory przewiduje prędkość kołową z dokładnością mniejszą niż jedna niepewność pomiarowa dla praktycznie każdego obserwowanego promienia.

Ale „doskonały” to nie to samo co „sprawdzony”

χ²/dof = 0,475 może również oznaczać, że niepewności pomiarowe _σi są nieco zawyżone. Gdyby prawdziwe błędy były o 30% mniejsze, ten sam model dałby χ²/dof ≈ 1,0.

Nie możemy odróżnić „model jest bardzo dobry” od „błędy są nieco hojne” na podstawie samego χ². Jest to standardowa niejednoznaczność statystyczna.

The Residuals – Punkt po punkcie

Surowa liczba 0,475 ukrywa informacje. Spojrzenie na poszczególne podciągnięcia ujawnia strukturę dopasowania: które punkty model przybija, a gdzie ma trudności.

Wykres ciągniony: \((V_{\mathrm{model}}-V_{\mathrm{obs}})/\sigma_i\) dla każdego punktu danych Gaia 2024.

R (kpc)	Vobs	σ	Vmodel	Pozostały	Proszę ciągnąć
Ładowanie reszt…

15 / 16 punktów w granicach 1σ

Każdy punkt danych z wyjątkiem R = 27,3 kpc ma |pull| < 1,0. W dobrze skalibrowanym modelu z błędami gaussowskimi spodziewamy się około 68% punktów w granicach 1σ - tutaj mamy 94%.

Sugeruje to, że model jest doskonały lub słupki błędów w punktach wewnętrznych są nieco zbyt duże.

Wartość odstająca przy R = 27,3 kpc – pull = +1,53σ

Model przewiduje _Vc = 195,3 km/s, podczas gdy Gaia mierzy 173 ± 17 km/s.

Rozbieżność wynosi 22,3 km/s, czyli 1,53σ. Nie jest to statystycznie alarmujące, ale ma znaczenie fizyczne: model spada zbyt wolno przy największym promieniu.

Najlepiej dopasowany parametr K = 0,01349

Stała sprzężenia K = 0,01349 kpc-¹ to amplituda pola ciemnej masy BeeTheory generowana na jednostkę masy widzialnej.

Została ona wyznaczona poprzez minimalizację χ² w 16 punktach danych Gaia z α = 0.0744 kpc-¹ ustalonym z połączonego dopasowania Newby + Gaia.

\(K=0.01349\,\mathrm{kpc}^{-1}\) \(\lambda=K\ell^2=\frac{K}{\alpha^2}=\frac{0.01349}{(0.0744)^2}=2.44\) \(\rho_{\mathrm{dark}}(R_\odot=8\,\mathrm{kpc})=0.37\,\mathrm{GeV/cm^3}\)

Obserwowana lokalna gęstość ciemnej materii wynosi:

\(\rho_{\mathrm{obs}}(R_\odot)=0.39\pm0.03\,\mathrm{GeV/cm^3}\)

Bezwymiarowe sprzężenie λ = 2,44 mieści się w zakresie λ ∈ [2, 7] obserwowanym w dopasowaniach BeeTheory, w tym w porównaniach jednoskładnikowych, wieloskładnikowych i atomowych H₂.

Ta uniwersalność – ten sam rząd wielkości od skali femtometrów do skali kiloparseków – jest najsilniejszym wewnętrznym sprawdzianem spójności ram BeeTheory.

Dlaczego K jest tutaj mniejsze niż w dopasowaniach jednodyskowych?

Poprzednie dopasowania wykorzystujące tylko cienki dysk jako źródło dały_Ksingle = 0,038 kpc-¹. Przy wszystkich sześciu galaktycznych składnikach przyczyniających się do ciemnego pola, całkowite źródło barionowe jest 2,8 razy większe.

Dlatego K musi być proporcjonalnie mniejsze, aby wytworzyć taką samą ciemną masę.

\(K_{\mathrm{multi}}\approx\frac{K_{\mathrm{single}}}{2.8}=\frac{0.038}{2.8}=0.0136\,\mathrm{kpc}^{-1}\)

To jest dokładnie to, co daje dopasowanie. Jest to sprawdzenie spójności, a nie zbieg okoliczności.

Dlaczego 27,3 kpc jest trudne dla każdego modelu?

Najbardziej oddalony punkt danych Gaia 2024, _Vc(27.3 kpc) = 173 ± 17 km/s, jest trudny nie tylko dla BeeTheory. Jest to najtrudniejszy punkt dla każdego modelu ciemnej materii zastosowanego do krzywej rotacji Drogi Mlecznej.

Model	Przewidywane Vc(27,3 kpc)	Pozostały	Proszę ciągnąć	χ²/dof
BeeTheory	195,3 km/s	+22.3	+1.53σ	0.475
Profil NFW	~198 km/s	+25	+1.5σ	0.44
Einasto α = 0,91	~191 km/s	+18	+1.1σ	0.38
Aureola izotermiczna	~218 km/s	+45	+2.6σ	2.6

Żaden standardowy model dwuparametrowy nie odtwarza _Vc(27.3 kpc) = 173 km/s z dokładnością do 1σ.

Istnieją trzy prawdopodobne wyjaśnienia.

Sam pomiar: R = 27,3 kpc to najbardziej odległy punkt, z najmniejszą liczbą znaczników kinematycznych i największą niepewnością systematyczną. Gaia DR4 może przesunąć tę wartość.
Wkład dysku gazowego: Dysk HI rozciąga się do około 25 kpc i przyczynia się do masy barionowej na dużych promieniach. Uwzględnienie go jako osobnego składnika mogłoby spowodować stromy spadek.
Mniejsza długość koherencji: α = 0,12 kpc-¹ lub ℓ = 8,3 kpc, lepiej pasuje do najbardziej zewnętrznego punktu, ale pogarsza wewnętrzne plateau.

Jak wyglądałby naprawdę solidny wynik?

Obecne dopasowanie jest statystycznie dobre. Ale dobre w modelu z 2 wolnymi parametrami i 16 punktami danych to nie to samo, co udowodnione.

0.475

Prąd χ²/dof, 16 punktów, 2 parametry

< 0.8

Cel: dobre dopasowanie do rozszerzonego zbioru danych

~1.0

Idealny: skalibrowany w próbkach galaktyk

Wymagania dotyczące solidnej walidacji

Dobre χ²/dof na zestawie danych treningowych: Osiągnięto. χ²/dof = 0,475 na Gaia 2024.
Prawidłowa lokalna gęstość ciemnej materii: Osiągnięto. 0,37 GeV/cm³ przewidywane vs 0,39 ± 0,03 obserwowane.
Spójne bezwymiarowe sprzężenie: Osiągnięto. λ = 2,44 w porównaniu z λ ≈ 3,5 z kalibracji H₂.
! Proszę odtworzyć najbardziej oddalony punkt Gaia z dokładnością do 1σ: Nie w pełni osiągnięte. Pozostałość przy R = 27,3 kpc wynosi +1,53σ.
Ślepe przewidywanie dotyczące innych galaktyk: Jeszcze nie zrobione. Katalog SPARC zapewnia naturalny test.
✗ Wyprowadzenie ℓ z pierwszych zasad: Jeszcze nie zrobione. ℓ jest obecnie dopasowany, a nie wyprowadzony.
Walidacja w skali klastra: Nie została jeszcze przeprowadzona. Soczewkowanie gromady kulistej jest kluczowym testem.
→ Rozszerzona krzywa rotacji poza 30 kpc: Gaia DR4 powinna zapewnić krytyczny test w najbliższej przyszłości.

Uczciwy status naukowy

BeeTheory w swojej obecnej formie jest fizycznie umotywowaną strukturą fenomenologiczną, która pasuje do krzywej rotacji Drogi Mlecznej, a także do najlepszych empirycznych profili ciemnej materii, z tą zaletą, że ma fizyczny mechanizm.

Nie jest to jeszcze teoria w ścisłym sensie, ponieważ K i ℓ są raczej dopasowane niż wyprowadzone.

Droga do pełnej teorii jest jasna: wyprowadzić ℓ = f(_Lsource) z postulatu masy falowej, przetestować uniwersalnie na SPARC i dokonać ślepej prognozy dla Gaia DR4.

Dane: Ou, X. et al., MNRAS 528, 2024. Model: BeeTheory v2, Dutertre 2023, dopasowanie wieloskładnikowe. Referencja SPARC: Lelli, McGaugh, Schombert, AJ 152, 2016. Bullet Cluster: Clowe et al., ApJL 648, 2006.