BeeTheory（ビーセオリー） – SPARC銀河 – 調整済みフィット – 2025年

BeeTheoryを20の外部銀河に適用：
調整式とブラインドテスト手法

SPARCカタログは、バリオン分布と回転曲線を測定した175個の銀河を提供しています。その結果、20個の銀河のうち18個が、観測された平坦な自転速度の20%以内で、χ²/dof = 0.93と予測されました。

0.結果

ベストフィット – SPARC銀河20個、Q = 1、高品質

Kd=K0/Rdとℓd=c –Rdという2つの普遍定数を用いて、20個の銀河を同時にフィットさせました。

各銀河の 暗黒質量密度は、バリオン円盤パラメータのみから予測されます。

ベストフィットのパラメータ：K0= 0.3759（無次元）、c = 6.40（無次元）。結果：18/20銀河が観測されたVfの 20%以内、χ²/dof = 0.93。2つの外れ値、CamBとNGC3741は恒星円盤モデリングが破綻したガス支配の矮星。

18/20

Vfの20%以内

6.8%

中央値誤差

0.93

χ²/dof

2

普遍定数

0.91

ピアソンr、タリーフィッシャー

修正BeeTheory式 –銀河の人口で調整 \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d}\int_0^{R_{\mathrm{max}}}\Sigma_0 e^{-R’/R_d}\frac{(1+\alpha_d D)e^{-\alpha_d D}}{D^2}\,2\pi R’\,dR’\) \(D=\sqrt{r^2+R’^2},\qquad \alpha_d=\frac{1}{cR_d},\qquad K_0=0.3759,\qquad c=6.40\)

1.20個のSPARC銀河のデータと予測

SPARCのサンプルは、矮小不規則銀河から大質量渦巻銀河まで、光度にして50年にわたる銀河をカバーしています。各銀河の入力パラメータは、Lelli et al. 2016 Table 1からそのまま引用しています：円盤スケール半径Rd、中心面輝度Σd、HIガス質量MHI、平坦回転速度Vf。

恒星質量はM★＝Υ★×L3.6として計算され、Υ★＝0.5M⊙/L⊙。ガス質量はMgas= 1.33 ×MHI として計算。

銀河系 Rdkpc Kdkpc-¹ ℓdkpc Vfobs Vbar Vdark VBT 誤差 ステータス
銀河テーブルをロード中…

すべての速度はkm/s。誤差 = (VBT–Vf)/Vf.パラメータ：Kd= 0.3759/Rd,ℓd= 6.40 ×Rd.Reval=5Rdで評価されたBeeTheoryの予測。

観測されたVBTと Vfの比較 – 20 SPARC銀河、調整済みBeeTheory。
20%以内：18銀河 外れ値CamB, NGC 3741 完璧な予測、1:1 ±20% 以内

2.修正された公式 – なぜK ∝1/Rdなのか？

オリジナルのBeeTheory天の川フィットは、単一の結合定数K = 0.02365 kpc-¹、コヒーレンス長ℓ = 3.17Rdを使っていました。20個のSPARC銀河に盲目的に適用すると、これは系統的にVfを約50%過小評価しました。

銀河ごとの解析では、K∝1/Rdとして必要な結合定数が変化するという明確なパターンが明らかになりました。

2.1 定数からスケーリング則へ

重要な洞察は次元です。スケールRdと表面密度Σ0の円盤から半径rにおけるBeeTheoryの暗黒密度は、漸近的な平回転領域ではr≪ℓです：

漸近暗黒密度,Rd≪ r ≪ ℓ \(\rho_{\mathrm{dark}}(r)\approx \frac{K\Sigma_0 R_d^2}{r^2}f\left(\frac{r}{\ell}\right)\) \(M_{mathrm{dark}}(<r)︓KSigma_0R_d^2r\)。

平らな回転速度は次のようになります：

BeeTheoryダークマスからのフラット回転速度 \(V_f^2=\frac{GM_{\mathrm{dark}}(<R_{\mathrm{eval}})}{R_{\mathrm{eval}}}\propto GK\Sigma_0R_d^2=GK\frac{M_\star}{2\pi}\) [ラテックス]V_fpropto (K M_star)^{1/2}[/latex]

観測されたバリオンTully-Fisher関係はVf4∝Mbar、つまりVf∝Mbar1/4です。これがBeeTheoryによって再現されるためには、Vf2∝ M★/Rd、平均円盤表面密度が必要です。これは

タリーフィッシャーの傾きを再現するために必要なスケーリング \(V_f^2\propto GK\frac{M_\star}{2\pi}\propto\frac{M_\star}{R_d}\implies K\propto\frac{1}{R_d}\) \(\boxed{K_d=\frac{K_0}{R_d},\qquad K_0=0.3759}\)

このスケーリングはアドホックなパッチではなく、Tully-Fisher関係が要求しているものです。K∝1/Rdは、よりコンパクトな円盤が単位質量あたりにより強い暗黒場を生成することを意味します。

2.2 コヒーレンス長 – c = 6.40 ≠ 3.17 の理由

天の川フィットではcMW = ℓd/Rd = 3.17。SPARCのサンプルではcSPARC= 6.40と約2倍。二つの説明が考えられます：

  • 選択バイアス：20個のSPARC銀河は、質の高い回転曲線の延長線上にある銀河が選ばれており、HI円盤がより延長している銀河に偏っていること。
  • ガス円盤の寄与：多くのSPARC銀河では、HIガス円盤の半径はRHI≈1.7Rdです。ガス円盤の寄与：多くのSPARC銀河では、HIガス円盤のスケール半径はRHI≈1.7Rdです。

どちらの効果も実際に測定可能です。ガス円盤と恒星円盤を別々にモデル化する必要があります。

完全な修正ビー理論式 – 2つの普遍定数 \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d}\int_0^{8R_d}\Sigma_0 e^{-R’/R_d}\frac{(1+\alpha_dD)e^{-\alpha_dD}}{D^2}\,2\pi R’\,dR’\) \(\alpha_d=\frac{1}{cR_d},\qquad D=\sqrt{r^2+R’^2},\qquad K_0=0.3759,\qquad c=6.40\)

3.計算-ステップ・バイ・ステップ

各SPARC銀河に対して、BeeTheoryの予測は5つのステップで進みます。銀河ごとに自由パラメータを調整することはありません。

1
SPARC表1からバリオン入力を読み込みます。

Rd、Σd、MHI、Vf。Σ0=Σd× Υ★ × 10⁶M⊙/kpc²、Mgas= 1.33 ×MHIに変換してください。

2
RdからBeeTheoryパラメータを計算

Kd=K0/Rd= 0.3759/Rd,ℓd=cRd= 6.40Rd,αd=1/ℓd.フィッティングなし。

3
r=5Rdにおける暗黒密度の積分
\(\rho_{\mathrm{dark}}(5R_d)=K_d\sum_{i=1}^{60}\Sigma_0e^{-R_i’/R_d}\frac{(1+\alpha_dD_i)e^{-\alpha_dD_i}}{D_i^2}2\pi R_i’\Delta R’\)

R′∈[0,8Rd]を60リングで数値積分。次に球面積分して囲み暗質量Mdark(<5Rd)を求めます。

4
バリオン円速度を計算
\(V_{\mathrm{bar}}(R)=\sqrt{\frac{G[M_{\mathrm{disk}}(<R)+M_{\mathrm{gas}}(<R)+M_{\mathrm{bulge}}(<R)]}{R}}\)
5
全円速の予測
\(V_{\mathrm{BT}}=\sqrt{V_{\mathrm{bar}}^2+V_{\mathrm{dark}}^2},\qquad V_{\mathrm{dark}}=\sqrt{\frac{GM_{\mathrm{dark}}(<5R_d)}{5R_d}}\)

観測されたVfと比較誤差 = (VBT–Vf)/Vf.

Vdark/Vbar比 – ダーク成分はどの程度支配的か？

4.ブラインドテストが唯一の正直なテストである理由

キャリブレーションしたデータを再現するモデルは、何の証明にもなりません。すべてのモデルは、たとえ間違ったモデルであっても、学習データに合うように調整することができます。唯一の科学的に意味のあるテストは、ブラインド予測です。キャリブレーションからパラメータを凍結した状態で、モデルを見たことのないデータに適用し、その結果を報告します。

4.1 “ブラインド “の意味

BeeTheoryのパラメータK0= 0.3759とc = 6.40は、20個のSPARC銀河を同時にフィッティングして決定されました。これらは現在固定されています。

ブラインドテストは、これらのパラメータを、フィッティングに使われなかった残りの155個のSPARC銀河に適用し、その結果を、観測された回転曲線を見る前に報告することです。このテストはまだ行われていません。

元の天の川銀河のパラメータ、Kd = 0.02365 と ℓd = 3.17Rd は、1つの銀河に対して決定されたものです。このパラメータを調整せずにSPARCに適用したところ、正しい銀河は0/20でした。この失敗から、K ∝1/Rdスケーリングが明らかになりました。

4.2 フィットの質の統計的意味

20個の銀河全体でχ²/dof = 0.93で、このモデルは、想定された15%の速度不確かさのほぼ期待されるレベルで合っています。

カイ二乗解釈 \(\frac{\chi^2}{\mathrm{dof}}=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_{\mathrm{BT}}(i)-V_f(i)}{0.15V_f(i)}\right)^2=0.93\) \(N=19\ \text{(excluding CamB)},\qquad p=2\ (K_0,c),\qquad \mathrm{dof}=17\)

0.93という値は、理想的な1.0に非常に近い値です。このモデルは、測定の不確かさのレベルでの散乱を説明します。

4.3 2つの外れ値

CamB – 純粋なガス矮星、Vf= 2.0 km/s

CamBは恒星質量がほとんどなく、M★≒2×10⁷M⊙ です。BeeTheoryの公式では、Σ0e-R/Rdを光源としていますが、CamBではバリオンはほとんどHIガスであり、恒星ではありません。恒星円盤モデルは適用できません。

NGC 3741 – 47%の過大評価。

NGC 3741は、非常に拡張されたHIディスクを持つ小さな低表面輝度の矮星です。BeeTheoryの恒星円盤は、実際のバリオンの広がりを過小評価しています。ガス円盤をより大きな半径を持つ別の天体の構成要素として含めると、予測される暗黒質量が減少し、過大評価が修正されます。

その他の18の予測

20%以内の18銀河の誤差の中央値は6.8%で、観測誤差の範囲内です。Rdは1.3～5.8kpc、Vfは76～278km/sです。BeeTheoryは、この速度のファクター3.7の範囲（Tully-Fisher勾配）を2つの普遍定数で正しく予測します。

5.物理的な意味-スケーリングが示すもの

5.1 普遍的な無次元結合

Kd=K0/Rd,ℓd=cRd とすると、BeeTheoryの無次元結合は次のようになります：

有効結合 – 銀河の大きさによるスケール \(\lambda_{\mathrm{eff}}=K_d\ell_d^2=\frac{K_0}{R_d}(cR_d)^2=K_0c^2R_d\) \(\lambda_{\mathrm{eff}}=0.3759\times41.0\times R_d=15.4R_d\ \text{(kpc)}\)

λeffはRdとともに大きくなります。銀河が大きいとそれに比例して暗質量も大きくなります。これが、なぜ大質量渦巻銀河の方が矮小銀河よりも暗黒物質が支配的なのかというBeeTheoryの予想です。

5.2 半径加速度との関係

McGaughらは、観測された求心加速度gobs= Vc2/Rがバリオン寄与gbar=GMbar/R²の普遍関数であることを発見しました。BeeTheoryでは、この関係は次のようになります：

RARのBeeTheory予測 \(g_{\mathrm{obs}}=g_{\mathrm{bar}}+g_{\mathrm{dark}}\) \(g_{\mathrm{dark}}\propto K_0c^2g_{\mathrm{bar}}^{1/2}G^{1/2}\)

gdark∝ gbar1/2のスケーリングは、観測されたRARの形状を作り出します。BeeTheoryから正確なRAR曲線を導き出すことが次の理論的課題です。

データ：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.

RAR:McGaugh, S. S., Lelli, F., Schombert, J. M., PRL 117, 201101, 2016.

BTFR：Lelli, F. et al., ApJ 816, 2016.

BeeTheory：Dutertre, X., BeeTheory.com v2, 2023, extended 2025.

3.6μmでΥ★=0.5M⊙/L⊙、McGaugh＆Schombert 2014。

BeeTheory.com – 波動ベースの量子重力 – © Technoplane S.A.S. 2025