BeeTheory – Grondslagen – Technische noot X

Anatomie van de residuen:
Een lineaire trend met schijfgrootte

De blinde test met 94 sterrenstelsels van noot IX liet een systematische residuele trend met de schijfgrootte zien. Deze notitie karakteriseert die trend kwantitatief, isoleert de grootste afwijkingen aan elke kant en identificeert de structurele oorsprong van de spreiding.

1. Het resultaat eerst

Een lineair residu, twee tegengestelde populaties

De voorspellingsfout schaalt lineair met de schaallengte van de schijf: $\text{error},(\%) \approx -31.7 + 12.8$,R_d$, met Pearson correlatie $r = +0.75$. De lijn kruist nul bij $R_d = 2.48$ kpc, in wezen de schijfgrootte van de Melkweg die de kalibratie verankerde. De twee uitersten van deze regressie komen overeen met twee fysiek verschillende populaties van uitschieters: grote massieve spiralen (te hoge voorspelling) aan de ene kant, compacte dwergen (te lage voorspelling) aan de andere kant.

2. Het residu is lineair in $R_d$

Door de voorspellingsfout uit te zetten tegen $R_d$, met elk punt gekleurd door het Hubble-type, wordt de lineariteit van de trend onmiddellijk zichtbaar. De rode lijn is de lineaire regressie van de fout op $R_d$ over alle 94 blinde sterrenstelsels.

Voorspellingsfout t.o.v. schijfgrootte – lineair in $R_d$, gekleurd door Hubble-type 0.30.5123510 -50%-25%+0%+25%+50%+75% → Over-voorspeld gebied ±30% band ← Onder-voorspeld gebied Schijfschaallengte $R_d$ (kpc) – logoschaal Voorspellingsfout (%) F579-V1 (Sd-Im) Rd=3.20 Vf=105 err=-12%F583-1 (Sd-Im) Rd=1,80 Vf=83 err=-30%F583-4 (Sd-Im) Rd=1,40 Vf=67 err=-29%IC2574 (Sd-Im) Rd=2.80 Vf=69 err=+35%KK98-251 (Sd-Im) Rd=0,30 Vf=17 err=-23%M33 (Sc-Scd) Rd=1,40 Vf=100 err=-2%NGC0055 (Sd-Im) Rd=1,80 Vf=87 err=-1%NGC0100 (Sc-Scd) Rd=2.30 Vf=83 err=+13%NGC0247 (Sd-Im) Rd=2,40 Vf=90 err=+20%NGC0289 (Sb-Sbc) Rd=3,50 Vf=155 err=+29%NGC0300 (Sd-Im) Rd=1,50 Vf=76 err=+0%NGC0801 (Sc-Scd) Rd=5,80 Vf=208 err=+57%NGC0891 (Sb-Sbc) Rd=4.10 Vf=212 err=+7%NGC0925 (Sd-Im) Rd=3,10 Vf=105 err=+48%NGC1003 (Sc-Scd) Rd=2.80 Vf=115 err=+12%NGC1090 (Sb-Sbc) Rd=3,80 Vf=170 err=+17%NGC1705 (S0-Sa) Rd=0,60 Vf=54 err=-19%NGC2366 (Sd-Im) Rd=1,30 Vf=55 err=+14%NGC2403 (Sc-Scd) Rd=1,80 Vf=131 err=-4%NGC2683 (Sb-Sbc) Rd=2,90 Vf=175 err=+15%NGC2903 (Sb-Sbc) Rd=2,60 Vf=184 err=-0%NGC2915 (Sd-Im) Rd=0,50 Vf=85 err=-38%NGC2955 (Sb-Sbc) Rd=5,50 Vf=266 err=+53%NGC2976 (Sc-Scd) Rd=0,75 Vf=80 err=-37%NGC3109 (Sd-Im) Rd=1,40 Vf=68 err=-25%NGC3521 (Sb-Sbc) Rd=2,80 Vf=225 err=+5%NGC3621 (Sd-Im) Rd=2.10 Vf=149 err=+28%NGC3726 (Sc-Scd) Rd=3.00 Vf=152 err=+19%NGC3741 (Sd-Im) Rd=0,68 Vf=51 err=+26%NGC3769 (Sc-Scd) Rd=2.80 Vf=112 err=+25%NGC3877 (Sc-Scd) Rd=2.70 Vf=163 err=+12%NGC3893 (Sc-Scd) Rd=2.80 Vf=159 err=+14%NGC3949 (Sb-Sbc) Rd=1,40 Vf=125 err=-21%NGC3953 (Sb-Sbc) Rd=3,50 Vf=200 err=-6%NGC3972 (Sc-Scd) Rd=1,60 Vf=135 err=-27%NGC3992 (Sb-Sbc) Rd=3.80 Vf=242 err=-15%NGC4010 (Sc-Scd) Rd=1,80 Vf=128 err=-14%NGC4013 (Sc-Scd) Rd=2.20 Vf=185 err=+8%NGC4051 (Sb-Sbc) Rd=1,90 Vf=110 err=+3%NGC4085 (Sc-Scd) Rd=1,20 Vf=135 err=-41%NGC4088 (Sb-Sbc) Rd=1,90 Vf=175 err=-27%NGC4100 (Sb-Sbc) Rd=1,80 Vf=162 err=-28%NGC4138 (S0-Sa) Rd=1,30 Vf=150 err=-44%NGC4157 (Sb-Sbc) Rd=2,60 Vf=185 err=-1%NGC4183 (Sc-Scd) Rd=1,60 Vf=110 err=-36%NGC4214 (Sd-Im) Rd=0,50 Vf=68 err=-27%NGC4217 (Sb-Sbc) Rd=2,80 Vf=180 err=+5%NGC4389 (Sb-Sbc) Rd=1,20 Vf=110 err=-43%NGC4559 (Sc-Scd) Rd=3.20 Vf=123 err=+28%NGC5005 (Sb-Sbc) Rd=3.00 Vf=260 err=-8%NGC5033 (Sc-Scd) Rd=4.50 Vf=195 err=+44%NGC5055 (Sb-Sbc) Rd=3,50 Vf=180 err=+32%NGC5371 (Sb-Sbc) Rd=3.80 Vf=225 err=+14%NGC5585 (Sd-Im) Rd=1,50 Vf=87 err=-5%NGC5907 (Sc-Scd) Rd=4,20 Vf=210 err=+32%NGC5985 (Sb-Sbc) Rd=4,50 Vf=295 err=+2%NGC6015 (Sc-Scd) Rd=2,40 Vf=142 err=+6%NGC6195 (Sb-Sbc) Rd=5,20 Vf=260 err=+46%NGC6503 (Sc-Scd) Rd=2.40 Vf=121 err=+39%NGC6674 (Sb-Sbc) Rd=5,50 Vf=260 err=+46%NGC6789 (Sd-Im) Rd=0,30 Vf=60 err=-63%NGC6946 (Sc-Scd) Rd=2,60 Vf=180 err=+10%NGC7331 (Sb-Sbc) Rd=3,20 Vf=265 err=+4%NGC7793 (Sd-Im) Rd=1,80 Vf=118 err=-1%UGC00128 (Sd-Im) Rd=7,50 Vf=135 err=+80%UGC02259 (Sd-Im) Rd=1,60 Vf=90 err=+0%UGC02487 (S0-Sa) Rd=7,50 Vf=330 err=+41%UGC02885 (Sc-Scd) Rd=8.50 Vf=290 err=+52%UGC05716 (Sd-Im) Rd=2,00 Vf=75 err=-6%UGC05721 (Sd-Im) Rd=1,20 Vf=85 err=+0%UGC05750 (Sd-Im) Rd=4,50 Vf=80 err=+38%UGC05764 (Sd-Im) Rd=0,40 Vf=57 err=-46%UGC05829 (Sd-Im) Rd=1,60 Vf=69 err=-10%UGC06399 (Sd-Im) Rd=2,50 Vf=89 err=+11%UGC06446 (Sd-Im) Rd=1,80 Vf=87 err=+6%UGC06614 (S0-Sa) Rd=4,50 Vf=200 err=+19%UGC06628 (Sd-Im) Rd=2.50 Vf=75 err=+7%UGC06667 (Sd-Im) Rd=2.50 Vf=90 err=+16%UGC06917 (Sd-Im) Rd=2.50 Vf=110 err=-12%UGC06983 (Sc-Scd) Rd=2.50 Vf=113 err=+4%UGC07125 (Sd-Im) Rd=4.50 Vf=75 err=+37%UGC07151 (Sc-Scd) Rd=1,30 Vf=82 err=-32%UGC07261 (Sd-Im) Rd=1,10 Vf=72 err=-29%UGC07399 (Sd-Im) Rd=1,40 Vf=93 err=-19%UGC07690 (Sd-Im) Rd=0,70 Vf=62 err=-29%UGC08286 (Sc-Scd) Rd=1,30 Vf=84 err=-4%UGC08490 (Sd-Im) Rd=0,65 Vf=80 err=-29%UGC08550 (Sd-Im) Rd=1,50 Vf=67 err=-17%UGC09037 (Sc-Scd) Rd=3.50 Vf=160 err=-8%UGC11455 (Sc-Scd) Rd=5.50 Vf=275 err=-28%UGC11557 (Sd-Im) Rd=3.00 Vf=90 err=+1%UGC11820 (Sd-Im) Rd=4,50 Vf=90 err=+26%UGCA281 (Sd-Im) Rd=0,50 Vf=40 err=-36%UGCA442 (Sd-Im) Rd=1,00 Vf=57 err=-44% err ≈ -31,7 +12,8-Rd nul bij Rd=2,48 kpc Pearson correlatier = +0.749 S0-Sa (T=0-2)Sb-Sbc (T=3-4)Sc-Scd (T=5-6)Sd-Im (T=7-10)
94 blinde sterrenstelsels uitgezet tegen schijfgrootte, gekleurd door Hubble-type. De rode lijn is de lineaire regressie van de fout op $R_d$. Deze kruist nul bij $R_d = 2,48$ kpc – in wezen de schijfgrootte die de oorspronkelijke kalibratie verankerde.

Fout als functie van schijfgrootte

$${fout},(\%) \; -31.7 \; +; 12.8 \times R_d \,[\text{kpc}]$$.

Lineaire fit op 94 blinde melkwegstelsels, Pearson $r = +0.75$, RMSE van residuen $= 18.4%$.

Vergelijking van functionele vormen

Er werden verschillende alternatieve parametrisaties vergeleken. De lineaire vorm is statistisch niet te onderscheiden van log- en vierkantswortelalternatieven:

ModelPeer $r$RMSEOpmerking
$text{err} = a + b (lineair)$+0.749$$18.4\%$Schoonste analytische vorm
$tekst{err} = a + b,\log_{10}R_d$$+0.748$$18.4\%$Statistisch gelijkwaardig
$tekst{err} = a + b,qrt{R_d}$$+0.768$$17.7\%$Marginaal beter, geen echte winst
$tekst{err} = a + b, R_d + c, R_d^2$$17.8\%$Kwadratische term zeer klein ($c ≤ -1.1$)

De lineaire vorm wordt daarom aangenomen als de eenvoudigste getrouwe beschrijving van de gegevens.

Hubble-type verdeling langs de lijn

Hubble-klasse $N$ Mediaan $R_d$ (kpc) Mediaanfout Positie
S0-Sa (vroeg type)42.9$+0.0\%$Midden, bij de nuldoorgang
Sb-Sbc (intermediair)233.2$+3.9\%$Rechts van het midden; staart in het overvoorspelde gebied
Sc-Scd (late spiraal)272.5$+7.7\%$Verspreid over het diagram
Sd-Im (dwerg / onregelmatig)401.6$-3.2\%$Links; staart in het ondergevoorspelde gebied

Het kleurenpatroon in de figuur is niet onafhankelijk van de lineaire trend – het is hetzelfde patroon dat door de morfologie-as heen te zien is. De Hubble-reeks in schijfstelsels correleert met de schijfgrootte: laat-type dwergen zijn overwegend compact, tussenliggende spiralen zijn overwegend groot. Elke kleur ligt daarom op een ander stuk van de regressielijn, met Sd-Im aan de linkerkant, Sc-Scd in het midden en Sb-Sbc aan de rechterkant.

Een structureel residu, geen willekeurige ruis

Een spreiding die lineair afhangt van één fysische parameter, en nul kruist op het kalibratiepunt, is de signatuur van een ontbrekende additieve constante in één van de relaties van het model, niet van willekeurige waarnemingsspreiding. De afwijking kan gecorrigeerd worden: het kan geabsorbeerd worden door een enkele extra vrijheidsgraad in de coherentie-lengtewet.

3. De tien meest overvoorspelde sterrenstelsels

Dit zijn de sterrenstelsels waarvoor BeeTheory een hogere vlakke rotatiesnelheid voorspelt dan waargenomen. Gesorteerd op de grootte van het residu:

GalaxyHubble-type$R_d$ (kpc)$M_\star/10^{10}$$f_text{gas}$$Sigma_d$$V_f$$V_text{tot}$Fout
UGC00128Sd-Im7.501.060.3960135243+80.0%
NGC0801Sb-Sbc5.802.010.32190208326+56.6%
NGC2955Sb-Sbc5.503.990.23420266406+52.7%
UGC02885Sc-Scd8.503.400.41150290441+52.0%
NGC0925Sc-Scd3.100.220.7572105155+48.0%
NGC6195Sb-Sbc5.203.400.26400260380+46.3%
NGC6674Sb-Sbc5.503.330.29350260380+46.2%
NGC5033Sb-Sbc4.501.270.46200195280+43.7%
UGC02487S0-Sa7.505.300.23300330465+40.8%
NGC6503Sc-Scd2.400.380.55210121168+38.9%
EigendomMediaanwaardeBereikVergelijking
$R_d$4,5 kpc2.4 – 8.5$2 keer$ groter dan mediaan
$M_star$$1,3 maal 10^{10},M_odot$.$2,2 \times 10^{9}$ – $5,3 \times 10^{10}$$8 keer$ massiever
$f_text{gas}$$0.41$$0.23$ – $0.87$Onder de mediaan (0,64)
Hubble $T$$5$ (Sbc)$1$ – $8$Geconcentreerd in tussenliggende spiralen
$V_f$$195$ km/s$69$ – $330$Snelste rotators in de steekproef

Profiel van de overvoorspelde groep

Grote, massieve, middentype spiralen. Deze sterrenstelsels bevinden zich aan de rechterkant van de regressielijn, ruim boven de nuldoorgang. De coherentie-lengtewet van het model $\ell = c_text{disk},R_d$ produceert waarden van $\ell$ boven 20 kpc in dit regime, en genereert meer golfveldmassa dan de waargenomen rotatie vereist.

4. De tien meest onder-voorspelde sterrenstelsels

Dit zijn de sterrenstelsels waarvoor BeeTheory een lagere vlakke rotatiesnelheid voorspelt dan waargenomen. Gesorteerd op de grootte van het residu:

GalaxyHubble-type$R_d$ (kpc)$M_\star/10^{10}$$f_text{gas}$$Sigma_d$$V_f$$V_text{tot}$Fout
NGC6789Sd-Im0.300.010.532506022-63.0%
UGC05764Sd-Im0.400.000.86805731-45.6%
UGCA442Sd-Im1.000.000.85155732-44.2%
NGC4138S0-Sa1.300.130.3325015085-43.6%
NGC4389Sb-Sbc1.200.070.3715011062-43.4%
NGC4085Sb-Sbc1.200.090.4220013579-41.1%
NGC2915Sd-Im0.500.010.841608553-38.2%
NGC2976Sb-Sbc0.750.040.292208050-37.4%
NGC4183Sc-Scd1.600.030.814011070-36.3%
UGCA281Sd-Im0.500.010.63804026-36.1%
EigendomMediaanwaardeBereikVergelijking
$R_d$1,1 kpc0.30 – 1.80$2 keer$ kleiner dan mediaan
$M_star$$2,7 maal 10^{8},M_odot$.$4 \times 10^{7}$ – $1,3 \times 10^{9}$$6 keer$ minder massief
$f_text{gas}$$0.58$$0.29$ – $0.86$Onder de mediaan (0,64)
Hubble $T$$8$ (Sd)$1$ – $10$Geconcentreerd in laat-type dwergen
$V_f$$82$ km/s$40$ – $150$Trage rotators

Profiel van de onder-voorspelde groep

Compacte, lichte dwergen en kleine spiralen. Deze sterrenstelsels bevinden zich aan de linkerkant van de regressielijn, ruim onder de nuldoorgang. De coherentie-lengte wet $\ell = c_\text{disk},R_d$ produceert $\ell$ van orde $1$-$3$ kpc in dit regime, mogelijk te kort om de volledige omvang van het golfveld te verzamelen.

5. Zij-aan-zij vergelijking van de drie groepen

Eigendom (mediaan) Te hoog voorspeld
(err > +30%, $N = 15$)
Goed voorspeld
(|err| ≤ 30%, $N = 67$)
Onder-voorspeld
(err < -30%, $N = 12$)
$R_d$ (kpc)4.52.41.1
$M_\star / 10^{10}$1.270.150.027
$M_text{gas} / 10^{10}$0.930.270.04
$f_text{gas}$0.410.640.58
$Sigma_d$200140115
Hubble $T$5 (Sbc)6 (Sc)8 (Sd)
$V_f$ (km/s)19511382

Elke eigenschap varieert monotoon van links naar rechts. De over-voorspelde groep is groter, massiever, meer ster-gedomineerd en sneller roterend; de onder-voorspelde groep is kleiner, lichter, gasrijker en langzamer; de goed-voorspelde meerderheid zit er tussenin. De Melkweg ($R_d = 2,6$ kpc, $V_f ca. 230$ km/s) valt van nature binnen het goed voorspelde regime waarin de kalibratie verankerd was.

6. Interpretatie

Het model heeft een enkele koppelingsparameter $lambda$ en drie universele geometrische constanten $(c_text{disk}, c_text{sph}, c_text{arm})$. Deze werden bepaald op een melkwegstelsel van gemiddelde grootte (de Melkweg, $R_d = 2,6$ kpc) en gevalideerd op tweeëntwintig melkwegstelsels van vergelijkbare grootte. De blinde test van Noot IX laat zien dat ze redelijk goed generaliseren, maar met een residu dat lineair met de schijfgrootte verschuift.

Een affiene correctie is voldoende

De lineariteit van het residu in $R_d$ – die goed past bij een enkele rechte lijn die nul kruist bij $R_d = 2,48$ kpc – is het kenmerk van een ontbrekende additieve offset in de coherentie-lengte relatie. De huidige wet van $\ell = c_\text{disk}\,R_d$ koppelt de golf-coherentielengte strikt evenredig aan de schijfschaal. Vervanging door een affiene relatie $\ell = c_text{disk}(R_d – R_0)$, waarbij $R_0$ een kleine offset is van ongeveer $2.5$ kpc, zou een rest opleveren die verdwijnt bij het ijkpunt en lineair groeit aan beide kanten – precies het waargenomen patroon.

De goed voorspelde meerderheid is in grote lijnen representatief

Twee derde van de steekproef valt in de goed voorspelde band. Deze 67 sterrenstelsels bestrijken het volledige Hubble-typenbereik en een factor 100$ in stellaire massa. Het geldigheidsgebied van het model is niet smal: het omvat het grootste deel van de SPARC-populatie, met afwijkingen die geconcentreerd zijn aan de twee uitersten van de schijfgrootte, precies zoals een lineair $R_d$-afhankelijk residu zou opleveren.

7. Samenvatting

1. De voorspellingsfout van de blinde test met 94 sterrenstelsels volgt een duidelijke lineaire trend met de lengte van de schijf: $text{error}(\%) \approx -31.7 + 12.8\,R_d$, met Pearson $r = +0.75$ en RMSE van residuen $= 18.4\%$.

2. De lineaire regressie kruist nul bij $R_d = 2,48$ kpc, in wezen de schijfgrootte van de Melkweg die de kalibratie verankerde. De twee uiteinden van de lijn komen overeen met twee fysiek verschillende populaties van uitschieters.

3. De 15 sterrenstelsels die meer dan $+30%$ te hoog voorspeld zijn, zijn grote, massieve spiralen van het tussentype: mediaan $R_d = 4,5$ kpc, $M_star approx 10^{10},M_odot$, $V_f approx 200$ km/s.

4. De 12 sterrenstelsels die meer dan $-30%$ te laag voorspeld zijn, zijn compacte dwergen met een lage massa: mediaan $R_d = 1,1$ kpc, $M_star approx 3 \times 10^{8},M_odot$, $V_f approx 80$ km/s.

5. De afwijking kan worden geabsorbeerd door een affiene correctie op de coherentie-lengte wet, $ell = c_text{disk}(R_d – R_0)$, met $R_0 approx 2.5$ kpc – waarbij een enkele nieuwe constante wordt geïntroduceerd.


Referenties. 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). – de Vaucouleurs, G. et al. – Derde referentiecatalogus van heldere sterrenstelsels, Springer (1991). – McGaugh, S. S. – De derde wet van galactische rotatie, Galaxies 2, 601 (2014). – Dutertre, X. – Bijentheorie™: Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).

BeeTheory.com – Op golven gebaseerde kwantumzwaartekracht – SPARC-residuen – © Technoplane S.A.S. 2026