BeeTheory – Statistisk analys – 2025

Läsa av passformen

χ²/dof = 0,475

Ett tal mellan 0 och ∞ som talar om för dig om din modell passar data bra, för bra eller inte alls. Här är vad det betyder i BeeTheory galaktiska simulering – och vad vi skulle behöva för att hävda ett verkligt solidt resultat.

Vad är χ²/dof?

När en modell förutsäger värden _Vmodel(_Ri) och data ger observerade värden _Vobs(_Ri) med mätosäkerheter _σi, är det reducerade chi-två:

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

där N är antalet datapunkter och p är antalet fria parametrar.

Här: N = 16 Gaia 2024-punkter, p = 2 anpassade parametrar, K och α, så dof = 14.

Varje term i summan är en pull: restvärdet uttryckt i enheter av mätosäkerheten.

En dragning på 0,5 betyder att modellen är av med en halv sigma – utmärkt. En dragning på 2,0 innebär en avvikelse på två sigma – värt att undersöka.

\(\frac{\chi^2}{\mathrm{dof}}\approx1\quad\Rightarrow\quad\text{modellen passar data på nivån för deras osäkerheter}\) \(\frac{\chi^2}{\mathrm{dof}}\ll1\quad\Rightarrow\quad\text{modellen passar för bra: osäkerheterna kan vara överskattade, eller överanpassning}\) \(\frac{\chi^2}{\mathrm{dof}}\gg1\quad\Rightarrow\quad\text{modellen är fel, eller osäkerheterna är underskattade}\)

Vårt nummer: 0,475

0.475

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

Statistiskt sett utmärkt

Modellen är inte fel

Osäkerheterna kan vara något större

χ²/dof-skala

0,475 ← oss

1,0 ideal

00.511.522.53+

< 0.3
överanpassad

0.3-0.7
utmärkt

0.7-1.3
bra

1.3-2.0
acceptabel

> 2.0
dålig

Ett värde på 0,475 ligger i den utmärkta zonen. Modellen har inte fel – den genomsnittliga residualen är bara 0,69σ för alla 16 datapunkterna.

Fysikaliskt innebär detta att BeeTheory förutsäger den cirkulära hastigheten med mindre än en mätosäkerhet vid praktiskt taget varje observerad radie.

Men ”utmärkt” är inte samma sak som ”bevisat”

χ²/dof = 0,475 kan också betyda att mätosäkerheterna _σi är något överskattade. Om de verkliga felen var 30% mindre skulle samma modell ge χ²/dof ≈ 1,0.

Vi kan inte skilja mellan ”modellen är mycket bra” och ”felen är något generösa” från χ² ensam. Detta är en standard statistisk tvetydighet.

The Residuals – Punkt för punkt

Det obearbetade talet 0,475 döljer information. När man tittar på enskilda dragningar avslöjas strukturen i passformen: vilka punkter modellen spikar och var den kämpar.

Pull plot: \((V_{\mathrm{modell}}-V_{\mathrm{obs}})/\sigma_i\) för varje Gaia 2024-datapunkt

R (kpc)	Vobs	σ	Vmodell	Återstående	Dra
Laddar residualer…

15 / 16 poäng inom 1σ

Varje datapunkt utom R = 27,3 kpc har |pull| < 1,0. I en välkalibrerad modell med gaussiska fel förväntar vi oss cirka 68% av punkterna inom 1σ - här har vi 94%.

Detta tyder på att modellen är utmärkt eller att felstaplarna på de inre punkterna är något för stora.

Utropstecknet vid R = 27,3 kpc – pull = +1,53σ

Modellen förutspår _Vc = 195,3 km/s, medan Gaia mäter 173 ± 17 km/s.

Avvikelsen är 22,3 km/s, eller 1,53σ. Detta är inte statistiskt alarmerande, men det är fysiskt signifikant: modellen avtar för långsamt vid den största radien.

Den bäst anpassade parametern K = 0,01349

Kopplingskonstanten K = 0,01349 kpc-¹ är amplituden för BeeTheory-fältet för mörk massa som genereras per enhet synlig massa.

Den bestämdes genom att minimera χ² över de 16 Gaia-datapunkterna med α = 0,0744 kpc-¹ fixerad från den kombinerade Newby + Gaia-anpassningen.

\(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 observerade lokala tätheten av mörk materia är:

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

Den dimensionslösa kopplingen λ = 2,44 ligger i intervallet λ ∈ [2, 7] som observerats i BeeTheory-anpassningar, inklusive enkomponent-, flerkomponent- och atomära H₂-jämförelser.

Denna universalitet – samma storleksordning från femtometersskalan till kiloparsecskalan – är den starkaste interna konsistenskontrollen av BeeTheory-ramverket.

Varför K är mindre här än i enstaka diskpassningar

Tidigare anpassningar där endast den tunna skivan användes som källa gav _Ksingle = 0,038 kpc-¹. Med alla sex galaktiska komponenter som bidrar till det mörka fältet är den totala baryoniska källan 2,8× större.

Därför måste K vara proportionellt mindre för att ge samma mörka massa.

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

Detta är exakt vad passformen ger. Det är en konsekvenskontroll, inte en tillfällighet.

Varför 27,3 kpc är svårt för alla modeller

Den yttersta datapunkten för Gaia 2024, _Vc(27,3 kpc) = 173 ± 17 km/s, är inte bara svår för BeeTheory. Det är den svåraste punkten för alla modeller för mörk materia som tillämpas på Vintergatans rotationskurva.

Modell	Vc(27,3 kpc) förutspådd	Återstående	Dra	χ²/dof
BeeTheory	195,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 standardmodell med två parametrar återger _Vc(27,3 kpc) = 173 km/s inom 1σ.

Det finns tre rimliga förklaringar.

Själva mätningen: R = 27,3 kpc är den mest avlägsna punkten, med de minsta kinematiska spårämnena och de största systematiska osäkerheterna. Gaia DR4 kan förändra värdet.
Ett bidrag från en gasskiva: HI-disken sträcker sig till cirka 25 kpc och bidrar till den baryoniska massan vid stora radier. Om den inkluderas som en separat komponent kan nedgången bli brantare.
En mindre koherenslängd: α = 0,12 kpc-¹, eller ℓ = 8,3 kpc, passar den yttersta punkten bättre men försämrar den inre platån.

Hur ett verkligt solitt resultat skulle se ut

Den nuvarande anpassningen är statistiskt sett bra. Men bra i en modell med 2 fria parametrar och 16 datapunkter är inte samma sak som bevisat.

0.475

Aktuell χ²/dof, 16 poäng, 2 parametrar

< 0.8

Mål: bra passform på utökat dataset

~1.0

Idealisk: kalibrerad över galaxprover

Krav för en solid validering

✓ Bra χ²/dof på träningsdatasetet: Uppnått. χ²/dof = 0,475 på Gaia 2024.
✓ Korrekt lokal densitet för mörk materia: Uppnådd. 0,37 GeV/cm³ förutspådd jämfört med 0,39 ± 0,03 observerad.
✓ Konsistent dimensionslös koppling: Uppnådd. λ = 2,44 jämfört med λ ≈ 3,5 från H₂-kalibrering.
! Reproducera den yttersta Gaia-punkten inom 1σ: Inte helt uppnått. Restvärdet vid R = 27,3 kpc är +1,53σ.
✗ Blind förutsägelse om andra galaxer: Ännu inte gjort. SPARC-katalogen utgör det naturliga testet.
✗ Härledning av ℓ från första principen: Ännu inte gjort. ℓ är för närvarande monterad, inte härledd.
✗ Validering i klusterskala: Än är det inte klart. Bullet Cluster-linser är ett viktigt test.
→ Utökad rotationskurva bortom 30 kpc: Gaia DR4 bör ge ett kritiskt test på kort sikt.

Den ärliga vetenskapliga statusen

BeeTheory i sin nuvarande form är ett fysikaliskt motiverat fenomenologiskt ramverk som passar Vintergatans rotationskurva såväl som de bästa empiriska profilerna för mörk materia, med fördelen att ha en fysikalisk mekanism.

Det är ännu inte en teori i strikt mening, eftersom K och ℓ är anpassade snarare än härledda.

Vägen till en fullständig teori är tydlig: härled ℓ = f(_Lsource) från våg-massapostulatet, testa universellt på SPARC och gör en blind förutsägelse för Gaia DR4.

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