BeeTheory – SPARC Galaxies – Justerad anpassning – 2025

BeeTheory tillämpad på 20 externa galaxer:
Justerad formel och blind testmetodik

SPARC-katalogen innehåller 175 galaxer med uppmätta baryoniska profiler och rotationskurvor. Vi tillämpar BeeTheory-ekvationen för mörk massa – justerar dess skalningslag för att matcha galaxpopulationen – och rapporterar resultatet: 18 av 20 galaxer förutspås inom 20% av deras observerade platta rotationshastighet, med χ²/dof = 0,93.

0. Resultat – först uttalat

Bästa anpassning – 20 SPARC-galaxer, Q = 1, hög kvalitet

Med den modifierade BeeTheory-formeln Kd = K0/Rd och ℓd = c –Rd, två universella konstanter som passar alla 20 galaxer samtidigt.

Den mörka massans densitet i varje galax förutsägs enbart utifrån dess baryoniska diskparametrar – ingen inställning per galax.

Parametrar med bästa passform: K0 = 0,3759, dimensionslös, och c = 6,40, dimensionslös. Resultat: 18/20 galaxer inom 20% av observeradVf, χ²/dof = 0,93. Två avvikande galaxer, CamB och NGC 3741, är gasdominerade dvärgar där modelleringen av stjärnskivan inte fungerar.

18/20

Inom 20 % avVf

6.8%

Medianfel

0.93

χ²/dof

2

Universella konstanter

0.91

Pearson r, Tully-Fisher

Modifierad BeeTheory-formel – justerad för galaxpopulation \(\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-galaxerna – data och förutsägelser

SPARC-urvalet omfattar galaxer som spänner över fem decennier i luminositet, från irreguljära dvärggalaxer till massiva spiraler. För varje galax hämtas ingångsparametrarna direkt från Lelli et al. 2016 Tabell 1: diskens skalradieRd, den centrala ytans ljusstyrka Σd, HI-gasmassanMHI och den plana rotationshastighetenVf.

Stellarmassan beräknas som M★ = Υ★ × L3,6, med Υ★ = 0,5 M⊙/L⊙. Gasmassan är beräknad som Mgas = 1,33 ×MHI.

Galax Rd kpc Kd kpc-¹ ℓd kpc Vf obs Vbar Vdark VBT fel Status
Laddar galaxtabell…

Alla hastigheter i km/s. Fel = (VBTVf)/Vf. Parametrar: Kd = 0,3759/Rd, ℓd = 6,40 ×Rd. BeeTheory-förutsägelse utvärderad vidReval = 5Rd.

VBT vsVf observerad – 20 SPARC-galaxer, justerad BeeTheory
Inom 20%, 18 galaxer Avvikande värden: CamB, NGC 3741 Perfekt förutsägelse, 1:1 ±20% inom skalet

2. Den modifierade formeln – varför K ∝ 1/Rd

Den ursprungliga BeeTheory-anpassningen av Vintergatan använde en enda kopplingskonstant K = 0,02365 kpc-¹ med koherenslängden ℓ = 3,17Rd. När detta tillämpades blint på 20 SPARC-galaxer underskattadesVf systematiskt med cirka 50%.

Analysen per galax visade ett tydligt mönster: den erforderliga kopplingskonstanten varierar som K ∝ 1/Rd.

2.1 Från en konstant till en skalningslag

Den viktigaste insikten är dimensionell. Den mörka tätheten enligt BeeTheory vid radien r från en skiva med skalan Rd och yttätheten Σ0 är, i den asymptotiska regimen med platt rotation r ≪ ℓ:

Asymptotisk mörk densitet,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 plana rotationshastigheten skalar då som:

Platt rotationshastighet från BeeTheory mörk massa \(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}\)

Det observerade baryoniska Tully-Fisher-förhållandet anger Vf4Mbar, vilket betyderVfMbar1/4. För att detta ska kunna återges av BeeTheory behöver vi Vf2M★/Rd, den genomsnittliga ytdensiteten för skivan. Detta kräver:

Nödvändig skalning för att återge Tully-Fisher-sluttningen \(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}\)

Denna skalning är inte en ad hoc-lapp – det är vad Tully-Fisher-relationen kräver. En koppling K ∝ 1/Rd innebär att mer kompakta skivor genererar ett starkare mörkt fält per massaenhet.

2.2 Koherenslängden – varför c = 6,40 ≠ 3,17

Vintergatans passning gav cMW = ℓd/Rd = 3,17. SPARC-urvalet ger cSPARC = 6,40, ungefär dubbelt så stort. Två förklaringar är möjliga:

Båda effekterna är verkliga och mätbara. För att fastställa det definitiva värdet på c krävs att gas- och stjärnskivorna modelleras separat.

Komplett modifierad BeeTheory-formel – två universalkonstanter \(\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. Beräkningen – steg för steg

För varje SPARC-galax sker BeeTheory-förutsägelsen i fem steg. Inga fria parametrar justeras per galax.

1
Läs baryoniska indata från SPARC Tabell 1

Rd, Σd,MHI ochVf. Konvertera Σ0 = Σd × Υ★ × 10⁶ M⊙/kpc², och Mgas = 1,33 ×MHI.

2
Beräkna BeeTheory-parametrar frånRd

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

3
Integrera den mörka densiteten vid 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 ringar, R′ ∈ [0, 8Rd]. Integrera sedan sfäriskt för att få den inneslutna mörka massan Mdark(<5Rd).

4
Beräkna baryonisk cirkulär hastighet
\(V_{\mathrm{bar}}(R)=\sqrt{\frac{G[M_{\mathrm{disk}}(<R)+M_{\mathrm{gas}}(<R)+M_{\mathrm{bulge}}(<R)]}{R}}\)
5
Förutsäg total cirkulär hastighet
\(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}}\)

Jämför med observeradVf. Fel = (VBTVf)/Vf.

Förhållandet Vdark/Vbar – Hur dominerande är den mörka komponenten?

4. Varför blind testning är det enda ärliga testet

En modell som reproducerar de data som den kalibrerades på bevisar ingenting. Varje modell, även en felaktig sådan, kan justeras för att passa dess träningsdata. Det enda vetenskapligt meningsfulla testet är en blind förutsägelse: tillämpa modellen på data som den aldrig har sett, med parametrar som är frusna från kalibreringen, och rapportera resultatet – oavsett vad det är.

4.1 Vad ”blind” betyder här

BeeTheory-parametrarna K0 = 0,3759 och c = 6,40 bestämdes genom att anpassa de 20 SPARC-galaxerna samtidigt. De är nu fixerade.

Blindtestet skulle vara: att tillämpa dessa parametrar på de återstående 155 SPARC-galaxerna, som inte användes i anpassningen, och rapportera resultatet innan man tittar på deras observerade rotationskurvor. Detta test har ännu inte utförts – det är nästa steg.

De ursprungliga Vintergatsparametrarna, Kd = 0,02365 och ℓd = 3,17Rd, bestämdes på en enda galax. Om man applicerade dem på SPARC utan justering fick man 0/20 galaxer korrekta – ett ärligt och viktigt misslyckande. Detta misslyckande avslöjade skalningen K ∝ 1/Rd.

4.2 Statistisk innebörd av anpassningskvalitet

Med χ²/dof = 0,93 för 20 galaxer passar modellen in på ungefär den förväntade nivån för de antagna hastighetsosäkerheterna på 15 %.

Chi-squared tolkning \(\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{(exklusive CamB)},\qquad p=2\ (K_0,c),\qquad \mathrm{dof}=17\)

Ett värde på 0,93 är mycket nära idealet 1,0. Modellen tar hänsyn till spridningen på nivån för mätosäkerhet.

4.3 De två extremvärdena

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

CamB har nästan ingen stjärnmassa, M★ ≈ 2×10⁷ M⊙. I BeeTheory-formeln används Σ0e-R/Rd som källa – men i CamB består baryonerna nästan uteslutande av HI-gas, inte av stjärnor. Stjärndiskmodellen är inte tillämpbar.

NGC 3741 – överskattad med 47%.

NGC 3741 är en liten dvärg med låg ytljusstyrka och en mycket utsträckt HI-skiva. BeeTheory-källan, stjärndisken, underskattar den faktiska baryoniska omfattningen. Om gasskivan inkluderas som en separat källkomponent med större skalradie skulle den beräknade mörka massan minska och överskattningen korrigeras.

De övriga 18 – äkta förutsägelser

För de 18 galaxerna inom 20% är medianfelet 6,8%, väl inom observationsosäkerheterna. Dessa spänner överRd från 1,3 till 5,8 kpc ochVf från 76 till 278 km/s. BeeTheory förutsäger korrekt denna faktor 3,7 i hastighetsintervall – Tully-Fisher-lutningen – med två universella konstanter.

5. Fysikalisk innebörd – vad skalningen avslöjar

5.1 Den universella dimensionslösa kopplingen

Med Kd = K0/Rd och ℓd = cRd är den dimensionslösa BeeTheory-kopplingen:

Effektiv koppling – skalar med galaxstorlek \(\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 växer medRd. Större galaxer genererar proportionellt mer mörk massa. Detta är BeeTheory-förutsägelsen för varför massiva spiraler är mer dominerade av mörk materia än dvärgar.

5.2 Koppling till det radiella accelerationsrelationen

McGaugh et al. fann att den observerade centripetalaccelerationen gobs = Vc2/R är en universell funktion av det baryoniska bidraget gbar = GMbar/R². I BeeTheory uppstår detta förhållande därför att:

BeeTheory-förutsägelse av RAR \(g_{\mathrm{obs}}=g_{\mathrm{bar}}+g_{\mathrm{dark}}\) \(g_{\mathrm{dark}}\propto K_0c^2g_{\mathrm{bar}}^{1/2}G^{1/2}\)

Skalningen gdarkgbar1/2 ger den observerade RAR-formen. Att härleda den exakta RAR-kurvan från BeeTheory är den omedelbara nästa teoretiska uppgiften.

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, utökad 2025.

Massa-till-ljus: Υ★ = 0,5 M⊙/L⊙ vid 3,6 μm, McGaugh & Schombert 2014.