BeeTheory – SPARC Galaxies – Justeret tilpasning – 2025

BeeTheory anvendt på 20 eksterne galakser:
Justeret formel og blindtestmetode

SPARC-kataloget indeholder 175 galakser med målte baryoniske profiler og rotationskurver. Vi anvender BeeTheorys ligning for mørk masse – og justerer dens skaleringslov, så den passer til galaksepopulationen – og rapporterer resultatet: 18 ud af 20 galakser forudsiges inden for 20 % af deres observerede flade rotationshastighed, med χ²/dof = 0,93.

0. Resultater – sagt først

Bedste tilpasning – 20 SPARC-galakser, Q = 1, høj kvalitet

Med den modificerede BeeTheory-formel Kd = K0/Rd og ℓd = c –Rd passer to universelle konstanter til alle 20 galakser samtidigt.

Den mørke massetæthed i hver galakse forudsiges alene ud fra dens baryoniske diskparametre – ingen tuning pr. galakse.

De bedst tilpassede parametre: K0 = 0,3759, dimensionsløs, og c = 6,40, dimensionsløs. Resultat: 18/20 galakser inden for 20 % af den observeredeVf, χ²/dof = 0,93. To afvigere, CamB og NGC 3741, er gasdominerede dværge, hvor modelleringen af stjerneskiven bryder sammen.

18/20

Inden for 20 % afVf

6.8%

Medianfejl

0.93

χ²/dof

2

Universelle konstanter

0.91

Pearson r, Tully-Fisher

Modificeret BeeTheory-formel – justeret for galaksepopulation \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d}\int_0^{R_{\mathrm{max}}}\Sigma_0 e^{-R’/R_d}\frac{(1+\alpha_d D)e^{-\alpha_d D}}{D^2}\,2\pi R’\,dR’\) \(D=\sqrt{r^2+R’^2},\qquad \alpha_d=\frac{1}{cR_d},\qquad K_0=0.3759,\qquad c=6.40\)

1. De 20 SPARC-galakser – data og forudsigelser

SPARC-samlingen dækker galakser, der spænder over fem årtier i luminositet, fra irregulære dværggalakser til massive spiraler. For hver galakse er inputparametrene taget direkte fra Lelli et al. 2016 Tabel 1: diskens skalaradiusRd, den centrale overfladelysstyrke Σd, HI-gasmassenMHI og den flade rotationshastighedVf.

Stjernemassen er beregnet som M★ = Υ★ × L3.6, med Υ★ = 0,5 M⊙/L⊙. Gasmassen er beregnet som Mgas = 1,33 ×MHI.

Galakse Rd kpc Kd kpc-¹ ℓd kpc Vf obs Vbar Vdark VBT Fejl Status
Indlæser galakse-tabel…

Alle hastigheder i km/s. Fejl = (VBTVf)/Vf. Parametre: Kd = 0,3759/Rd, ℓd = 6,40 ×Rd. BeeTheory-forudsigelse evalueret vedReval = 5Rd.

VBT vsVf observeret – 20 SPARC-galakser, justeret BeeTheory
Inden for 20 %, 18 galakser Afvigelser: CamB, NGC 3741 Perfekt forudsigelse, 1:1 ±20% indhylling

2. Den modificerede formel – hvorfor K ∝ 1/Rd

Den oprindelige BeeTheory-mælkevejstilpasning brugte en enkelt koblingskonstant K = 0,02365 kpc-¹ med kohærenslængde ℓ = 3,17Rd. Når den blev anvendt blindt på 20 SPARC-galakser, undervurderede den systematiskVf med ca. 50 %.

Analysen pr. galakse afslørede et klart mønster: Den nødvendige koblingskonstant varierer som K ∝ 1/Rd.

2.1 Fra en konstant til en skaleringslov

Den vigtigste indsigt er dimensionel. Den mørke tæthed i BeeTheory ved radius r fra en skive med skalaen Rd og overfladetætheden Σ0 er i det asymptotiske regime med flad rotation r ≪ ℓ:

Asymptotisk mørk tæthed,Rd ≪ r ≪ ℓ \(\rho_{\mathrm{dark}}(r)\approx \frac{K\Sigma_0 R_d^2}{r^2}f\left(\frac{r}{\ell}\right)\) \(M_{\mathrm{dark}}(<r)\propto K\Sigma_0R_d^2r\)

Den flade rotationshastighed skaleres derefter som:

Flad rotationshastighed fra BeeTheory mørk masse \(V_f^2=\frac{GM_{\mathrm{dark}}(<R_{\mathrm{eval}})}{R_{\mathrm{eval}}}\propto GK\Sigma_0R_d^2=GK\frac{M_\star}{2\pi}\) \(V_f\propto (K M_\star)^{1/2}\)

Den observerede baryoniske Tully-Fisher-relation siger Vf4Mbar, hvilket betyderVfMbar1/4. For at dette kan gengives af BeeTheory, har vi brug for Vf2M★/Rd, den gennemsnitlige overfladetæthed på disken. Dette kræver:

Nødvendig skalering for at gengive Tully-Fisher-hældningen \(V_f^2\propto GK\frac{M_\star}{2\pi}\propto\frac{M_\star}{R_d}\implies K\propto\frac{1}{R_d}\) \(\boxed{K_d=\frac{K_0}{R_d},\qquad K_0=0.3759}\)

Denne skalering er ikke en ad hoc-løsning – det er, hvad Tully-Fisher-relationen kræver. En kobling K ∝ 1/Rd betyder, at mere kompakte skiver genererer et stærkere mørkt felt pr. masseenhed.

2.2 Kohærenslængden – hvorfor c = 6,40 ≠ 3,17

Mælkevejstilpasningen gav cMW = ℓd/Rd = 3,17. SPARC-prøven giver cSPARC = 6,40, cirka dobbelt så stor. Der er to mulige forklaringer:

Begge effekter er reelle og målbare. Den endelige værdi af c kræver, at gas- og stjerneskiver modelleres separat.

Komplet modificeret BeeTheory-formel – to universelle konstanter \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d}\int_0^{8R_d}\Sigma_0 e^{-R’/R_d}\frac{(1+\alpha_dD)e^{-\alpha_dD}}{D^2}\,2\pi R’\,dR’\) \(\alpha_d=\frac{1}{cR_d},\qquad D=\sqrt{r^2+R’^2},\qquad K_0=0.3759,\qquad c=6.40\)

3. Beregningen – trin for trin

For hver SPARC-galakse foregår BeeTheory-forudsigelsen i fem trin. Ingen frie parametre justeres pr. galakse.

1
Læs baryoniske input fra SPARC tabel 1

Rd, Σd,MHI ogVf. Konverter Σ0 = Σd × Υ★ × 10⁶ M⊙/kpc², og Mgas = 1,33 ×MHI.

2
Beregn BeeTheory-parametre fraRd

Kd = K0/Rd = 0,3759/Rd, ℓd = cRd = 6,40Rd, og αd = 1/ℓd. Ingen tilpasning.

3
Integrer den mørke tæthed ved r = 5Rd
\(\rho_{\mathrm{dark}}(5R_d)=K_d\sum_{i=1}^{60}\Sigma_0e^{-R_i’/R_d}\frac{(1+\alpha_dD_i)e^{-\alpha_dD_i}}{D_i^2}2\pi R_i’\Delta R’\)

Numerisk integration med 60 ringe, R′ ∈ [0, 8Rd]. Integrer derefter sfærisk for at få den indesluttede mørke masse Mdark(<5Rd).

4
Beregn baryonisk cirkulær hastighed
\(V_{\mathrm{bar}}(R)=\sqrt{\frac{G[M_{\mathrm{disk}}(<R)+M_{\mathrm{gas}}(<R)+M_{\mathrm{bulge}}(<R)]}{R}}\)
5
Forudsig den samlede cirkulære hastighed
\(V_{\mathrm{BT}}=\sqrt{V_{\mathrm{bar}}^2+V_{\mathrm{dark}}^2},\qquad V_{\mathrm{dark}}=\sqrt{\frac{GM_{\mathrm{dark}}(<5R_d)}{5R_d}}\)

Sammenlign med observeretVf. Fejl = (VBTVf)/Vf.

Vdark/Vbar-forhold – hvor dominerende er den mørke komponent?

4. Hvorfor blindtest er den eneste ærlige test

En model, der gengiver de data, den blev kalibreret på, beviser intet. Alle modeller, selv de forkerte, kan justeres, så de passer til træningsdataene. Den eneste videnskabeligt meningsfulde test er en blind forudsigelse: Anvend modellen på data, den aldrig har set, med parametre, der er fastfrosset fra kalibreringen, og rapporter resultatet – uanset hvad det er.

4.1 Hvad “blind” betyder her

BeeTheory-parametrene K0 = 0,3759 og c = 6,40 blev bestemt ved at tilpasse de 20 SPARC-galakser samtidigt. De er nu faste.

Blindtesten ville være: at anvende disse parametre på de resterende 155 SPARC-galakser, som ikke blev brugt i tilpasningen, og rapportere resultatet, før man ser på deres observerede rotationskurver. Denne test er endnu ikke udført – det er næste skridt.

De oprindelige Mælkevejsparametre, Kd = 0,02365 og ℓd = 3,17Rd, blev bestemt på en enkelt galakse. At anvende dem på SPARC uden justering gav 0/20 korrekte galakser – en ærlig og vigtig fejl. Denne fejl afslørede K ∝ 1/Rd-skaleringen.

4.2 Statistisk betydning af tilpasningskvalitet

Med χ²/dof = 0,93 på tværs af 20 galakser passer modellen nogenlunde til det forventede niveau for de antagne hastighedsusikkerheder på 15 %.

Chi-kvadrat fortolkning \(\frac{\chi^2}{\mathrm{dof}}=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_{\mathrm{BT}}(i)-V_f(i)}{0.15V_f(i)}\right)^2=0.93\) \(N=19\ \text{(undtagen CamB)},\qquad p=2\ (K_0,c),\qquad \mathrm{dof}=17\)

En værdi på 0,93 er meget tæt på idealet 1,0. Modellen tager højde for spredningen på niveauet for måleusikkerhed.

4.3 De to afvigende værdier

CamB – ren gasdværg,Vf = 2,0 km/s

CamB har næsten ingen stjernemasse, M★ ≈ 2×10⁷ M⊙. BeeTheory-formlen bruger Σ0e-R/Rd som kilde – men i CamB er baryonerne næsten udelukkende HI-gas, ikke stjerner. Stjerneskivemodellen kan ikke anvendes.

NGC 3741 – overvurderet med 47 %.

NGC 3741 er en lille dværg med lav overfladelysstyrke og en meget udstrakt HI-skive. BeeTheory-kilden, stjerneskiven, undervurderer den faktiske baryoniske udstrækning. Hvis man inkluderede gasskiven som en separat kildekomponent med større skalaradius, ville det reducere den forudsagte mørke masse og korrigere overvurderingen.

De andre 18 – ægte forudsigelser

For de 18 galakser inden for 20 % er medianfejlen 6,8 %, hvilket ligger inden for observationsusikkerhederne. Disse spænder overRd fra 1,3 til 5,8 kpc ogVf fra 76 til 278 km/s. BeeTheory forudsiger korrekt denne faktor 3,7 i hastighed – Tully-Fisher-hældningen – med to universelle konstanter.

5. Fysisk betydning – hvad skaleringen afslører

5.1 Den universelle dimensionsløse kobling

Med Kd = K0/Rd og ℓd = cRd er den dimensionsløse BeeTheory-kobling:

Effektiv kobling – skalerer med galaksestørrelsen \(\lambda_{\mathrm{eff}}=K_d\ell_d^2=\frac{K_0}{R_d}(cR_d)^2=K_0c^2R_d\) \(\lambda_{\mathrm{eff}}=0.3759\times41.0\times R_d=15.4R_d\ \text{(kpc)}\)

λeff vokser medRd. Større galakser genererer forholdsmæssigt mere mørk masse. Dette er bi-teoriens forudsigelse af, hvorfor massive spiraler er mere domineret af mørkt stof end dværge.

5.2 Forbindelse til den radiale accelerationsrelation

McGaugh et al. fandt, at den observerede centripetalacceleration gobs = Vc2/R er en universel funktion af det baryoniske bidrag gbar = GMbar/R². I BeeTheory opstår denne relation, fordi:

BeeTheory-forudsigelse af RAR \(g_{\mathrm{obs}}=g_{\mathrm{bar}}+g_{\mathrm{dark}}\) \(g_{\mathrm{dark}}\propto K_0c^2g_{\mathrm{bar}}^{1/2}G^{1/2}\)

Skaleringen gdarkgbar1/2 giver den observerede RAR-form. At udlede den nøjagtige RAR-kurve fra BeeTheory er den umiddelbare næste teoretiske opgave.

Data: Lelli, F., McGaugh, S. S., Schombert, J. M., SPARC: Mass Models for 175 Disk Galaxies with Spitzer Photometry and Accurate Rotation Curves, AJ 152, 157, 2016.

RAR: McGaugh, S. S., Lelli, F., Schombert, J. M., PRL 117, 201101, 2016.

BTFR: Lelli, F. et al, ApJ 816, 2016.

BeeTheory: Dutertre, X., BeeTheory.com v2, 2023, udvidet 2025.

Masse-til-lys: Υ★ = 0,5 M⊙/L⊙ ved 3,6 μm, McGaugh & Schombert 2014.