ビー理論暗黒質量モデルの数値シミュレーション

BeeTheory暗黒質量シミュレーションの背後にある数値積分、バリオン速度分解、χ²最小化、実装の選択についての完全で再現可能な説明。

このテクニカルページでは、BeeTheoryによる天の川の隠れた質量の数値シミュレーションの再現方法を説明します。観測データ、バリオンモデル、波動ベースの密度方程式、数値積分、フィッティング方法、収束テスト、リファレンスコードについて説明しています。

目的は簡単で、目に見える天の川円盤から始め、BeeTheoryの波動密度モデルを適用し、有効な隠れた質量を計算し、その結果得られる円速曲線をガイア時代の観測結果と比較することです。

0.私たちは何を計算し、なぜ計算するのか

天の川の自転曲線は、銀河中心からの距離Rの関数としての星の円速度Vc(R)です。私たちが直接見ることのできる全質量分布よりも、はるかに正確に測定されています。

観測された速度と、目に見えるバリオン物質が予測する速度との間にある欠損が、隠れた質量の問題です。標準モデルでは、目に見えない粒子ハロー（通常は冷たい暗黒物質）が存在します。BeeTheoryは別の解釈を提案します。目に見える質量の要素はすべて、3次元で指数関数的に減衰する波動場を放射しており、蓄積された場は隠れた質量のように振る舞います。

シミュレーションは3つのことを行います：

1. Freemanの解析的円盤公式を用いて、可視円盤とバルジから^Vcbar(R)を計算します。
2. ビーセオリーの密度ρ_dark(r; ℓ, λ)を各半径で数値積分し、^VcDM(R)に変換します。
3. ガイア時代の天の川回転曲線に対してχ²(ℓ, λ)を最小化し、代表的なベストフィットのパラメータを見つけます。

シミュレーションは再現できるように設計されています。特別な天体物理学ライブラリがなくても、JavaScriptやPythonで実行できます。

表記方法

• Rは円盤面における円筒状の銀河中心半径で、単位はkpc。
• rは球形の銀河中心半径。
• zは円盤からの高さ。
• ℓはBeeTheoryのコヒーレンス長で、単位はkpc。
• λは無次元波動質量結合。
• Gはkpc km² s-² M⊙-¹の単位で使われています。

1.観測データ – ガイア2024

データセットはガイア時代の天の川回転曲線に基づいています。半径R、円速度_Vc、不確かさσを使用しています。オリジナルの技術版では、R = 4-27.3 kpcをカバーする16個のデータが使われています。

重要な観測事実は

• _Vc(R⊙=8kpc)≒230km/s、太陽の軌道に近いアンカーポイント。
• _Vcは約5〜20kpcでほぼ平坦。
• _Vcは円盤の外側で減少し、27.3kpc付近で173±17km/sに達します。

このためには、中間半径で強く上昇し、半径が大きくなると弱くなる有効な隠れた質量分布が必要です。

銀河系の一番内側は、中心棒と非円運動が支配的なので除外します。その領域では、単純化した軸対称モデルは信頼できません。

3.ビー理論の暗黒密度方程式

3.1 物理的仮定

銀河円盤の位置r′にあるすべての可視質量要素dVは、重力波磁場を生成し、磁場点rに有効質量密度を作ります：

任意の点(R,z)における全暗黒密度は、すべての円盤要素の重ね合わせです。光源は円盤なので、体積密度の項は表面密度の項になります：

完全な二重積分は

3.3 漸近極限の解析的検証

_Rd≪ r ≪ ℓ の場合、式は単純化されます。円盤は半径よりずっと小さく、コヒーレンス長はやはり半径よりずっと大きい。

だって：

この漸近密度はρ ∝ r-²として振る舞い、M( [latx] \rho(r)｟propto r^{-2}Longrightarrow M(<r)｟propto r Longrightarrow V_c=sqrt{prac{GM(<r)}{r}} approx\mathrm{constant}[/latex].

代表的な値を差し込むと、観測された天の川の値に近い局所密度が得られます：

5. χ² 最小化とパラメータ・フィッティング

5.1 目的関数

フィットはカイ2乗を最小化します：

N = 16データポイント、p = 2自由パラメータ、ℓとλ。

5.2 2パスグリッド探索

ランドスケープがλとℓの間に長くカーブした縮退の谷を持つため、勾配降下ではなくグリッド探索が用いられます。

• パス1:ℓとλの粗いグリッド。
• パス2:粗い最小値を中心とした局所的な絞り込み。

代表的なベストフィット領域は

5.3 退廃の谷

χ²ランドスケープは対称的なボウルではありません。これは、先行密度正規化が結合強度に強く依存し、平坦回転領域ではコヒーレンス長に弱く依存するだけだからです。

回転曲線の外側の減少は、ℓ値が小さいほど実効密度が早く抑えられるため、この縮退を破ります。

球状星団、恒星ストリーム、ハロー星、衛星銀河など、30kpcを超えるデータは、ℓをより厳密に制約するために不可欠です。

7.完全参照コード

以下のJavaScriptリファレンス実装は、主なシミュレーションロジックを再現しています。これは技術的な検証を目的としており、カスタムスクリプトが許可されていない限り、標準的なWordPressのコンテンツブロック内に直接置かず、適切なスクリプト環境に置く必要があります。

// 物理定数
定数 G = 4.302e-3; // kpc km² s-² M☉-¹.
const Sig0 = 800.0; // M☉ pc-².
const Rd = 2.6; // kpc
const Mdisk = 3.5e10; // M☉.
const Mbulge = 1.2e10; // M☉.
const CONV_RHO = (1.989e30 * 5.61e26) / (3.086e21)**3；

// ガイア時代の回転データ
const OBS_R = [4,5,6,7,8,9,10,11,12,14,16,18,20,22,24,27.3]；
const OBS_V = [220,228,232,231,230,229,228,227,226,224,222,219,215,208,200,173]；
const OBS_ERR = [10,8,7,7,6,6,7,8,9,10,11,13,17]；

// バリオン速度プレースホルダ
関数 vBaryonic(R) {
  // 完全な実装では、Freemanの円盤の公式を使用します。
  // 修正Bessel関数I0, I1, K0, K1を使用します.
  const vb2 = G * Mbulge / Math.max(R, 0.2)；
  return Math.sqrt(Math.max(0, vb2))；
}

// ビー理論の暗黒密度
関数 rhoDark(r, ell, lam) { // ビー理論の暗黒密度
  const N = 40；
  const dRp = 30.0 / N；
  let sum = 0；

  for (let i = 0; i < N; i++) { .
    const Rp = (i + 0.5) * dRp；
    const Sig = Sig0 * Math.exp(-Rp / Rd)；
    const K = (2 * Math.PI * ell / r)
      * Math.sinh(Math.min(r / ell, 30))
      * Math.exp(-Math.min((r + Rp) / ell, 30))；

    sum += Sig * Rp * K * dRp；
  }

  return (lam / ell) * sum；
}

// 囲みダークマス
関数 enclosedDarkMass(R, ell, lam) { // 密閉されたダークマス
  const Nr = 30；
  const dr = R / Nr；
  M = 0；

  for (let j = 0; j < Nr; j++) { .
    const rj = (j + 0.5) * dr；
    M += 4 * Math.PI * rj * rhoDark(rj, ell, lam) * dr；
  }

  M を返します；
}

// 暗黒速度と全速度
関数vDM(R, ell, lam) { // 暗黒・全速度
  return Math.sqrt(Math.max(0, G * enclosedDarkMass(R, ell, lam) / R))；
}

関数vTotal(R, ell, lam) { } }.
  const vb = vBaryonic(R)；
  const vd = vDM(R, ell, lam)；
  return Math.sqrt(vb * vb + vd * vd)；
}

// カイ二乗
関数 chiSquared(ell, lam) { } // カイ二乗
  let s = 0；

  for (let i = 0; i < OBS_R.length; i++) { .
    const dv = (vTotal(OBS_R[i], ell, lam) - OBS_V[i]) / OBS_ERR[i]；
    s += dv * dv；
  }

  return s / (OBS_R.length - 2)；
}

完全版には、フリーマン円盤速度に対する修正ベッセル関数I₀、I₁、_K0、_K1の正確な実装が含まれているはずです。

参考文献

• アブラモウィッツ、M.、ステガン、I.A.-数学関数ハンドブック、ドーバー、1972年。
• Bovy, J., Rix, H.-W.-A Direct Dynamical Measurement of theMilky Way'sDisk Surface Density Profile, ApJ 779, 115, 2013.
• 渦巻銀河とS0銀河の円盤について, ApJ 160, 811, 1970.
• McMillan, P. J. -Themass distribution and gravitational potential of Milky Way, MNRAS 465, 76, 2017.
• Ou,X,Eilers,A.-C.,Necib,L.,Frebel,A.-円速度曲線から推測される天の川の暗黒物質プロファイル, MNRAS 528, 693, 2024.
• Pato, M., Iocco, F., Bertone, G. -Dynamical constraints on dark matter distribution in Milky Way, JCAP 12, 001, 2015.
• Portail, M. et al. -Dynamical modelling of Galactic bulge and bar, MNRAS 465, 1621, 2017.