【日本語で解説】Jeans方程式の導出とその応用-Shinoryo's Weblog

Chuk YongによるPixabay(https://pixabay.com/)からの画像

無衝突Boltzmann方程式(collisionless Boltzmann equation)とは, あるポテンシャル下で運動する, 恒星系(必ずしもそれに限らないが)の位相空間 $(r, v)$ 上の分布関数 $f (r, v)$ の時間発展を与える式である. Jeans方程式(Jeans equation)とは, 無衝突Boltzmann方程式を速度空間で平均化した方程式である. ここでは, 無衝突Boltzmann方程式について紹介し, 無衝突Boltzmann方程式からJeans方程式を導く.

本記事のPDFver.はこちら $\Rightarrow$ https://drive.google.com/file/d/1j5hdJ3538rjs70xiZWCU6WnzYJGOcp51/view?usp=sharing

目次[非表示]

1恒星系の分布関数
2無衝突Boltzmann方程式
3Jeans方程式
4球対称Jeans方程式
- 4.1球対称Jeans方程式(微分形)
- 4.2球対称Jeans方程式(積分形)
5ビリアル定理
6参考文献
7脚注

恒星系の分布関数

「恒星系の分布関数」 $f$ とは, 一般には以下のように定義される.

$f (r, v, t) d^{3} r d^{3} v$ ：ある時刻 $t$ において, 位相空間のある位置 $(r, v)$ における微小体積要素 $d^{3} r d^{3} v$ における存在確率

これを用いると, ある時刻 $t$ において, 空間位置 $r$ における数密度 $n (r, t)$ は

\begin{aligned} (1) & n (r, t) & = \int f (r, v, t) d^{3} v \end{aligned}

と書ける. また, 各方向の平均速度 $⟨ v_{i} ⟩$ は,

\begin{aligned} (2) & ⟨ v_{i} ⟩ & = \frac{1}{n} \int v_{i} f d^{3} v \end{aligned}

と書ける.

無衝突Boltzmann方程式

無衝突系(collisionless)：位相空間の各位置において, 恒星が新たに生まれたりあるいは死んだりして数が変化せず, また恒星同士が衝突しない系

無衝突系であれば, 恒星系の分布関数は位相空間における連続の式に従う：

\begin{aligned} (3) & \frac{\partial f}{\partial t} + \sum_{i = 1}^{3} \frac{\partial (f {\dot{x}}_{i})}{\partial x_{i}} + \sum_{i = 1}^{3} \frac{\partial (f {\dot{v}}_{i})}{\partial v_{i}} & = 0 . \end{aligned}

ここで,

${\dot{x}}_{i} = v_{i}$ と $x_{i}$ は互いに独立
運動方程式より ${\dot{v}}_{i}$ は $x_{i}$ だけの関数

であるから, 式 $(3)$ より

\begin{aligned} (4) & \frac{\partial f}{\partial t} + v \cdot \nabla f + \dot{v} \cdot \frac{\partial f}{\partial v} & = 0 \end{aligned}

となり, 運動方程式

\begin{aligned} (5) & \frac{d v}{d t} & = - \nabla Φ \end{aligned}

より

\begin{array}{r} (6) & \frac{\partial f}{\partial t} + v \cdot \nabla f - \nabla Φ \cdot \frac{\partial f}{\partial v} = 0 \end{array}

となる(無衝突Boltzmann方程式(collisionless Boltzmann equation)).

Jeans方程式

速度空間上での積分

式 $(6)$ の両辺を速度空間上で積分する. 各項計算すると, 第1項は

\begin{aligned} (7) & \int \frac{\partial f}{\partial t} d^{3} v & = \frac{\partial}{\partial t} \int f d^{3} v = \frac{\partial n}{\partial t} \end{aligned}

となり, 第2項は

\begin{aligned} \int v_{i} \frac{\partial f}{\partial x_{i}} d^{3} v & = \int \frac{\partial (f v_{i})}{\partial x_{i}} d^{3} v \\ = \frac{\partial}{\partial x_{i}} \int f v_{i} d^{3} v \\ (8) & = \frac{\partial (n ⟨ v_{i} ⟩)}{\partial x_{i}} \end{aligned}

となる. 最後の項は,

$Φ$ は $v$ には依存しないこと
$v \to \infty$ で十分速く $f \to 0$ になること

を利用すれば, Gaussの発散定理より表面項のみが残り,

\begin{array}{r} (9) & \int \frac{\partial Φ}{\partial x_{i}} \frac{\partial f}{\partial v_{i}} d^{3} v = \int \frac{\partial}{\partial v_{i}} (f \frac{\partial Φ}{\partial x_{i}}) d^{3} v = 0 \end{array}

となる. 以上から,

\begin{aligned} (10) & \frac{\partial n}{\partial t} + \sum_{i = 1}^{3} \frac{\partial (n ⟨ v_{i} ⟩)}{\partial x_{i}} & = 0 \end{aligned}

となる.

$v_{i}$ のモーメントをかけて積分

式 $(6)$ に $v_{i}$ のモーメントをかけて積分する. 各項計算すると, 第1項は

\begin{aligned} (11) & \int v_{j} \frac{\partial f}{\partial t} d^{3} v & = \int \frac{\partial (f v_{j})}{\partial t} d^{3} v = \frac{\partial}{\partial t} \int f v_{j} d^{3} v = \frac{\partial (n ⟨ v_{j} ⟩)}{\partial t} \end{aligned}

となり, 第2項は

\begin{aligned} \int v_{i} v_{j} \frac{\partial f}{\partial x_{i}} d^{3} v & = \int \frac{\partial (f v_{i} v_{j})}{\partial x_{i}} d^{3} v \\ = \frac{\partial}{\partial x_{i}} \int f v_{i} v_{j} d^{3} v \\ (12) & = \frac{\partial (n ⟨ v_{i} v_{j} ⟩)}{\partial x_{i}} \end{aligned}

となる. 最後の項は,

$Φ$ は $v$ には依存しないこと
$v \to \infty$ で十分速く $f \to 0$ になること
$\partial v_{j} / \partial v_{i} = δ_{i j}$

を利用すれば,

\begin{aligned} \int v_{j} \frac{\partial Φ}{\partial x_{i}} \frac{\partial f}{\partial v_{i}} d^{3} v & = \frac{\partial Φ}{\partial x_{i}} \int v_{j} \frac{\partial f}{\partial v_{i}} d^{3} v \\ = \frac{\partial Φ}{\partial x_{i}} \int (\frac{\partial (f v_{j})}{\partial v_{i}} - f \frac{\partial v_{j}}{\partial v_{i}}) d^{3} v \\ = - \frac{\partial Φ}{\partial x_{i}} n δ_{i j} \\ (13) & = - n \frac{\partial Φ}{\partial x_{j}} \end{aligned}

となる. 以上から,

\begin{aligned} (14) & \frac{\partial (n ⟨ v_{j} ⟩)}{\partial t} + \sum_{i = 1}^{3} \frac{\partial (n ⟨ v_{i} v_{j} ⟩)}{\partial x_{i}} + n \frac{\partial Φ}{\partial x_{j}} & = 0, \end{aligned}

つまり

\begin{aligned} (15) & \frac{\partial n}{\partial t} ⟨ v_{j} ⟩ + n \frac{\partial ⟨ v_{j} ⟩}{\partial t} + \sum_{i = 1}^{3} \frac{\partial (n ⟨ v_{i} v_{j} ⟩)}{\partial x_{i}} + n \frac{\partial Φ}{\partial x_{j}} & = 0, \end{aligned}

Jeans方程式

式 $(10)$ , $(15)$ から $\partial n / \partial t$ を消去すれば,

\begin{aligned} (16) & - ⟨ v_{j} ⟩ \sum_{i = 1}^{3} \frac{\partial (n ⟨ v_{i} ⟩)}{\partial x_{i}} + n \frac{\partial ⟨ v_{j} ⟩}{\partial t} + \sum_{i = 1}^{3} \frac{\partial (n ⟨ v_{i} v_{j} ⟩)}{\partial x_{i}} + n \frac{\partial Φ}{\partial x_{j}} & = 0 \end{aligned}

となる. 速度分散 $σ_{i j}$ を

\begin{aligned} (17) & σ_{i j}^{2} & = ⟨ (v_{i} - ⟨ v_{i} ⟩) (v_{j} - ⟨ v_{j} ⟩) ⟩ = ⟨ v_{i} v_{j} ⟩ - ⟨ v_{i} ⟩ ⟨ v_{j} ⟩ \end{aligned}

のように導入すれば, 整理することで

\begin{array}{r} (18) & n \frac{\partial ⟨ v_{j} ⟩}{\partial t} + \sum_{i = 1}^{3} n ⟨ v_{i} ⟩ \frac{\partial ⟨ v_{j} ⟩}{\partial x_{i}} = - n \frac{\partial Φ}{\partial x_{j}} - \sum_{i = 1}^{3} \frac{\partial (n σ_{i j}^{2})}{\partial x_{i}} \end{array}

を得る(Jeans方程式(Jeans equation)).

球対称Jeans方程式

球対称Jeans方程式(微分形)

ここでは, 球対称なJeans方程式を考えていく.

式 $(6)$ を球座標で記述すると,

\begin{aligned} (19) & \frac{\partial f}{\partial t} + v_{r} \frac{\partial f}{\partial r} + \frac{v_{θ}}{r} \frac{\partial f}{\partial θ} + \frac{v_{ϕ}}{r \sin θ} \frac{\partial f}{\partial ϕ} + {\dot{v}}_{r} \frac{\partial f}{\partial v_{r}} + {\dot{v}}_{θ} \frac{\partial f}{\partial v_{θ}} + {\dot{v}}_{ϕ} \frac{\partial f}{\partial v_{ϕ}} & = 0 \end{aligned}

となる. ここで, 座標微分と速度成分の関係は

\begin{array}{r} (20) & {\begin{aligned} v_{r} & = \dot{r}, \\ v_{θ} & = r \dot{θ}, \\ v_{ϕ} & = r \sin θ \cdot \dot{ϕ} \end{aligned} \end{array}

であり, また速度成分の微分は

\begin{array}{r} (21) & {\begin{aligned} {\dot{v}}_{r} & = \frac{v_{θ}^{2} + v_{ϕ}^{2}}{r} - \frac{d Φ}{d r}, \\ {\dot{v}}_{θ} & = \frac{v_{ϕ}^{2} \cot θ - v_{r} v_{θ}}{r}, \\ {\dot{v}}_{ϕ} & = - \frac{v_{ϕ} v_{r} + v_{ϕ} v_{θ} \cot θ}{r} \end{aligned} \end{array}

となることに注意せよ( $Φ$ は $r$ にしか依存しないと仮定している).

式 $(19)$ をそのまま速度空間上で積分する. 式 $(19)$ をそのまま速度空間上で積分した結果は,

$v \to \infty$ で十分速く $f \to 0$ になること

を利用すれば,

\begin{aligned} (22) & \frac{\partial n}{\partial t} + \frac{\partial (n ⟨ v_{r} ⟩)}{\partial r} + \frac{1}{r} \frac{\partial (n ⟨ v_{θ} ⟩)}{\partial θ} + \frac{1}{r \sin θ} \frac{\partial (n ⟨ v_{ϕ} ⟩)}{\partial ϕ} + n \frac{2 ⟨ v_{r} ⟩}{r} + n \frac{⟨ v_{θ} ⟩ \cot θ}{r} & = 0 \end{aligned}

となる.

式 $(19)$ に $v_{r}$ をかけて, 速度空間上で積分する. 式 $(22)$ を考慮すれば,

\begin{aligned} n \frac{\partial ⟨ v_{r} ⟩}{\partial t} + n (⟨ v_{r} ⟩ \frac{\partial ⟨ v_{r} ⟩}{\partial r} + \frac{⟨ v_{θ} ⟩}{r} \frac{\partial ⟨ v_{r} ⟩}{\partial θ} + \frac{⟨ v_{ϕ} ⟩}{r \sin θ} \frac{\partial ⟨ v_{r} ⟩}{\partial ϕ}) \\ + \frac{\partial (n σ_{r r}^{2})}{\partial r} + \frac{1}{r} \frac{\partial (n σ_{r θ}^{2})}{\partial θ} + \frac{1}{r \sin θ} \frac{\partial (n σ_{r ϕ}^{2})}{\partial ϕ} \\ + \frac{n}{r} {2 σ_{r r}^{2} - (σ_{θ θ}^{2} + σ_{ϕ ϕ}^{2} + {⟨ v_{θ} ⟩}^{2} + {⟨ v_{ϕ} ⟩}^{2}) + σ_{r θ}^{2} \cot θ} \\ (23) & = - n \frac{d Φ}{d r} \end{aligned}

となる.

式 $(23)$ に, さまざまな仮定を加えていく.

分布関数の球面対称性を仮定する.
定常状態を考え, $\partial / \partial t = 0$
平均からの差を考え, $⟨ v_{r} ⟩ = ⟨ v_{θ} ⟩ = ⟨ v_{ϕ} ⟩ = 0$
視線速度と接線速度には特に相関関係はなく, $σ_{r θ}^{2} = σ_{r ϕ}^{2} = 0$
$σ_{r r}^{2} = σ_{r}^{2}$ , $σ_{θ θ}^{2} = σ_{θ}^{2}$ , $σ_{ϕ ϕ}^{2} = σ_{ϕ}^{2}$ と書く

これらを用いれば,

\begin{aligned} (24) & \frac{d (n σ_{r}^{2})}{d r} + \frac{n}{r} (2 σ_{r}^{2} - σ_{θ}^{2} - σ_{ϕ}^{2}) & = - n \frac{d Φ}{d r} \\ (25) & ∴ \frac{1}{n} \frac{d (n σ_{r}^{2})}{d r} + \frac{2 σ_{r}^{2} - σ_{θ}^{2} - σ_{ϕ}^{2}}{r} & = - \frac{d Φ}{d r} \end{aligned}

となる. 異方性パラメータ $β$ を

\begin{aligned} (26) & β & = 1 - \frac{σ_{θ}^{2} + σ_{ϕ}^{2}}{2 σ_{r}^{2}} \end{aligned}

とすれば,

\begin{array}{r} (27) & \frac{1}{n (r)} \frac{d (n (r) σ_{r}^{2} (r))}{d r} + \frac{2 β (r) σ_{r}^{2} (r)}{r} = - \frac{d Φ (r)}{d r} \end{array}

を得る. 式 $(27)$ を, 球面対称性を仮定したJeans方程式という.

球対称Jeans方程式(積分形)

式 $(27)$ は $n (r) σ_{r}^{2} (r)$ に関する非同次1階線形微分方程式である. まず同次方程式

\begin{array}{r} (28) & \frac{d (n (r) σ_{r}^{2} (r))}{d r} + \frac{2 β_{ani} (r)}{r} n (r) σ_{r}^{2} (r) = 0 \end{array}

を考えると, その解は

\begin{aligned} (29) & n (r) σ_{r}^{2} (r) & = C \exp [- \int_{\infty}^{r} d r^{'} \frac{2 β_{ani} (r^{'})}{r^{'}}] = C \exp [\int_{r}^{\infty} d r^{'} \frac{2 β_{ani} (r^{'})}{r^{'}}] \end{aligned}

となる( $C$ ：任意定数). 次に, 非同次方程式

\begin{array}{r} (30) & \frac{d (n (r) σ_{r}^{2} (r))}{d r} + \frac{2 β_{ani} (r)}{r} n (r) σ_{r}^{2} (r) = - n (r) \frac{d Φ (r)}{d r} \end{array}

を考えるために, 式 $(29)$ の任意定数 $C$ を関数 $C (r)$ に置き換える(定数変化法). 式 $(29)$ の任意定数 $C$ を関数 $C (r)$ に置き換えたものを式 $(30)$ に代入すると,

\begin{aligned} (31) & \frac{d C (r)}{d r} & = - n (r) \frac{d Φ (r)}{d r} \exp [- \int_{r}^{\infty} d r^{'} \frac{2 β_{ani} (r^{'})}{r^{'}}] \end{aligned}

となり, これを積分すると

\begin{aligned} C (r) & = - \int_{\infty}^{r} d r^{''} n (r^{''}) \frac{d Φ (r^{''})}{d r^{''}} \exp [- \int_{r^{''}}^{\infty} d r^{'} \frac{2 β_{ani} (r^{'})}{r^{'}}] \\ (32) & = \int_{r}^{\infty} d r^{''} n (r^{''}) \frac{d Φ (r^{''})}{d r^{''}} \exp [- \int_{r^{''}}^{\infty} d r^{'} \frac{2 β_{ani} (r^{'})}{r^{'}}] \end{aligned}

となる. ゆえに, 非同次方程式 $(30)$ の解は

\begin{aligned} n (r) σ_{r}^{2} (r) & = \int_{r}^{\infty} d r^{''} n (r^{''}) \frac{d Φ (r^{''})}{d r^{''}} \exp [- \int_{r^{''}}^{\infty} d r^{'} \frac{2 β_{ani} (r^{'})}{r^{'}}] \exp [- \int_{\infty}^{r} d r^{'} \frac{2 β_{ani} (r^{'})}{r^{'}}] \\ = \int_{r}^{\infty} d r^{''} n (r^{''}) \frac{d Φ (r^{''})}{d r^{''}} \exp [- \int_{r^{''}}^{r} d r^{'} \frac{2 β_{ani} (r^{'})}{r^{'}}] \\ = \int_{r}^{\infty} d r^{'} n (r^{'}) \frac{d Φ (r^{'})}{d r^{'}} \exp [- \int_{r^{'}}^{r} d r^{''} \frac{2 β_{ani} (r^{''})}{r^{''}}] \\ (33) & = \int_{r}^{\infty} d r^{'} n (r^{'}) \frac{d Φ (r^{'})}{d r^{'}} \exp [\int_{r}^{r^{'}} d r^{''} \frac{2 β_{ani} (r^{''})}{r^{''}}] \end{aligned}

となる. したがって, 結果としてJeans方程式の積分形

\begin{array}{r} (34) & σ_{r}^{2} (r) = \frac{1}{n (r)} \int_{r}^{\infty} d r^{'} n (r^{'}) \frac{d Φ (r^{'})}{d r^{'}} \exp [\int_{r}^{r^{'}} d r^{''} \frac{2 β_{ani} (r^{''})}{r^{''}}] \end{array}

が得られる.

このような表式であれば, 速度分散 $σ_{r}^{2}$ の分布から右辺の重力場 $Φ (r)$ を求めることができることが, より明示的になる.

ビリアル定理

Jeans方程式をさらに座標空間でモーメント積分をとることによって, 位置によらないその恒星系の大局的な力学関係を調べることができる.

テンソルビリアル定理

式 $(14)$ の両辺に $m x_{k}$ をかけて( $m$ ：恒星の典型的な質量, これをかけて $n (r)$ を質量密度 $ρ (r)$ にする)空間全体で積分すると,

\begin{aligned} (35) & \int x_{k} \frac{\partial ρ ⟨ v_{j} ⟩}{\partial t} d^{3} x + \int x_{k} \sum_{i = 1}^{3} \frac{\partial}{\partial x_{i}} (ρ ⟨ v_{i} v_{j} ⟩) d^{3} x + \int ρ x_{k} \frac{\partial Φ}{\partial x_{j}} d^{3} x & = 0 \end{aligned}

が得られる. 以下では, $x \to \infty$ で $ρ \to 0$ の場合を考える. 式 $(35)$ の第2項は,

\begin{aligned} \int x_{k} \sum_{i = 1}^{3} \frac{\partial}{\partial x_{i}} (ρ ⟨ v_{i} v_{j} ⟩) d^{3} x & = \int (\sum_{i = 1}^{3} \frac{\partial}{\partial x_{i}} (ρ ⟨ v_{i} v_{j} ⟩ x_{k}) - \sum_{i = 1}^{3} ρ ⟨ v_{i} v_{j} ⟩ \frac{\partial x_{k}}{\partial x_{i}}) d^{3} x \\ (36) & = - \int ρ ⟨ v_{j} v_{k} ⟩ d^{3} x \end{aligned}

となる. ゆえに,

$K_{j k}$ ：運動エネルギーテンソル

\begin{array}{r} (37) & K_{j k} = \frac{1}{2} \int ρ ⟨ v_{j} v_{k} ⟩ d^{3} x \end{array}

$W_{j k}$ ：ポテンシャルエネルギーテンソル

\begin{array}{r} (38) & W_{j k} = - \int ρ x_{k} \frac{\partial Φ}{\partial x_{j}} d^{3} x \end{array}

を用いると, 式 $(35)$ は

\begin{aligned} (39) & \int x_{k} \frac{\partial ρ ⟨ v_{j} ⟩}{\partial t} d^{3} x & = 2 K_{j k} + W_{j k} \end{aligned}

と書ける. 式 $(39)$ の右辺は $i$ と $j$ の入れ替えに関して対称である^*1から, $i$ と $j$ を入れ替えたものとの和をとって2で割ると,

\begin{aligned} \frac{1}{2} (\int x_{k} \frac{\partial ρ ⟨ v_{j} ⟩}{\partial t} d^{3} x + \int x_{j} \frac{\partial ρ ⟨ v_{k} ⟩}{\partial t} d^{3} x) & = 2 K_{j k} + W_{j k}, \\ (40) & \frac{1}{2} \frac{d}{d t} (\int ρ (⟨ v_{j} ⟩ x_{k} + x_{j} ⟨ v_{k} ⟩) d^{3} x) & = 2 K_{j k} + W_{j k} \end{aligned}

となる. 一方,

$I_{j k}$ ：慣性モーメントテンソル

\begin{array}{r} (41) & I_{j k} = \int ρ x_{j} x_{k} d^{3} x \end{array}

を定義すれば, 式 $(10)$ と発散定理を用いると,

\begin{aligned} \frac{d I_{j k}}{d t} & = \frac{d}{d t} \int ρ x_{j} x_{k} d^{3} x \\ = \int \frac{\partial ρ}{\partial t} x_{j} x_{k} d^{3} x \\ = - \sum_{i = 1}^{3} \int \frac{\partial ρ ⟨ v_{i} ⟩}{\partial x_{i}} x_{j} x_{k} d^{3} x \\ = - \sum_{i = 1}^{3} \int \frac{\partial (ρ ⟨ v_{i} ⟩ x_{j} x_{k})}{\partial x_{i}} d^{3} x + \sum_{i = 1}^{3} \int ρ ⟨ v_{i} ⟩ \frac{\partial (x_{j} x_{k})}{\partial x_{i}} d^{3} x \\ = \sum_{i = 1}^{3} \int ρ ⟨ v_{i} ⟩ \frac{\partial (x_{j} x_{k})}{\partial x_{i}} \\ (42) & = \int ρ (⟨ v_{j} ⟩ x_{k} + x_{j} ⟨ v_{k} ⟩) d^{3} x \end{aligned}

であるから,

\begin{aligned} (43) & \frac{1}{2} \frac{d}{d t} \frac{d I_{j k}}{d t} & = 2 K_{j k} + W_{j k} \end{aligned}

が成り立つ. ゆえに, 最終的には

\begin{array}{r} (44) & \frac{1}{2} \frac{d^{2} I_{j k}}{d t^{2}} = 2 K_{j k} + W_{j k} \end{array}

というテンソルビリアル定理(Tensor Virial Theorem)を表す式に帰着できる.

スカラービリアル定理

系が力学平衡にあるときは, 式 $(44)$ の左辺は $0$ になり,

\begin{aligned} 2 K_{j k} + W_{j k} & = 0 \end{aligned}

となる. この両辺のトレースをとって,

$K$ ：系全体の運動エネルギー

\begin{aligned} (45) & K & = \sum_{j} \frac{1}{2} \int ρ ⟨ v_{j} v_{j} ⟩ d^{3} x \end{aligned}

$W$ ：系全体のポテンシャルエネルギー

\begin{aligned} (46) & W & = - \sum_{j} \int ρ x_{j} \frac{\partial Φ}{\partial x_{j}} d^{3} x \end{aligned}

を定義すれば,

\begin{array}{r} (47) & 2 K + W = 0 \end{array}

が得られる(スカラービリアル定理(Scalar Virial Theorem)). 系の全エネルギーを

\begin{aligned} (48) & E & = K + W \end{aligned}

とすると,

\begin{array}{r} (49) & E = - K = \frac{W}{2} \end{array}

が得られる. すなわち,

系からエネルギーを奪うと, 系の運動エネルギーが増える(系の「比熱」が負)
系のポテンシャルエネルギーは, 系の全エネルギーの半分である

が成り立つ.

スカラービリアル半径の応用

系全体の運動エネルギー $K$ は, 系全体の質量

\begin{aligned} (50) & M & = \int ρ d^{3} x \end{aligned}

と系全体の平均的な速度分散

\begin{aligned} (51) & ⟨ v^{2} ⟩ & = \frac{1}{M} \sum_{j} \int ρ ⟨ v_{j}^{2} ⟩ d^{3} x \end{aligned}

を定義すれば,

\begin{aligned} (52) & K & = \frac{1}{2} M ⟨ v^{2} ⟩ \end{aligned}

と書けることに注意せよ. また, 系全体のポテンシャルエネルギー $W$ は, ビリアル半径と呼ばれる半径 $r_{v}$ を

\begin{aligned} (53) & \frac{1}{r_{v}} & = - \frac{1}{M^{2}} \int \frac{ρ (x) ρ (x^{'})}{| x - x^{'} |} d^{3} x d x^{'} \end{aligned}

のように定義すれば,

\begin{aligned} (54) & W & = - \frac{G M^{2}}{2 r_{v}} = - \frac{G M^{2}}{r_{g}} \end{aligned}

と書けることに注意せよ( $r_{g}$ は重力半径(gravitational radius)と呼ばれる).

式 $(47)$ , $(52)$ , $(54)$ より,

\begin{array}{r} (55) & ⟨ v^{2} ⟩ = \frac{G M}{r_{g}} \end{array}

と書ける. ゆえに, 系全体の平均的な速度分散 $⟨ v^{2} ⟩$ と重力半径 $r_{g}$ が観測からわかれば, 系全体の質量 $M$ を求めることができる.

重力半径 $r_{g}$ には曖昧さが残るが, 典型的には「質量の半分が入っている半径」(射影はされていない) $r_{h}$ との関係が $r_{h} ≃ r_{g}$ であるという経験則が知られており,

\begin{aligned} (56) & ⟨ v^{2} ⟩ & = 0.4 \frac{G M}{r_{h}} \end{aligned}

となるそうである. 例えば, 密度一様で半径が $2^{1 / 3} r_{h}$ となるような^*2仮想球を考えると,

\begin{aligned} (57) & ⟨ v^{2} ⟩ & = \frac{3 G M}{5 \times 2^{1 / 3} r_{h}} = 0.48 \frac{G M}{r_{h}} \end{aligned}

となる^*3. また, 球状星団の分布をよく説明するPlummerモデル

\begin{aligned} (58) & ρ (r) & = \frac{3}{4 π} \frac{M a^{2}}{{(r^{2} + a^{2})}^{5 / 2}} \end{aligned}

を仮定する( $a > 0$ はPlummer半径と呼ばれるスケール半径)と,

\begin{aligned} (59) & ⟨ v^{2} ⟩ & = \frac{3 π G M}{32 \sqrt{2^{2 / 3} - 1} r_{h}} ≃ 0.38 \frac{G M}{r_{h}} \end{aligned}

となる^*4.

参考文献

Binney, James; Tremaine, Scott. Galactic dynamics. 2nd ed., Princeton, Princeton University Press, 2008, 885p., (Princeton series in astrophysics), ISBN 978-0-691-13026-2.

Princeton University Press

Galactic Dynamics (Princeton Series in Astrophysics)

23369329

Amazonで見る楽天市場で見る Yahoo!ショッピングで見る

Schechter, Paul. 8.902 Astrophysics II. Fall 2004. Massachusetts Institute of Technology: MIT OpenCourseWare, https://ocw.mit.edu. License: Creative Commons BY-NC-SA.

千葉柾司. 銀河考古学. 日本評論社, 2015, p.44-50, (新天文学ライブラリー, 2), ISBN 978-4-535-60741-5.

銀河考古学新天文学ライブラリー

Amazonで見る楽天市場で見る Yahoo!ショッピングで見る

牧野淳一郎. “理論天体物理学特論I 講義資料(3回目)”. 2003-11-09. http://jun-makino.sakura.ne.jp/kougi/stellar_dynamics_2004/note3/note3-e.html, (参照 2022-05-03).

脚注

*1： $W_{j k}$ が対称テンソルであるというのは, 自己重力系であるということから分かる. 自己重力系であるという仮定から, 重力ポテンシャルは,

\begin{aligned} (60) & Φ (x) & = - \int \frac{G ρ (x^{'})}{| x - x^{'} |} d^{3} x^{'}, \end{aligned}

となり, その勾配は

\begin{aligned} (61) & \nabla_{x} Φ (x) & = \int \frac{G ρ (x^{'}) (x - x^{'})}{{| x - x^{'} |}^{3}} d^{3} x^{'} \end{aligned}

となる. これを用いて, $W_{j k}$ は

\begin{aligned} (62) & W_{j k} & = - \iint \frac{G ρ (x) ρ (x^{'}) x_{k} (x_{j} - x_{j}^{'})}{{| x - x^{'} |}^{3}} d^{3} x d^{3} x^{'} \end{aligned}

となる. 積分変数 $x$ と $x^{'}$ を入れ替えたものを加えて $2$ で割ると,

\begin{aligned} W_{j k} & = - \frac{1}{2} (\iint \frac{G ρ (x) ρ (x^{'}) x_{k} (x_{j} - x_{j}^{'})}{{| x - x^{'} |}^{3}} d^{3} x d^{3} x^{'} \\ + \iint \frac{G ρ (x^{'}) ρ (x) x_{k}^{'} (x_{j}^{'} - x_{j})}{{| x^{'} - x |}^{3}} d^{3} x d^{3} x^{'}) \\ (63) & = - \frac{1}{2} \iint \frac{G ρ (x) ρ (x^{'}) (x_{j} - x_{j}^{'}) (x_{k} - x_{k}^{'})}{{| x - x^{'} |}^{3}} d^{3} x d^{3} x^{'} \end{aligned}

となり, $W_{j k}$ が対称テンソルであることが分かる. 以降の議論においては, この系が自己重力系であることを仮定していることに注意せよ.

*2： $r_{h}$ が「質量の半分が入っている半径」となるように定めた.

*3：密度一様な球内での重力ポテンシャルは

\begin{aligned} (64) & Φ (r) & = \frac{2 π G ρ}{3} r^{2} - 2^{4 / 3} π G ρ r_{h} \end{aligned}

となるから, $W$ を計算すると

\begin{aligned} W & = - \int_{0}^{2^{1 / 3} r_{h}} ρ r \frac{d Φ (r)}{d r} 4 π r^{2} d r \\ = - \frac{16 π^{2} G ρ^{2}}{3} \int_{0}^{2^{1 / 3} r_{h}} r^{4} d r \\ = - \frac{16 π^{2} G ρ^{2} 2^{5 / 3} r_{h}^{5}}{15} \\ (65) & = - \frac{3 G M^{2}}{5 \times 2^{1 / 3} r_{h}} \end{aligned}

となる.

*4：Plummerモデルに従う恒星系が作り出す重力ポテンシャルは,

\begin{aligned} (66) & Φ (r) & = - \frac{G M}{\sqrt{r^{2} + a^{2}}} \end{aligned}

となるから, $W$ を計算すると

\begin{aligned} W & = - \int_{0}^{a} ρ (r) r \frac{d Φ (r)}{d r} 4 π r^{2} d r \\ (67) & = - \frac{3 π G M^{2}}{32 a} \end{aligned}

となる. また, 半径 $r$ 以内の質量 $M (r)$ は

\begin{aligned} (68) & M (r) & = \frac{M r^{3}}{{(r^{2} + a^{2})}^{3 / 2}} \end{aligned}

であるから, 「質量の半分が入っている半径」 $r_{h}$ と $a$ の関係は

\begin{aligned} (69) & a & = \sqrt{2^{2 / 3} - 1} r \end{aligned}

となる. ゆえに

\begin{aligned} (70) & W & = - \frac{3 π G M^{2}}{32 \sqrt{2^{2 / 3} - 1} r} \end{aligned}

のように計算できる.

【日本語で解説】Jeans方程式の導出とその応用

恒星系の分布関数

無衝突Boltzmann方程式