BeeTheory – Fundamenter – Teknisk note XII
Formalisering:
Den galaktiske skala BeeTheory-beregning
Dette notat formaliserer BeeTheory-rammen, som den anvendes på en skivegalakse. Det specificerer de observerede input, den geometriske nedbrydning af den baryoniske fordeling, integralligningerne, der definerer bølgefeltet for hver komponent, og kæden af operationer, der giver den forudsagte rotationskurve. Proceduren er strengt ensrettet: Den observerede baryoniske struktur bestemmer bølgefeltet, som bestemmer rotationskurven – aldrig omvendt.
1. Beregningen i et diagram
En ensrettet kæde
Observeret fotometri $\;\longrightarrow\;$ Baryonisk nedbrydning $(\rho_\text{bar})$.
$\stor\nedadgående pil$
Bølgefeltkonvolution $\;\longrightarrow\;$ Bølgetæthed $(\rho_\text{wave}) $\big\down$
$\stor\nedadgående pil$
Masseintegration $\;\longrightarrow\;$ Indesluttet bølgemasse $(M_\text{wave})$.
$\big\nedadgående pil$
Newtonsk relation $\;\longrightarrow\;$ Forudsagt rotationskurve $(V_c)$.
Intet trin er inverteret. Rotationskurven $V_c(R)$ bruges aldrig som input.
2. Observationsmæssige input
For hver galakse kræver beregningen fem offentliggjorte observationer. Det er de eneste galaksespecifikke størrelser; alt andet beregnes ud fra dem. Der foretages ingen tilpasning til rotationskurven på dette tidspunkt.
| Symbol | Mængde | Kilde |
|---|---|---|
| $T$ | Hubbles morfologiske type | Katalog (de Vaucouleurs et al. 1991, SPARC) |
| $R_d$ | Stjerneskivens skalalængde (kpc) | Spitzer 3,6 µm fotometri (SPARC) |
| $\Sigma_d$ | Overfladens lysstyrke på den centrale skive ($L_\odot/\text{pc}^2$) | Spitzer 3,6 µm fotometri (SPARC) |
| $M_\tekst{HI}$ | Samlet atomar brintmasse ($M_odot$) | 21-cm radioobservationer (SPARC) |
| $\Upsilon_\stjerne$ | Stjernernes masse-til-lys-forhold ved 3,6 µm | Fast universel: $0.5\,M_\odot/L_\odot$ (McGaugh 2014) |
For Mælkevejen er $R_d$, $Sigma_d$ og $M_text{HI}$ erstattet af de tilsvarende værdier, der er bestemt ud fra interne stjerneundersøgelser (Bovy & Rix 2013) og 21-cm-kort. Den samme inputvektor med fem mængder bruges.
3. Baryonisk nedbrydning – fem geometriske komponenter
Ud fra de fem observationsinput opdeles den baryoniske masse i fem forskellige geometriske komponenter. Hver komponent har sin egen tæthedsprofil og karakteristiske skala.
3.1 Samlet stjerne- og gasmasse
$$M_\star \;=\; 2\pi\,R_d^2\,\Sigma_d\,\Upsilon_\star$$
$$M_\tekst{gas} \;=\; 1.33\,M_\text{HI} \qquad \text{(He-korrektion; Arnett 1996)}$$
3.2 Komponenternes masse og skala
| Komponent | Masse | Skala | Aktivering |
|---|---|---|---|
| Udbuling | $M_b = 0,20\,M_\star$ | $r_b = \max(0.5\,R_d,\,0.3\text{ kpc})$. | Hvis $T \leq 4$ |
| Tynd disk | $M_\text{thin} = 0.75\,(M_\star – M_b)$. | $R_d$ | Altid |
| Tyk disk | $M_\text{thick} = 0.25\,(M_\star – M_b)$. | $1.5\,R_d$ | Altid |
| Gasring | $M_\text{gas} = 1.33\,M_\text{HI}$. | $R_g = 1,7\,R_d$ (Broeils & Rhee 1997) | Altid |
| Spiralformede arme | $M_\text{arm} = 0.10\,M_\text{thin}$ (effektiv) | $R_d$ (følger en tynd disk) | Altid |
3.3 Tæthedsprofiler
Udbuling (3D Hernquist)
$$\rho_b(r) \;=\; \frac{M_b\,r_b}{2\pi\,r\,(r + r_b)^3}$$
Tynde og tykke stjerneskiver (2D-eksponentiel)
$$\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-eksponentiel med centralt hul)
$$\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$$
Spiralarmsoverskud (2D, følger tynd skive)
$$\Sigma_\text{arm}(R) \;=\; 0.10\;\Sigma_\text{thin}(R)$$
4. Bølgekernen
Hvert baryonisk masseelement genererer et BeeTheory-bølgefelt. Feltet i et punkt $vec{r}$ produceret af et kildeelement i $vec{r},’$ adskilt af $D = |vec{r} – vec{r},’|$ styres af Yukawa-kernen, der er afledt af den regulariserede bølgefunktion i note I:
BeeTheory-bølgekerne
$$\mathcal{K}_i(D) \;=\; K_0\,\frac{(1 + \alpha_i\,D)\,e^{-\alpha_i\,D}}{D^2}, \qquad \alpha_i \;=\; \frac{1}{\ell_i}$$
Her er $K_0$ den universelle bølgemasseamplitude (et enkelt dimensionsløst tal), og $\ell_i$ er kohærenslængden for komponent $i$. Kernen koder for en kvasi-newtonsk $1/D^2$-opførsel ved korte afstande, moduleret af en eksponentiel afskærmning ved skalaer over $\ell_i$. Formen $(1 + \alpha D)\,e^{-\alpha D}$ sikrer kontinuitet og endelig total indesluttet masse i det uendelige.
4.1 Komponenternes kohærenslængder
Kohærenslængden for hver komponent bestemmes af dens naturlige geometriske skala, ganget med en dimensionsløs konstant, der er specifik for dens dimensionalitet:
| Komponent | Kohærenslængde | Geometrisk konstant |
|---|---|---|
| Udbuling (3D-kugle) | $\ell_b = c_\text{sph}\,r_b$ | $c_\tekst{sph}$ |
| Tynd disk (2D) | $\ell_\tekst{tynd} = c_\tekst{disk}\,R_d$ | $c_\tekst{disk}$ |
| Tyk disk (2D) | $\ell_\text{thick} = c_\text{disk}\,(1.5\,R_d)$. | $c_\tekst{disk}$ |
| Gasring (2D) | $\ell_\tekst{gas} = c_\tekst{disk}\,R_g$$. | $c_\tekst{disk}$ |
| Spiralarme (2D, azimutalt koncentreret) | $\ell_\text{arm} = c_\text{arm}\,R_d$ | $c_\tekst{arm}$ |
De tre geometriske konstanter $(c_\text{sph},\,c_\text{disk},\,c_\text{arm})$ er universelle – de varierer ikke fra galakse til galakse. Sammen med den globale bølgemasseamplitude $K_0$ og bølgefeltkoblingen $\lambda$ udgør de det fulde sæt af parametre på teoriniveau.
5. Bølgefeltkonvolution – integralligninger pr. komponent
Bølgefeltets tæthed ved en position $\vec{r}$ er sammenfoldningen af den baryoniske kildefordeling med bølgekernen. For et galaktisk symmetrisk system (aksialsymmetrisk, monopolær tilnærmelse) bidrager hver baryonisk komponent additivt:
Samlet bølgefelt-tæthed ved radius $r$
$$\rho_\text{wave}(r) \;=\; \lambda \;\sum_{i \in \{\text{tynd, tyk, gas, arm, bule}\}} \rho_\text{wave}^{(i)}(r)$$
De fem integraler er skrevet nedenfor, et pr. komponent. Hvert integral omdanner en baryonisk massefordeling til en bølgefeltmassefordeling i det samme rumlige punkt.
5.1 Bulge – 3D-skalintegration
$$\rho_\tekst{bølge}^{(b)}(r) \;=\; \int_0^{r_\tekst{max}} \rho_b(r’)\;\mathcal{K}_b\!\left(\sqrt{r^2 + r’^2}\right)\;4\pi r’^2\,dr’$$$\rho_\text{wave}^{(b)}(r’)\;=\; \int_0^{r_text{max}}
Integrationen sker over koncentriske sfæriske skaller med radius $r’$. Feltpunktet ved radius $r$ fra centrum ser hver skal ved en effektiv adskillelse $D = \sqrt{r^2 + r’^2}$ i den monopolare tilnærmelse. Integrationen strækker sig ud til $r_\text{max} = 6\,r_b$, hvor udbulingstætheden er numerisk ubetydelig.
5.2 Tynd disk – 2D-ringintegration
$$\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’$$$\rho_text{wave}^{(\text{thin}}(r’)\;=\; \int_0^{R_text{thin}}(r’)
Skiven nedbrydes i koncentriske ringe med radius $R’$ og infinitesimal bredde $dR’$, som hver bærer overflademasse $\Sigma_\text{thin}(R’)\,2\pi R’\,dR’$. Den samme monopolære tilnærmelse gælder: Bølgefeltet ved radius $r$ fra centrum modtager bidrag fra hver ring ved effektiv adskillelse $D = \sqrt{r^2 + R’^2}$. Integrationsområdet er $R_\text{max} = 8\,R_d$.
5.3 Tyk disk – 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 af den tynde disk, med $\Sigma_\text{thick}(R’)$ som kildetæthed og en kerneparameter $\alpha_\text{thick} = 1/(c_\text{disk}\,\cdot 1.5\,R_d)$. Den bredere radiale udstrækning af den tykke disk resulterer i et lidt bredere bølgekohærensområde.
5.4 Gasring – 2D-ringintegration med central udtømning
$$\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’$$
Gasfordelingen har et centralt hul, der fanges af den radiale afskærmning ved $R_\text{hole} = 0,5\,R_g$ i den nedre integrationsgrænse. Uden for denne afgrænsning strækker gassen sig længere end stjerneskiven; dette afspejles i den større karakteristiske skala $R_g = 1,7\,R_d$, som indgår i kohærenslængden $\ell_\text{gas} = c_\text{disk}\,R_g$.
5.5 Spiralarmsoverskud – 2D-ringintegration med reduceret amplitude
$$\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’$$
Spiralarmene behandles som en aksialt gennemsnitlig forstærkning af den tynde skives overfladetæthed på niveauet $10\%$, med deres egen kohærenslængde $\ell_\text{arm} = c_\text{arm}\,R_d$. Kernen er derfor smallere end den tynde skives kerne, hvilket afspejler den azimutale koncentration af spiralstrukturen.
6. Indesluttet bølgemasse og forudsagt rotationskurve
Når den samlede bølgefeltstæthed $\rho_\text{wave}(r)$ er kendt, fås den indesluttede bølgefeltmasse inden for en kugle med radius $R$ ved radial integration:
Masse af lukket bølgefelt
$$M_\text{wave}(R) \;=\; \int_0^{R} 4\pi\,r^2\,\rho_\text{wave}(r)\,dr$$
Den forudsagte cirkulære hastighed ved radius $R$ følger derefter af den newtonske relation, der kombinerer de baryoniske og bølgefeltets bidrag i kvadratur:
Forudsagt cirkulær hastighed
$$V_c^2(R) \;=\; V_\text{bar}^2(R) \;+\; \frac{G\,M_\text{wave}(R)}{R}$$
Den baryoniske hastighed $V_\text{bar}(R)$ er i sig selv den kvadratiske sum af bidrag fra de fire disklignende komponenter (Freeman 1970-formel for hver eksponentiel profil) og udbulingen (Hernquist formel for indesluttet masse):
$$V_\text{bar}^2(R) \;=\; V_\text{bulge}^2 + V_\text{thin}^2 + V_\text{thick}^2 + V_\text{gas}^2$$
hvor hver $V_i(R)$ er den standard newtonske cirkulære hastighed for den tilsvarende massefordeling.
7. Parametre på teoriniveau
Den komplette BeeTheory-ramme, som anvendes på galakser, indeholder fem parametre på teoriniveau. Disse er universelle: De varierer ikke fra galakse til galakse.
| Symbol | Betydning | Rolle |
|---|---|---|
| $K_0$ | Bølgemassens amplitude | Indstiller den dimensionsløse skala for bølgekernen |
| $c_\tekst{sph}$ | Geometrisk 3D-konstant | Forholdet $\ell/r_\text{scale}$ for sfæriske kilder (bulge) |
| $c_\tekst{disk}$ | 2D geometrisk konstant | Forholdet $\ell/R_\text{scale}$ for disk- og ringkilder |
| $c_\tekst{arm}$ | Spiralens geometriske konstant | Forholdet $\ell/R_d$ for det azimutalt koncentrerede armoverskud |
| $\lambda$ | Global kobling af bølgefelt | Skalerer den samlede bølgefelt-tæthed |
Parametrenes universalitet
Alle fem parametre er globale. De samme numeriske værdier gælder for Mælkevejen, for irregulære dværggalakser og for massive spiraler. Den galaksespecifikke information kommer kun ind gennem de fem observationsinput $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$. Modellen indeholder ikke nogen parameter, der kan indstilles pr. galakse.
8. Beregningens ensrettede karakter
En åben kæde – ingen feedback
Hele beregningen flyder fra input til output, i én retning. De fotometriske og 21-cm-observationer bestemmer den baryoniske nedbrydning. Den baryoniske nedbrydning bestemmer bølgefeltets tæthed. Bølgefeltets tæthed bestemmer den indesluttede bølgemasse. Den indesluttede bølgemasse bestemmer den forudsagte rotationskurve. På intet tidspunkt påvirker rotationskurven noget tidligere trin i beregningen.
Denne ensrettethed har tre vigtige konsekvenser.
(a) Når de fem parametre på teoriniveau er fastsat, er rotationskurven en streng forudsigelse, ikke en tilpasning. Sammenligning med den observerede rotationskurve er en test, ikke en kalibrering.
(b) Modellen har ingen mekanisme til galakse-for-galakse-justering. Enhver ændring af forudsigelsen af rotationskurven skal komme fra en ændring af inputvektoren $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$ eller fra en ændring af de universelle parametre på teoriniveau $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$.
(c) At kalibrere $lambda$ på en referencegalakse er ikke det samme som at tilpasse den til galaksens rotationskurve. Kalibreringen bestemmer et enkelt globalt tal; rotationskurven ved alle andre radier af referencegalaksen og rotationskurverne for alle andre galakser er derefter strenge forudsigelser af den kalibrerede ramme.
9. Den centrale overfladetætheds rolle (Note XI revision)
Diagnosen i note XI identificerede, at den resterende forudsigelsesfejl korrelerer stærkt med den centrale baryoniske overfladetæthed $\Sigma_d$, uafhængigt af diskens skalalængde $R_d$. Formaliseringen præsenteret ovenfor er versionen af modellen, før denne opdagelse blev indarbejdet – den bruger kun $R_d$ i kohærenslængdeudtrykkene $\ell_i = c_i\,R_d$.
Hvor kommer forfiningen ind
I den raffinerede model vil kohærenslængderne $\ell_i$ afhænge af både $R_d$ og $\Sigma_d$, hvilket erstatter den strenge lineære relation $\ell_i = c_i\,R_d$ med en funktion $\ell_i = c_i\,R_d\,\phi(\Sigma_d/\Sigma_\text{ref})$, der absorberer det restbeløb, der blev identificeret i note XI. Den funktionelle form af $phi$ og dens parametre vil blive bestemt i efterfølgende noter, først på kalibreringssættet med 22 galakser og derefter valideret ved blind forudsigelse på den resterende SPARC-prøve.
Beregningens ensrettede struktur bevares ved denne forbedring: $\Sigma_d$ er et observationsinput, de modificerede kohærenslængder indgår i de samme sammenfoldningsintegraler, og rotationskurven fremkommer som før. Der er kun tilføjet et operationelt link – afhængigheden af $\ell_i$ af en anden observerbar størrelse.
10. Sammenfatning af metoden
1. Input. Fem observationer pr. galakse: Hubble-type $T$, diskskala $R_d$, overfladelysstyrke $\Sigma_d$, HI-masse $M_\text{HI}$ og det universelle forhold mellem stjernemasse og lys $\Upsilon_\star$.
2. Baryonisk nedbrydning. Fem komponenter: bulge (hvis $T \leq 4$), tynd skive, tyk skive, gasring, spiralarmsoverskud. Hver har en analytisk tæthedsprofil.
3. Bølgekerne. Universel form af Yukawa-typen $\mathcal{K}_i(D) = K_0\,(1 + \alpha_i D)\,e^{-\alpha_i D}/D^2$ med kohærenslængde $\ell_i = c_i\,R_\text{scale}$ bestemt af den geometriske udstrækning af hver komponent.
4. Konvolution. Hver komponent genererer en bølgefeltstæthed via et endimensionalt integral over ringe (2D-komponenter) eller skaller (3D-bule). Den samlede bølgefeltstæthed er summen af de fem komponenter, skaleret med den globale kobling $\lambda$.
5. Output. Den inkluderede bølgemasse $M_\text{wave}(R)$ integreres og kombineres med den baryoniske hastighed $V_\text{bar}(R)$ for at give den forudsagte rotationskurve $V_c(R)$.
6. Parametre på teoriniveau. $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$ – universel, ingen tuning pr. galakse. En forbedring under undersøgelse vil tilføje en afhængighed af $\Sigma_d$.
7. Retning. Inputs → baryoner → bølgefelt → rotationskurve. Ingen feedback. Rotationskurven er en forudsigelse, ikke en tilpasning.
Referencer. 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 – Bølgebaseret kvantegravitation – Galaktisk metodologi – © Technoplane S.A.S. 2026