BeeTheory – Tilastollinen analyysi – 2025

Sovituksen lukeminen

χ²/dof = 0,475

0:n ja ∞:n välillä oleva luku, joka kertoo, sopiiko mallisi aineistoon hyvin, liian hyvin vai ei lainkaan. Seuraavassa kerrotaan, mitä se tarkoittaa BeeTheoryn galaktisessa simulaatiossa – ja mitä tarvitsisimme, jotta voisimme väittää, että tulos olisi todella luotettava.

Mikä on χ²/dof?

Kun malli ennustaa arvot _Vmodel(_Ri) ja data antaa havaitut arvot _Vobs(_Ri) mittausepävarmuuksilla _σi, redusoitu chi-neliö on:

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

jossa N on datapisteiden lukumäärä ja p on vapaiden parametrien lukumäärä.

Tässä: p = 2 sovitettua parametria, K ja α, joten dof = 14.

Summan jokainen termi on veto: jäännös ilmaistuna mittausepävarmuuden yksiköissä.

Vetovoima 0,5 tarkoittaa, että malli poikkeaa puoli sigmaa – erinomainen. Veto 2,0 tarkoittaa kahden sigman poikkeamaa – kannattaa tutkia.

\(\frac{\chi^2}{\mathrm{dof}}\approx1\quad\Rightarrow\quad\text{malli sopii dataan niiden epävarmuuksien tasolla}\) \(\frac{\chi^2}{\mathrm{dof}}\ll1\quad\Rightarrow\quad\text{malli sopii liian hyvin: epävarmuudet saattavat olla yliarvioituja eli overfitting}\) \(\frac{\chi^2}{\mathrm{dof}}\gg1\quad\Rightarrow\quad\text{malli on väärä tai epävarmuudet aliarvioidaan}\)

Meidän numeromme: 0.475

0.475

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

Tilastollisesti erinomainen

Malli ei ole väärässä

Epävarmuudet voivat olla hieman suuria

χ²/dof asteikko

0,475 ← us

1.0 ihanteellinen

00.511.522.53+

< 0.3
over-fit

0.3-0.7
erinomainen

0.7-1.3
hyvä

1.3-2.0
hyväksyttävä

> 2.0
huono

Arvo 0,475 sijoittuu erinomaiselle alueelle. Malli ei ole väärässä – keskimääräinen jäännös on vain 0,69σ kaikissa 16 datapisteessä.

Fyysisesti tämä tarkoittaa, että BeeTheory ennustaa ympyränopeuden alle yhden mittausepävarmuuden tarkkuudella lähes jokaisella havaitulla säteellä.

Mutta ”erinomainen” ei ole sama kuin ”todistettu”.

χ²/dof = 0,475 voi myös tarkoittaa, että mittausepävarmuudet _σi on hieman yliarvioitu. Jos todelliset virheet olisivat 30 % pienemmät, sama malli antaisi χ²/dof ≈ 1,0.

Pelkän χ²:n perusteella emme voi erottaa ”malli on erittäin hyvä” ja ”virheet ovat hieman anteliaita”. Tämä on tavanomainen tilastollinen epäselvyys.

The Residuals – Kohta kohdalta

Raakaluku 0,475 kätkee sisäänsä tietoa. Yksittäisten vetojen tarkastelu paljastaa sovituksen rakenteen: mitkä kohdat malli onnistuu ja missä se on vaikeuksissa.

Pull plot: \((V_{\mathrm{model}}-V_{\mathrm{obs}})/\sigma_i\) kullekin Gaia 2024 datapisteelle.

R (kpc)	Vobs	σ	Vmalli	Jäännös	Vedä
Jäännösten lataaminen…

15 / 16 pistettä 1σ:n sisällä

Jokaisen datapisteen paitsi R = 27,3 kpc:n kohdalla |pull| < 1,0. Hyvin kalibroidussa mallissa, jossa on Gaussin virheet, odotamme, että noin 68 % pisteistä on 1σ:n sisällä - tässä tapauksessa 94 %.

Tämä viittaa siihen, että malli on erinomainen tai sisäpisteiden virhepalkit ovat hieman liian suuret.

Poikkeus R = 27,3 kpc:n kohdalla – veto = +1,53σ.

Malli ennustaa _Vc = 195,3 km/s, kun taas Gaia mittaa 173 ± 17 km/s.

Ero on 22,3 km/s eli 1,53σ. Tämä ei ole tilastollisesti hälyttävää, mutta se on fysikaalisesti merkittävää: malli laskee liian hitaasti suurimmalla säteellä.

Parhaiten sopiva parametri K = 0,01349.

Kytkentävakio K = 0,01349 kpc-¹ on BeeTheoryn pimeän massakentän amplitudi, joka syntyy näkyvää massaa kohti.

Se määritettiin minimoimalla χ² 16 Gaia-datapisteen yli siten, että α = 0,0744 kpc-¹ on korjattu yhdistetystä Newby- ja Gaia-sovituksesta.

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

Havaittu paikallinen pimeän aineen tiheys on:

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

Mitaton kytkentä λ = 2,44 sijoittuu alueelle λ ∈ [2, 7], joka on havaittu kaikissa BeeTheory-sovituksissa, mukaan lukien yksikomponenttiset, monikomponenttiset ja atomiset H₂-vertailut.

Tämä yleispätevyys – sama suuruusluokka femtometrin mittakaavasta kiloparsec-mittakaavaan – on BeeTheory-kehyksen vahvin sisäisen johdonmukaisuuden tarkistus.

Miksi K on tässä pienempi kuin yhden levyn sovituksissa.

Aiemmat sovitukset, joissa lähteenä käytettiin vain ohutta kiekkoa, antoivat_Ksingle = 0,038 kpc-¹. Kun kaikki kuusi galaktista komponenttia osallistuvat pimeään kenttään, baryoninen kokonaislähde on 2,8 kertaa suurempi.

Siksi K:n on oltava suhteessa pienempi, jotta saadaan sama pimeä massa.

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

Tämä on juuri se, mitä sovitus antaa. Se on johdonmukaisuuden tarkistus, ei sattuma.

Miksi 27,3 kpc on vaikea jokaiselle mallille?

Kaikkein uloin Gaia 2024:n datapiste, _Vc(27,3 kpc) = 173 ± 17 km/s, ei ole vain BeeTheoryn kannalta vaikea. Se on vaikein piste jokaiselle Linnunradan kiertokäyrään sovelletulle pimeän aineen mallille.

Malli	Vc(27.3 kpc) ennustettu	Jäännös	Vedä	χ²/dof
BeeTheory	195,3 km/s	+22.3	+1.53σ	0.475
NFW-profiili	~198 km/s	+25	+1.5σ	0.44
Einasto α = 0,91	~191 km/s	+18	+1.1σ	0.38
Isoterminen halo	~218 km/s	+45	+2.6σ	2.6

Mikään kahden parametrin standardimalli ei toista _Vc(27,3 kpc) = 173 km/s 1σ:n tarkkuudella.

Tähän on kolme uskottavaa selitystä.

Itse mittaus: R = 27,3 kpc on kaukaisin piste, jossa on vähiten kinemaattisia merkkiaineita ja suurimmat systemaattiset epävarmuudet. Gaia DR4 saattaa siirtää arvoa.
Kaasukiekon osuus: HI-kiekko ulottuu noin 25 kpc:n etäisyydelle ja vaikuttaa baryoniseen massaan suurilla säteillä. Sen ottaminen mukaan erillisenä komponenttina voisi jyrkentää laskua.
Pienempi koherenssin pituus α = 0,12 kpc-¹ eli ℓ = 8,3 kpc sopii paremmin uloimpaan pisteeseen, mutta huonontaa sisemmän tasangon tilannetta.

Miltä aidosti vankka tulos näyttäisi?

Nykyinen sovitus on tilastollisesti hyvä. Mutta hyvä malli, jossa on 2 vapaata parametria ja 16 datapistettä, ei ole sama kuin todistettu.

0.475

Nykyinen χ²/dof, 16 pistettä, 2 muuttujaa

< 0.8

Tavoite: hyvä sovitus laajennettuun aineistoon

~1.0

Ihanteellinen: kalibroitu eri galaksinäytteissä.

Vaatimukset vankalle validoinnille

✓ Hyvä χ²/dof koulutustietokannassa: Saavutettu. χ²/dof = 0,475 Gaia 2024 -aineistossa.
✓ Oikea paikallinen pimeän aineen tiheys: Saavutettu. 0,37 GeV/cm³ ennustettu vs. 0,39 ± 0,03 havaittu.
✓ Johdonmukainen dimensioton kytkentä: λ = 2,44 verrattuna H₂-kalibroinnin λ ≈ 3,5:een.
! Toista uloin Gaia-piste 1σ:n tarkkuudella: Ei täysin saavutettu. Jäännös R = 27,3 kpc:n kohdalla on +1,53σ.
✗ Sokea ennustus muista galakseista: Ei vielä tehty. SPARC-katalogi tarjoaa luonnollisen testin.
✗ ℓ:n johtaminen ensimmäisistä periaatteista: Ei vielä tehty. ℓ on tällä hetkellä sovitettu, ei johdettu.
✗ Klusterin mittakaavan validointi: Ei vielä tehty. Bullet Cluster -linssit ovat keskeinen testi.
→ Laajennettu kiertokäyrä 30 kpc:n ulkopuolella: Gaia DR4:n pitäisi tarjota kriittinen testi lähitulevaisuudessa.

Rehellinen tieteellinen asema

BeeTeoria nykyisessä muodossaan on fysikaalisesti motivoitu fenomenologinen kehys, joka sopii Linnunradan pyörimisliikkeeseen sekä parhaisiin empiirisiin pimeän aineen profiileihin, ja sen etuna on fysikaalinen mekanismi.

Se ei ole vielä teoria varsinaisessa mielessä, koska K ja ℓ on pikemminkin sovitettu kuin johdettu.

Tie täydelliseen teoriaan on selvä: johdetaan ℓ = f(_Lsource) aaltomassapostulaatista, testataan universaalisti SPARC:lla ja tehdään sokea ennuste Gaia DR4:lle.

Tiedot: MNRAS 528, 2024. Malli: Dutertre 2023, monikomponenttisovitus. SPARC-viite: Lelli, McGaugh, Schombert, AJ 152, 2016. Bullet Cluster: Clowe et al., ApJL 648, 2006.