BeeTheory – Grundlagen – Technischer Hinweis XII
Formalisierung:
Die Berechnung der Bienentheorie auf galaktischer Skala
Diese Notiz formalisiert den Rahmen der BeeTheory, wie er auf eine Scheibengalaxie angewendet wird. Sie spezifiziert die Beobachtungsdaten, die geometrische Zerlegung der baryonischen Verteilung, die Integralgleichungen, die das Wellenfeld für jede Komponente definieren, und die Kette von Operationen, die die vorhergesagte Rotationskurve ergibt. Das Verfahren ist streng unidirektional: Die beobachtete baryonische Struktur bestimmt das Wellenfeld, das wiederum die Rotationskurve bestimmt – niemals umgekehrt.
1. Die Berechnung in einem Diagramm
Eine unidirektionale Kette
Beobachtete Photometrie $\;\longrightarrow\;$ Baryonenzerlegung $(\rho_\text{bar})$
$\big\downarrow$
Wellenfeldfaltung $\;\longrightarrow\;$ Wellendichte $(\rho_\text{wave})$
$\big\downarrow$
Massenintegration $\;\longrightarrow\;$ Eingeschlossene Wellenmasse $(M_\text{wave})$
$\big\downarrow$
Newtonsche Beziehung $\;\longrightarrow\;$ Vorhergesagte Rotationskurve $(V_c)$
Es wird kein Schritt invertiert. Die Rotationskurve $V_c(R)$ wird nie als Eingabe verwendet.
2. Eingaben durch Beobachtung
Für jede Galaxie erfordert die Berechnung fünf veröffentlichte Beobachtungsgrößen. Dies sind die einzigen galaxienspezifischen Größen; alles andere wird aus ihnen errechnet. In diesem Stadium wird keine Anpassung an die Rotationskurve vorgenommen.
| Symbol | Menge | Quelle |
|---|---|---|
| $T$ | Hubble morphologischer Typ | Katalog (de Vaucouleurs et al. 1991, SPARC) |
| $R_d$ | Skalenlänge der stellaren Scheibe (kpc) | Spitzer 3,6 µm Photometrie (SPARC) |
| $\Sigma_d$ | Helligkeit der zentralen Scheibenoberfläche ($L_\odot/\text{pc}^2$) | Spitzer 3,6 µm Photometrie (SPARC) |
| $M_\text{HI}$ | Gesamte atomare Wasserstoffmasse ($M_odot$) | 21-cm-Radiobeobachtungen (SPARC) |
| $\Upsilon_\star$ | Stellares Masse-zu-Licht-Verhältnis bei 3,6 µm | Festes Universal: $0,5\,M_\odot/L_\odot$ (McGaugh 2014) |
Für die Milchstraße werden $R_d$, $Sigma_d$ und $M_text{HI}$ durch die entsprechenden Werte ersetzt, die aus internen stellaren Durchmusterungen (Bovy & Rix 2013) und 21-cm-Karten ermittelt wurden. Es wird derselbe Fünf-Quantitäten-Eingangsvektor verwendet.
3. Baryonische Zersetzung – fünf geometrische Komponenten
Aus den fünf Beobachtungsdaten wird die baryonische Masse in fünf verschiedene geometrische Komponenten aufgeteilt. Jede Komponente hat ihr eigenes Dichteprofil und ihre eigene charakteristische Skala.
3.1 Gesamte Stern- und Gasmassen
$$M_\star \;=\; 2\pi\,R_d^2\,\Sigma_d\,\Upsilon_\star$$
$$M_\text{gas} \;=\; 1.33\,M_\text{HI} \qquad \text{(He-Korrektur; Arnett 1996)}$$
3.2 Massen und Maßstäbe der Komponenten
| Komponente | Masse | Skala | Aktivierung |
|---|---|---|---|
| Wulst | $M_b = 0,20\,M_\star$ | $r_b = \max(0.5\,R_d,\,0.3\text{ kpc})$ | Wenn $T \leq 4$ |
| Dünne Scheibe | $M_\text{thin} = 0.75\,(M_\star – M_b)$ | $R_d$ | Immer |
| Dicke Scheibe | $M_\text{dick} = 0.25\,(M_\star – M_b)$ | $1.5\,R_d$ | Immer |
| Gasring | $M_\text{gas} = 1.33\,M_\text{HI}$ | $R_g = 1,7\,R_d$ (Broeils & Rhee 1997) | Immer |
| Spiralförmige Arme | $M_\text{arm} = 0.10\,M_\text{dünn}$ (effektiv) | $R_d$ (folgt dünne Scheibe) | Immer |
3.3 Dichteprofile
Bulge (3D Hernquist)
$$\rho_b(r) \;=\; \frac{M_b\,r_b}{2\pi\,r\,(r + r_b)^3}$$
Dünne und dicke stellare Scheiben (2D-Exponential)
$$\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-Exponential mit zentralem Loch)
$$$Sigma_\text{gas}(R) \;=\; \frac{M_\text{gas}}{2\pi\,R_g^2}\,\exp\!\links(-\frac{R_\text{hole}}{R} – \frac{R}{R_g}\rechts), \quad R_\text{hole} = 0.5\,R_g$$
Spiralarm-Überschuss (2D, folgt dünner Scheibe)
$$\Sigma_\text{arm}(R) \;=\; 0.10\;\Sigma_\text{dünn}(R)$$
4. Der Wellenkern
Jedes baryonische Massenelement erzeugt ein BeeTheory-Wellenfeld. Das Feld an einem Punkt $vec{r}$, das von einem Quellenelement an $vec{r},’$ erzeugt wird, das durch $D = |vec{r} – vec{r},’|$ getrennt ist, unterliegt dem Yukawa-Kernel, der aus der regularisierten Wellenfunktion von Anmerkung I abgeleitet ist:
Bienentheorie-Wellenkern
$$\mathcal{K}_i(D) \;=\; K_0\,\frac{(1 + \alpha_i\,D)\,e^{-\alpha_i\,D}}{D^2}, \qquad \alpha_i \;=\; \frac{1}{\ell_i}$$
Dabei ist $K_0$ die universelle Wellen-Masse-Amplitude (eine einzelne dimensionslose Zahl) und $\ell_i$ die Kohärenzlänge der Komponente $i$. Der Kernel kodiert ein quasi-newtonsches $1/D^2$ Verhalten bei kurzen Abständen, das bei Skalen jenseits von $\ell_i$ durch einen exponentiellen Cutoff moduliert wird. Die Form $(1 + \alpha D)\,e^{-\alpha D}$ gewährleistet Kontinuität und eine endliche eingeschlossene Gesamtmasse bei Unendlichkeit.
4.1 Kohärenzlängen der Komponenten
Die Kohärenzlänge jeder Komponente wird durch ihren natürlichen geometrischen Maßstab bestimmt, multipliziert mit einer dimensionslosen Konstante, die für ihre Dimensionalität spezifisch ist:
| Komponente | Länge der Kohärenz | Geometrische Konstante |
|---|---|---|
| Bulge (3D-Kugel) | $\ell_b = c_\text{sph}\,r_b$ | $c_\text{sph}$ |
| Dünne Scheibe (2D) | $\ell_\text{thin} = c_\text{disk}\,R_d$ | $c_\text{disk}$ |
| Dicke Scheibe (2D) | $\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}$ |
| Spiralförmige Arme (2D, azimutal konzentriert) | $\ell_\text{arm} = c_\text{arm}\,R_d$ | $c_\text{arm}$ |
Die drei geometrischen Konstanten $(c_\text{sph},\,c_\text{disk},\,c_\text{arm})$ sind universell – sie variieren nicht von Galaxie zu Galaxie. Zusammen mit der globalen Wellenmassenamplitude $K_0$ und der Wellenfeldkopplung $\lambda$ bilden sie den vollständigen Satz von Parametern auf Theorieebene.
5. Wellenfeldfaltung – Integralgleichungen pro Komponente
Die Wellenfelddichte an einer Position $\vec{r}$ ist die Faltung der baryonischen Quellenverteilung mit dem Wellenkern. Für ein galaktisch symmetrisches System (achsensymmetrisch, monopolare Näherung) trägt jede baryonische Komponente additiv bei:
Gesamtwellenfelddichte am Radius $r$
$$$rho_\text{wave}(r) \;=\; \lambda \;\sum_{i \in \{\text{dünn, dick, Gas, Arm, Wulst}\} \rho_\text{wave}^{(i)}(r)$$
Die fünf Integrale sind unten aufgeführt, eines pro Komponente. Jedes Integral wandelt eine baryonische Massenverteilung in eine Wellenfeld-Massenverteilung an demselben räumlichen Punkt um.
5.1 Bulge – Integration der 3D-Schale
$$$\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’$$
Die Integration erfolgt über konzentrische Kugelschalen mit dem Radius $r’$. Der Feldpunkt am Radius $r$ vom Zentrum sieht jede Schale in einem effektiven Abstand $D = \sqrt{r^2 + r’^2}$ in der monopolaren Approximation. Die Integration erstreckt sich bis $r_\text{max} = 6\,r_b$, jenseits dessen die Bulge-Dichte numerisch vernachlässigbar ist.
5.2 Dünne Scheibe – 2D-Ring-Integration
$$\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’$$
Die Scheibe wird in konzentrische Ringe mit Radius $R’$ und infinitesimaler Breite $dR’$ zerlegt, von denen jeder eine Oberflächenmasse $\Sigma_\text{thin}(R‘)\,2\pi R’\,dR’$ trägt. Es gilt die gleiche monopolare Näherung: das Wellenfeld am Radius $r$ vom Zentrum erhält Beiträge von jedem Ring im effektiven Abstand $D = \sqrt{r^2 + R’^2}$. Der Integrationsbereich ist $R_\text{max} = 8\,R_d$.
5.3 Dicke Scheibe – 2D-Ring-Integration
$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{dick}(R‘)\;\mathcal{K}_\text{dick}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
Identisch mit der Integration der dünnen Scheibe, mit $\Sigma_\text{thick}(R‘)$ als Quelldichte und einem Kernelparameter $\alpha_\text{thick} = 1/(c_\text{disk}\,\cdot 1.5\,R_d)$. Die größere radiale Ausdehnung der dicken Scheibe führt zu einem etwas breiteren Wellenkohärenzbereich.
5.4 Gasring – 2D-Ringintegration mit zentraler Verarmung
$$\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’$$
Die Gasverteilung hat ein zentrales Loch, das durch den radialen Cutoff bei $R_\text{hole} = 0.5\,R_g$ in der unteren Grenze der Integration erfasst wird. Außerhalb dieser Grenze erstreckt sich das Gas weiter als die stellare Scheibe; dies spiegelt sich in der größeren charakteristischen Skala $R_g = 1,7\,R_d$ wider, die in die Kohärenzlänge $\ell_\text{gas} = c_\text{disk}\,R_g$ eingeht.
5.5 Spiralarm-Exzess – 2D-Ring-Integration mit reduzierter 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’$$
Die Spiralarme werden als axial gemittelte Verstärkung der Oberflächendichte der dünnen Scheibe auf dem Niveau von $10\%$ behandelt, mit ihrer eigenen Kohärenzlänge $\ell_\text{arm} = c_\text{arm}\,R_d$. Der Kernel ist also schmaler als der Kernel der dünnen Scheibe und spiegelt die azimutale Konzentration der Spiralstruktur wider.
6. Eingeschlossene Wellenmasse und vorhergesagte Rotationskurve
Sobald die gesamte Wellenfelddichte $\rho_\text{wave}(r)$ bekannt ist, erhält man durch radiale Integration die eingeschlossene Wellenfeldmasse innerhalb einer Kugel mit Radius $R$:
Eingeschlossene Wellenfeldmasse
$$M_\text{Welle}(R) \;=\; \int_0^{R} 4\pi\,r^2\,\rho_\text{Welle}(r)\,dr$$
Die vorhergesagte kreisförmige Geschwindigkeit bei Radius $R$ ergibt sich dann aus der Newtonschen Beziehung, die die baryonischen und Wellenfeldbeiträge in Quadratur kombiniert:
Vorausgesagte Kreisgeschwindigkeit
$$V_c^2(R) \;=\; V_\text{bar}^2(R) \;+\; \frac{G\,M_\text{wave}(R)}{R}$$
Die baryonische Geschwindigkeit $V_\text{bar}(R)$ ist ihrerseits die quadratische Summe der Beiträge der vier scheibenförmigen Komponenten (Freeman-Formel von 1970 für jedes Exponentialprofil) und des Bulge (Hernquist-Formel für eingeschlossene Masse):
$$V_\text{bar}^2(R) \;=\; V_\text{bulge}^2 + V_\text{thin}^2 + V_\text{thick}^2 + V_\text{gas}^2$$
wobei jedes $V_i(R)$ die Newtonsche Standard-Kreisgeschwindigkeit der entsprechenden Massenverteilung ist.
7. Parameter auf Theorieebene
Der vollständige BeeTheory-Rahmen, wie er auf Galaxien angewandt wird, enthält fünf Parameter auf der Theorieebene. Diese sind universell: Sie variieren nicht von Galaxie zu Galaxie.
| Symbol | Bedeutung | Rolle |
|---|---|---|
| $K_0$ | Amplitude der Wellenmasse | Legt die dimensionslose Skala des Wellenkerns fest |
| $c_\text{sph}$ | Geometrische 3D-Konstante | Verhältnis $\ell/r_\text{scale}$ für sphärische Quellen (Bulge) |
| $c_\text{disk}$ | 2D geometrische Konstante | Verhältnis $\ell/R_\text{scale}$ für Scheiben- und Ringquellen |
| $c_\text{arm}$ | Geometrische Konstante der Spirale | Verhältnis $\ell/R_d$ für den azimutal konzentrierten Armüberschuss |
| $\lambda$ | Globale Wellenfeldkopplung | Skaliert die gesamte Wellenfelddichte |
Universalität der Parameter
Alle fünf Parameter sind global. Die gleichen numerischen Werte gelten für die Milchstraße, für zwergförmige irreguläre Galaxien und für massive Spiralen. Die galaxiespezifischen Informationen fließen nur über die fünf Beobachtungsdaten $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$ ein. Das Modell enthält keine pro Galaxie einstellbaren Parameter.
8. Die unidirektionale Natur der Berechnung
Eine offene Kette – keine Rückkopplung
Die gesamte Berechnung fließt von den Eingaben zu den Ausgaben, in eine Richtung. Die photometrischen und 21-cm-Beobachtungen bestimmen die baryonische Zersetzung. Die baryonische Zersetzung bestimmt die Wellenfelddichte. Die Wellenfelddichte bestimmt die eingeschlossene Wellenmasse. Die eingeschlossene Wellenmasse bestimmt die vorhergesagte Rotationskurve. Zu keinem Zeitpunkt beeinflusst die Rotationskurve einen früheren Schritt der Berechnung.
Diese Unidirektionalität hat drei wichtige Konsequenzen.
(a) Sobald die fünf Parameter auf Theorieebene festgelegt sind, ist die Rotationskurve eine strenge Vorhersage und keine Anpassung. Der Vergleich mit der beobachteten Rotationskurve ist ein Test, keine Kalibrierung.
(b) Das Modell hat keinen Mechanismus für eine galaxienweise Anpassung. Jede Änderung der Rotationskurvenvorhersage muss durch eine Änderung des Eingabevektors $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$ oder durch eine Änderung der universellen Theorieparameter $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$ erfolgen.
(c) Die Kalibrierung von $lambda$ an einer Referenzgalaxie ist nicht dasselbe wie die Anpassung an die Rotationskurve dieser Galaxie. Die Kalibrierung bestimmt eine einzige globale Zahl; die Rotationskurve bei allen anderen Radien der Referenzgalaxie und die Rotationskurven aller anderen Galaxien sind dann strenge Vorhersagen des kalibrierten Rahmens.
9. Die Rolle der zentralen Oberflächendichte (Anmerkung XI – Überarbeitung)
Die Diagnose in Anmerkung XI hat gezeigt, dass der restliche Vorhersagefehler stark mit der zentralen baryonischen Oberflächendichte $\Sigma_d$ korreliert, unabhängig von der Skalenlänge der Scheibe $R_d$. Die oben vorgestellte Formalisierung ist die Version des Modells , bevor diese Erkenntnis eingearbeitet wurde – sie verwendet nur $R_d$ in den Kohärenzlängenausdrücken $\ell_i = c_i\,R_d$.
Wo die Verfeinerung einsetzt
Im verfeinerten Modell werden die Kohärenzlängen $\ell_i$ sowohl von $R_d$ als auch von $\Sigma_d$ abhängen, wobei die streng lineare Beziehung $\ell_i = c_i\,R_d$ durch eine Funktion $\ell_i = c_i\,R_d\,\phi(\Sigma_d/\Sigma_\{ref})$ ersetzt wird, die den in Anmerkung XI identifizierten Rest absorbiert. Die funktionale Form von $\phi$ und ihre Parameter werden in den folgenden Anmerkungen bestimmt, zunächst anhand des Kalibrierungssatzes von 22 Galaxien und dann durch blinde Vorhersage anhand der verbleibenden SPARC-Stichprobe validiert.
Die unidirektionale Struktur der Berechnung wird durch diese Verfeinerung beibehalten: $\Sigma_d$ ist eine Beobachtungseingabe, die modifizierten Kohärenzlängen fließen in dieselben Faltungsintegrale ein, und die Rotationskurve entsteht wie zuvor. Es wird nur ein operatives Glied hinzugefügt – die Abhängigkeit von $\ell_i$ von einer zweiten Beobachtungsgröße.
10. Zusammenfassung der Methodik
1. Eingaben. Fünf Observablen pro Galaxie: Hubble-Typ $T$, Scheibengröße $R_d$, Oberflächenhelligkeit $\Sigma_d$, HI-Masse $M_\text{HI}$ und das universelle stellare Masse-zu-Licht-Verhältnis $\Upsilon_\star$.
2. Baryonische Zersetzung. Fünf Komponenten: Bulge (wenn $T \leq 4$), dünne Scheibe, dicke Scheibe, Gasring, Spiralarmüberschuss. Jede Komponente hat ein analytisches Dichteprofil.
3. Wellenkern. Universelle Yukawa-Form $\mathcal{K}_i(D) = K_0\,(1 + \alpha_i D)\,e^{-\alpha_i D}/D^2$ mit Kohärenzlänge $\ell_i = c_i\,R_\text{scale}$, die durch die geometrische Ausdehnung jeder Komponente bestimmt wird.
4. Faltung. Jede Komponente erzeugt eine Wellenfelddichte über ein eindimensionales Integral über Ringe (2D-Komponenten) oder Schalen (3D-Wulst). Die gesamte Wellenfelddichte ist die Summe der fünf Komponenten, skaliert durch die globale Kopplung $\lambda$.
5. Ausgabe. Die eingeschlossene Wellenmasse $M_\text{wave}(R)$ wird integriert und mit der baryonischen Geschwindigkeit $V_\text{bar}(R)$ kombiniert, um die vorhergesagte Rotationskurve $V_c(R)$ zu erhalten.
6. Parameter auf Theorieebene. $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$ – universal, keine Abstimmung pro Galaxie. Eine untersuchte Verfeinerung wird eine Abhängigkeit von $\Sigma_d$ hinzufügen.
7. Richtung. Eingaben → Baryonen → Wellenfeld → Rotationskurve. Keine Rückkopplung. Die Rotationskurve ist eine Vorhersage, keine Anpassung.
Referenzen. Lelli, F., McGaugh, S. S., Schombert, J. M. – SPARC: Massenmodelle für 175 Scheibengalaxien mit Spitzer-Photometrie und präzisen Rotationskurven, AJ 152, 157 (2016). – Freeman, K. C. – Über die Scheiben von Spiral- und S0-Galaxien, ApJ 160, 811 (1970). – Hernquist, L. – Ein analytisches Modell für sphärische Galaxien und Bulges, ApJ 356, 359 (1990). – Broeils, A. H., Rhee, M.-H. – Kurze 21-cm WSRT-Beobachtungen von Spiralgalaxien und irregulären Galaxien, A&A 324, 877 (1997). – McGaugh, S. S. – Das dritte Gesetz der galaktischen Rotation, Galaxies 2, 601 (2014). – Bovy, J., Rix, H.-W. – Eine direkte dynamische Messung des Dichteprofils der Scheibenoberfläche der Milchstraße, ApJ 779, 115 (2013). – Arnett, D. – Supernovae und Nukleosynthese, Princeton (1996). – Dutertre, X. – Bee Theory™: Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).
BeeTheory.com – Wellenbasierte Quantengravitation – Galaktische Methodik – © Technoplane S.A.S. 2026