BeeTheory – Foundations – Uwaga techniczna XVII
Pięć elementów geometrycznych:
Kompletny wykaz parametrów
Przed rozszerzeniem ram BeeTheory na duże próbki galaktyk, niniejsza notatka konsoliduje warstwę modelowania: dla każdego z pięciu komponentów geometrycznych używanych do opisu galaktyki dyskowej, wyraźnie wymienia wymagane parametry, profil gęstości, długość koherencji pola falowego i geometrię integracji. Jest to specyfikacja operacyjna, która napędza każde obliczenie BeeTheory, począwszy od Notatki VII.
1. Najpierw wynik – w skrócie
Na galaktykę: 5 danych obserwacyjnych → 5 składników barionowych → pole falowe
Każda galaktyka jest opisana przez pięć danych obserwacyjnych, które napędzają pięcioskładnikowy rozkład barionowy: wybrzuszenie (3D), cienki dysk (2D), gruby dysk (2D), pierścień gazowy (2D z centralną dziurą) i nadmiar ramienia spiralnego (2D, węższe jądro). Wraz z czterema uniwersalnymi parametrami teorii $(K_0, c_\text{sph}, c_\text{disk}, c_\text{arm})$ i jednym globalnym sprzężeniem $\lambda$, w pełni określa to obliczenia pola falowego.
Parametry całkowite: 5 danych obserwacyjnych + do 18 wyprowadzonych parametrów składowych + 5 parametrów teorii uniwersalnej. Poza tymi parametrami nie ma możliwości dostosowania dla poszczególnych galaktyk.
2. Dane obserwacyjne (na galaktykę)
| Symbol | Ilość | Źródło |
|---|---|---|
| $T$ | Typ morfologiczny Hubble’a | Katalog (de Vaucouleurs, SPARC) |
| $R_d$ | Długość w skali dysku gwiezdnego (kpc) | Fotometria Spitzera 3,6 µm |
| $\Sigma_d$ | Jasność powierzchni dysku centralnego ($L_\odot/\text{pc}^2$) | Fotometria Spitzera 3,6 µm |
| $M_\text{HI}$. | Całkowita masa atomowa wodoru ($M_odot$) | 21-cm obserwacje radiowe |
| $\Upsilon_\star$ | Stosunek masy gwiazdowej do światła przy 3,6 µm | Stały uniwersalny: 0,5\,M_\odot/L_\odot$ (McGaugh 2014) |
Na podstawie tych danych wejściowych obliczane są jednorazowo dwie zintegrowane wielkości masy:
$$M_\star \;=\; 2\pi\,R_d^2\,\Sigma_d\,\Upsilon_\star \qquad\text{(masa gwiazdowa)}$ $$
$$M_\text{gas} \;=\; 1.33\,M_\text{HI} \qquad\text{(masa gazu, poprawka He)}$$
3. Komponent 1 – Wybrzuszenie (3D Hernquist)
Wybrzuszenie to trójwymiarowa sferyczna koncentracja w centrum galaktyki. Jest ono aktywne tylko w przypadku galaktyk wczesnego i średniego typu. W galaktykach spiralnych i nieregularnych późnego typu wybrzuszenie nie występuje.
Aktywacja: $T \leq 4$ (spirale S0, Sa, Sb, Sbc). Wyłączone dla $T \geq 5$ (Sc, Sd, Im).
| Parametr | Symbol | Formuła |
|---|---|---|
| Masa wybrzuszenia | $M_b$ | $0.20 \cdot M_\star$ |
| Promień skali | $r_b$ | $\max(0.5\,R_d,\;0.3\text{ kpc})$ |
| Długość koherencji | $ell_b$ | $c_\text{sph} \cdot r_b$ |
Profil gęstości
$$\rho_b(r) \;=\; \frac{M_b\,r_b}{2\pi\,r\,(r + r_b)^3}$$
Całkowanie pola falowego – powłoki sferyczne
$$\rho_\text{wave}^{(b)}(r) \;=\; \int_0^{6r_b} \rho_b(r’) \cdot K_0\frac{(1 + \alpha_b D)\,e^{-\alpha_b D}}{D^2} \cdot 4\pi r’^2\,dr’, \quad D = \sqrt{r^2 + r’^2}, \quad \alpha_b = 1/\ell_b$$
Liczba parametrów: 3 ($M_b$, $r_b$, $\ell_b$) po aktywacji, 0 w przeciwnym razie.
4. Składnik 2 – Cienki dysk gwiezdny (2D wykładniczy)
Cienki dysk zawiera większość masy gwiezdnej, która nie znajduje się w wybrzuszeniu. Jest to geometrycznie najcieńszy komponent gwiezdny, o najmniejszej rozciągłości pionowej. Zawsze aktywny.
| Parametr | Symbol | Formuła |
|---|---|---|
| Masa cienkiego dysku | $M_\text{thin}$. | $0.75 \cdot (M_\star – M_b) $ |
| Długość skali | $R_d$ | Obserwowane (dane wejściowe) |
| Długość koherencji | $\ell_\text{thin}$. | $c_\text{disk} \cdot R_d$ |
Profil gęstości i integracja
$$\Sigma_\text{thin}(R) \;=\; \frac{M_\text{thin}}{2\pi\,R_d^2}\,e^{-R/R_d}$$
$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{8R_d} \Sigma_\text{thin}(R’) \cdot K_0\frac{(1 + \alpha_\text{thin} D)\,e^{-\alpha_\text{thin} D}}{D^2} \cdot 2\pi R’\,dR’$$.
Liczba parametrów: 3 ($M_\text{thin}$, $R_d$, $\ell_\text{thin}$). Całkowanie po koncentrycznych pierścieniach $R’$.
5. Składnik 3 – Gruby dysk gwiezdny (2D wykładniczy, szerszy)
Gruby dysk składa się ze starszych, dynamicznie cieplejszych gwiazd rozmieszczonych w szerszej skali radialnej niż cienki dysk. Zawsze aktywny. Niesie 25% masy gwiezdnej innej niż kulista.
| Parametr | Symbol | Formuła |
|---|---|---|
| Masa grubego dysku | $M_\text{thick}$. | $0.25 \cdot (M_\star – M_b) $ |
| Długość skali | $R_\text{thick}$. | $1.5 \cdot R_d$ |
| Długość koherencji | $\ell_\text{thick}$. | $c_\text{disk} \cdot R_\text{thick} = 1.5\,c_\text{disk}\,R_d$ |
Profil gęstości i integracja
$$\Sigma_\text{thick}(R) \;=\; \frac{M_\text{thick}}{2\pi\,R_\text{thick}^2}\,e^{-R/R_\text{thick}}$$
$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{8R_\text{thick}} \Sigma_\text{thick}(R’) \cdot K_0\frac{(1 + \alpha_\text{thick} D)\,e^{-\alpha_\text{thick} D}}{D^2} \cdot 2\pi R’\,dR’$$
Liczba parametrów: 3 ($M_\text{thick}$, $R_\text{thick}$, $\ell_\text{thick}$). Taka sama geometria pierścienia jak w przypadku cienkiego dysku.
6. Składnik 4 – Pierścień gazowy (HI + He, 2D z centralnym otworem)
Neutralny gaz atomowy galaktyki (z poprawką na hel) jest rozłożony na szerszą skalę niż dysk gwiazdowy i jest centralnie zubożony. Jest to najbardziej rozległy składnik barionowy, rozciągający się znacznie poza dysk optyczny.
| Parametr | Symbol | Formuła |
|---|---|---|
| Masa gazu | $M_\text{gas}$. | $1.33 \cdot M_\text{HI}$ |
| Długość skali gazu | $R_g$ | $1.7 \cdot R_d$ (Broeils & Rhee 1997) |
| Promień otworu centralnego | $R_\text{hole}$ | $0.5 \cdot R_g$ |
| Długość koherencji | $\ell_\text{gas}$ | $c_\text{disk} \cdot R_g = 1.7\,c_\text{disk}\,R_d$ |
Profil gęstości i integracja
$$\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)$$
$$\rho_\text{wave}^{(\text{gas})}(r) \;=\; \int_{R_\text{hole}}^{8R_g} \Sigma_\text{gas}(R’) \cdot K_0\frac{(1 + \alpha_\text{gas} D)\,e^{-\alpha_\text{gas} D}}{D^2} \cdot 2\pi R’\,dR’$$
Forma podwójnie wykładnicza wychwytuje zarówno centralne zubożenie (człon $-R_\text{hole}/R$ tłumi profil przy małym $R$, gdzie neutralny wodór jest zazwyczaj fotojonizowany lub w postaci molekularnej), jak i zewnętrzny spadek (człon $-R/R_g$). Dolna granica całkowania zaczyna się od $R_\text{hole}$, gdzie profil staje się nieistotny.
Liczba parametrów: 4 ($M_\text{gas}$, $R_g$, $R_\text{hole}$, $\ell_\text{gas}$). Całkowanie pierścieniowe z obciętym promieniem wewnętrznym.
7. Komponent 5 – Nadmiar ramienia spiralnego (2D, węższe jądro)
Ramiona spiralne są azymutalną modulacją gęstości powierzchniowej cienkiego dysku. Są one traktowane, w przybliżeniu osiowo-symetrycznego monopolu BeeTheory, jako efektywne jednorodne wzmocnienie profilu cienkiego dysku na poziomie 10%, ale z wyraźną długością koherencji, która odzwierciedla węższy zakres kątowy struktury ramienia w porównaniu z gładkim dyskiem.
| Parametr | Symbol | Formuła |
|---|---|---|
| Masa efektywna ramienia | $M_\text{arm}$ | $0.10 \cdot M_\text{thin}$ |
| Skala promieniowa | $R_d$ | Podąża za cienkim dyskiem |
| Długość koherencji | $\ell_\text{arm}$. | $c_\text{arm} \cdot R_d$ (węższy niż $\ell_\text{thin}$) |
Profil gęstości i integracja
$$\Sigma_\text{arm}(R) \;=\; 0.10 \cdot \Sigma_\text{thin}(R) \;=\; \frac{0.10\,M_\text{thin}}{2\pi\,R_d^2}\,e^{-R/R_d}$$
$$\rho_\text{wave}^{(\text{arm})}(r) \;=\; \int_0^{8R_d} \Sigma_\text{arm}(R’) \cdot K_0\frac{(1 + \alpha_\text{arm} D)\,e^{-\alpha_\text{arm} D}}{D^2} \cdot 2\pi R’\,dR’$$
Ponieważ $c_\text{arm} < c_\text{disk}$, jądro ramienia spiralnego jest bardziej zlokalizowane niż jądro cienkiego dysku - pole jest wzmocnione w niewielkich odległościach, ale wykładniczo tłumione powyżej kilku kpc. Odzwierciedla to fakt, że prawdziwe ramiona spiralne wytwarzają intensywne lokalne cechy grawitacyjne, ale nie rozszerzają spójności na cały dysk.
Liczba parametrów: 3 ($M_\text{arm}$, $R_d$, $\ell_\text{arm}$). Taka sama geometria pierścienia jak w przypadku cienkiego dysku.
8. Tabela podsumowująca – wszystkie komponenty jednocześnie
| # | Komponent | Geometria | Masa | Skala promieniowa | Spójność $ell$ | Aktywacja | Parametry |
|---|---|---|---|---|---|---|---|
| 1 | Wybrzuszenie | Kula 3D Hernquista | $0.20\,M_\star$ | $r_b = \max(0.5R_d,\,0.3)$. | $c_\text{sph}\,r_b$. | $T \leq 4$ | 3 |
| 2 | Cienki dysk | Wykładniczy 2D | $0.75\,(M_\star – M_b)$. | $R_d$ | $c_\text{disk}\,R_d$ | Zawsze | 3 |
| 3 | Gruby dysk | Wykładniczy 2D | $0.25\,(M_\star – M_b)$. | $1.5\,R_d$ | $1.5\,c_\text{disk}\,R_d$ | Zawsze | 3 |
| 4 | Pierścień gazowy | Eksp. 2D z centralnym otworem | $1.33\,M_\text{HI}$. | $1.7\,R_d$, $R_\text{hole} = 0.85\,R_d$. | $1.7\,c_\text{disk}\,R_d$ | Zawsze | 4 |
| 5 | Ramiona spiralne | Nadwyżka azymutalna 2D | $0.10\,M_\text{thin}$ | $R_d$ (podąża za cienkim) | $c_\text{arm}\,R_d$ | Zawsze | 3 |
9. Parametry teorii uniwersalnej (identyczne dla wszystkich galaktyk)
Pięć liczb określa jądro falowe BeeTheory. Są one uniwersalne – te same wartości mają zastosowanie do Drogi Mlecznej, karłów i masywnych spirali. Nie różnią się one w zależności od galaktyki i są wyznaczane raz na próbce kalibracyjnej.
| Symbol | Wartość | Rola |
|---|---|---|
| $K_0$ | $0.3759$ | Amplituda masy falowej – ustawia bezwymiarową skalę jądra |
| $c_\text{sph}$. | $0.41$ | Współczynnik koherencji 3D: $\ell_b / r_b$ dla źródeł sferycznych (bulge) |
| $c_\text{disk}$ | $3.17$ | Współczynnik koherencji 2D: $\ell / R_\text{scale}$ dla dysków i pierścienia gazowego |
| $c_\text{arm}$ | $2.0$ | Spiralny współczynnik koherencji: węższe jądro dla modulacji ramienia |
| $\lambda$ | $0.4957$ | Globalne sprzężenie pola falowego (skaluje całkowitą gęstość fal) |
Samo jądro fali, identyczne dla każdego komponentu, jest następujące:
$$\mathcal{K}(D) \;=\; K_0 \cdot \frac{(1 + \alpha D)\,e^{-\alpha D}}{D^2}, \qquad \alpha = 1/\ell$$
10. Od komponentów do krzywej rotacji
Całkowita gęstość pola falowego w promieniu $r$ jest sumą pięciu składowych, skalowanych przez globalne sprzężenie:
$$\rho_\text{wave}(r) \;=\; \lambda \cdot \!\!\sum_{i \in \{b,\text{thin},\text{thick},\text{gas},\text{arm}}\!\!\rho_\text{wave}^{(i)}(r)$$
Poniżej przedstawiono masę fali zamkniętej i przewidywaną prędkość kołową:
$$M_\text{wave}(R) \;=\; \int_0^R 4\pi r^2 \rho_\text{wave}(r)\,dr$$
$$V_c^2(R) \;=\; V_\text{bar}^2(R) \;+\; \frac{G\,M_\text{wave}(R)}{R}$$
gdzie $V_\text{bar}(R)$ jest newtonowską prędkością kołową widocznych barionów (wzór Freemana 1970 dla każdego wykładniczego składnika dysku, masa zamknięta Hernquista dla wybrzuszenia, wszystkie połączone w kwadraturze).
11. Księgowanie parametrów – podsumowanie
Na galaktykę, co wchodzi i co jest wyprowadzane
Dane obserwacyjne (na galaktykę): 5 wielkości ($T$, $R_d$, $\Sigma_d$, $M_\text{HI}$, $\Upsilon_\star$).
Pochodne parametry składowe (na galaktykę): 13 jeśli $T > 4$ (bez wybrzuszenia), 16 jeśli $T \leq 4$ (z wybrzuszeniem). Wszystkie obliczone na podstawie 5 powyższych danych wejściowych za pomocą wzorów deterministycznych.
Parametry teorii uniwersalnej: 5 liczb ($K_0$, $c_\text{sph}$, $c_\text{disk}$, $c_\text{arm}$, $\lambda$). Identyczne dla każdej galaktyki.
Darmowe parametry dopasowania dla każdej galaktyki: $\mathbf{0}$. Model nie ma dopasowania galaktyka po galaktyce.
12. Podsumowanie
1. Każda galaktyka jest opisana przez pięć elementów geometrycznych: wybrzuszenie (3D Hernquist, opcjonalnie), cienki dysk gwiezdny, gruby dysk gwiezdny, pierścień gazowy z centralną dziurą i nadmiar ramienia spiralnego (wszystkie cztery ostatnie są wykładnikami 2D).
2. Wybrzuszenie jest aktywowane tylko dla $T \leq 4$ (od S0 do Sbc). Cztery komponenty 2D są zawsze obecne.
3. Każdy komponent wnosi oddzielną całkę do gęstości pola falowego: całkowanie powłoki sferycznej dla wybrzuszenia, całkowanie pierścienia dla czterech komponentów 2D.
4. Liczba wyprowadzonych parametrów na galaktykę wynosi co najwyżej 16 (z wybrzuszeniem) lub 13 (bez wybrzuszenia), wszystkie obliczone deterministycznie na podstawie 5 danych obserwacyjnych.
5. Pięć uniwersalnych parametrów teorii $(K_0, c_text{sph}, c_text{disk}, c_text{arm}, lambda)$ jest identycznych dla wszystkich galaktyk – nie są one dostosowywane dla każdej galaktyki.
6. W modelu nie istnieje żaden parametr wolny dla galaktyki. Po ustaleniu $lambda$ na próbce kalibracyjnej, wszystkie kolejne krzywe rotacji są czystymi przewidywaniami.
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). – Hernquist, L. – An analytical model for spherical galaxies and bulges, ApJ 356, 359 (1990). Profil gęstości wybrzuszenia. – Freeman, K. C. – On the disks of spiral and S0 galaxies, ApJ 160, 811 (1970). Wykładnicza prędkość kołowa dysku. – Broeils, A. H., Rhee, M.-H. – Short 21-cm WSRT observations of spiral and irregular galaxies, A&A 324, 877 (1997). Stosunek skali gazu do dysku gwiezdnego. – McGaugh, S. S. – Trzecie prawo rotacji galaktyk, Galaxies 2, 601 (2014). $\Upsilon_\star$ przy 3,6 µm. – Bland-Hawthorn, J., Gerhard, O. – The Galaxy in Context, ARA&A 54, 529 (2016). Rozkład strukturalny Drogi Mlecznej. – Dutertre, X. – Bee Theory™: Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).
BeeTheory.com – Kwantowa grawitacja oparta na falach – Warstwa modelująca – © Technoplane S.A.S. 2026