BeeTheory – Statistisk analyse – 2025

Læsning af pasformen

χ²/dof = 0,475

Et tal mellem 0 og ∞, der fortæller dig, om din model passer godt, for godt eller slet ikke til dataene. Her er, hvad det betyder i BeeTheorys galaktiske simulering – og hvad vi skal bruge for at få et virkelig solidt resultat.

Hvad er χ²/dof?

Når en model forudsiger værdier Vmodel(Ri), og dataene giver observerede værdier Vobs(Ri) med måleusikkerheder σi, er det reducerede chi-kvadrat:

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

hvor N er antallet af datapunkter, og p er antallet af frie parametre.

I dette tilfælde: N = 16 Gaia 2024-punkter, p = 2 tilpassede parametre, K og α, så dof = 14.

Hvert udtryk i summen er et pull: restværdien udtrykt i enheder af måleusikkerheden.

Et pull på 0,5 betyder, at modellen afviger med en halv sigma – fremragende. Et pull på 2,0 betyder en uoverensstemmelse på to sigma – det er værd at undersøge.

\(\frac{\chi^2}{\mathrm{dof}}\approx1\quad\Rightarrow\quad\text{modellen passer til dataene på niveau med deres usikkerheder}\) \(\frac{\chi^2}{\mathrm{dof}}\ll1\quad\Rightarrow\quad\text{modellen passer for godt: usikkerhederne kan være overvurderet eller overtilpasning}\) \(\frac{\chi^2}{\mathrm{dof}}\gg1\quad\Rightarrow\quad\text{modellen er forkert, eller usikkerhederne er undervurderet}\)

Vores nummer: 0,475

0.475

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

χ²/dof-skala

0,475 ← us
1,0 ideel
00.511.522.53+
< 0.3
overtilpasset
0.3-0.7
fremragende
0.7-1.3
god
1.3-2.0
acceptabel
> 2.0
dårlig

En værdi på 0,475 ligger i den fremragende zone. Modellen tager ikke fejl – den gennemsnitlige residual er kun 0,69σ på tværs af alle 16 datapunkter.

Fysisk betyder det, at BeeTheory forudsiger den cirkulære hastighed med mindre end én måleusikkerhed ved stort set alle observerede radier.

Men “fremragende” er ikke det samme som “bevist”

χ²/dof = 0,475 kan også betyde, at måleusikkerhederne σi er lidt overvurderede. Hvis de sande fejl var 30 % mindre, ville den samme model give χ²/dof ≈ 1,0.

Vi kan ikke skelne mellem “modellen er meget god” og “fejlene er lidt for store” ud fra χ² alene. Dette er en standard statistisk tvetydighed.

The Residuals – punkt for punkt

Det rå tal 0,475 skjuler information. Når man ser på individuelle træk, afsløres strukturen i tilpasningen: hvilke punkter modellen klarer, og hvor den kæmper.

Pull-plot: \((V_{\mathrm{model}}-V_{\mathrm{obs}})/\sigma_i\) for hvert Gaia 2024-datapunkt

R (kpc)VobsσVmodelResterendeTræk
Indlæsning af residualer…

15 / 16 point inden for 1σ

Alle datapunkter undtagen R = 27,3 kpc har |pull| < 1,0. I en velkalibreret model med gaussiske fejl forventer vi ca. 68 % af punkterne inden for 1σ - her har vi 94 %.

Det tyder på, at modellen er fremragende, eller at fejlbjælkerne på de indre punkter er lidt for store.

Afvigelsen ved R = 27,3 kpc – pull = +1,53σ

Modellen forudsiger Vc = 195,3 km/s, mens Gaia måler 173 ± 17 km/s.

Afvigelsen er 22,3 km/s eller 1,53σ. Det er ikke statistisk alarmerende, men det er fysisk signifikant: Modellen falder for langsomt ved den største radius.

Den bedst tilpassede parameter K = 0,01349

Koblingskonstanten K = 0,01349 kpc-¹ er amplituden af det mørke BeeTheory-massefelt, der genereres pr. enhed synlig masse.

Den blev bestemt ved at minimere χ² over de 16 Gaia-datapunkter med α = 0,0744 kpc-¹ fastsat ud fra den kombinerede Newby + Gaia-tilpasning.

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

Den observerede lokale tæthed af mørkt stof er:

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

Den dimensionsløse kobling λ = 2,44 ligger i intervallet λ ∈ [2, 7], der er observeret på tværs af BeeTheory-tilpasninger, herunder enkeltkomponent-, multikomponent- og atomare H₂-sammenligninger.

Denne universalitet – samme størrelsesorden fra femtometerskalaer til kiloparsecskalaer – er den stærkeste interne konsistenskontrol af BeeTheory-rammen.

Hvorfor K er mindre her end i single-disk fits

Tidligere tilpasninger, der kun brugte den tynde skive som kilde, gav Ksingle = 0,038 kpc-¹. Med alle seks galaktiske komponenter, der bidrager til det mørke felt, er den samlede baryoniske kilde 2,8× større.

Derfor skal K være proportionalt mindre for at producere den samme mørke masse.

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

Det er præcis, hvad tilpasningen giver. Det er et konsistenstjek, ikke en tilfældighed.

Hvorfor 27,3 kpc er svært for alle modeller

Det yderste datapunkt i Gaia 2024, Vc(27,3 kpc) = 173 ± 17 km/s, er ikke kun vanskeligt for BeeTheory. Det er det sværeste punkt for alle modeller for mørkt stof, der anvendes på Mælkevejens rotationskurve.

ModelVc(27,3 kpc) forudsagtResterendeTrækχ²/dof
BeeTheory195,3 km/s+22.3+1.53σ0.475
NFW-profil~198 km/s+25+1.5σ0.44
Einasto α = 0,91~191 km/s+18+1.1σ0.38
Isotermisk halo~218 km/s+45+2.6σ2.6

Ingen standardmodel med to parametre gengiver Vc(27,3 kpc) = 173 km/s inden for 1σ.

Der er tre plausible forklaringer.

  • Selve målingen: R = 27,3 kpc er det fjerneste punkt med de færreste kinematiske sporstoffer og de største systematiske usikkerheder. Gaia DR4 kan flytte værdien.
  • Et bidrag fra gasskiven: HI-disken strækker sig til omkring 25 kpc og bidrager til den baryoniske masse ved store radier. Hvis man inkluderer den som en separat komponent, kan faldet blive stejlere.
  • En mindre kohærenslængde: α = 0,12 kpc-¹, eller ℓ = 8,3 kpc, passer bedre til det yderste punkt, men forværrer det indre plateau.

Sådan ville et virkelig solidt resultat se ud

Den nuværende tilpasning er statistisk god. Men godt i en model med 2 frie parametre og 16 datapunkter er ikke det samme som bevist.

Krav til en solid validering

  • ✓ God χ²/dof på træningsdatasættet: Opnået. χ²/dof = 0,475 på Gaia 2024.
  • ✓ Korrekt lokal tæthed af mørkt stof: Opnået. 0,37 GeV/cm³ forudsagt mod 0,39 ± 0,03 observeret.
  • ✓ Konsistent dimensionsløs kobling: Opnået. λ = 2,44 sammenlignet med λ ≈ 3,5 fra H₂-kalibrering.
  • ! Gengiv det yderste Gaia-punkt inden for 1σ: Ikke fuldt ud opnået. Restværdien ved R = 27,3 kpc er +1,53σ.
  • ✗ Blind forudsigelse af andre galakser: Det er ikke gjort endnu. SPARC-kataloget er den naturlige test.
  • ✗ Udledning af ℓ fra de første principper: Endnu ikke gjort. ℓ er i øjeblikket tilpasset, ikke udledt.
  • ✗ Validering i klyngeskala: Det er ikke gjort endnu. Bullet Cluster-linser er en vigtig test.
  • → Udvidet rotationskurve ud over 30 kpc: Gaia DR4 bør give en kritisk test på kort sigt.

Den ærlige videnskabelige status

BeeTheory i sin nuværende form er en fysisk motiveret fænomenologisk ramme, som passer til Mælkevejens rotationskurve og de bedste empiriske profiler af mørkt stof, og som har den fordel, at den har en fysisk mekanisme.

Det er endnu ikke en teori i streng forstand, fordi K og ℓ er tilpasset snarere end udledt.

Vejen til en fuld teori er klar: udled ℓ = f(Lsource) fra bølgemassepostulatet, test universelt på SPARC, og lav en blind forudsigelse for Gaia DR4.

Data: Ou, X. et al, MNRAS 528, 2024. Model: BeeTheory v2, Dutertre 2023, multikomponenttilpasning. SPARC-reference: Lelli, McGaugh, Schombert, AJ 152, 2016. Bullet Cluster: Clowe et al, ApJL 648, 2006.

BeeTheory.com – Bølgebaseret kvantegravitation – © Technoplane S.A.S. 2025