ダークマターと銀河回転曲線問題の数学

Q: ダークマターと銀河回転曲線問題の数学は社会やプログラミングでどのように使われていますか？

NFWプロファイルによる銀河回転曲線のフィット計算コード: 宇宙論シミュレーションで標準的に用いられるNavarro-Frenk-White (NFW) 暗黒物質ハロープロファイルから銀河回転速度曲線を計算するPythonコードです。

外側の星が遅くならない謎 — ニュートン力学の破綻と未知の暗黒物質ハロー

この単元の発見をシェア

3行要約（Key Takeaways / AIOサマリー）

1: 1970年代にヴェラ・ルービンが渦巻銀河の観測から発見した、現代宇宙論最大の謎『銀河回転曲線問題』
2: ニュートン力学とケプラーの法則によれば、銀河の外縁部の星は中心から離れるほど $v \propto 1/\sqrt{r}$ で減速するはずだが、実際には外側も猛スピードのまま平坦（ $v \approx \text{const}$ ）に回転していた
3: 目に見える星やガスの引力だけでは銀河が遠心力で四散してしまうため、光を発しない未知の質量『ダークマター（暗黒物質）』が銀河全体を球状に包み込んでいる（銀河総質量の約85%を占める）ことが数学的に確定した

なぜこれを学ぶのか？ — 現実社会とテクノロジーでの価値

Purpose & Motivation

なぜ生まれたのか？

『数式による理論予想と、現実の観測データが決定的に食い違ったとき、科学はどう進歩するのか』という物理学史の最高のドラマを体験できます。万有引力と等速円運動の向心力 $v = \sqrt{GM(r)/r}$ という基本公式から、宇宙の85%が未知の物質であるという驚愕の結論に至る論理的道筋を学べます。

身につく武器・思考力

球対称質量分布の積分（ガウスの法則・ニュートンの球殻定理）の直感理解、質量プロファイル $M(r) \propto r$ のモデリング力、科学的データ検証の思考法

実社会・ITでの活躍現場

•すばる望遠鏡のHSC（超広視野主焦点カメラ）による重力レンズ効果を用いた全天ダークマター地図作成
•XMASSやXENONnT実験など、地下深部検出器によるダークマター直接検知（WIMP粒子探索）
•銀河団衝突（弾丸銀河団）のX線高温ガスと重力ポテンシャルの分離解析
•超巨大宇宙シミュレーション（IllustrisTNG等）における銀河形成・フィラメント構造の進化計算

Layer 1 — 直感・シミュレーター

『観測現実（暗黒物質ハローあり）』と『理論予想（目に見える星のみ）』を切り替えてみてください。理論予想では外側の星がケプラー減速して極端に遅くなり、右側のグラフの点線に沿って落ちていきます。観測現実では青紫色の暗黒物質ハローの重力によって、外縁部までフラットな高速度を保ち続ける様子が確認できます。

外縁部も高速回転（平坦曲線）

ヴェラ・ルービンの大発見：平坦な銀河回転曲線

暗黒物質（ダークマター）の推定質量比率85% (通常約85%)

直感チェック！ミニクイズ

銀河の回転曲線が外縁部まで平坦（v(r) = 一定）であるとき、内包質量 M(r) と距離 r の関係はどうなるか？

Layer 2 — 定義・受験・証明

定義

**銀河回転曲線問題（Galaxy Rotation Problem）**とは、円盤銀河における恒星や星間ガスの円軌道公転速度 $v(r)$ が、中心からの距離 $r$ が大きくなってもニュートン力学の予想（ $v \propto r^{-1/2}$ ）のように低下せず、ほぼ平坦（ $v(r) \approx \text{const}$ ）な値を維持する観測的事実である。

公式・数理モデル

\text{等速円運動の運動方程式: } \frac{v(r)^2}{r} = \frac{G M(r)}{r^2} \implies v(r) = \sqrt{\frac{G M(r)}{r}} \text{バリオン（可視光）物質のみ: } M(r) \approx M_{\text{core}} = \text{一定} \implies v(r) \propto \frac{1}{\sqrt{r}} \text{観測事実: } v(r) = v_0 (\text{一定}) \implies M(r) = \frac{v_0^2}{G} r \implies \rho_{\text{DM}}(r) \propto \frac{1}{r^2}

証明・導出

**平坦な回転曲線から導かれるダークマター密度分布（球殻定理）** 半径 $r$ の軌道を回る星に対して、ニュートンの球殻定理より、星より内側にある質量 $M(r)$ のみが原点からの中心力として作用する：

\frac{m v^2}{r} = G \frac{M(r) m}{r^2} \implies v(r)^2 = \frac{G M(r)}{r}

もし質量が目に見える星（銀河中心部に密集）だけであれば、外側（ $r > R_{\text{disk}}$ ）では $M(r) \approx M_{\text{total}}$ となり、速度は：

v(r) = \sqrt{\frac{G M_{\text{total}}}{r}} \propto r^{-1/2}

と急激に減少するはずである。

しかし、観測事実は外縁部まで $v(r) = v_0$ （一定）である。これを方程式に代入すると：

v_0^2 = \frac{G M(r)}{r} \implies M(r) = \frac{v_0^2}{G} r

質量 $M(r)$ が中心からの距離 $r$ に完全に比例して無限に増え続けなければならない。球対称密度分布 $\rho(r)$ と内包質量の関係 $M(r) = \int_0^r 4\pi r'^2 \rho(r') dr'$ を両辺 $r$ で微分すると：

\frac{dM}{dr} = 4\pi r^2 \rho(r) = \frac{v_0^2}{G} \implies \rho(r) = \frac{v_0^2}{4\pi G} \frac{1}{r^2}

すなわち、銀河の外側には可視光では見えない質量が $\rho \propto 1/r^2$ の密度で広大に分布している（等温ダークマターハロー）ことが数学的に厳密に導かれる（Q.E.D.）。

押さえておきたいポイント

可視光で見える星やガスの質量は銀河全体のわずか15%程度にすぎず、残りの約85%は未解明のダークマターである
暗黒物質の存在を仮定しない別の解決策として、長距離で重力法則が修正されるとする『修正ニュートン力学（MOND）』も提唱されている
重力レンズ効果や宇宙背景放射（CMB）の観測からも、冷たい暗黒物質（CDM）モデルが圧倒的に支持されている

参照基準・学習指導要領との対応

本トピックの内容・数理モデルは、文部科学省学習指導要領（算数・数学・情報）および国際標準カリキュラムに準拠し、エンジニアリング現場の実装事例と照らし合わせてMath World Lab編集部により検証・作成されています。

Layer 3 — 工学応用・現場実装

NFWプロファイルによる銀河回転曲線のフィット計算コード

宇宙論シミュレーションで標準的に用いられるNavarro-Frenk-White (NFW) 暗黒物質ハロープロファイルから銀河回転速度曲線を計算するPythonコードです。

1import numpy as np
2
3# 物理定数
4G = 4.30091e-6 # (km/s)^2 * kpc / M_sun
5
6# 天の川銀河の代表パラメータ
7# 1. 円盤バリオン質量 (M_disk = 5e10 M_sun, スケール長 R_d = 3.0 kpc)
8def v_disk(r, M_d=5.0e10, R_d=3.0):
9    y = r / (2.0 * R_d)
10    # 簡易近似プロファイル
11    return np.sqrt(G * M_d * (r**2) / ((r**2 + R_d**2)**1.5))
12
13# 2. NFWダークマターハロー (M_vir = 1e12 M_sun, 濃度 c=12, R_vir = 200 kpc)
14def v_halo_nfw(r, v200=150.0, c=12.0, R_vir=200.0):
15    rs = R_vir / c
16    x = r / rs
17    # NFW内包質量関数
18    f_c = np.log(1.0 + c) - c / (1.0 + c)
19    f_x = np.log(1.0 + x) - x / (1.0 + x)
20    m_enc = (f_x / f_c) * 1.0e12
21    return np.sqrt(G * m_enc / r)
22
23radii_kpc = np.array([2.0, 5.0, 8.5, 15.0, 25.0, 40.0]) # 太陽系は 8.5 kpc
24print(f"{'距離 (kpc)':<10} | {'星のみ (km/s)':<15} | {'DM合算 (km/s)':<15}")
25print("-" * 45)
26
27for r in radii_kpc:
28    vd = v_disk(r)
29    vh = v_halo_nfw(r)
30    v_total = np.sqrt(vd**2 + vh**2)
31    print(f"{r:<10.1f} | {vd:<15.1f} | {v_total:<15.1f}")