BeeTheory – Grundlagen – Technischer Hinweis XI
Identifizieren Sie den fehlenden Parameter:
Schritt 1 – Systematische Korrelationsanalyse
Bevor wir das Modell ändern, stellen wir in dieser Notiz fest, welcher beobachtbare Parameter den Restfehler am besten vorhersagt. Wir arbeiten mit dem 22-Galaxien-Kalibrierungsset aus Anmerkung VIII und testen die Korrelation des Vorhersagefehlers mit jeder physikalisch bedeutsamen Variable und dann mit jeder bivariaten Kombination, um genau zu ermitteln, was das aktuelle Modell ausgelassen hat.
1. Das Ergebnis zuerst
Der fehlende Parameter ist die zentrale Oberflächendichte
Die zentrale baryonische Oberflächendichte $\Sigma_d$ hat die stärkste nicht-triviale Korrelation mit dem Vorhersagefehler: $r = +0,62$, $R^2 = 0,39$ für sich allein.
Die Kombination von $\Sigma_d$ mit der Scheibengröße $R_d$ in einem bivariaten Modell erklärt $R^2 = 0,43$ der Restvarianz, verglichen mit $R^2 = 0,07$ mit $R_d$ allein. Der RMS-Rest fällt von $19,5\%$ auf $14,9\%$.
Nachdem sowohl $R_d$ als auch $\Sigma_d$ absorbiert wurden, trägt keine zusätzliche physikalische Beobachtung Informationen über den Restwert.
2. Methode
Wenn wir mit dem Kalibrierungssatz von 22 Galaxien arbeiten (Anmerkung VIII), haben wir für jede Galaxie den Vorhersagefehler $text{err} = (V_text{tot} – V_f)/V_f$ und eine Liste von messbaren physikalischen Parametern. Wir berechnen die Pearson und Spearman Korrelationen zwischen dem Fehler und jeder Kandidatenvariablen und testen dann bivariate Regressionen der Form:
$$\text{err}(\%) \;=\; a \cdot R_d \;+\; b \cdot X \;+\; c$$
wobei $X$ jede Kandidatenvariable ist. Die beste $X$ ist diejenige, die die erklärte Varianz $R^2$ für die 22 Galaxien maximiert. Selbstreferenzielle Variablen – also solche, die aus der Modellausgabe abgeleitet sind, wie $V_\text{wave}$ oder $V_\text{tot}$ – werden von der Suche ausgeschlossen, da ihre Korrelation mit dem Fehler tautologisch ist.
3. Univariate Korrelationen
Die 24 getesteten Kandidatenvariablen, geordnet nach absoluter Pearson-Korrelation mit dem Fehler. Gold schattierte Zeilen sind Variablen, die aus dem Modell selbst abgeleitet sind (tautologisch); rot schattierte Zeilen sind echte physikalische Beobachtungsgrößen mit $|r| > 0,5$.
| Variabel | Beschreibung | Birnbaum $r$ | $p$-Wert | Bedeutung |
|---|---|---|---|---|
| Vw_über_Vf | Verhältnis Vw / Vf | +0.974 | 0.0000 | ★★★ |
| V_dynamisch | V_dyn = √(GM_bar/Rd) | +0.632 | 0.0021 | ★★★ |
| log_Sigma_d | log₁₀(Σ_d) | +0.622 | 0.0026 | ★★★ |
| M_gas | Gasmasse (M_sun) | +0.609 | 0.0034 | ★★★ |
| M_HI | HI-Masse (M_sun) | +0.609 | 0.0034 | ★★★ |
| T | Hubble-Typ | -0.585 | 0.0053 | ★★ |
| Vbar | Baryonisch Vbar (km/s) | +0.582 | 0.0057 | ★★ |
| M_bar_over_Rd2 | M_bar / Rd² | +0.559 | 0.0084 | ★★ |
| Vtot | Vorausgesagte Vtot (km/s) | +0.555 | 0.0090 | ★★ |
| Vw | Welle Vw (km/s) | +0.550 | 0.0098 | ★★ |
| Vbar_über_Vf | Verhältnis Vbar / Vf | +0.519 | 0.0158 | ★★ |
| log_M_gas | log₁₀(M_gas) | +0.506 | 0.0193 | ★★ |
| log_M_bar | log₁₀(M_bar) | +0.505 | 0.0196 | ★★ |
| M_bar | Baryonische Masse (M_sun) | +0.498 | 0.0214 | ★★ |
| log_M_star | log₁₀(M_star) | +0.449 | 0.0414 | ★★ |
| Sigma_d | Oberflächendichte (L/pc²) | +0.426 | 0.0544 | ★★ |
| M_stern_über_Rd2 | M_star / Rd² | +0.426 | 0.0544 | ★★ |
| M_star | Stellare Masse (M_sun) | +0.389 | 0.0815 | ★ |
Lesen Sie die Tabelle
Die höchste einzelne Korrelation ist $V_\text{wave}/V_f = +0,974$. Das ist tautologisch: Konstruktionsbedingt skaliert der Fehler direkt mit $V_\text{wave}$, so dass diese Variable einfach die Struktur der Vorhersageformel widerspiegelt und nicht einen externen physikalischen Antrieb.
Unter den echten physikalischen Observablen sind die höchsten Korrelationen $\log(\Sigma_d) = +0.622$, $V_\text{dynamisch} = +0.632$, $M_\text{gas} = +0.609$ und Hubble-Typ $T = -0.585$. Diese vier Signale stehen in einem physikalischen Zusammenhang: dichte Scheiben sind tendenziell massereicher, von früherem Typ und haben eine höhere baryonische dynamische Geschwindigkeit. Die Frage ist, welches der grundlegende Treiber ist.
4. Herausfiltern der redundanten Variablen
Mehrere der am stärksten korrelierten Variablen sind selbst stark mit $R_d$ korreliert, der Variable, von der bereits bekannt ist, dass sie den Fehler verursacht. Die Frage ist, welche davon unabhängige Informationen enthält.
| Variabel | Korrelation mit $R_d$ | Status |
|---|---|---|
| $\log(M_\star)$ | $r = +0.88$ | Redundant mit $R_d$ |
| $\log(M_\text{bar})$ | $r = +0.87$ | Redundant mit $R_d$ |
| $\log(M_\text{gas})$ | $r = +0.86$ | Redundant mit $R_d$ |
| Hubble-Typ $T$ | $r = -0.66$ | Teilweise redundant |
| $V_\text{dynamisch}$ | $r = +0.50$ | Teilweise unabhängig |
| $M_\text{bar}/R_d^2$ | $r = -0.19$ | Unabhängig |
| $\log(\Sigma_d)$ | $r = +0.10$ | Unabhängig |
Die Massen korrelieren nahezu perfekt mit $R_d$: eine größere Scheibe enthält einfach mehr baryonisches Material. Diese Variablen tragen also im Wesentlichen die gleiche Information wie $R_d$ selbst. Im Gegensatz dazu sind $\Sigma_d$ (zentrale Oberflächendichte) und $M_\text{bar}/R_d^2$ (mittlere baryonische Oberflächendichte) in dieser Stichprobe fast orthogonal zu $R_d$: Sie erfassen die strukturelle Eigenschaft „wie kompakt die Materie ist“, unabhängig davon, „wie ausgedehnt die Scheibe ist“.
5. Fehler versus Oberflächendichte – Visualisierung
Aufzeichnung des Fehlers gegen $\log_{10}(\Sigma_d)$ allein, eingefärbt nach Hubble-Typ:
Der Trend ist eindeutig und monoton: Galaxien mit einer höheren zentralen Oberflächendichte werden von der BeeTheory systematisch überschätzt, während diffuse Scheiben mit geringer Dichte unterschätzt werden. Die Anpassungsneigung von $33$ Prozentpunkten pro Dekade von $\Sigma_d$ stimmt über den gesamten Bereich von 15 bis 605 $L_\odot/\text{pc}^2$ gut mit den Daten überein.
6. Bivariate Modelle – Vergleich
Wenn Sie $R_d$ zu jeder Kandidatenvariable hinzufügen, ergibt sich eine klarere Rangfolge. Die folgende Tabelle zeigt die erklärte Varianz $R^2$, wenn $R_d$ mit jeder zweiten Variable gepaart wird (tautologische Kombinationen ausgeschlossen):
| Bivariates Modell | $R^2$ | RMS-Restwert | Anmerkungen |
|---|---|---|---|
| $\text{err} = a R_d + c$ (univariate Basislinie) | 0.074 | $19.5\%$ | Referenz, keine zweite Variable |
| $\text{err} = a R_d + b f_\text{gas} + c$ | 0.101 | $19.3\%$ | Vernachlässigbare Verbesserung |
| $\text{err} = a R_d + b \log M_\star + c$ | 0.272 | $17.3\%$ | – |
| $\text{err} = a R_d + b V_\text{bar} + c$ | 0.345 | $16.4\%$ | – |
| $\text{err} = a R_d + b \log M_\text{gas} + c$ | 0.359 | $16.3\%$ | – |
| $\text{err} = a R_d + b T + c$ | 0.367 | $16.2\%$ | – |
| $\text{err} = a R_d + b \log M_\text{bar} + c$ | 0.373 | $16.1\%$ | – |
| $\text{err} = a R_d + b\,V_\text{dynamisch} + c$ | 0.402 | $15.7\%$ | Stark |
| $\text{err} = a R_d + b \log\Sigma_d + c$ | 0.430 | $15.3\%$ | Unabhängig von $R_d$ |
| $\text{err} = a R_d + b (M_\text{bar}/R_d^2) + c$ | 0.459 | $14.9\%$ | Bestes nicht-tautologisches Modell |
Das beste bivariate Modell
$$\text{err}(\%) \;=\; a\,R_d \;+\; b\,\frac{M_\text{bar}}{R_d^2} \;+\; c, \qquad R^2 = 0.46$$
Die Variable $M_\text{bar}/R_d^2$ ist die mittlere baryonische Oberflächendichte der Scheibe, $\langle \Sigma_\text{bar} \rangle = M_\text{bar}/(\pi R_d^2)$. Sie enthält Informationen darüber, wie kompakt die sichtbare Materie ist, unabhängig davon, wie groß die Scheibe ist. Dies ist die Variable, die die BeeTheory derzeit nicht berücksichtigt.
7. Abschlussprüfung – was bleibt übrig, nachdem $R_d$ und $\Sigma_d$ berücksichtigt worden sind
Wenn $R_d$ und $\log \Sigma_d$ zusammen den strukturellen Defekt erfassen, sollte das Residuum der bivariaten Anpassung mit jeder physikalischen Beobachtung unkorreliert sein. Dies zu testen ist die formale Abschlussprüfung:
| Variabel | Korrelation mit Residuum | Status |
|---|---|---|
| $R_d$ | $+0.00$ | Nach Konstruktion |
| $\log \Sigma_d$ | $+0.00$ | Nach Konstruktion |
| $\log M_\star$ | $-0.05$ | Absorbiert |
| $\log M_\text{bar}$ | $+0.07$ | Absorbiert |
| $\log M_\text{gas}$ | $+0.14$ | Absorbiert |
| Hubble-Typ $T$ | $-0.04$ | Absorbiert |
| $V_\text{dynamisch}$ | $+0.08$ | Absorbiert |
| $V_\text{bar}$ | $+0.05$ | Absorbiert |
| $f_\text{gas}$ | $+0.28$ | Geringfügig; unter der Signifikanzgrenze |
Nach Berücksichtigung von $R_d$ und $\log \Sigma_d$ bleibt keine physikalische Beobachtungsgröße mehr signifikant mit dem Restfehler korreliert. Die strukturelle Information im Fehler wurde von diesen beiden Variablen vollständig erfasst. Die verbleibende RMS-Streuung von $15%$ steht im Einklang mit der Beobachtungsunsicherheit bei den SPARC-Eingabeparametern und mit der intrinsischen Variabilität von Galaxie zu Galaxie, die von keinem dieser aggregierten Deskriptoren erfasst wird.
8. Physikalische Interpretation
Das aktuelle BeeTheory-Modell verwendet die Skalenlänge der Scheibe $R_d$ an zwei Stellen: als räumliche Skala der baryonischen Verteilung (das Exponentialprofil $Sigma propto e^{-R/R_d}$) und als Kohärenzlänge des Wellenkerns ($ell = c_text{disk},R_d$). Die Amplitude des baryonischen Profils $\Sigma_0$ ist implizit und so skaliert, dass sie nach der Integration die korrekte Sternmasse ergibt.
Was die Oberflächendichte physikalisch darstellt
Die mittlere baryonische Oberflächendichte $langle Sigma_text{bar} rangle = M_text{bar}/(pi R_d^2)$ ist die Masse pro Flächeneinheit der Scheibe. Zwei Galaxien mit gleichem $R_d$ aber unterschiedlichem $\Sigma_d$ haben die gleiche geometrische Ausdehnung, aber unterschiedliche Mengen an eingepackter Materie. Das aktuelle Modell behandelt nur die geometrische Ausdehnung ($R_d$) als relevant für die Wellenkohärenzlänge und ignoriert, wie konzentriert die Materie ist. Dies ist genau der Parameter, den die Restanalyse als fehlend identifiziert.
Die Richtung des Effekts
Die Korrelation ist positiv: der Fehler wächst mit der Oberflächendichte. Das bedeutet, dass bei festem $R_d$ dichtere Scheiben vom Modell überschätzt werden – das Wellenfeld ist relativ zur Rotationskurve zu stark. Umgekehrt sagt das Modell für einen gegebenen $R_d$ diffuse Scheiben mit geringer Dichte zu wenig voraus. Eine plausible physikalische Interpretation: Die Wellenkohärenzlänge sollte nicht nur von der geometrischen Ausdehnung der Quelle abhängen, sondern auch von ihrer Konzentration, wobei dichtere Materie eine stärker lokalisierte Wellenantwort erzeugt. Dies würde natürlich die Amplitude des Wellenfeldes in Scheiben mit hohem $\Sigma$ unterdrücken und sie in solchen mit niedrigem $\Sigma$ verstärken.
9. Zusammenfassung von Schritt 1
1. Bei dem Kalibrierungsset mit 22 Galaxien korreliert der Vorhersagefehler am stärksten mit der zentralen Oberflächendichte $\Sigma_d$ ($r = +0,62$) unter den echten physikalischen Observablen.
2. Andere Variablen, die auf den ersten Blick stark korreliert erscheinen (Sternmasse, Gasmasse, baryonische Masse), erweisen sich als hochgradig redundant mit $R_d$ (Korrelationen $\geq 0,86$ mit $R_d$) und enthalten daher wenig neue Informationen.
3. Das beste nicht-tautologische bivariate Modell ist $\text{err} = a\,R_d + b\,(M_\text{bar}/R_d^2) + c$, mit $R^2 = 0,46$ und RMS-Rest $14,9\%$. Die zweite Variable ist die mittlere baryonische Oberflächendichte der Scheibe.
4. Nach Berücksichtigung von $R_d$ und $\Sigma_d$ behält keine andere beobachtbare Größe eine signifikante Korrelation mit dem Residuum. Die Diagnose ist abgeschlossen.
5. Der fehlende Parameter ist identifiziert: Das aktuelle BeeTheory-Modell berücksichtigt die geometrische Ausdehnung der baryonischen Verteilung ($R_d$), aber nicht ihre Oberflächendichte ($\Sigma_d$). Der nächste Schritt besteht darin, $\Sigma_d$ als zweite Eingabe für die Wellenkohärenzlänge einzubeziehen und dann das Modell auf die 22-Galaxien-Menge zu übertragen.
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). – Pearson, K. – Mathematische Beiträge zur Theorie der Evolution III, Phil. Trans. R. Soc. A 187, 253 (1896). Korrelationskoeffizient. – Dutertre, X. – Bee Theory™: Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).
BeeTheory.com – Wellenbasierte Quantengravitation – Diagnoseschritt 1 – © Technoplane S.A.S. 2026