重力多体系

万有引力を及ぼしあう $N$ 個の天体系において，特定の天体が受ける力の記述。OKWaveより。

ポテンシャルエネルギーと力をつなぐ計算だから，決して難しい問題ではないが，エレガントに記述するには工夫がいる。

【問題】

万有引力を及ぼしあう $N$ 個の天体系において， $i$ 番目の天体が他の天体から受ける力は

$\boldsymbol{F}_i = -\sum_{j=1}^N \frac{GM_iM_j(\boldsymbol{R}_i-\boldsymbol{R}_j)}{|\boldsymbol{R}_i-\boldsymbol{R}_j|^3}$ 　…(1)

となる。ただし， $M_i,\boldsymbol{R}_i=(x_i,y_i,z_i)$ は $i$ 番目の天体の質量および位置である。
ここで関数

$V(\boldsymbol{R}_1,\boldsymbol{R}_2,\cdots,\boldsymbol{R}_N) = -\sum_{i<j}\frac{GM_iM_j}{|\boldsymbol{R}_i-\boldsymbol{R}_j|}$

を定義すると，(1)は

$\boldsymbol{F}_i = -\frac{\partial}{\partial \boldsymbol{R}_i}V(\boldsymbol{R}_1,\boldsymbol{R}_2,\cdots,\boldsymbol{R}_N) = -\left(\frac{\partial V}{\partial x_i},\frac{\partial V}{\partial y_i},\frac{\partial V}{\partial z_i}\right)$

と書けることを示せ。

【解答】

$|\boldsymbol{R}_i-\boldsymbol{R}_j| = \sqrt{(x_i-x_j)^2+(y_i-y_j)^2+(z_i-z_j)^2}$

より，

$\frac{\partial}{\partial x_i}\frac{1}{|\boldsymbol{R}_i-\boldsymbol{R}_j|} = -\frac{x_i-x_j}{|\boldsymbol{R}_i-\boldsymbol{R}_j|^3}\qquad,\qquad\frac{\partial}{\partial x_j}\frac{1}{|\boldsymbol{R}_i-\boldsymbol{R}_j|} = \frac{x_i-x_j}{|\boldsymbol{R}_i-\boldsymbol{R}_j|^3}$

面倒なので，

$f(i,j) = -\frac{GM_iM_j(x_i - x_j)}{|\boldsymbol{R}_i-\boldsymbol{R}_j|^3} = -f(j,i)$

と書くと，力 $\boldsymbol{F}_i$ の $x$ 成分は

$F_{ix} = \sum_jf(i,j)$

となるが，これを示せばよい。

$F_{ix} = -\frac{\partial V}{\partial x_i} = \frac{\partial}{\partial x_i}\sum_{i<j}\frac{GM_iM_j}{|\boldsymbol{R}_i-\boldsymbol{R}_j|}$
　　　 $= \sum_{k<j} \left[ \frac{\partial x_k}{\partial x_i}f(k,j) - \frac{\partial x_j}{\partial x_i}f(k,j) \right]$
　　　 $= \sum_{j(>i)}f(i,j) - \sum_{k(<i)}f(k,i)$ 　　※　 $i$ は固定されていることに注意
　　　 $= \sum_{j(>i)}f(i,j) + \sum_{k(<i)}f(i,k)$
　　　 $= \sum_jf(i,j)$ 　　※　 $f(i,i)=0$

他の成分についても同様。（証明終）

※　 $\frac{\partial x_k}{\partial x_i} = \delta_{ki} = \Big\{\begin{matrix}1\quad(k=i)\\0\quad(k\ne i)\end{matrix}$

タグ：

+ タグ編集

「重力多体系」をウィキ内検索

最終更新：2010年11月13日 22:38

科学のおもちゃ箱 @wiki

重力多体系

重力多体系

メニュー

更新履歴