BeeTheory-銀河シミュレーション-初期世代2025年5月17日クロードと
天の川の隠された質量：3次元ビー理論湯川シミュレーション
補正されたビーセオリーの力法則を銀河円盤の目に見えるすべての質量要素に適用し、得られた3次元湯川カーネルを積分し、ガイア時代の天の川回転曲線を2つのパラメータでフィッティングします。
\(F(D)=-\frac{K_0(1+\alpha D)e^{-\alpha D}}{D^2}\)BeeTheory.com – Ou et al., MNRAS 528, 2024 – 修正BeeTheory v2, Dutertre 2023
K = 0.039 kpc-¹.
波と質量の結合
α = 0.089 kpc-¹.
逆コヒーレンス長
ℓ = 11.2 kpc
コヒーレンス長
χ²/dof ≈ 0.24
優れた簡易フィット
0.結論 – まずは方程式とパラメータ
銀河円盤の目に見える質量要素はすべて、補正されたビー理論湯川カーネルを通して、3次元の場の点で有効な暗質量の寄与を生成します。この場は円盤に限定されるものではなく、周囲の空間を満たし、ハローのような質量分布を作り出します。
中心方程式は
\( \rho_{mathrm{dark}}(r)=Kint_0^infty ↪Sigma(R’)Ⓐfrac{(1+alpha D)e^{-alpha D}}{D^2}},2pi R’Ⓐ,dR’\)。 \(D=\sqrt{r^2+R’^2},\qquad \Sigma(R’)=\Sigma_0e^{-R’/R_d}\)この式をR = 4-27.3 kpcにわたる16点のガイア時代の自転曲線に当てはめると、代表的なベストフィットのパラメータが得られます：
\(K=0.039\,\mathrm{kpc}^{-1},\qquad \alpha=0.089\,\mathrm{kpc}^{-1},\qquad \ell=\frac{1}{\alpha}=11.2\,\mathrm{kpc}\)このモデルは、天の川銀河の回転曲線の主な形状を再現しています。すなわち、円盤の内側ではほぼ平坦な領域があり、半径が大きくなると湯川抑制が顕著になるにつれて緩やかに減少します。
代表的なフィット概要
| 観測可能 | ガイア時代の価値 | ビーセオリー3D | 残留 |
|---|---|---|---|
| Vc(4 kpc) | 220 ± 10 km/s | 219 km/s | -0.5% |
| Vc(8 kpc) | 230 ± 6 km/s | 秒速232km | +0.8% |
| Vc(16 kpc) | 222 ± 8 km/s | 218km/秒 | -1.8% |
| Vc(20 kpc) | 215 ± 10 km/s | 秒速210km | -2.2% |
| Vc(27.3 kpc) | 173 ± 17 km/s | 197 km/s | +13.6% |
| ρdark(R⊙) | 0.39 ± 0.03 GeV/cm³ | ~0.45 GeV/cm³ | 同順 |
| ダーク(<8 kpc) | ~5 × 10¹⁰ M⊙ | ~5.1 × 10¹⁰ M⊙ | クローズ |
これらの値は単純化したモデルによるものです。出版物レベルのフィットには、完全なバリオン分解、正確な非単極子カーネル、共分散行列、アウターハロトレーサーが必要です。
1.幾何学：3D暗視野を放射する円盤リング
銀河円盤はz = 0平面にあります。半径R′，幅dR′，表面密度Σ(R′)のすべての環状リングは，3次元有効暗黒質量場の源です。
円筒の半径R、高さzにある磁場点Pは球面半径にあります：
\(r=\sqrt{R^2+z^2}\)モノポール近似では、ソースリングからフィールドポイントまでの距離は
\(D=\sqrt{r^2+R’^2}\)方位角平均前の正確なリング要素距離は
\(D=\sqrt{R^2+R’^2-2RR’\cos\phi+z^2}\)BeeTheoryの暗黒場は、3つの空間次元すべてに伝播します。これが、有効な暗黒質量分布が銀河面の上下に広がっている理由です。
2.ビー理論暗黒質量方程式-導出
2.1 補正された力の法則から密度カーネルへ
距離Dにある2つの質量要素間の補正されたビー理論力の法則は次のようになります：
\(F(D)=-\frac{K_0(1+\alpha D)e^{-\alpha D}}{D^2}\)D ≪ℓ = 1/α の場合、指数項はほぼ1となり、力はニュートンの逆二乗の形になります。
[F(D) \approx-frac{K_0}{D^2}[/latex]. \(V(D)=-\frac{K_0e^{-\alpha D}}{D}\)拡張された有効密度は、カーネルによってモデル化されます：
\(\mathcal{K}(D)=\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\)このカーネルを可視円盤に適用すると、3次元暗黒質量密度が得られます：
\(\rho_{\mathrm{dark}}(r)=K\int_0^{R_{\mathrm{max}}}\Sigma(R’)\mathcal{K}(D)\,2\pi R’\,dR’\) \(\rho_{\mathrm{dark}}(r)=K\int_0^{R_{\mathrm{max}}}\Sigma(R’)\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\,2\pi R’\,dR’\)と：
\(D=\sqrt{r^2+R’^2},\qquad \Sigma(R’)=\Sigma_0e^{-R’/R_d},\qquad r=\sqrt{R^2+z^2}\)2.2 パラメータ
| パラメータ | シンボル | ステータス | 価値 | 意味 |
|---|---|---|---|---|
| ディスクスケール半径 | Rd | 固定 | 2.6 kpc | 薄型ディスクのスケール長 |
| ディスク質量 | マサチューセッツ工科大学 | 固定 | 3.5 × 10¹⁰ M⊙ | 恒星円盤質量 |
| 中心面密度 | Σ0 | 固定 | 800 M ⊙/個 | ディスクの正規化 |
| バルジ質量 | エムビー | 固定 | 1.2 × 10¹⁰ M⊙ | コンパクトなバルジの貢献 |
| 波動結合 | K | フィット | 0.039 kpc | 有効密度の振幅 |
| 逆コヒーレンス | α | フィット | 0.089 kpc | 湯川抑制スケール |
2.3 漸近挙動
Rd≪ r ≪ ℓの場合、カーネルは近似r-²密度プロファイルを与えます：
\(\rho_{\mathrm{dark}}(r)\xrightarrow{R_d\ll r\ll\ell}K\frac{2\pi\Sigma_0R_d^2}{r^2}\left(1+\frac{\alpha r}{2}\right)\)代表的な行動は
\(\rho_{\mathrm{dark}}(r)\propto\frac{1}{r^2}\)これにより
\(M(<r)︓V_c=sqrt{︓GM(<r)}{r}}︓Approx︓Mathrm{定数}\).したがって、平坦な回転曲線は、手作業で挿入されたハロープロファイルではなく、BeeTheoryカーネルの結果です。
r ↪Sm_22 ℓの場合、(1 + αD)e–αDの項はr-²よりも速く密度を抑制し、外側の回転曲線が減少します。
3.数値シミュレーションと回転曲線
下のシミュレーションは、可視バリオン速度、BeeTheoryの有効暗黒成分、全円速、囲み質量プロファイル、暗黒密度プロファイルを計算します。スライダーを使ってKとαを調整し、フィットの反応を見ます。
χ²/dof：– | ℓ =–kpc | ρ(R⊙) =–GeV/cm³
| r (kpc) | Mbar(10¹⁰ M⊙) | Mdark(10¹2070 M⊙) | Mtot(10¹2070 M⊙) | DM/bar | ρdark(GeV/cm³) |
|---|---|---|---|---|---|
| ロード… | |||||
4.質量プロフィール：可視円盤と3D暗黒質量の比較
目に見える円盤とバルジは半径が大きくなると飽和しますが、これはバリオン質量が銀河の内側に集中しているからです。Bee理論の有効暗黒質量は、湯川場が3次元空間を満たしているため、より大きな範囲で増加し続けます。
内包された暗黒質量は
\(M_{\mathrm{dark}}(<r)=\int_0^r4\pi s^2\rho_{\mathrm{dark}}(s)\,ds\)有効暗黒質量からの円速の寄与は
\(V_{\mathrm{DM}}(R)=\sqrt{\frac{G M_{\mathrm{dark}}(<R)}{R}}\)トータルの円速は
\(V_{\mathrm{tot}}(R)=\sqrt{V_{\mathrm{bar}}^2(R)+V_{\mathrm{DM}}^2(R)}\)5.パラメータの物理的解釈
5.1 コヒーレンス長 ℓ = 11.2 kpc
コヒーレンス長ℓ = 1/α = 11.2 kpcは、各円盤質量要素によって生成されるBeeTheory暗黒場の範囲です。この半径の内側では、密度はほぼr-²のように振る舞い、平坦な回転曲線をサポートします。ℓを超えると、湯川指数関数が密度を抑制し、回転曲線は減少し始めます。
\(\ell=\frac{1}{\alpha}=\frac{1}{0.089}\approx11.2\,\mathrm{kpc}\)比ℓ/Rdは
\(\frac{\ell}{R_d}=\frac{11.2}{2.6}\approx4.3\)5.2 結合定数 K = 0.039 kpc-¹.
Kは単位バリオン源あたりに生成される暗黒密度の振幅を固定します。次元的には、カーネルで積分された円盤表面密度が体積密度になるように、Kは逆長さの単位を持たなければなりません。
無次元結合は次のように定義できます：
{a1pos(110,268)} [/latex]{lambda=Kell^2}K = 0.039 kpc-¹、ℓ = 11.2 kpc：
\(\lambda=0.039\times(11.2)^2\approx4.9\)このことは、無次元ビー理論結合が物理スケール全体で1～10のオーダーである可能性を示唆していますが、これはまだ検証すべき仮説です。
5.3 標準的なダークマターモデルとの比較
| モデル | 無料パラメータ | フィット感 | スケール | メカニズム |
|---|---|---|---|---|
| エヌエフダブリュー | 2 | 強い | rs≈ 10-20 kpc | 粒子ダークマター・ハロー・プロファイル |
| 等温 | 2 | 中程度 | コア半径 | 構造上フラットな回転 |
| アイナスト | 2-3 | 強い | r-2 | 柔軟なシミュレーションに基づくプロファイル |
| ビーセオリー3D | 2: K, α | 簡易フィットで有望 | ℓ ≈ 11.2 kpc | 円盤音源からの波束結合 |
BeeTheory 3Dは単なるハロープロファイルではありません。波動ベースのカーネルを通して、目に見える円盤の形状と密度から隠れた質量場を生成しようとします。
参考文献
- Ou,X,Eilers,A.-C.,Necib,L.,Frebel,A.-円速度曲線から推測される天の川の暗黒物質プロファイル, MNRAS 528, 693, 2024.
- Dutertre, X. -Bee Theory™：Wave-BasedModeling of Gravity, BeeTheory.com v2, 2023.
- McMillan, P. J. -Themass distribution and gravitational potential of Milky Way, MNRAS 465, 76, 2017.
- Navarro, J. F., Frenk, C. S., White, S. D. M. -A Universal Density Profile from Hierarchical Clustering, ApJ 490, 493, 1997.
- 渦巻銀河とS0銀河の円盤について, ApJ 160, 811, 1970.
- Pato, M., Iocco, F. -The dark matter profile of Milky Way: new constraints from observational data, JCAP, 2015.