リープ・フロッグ法

リープ・フロッグ法は、微分方程式の数値積分法 (常微分方程式の数値解法) の一種、2次のシンプレクティック数値積分法である^[1]。リープ・フロッグ法は、

${\displaystyle {\ddot {x}}={\frac {d^{2}x}{dt^{2}}}=F(x)}$

または

${\displaystyle {\dot {v}}={\frac {dv}{dt}}=F(x),\quad {\dot {x}}={\frac {dx}{dt}}=v}$

という形式の微分方程式を解く際に用いられ、特に、古典力学における力学系の計算で重要である。

リープ・フロッグ法における時間積分は、互いの上を蛙跳び (Leap-frog) するように位置 ${\displaystyle x(t)}$ と速度 ${\displaystyle v(t)={\dot {x}}(t)}$ をずらして時間発展させるのが特徴である。リープ・フロッグ法は、オイラー法が一次精度であるのとは対照的に、二次精度の数値積分法である。また、オイラー法とは異なり、時間幅が定数 ${\displaystyle \Delta t}$ でありなおかつ ${\displaystyle \Delta t\leq 2/\omega }$ なる周期運動で安定となる^[2]。

リープ・フロッグ法では、以下の式で位置と速度を更新する。

${\displaystyle {\begin{aligned}x_{i}&=x_{i-1}+v_{i-1/2}\Delta t,\\a_{i}&=F(x_{i}),\\v_{i+1/2}&=v_{i-1/2}+a_{i}\Delta t.\end{aligned}}}$

ここで、${\displaystyle x_{i}}$ は ${\displaystyle i}$ ステップ目での位置で、${\displaystyle v_{i+1/2}}$ は ${\displaystyle i+1/2}$ ステップ目の速度、${\displaystyle a_{i}}$ は ${\displaystyle i}$ ステップ目の加速度である。${\displaystyle \Delta t}$ は時間ステップの大きさである。これらの式は、半整数ステップを消去することによって以下のような整数ステップのみの式で表現することができる^[3]。

${\displaystyle {\begin{aligned}x_{i+1}&=x_{i}+v_{i}\Delta t+{\frac {1}{2}}a_{i}\Delta t^{2},\\v_{i+1}&=v_{i}+{\frac {1}{2}}(a_{i}+a_{i+1})\Delta t.\end{aligned}}}$

ただし、この形式では、時間ステップ ${\displaystyle \Delta t}$ が一定値でない限り安定ではない^[4]。

リープ・フロッグ法は、加速度が速度に非依存である、重力計算に用いられることが多い^[5]。なお、重力計算にはルンゲ・クッタ法のような高次精度の数値積分法もよく用いられている。

力学系のシミュレーションに際して、リープ・フロッグ法にはいくつか利点がある。一つ目は、時間可逆性である。これは、${\displaystyle n}$ 段時間積分したのち、時間を逆向きに ${\displaystyle n}$ 段数値積分すると、初期位置に戻るという性質である。二つ目は、シンプレクティック性であり、これはエネルギー保存性を意味している。この性質は、軌道力学において有用である。4次ルンゲ・クッタ法のような他の多くの数値積分法は、エネルギーが保存せず、時間とともに誤差がどんどん増大してしまう。時間可逆性やシンプレクティック性から、リープ・フロッグ法はハミルトン・モンテカルロ法にも用いられている。ハミルトン・モンテカルロ法は、直接サンプリングが困難な確率分布からランダムサンプルを得るための手法である^[6]。

関連項目

• オイラー法
• ルンゲ・クッタ法
• ハミルトン力学
• 多体問題
• N体シミュレーション
• シンプレクティック数値積分法

参考文献