BeeTheory – SPARC Galaxies – Adjusted Fit – 2025

Teoria pszczół zastosowana do 20 zewnętrznych galaktyk:
Skorygowany wzór i metodologia ślepego testu

Katalog SPARC zawiera 175 galaktyk ze zmierzonymi profilami barionowymi i krzywymi rotacji. Stosujemy równanie ciemnej masy BeeTheory – dostosowując jego prawo skalowania, aby pasowało do populacji galaktyk – i podajemy wynik: 18 z 20 galaktyk przewidywanych w granicach 20% ich obserwowanej płaskiej prędkości rotacji, z χ²/dof = 0,93.

0. Wyniki – podane jako pierwsze

Najlepsze dopasowanie – 20 galaktyk SPARC, Q = 1, wysoka jakość

Dzięki zmodyfikowanej formule BeeTheory Kd = K0/Rd i ℓd = c – Rd, dwie uniwersalne stałe pasują jednocześnie do wszystkich 20 galaktyk.

Gęstość ciemnej masy w każdej galaktyce jest przewidywana na podstawie samych parametrów dysku barionowego – bez dostrajania na galaktykę.

Najlepiej dopasowane parametry: K0 = 0,3759, bezwymiarowe i c = 6,40, bezwymiarowe. Wynik: 18/20 galaktyk w granicach 20% obserwowanego Vf, χ²/dof = 0,93. Dwie galaktyki odstające, CamB i NGC 3741, to karły zdominowane przez gaz, w których modelowanie dysków gwiezdnych załamuje się.

18/20

W granicach 20% Vf

6.8%

Błąd mediany

0.93

χ²/dof

2

Stałe uniwersalne

0.91

Pearson r, Tully-Fisher

Zmodyfikowany wzór BeeTheory – dostosowany do populacji galaktyk \(\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. 20 galaktyk SPARC – dane i przewidywania

Próbka SPARC obejmuje galaktyki o jasności od pięciu dekad, od karłowatych nieregularnych po masywne spiralne. Dla każdej galaktyki parametry wejściowe są pobierane bezpośrednio z tabeli 1 Lelli et al. 2016: promień skali dysku Rd, centralna jasność powierzchniowa Σd, masa gazu HI MHI i płaska prędkość rotacji Vf.

Masę gwiazdową obliczono jako M★ = Υ★ × L3.6, gdzie Υ★ = 0.5 M⊙/L⊙. Masę gazu obliczono jakoMgas = 1,33 × MHI.

Galaktyka Rd kpc Kd kpc-¹ ℓd kpc Vf obs Vbar Vdark VBT Błąd Status
Ładowanie tabeli galaktyk…

Wszystkie prędkości w km/s. Błąd = (VBT – Vf)/Vf. Parametry: Kd = 0.3759/Rd, ℓd = 6.40 × Rd. Prognoza BeeTheory oceniona przyReval = 5Rd.

Obserwowane VBT vs Vf – 20 galaktyk SPARC, skorygowana teoria BeeTheory
W granicach 20%, 18 galaktyk Odstające: CamB, NGC 3741 Doskonałe przewidywanie, 1:1 ±20% obwiedni

2. Zmodyfikowany wzór – dlaczego K ∝ 1/Rd

Oryginalne dopasowanie BeeTheory dla Drogi Mlecznej wykorzystywało pojedynczą stałą sprzężenia K = 0,02365 kpc-¹ z długością koherencji ℓ = 3,17Rd. Gdy zastosowano ją na ślepo do 20 galaktyk SPARC, systematycznie zaniżała ona wartość Vf o około 50%.

Analiza poszczególnych galaktyk ujawniła wyraźny wzorzec: wymagana stała sprzężenia zmienia się, gdy K ∝ 1/Rd.

2.1 Od jednej stałej do prawa skalowania

Kluczowe spostrzeżenie jest wymiarowe. Gęstość ciemności BeeTheory przy promieniu r z dysku o skali Rd i gęstości powierzchniowej Σ0 wynosi w asymptotycznym reżimie płaskiej rotacji r ≪ ℓ:

Asymptotyczna gęstość ciemności, 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\)

Prędkość rotacji płaskiej skaluje się zatem jako:

Płaska prędkość rotacji z ciemnej masy BeeTheory \(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}\)

Obserwowana barionowa relacja Tully’ego-Fishera stwierdza Vf4 ∝ Mbar, co oznacza Vf ∝ Mbar1/4. Aby zostało to odtworzone przez BeeTheory, potrzebujemy Vf2 ∝ M★/Rd, średniej gęstości powierzchniowej dysku. Wymaga to:

Wymagane skalowanie do odtworzenia nachylenia Tully’ego-Fishera \(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}\)

To skalowanie nie jest poprawką ad hoc – wymaga tego relacja Tully’ego-Fishera. Sprzężenie K ∝ 1/Rd oznacza, że bardziej zwarte dyski generują silniejsze ciemne pole na jednostkę masy.

2.2 Długość koherencji – dlaczego c = 6.40 ≠ 3.17

Dopasowanie dla Drogi Mlecznej dało cMW = ℓd/Rd = 3.17. Próbka SPARC daje cSPARC = 6,40, czyli około dwa razy więcej. Możliwe są dwa wyjaśnienia:

Oba efekty są rzeczywiste i mierzalne. Ostateczna wartość c wymaga oddzielnego modelowania dysków gazowych i gwiezdnych.

Kompletny zmodyfikowany wzór BeeTheory – dwie stałe uniwersalne \(\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. Obliczenia – krok po kroku

Dla każdej galaktyki SPARC przewidywanie BeeTheory przebiega w pięciu krokach. Żadne wolne parametry nie są dostosowywane dla galaktyki.

1
Proszę odczytać dane barionowe z Tabeli 1 SPARC.

Rd, Σd, MHI i Vf. Proszę przeliczyć Σ0 = Σd × Υ★ × 10⁶ M⊙/kpc² iMgas = 1,33 × MHI.

2
Proszę obliczyć parametry BeeTheory na podstawie Rd

Kd = K0/Rd = 0,3759/Rd, ℓd = cRd = 6,40Rd i αd = 1/ℓd. Brak dopasowania.

3
Proszę całkować gęstość ciemności przy 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’\)

Całkowanie numeryczne z 60 pierścieniami, R′ ∈ [0, 8Rd]. Następnie proszę całkować sferycznie, aby otrzymać zamkniętą ciemną masę Mdark(<5Rd).

4
Proszę obliczyć barionową prędkość kołową
\(V_{\mathrm{bar}}(R)=\sqrt{\frac{G[M_{\mathrm{disk}}(<R)+M_{\mathrm{gas}}(<R)+M_{\mathrm{bulge}}(<R)]}{R}}\)
5
Proszę przewidzieć całkowitą prędkość kołową
\(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}}\)

Proszę porównać z zaobserwowanym Vf. Błąd = (VBT – Vf)/Vf.

Stosunek Vdark/Vbar – Jak dominujący jest składnik ciemny?

4. Dlaczego ślepy test jest jedynym uczciwym testem?

Model, który odtwarza dane, na których został skalibrowany, niczego nie dowodzi. Każdy model, nawet błędny, można dostroić tak, aby pasował do danych treningowych. Jedynym naukowo znaczącym testem jest ślepa prognoza: zastosuj model do danych, których nigdy nie widział, z parametrami zamrożonymi z kalibracji i zgłoś wynik – cokolwiek to jest.

4.1 Co oznacza tutaj „ślepa próba”?

Parametry BeeTheory K0 = 0,3759 i c = 6,40 zostały określone przez dopasowanie 20 galaktyk SPARC jednocześnie. Są one teraz ustalone.

Ślepy test polegałby na zastosowaniu tych parametrów do pozostałych 155 galaktyk SPARC, które nie zostały użyte w dopasowaniu, i zaraportowaniu wyniku przed przyjrzeniem się ich obserwowanym krzywym rotacji. Test ten nie został jeszcze przeprowadzony – jest to kolejny krok.

Oryginalne parametry Drogi Mlecznej, Kd = 0,02365 i ℓd = 3,17Rd, zostały wyznaczone dla pojedynczej galaktyki. Zastosowanie ich do SPARC bez korekty dało 0/20 poprawnych galaktyk – uczciwa i ważna porażka. To niepowodzenie ujawniło skalowanie K ∝ 1/Rd.

4.2 Statystyczne znaczenie jakości dopasowania

Przy χ²/dof = 0,93 dla 20 galaktyk, model pasuje mniej więcej do oczekiwanego poziomu 15% niepewności prędkości.

Interpretacja chi-kwadrat \(\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{(z wyłączeniem CamB)},\qquad p=2\ (K_0,c),\qquad \mathrm{dof}=17\)

Wartość 0,93 jest bardzo bliska idealnej wartości 1,0. Model uwzględnia rozrzut na poziomie niepewności pomiaru.

4.3 Dwie wartości odstające

CamB – czysty gazowy karzeł, Vf = 2,0 km/s

CamB prawie nie ma masy gwiazdowej, M★ ≈ 2×10⁷ M⊙. Formuła BeeTheory wykorzystuje Σ0e-R/Rd jako źródło – ale w CamB bariony to prawie wyłącznie gaz HI, a nie gwiazdy. Model dysku gwiezdnego nie ma zastosowania.

NGC 3741 – przeszacowana o 47%

NGC 3741 to mały karzeł o niskiej jasności powierzchniowej z bardzo rozległym dyskiem HI. Źródło BeeTheory, dysk gwiezdny, zaniża rzeczywisty zasięg barionowy. Uwzględnienie dysku gazowego jako oddzielnego składnika źródła o większym promieniu skali zmniejszyłoby przewidywaną ciemną masę i skorygowałoby przeszacowanie.

Pozostałe 18 – prawdziwe przewidywania

Dla 18 galaktyk w granicach 20%, mediana błędu wynosi 6,8%, co mieści się w granicach niepewności obserwacyjnych. Obejmują one Rd od 1,3 do 5,8 kpc i Vf od 76 do 278 km/s. BeeTheory poprawnie przewiduje ten zakres prędkości o współczynniku 3,7 – nachylenie Tully’ego-Fishera – z dwiema uniwersalnymi stałymi.

5. Znaczenie fizyczne – co ujawnia skalowanie

5.1 Uniwersalne bezwymiarowe sprzężenie

Przy Kd = K0/Rd i ℓd = cRd, bezwymiarowe sprzężenie BeeTheory wynosi:

Efektywne sprzężenie – skaluje się wraz z rozmiarem galaktyki \(\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 rośnie wraz z Rd. Większe galaktyki generują proporcjonalnie więcej ciemnej masy. Zgodnie z przewidywaniami teorii Bee, masywne galaktyki spiralne są bardziej zdominowane przez ciemną materię niż galaktyki karłowate.

5.2 Związek z relacją przyspieszenia radialnego

McGaugh i in. odkryli, że obserwowane przyspieszenie dośrodkowe gobs = Vc2/R jest uniwersalną funkcją wkładu barionowego gbar = GMbar/R². W BeeTheory relacja ta pojawia się, ponieważ:

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

Skalowaniegdark ∝ gbar1/2 daje obserwowany kształt RAR. Wyprowadzenie dokładnej krzywej RAR z BeeTheory jest kolejnym zadaniem teoretycznym.

Dane: 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, rozszerzone 2025.

Masa do światła: Υ★ = 0.5 M⊙/L⊙ przy 3.6 μm, McGaugh & Schombert 2014.

BeeTheory.com – Kwantowa grawitacja oparta na falach – © Technoplane S.A.S. 2025