Bijentheorie – SPARC melkwegstelsels – Aangepaste pasvorm – 2025

BeeTheory toegepast op 20 externe sterrenstelsels:
Aangepaste formule en blinde testmethodologie

De SPARC-catalogus bevat 175 sterrenstelsels met gemeten baryonische profielen en rotatiecurves. We passen de BeeTheory-vergelijking voor donkere massa toe – waarbij we de schaalwet aanpassen aan de melkwegpopulatie – en rapporteren het resultaat: 18 van de 20 melkwegstelsels worden voorspeld binnen 20% van hun waargenomen vlakke rotatiesnelheid, met χ²/dof = 0,93.

0. Resultaten – Eerst vermeld

Beste pasvorm – 20 SPARC melkwegstelsels, Q = 1, hoge kwaliteit

Met de aangepaste BeeTheory-formule Kd = K0/Rd en ℓd = c – Rd passen twee universele constanten bij alle 20 melkwegstelsels tegelijk.

De donkere massadichtheid van elk sterrenstelsel wordt alleen voorspeld op basis van de baryonische schijfparameters – geen afstemming per sterrenstelsel.

Best-fit parameters: K0 = 0,3759, dimensieloos, en c = 6,40, dimensieloos. Resultaat: 18/20 melkwegstelsels binnen 20% van de waargenomen Vf, χ²/dof = 0,93. Twee uitbijters, CamB en NGC 3741, zijn door gas gedomineerde dwergen waar de modellering van de stellaire schijf niet werkt.

18/20

Binnen 20% van Vf

6.8%

Mediaanfout

0.93

χ²/dof

2

Universele constanten

0.91

Pearson r, Tully-Fisher

Aangepaste BeeTheory-formule – aangepast voor melkwegpopulatie \(\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 melkwegstelsels – gegevens en voorspellingen

Het SPARC-monster omvat melkwegstelsels van vijf decennia in helderheid, van dwergonregelmatige tot massieve spiralen. Voor elk sterrenstelsel zijn de invoerparameters rechtstreeks overgenomen uit tabel 1 van Lelli et al. 2016: schijfschaalradius Rd, centrale oppervlaktehelderheid Σd, HI-gasmassaMHI en vlakke rotatiesnelheid Vf.

De stellaire massa wordt berekend als M★ = Υ★ × L3.6, met Υ★ = 0,5 M⊙/L⊙. De gasmassa wordt berekend als Mgas = 1,33 ×MHI.

Melkweg Rd kpc Kd kpc-¹ ℓd kpc Vf obs Vbar Vdark VBT Fout Status
Sterrenstelsel tabel wordt geladen…

Alle snelheden in km/s. Fout = (VBTVf)/Vf. Parameters: Kd = 0,3759/Rd, ℓd = 6,40 × Rd. BeeTheorie voorspelling geëvalueerd bijReval = 5Rd.

Waargenomen VBT vs Vf – 20 SPARC melkwegstelsels, aangepaste BeeTheory
Binnen 20%, 18 sterrenstelsels Uitschieters: CamB, NGC 3741 Perfecte voorspelling, 1:1 ±20% omhullende

2. De gewijzigde formule – Waarom K ∝ 1/Rd

De originele BeeTheory Melkweg fit gebruikte een enkele koppelingsconstante K = 0,02365 kpc-¹ met coherentielengte ℓ = 3,17Rd. Bij blinde toepassing op 20 SPARC melkwegstelsels werd Vf systematisch met ongeveer 50% onderschat.

De analyse per melkwegstelsel liet een duidelijk patroon zien: de vereiste koppelingsconstante varieert als K ∝ 1/Rd.

2.1 Van één constante naar een schaalwet

Het belangrijkste inzicht is dimensionaal. De BeeTheory donkere dichtheid bij straal r van een schijf met schaal Rd en oppervlaktedichtheid Σ0 is, in het asymptotische vlakke-rotatieregime r ≪ ℓ:

Asymptotische donkere dichtheid, 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\)

De vlakke rotatiesnelheid schaalt dan als:

Vlakke rotatiesnelheid van BeeTheory donkere 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}\)

De waargenomen baryonische Tully-Fisherrelatie geeft Vf4Mbar, wat betekent dat VfMbar1/4. Om dit door de BeeTheory te kunnen reproduceren, hebben we Vf2M★/Rd nodig, de gemiddelde oppervlaktedichtheid van de schijf. Dit vereist:

Vereiste schaling om de Tully-Fisher-helling te reproduceren \(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}\)

Deze schaling is geen ad-hoc patch – het is wat de Tully-Fisher relatie vereist. Een koppeling K ∝ 1/Rd betekent dat compactere schijven een sterker donker veld per massa-eenheid opwekken.

2.2 De coherentielengte – Waarom c = 6,40 ≠ 3,17

De Melkweg fit gaf cMW = ℓd/Rd = 3,17. Het SPARC monster geeft cSPARC = 6,40, ongeveer twee keer zo groot. Er zijn twee verklaringen mogelijk:

Beide effecten zijn reëel en meetbaar. Voor de definitieve waarde van c moeten gas- en stellaire schijven afzonderlijk worden gemodelleerd.

Volledige aangepaste BeeTheory-formule – twee universele constanten \(\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. De berekening – Stap voor stap

Voor elk SPARC sterrenstelsel verloopt de BeeTheory voorspelling in vijf stappen. Per melkwegstelsel worden geen vrije parameters aangepast.

1
Baryonische inputs uit SPARC Tabel 1 lezen

Rd, Σd,MHI en Vf. Converteer Σ0 = Σd × Υ★ × 10⁶ M⊙/kpc², en Mgas = 1,33 ×MHI.

2
Bereken BeeTheory-parameters van Rd

Kd = K0/Rd = 0,3759/Rd, ℓd = cRd = 6,40Rd, en αd = 1/ℓd. Geen aanpassing.

3
Integreer de donkere dichtheid op 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’\)

Numerieke integratie met 60 ringen, R′ ∈ [0, 8Rd]. Integreer dan sferisch om ingesloten donkere massa Mdark(<5Rd) te krijgen.

4
Bereken de baryonische cirkelsnelheid
\(V_{\mathrm{bar}}(R)=\sqrt{\frac{G[M_{\mathrm{disk}}(<R)+M_{\mathrm{gas}}(<R)+M_{\mathrm{bulge}}(<R)]}{R}}\)
5
Voorspel de totale cirkelsnelheid
\(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}}\)

Vergelijk met de waargenomen Vf. Fout = (VBTVf)/Vf.

Vdark/Vbar-verhouding – Hoe dominant is de donkere component?

4. Waarom blind testen de enige eerlijke test is

Een model dat de gegevens reproduceert waarop het gekalibreerd is, bewijst niets. Elk model, zelfs een fout model, kan getuned worden om aan de trainingsgegevens te voldoen. De enige wetenschappelijk zinvolle test is een blinde voorspelling: pas het model toe op gegevens die het nog nooit heeft gezien, met parameters die bevroren zijn van de kalibratie, en rapporteer het resultaat – wat het ook is.

4.1 Wat “Blind” hier betekent

De BeeTheory parameters K0 = 0,3759 en c = 6,40 werden bepaald door de 20 SPARC melkwegstelsels tegelijkertijd te passen. Ze liggen nu vast.

De blinde test zou zijn: pas deze parameters toe op de resterende 155 SPARC melkwegstelsels, die niet in de fit zijn gebruikt, en rapporteer het resultaat voordat u naar hun waargenomen rotatiecurves kijkt. Deze test is nog niet uitgevoerd – het is de volgende stap.

De oorspronkelijke Melkwegparameters, Kd = 0,02365 en ℓd = 3,17Rd, werden bepaald op één enkel sterrenstelsel. Ze zonder aanpassing toepassen op SPARC gaf 0/20 melkwegstelsels correct – een eerlijke en belangrijke mislukking. Deze mislukking onthulde de K ∝ 1/Rd schaling.

4.2 Statistische betekenis van aanpaskwaliteit

Met χ²/dof = 0,93 voor 20 melkwegstelsels past het model ongeveer op het verwachte niveau van de aangenomen 15% snelheidsonzekerheden.

Chi-kwadraat interpretatie \(\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{(exclusief CamB)},\quad p=2{(K_0,c),\quad \mathrm{dof}=17\)

Een waarde van 0,93 ligt heel dicht bij de ideale 1,0. Het model houdt rekening met de spreiding op het niveau van de meetonzekerheid.

4.3 De twee uitschieters

CamB – zuivere gasdwerg, Vf = 2,0 km/s

CamB heeft bijna geen stellaire massa, M★ ≈ 2×10⁷ M⊙. De BeeTheory-formule gebruikt Σ0e-R/Rd als bron – maar in CamB zijn de baryonen bijna helemaal HI-gas, geen sterren. Het stellaire schijfmodel is niet van toepassing.

NGC 3741 – 47% overschat

NGC 3741 is een kleine dwerg met een lage oppervlaktehelderheid en een zeer uitgebreide HI-schijf. De BeeTheory-bron, de stellaire schijf, onderschat de werkelijke baryonische omvang. Het opnemen van de gasschijf als een aparte broncomponent met een grotere schaalradius zou de voorspelde donkere massa verminderen en de overschatting corrigeren.

De andere 18 – echte voorspellingen

Voor de 18 sterrenstelsels binnen 20% is de mediane fout 6,8%, ruim binnen de waarnemingsonzekerheden. Deze variëren Rd van 1,3 tot 5,8 kpc en Vf van 76 tot 278 km/s. BeeTheory voorspelt dit snelheidsbereik met een factor 3,7 – de Tully-Fisher-helling – correct met twee universele constanten.

5. Fysische betekenis – Wat de schaling onthult

5.1 De universele dimensieloze koppeling

Met Kd = K0/Rd en ℓd = cRd is de dimensieloze BeeTheory koppeling:

Effectieve koppeling – schaalt met de grootte van het melkwegstelsel \(\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 groeit met Rd. Grotere sterrenstelsels genereren verhoudingsgewijs meer donkere massa. Dit is de voorspelling van de BeeTheory waarom massieve spiralen meer donkere materie bevatten dan dwergen.

5.2 Verband met de radiale versnellingsrelatie

McGaugh et al. vonden dat de waargenomen centripetale versnellinggobs = Vc2/R een universele functie is van de baryonische bijdrage gbar = GMbar/R². In de Bijentheorie ontstaat deze relatie omdat:

BeeTheorie voorspelling van de RAR \(g_{\mathrm{obs}}=g_{\mathrm{bar}}+g_{\mathrm{dark}}\) \(g_{\mathrm{dark}}\propto K_0c^2g_{\mathrm{bar}}^{1/2}G^{1/2}\)

De gdarkgbar1/2 schaling produceert de waargenomen RAR-vorm. Het afleiden van de exacte RAR-curve uit BeeTheory is de eerstvolgende theoretische taak.

Gegevens: 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.

Bijentheorie: Dutertre, X., BeeTheory.com v2, 2023, uitgebreid 2025.

Lichtmassa: Υ★ = 0,5 M⊙/L⊙ bij 3,6 μm, McGaugh & Schombert 2014.