BeeTheory – Foundations – Uwaga techniczna XII
Formalizacja:
Obliczenia teorii pszczół w skali galaktycznej
Niniejsza notatka formalizuje teorię BeeTheory zastosowaną do galaktyki dyskowej. Określa ona dane obserwacyjne, geometryczną dekompozycję rozkładu barionowego, równania całkowe definiujące pole falowe dla każdego składnika oraz łańcuch operacji, które dają przewidywaną krzywą rotacji. Procedura jest ściśle jednokierunkowa: obserwowana struktura barionowa określa pole falowe, które określa krzywą rotacji – nigdy odwrotnie.
1. Obliczenia na jednym wykresie
Łańcuch jednokierunkowy
Obserwowana fotometria $\;\longrightarrow\;$ Rozkład barionowy $(\rho_\text{bar})$
$\big\downarrow$
Sprzężenie pola falowego $\;\longrightarrow\;$ Gęstość falowa $(\rho_\text{wave})$
$\big\downarrow$
Całkowanie masy $\;\longrightarrow\;$ Zamknięta masa fali $(M_\text{wave})$
$\big\downarrow$
Relacja newtonowska $\;\longrightarrow\;$ Przewidywana krzywa rotacji $(V_c)$
Żaden krok nie jest odwracany. Krzywa rotacji $V_c(R)$ nigdy nie jest używana jako dane wejściowe.
2. Dane wejściowe z obserwacji
Dla każdej galaktyki obliczenia wymagają pięciu opublikowanych obserwabli. Są to jedyne wielkości specyficzne dla galaktyki; wszystko inne jest obliczane na ich podstawie. Na tym etapie nie jest wykonywane dopasowanie do krzywej rotacji.
| Symbol | Ilość | Źródło |
|---|---|---|
| $T$ | Typ morfologiczny Hubble’a | Katalog (de Vaucouleurs et al. 1991, SPARC) |
| $R_d$ | Długość w skali dysku gwiezdnego (kpc) | Fotometria Spitzera 3,6 µm (SPARC) |
| $\Sigma_d$ | Jasność powierzchni dysku centralnego ($L_\odot/\text{pc}^2$) | Fotometria Spitzera 3,6 µm (SPARC) |
| $M_\text{HI}$. | Całkowita masa atomowa wodoru ($M_odot$) | 21-cm obserwacje radiowe (SPARC) |
| $\Upsilon_\star$ | Stosunek masy gwiazdowej do światła przy 3,6 µm | Stały uniwersalny: 0,5\,M_\odot/L_\odot$ (McGaugh 2014) |
Dla Drogi Mlecznej, $R_d$, $Sigma_d$ i $M_text{HI}$ zostały zastąpione analogicznymi wartościami wyznaczonymi z wewnętrznych przeglądów gwiazd (Bovy & Rix 2013) i map 21-cm. Używany jest ten sam pięcioilościowy wektor wejściowy.
3. Dekompozycja barionowa – pięć składników geometrycznych
Na podstawie pięciu danych obserwacyjnych masa barionowa jest dzielona na pięć odrębnych składników geometrycznych. Każdy składnik ma swój własny profil gęstości i charakterystyczną skalę.
3.1 Całkowita masa gwiazd i gazu
$$M_\star \;=\; 2\pi\, R_d^2\, \Sigma_d\, \Upsilon_\star$$.
$$M_\text{gas} \;=\; 1.33\,M_\text{HI} \qquad \text{(korekta He; Arnett 1996)}$$
3.2 Masy i wagi składników
| Komponent | Masa | Skala | Aktywacja |
|---|---|---|---|
| Wybrzuszenie | $M_b = 0,20\,M_\star$. | $r_b = \max(0.5\,R_d,\,0.3\text{kpc})$. | Jeśli $T \leq 4$ |
| Cienki dysk | $M_\text{thin} = 0.75\,(M_\star – M_b)$. | $R_d$ | Zawsze |
| Gruby dysk | $M_\text{thick} = 0.25\,(M_\star – M_b)$. | $1.5\,R_d$ | Zawsze |
| Pierścień gazowy | $M_\text{gas} = 1,33\,M_\text{HI}$. | $R_g = 1,7\,R_d$ (Broeils & Rhee 1997) | Zawsze |
| Ramiona spiralne | $M_\text{arm} = 0,10\,M_\text{thin}$ (efektywnie) | $R_d$ (podąża za cienkim dyskiem) | Zawsze |
3.3 Profile gęstości
Wybrzuszenie (3D Hernquista)
$$\rho_b(r) \;=\; \frac{M_b\,r_b}{2\pi\,r\,(r + r_b)^3}$$
Cienkie i grube dyski gwiezdne (2D wykładnicze)
$$\Sigma_\text{thin}(R) \;=\; \frac{M_\text{thin}}{2\pi\,R_d^2}\,e^{-R/R_d}$$
$$\Sigma_\text{thick}(R) \;=\; \frac{M_\text{thick}}{2\pi\,(1.5\,R_d)^2}\,e^{-R/(1.5R_d)}$$
Pierścień gazowy (wykładniczy 2D z centralnym otworem)
$$\Sigma_\text{gas}(R) \;=\; \frac{M_\text{gas}}{2\pi\,R_g^2}\,\exp\!\left(-\frac{R_\text{hole}}{R} – \frac{R}{R_g}\right), \quad R_\text{hole} = 0.5\,R_g$$
Nadmiar ramienia spiralnego (2D, podąża za cienkim dyskiem)
$$\Sigma_\text{arm}(R) \;=\; 0.10\;\Sigma_\text{thin}(R)$$
4. Jądro fali
Każdy element masy barionowej generuje pole falowe BeeTheory. Pole w punkcie $vec{r}$ wytwarzane przez element źródłowy w punkcie $vec{r},’$ oddzielonym przez $D = |vec{r} – vec{r},’|$ jest regulowane przez jądro typu Yukawy wyprowadzone z uregulowanej funkcji falowej z Uwagi I:
Jądro fali teorii pszczół
$$\mathcal{K}_i(D) \;=\; K_0\,\frac{(1 + \alpha_i\,D)\,e^{-\alpha_i\,D}}{D^2}, \qquad \alpha_i \;=\; \frac{1}{\ell_i}$$.
Tutaj $K_0$ jest uniwersalną amplitudą fali-masy (pojedyncza liczba bezwymiarowa), a $\ell_i$ jest długością koherencji składnika $i$. Jądro koduje quasi-newtonowskie zachowanie $1/D^2$ w małych odległościach, modulowane przez wykładnicze odcięcie w skalach powyżej $\ell_i$. Postać $(1 + \alpha D)\,e^{-\alpha D}$ zapewnia ciągłość i skończoną całkowitą masę zamkniętą w nieskończoności.
4.1 Długości spójności komponentów
Długość koherencji każdego komponentu jest ustalana przez jego naturalną skalę geometryczną, pomnożoną przez bezwymiarową stałą specyficzną dla jego wymiarowości:
| Komponent | Długość koherencji | Stała geometryczna |
|---|---|---|
| Wybrzuszenie (sfera 3D) | $\ell_b = c_\text{sph}\,r_b$ | $c_\text{sph}$. |
| Cienki dysk (2D) | $\ell_\text{thin} = c_\text{disk}\,R_d$ | $c_\text{disk}$ |
| Gruby dysk (2D) | $\ell_\text{thick} = c_\text{disk}\,(1.5\,R_d)$ | $c_\text{disk}$ |
| Pierścień gazowy (2D) | $\ell_\text{gas} = c_\text{disk}\,R_g$ | $c_\text{disk}$ |
| Ramiona spiralne (2D, skoncentrowane azymutalnie) | $\ell_\text{arm} = c_\text{arm}\,R_d$ | $c_\text{arm}$ |
Trzy stałe geometryczne $(c_\text{sph},\,c_\text{disk},\,c_\text{arm})$ są uniwersalne – nie różnią się w zależności od galaktyki. Wraz z globalną amplitudą masy falowej $K_0$ i sprzężeniem pola falowego $\lambda$ stanowią one pełny zestaw parametrów na poziomie teorii.
5. Sprzężenie pola falowego – równania całkowe na składową
Gęstość pola falowego w pozycji $\vec{r}$ jest splotem rozkładu źródła barionowego z jądrem falowym. W przypadku galaktycznie symetrycznego układu (osiowo symetrycznego, w przybliżeniu monopolarnego), każdy składnik barionowy wnosi addytywny wkład:
Całkowita gęstość pola falowego w promieniu $r$
$$\rho_\text{wave}(r) \;=\; \lambda \;\\sum_{i \in \{\text{cienki, gruby, gaz, ramię, wybrzuszenie}}} \rho_\text{wave}^{(i)}(r)$$.
Poniżej przedstawiono pięć całek, po jednej na składnik. Każda całka przekształca rozkład masy barionowej w rozkład masy pola falowego w tym samym punkcie przestrzennym.
5.1 Wybrzuszenie – integracja powłoki 3D
$$\rho_\text{wave}^{(b)}(r) \;=\; \int_0^{r_\text{max}} \rho_b(r’)\;\mathcal{K}_b\!\left(\sqrt{r^2 + r’^2}\right)\;4\pi r’^2\,dr’$$
Całkowanie odbywa się po koncentrycznych powłokach sferycznych o promieniu $r’$. Punkt pola w promieniu $r$ od środka widzi każdą powłokę w efektywnej separacji $D = \sqrt{r^2 + r’^2}$ w przybliżeniu monopolarnym. Całkowanie rozciąga się do $r_\text{max} = 6\,r_b$, poza którym gęstość wybrzuszenia jest numerycznie zaniedbywalna.
5.2 Cienki dysk – integracja pierścienia 2D
$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{thin}(R’) \;\mathcal{K}_\text{thin}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
Dysk jest podzielony na koncentryczne pierścienie o promieniu $R’$ i nieskończenie małej szerokości $dR’$, z których każdy ma masę powierzchniową $\Sigma_\text{thin}(R’)\,2\pi R’\,dR’$. Obowiązuje to samo przybliżenie monopolarne: pole falowe w promieniu $r$ od środka otrzymuje wkład od każdego pierścienia w efektywnej separacji $D = \sqrt{r^2 + R’^2}$. Zakres całkowania wynosi $R_\text{max} = 8\,R_d$.
5.3 Gruby dysk – integracja pierścienia 2D
$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{thick}(R’) \;\mathcal{K}_\text{thick}\!\left(\sqrt{r^2 + R’^2}\right) \;2\pi R’\,dR’$$
Identycznie jak w przypadku całkowania cienkiego dysku, z $\Sigma_\text{thick}(R’)$ jako gęstością źródła i parametrem jądra $\alpha_\text{thick} = 1/(c_\text{disk}\,\cdot 1.5\,R_d)$. Szerszy zasięg radialny grubego dysku skutkuje nieco szerszym zakresem koherencji fal.
5.4 Pierścień gazowy – integracja pierścienia 2D z centralnym zubożeniem
$$\rho_\text{wave}^{(\text{gas})}(r) \;=\; \int_{R_\text{hole}}^{R_\text{max}} \Sigma_\text{gas}(R’)\;\mathcal{K}_\text{gas}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
Rozkład gazu ma centralną dziurę, uchwyconą przez radialne odcięcie przy $R_\text{hole} = 0.5\,R_g$ w dolnej granicy całkowania. Poza tym odcięciem gaz rozciąga się dalej niż dysk gwiezdny; znajduje to odzwierciedlenie w większej skali charakterystycznej $R_g = 1.7\,R_d$, która wpływa na długość koherencji $\ell_\text{gas} = c_\text{disk}\,R_g$.
5.5 Nadmiar ramienia spiralnego – integracja pierścienia 2D ze zmniejszoną amplitudą
$$\rho_\text{wave}^{(\text{arm})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{arm}(R’)\;\mathcal{K}_\text{arm}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
Ramiona spiralne są traktowane jako uśrednione osiowo wzmocnienie gęstości powierzchniowej cienkiego dysku na poziomie $10\%$, z ich własną długością koherencji $\ell_\text{arm} = c_\text{arm}\,R_d$. Jądro jest zatem węższe niż jądro cienkiego dysku, odzwierciedlając azymutalną koncentrację struktury spiralnej.
6. Zamknięta masa falowa i przewidywana krzywa rotacji
Gdy znana jest całkowita gęstość pola falowego $\rho_\text{wave}(r)$, zamkniętą masę pola falowego w kuli o promieniu $R$ uzyskuje się przez całkowanie promieniowe:
Zamknięta masa pola falowego
$$M_\text{wave}(R) \;=\; \int_0^{R} 4\pi\,r^2\,\rho_\text{wave}(r)\,dr$$
Przewidywana prędkość kołowa w promieniu $R$ wynika następnie z relacji Newtona, łącząc wkłady pola barionowego i falowego w kwadraturze:
Przewidywana prędkość kołowa
$$V_c^2(R) \;=\; V_\text{bar}^2(R) \;+\; \frac{G\,M_\text{wave}(R)}{R}$$
Prędkość barionowa $V_\text{bar}(R)$ jest sama w sobie kwadratową sumą wkładów z czterech składników dysko-podobnych (wzór Freemana 1970 dla każdego profilu wykładniczego) i wybrzuszenia (wzór Hernquista na masę zamkniętą):
$$V_\text{bar}^2(R) \;=\; V_\text{bulge}^2 + V_\text{thin}^2 + V_\text{thick}^2 + V_\text{gas}^2$$.
gdzie każde $V_i(R)$ jest standardową prędkością kołową Newtona odpowiedniego rozkładu masy.
7. Parametry na poziomie teorii
Kompletna struktura BeeTheory, zastosowana do galaktyk, zawiera pięć parametrów na poziomie teorii. Są one uniwersalne: nie różnią się w zależności od galaktyki.
| Symbol | Znaczenie | Rola |
|---|---|---|
| $K_0$ | Amplituda masy falowej | Ustawia bezwymiarową skalę jądra fali |
| $c_\text{sph}$. | Stała geometryczna 3D | Stosunek $\ell/r_\text{scale}$ dla źródeł sferycznych (wybrzuszenie) |
| $c_\text{disk}$ | Stała geometryczna 2D | Stosunek $\ell/R_\text{scale}$ dla źródeł dyskowych i pierścieniowych |
| $c_\text{arm}$ | Stała geometryczna spirali | Stosunek $\ell/R_d$ dla azymutalnie skoncentrowanego nadmiaru ramienia |
| $\lambda$ | Globalne sprzężenie pola falowego | Skaluje całkowitą gęstość pola falowego |
Uniwersalność parametrów
Wszystkie pięć parametrów ma charakter globalny. Te same wartości liczbowe mają zastosowanie do Drogi Mlecznej, nieregularnych karłów i masywnych spirali. Informacje specyficzne dla galaktyki są wprowadzane tylko poprzez pięć danych obserwacyjnych $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star) $. Model nie zawiera żadnego przestrajalnego parametru na galaktykę.
8. Jednokierunkowy charakter obliczeń
Otwarty łańcuch – brak sprzężenia zwrotnego
Całe obliczenie przepływa od wejść do wyjść, w jednym kierunku. Obserwacje fotometryczne i 21-cm określają rozkład barionowy. Rozkład barionowy określa gęstość pola falowego. Gęstość pola falowego określa masę fali zamkniętej. Zamknięta masa falowa określa przewidywaną krzywą rotacji. W żadnym momencie krzywa rotacji nie wpływa na żaden wcześniejszy etap obliczeń.
Ta jednokierunkowość ma trzy ważne konsekwencje.
(a) Po ustaleniu pięciu parametrów na poziomie teorii, krzywa rotacji jest ścisłym przewidywaniem, a nie dopasowaniem. Porównanie z obserwowaną krzywą rotacji jest testem, a nie kalibracją.
(b) Model nie ma mechanizmu dostosowywania galaktyka po galaktyce. Każda modyfikacja przewidywanej krzywej rotacji musi wynikać z modyfikacji wektora wejściowego $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$ lub ze zmiany parametrów na poziomie teorii uniwersalnej $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$.
(c) Kalibracja $lambda$ na galaktyce referencyjnej nie jest tym samym, co dopasowanie jej do krzywej rotacji tej galaktyki. Kalibracja określa pojedynczą liczbę globalną; krzywa rotacji na wszystkich innych promieniach galaktyki referencyjnej oraz krzywe rotacji wszystkich innych galaktyk są następnie ścisłymi przewidywaniami skalibrowanej struktury.
9. Rola centralnej gęstości powierzchniowej (Uwaga XI rewizja)
Diagnoza z Uwagi XI wykazała, że resztkowy błąd przewidywania silnie koreluje z centralną barionową gęstością powierzchniową $\Sigma_d$, niezależnie od długości skali dysku $R_d$. Przedstawiona powyżej formalizacja jest wersją modelu przed uwzględnieniem tego odkrycia – wykorzystuje ona tylko $R_d$ w wyrażeniach długości koherencji $\ell_i = c_i\,R_d$.
Gdzie pojawi się udoskonalenie
W udoskonalonym modelu, długości koherencji $\ell_i$ będą zależeć zarówno od $R_d$ jak i $\Sigma_d$, zastępując ścisłą liniową zależność $\ell_i = c_i\,R_d$ funkcją $\ell_i = c_i\,R_d\,\phi(\Sigma_d/\Sigma_\text{ref})$, która absorbuje resztki zidentyfikowane w Nocie XI. Postać funkcyjna $phi$ i jej parametry zostaną określone w kolejnych notatkach, najpierw na zbiorze kalibracyjnym 22 galaktyk, a następnie zweryfikowane przez ślepą prognozę na pozostałej próbce SPARC.
Jednokierunkowa struktura obliczeń jest zachowana przez to udoskonalenie: $\Sigma_d$ jest wejściem obserwacyjnym, zmodyfikowane długości koherencji zasilają te same całki splotowe, a krzywa rotacji pojawia się jak poprzednio. Dodano tylko jedno powiązanie operacyjne – zależność $\ell_i$ od drugiej obserwowalnej.
10. Podsumowanie metodologii
1. Dane wejściowe. Pięć obserwabli na galaktykę: typ Hubble’a $T$, skala dysku $R_d$, jasność powierzchniowa $\Sigma_d$, masa HI $M_\text{HI}$ oraz uniwersalny stosunek masy gwiazdowej do światła $\Upsilon_\star$.
2. Rozkład barionowy. Pięć składników: wybrzuszenie (jeśli $T \leq 4$), cienki dysk, gruby dysk, pierścień gazowy, nadmiar ramienia spiralnego. Każdy z nich posiada analityczny profil gęstości.
3. Jądro falowe. Uniwersalna postać typu Yukawy $\mathcal{K}_i(D) = K_0\,(1 + \alpha_i D)\,e^{-\alpha_i D}/D^2$ z długością koherencji $\ell_i = c_i\,R_\text{scale}$ określoną przez geometryczny zasięg każdego komponentu.
4. Konwolucja. Każdy komponent generuje gęstość pola falowego poprzez jednowymiarową całkę nad pierścieniami (komponenty 2D) lub powłokami (wybrzuszenie 3D). Całkowita gęstość pola falowego jest sumą pięciu komponentów, skalowaną przez globalne sprzężenie $\lambda$.
5. Wynik. Dołączona masa falowa $M_\text{wave}(R)$ jest całkowana i łączona z prędkością barionową $V_\text{bar}(R)$, aby uzyskać przewidywaną krzywą rotacji $V_c(R)$.
6. Parametry na poziomie teorii. $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$ – uniwersalny, bez dostrajania na galaktykę. Badane udoskonalenie doda zależność od $\Sigma_d$.
7. Kierunek. Wejścia → bariony → pole falowe → krzywa rotacji. Brak sprzężenia zwrotnego. Krzywa rotacji jest przewidywaniem, a nie dopasowaniem.
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). – Freeman, K. C. – On the disks of spiral and S0 galaxies, ApJ 160, 811 (1970). – Hernquist, L. – An analytical model for spherical galaxies and bulges, ApJ 356, 359 (1990). – Broeils, A. H., Rhee, M.-H. – Short 21-cm WSRT observations of spiral and irregular galaxies, A&A 324, 877 (1997). – McGaugh, S. S. – Trzecie prawo rotacji galaktyk, Galaxies 2, 601 (2014). – Bovy, J., Rix, H.-W. – A direct dynamical measurement of the Milky Way’s disk surface density profile, ApJ 779, 115 (2013). – Arnett, D. – Supernovae and Nucleosynthesis, Princeton (1996). – Dutertre, X. – Bee Theory™: Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).
BeeTheory.com – Kwantowa grawitacja oparta na falach – Metodologia galaktyczna – © Technoplane S.A.S. 2026