BeeTheory – Foundations – Teknisk anvisning XII
Formalisering:
BeeTheory-beräkning i galaktisk skala
Denna not formaliserar BeeTheory-ramverket när det tillämpas på en skivgalax. Här specificeras observationsdata, den geometriska uppdelningen av baryonfördelningen, integralekvationerna som definierar vågfältet för varje komponent och den kedja av operationer som ger den förutsagda rotationskurvan. Proceduren är strikt enkelriktad: den observerade baryonstrukturen bestämmer vågfältet, som bestämmer rotationskurvan – aldrig tvärtom.
1. Beräkningen i ett diagram
En enkelriktad kedja
Observerad fotometri $\;\longrightarrow\;$ Baryonisk sönderdelning $(\rho_\text{bar})$
$\stor\nedåtpil$
Vågfältskonvolution $\;\longrightarrow\;$ Vågtäthet $(\rho_\text{wave})$\big\downarrow
$\big\nedåtpil$
Massintegration $\;\longrightarrow\;$ Innesluten vågmassa $(M_\text{våg})$\big\downarrow
$\big\nedåtpil$
Newtonskt förhållande $\;\longrightarrow\;$ Förutsedd rotationskurva $(V_c)$
Inget steg är inverterat. Rotationskurvan $V_c(R)$ används aldrig som indata.
2. Observationella ingångar
För varje galax kräver beräkningen fem publicerade observabler. Dessa är de enda galaxspecifika storheterna; allt annat beräknas utifrån dem. Ingen anpassning mot rotationskurvan utförs i detta skede.
| Symbol | Kvantitet | Källa |
|---|---|---|
| $T$ | Hubbles morfologiska typ | Katalog (de Vaucouleurs m.fl. 1991, SPARC) |
| $R_d$ | Stellar disk skalans längd (kpc) | Spitzer 3,6 µm fotometri (SPARC) |
| $\Sigma_d$ | Ljusstyrka på den centrala skivans yta ($L_\odot/\text{pc}^2$) | Spitzer 3,6 µm fotometri (SPARC) |
| $M_\text{HI}$$ | Total atomär vätemassa ($M_odot$) | 21-cm radioobservationer (SPARC) |
| $\Upsilon_\stjärna$ | Stjärnans förhållande mellan massa och ljus vid 3,6 µm | Fast universal: $0.5\,M_\odot/L_\odot$ (McGaugh 2014) |
För Vintergatan ersätts $R_d$, $Sigma_d$ och $M_text{HI}$ med de analoga värden som bestämts från interna stjärnundersökningar (Bovy & Rix 2013) och 21-cm-kartor. Samma inmatningsvektor med fem kvantiteter används.
3. Baryonisk sönderdelning – fem geometriska komponenter
Utifrån de fem observationerna delas den baryoniska massan upp i fem distinkta geometriska komponenter. Varje komponent har sin egen densitetsprofil och karakteristiska skala.
3.1 Total stjärn- och gasmassa
$$M_\star \;=\; 2\pi\,R_d^2\,\Sigma_d\,\Upsilon_\star$$$
$$M_\text{gas} \;=\; 1.33\,M_\text{HI} \qquad \text{(He-korrektion; Arnett 1996)}$$$
3.2 Komponenternas massa och skalor
| Komponent | Massa | Skala | Aktivering |
|---|---|---|---|
| Utbuktning | $M_b = 0,20\,M_\stjärna$ | $r_b = \max(0,5\,R_d,\,0,3\text{ kpc})$ | Om $T \leq 4$ |
| Tunn disk | M_\text{thin} = 0,75\,(M_\star – M_b)$ $. | $R_d$ | Alltid |
| Tjock disk | M_\text{thick} = 0.25\,(M_\star – M_b)$ $. | $1.5\,R_d$ | Alltid |
| Gasring | $M_\text{gas} = 1,33\,M_\text{HI}$ | $R_g = 1,7\,R_d$ (Broeils & Rhee 1997) | Alltid |
| Spiralformade armar | $M_\text{arm} = 0,10\,M_\text{tunn}$ (effektiv) | $R_d$ (följer tunn skiva) | Alltid |
3.3 Densitetsprofiler
Utbuktning (3D Hernquist)
$$\rho_b(r) \;=\; \frac{M_b\,r_b}{2\pi\,r\,(r + r_b)^3}}$$$
Tunna och tjocka stjärnskivor (2D exponentiell)
$$\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)}$$
Gasring (2D exponentiell med centralt hål)
$$\Sigma_\text{gas}(R) \;=\; \frac{M_\text{gas}}{2\pi\,R_g^2}\,\exp\!\left(-\frac{R_\text{hål}}{R} – \frac{R}{R_g}\right), \kvad R_\text{hål} = 0,5\,R_g$$$
Spiralarmsöverskott (2D, följer tunn skiva)
$$\Sigma_\text{arm}(R) \;=\; 0,10\;\Sigma_\text{tunn}(R)$$$
4. Vågkärnan
Varje baryoniskt masselement genererar ett BeeTheory-vågfält. Fältet vid en punkt $vec{r}$ som produceras av ett källelement vid $vec{r},’$ åtskilt av $D = |vec{r} – vec{r},’|$ styrs av kärnan av Yukawa-typ som härleds från den regulariserade vågfunktionen i not I:
BeeTheory vågkärna
$$\mathcal{K}_i(D) \;=\; K_0\,\frac{(1 + \alpha_i\,D)\,e^{-\alpha_i\,D}}{D^2}, \qquad \alpha_i \;=\; \frac{1}{\ell_i}$$$
Här är $K_0$ den universella vågmassamamplituden (ett enda dimensionslöst tal) och $\ell_i$ är koherenslängden för komponent $i$. Kärnan kodar ett kvasi-Newtonskt $1/D^2$ beteende vid korta separationer, modulerat av en exponentiell cutoff vid skalor bortom $\ell_i$. Formen $(1 + \alpha D)\,e^{-\alpha D}$ säkerställer kontinuitet och ändlig total innesluten massa vid oändlighet.
4.1 Komponenternas koherenslängder
Varje komponents koherenslängd bestäms av dess naturliga geometriska skala, multiplicerad med en dimensionslös konstant som är specifik för dess dimensionalitet:
| Komponent | Koherenslängd | Geometrisk konstant |
|---|---|---|
| Utbuktning (3D-sfär) | $\ell_b = c_\text{sph}\,r_b$ | $c_\text{sph}$$ |
| Tunn skiva (2D) | $\ell_\text{tunn} = c_\text{disk}\,R_d$ | $c_\text{disk}$ |
| Tjock skiva (2D) | $\ell_\text{thick} = c_\text{disk}\,(1.5\,R_d)$\ell_\text{thick} = c_\text{disk}\,(1.5\,R_d)$\ell_\text{thick} = c_\text{disk}\,(1.5\,R_d)$ | $c_\text{disk}$ |
| Gasring (2D) | $\ell_\text{gas} = c_\text{disk}\,R_g$$ | $c_\text{disk}$ |
| Spiralarmar (2D, azimutalt koncentrerade) | $\ell_\text{arm} = c_\text{arm}\,R_d$$ | $c_\text{arm}$$ |
De tre geometriska konstanterna $(c_\text{sph},\,c_\text{disk},\,c_\text{arm})$ är universella – de varierar inte från galax till galax. Tillsammans med den globala vågmassamamplituden $K_0$ och vågfältskopplingen $\lambda$ utgör de den fullständiga uppsättningen parametrar på teorinivå.
5. Vågfältskonvolution – integralekvationer per komponent
Vågfältsdensiteten vid en position $\vec{r}$ är en sammanfogning av den baryoniska källfördelningen med vågkärnan. För ett galaktiskt symmetriskt system (axialsymmetriskt, monopolär approximation) bidrar varje baryonisk komponent additivt:
Total vågfältsdensitet vid radien $r
$$\rho_\text{våg}(r) \;=\; \lambda \;\sum_{i \in \{\text{tunn, tjock, gas, arm, utbuktning}\}} \rho_\text{våg}^{(i)}(r)$$$
De fem integralerna är skrivna nedan, en per komponent. Varje integral omvandlar en baryonisk massfördelning till en vågfältsmassfördelning vid samma rumsliga punkt.
5.1 Bulge – integration av 3D-skal
$$\rho_\text{våg}^{(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’$$$
Integrationen sker över koncentriska sfäriska skal med radien $r’$. Fältpunkten vid radien $r$ från centrum ser varje skal vid en effektiv separation $D = \sqrt{r^2 + r’^2}$ i den monopolära approximationen. Integrationen sträcker sig ut till $r_\text{max} = 6\,r_b$, bortom vilket utbuktningens densitet är numeriskt försumbar.
5.2 Integration av tunn skiva och 2D-ring
$$\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’$$$
Skivan är uppdelad i koncentriska ringar med radien $R’$ och den infinitesimala bredden $dR’$, var och en med ytmassan $\Sigma_\text{thin}(R’)\,2\pi R’\,dR’$. Samma monopolära approximation gäller: vågfältet vid radien $r$ från centrum får bidrag från varje ring vid effektiv separation $D = \sqrt{r^2 + R’^2}$. Integrationsområdet är $R_\text{max} = 8\,R_d$.
5.3 Tjock skiva – 2D-ringintegration
$$\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’$$$
Identisk med integrationen för den tunna skivan, med $\Sigma_\text{thick}(R’)$ som källdensitet och en kernelparameter $\alpha_\text{thick} = 1/(c_\text{disk}\,\cdot 1,5\,R_d)$. Den tjocka diskens bredare radiella utbredning resulterar i ett något bredare vågkoherensområde.
5.4 Gasring – 2D-ringintegration med central utarmning
$$\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’$$
Gasfördelningen har ett centralt hål, som fångas upp av den radiella avgränsningen vid $R_\text{hole} = 0,5\,R_g$ i den nedre integrationsgränsen. Utanför denna avgränsning sträcker sig gasen längre än stjärndisken; detta återspeglas i den större karakteristiska skalan $R_g = 1,7\,R_d$, som matas in i koherenslängden $\ell_\text{gas} = c_\text{disk}\,R_g$.
5.5 Spiralarmsöverskott – 2D-ringintegration med reducerad 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’$$
Spiralarmarna behandlas som en axiellt medelvärdesbildad förstärkning av den tunna skivans yttäthet på nivån $10\%$, med sin egen koherenslängd $\ell_\text{arm} = c_\text{arm}\,R_d$. Kärnan är därför smalare än kärnan för den tunna skivan, vilket återspeglar den azimutala koncentrationen av spiralstrukturen.
6. Inkapslad vågmassa och beräknad rotationskurva
När den totala vågfältsdensiteten $\rho_\text{wave}(r)$ är känd, erhålls den inneslutna vågfältsmassan inom en sfär med radien $R$ genom radiell integration:
Massa för inneslutet vågfält
$$M_\text{våg}(R) \;=\; \int_0^{R} 4\pi\,r^2\,\rho_\text{våg}(r)\,dr$$
Den förutsagda cirkulära hastigheten vid radien $R$ följer sedan av den newtonska relationen, som kombinerar de baryoniska och vågfältbidragen i kvadratur:
Förutsedd cirkulär hastighet
$$V_c^2(R) \;=\; V_\text{bar}^2(R) \;+\; \frac{G\,M_\text{våg}(R)}{R}$$$
Den baryoniska hastigheten $V_\text{bar}(R)$ är i sig själv den kvadratiska summan av bidrag från de fyra skivliknande komponenterna (Freeman 1970 formel för varje exponentiell profil) och utbuktningen (Hernquist formel för sluten massa):
$$V_\text{bar}^2(R) \;=\; V_\text{bulge}^2 + V_\text{thin}^2 + V_\text{thick}^2 + V_\text{gas}^2$$$
där varje $V_i(R)$ är den newtonska standardcirkelhastigheten för motsvarande massfördelning.
7. Parametrar på teorinivå
Det kompletta BeeTheory-ramverket, som tillämpas på galaxer, innehåller fem parametrar på teorinivå. Dessa är universella: de varierar inte från galax till galax.
| Symbol | Betydelse | Roll |
|---|---|---|
| $K_0$ | Vågmassans amplitud | Ställer in den dimensionslösa skalan för vågkärnan |
| $c_\text{sph}$$ | Geometrisk konstant i 3D | Förhållandet $\ell/r_\text{scale}$ för sfäriska källor (bulge) |
| $c_\text{disk}$ | 2D geometrisk konstant | Förhållandet $\ell/R_\text{scale}$ för disk- och ringkällor |
| $c_\text{arm}$$ | Spiral geometrisk konstant | Förhållandet $\ell/R_d$ för det azimutalt koncentrerade armöverskottet |
| $\lambda$$ | Global vågfältskoppling | Skalar den totala vågfältsdensiteten |
Parametrarnas allmängiltighet
Alla fem parametrarna är globala. Samma numeriska värden gäller för Vintergatan, för irreguljära dvärggalaxer och för massiva spiraler. Den galaxspecifika informationen kommer endast in genom de fem observationsdata $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$. Modellen innehåller inte någon avstämbar parameter per galax.
8. Beräkningens enkelriktade karaktär
En öppen kedja – ingen återkoppling
Hela beräkningen går från ingångar till utgångar, i en riktning. De fotometriska och 21-cm-observationerna bestämmer den baryoniska nedbrytningen. Den baryoniska nedbrytningen bestämmer vågfältsdensiteten. Vågfältsdensiteten bestämmer den inneslutna vågmassan. Den inneslutna vågmassan bestämmer den förutspådda rotationskurvan. Rotationskurvan påverkar inte vid något tillfälle något tidigare steg i beräkningen.
Denna ensidighet har tre viktiga konsekvenser.
(a) När de fem parametrarna på teorinivå är fastställda är rotationskurvan en strikt förutsägelse, inte en anpassning. Jämförelsen med den observerade rotationskurvan är ett test, inte en kalibrering.
(b) Modellen har ingen mekanism för galax-för-galax-justering. Varje modifiering av prediktionen av rotationskurvan måste komma från en modifiering av ingångsvektorn $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$ eller från en ändring av de universella parametrarna på teorinivå $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$.
(c) Att kalibrera $lambda$ på en referensgalax är inte samma sak som att anpassa den till den galaxens rotationskurva. Kalibreringen fastställer ett enda globalt tal; rotationskurvan vid alla andra radier för referensgalaxen, och rotationskurvorna för alla andra galaxer, är sedan strikta förutsägelser av det kalibrerade ramverket.
9. Den centrala ytdensitetens roll (revidering av not XI)
Diagnosen i not XI identifierade att det kvarvarande prediktionsfelet korrelerar starkt med den centrala baryoniska yttätheten $\Sigma_d$, oberoende av diskens skalängd $R_d$. Formaliseringen som presenteras ovan är versionen av modellen innan detta resultat införlivas – den använder endast $R_d$ i uttrycken för koherenslängd $\ell_i = c_i\,R_d$.
Var förfiningen kommer in
I den förfinade modellen kommer koherenslängderna $\ell_i$ att bero på både $R_d$ och $\Sigma_d$, vilket ersätter det strikt linjära förhållandet $\ell_i = c_i\,R_d$ med en funktion $\ell_i = c_i\,R_d\,\phi(\Sigma_d/\Sigma_\text{ref})$ som absorberar den restpost som identifierades i not XI. Den funktionella formen för $\phi$ och dess parametrar kommer att fastställas i efterföljande noter, först på kalibreringsuppsättningen med 22 galaxer och sedan valideras genom blind prediktion på det återstående SPARC-urvalet.
Beräkningens enkelriktade struktur bevaras genom denna förfining: $\Sigma_d$ är en observationsinmatning, de modifierade koherenslängderna matas in i samma konvolutionsintegraler och rotationskurvan framträder som tidigare. Endast en operationell länk läggs till – beroendet av $\ell_i$ av en andra observabel.
10. Sammanfattning av metodiken
1. Ingångar. Fem observabler per galax: Hubble-typ $T$, diskskala $R_d$, ytans ljusstyrka $\Sigma_d$, HI-massa $M_\text{HI}$ och det universella förhållandet mellan stjärnornas massa och ljus $\Upsilon_\star$.
2. Baryonisk nedbrytning. Fem komponenter: bulge (om $T \leq 4$), tunn skiva, tjock skiva, gasring, spiralarmsöverskott. Var och en har en analytisk densitetsprofil.
3. Vågkärnan. Universell form av Yukawa-typ $\mathcal{K}_i(D) = K_0\,(1 + \alpha_i D)\,e^{-\alpha_i D}/D^2$ med koherenslängd $\ell_i = c_i\,R_\text{scale}$ bestämd av den geometriska omfattningen av varje komponent.
4. Konvolution. Varje komponent genererar en vågfältsdensitet via en endimensionell integral över ringar (2D-komponenter) eller skal (3D-utbuktning). Den totala vågfältsdensiteten är summan av de fem komponenterna, skalad med den globala kopplingen $\lambda$.
5. Utdata. Den inneslutna vågmassan $M_\text{wave}(R)$ integreras och kombineras med den baryoniska hastigheten $V_\text{bar}(R)$ för att ge den förväntade rotationskurvan $V_c(R)$.
6. Parametrar på teorinivå. $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$ – universell, ingen inställning per galax. En förfining som studeras kommer att lägga till ett beroende av $\Sigma_d$.
7. Riktning. Ingångar → baryoner → vågfält → rotationskurva. Ingen återkoppling. Rotationskurvan är en förutsägelse, inte en passform.
Referenser. 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. – The third law of galactic rotation, 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 – Vågbaserad kvantgravitation – Galaktisk metodik – © Technoplane S.A.S. 2026