BeeTheory – Statistische analyse – 2025

De pasvorm aflezen

χ²/dof = 0,475

Een getal tussen 0 en ∞ dat aangeeft of uw model goed, te goed of helemaal niet bij de gegevens past. Hier is wat het betekent in de galactische simulatie van BeeTheory – en wat we nodig zouden hebben om een echt solide resultaat te claimen.

Wat is χ²/dof?

Wanneer een model waarden _Vmodel(_Ri) voorspelt en de gegevens waargenomen waarden_Vobs(_Ri) met meetonzekerheden _σi opleveren, is de gereduceerde chi-kwadraat:

\(\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\)

waarbij N het aantal gegevenspunten is en p het aantal vrije parameters.

Hier: N = 16 Gaia 2024 punten, p = 2 passende parameters, K en α, dus dof = 14.

Elke term in de som is een trek: het residu uitgedrukt in eenheden van de meetonzekerheid.

Een pull van 0,5 betekent dat het model er een halve sigma naast zit – uitstekend. Een pull van 2,0 betekent een afwijking van twee sigma’s – het onderzoeken waard.

\(\frac{\chi^2}{\mathrm{dof}}{het model past bij de gegevens op het niveau van hun onzekerheden}\). \(\frac{\chi^2}{\mathrm{dof}}\ll1quad\Rightarrowquad\text{model past te goed: onzekerheden kunnen overschat zijn, of overfitting}\) \(\frac{\chi^2}{\mathrm{dof}{\gg1}quadraat}{het model klopt niet, of de onzekerheden zijn onderschat}\).

Ons getal: 0,475

0.475

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

Statistisch uitstekend

Het model is niet verkeerd

De onzekerheden kunnen enigszins groot zijn

χ²/dof schaal

0,475 ← ons

1,0 ideaal

00.511.522.53+

< 0.3
over-fit

0.3-0.7
uitstekend

0.7-1.3
goed

1.3-2.0
acceptabel

> 2.0
slecht

Een waarde van 0,475 ligt in de uitstekende zone. Het model is niet verkeerd – het gemiddelde residu is slechts 0,69σ over alle 16 gegevenspunten.

Natuurkundig gezien betekent dit dat BeeTheory de cirkelsnelheid tot op minder dan één meetonzekerheid nauwkeurig voorspelt bij vrijwel elke waargenomen straal.

Maar “uitstekend” is niet hetzelfde als “bewezen”.

χ²/dof = 0,475 zou ook kunnen betekenen dat de meetonzekerheden _σi iets overschat zijn. Als de werkelijke fouten 30% kleiner zouden zijn, zou hetzelfde model χ²/dof ≈ 1,0 geven.

We kunnen op basis van χ² alleen geen onderscheid maken tussen “model is erg goed” en “fouten zijn enigszins genereus”. Dit is een standaard statistische dubbelzinnigheid.

The Residuals – Punt voor punt

Het ruwe getal 0,475 verbergt informatie. Als u naar individuele pulls kijkt, ziet u de structuur van de fit: op welke punten het model spijkers met koppen slaat, en waar het moeite heeft.

Trekgrafiek: \((V_{\mathrm{model}}-V_{\mathrm{obs}})/\sigma_i\) voor elk Gaia 2024-gegevenspunt

R (kpc)	Vobs	σ	Vmodel	Overblijvend	Trek
Residu’s laden…

15 / 16 punten binnen 1σ

Elk gegevenspunt behalve R = 27,3 kpc heeft |pull| < 1,0. In een goed gekalibreerd model met Gaussische fouten verwachten we ongeveer 68% van de punten binnen 1σ - hier is dat 94%.

Dit suggereert dat het model uitstekend is of dat de foutbalken op de binnenste punten iets te groot zijn.

De uitschieter op R = 27,3 kpc – pull = +1,53σ

Het model voorspelt _Vc = 195,3 km/s, terwijl Gaia 173 ± 17 km/s meet.

De afwijking is 22,3 km/s, of 1,53σ. Dit is statistisch niet alarmerend, maar fysisch wel significant: het model neemt te langzaam af bij de grootste straal.

De best passende parameter K = 0,01349

De koppelingsconstante K = 0,01349 kpc-¹ is de amplitude van het BeeTheory donkere massaveld dat per eenheid zichtbare massa wordt opgewekt.

Het werd bepaald door χ² te minimaliseren over de 16 Gaia-gegevenspunten met α = 0,0744 kpc-¹, vastgesteld op basis van de gecombineerde Newby + Gaia fit.

\(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}\)

De waargenomen lokale donkere materiedichtheid is:

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

De dimensieloze koppeling λ = 2,44 ligt in het bereik λ ∈ [2, 7] dat is waargenomen in de pastings van BeeTheory, inclusief vergelijkingen voor één component, meerdere componenten en atomair H₂.

Deze universaliteit – dezelfde orde van grootte van femtometerschalen tot kiloparsecschalen – is de sterkste interne consistentiecontrole van het BeeTheory-raamwerk.

Waarom K hier kleiner is dan in pasvormen met één schijf

Eerdere fits die alleen de dunne schijf als bron gebruikten, gaven _Ksingle = 0,038 kpc-¹. Met alle zes galactische componenten die bijdragen aan het donkere veld, is de totale baryonische bron 2,8× groter.

Daarom moet K proportioneel kleiner zijn om dezelfde donkere massa te produceren.

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

Dit is precies wat de pasvorm geeft. Het is een consistentiecontrole, geen toeval.

Waarom 27,3 kpc moeilijk is voor elk model

Het buitenste datapunt van Gaia 2024, _Vc(27.3 kpc) = 173 ± 17 km/s, is niet alleen moeilijk voor BeeTheory. Het is het moeilijkste punt voor elk donkere-materiemodel dat wordt toegepast op de rotatiecurve van de Melkweg.

Model	Vc(27.3 kpc) voorspeld	Overblijvend	Trek	χ²/dof
Bijentheorie	195,3 km/s	+22.3	+1.53σ	0.475
NFW profiel	~198 km/s	+25	+1.5σ	0.44
Einasto α = 0,91	~191 km/s	+18	+1.1σ	0.38
Isotherme halo	~218 km/s	+45	+2.6σ	2.6

Geen enkel standaard twee-parametermodel reproduceert _Vc(27,3 kpc) = 173 km/s binnen 1σ.

Er zijn drie plausibele verklaringen.

De meting zelf: R = 27,3 kpc is het verste punt, met de minste kinematische tracers en de grootste systematische onzekerheden. Gaia DR4 kan de waarde verschuiven.
Een bijdrage van de gasschijf: De HI-schijf strekt zich uit tot ongeveer 25 kpc en draagt bij aan de baryonische massa bij grote stralen. Als deze als aparte component wordt meegenomen, kan de afname steiler worden.
Een kleinere coherentielengte: α = 0,12 kpc-¹, of ℓ = 8,3 kpc, past beter bij het buitenste punt, maar verslechtert het binnenste plateau.

Hoe een echt solide resultaat eruit zou zien

De huidige fit is statistisch goed. Maar goed in een model met 2 vrije parameters en 16 gegevenspunten is niet hetzelfde als bewezen.

0.475

Huidige χ²/dof, 16 punten, 2 parameters

< 0.8

Doel: goede pasvorm op uitgebreide dataset

~1.0

Ideaal: gekalibreerd voor monsters van melkwegstelsels

Vereisten voor een solide validatie

✓ Goede χ²/dof op de trainingsdataset: Bereikt. χ²/dof = 0,475 op Gaia 2024.
✓ Correcte lokale donkere materie dichtheid: Bereikt. 0,37 GeV/cm³ voorspeld vs 0,39 ± 0,03 waargenomen.
✓ Consistente dimensieloze koppeling: Bereikt. λ = 2,44 vergeleken met λ ≈ 3,5 van H₂ kalibratie.
! Reproduceer het buitenste Gaia-punt binnen 1σ: Niet volledig bereikt. Het residu op R = 27,3 kpc is +1,53σ.
✗ Blinde voorspelling op andere sterrenstelsels: Nog niet gedaan. De SPARC-catalogus biedt de natuurlijke test.
✗ Afleiding van ℓ uit eerste principes: Nog niet gedaan. ℓ is momenteel gepast, niet afgeleid.
✗ Validatie op clusterschaal: Nog niet gedaan. Bullet Cluster lensing is een belangrijke test.
→ Uitgebreide rotatiecurve voorbij 30 kpc: Gaia DR4 moet op korte termijn een kritische test opleveren.

De eerlijke wetenschappelijke status

De BeeTheory in zijn huidige vorm is een fysisch gemotiveerd fenomenologisch raamwerk dat zowel past bij de rotatiecurve van de Melkweg als bij de beste empirische donkere materieprofielen, met als voordeel dat het een fysisch mechanisme heeft.

Het is nog geen theorie in de strikte zin, omdat K en ℓ eerder gepast dan afgeleid zijn.

De weg naar een volledige theorie is duidelijk: leid ℓ = f(_Lsource) af van het golfmassapostulaat, test universeel op SPARC, en doe een blinde voorspelling voor Gaia DR4.

Gegevens: Ou, X. et al., MNRAS 528, 2024. Model: BeeTheory v2, Dutertre 2023, multi-component fit. SPARC referentie: Lelli, McGaugh, Schombert, AJ 152, 2016. Kogel cluster: Clowe et al., ApJL 648, 2006.