BeeTheory – 基礎 – テクニカルノートXII
公式化
銀河系規模の蜂理論計算
このノートは、BeeTheoryの枠組みを円盤銀河に適用したものです。観測入力、バリオン分布の幾何学的分解、各成分の波動場を定義する積分方程式、そして予測される自転曲線を得るための一連の操作。観測されたバリオン構造が波動場を決定し、波動場が回転曲線を決定します。
1.1つのダイアグラムでの計算
一方向チェーン
観測された光度 $(Baryonic decomposition $(\rho_text{bar})$)
big
波動場畳み込み $( \rho_text`)
γ
質量積分$(M_text{wave})$ Enclosed wave mass $(M_text{wave})$
(M_text{wave})$。
ニュートン関係$(V_c)$ Predicted rotation curve $(V_c)$
ステップは反転せず。回転曲線$V_c(R)$を入力として使うことはありません。
2.観測入力
各銀河について、計算には5つの観測値が必要です。これらは銀河固有の量だけで、他はすべてこれらから計算されます。この段階では自転曲線に対するフィッティングは行いません。
| シンボル | 数量 | ソース |
|---|---|---|
| $T$ | ハッブル形態型 | カタログ（de Vaucouleurs et al.） |
| R_d$ | 恒星円盤のスケール長 (kpc) | スピッツァー3.6µm測光 (SPARC) |
| シグマ | 中心円盤表面の明るさ ($L_odot/text{pc}^2$) | スピッツァー3.6µm測光 (SPARC) |
| M_text{HI}$ | 全原子水素質量 ($M_odot$) | 21cm電波観測（SPARC） |
| Upsilon_star$ | 波長3.6µmにおける恒星の質量光度比 | 固定ユニバーサル：$0.5,M, L_odot$（マクゴー2014） |
天の川については、$R_d$, $Sigma_d$, $M_text{HI}$は、内部恒星サーベイ(Bovy & Rix 2013)と21cmマップから決定された類似の値で置き換えられています。同じ5量体入力ベクトルを使用。
3.バリオン分解 – 5つの幾何学的成分
5つの観測入力から、バリオン質量は5つの異なる幾何学的成分に分割されます。各成分は密度分布と特徴的なスケールを持ちます。
3.1 恒星とガスの総質量
M_star \;=;2pi, R_d^2, ˶Sigma_d˶, ˶Upsilon_star
M_text{gas}\M_text {HI}.\(He補正; Arnett 1996)}$.
3.2 コンポーネントの質量とスケール
| コンポーネント | 質量 | スケール | 活性化 |
|---|---|---|---|
| バルジ | M_b = 0.20, M_star$ | r_b = \max(0.5,R_d,Γ, 0.3text{ kpc})$. | もし $T |
| 薄型ディスク | M_text{thin} = 0.75,(M_static – M_b)$ | R_d$ | 常に |
| 厚いディスク | M_text{thick} = 0.25*,(M_static – M_b)$ | 1.5ドル | 常に |
| ガスリング | M_text{gas} = 1.33, M_text{HI}$ | R_g = 1.7，R_d$ (Broeils & Rhee 1997) | 常に |
| スパイラルアーム | M_text{arm} = 0.10、M_text{thin}$（有効） | R_d$ (薄い円盤に従う) | 常に |
3.3 密度プロファイル
バルジ (3D Hernquist)
rho_b(r)
薄い恒星円盤と厚い恒星円盤(2次元指数)
$$\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)}$$
ガスリング(中心穴のある2次元指数)
Sigma_text{gas}(R) ￤$$Sigma;=;￤frac{M_text{gas}}{2pi,R_g^2}}, ￤exprac!￤left(-￤R_text{hole}}{R} -￤frac{R}{R_g}￤right), ￤quad R_text{hole} = 0.5,R_g$.
螺旋腕過剰（2次元、薄い円盤に従う）
Sigma_text{arm}(R) Γ;=Γ;0.10;ΓSigma_text{thin}(R)$$.
4.波動カーネル
各バリオン質量要素はビー理論の波動場を生成します。D=｜vec{r}-｜vec{r},’|$で区切られた$vec{r},’$にある光源要素によって生成される点$vec{r}$における場は、注Ⅰの正則化された波動関数から導かれる湯川型カーネルによって支配されます：
ビー理論の波動カーネル
K_0, ￤frac{(1 + ￤αi)￤e^{-￤αi)￤}{D^2}, ￤qquad ￤frac{1}{ell_i}$.
ここで$K_0$は普遍的な波動質量振幅（無次元数）、$iell_i$は成分$i$のコヒーレンス長。K_0$は普遍的な波動質量振幅(無次元数)、$i$は成分$i$のコヒーレンス長。カーネルは短い距離では準ニュートン的な$1/D^2$の振る舞いをし、$i$を超えるスケールでは指数関数的なカットオフで変調されます。(1 + Γ D)Γ,e^{-Γ D}$の形は連続性を保証し、無限大での全封入質量は有限。
4.1 コンポーネントのコヒーレンス長
各成分のコヒーレンス長は、その自然な幾何学的スケールに、その次元に特有の無次元定数を掛けたものです：
| コンポーネント | コヒーレンス長 | 幾何定数 |
|---|---|---|
| バルジ（3D球体） | well_b = c_text{sph}, r_b$. | c_text{sph}$ |
| 薄型ディスク（2D） | c_ell_text{thin} = c_text{disk}, R_d$ | c_text{disk}$ |
| 厚いディスク（2D） | cell_text{thick} = c_text{disk},(1.5, R_d)$ | c_text{disk}$ |
| ガスリング（2D） | l_ell_text{gas} = c_text{disk}}, R_g$. | c_text{disk}$ |
| スパイラルアーム（2D、方位角集中型） | arm = c_text_arm | c_text{arm}$ |
つの幾何学定数$(c_text{sph},Γ,c_text{disk},Γ,c_text{arm})$は普遍的なもので、銀河によって変わることはありません。大域的な波動質量振幅$K_0$と波動場結合$lambda$と合わせて、理論レベルのパラメータ一式を構成します。
5.波動場の畳み込み – 成分ごとの積分方程式
ある位置$vec{r}$での波動場密度は、バリオン源分布と波動カーネルの畳み込みです。銀河対称系（軸対称、単極近似）では、各バリオン成分は加法的に寄与します：
半径$r$における全波動場密度
sum_{i\rho_text_wave}^{(i)}(r)$$。
以下に5つの積分を，各成分ごとに1つずつ示します。各積分は，同じ空間点におけるバリオン質量分布を波動場質量分布に変換します。
5.1 バルジ – 3Dシェルの統合
波動}^{(b)}(r) \;=; ￤int_0^{r_text{max}} ￤int_0^{r_text{max\rho_b(r’)￤(￤sqrt{r^2 + r’^2}￤right)￤pi r’^2,dr’$$
積分は半径$r’$の同心球殻上で行います。中心から半径$r$にある磁場点は、単極近似で$D = \sqrt{r^2 + r’^2}$の有効分離で各殻を見ます。この積分は$r_text{max} = 6まで拡張され、それ以上ではバルジ密度は数値的に無視できます。
5.2 薄い円盤と2Dリングの統合
$$\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)Γ;2pi R’Γ,dR’$$
円盤は半径$R’$、幅$dR’$の同心円に分解され、各リングは表面質量$Sigma_text{thin}(R’)。同じ単極近似が適用され、中心から半径$r$の波動場は、実効的な分離$D = \sqrt{r^2 + R’^2}$で各リングからの寄与を受けます。積分範囲は$R_text{max} = 8,R_d$ 。
5.3 厚いディスクと2Dリングの統合
$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}}\Sigma_\text{thick}(R’)\;\mathcal{K}_\text{thick}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
薄い円盤の積分と同じで、$Sigma_text{thick}(R’)$を線源密度、カーネルパラメータ$alpha_text{thick} = 1/(c_text{disk},\cdot 1.5,R_d)$ 。厚い円盤の半径方向の広がりが大きいため、波動コヒーレンスの範囲が若干広い。
5.4 ガスリング – 中心空乏を伴う2Dリング統合
$$\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’$$
ガス分布は中心に穴があり、積分の下限$R_text{hole} = 0.5,R_g$の半径方向のカットオフでとらえられます。このカットオフの外側では、ガスは恒星円盤よりも広がっています。これは、コヒーレンス長$ell_text{gas} = c_text{disk},R_g$に反映される大きな特性スケール$R_g = 1.7,R_d$に反映されます。
5.5 スパイラルアームの過剰 – 振幅を減らした2Dリング積分
$$\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’$$
渦巻きの腕は、薄い円盤の表面密度の軸方向に平均した増強として扱われ、そのコヒーレンス長は$ell_text{arm} = c_text{arm}、R_d$。したがって、このカーネルは薄い円盤のカーネルよりも狭く、らせん構造の方位角集中を反映しています。
6.封じ込められた波の質量と予測回転曲線
全波動場密度$rho_text{wave}(r)$が分かれば、半径$R$の球内の囲まれた波動場の質量を半径方向に積分して求めます：
閉じた波動場の質量
M_text{wave}(R) ˶;=; ˶int_0^{R} 4pi,r^2 ˶rho_text{wave}(r)˶,dr$
半径$R$での円周速度の予測は、バリオン寄与と波動場寄与を直交させてニュートンの関係から導かれます：
予測円速
V_c^2(R) Γ;=; V_text{bar}^2(R) Γ;+; Γfrac{G,M_text{wave}(R)}{R}$.
バリオン速度$V_text{bar}(R)$は，4つの円盤状成分（Freeman 1970の指数分布の公式）とバルジ（Hernquistの囲み質量の公式）からの寄与の2次和です：
V_text{bar}^2(R) \;=V_text{bulge}^2 + V_text{thin}^2 + V_text{thick}^2 + V_text{gas}^2$.
ここで、各$V_i(R)$は、対応する質量分布の標準ニュートン円速度です。
7.理論レベルのパラメータ
銀河に適用される完全なBeeTheoryのフレームワークには、理論レベルで5つのパラメータが含まれています。これらは普遍的なもので、銀河によって変わることはありません。
| シンボル | 意味 | 役割 |
|---|---|---|
| $K_0$ | 波質量の振幅 | 波動カーネルの無次元スケールを設定します。 |
| c_text{sph}$ | 3次元幾何学定数 | 球面光源(bulge)の比率 $ell/r_text |
| c_text{disk}$ | 2次元幾何定数 | 円盤光源とリング光源の比率$ell/R_text{scale}$。 |
| c_text{arm}$ | スパイラル幾何学定数 | Ratio $ell/R_d$ for azimuthally concentrated arm excess (方位角集中腕過剰の比率$ell/R_d$) |
| $lambda$ | グローバル波動場結合 | 全波動場密度のスケール |
パラメータの普遍性
5つのパラメータはすべてグローバルです。天の川銀河にも、矮小不規則銀河にも、大質量渦巻き銀河にも同じ数値が適用されます。銀河固有の情報は、5つの観測入力$(T,∕R_d,∕Sigma_d,∕M_text{HI},∕Upsilon_star)$ を通してのみ入ります。このモデルには銀河ごとに調整可能なパラメータは含まれていません。
8.計算の一方向性
オープンチェーン – フィードバックなし
計算全体は、入力から出力へと一方向に流れます。測光と21cmの観測がバリオン分解を決定します。バリオン分解が波動場密度を決定。波動場密度は包含波質量を決定します。波動質量が回転曲線の予測値を決定。自転曲線は、それ以前のどのステップの計算にも影響を与えません。
この一方向性は、3つの重要な結果をもたらします。
(a)5つの理論レベルのパラメータが固定されると、回転曲線はフィットではなく厳密な予測になります。観測された回転曲線との比較はテストであり、校正ではありません。
(b)このモデルには銀河ごとの調整のメカニズムがありません。回転曲線の予測値を修正するには、入力ベクトル$(T,˶,R_d,˶,Sigma_d,˶,M_text{HI},˶,Upsilon_star)$を修正するか、普遍理論レベルのパラメータ$(K_0,˶,c_text{sph},˶,c_text{disk},˶,c_text{arm},˶,lambda)$を変更する必要があります。
(c)$lambda$を参照銀河で較正することは、その銀河の自転曲線に当てはめることとは違います。基準銀河の他の半径での自転曲線や、他のすべての銀河の自転曲線は、校正された枠組みの厳密な予測になります。
9.中心面密度の役割（注XI改訂版）
残差予測誤差は、円盤のスケール長$R_d$とは無関係に、中心バリオン表面密度$Sigma_d$と強い相関があることが、注XIの診断で明らかになりました。上で紹介した公式化は、この発見を取り入れる前のモデルのバージョンで、コヒーレンス長の式に$R_d$だけを使っています$ell_i = c_i,R_d$.
リファインメントが入るところ
精密化モデルでは、コヒーレンス長$ell_i$は$R_d$と$Sigma_d$の両方に依存し、厳密な線形関係$ell_i = c_i,R_d$を、注XIで同定した残差を吸収する関数$ell_i = c_i,R_d,Γphi(ΓSigma_d/ΓSigma_text{ref})$に置き換えます。の関数形とそのパラメータは、次のノートで、まず22銀河の較正セットで決定し、次に残りのSPARCサンプルでブラインド予測により検証します。
Sigma_d$は観測入力であり、修正されたコヒーレンス長は同じ畳み込み積分に入り、回転曲線は以前と同じように現れます。操作上のリンクが一つ追加されただけで、$ell_i$の第二観測量への依存性です。
10.方法論のまとめ
1.入力。ハッブル型$T$、円盤スケール$R_d$、表面輝度$Sigma_d$、HI質量$M_text{HI}$、恒星質量光量比$Upsilon_star$。
2.バリオン分解。バルジ（$T$4$の場合）、薄い円盤、厚い円盤、ガス環、渦状腕過剰の5成分。それぞれ解析的密度分布。
3.波動カーネル。Universal Yukawa-type form $mathcal{K}_i(D) = K_0,(1 + \alpha_i D) \,e^{-αalpha_i D}/D^2$ with coherence length $ell_i = c_i,R_text{scale}$ （各成分の幾何学的広がりで決まる）。
4.畳み込み。各成分は、リング(2次元成分)またはシェル(3次元バルジ)上の1次元積分によって波動場密度を生成します。全波動場密度は、大域結合$lambda$でスケーリングされた5つの成分の和です。
5.出力。内包された波の質量$M_text{wave}(R)$を積分し、バリオン速度$V_text{bar}(R)$と組み合わせることで、予測される回転曲線$V_c(R)$が得られます。
6.理論レベルのパラメータ(K_0,Γ,c_text{sph},Γ,c_text{disk},Γ,c_text{arm},Γ,Γlambda)$ – universal, no per-galaxy tuning.検討中の改良で$Sigma_d$依存性が追加される予定。
7.方向。入力→バリオン→波動場→回転曲線。フィードバックなし。回転曲線はフィットではなく予測。
参考文献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. –渦巻銀河とS0銀河の円盤について, ApJ 160, 811 (1970).- Hernquist, L. –An analytical model for spherical galaxies and bulges, ApJ 356, 359 (1990).- ブロイユ, A. H., Rhee, M.-H.– 渦巻銀河と不規則銀河の 21cm WSRT による短波長観測, A&A 324, 877 (1997).- McGaugh, S. S. –The third law of galactic rotation, Galaxies 2, 601 (2014).- Bovy, J., Rix, H.-W. – 天の川銀河の円盤表面密度分布の動的直接測定, ApJ 779, 115 (2013).- Arnett, D. –Supernovae and Nucleosynthesis, Princeton (1996).- Dutertre, X. –Bee Theory™：Wave-BasedModeling of Gravity, v2, BeeTheory.com (2023)。
Bee理論.com – 波動量子重力 – 銀河の方法論 – © Technoplane S.A.S. 2026