BeeTheory – Foundations – Uwaga techniczna X

Anatomia reszt:
Liniowy trend z rozmiarem dysku

Ślepy test 94 galaktyk opisany w Nocie IX wykazał systematyczny trend resztkowy wraz z rozmiarem dysku. Niniejsza notatka charakteryzuje ten trend ilościowo, izoluje największe odchylenia po każdej stronie i identyfikuje strukturalne pochodzenie dyspersji.

1. Najpierw wynik

Liniowa reszta, dwie przeciwne populacje

Błąd predykcji skaluje się liniowo z długością skali dysku: $\text{error}\,(\%) \ około -31.7 + 12.8\,R_d$, z korelacją Pearsona $r = +0.75$. Linia przecina zero przy $R_d = 2,48$ kpc, zasadniczo rozmiarze dysku Drogi Mlecznej, który zakotwiczył kalibrację. Dwa krańce tej regresji odpowiadają dwóm fizycznie odrębnym populacjom odstającym: duże masywne spirale (nadmiernie przewidywane) na jednym końcu, zwarte karły (niedostatecznie przewidywane) na drugim.

2. Reszta jest liniowa w $R_d$

Wykres błędu predykcji względem $R_d$, z każdym punktem pokolorowanym według typu Hubble’a, natychmiast uwidacznia liniowość trendu. Czerwona linia to regresja liniowa błędu względem $R_d$ dla wszystkich 94 ślepych galaktyk.

Wykres 94 ślepych galaktyk w zależności od rozmiaru dysku, pokolorowanych według typu Hubble’a. Czerwona linia to regresja liniowa błędu na $R_d$. Przecina ona zero przy $R_d = 2,48$ kpc – zasadniczo jest to rozmiar dysku, który zakotwiczył oryginalną kalibrację.

Błąd jako funkcja rozmiaru dysku

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

Dopasowanie liniowe dla 94 ślepych galaktyk, Pearson $r = +0.75$, RMSE reszt $= 18.4\%$.

Porównanie form funkcjonalnych

Porównano kilka alternatywnych parametryzacji. Forma liniowa jest statystycznie nieodróżnialna od alternatyw logarytmicznych i kwadratowo-krotnych:

Model	Pearson $r$	RMSE	Komentarz
$\text{err} = a + b\,R_d$ (liniowy)	$+0.749$	$18.4\%$	Najczystsza forma analityczna
$\text{err} = a + b\,\log_{10}R_d$	$+0.748$	$18.4\%$	Statystycznie równoważny
$\text{err} = a + b\,\sqrt{R_d}$.	$+0.768$	$17.7\%$	Nieznacznie lepiej, bez realnych korzyści
$\text{err} = a + b\,R_d + c\,R_d^2$.	–	$17.8\%$	Wyrażenie kwadratowe jest bardzo małe ($c około -1,1$)

W związku z tym przyjęto formę liniową jako najprostszy wierny opis danych.

Rozkład typu Hubble’a wzdłuż linii

Klasa Hubble’a	$N$	Mediana $R_d$ (kpc)	Błąd mediany	Stanowisko
S0-Sa (wczesny typ)	4	2.9	$+0.0\%$	Środek, w pobliżu przejścia przez zero
Sb-Sbc (pośredni)	23	3.2	$+3.9\%$	Na prawo od środka; ogon w obszarze nadmiernie przewidywanym
Sc-Scd (późna spirala)	27	2.5	$+7.7\%$	Rozłożone na wykresie
Sd-Im (karzeł / nieregularny)	40	1.6	$-3.2\%$	Lewa strona; ogon w obszarze niedostatecznie przewidywanym

Wzór kolorów na rysunku nie jest niezależną sygnaturą od trendu liniowego – jest to ta sama sygnatura widziana przez oś morfologii. Sekwencja Hubble’a w galaktykach dyskowych koreluje z rozmiarem dysku: późne karły są przeważnie zwarte, pośrednie spirale są przeważnie duże. Każdy kolor znajduje się zatem wzdłuż innego odcinka linii regresji, z Sd-Im po lewej, Sc-Scd w centrum i Sb-Sbc po prawej.

Pozostałość strukturalna, a nie przypadkowy szum

Rozrzut, który zależy liniowo od pojedynczego parametru fizycznego i przekracza zero w punkcie kalibracji, jest sygnaturą brakującej stałej addytywnej w jednej z relacji modelu, a nie losowego rozrzutu obserwacyjnego. Odchylenie to można skorygować: może ono zostać wchłonięte przez pojedynczy dodatkowy stopień swobody w prawie długości koherencji.

3. Dziesięć galaktyk o najbardziej zawyżonych przewidywaniach

Są to galaktyki, dla których BeeTheory przewiduje płaską prędkość rotacji wyższą niż obserwowana. Posortowane według wielkości pozostałości:

Galaktyka	Typ Hubble’a	$R_d$ (kpc)	$M_\star/10^{10}$.	$f_\text{gas}$.	$\Sigma_d$	$V_f$	$V_\text{tot}$.	Błąd
UGC00128	Sd-Im	7.50	1.06	0.39	60	135	243	+80.0%
NGC0801	Sb-Sbc	5.80	2.01	0.32	190	208	326	+56.6%
NGC2955	Sb-Sbc	5.50	3.99	0.23	420	266	406	+52.7%
UGC02885	Sc-Scd	8.50	3.40	0.41	150	290	441	+52.0%
NGC0925	Sc-Scd	3.10	0.22	0.75	72	105	155	+48.0%
NGC6195	Sb-Sbc	5.20	3.40	0.26	400	260	380	+46.3%
NGC6674	Sb-Sbc	5.50	3.33	0.29	350	260	380	+46.2%
NGC5033	Sb-Sbc	4.50	1.27	0.46	200	195	280	+43.7%
UGC02487	S0-Sa	7.50	5.30	0.23	300	330	465	+40.8%
NGC6503	Sc-Scd	2.40	0.38	0.55	210	121	168	+38.9%

Własność	Wartość mediany	Zasięg	Porównanie
$R_d$	4,5 kpc	2.4 – 8.5	$2\razy$ większa niż mediana
$M_\star$	$1,3 razy 10^{10}\,M_\odot$.	2,2 razy 10^{9}$ – 5,3 razy 10^{10}$.	8 razy bardziej masywny
$f_\text{gas}$.	$0.41$	$0.23$ – $0.87$	Poniżej mediany (0,64)
Hubble $T$	5$ (Sbc)	$1$ – $8$	Skoncentrowany w spiralach pośrednich
$V_f$	$195$ km/s	$69$ – $330$	Najszybsze rotatory w próbce

Profil nadmiernie przewidywanej grupy

Duże, masywne galaktyki spiralne typu pośredniego. Galaktyki te znajdują się po prawej stronie linii regresji, znacznie powyżej przejścia przez zero. Prawo długości koherencji modelu $\ell = c_\text{disk}\,R_d$ wytwarza wartości $\ell$ powyżej 20 kpc w tym reżimie, generując więcej masy pola falowego niż wymaga tego obserwowana rotacja.

4. Dziesięć najbardziej niedoszacowanych galaktyk

Są to galaktyki, dla których BeeTheory przewiduje płaską prędkość rotacji niższą niż obserwowana. Posortowane według wielkości pozostałości:

Galaktyka	Typ Hubble’a	$R_d$ (kpc)	$M_\star/10^{10}$.	$f_\text{gas}$.	$\Sigma_d$	$V_f$	$V_\text{tot}$.	Błąd
NGC6789	Sd-Im	0.30	0.01	0.53	250	60	22	-63.0%
UGC05764	Sd-Im	0.40	0.00	0.86	80	57	31	-45.6%
UGCA442	Sd-Im	1.00	0.00	0.85	15	57	32	-44.2%
NGC4138	S0-Sa	1.30	0.13	0.33	250	150	85	-43.6%
NGC4389	Sb-Sbc	1.20	0.07	0.37	150	110	62	-43.4%
NGC4085	Sb-Sbc	1.20	0.09	0.42	200	135	79	-41.1%
NGC2915	Sd-Im	0.50	0.01	0.84	160	85	53	-38.2%
NGC2976	Sb-Sbc	0.75	0.04	0.29	220	80	50	-37.4%
NGC4183	Sc-Scd	1.60	0.03	0.81	40	110	70	-36.3%
UGCA281	Sd-Im	0.50	0.01	0.63	80	40	26	-36.1%

Własność	Wartość mediany	Zasięg	Porównanie
$R_d$	1.1 kpc	0.30 – 1.80	2 razy $ mniejsza niż mediana
$M_\star$	2,7 razy 10^{8}\,M_\odot$.	$4 \ razy 10^{7}$ – $1.3 \ razy 10^{9}$	6 razy mniej masywny
$f_\text{gas}$.	$0.58$	$0.29$ – $0.86$	Poniżej mediany (0,64)
Hubble $T$	8$ (Sd)	$1$ – $10$	Skoncentrowane w karłach późnego typu
$V_f$	$82$ km/s	$40$ – $150$	Wolne rotatory

Profil niedostatecznie przewidywanej grupy

Kompaktowe karły o niskiej masie i małe spirale. Galaktyki te znajdują się po lewej stronie linii regresji, znacznie poniżej przejścia przez zero. Prawo długości koherencji $\ell = c_\text{disk}\,R_d$ daje $\ell$ rzędu $1$-3$ kpc w tym reżimie, prawdopodobnie zbyt krótkie, aby zebrać pełny zakres pola falowego.

5. Porównanie trzech grup obok siebie

Nieruchomość (mediana)	Przewidywanie zawyżone (err > +30%, $N = 15$)	Dobrze przewidywane (\|err\| ≤ 30%, $N = 67$)	Zaniżone przewidywania (err < -30%, $N = 12$)
$R_d$ (kpc)	4.5	2.4	1.1
$M_\star / 10^{10}$.	1.27	0.15	0.027
$M_\text{gas} / 10^{10}$	0.93	0.27	0.04
$f_\text{gas}$.	0.41	0.64	0.58
$\Sigma_d$	200	140	115
Hubble $T$	5 (Sbc)	6 (Sc)	8 (Sd)
$V_f$ (km/s)	195	113	82

Każda właściwość zmienia się monotonicznie od lewej do prawej. Przewidywana grupa jest większa, bardziej masywna, bardziej zdominowana przez gwiazdy i szybciej się obraca; niedostatecznie przewidywana grupa jest mniejsza, lżejsza, bogatsza w gaz i wolniejsza; dobrze przewidywana większość znajduje się pomiędzy. Droga Mleczna ($R_d = 2,6$ kpc, $V_f około 230$ km/s) naturalnie mieści się w dobrze przewidywanym reżimie, w którym zakotwiczono kalibrację.

6. Interpretacja

Model posiada pojedynczy parametr sprzężenia $\lambda$ oraz trzy uniwersalne stałe geometryczne $(c_\text{disk}, c_\text{sph}, c_\text{arm})$. Zostały one wyznaczone na galaktyce o pośrednim rozmiarze ( Droga Mleczna, $R_d = 2,6$ kpc) i zweryfikowane na dwudziestu dwóch galaktykach o podobnym zakresie rozmiarów. Ślepy test z Uwagi IX pokazuje, że uogólniają się one dość dobrze, ale z pozostałością, która dryfuje liniowo wraz z rozmiarem dysku.

Korekta afiniczna jest wystarczająca

Liniowość reszt w $R_d$ – dobrze dopasowana przez pojedynczą linię prostą przecinającą zero przy $R_d = 2,48$ kpc – jest sygnaturą brakującego addytywnego przesunięcia w relacji długość koherencji. Obecne prawo $\ell = c_\text{disk}\,R_d$ wiąże długość koherencji fali ściśle proporcjonalnie do skali dysku. Zastąpienie go afiniczną zależnością $\ell = c_\text{disk}(R_d – R_0)$, gdzie $R_0$ jest małym przesunięciem około $2.5$ kpc, dałoby resztę, która znika w punkcie kalibracji i rośnie liniowo po obu stronach – dokładnie taki wzorzec, jaki zaobserwowano.

Dobrze przewidywana większość jest ogólnie reprezentatywna

Dwie trzecie próbki mieszczą się w dobrze przewidywanym zakresie. Te 67 galaktyk obejmuje pełen zakres typów Hubble’a i czynnik $\sim 100$ w masie gwiazdowej. Domena ważności modelu nie jest wąska: obejmuje większość populacji SPARC, z odchyleniami skoncentrowanymi na dwóch skrajnych rozmiarach dysku, dokładnie tak, jak wynikałoby to z liniowej pozostałości zależnej od $R_d$.

7. Podsumowanie

1. Błąd predykcji ślepego testu 94 galaktyk podąża za czystym liniowym trendem długości skali dysku: $\text{error}(\%) \approx -31.7 + 12.8\,R_d$, z Pearson $r = +0.75$ i RMSE reszt $= 18.4\%$.

2. Regresja liniowa przecina zero przy $R_d = 2,48$ kpc, zasadniczo rozmiarze dysku Drogi Mlecznej, który zakotwiczył kalibrację. Dwa końce linii odpowiadają dwóm fizycznie odrębnym populacjom odstającym.

3. 15 galaktyk, których przewidywania są zawyżone o więcej niż $+30\%$ to duże, masywne galaktyki spiralne średniego typu: mediana $R_d = 4.5$ kpc, $M_\star \approx 10^{10}\,M_\odot$, $V_f \approx 200$ km/s.

4. 12 galaktyk, których przewidywania są zaniżone o ponad $-30\%$, to zwarte karły o niskiej masie: mediana $R_d = 1,1$ kpc, $M_\star \ około 3 \ razy 10^{8}\,M_\odot$, $V_f \ około 80$ km/s.

5. Odchylenie jest absorbowalne przez afiniczną korektę prawa długości koherencji, $\ell = c_\text{disk}(R_d – R_0)$, z $R_0 \ około 2,5$ kpc – wprowadzając jedną nową stałą.

Referencje. 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. – Third Reference Catalogue of Bright Galaxies, Springer (1991). – McGaugh, S. S. – The third law of galactic rotation, Galaxies 2, 601 (2014). – Dutertre, X. – Bee Theory™: Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).

Anatomia reszt:Liniowy trend z rozmiarem dysku