BeeTheory – Grondslagen – Technische notitie XII
Formalisatie:
De bijentheoretische berekening op galactische schaal
Deze notitie formaliseert het BeeTheory-raamwerk zoals het wordt toegepast op een schijfstelsel. Het specificeert de waarnemingsinput, de geometrische decompositie van de baryonische distributie, de integraalvergelijkingen die het golfveld voor elke component definiëren, en de keten van operaties die de voorspelde rotatiecurve oplevert. De procedure is strikt eenrichtingsverkeer: de waargenomen baryonische structuur bepaalt het golfveld, dat de rotatiecurve bepaalt – nooit andersom.
1. De berekening in één diagram
Een eenrichtingsketting
Waargenomen fotometrie $longrightarrow;$ Baryonische ontbinding $(\rho_text{bar})$
$grote pijl omlaag
Golfveldconvolutie $;\longarrow;$ Golfdichtheid $(\rho_text{wave})$
$grote pijl omlaag
Massaintegratie $\;\longarrow;$ingesloten golfmassa $(M_text{wave})$
$grote pijl omlaag
Voorspelde rotatiecurve $(V_c)$.
Er wordt geen stap omgekeerd. De rotatiecurve $V_c(R)$ wordt nooit als invoer gebruikt.
2. Waarnemingsinput
Voor elk sterrenstelsel heeft de berekening vijf gepubliceerde waarneembaarheden nodig. Dit zijn de enige melkwegspecifieke grootheden; al de rest wordt daaruit berekend. Er wordt in dit stadium niet tegen de rotatiecurve gepast.
| Symbool | Hoeveelheid | Bron |
|---|---|---|
| $T$ | Hubble morfologisch type | Catalogus (de Vaucouleurs et al. 1991, SPARC) |
| $R_d$ | Schaallengte stellaire schijf (kpc) | Spitzer 3,6 µm fotometrie (SPARC) |
| $Sigma_d$ | Helderheid oppervlak centrale schijf ($L_\odot/{pc}^2$) | Spitzer 3,6 µm fotometrie (SPARC) |
| $M_text{HI}$ | Totale atomaire waterstofmassa ($M_odot$) | 21-cm radiowaarnemingen (SPARC) |
| $Upsilon_ster | Verhouding stellaire massa/licht bij 3,6 µm | Vast universeel: $0.5,M_\odot/L_\odot$ (McGaugh 2014) |
Voor de Melkweg zijn $R_d$, $Sigma_d$ en $M_text{HI}$ vervangen door de analoge waarden die bepaald zijn uit interne stellaire onderzoeken (Bovy & Rix 2013) en 21-cm-kaarten. Dezelfde invoervector met vijf grootheden wordt gebruikt.
3. Baryonische ontbinding – vijf geometrische componenten
Op basis van de vijf waarnemingsgegevens wordt de baryonische massa verdeeld in vijf verschillende geometrische componenten. Elke component heeft zijn eigen dichtheidsprofiel en karakteristieke schaal.
3.1 Totale stellaire en gasmassa’s
$$$M___sterren_$; 2_pi,R_d^2,\Sigma_d,\Upsilon__sterren$$
$$M_text{gas} \1,33, M_text{HI} \quad {(He correctie; Arnett 1996)}$
3.2 Componentenmassa’s en -schalen
| Component | Massa | Schaal | Activering |
|---|---|---|---|
| Bulge | $M_b = 0,20,M_star | $r_b = \max(0,5,R_d,\,0,3{ kpc})$ | Als $T ≤ 4$ |
| Dunne schijf | $M_text{dun} = 0.75,(M_star – M_b)$ | $R_d$ | Altijd |
| Dikke schijf | $M_text{dikte} = 0.25,(M_star – M_b)$ | $1,5,R_d$ | Altijd |
| Gasring | $M_text{gas} = 1,33,M_text{HI}$ | $R_g = 1.7,R_d$ (Broeils & Rhee 1997) | Altijd |
| Spiraalarmen | $M_text{arm} = 0,10,M_text{thin}$ (effectief) | $R_d$ (volgt dunne schijf) | Altijd |
3.3 Dichtheidsprofielen
Uitstulping (3D Hernquist)
$$\rho_b(r) \;=; frac{M_b,r_b}{2\pi,r,(r + r_b)^3}$
Dunne en dikke sterrenschijven (2D exponentieel)
$$\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 exponentieel met centraal gat)
$$Sigma_text{gas}(R) \;=; \frac{M_text{gas}{2{M_text{gas}{2{M_text{gas}}{2{M_text{gas}{2{M_text{gas}}{2{M_text{gas}{2}{R_g^2}{R_g}{R_g}{R_g}{R_g}{R_g}{R_g}{R}{R_g}{R}{R_g}{R}{R}{R_g}{R}{R}{R_g}{R}{R}{R_g}{R}{R}{R}{R}{R}{R}{R}{R}{R}=0.5}$$
Spiraalarmovermaat (2D, volgt dunne schijf)
$$ \Sigma_text{arm}(R) \;=; 0.10;\Sigma_text{dun}(R)$$
4. De golfkernel
Elk baryonisch massa-element genereert een BeeTheory-golfveld. Het veld in een punt $vec{r}$ geproduceerd door een bronelement op $vec{r},’$ gescheiden door $D = |vec{r} – vec{r},’|$ wordt beheerst door de Yukawa-type kernel afgeleid van de geregulariseerde golffunctie van Noot I:
Bijentheorie golfkernel
K_0,frac{(1 + \alpha_i,D)^{- \alpha_i,D}{D^2}, \kwadraat \alpha_i $;=; \frac{1}{\ell_i}}$.
Hier is $K_0$ de universele golf-massa-amplitude (een enkel dimensieloos getal) en $\ell_i$ de coherentielengte van component $i$. De kernel codeert een quasi-Newtoniaans $1/D^2$ gedrag op korte afstanden, gemoduleerd door een exponentiële cutoff op schalen voorbij $\ell_i$. De vorm $(1 + \alpha D)\,e^{-\alpha D}$ zorgt voor continuïteit en eindige totale ingesloten massa op oneindig.
4.1 Component coherentielengtes
De coherentielengte van elke component wordt bepaald door zijn natuurlijke geometrische schaal, vermenigvuldigd met een dimensieloze constante die specifiek is voor zijn dimensionaliteit:
| Component | Coherentielengte | Geometrische constante |
|---|---|---|
| Uitstulping (3D-bol) | $r_b = c_tekst{sph}r_b$ | $c_text{sph}$ |
| Dunne schijf (2D) | $ cell_text{thin} = c_text{disk}},R_d$ | $c_text{disk}$ |
| Dikke schijf (2D) | $ell_tekst{dikte} = c_tekst{schijf},(1.5},R_d)$ | $c_text{disk}$ |
| Gasring (2D) | $r_tekst{gas} = c_tekst{schijf}r,R_g$ | $c_text{disk}$ |
| Spiraalarmen (2D, azimuthaal geconcentreerd) | $_tekst{arm} = c_tekst{arm}, R_d$ | $c_text{arm}$ |
De drie geometrische constanten $(c_text{sph},¦,c_text{disk},¦,c_text{arm})$ zijn universeel – ze variëren niet van melkwegstelsel tot melkwegstelsel. Samen met de globale golfmassa-amplitude $K_0$ en de golfveldkoppeling $\lambda$ vormen ze de volledige set parameters op theorieniveau.
5. Golfveldconvolutie – integraalvergelijkingen per component
De golfvelddichtheid op een positie $vec{r}$ is de convolutie van de baryonische bronverdeling met de golfkernel. Voor een galactisch symmetrisch systeem (axiaal symmetrisch, monopolaire benadering) draagt elke baryonische component additief bij:
Totale golfvelddichtheid bij een straal van $r$
$$rho_text{wave}(r) \;=; \lambda \;sum_{i \in \{{dun, dik, gas, arm, uitstulping}}}} \rho_text{wave}^{(i)}(r)$$
De vijf integralen staan hieronder, één per component. Elke integraal zet een baryonische massaverdeling om in een golfveldmassaverdeling op hetzelfde ruimtelijke punt.
5.1 Bulge – 3D schelpintegratie
$$ \rho_{wave}^{(b)}(r) \;=; \int_0^{r_{text{max}} \rho_b(r’)\;\mathcal{K}_b!\links(\sqrt{r^2 + r’^2})\;4\pi r’^2,dr’$$
De integratie gaat over concentrische sferische schillen met straal $r’$. Het veldpunt op straal $r$ van het centrum ziet elke schil op een effectieve afstand $D = \sqrt{r^2 + r’^2}$ in de monopolaire benadering. De integratie strekt zich uit tot $r_text{max} = 6\,r_b$, waarboven de uitstulpingsdichtheid numeriek verwaarloosbaar is.
5.2 Dunne schijf – 2D ringintegratie
$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_{text{thin}(R’)º;\mathcal{K}_{text{thin}{\link(\sqrt{r^2 + R’^2}}};\2pi R’º,dR’$$
De schijf wordt ontleed in concentrische ringen met straal $R’$ en infinitesimale breedte $dR’$, elk met oppervlaktemassa $\Sigma_text{thin}(R’)\pi R’\,dR’$. Dezelfde monopolaire benadering is van toepassing: het golfveld op straal $r$ van het centrum ontvangt bijdragen van elke ring op effectieve afstand $D = \sqrt{r^2 + R’^2}$. Het integratiebereik is $R_text{max} = 8\,R_d$.
5.3 Dikke schijf – 2D ringintegratie
$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_{text{dikte}(R’);\mathcal{K}_{text{dikte}!\links(\sqrt{r^2 + R’^2})\rechts);2\pi R’,dR’$$
Identiek aan de integratie van de dunne schijf, met $\Sigma_{dikte}(R’)$ als de brondichtheid en een kernelparameter $\alpha_{dikte} = 1/(c_\text{disk},\cdot 1.5,R_d)$. De grotere radiale omvang van de dikke schijf resulteert in een iets breder golfcoherentiebereik.
5.4 Gasring – 2D ringintegratie met centrale uitputting
$$\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’$$
De gasdistributie heeft een centraal gat, gevangen door de radiale cutoff op $R_text{hole} = 0,5,R_g$ in de ondergrens van de integratie. Buiten deze cutoff strekt het gas zich verder uit dan de stellaire schijf; dit wordt weerspiegeld in de grotere karakteristieke schaal $R_g = 1,7 R_d$, die wordt ingevoerd in de coherentielengte $\_text{gas} = c_text{disk},R_g$.
5.5 Spiraalarmovermaat – 2D ringintegratie met gereduceerde 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’$$
De spiraalarmen worden behandeld als een axiaal gemiddelde versterking van de oppervlaktedichtheid van de dunne schijf op het niveau van $10%$, met hun eigen coherentielengte $\ell_text{arm} = c_text{arm},R_d$. De kernel is daarom smaller dan de kernel van de dunne schijf, wat de azimutale concentratie van de spiraalstructuur weerspiegelt.
6. Ingesloten golfmassa en voorspelde rotatiecurve
Zodra de totale golfvelddichtheid $\rho_text{wave}(r)$ bekend is, wordt de ingesloten golfveldmassa binnen een bol met straal $R$ verkregen door radiale integratie:
Ingesloten golfveldmassa
$$M_text{wave}(R) \;=; \int_0^{R} 4\pi,r^2,\rho_text{wave}(r)\,dr$$
De voorspelde cirkelsnelheid bij een straal van $R$ volgt dan uit de Newtoniaanse relatie, waarbij de baryonische en golfveldbijdragen in kwadratuur worden gecombineerd:
Voorspelde cirkelsnelheid
$$V_c^2(R) ^;=; V_\text{bar}^2(R) ^;+; \frac{G,M_\text{wave}(R)}{R}$
De baryonische snelheid $V_\text{bar}(R)$ is zelf de kwadratische som van bijdragen van de vier schijfachtige componenten (Freeman 1970-formule voor elk exponentieel profiel) en de uitstulping (Hernquist-formule voor ingesloten massa):
$$V_\text{bar}^2(R) \;=; V_\text{bulge}^2 + V_\text{thin}^2 + V_\text{thick}^2 + V_\text{gas}^2$
waarbij elke $V_i(R)$ de standaard Newtoniaanse cirkelsnelheid van de corresponderende massaverdeling is.
7. Parameters op theoretisch niveau
Het volledige BeeTheory raamwerk, zoals toegepast op melkwegstelsels, bevat vijf parameters op theorieniveau. Deze zijn universeel: ze variëren niet van melkwegstelsel tot melkwegstelsel.
| Symbool | Betekenis | Rol |
|---|---|---|
| $K_0$ | Golfmassa-amplitude | Stelt de dimensieloze schaal van de golfkern in |
| $c_text{sph}$ | 3D geometrische constante | Verhouding $\ell/r_text{scale}$ voor sferische bronnen (uitstulping) |
| $c_text{disk}$ | 2D geometrische constante | Verhouding $R_text{scale}$ voor schijf- en ringbronnen |
| $c_text{arm}$ | Geometrische constante van spiraal | Verhouding $\ell/R_d$ voor het azimutaal geconcentreerde armoverschot |
| $lambda$ | Globale golfveldkoppeling | Schaalt de totale golfvelddichtheid |
Universaliteit van de parameters
Alle vijf parameters zijn globaal. Dezelfde numerieke waarden gelden voor de Melkweg, voor dwergonregelmatige stelsels en voor massieve spiralen. De melkwegspecifieke informatie komt alleen binnen via de vijf waarnemingsingangen $(T,\,R_d,\Sigma_d,\,M_text{HI},\,\Upsilon_star)$. Het model bevat geen per melkwegstelsel instelbare parameter.
8. Het eenrichtingskarakter van de berekening
Een open keten – geen terugkoppeling
De hele berekening loopt van input naar output, in één richting. De fotometrische en 21-cm waarnemingen bepalen de baryonische decompositie. De baryonische decompositie bepaalt de golfvelddichtheid. De golfvelddichtheid bepaalt de ingesloten golfmassa. De ingesloten golfmassa bepaalt de voorspelde rotatiecurve. Op geen enkel moment beïnvloedt de rotatiecurve een eerdere stap van de berekening.
Deze unidirectionaliteit heeft drie belangrijke gevolgen.
(a) Zodra de vijf parameters op theorieniveau vastliggen, is de rotatiecurve een strikte voorspelling, geen fit. Vergelijking met de waargenomen rotatiecurve is een test, geen kalibratie.
(b) Het model heeft geen mechanisme voor aanpassing per melkwegstelsel. Elke verandering van de voorspelling van de rotatiecurve moet komen van een verandering van de ingangsvector $(T,\,R_d,\Sigma_d,\,M_text{HI},\,\Upsilon_ster)$ of van een verandering van de universele parameters op theoretisch niveau $(K_0,\,c_text{sph},\,c_text{schijf},\,c_text{arm},\lambda)$.
(c) Het ijken van $lambda$ op een referentiestelsel is niet hetzelfde als het passen op de rotatiecurve van dat stelsel. De ijking bepaalt een enkel globaal getal; de rotatiecurve bij alle andere stralen van het referentiestelsel, en de rotatiecurves van alle andere melkwegstelsels, zijn dan strikte voorspellingen van het gekalibreerde raamwerk.
9. De rol van centrale oppervlaktedichtheid (Nota XI herziening)
De diagnostiek van Noot XI stelde vast dat de resterende voorspellingsfout sterk correleert met de centrale baryonische oppervlaktedichtheid $\Sigma_d$, onafhankelijk van de schijfschaallengte $R_d$. De hierboven gepresenteerde formalisatie is de versie van het model voordat deze bevinding is opgenomen – het gebruikt alleen $R_d$ in de coherentielengte-uitdrukkingen $\ell_i = c_i,R_d$.
Waar de verfijning binnenkomt
In het verfijnde model hangen de coherentielengtes $\ell_i$ af van zowel $R_d$ als $\Sigma_d$, waarbij de strikte lineaire relatie $\ell_i = c_i,R_d$ vervangen wordt door een functie $\ell_i = c_i,R_d,\phi(\Sigma_d/\Sigma_text{ref})$ die het residu absorbeert dat in opmerking XI geïdentificeerd is. De functionele vorm van $phi$ en zijn parameters zullen in volgende aantekeningen bepaald worden, eerst op de kalibratieset van 22 melkwegstelsels, daarna gevalideerd door blinde voorspelling op de resterende SPARC-monsters.
De eenrichtingsstructuur van de berekening blijft behouden door deze verfijning: $\Sigma_d$ is een waarnemingsinput, de gewijzigde coherentielengtes voeden dezelfde convolutie-integralen, en de rotatiecurve ontstaat als voorheen. Er wordt slechts één operationele schakel toegevoegd – de afhankelijkheid van $\ell_i$ van een tweede waarneembaar gegeven.
10. Samenvatting van de methodologie
1. Invoeren. Vijf observabelen per melkwegstelsel: Hubble-type $T$, schijfschaal $R_d$, oppervlaktehelderheid $\Sigma_d$, HI-massa $M_text{HI}$ en de universele stellaire massa/lichtverhouding $Upsilon_star$.
2. Baryonische ontbinding. Vijf componenten: uitstulping (als $T ≤ 4$), dunne schijf, dikke schijf, gasring, spiraalarmoverschot. Elk heeft een analytisch dichtheidsprofiel.
3. Golfkernel. Universele Yukawa-type vorm $\mathcal{K}_i(D) = K_0,(1 + \alpha_i D)\,e^{-\alpha_i D}/D^2$ met coherentielengte $\ell_i = c_i,R_text{scale}$ bepaald door de geometrische omvang van elke component.
4. Convolutie. Elke component genereert een golfvelddichtheid via een eendimensionale integraal over ringen (2D componenten) of schelpen (3D uitstulping). De totale golfvelddichtheid is de som van de vijf componenten, geschaald door de globale koppeling $lambda$.
5. Output. De ingesloten golfmassa $M_\text{wave}(R)$ wordt geïntegreerd en gecombineerd met de baryonische snelheid $V_\text{bar}(R)$ om de voorspelde rotatiecurve $V_c(R)$ te verkrijgen.
6. Parameters op theoretisch niveau. $(K_0,\,c_text{sph},\,c_text{disk},\,c_text{arm},\,\lambda)$ – universeel, geen afstemming per melkwegstelsel. Een bestudeerde verfijning zal een afhankelijkheid van $Sigma_d$ toevoegen.
7. Richting. Ingangen → baryonen → golfveld → rotatiecurve. Geen terugkoppeling. De rotatiecurve is een voorspelling, geen fit.
Referenties. 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. – Een analytisch model voor bolvormige sterrenstelsels en bulges, ApJ 356, 359 (1990). – Broeils, A. H., Rhee, M.-H. – Korte 21-cm WSRT-waarnemingen van spiraalstelsels en onregelmatige sterrenstelsels, A&A 324, 877 (1997). – McGaugh, S. S. – De derde wet van galactische rotatie, Galaxies 2, 601 (2014). – Bovy, J., Rix, H.-W. – Een directe dynamische meting van het dichtheidsprofiel van het schijfoppervlak van de Melkweg, 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 – Op golven gebaseerde kwantumzwaartekracht – Galactische methodologie – © Technoplane S.A.S. 2026