BeeTheory – 基礎 – テクニカルノート XI

欠落したパラメータの特定
ステップ1 – 系統的相関分析

このノートでは、モデルを修正する前に、どの観測可能なパラメータが最もよく残差 を予測するかを診断します。注Ⅷの22銀河の較正セットを使って、物理的に意味のあるすべての変数と予測誤差の相関をテストし、次にすべての二変数の組み合わせと相関をテストして、現在のモデルが何を省略しているかを厳密に特定します。

1.結果はまず

欠けているパラメータは中心表面密度

中心バリオン表面密度$Sigma_d$は、予測誤差と最も強い非自明な相関があります:$r = +0.62$、$R^2 = 0.39$。

R_d$単独では$R^2 = 0.07$であるのに対し、$Sigma_d$とディスクサイズ$R_d$を組み合わせた二変量モデルでは$R^2 = 0.43$の残差分散を説明。RMS残差は$19.5%から$14.9%に減少。

R_d$と$Sigma_d$の両方を吸収した後は、残差に関する情報を伝える追加の物理観測量はありません。

2.方法

22銀河の較正セット(注VIII)を使って、各銀河の予測誤差$text{err} = (V_text{tot} – V_f)/V_f$と測定可能な物理パラメータのリストを用意します。誤差と各候補変数のピアソン相関とスピアマン相関を計算し、二変量回帰を検定します:

err}(\;==)a♪R_d♪+♪b♪X♪+♪c$$。

ここで、$X$はそれぞれの候補変数です。X$は22個の銀河の説明分散$R^2$を最大にするものが最適です。V_text{wave}$や$V_text{tot}$のような自己言及変数(モデルの出力から得られる変数)は、誤差との相関がトートロジーなので、探索から除外されます。

3.一変量相関

テストされた24の候補変数。誤差との絶対ピアソン相関でランク付け。金色の網掛けをした行はモデル自身から導かれた変数(トートロジー)、赤の網掛けをした行は$|r| > 0.5$の本物の物理観測値。

可変 説明 ピアソン $r$ p$値 意義
Vw_over_VfVw/Vf比+0.9740.0000★★★
V_dynamicalV_dyn = √(GM_bar/Rd)+0.6320.0021★★★
log_Sigma_dlog₁₀(Σ_d)+0.6220.0026★★★
M_ガスガス質量 (M_sun)+0.6090.0034★★★
M_HIHI質量 (M_sun)+0.6090.0034★★★
Tハッブル型-0.5850.0053★★
バールバリオンVbar (km/s)+0.5820.0057★★
M_bar_over_Rd2M_bar / Rd²+0.5590.0084★★
ボルト予想Vtot (km/s)+0.5550.0090★★
ブイ波の速度(km/s)+0.5500.0098★★
Vbar_over_VfVbar/Vf比+0.5190.0158★★
ログ_M_ガスlog₁₀(M_gas)+0.5060.0193★★
ログ_M_バーlog₁₀(M_bar)+0.5050.0196★★
M_barバリオン質量 (M_sun)+0.4980.0214★★
ログ_M_スターlog₁₀(M_star)+0.4490.0414★★
シグマ表面密度 (L/pc²)+0.4260.0544★★
M_star_over_Rd2M_star / Rd²+0.4260.0544★★
M_star恒星質量 (M_sun)+0.3890.0815

表を読むと

最も高い相関は$V_text{wave}/V_f = +0.974$です。これはトートロジーです。構造上、誤差は$V_text{wave}$で直接スケールするので、この変数は単に予測式の構造を反映しているだけで、外部の物理的なドライバーではありません。

本物の物理観測量の中で最も相関が高いのは、$log(˶‾᷅˵)=+0.622$、$V_text{dynamical} = +0.632$、$M_text{gas} = +0.609$、ハッブル型$T = -0.585$です。これらの4つのシグナルは物理的に関連しています。密度が高い円盤は質量が大きく、型が早く、バリオン動力学速度が速い傾向があります。問題は、どれが根本的な原動力なのかということです。

4.冗長変数のフィルタリング

相関の高い変数のうちのいくつかは、それ自体が$R_d$と強い相関を持っています。問題は、どれが独立した情報を持っているかということです。

可変 R_d$との相関 ステータス
log(M_star)$$r = +0.88$R_d$と冗長
log(M_text{bar})$.$r = +0.87$R_d$と冗長
log(M_text{gas})$$r = +0.86$R_d$と冗長
ハッブル型 $T$$r = -0.66$一部冗長
V_text{dynamical}$$r = +0.50$一部独立
M_text{bar}/R_d^2$$r = -0.19$独立系
log(˶´⚰︎`˵)$.$r = +0.10$独立系

質量は$R_d$とほぼ完全に相関しています。したがって、これらの変数は$R_d$そのものと本質的に同じ情報を持っています。これに対して、$Sigma_d$(中心表面密度)と$M_text{bar}/R_d^2$(平均バリオン表面密度)は、このサンプルでは$R_d$とほとんど直交しています:これらは、「円盤がどれだけ広がっているか」とは関係なく、「物質がどれだけコンパクトか」という構造的性質をとらえています。

5.表面密度に対する誤差 – 可視化

log_{10}( \Sigma_d)$ だけに対する誤差をハッブル型で色分けしてプロット:

予測誤差と中心面密度$Sigma_d$ – 22個の較正銀河 101001000 -50%-25%+0%+25%+50% 中心面密度 $Sigma_d$ ($L_Lodot/$pc$^2$) – log scale 予測誤差 err = 33.0-log(Σd) – 60ピアソン r = 0.622, R² = 0.387 D631-7 Σd=115 err=-12%.DDO064 Σd=120 誤差=+13DDO064 Σd=120 err=+13% DDO154 Σd=45 err=+4DDO161 Σd=35 err=+11DDO168 Σd=100 err=-21DDO170 Σd=25 err=+15ESO116-G012 Σd=115 err=+14ESO444-G084 Σd=60 err=+6F561-1 Σd=30 err=-15F563-1 Σd=20 err=-18F563-V1 Σd=25 err=-36F563-V2 Σd=30 err=-27F565-V2 Σd=18 err=-39% F567-2 Σd=30 err=-27F565-V2 Σd=18 err=-39% F567-2 Σd=15 err=-23F568-1 Σd=40 err=-8F568-3 Σd=35 err=-12F568-V1 Σd=20 err=-24F571-8 Σd=50 err=+14% F574-1 Σd=40 err=-8F574-1 Σd=30 err=-2% NGC2841 Σd=35 err=-12F574-1 Σd=30 err=-2% NGC2841 Σd=605 err=+18NGC3198 Σd=153 err=+43% Σd=30 err=-2% F574-1 S0-SaSb-SbcSc-ScdSd-Im
一変量フィット誤差$( \%) = 33log( \Sigma_d) – 60$, Pearson $r = 0.62$, $R^2 = 0.39$。

中心面密度の高い銀河はBeeTheoryによって系統的に過大に予測され、拡散した低密度の円盤は過小に予測されます。シグマ$d$の10年あたり$33$パーセントポイントというフィットの傾きは、15から605$L_odot/text{pc}^2$の全範囲でデータとよく一致します。

6.二変量モデル – 比較

各候補変数に$R_d$を加えると、より明確なランキングが得られます。下の表は、$R_d$を各第2変数とペアにしたときの説明される分散$R^2$です(同語反復の組み合わせは除く):

二変量モデル $R^2$ 実効残差 備考
text{err} = a R_d + c$ (一変量ベースライン)0.074$19.5\%$参照、2番目の変数なし
a R_d + b f_text{gas}+ c$0.101$19.3\%$ごくわずかな改善
text{err} = a R_d + b ୧log M_interstar + c$.0.272$17.3\%$
$text{err} = a R_d + b V_text{bar}+ c$0.345$16.4\%$
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{dynamical} + c$.0.402$15.7\%$強い
text{err} = a R_d + b Γ + c$.0.430$15.3\%$R_d$に依存せず
text{err} = a R_d + b (M_text{bar}/R_d^2) + c$.0.459$14.9\%$最高の非同語モデル

最良の2変量モデル

$$\text{err}(\%) \;=\; a\,R_d \;+\; b\,\frac{M_\text{bar}}{R_d^2}\c, \qquad R^2 = 0.46$$.

変数$M_text{bar}/R_d^2$は円盤の平均バリオン表面密度、$langle\rangle = M_text_bar}/(Γpi R_d^2)$。これは、円盤の大きさとは関係なく、目に見える物質がどれだけコンパクトかという情報を持っています。これはBeeTheoryが現在説明できない変数です。

7.閉包チェック – $R_d$ と $Sigma_d$ を考慮した後に残るもの。

もし$R_d$と$log \Sigma_d$が一緒に構造的欠陥を捉えているならば、二変量フィットの残差はあらゆる物理的観測量と相関がないはずです。これをテストするのが正式な閉包チェックです:

可変 残差との相関 ステータス
R_d$$+0.00$工事別
log$+0.00$工事別
log M_star$$-0.05$吸収
log M_text{bar}$$+0.07$吸収
log M_text{gas}$$+0.14$吸収
ハッブル型 $T$$-0.04$吸収
V_text{dynamical}$$+0.08$吸収
V_text{bar}$$+0.05$吸収
f_text{gas}$。$+0.28$限界;有意水準以下

R_d$と$log ↪Sigma_d$ を考慮した後では、どの物理観測変数も残差誤差と有意な相関を保ちません。誤差の構造情報はこの2つの変数によって完全に捕捉されました。残りの$15%$のRMSのばらつきは、SPARCの入力パラメータに関する観測の不確かさと、これらの総体的な記述子のどれにも捕らえられない内在的な銀河間の変動と一致します。

8.物理的解釈

現在のBeeTheoryモデルでは、円盤のスケール長$R_d$をバリオン分布の空間スケール(指数分布$Sigma propto e^{-R/R_d}$)と波動カーネルのコヒーレンス長($ell = c_text{disk},R_d$)の2箇所で使っています。バリオン分布$Sigma_0$の振幅は暗黙的で、積分すると正しい恒星質量になるようにスケーリングされます。

表面密度が物理的に表すもの

平均バリオン表面密度 $langle Sigma_text{bar} rangle = M_text{bar}/(pi R_d^2)$ は、円盤の単位面積あたりの質量です。R_d$が同じで$Sigma_d$が異なる2つの銀河は、幾何学的な広がりは同じですが、詰まっている物質の量が異なります。現在のモデルでは、幾何学的な広がり($R_d$)だけが波動コヒーレンス長に関係しており、物質がどれだけ集中しているかは無視されています。これこそが、残差分析で欠落しているパラメータなのです。

効果の方向

この相関は正の相関です。つまり、$R_d$を固定した場合、密度の高い円盤はモデルによって過剰に予測されます。逆に、$R_d$が一定であれば、低密度の円盤は過小評価されます。もっともらしい物理的解釈:波のコヒーレンス長は、光源の幾何学的な広がりだけでなく、その濃度にも依存するはずです。これは当然、高シグマ円盤では波動場の振幅を抑制し、低シグマ円盤では増幅します。

9.ステップ1のまとめ

1.22銀河の較正セットにおいて、予測誤差は、真の物理観測量のうち、中心表面密度$Sigma_d$($r = +0.62$)と最も強く相関しています。

2.最初に強い相関があると思われた他の変数(恒星質量、ガス質量、バリオン質量)は、$R_d$と非常に冗長($R_d$との相関は$0.86$)であることがわかり、新しい情報はほとんどありません。

3.最良の非同語2変量モデルは$text{err} = a,R_d + b,(M_text“/bar}/R_d^2) + c$で、$R^2 = 0.46$、RMS残差$14.9%$。つ目の変数は円盤の平均バリオン表面密度。

4.R_d$と$Sigma_d$を考慮しても、他の観測量が残差と有意な相関を保つことはありません。診断終了。

5.現在のBeeTheoryモデルは、バリオン分布の幾何学的な広がり($R_d$)は説明しますが、表面密度($Sigma_d$)は説明しません。次のステップは、$Sigma_d$を波動コヒーレンス長への2番目の入力として組み込み、22銀河セットでモデルを再フィットすることです。


参考文献Lelli, F., McGaugh, S. S., Schombert, J. M. –SPARC: Mass Models for 175 DiskGalaxies with Spitzer Photometry and Accurate RotationCurves, AJ 152, 157 (2016).- Pearson, K. –Mathematical contributions to the theory of evolution III, Phil.Trans.R. Soc. A 187, 253 (1896)。相関係数。- Dutertre, X. –Bee Theory™:Wave-BasedModeling of Gravity, v2, BeeTheory.com (2023).

BeeTheory.com – 波動量子重力 – 診断ステップ1 – © Technoplane S.A.S. 2026