三体問題シミュレーター

カオス的な重力、双子軌道の発散、数値的安定性

インタラクティブシミュレーションを読み込んでいます...

なぜ更新するたびに同じ結果になるのか

三体問題はカオス的ですが、ランダムではありません。カオスとは、初期条件のわずかな違いが指数関数的に増大することを意味しますが、ページを更新しても何も変わりません。同じ初期位置、同じ初速度、同じdt、同じアルゴリズムなら、毎回同じ結果になります。決定論的カオスと計算の再現性は、まったく別の概念です。カオスを実際に見るには、摂動を0.000001に設定し、破線で示される双子軌道を観察してください。両方の軌道は10⁻⁶だけ離れた位置から始まり、やがてまったく異なる経路へと発散していきます。その指数関数的な分離こそがカオスであり、リアプノフ指数が働いている証拠です。軌道の下にある発散グラフは、その分離量のlog₁₀を時間の経過とともにプロットしたもので、直線的に上昇していれば指数関数的な成長を意味します。

なぜ三体が難しいのか

互いに引き合う二つの天体は、単純で永遠に繰り返す楕円を描きます。この場合はニュートンが紙と鉛筆で解きました。ところが三つ目の天体を加えると、きれいな公式は消え去ります。誰がどこへ行き着くかを表す一般的な方程式は存在せず、唯一の方法は、このツールのように微小な dt ずつシミュレーションを進めることだけです。figure-8 プリセットを押すと、まれな例外の一つ、三つの等しい質量が単一のループ状の軌道を永遠に追いかけ合う様子が見られます。

存在するのに役に立たない解

三体問題は「解けない」と言われがちですが、1912年にカール・スンドマンは t1/3 のべき級数が収束し、あらゆる時刻の天体の位置を与えることを証明しました(三体衝突の初期条件を除く)。ただし収束が極端に遅く、天文学的な精度を得るには 108000000 項ものオーダーが必要で、どんなコンピュータでも到底不可能です。形式的な解は確かに存在しますが、ここでのように一歩ずつシミュレーションを進めることが、唯一実用的な手段なのです。

全プロセスの詳細解説

  1. 初期条件から求める3つの等質量の8の字軌道のエネルギー 7 ステップ

    8の字コレオグラフィーは3つの等質量体を単一の閉曲線上に周回させます。初期条件のみから t = 0 でのエネルギーを求めなさい。 これは8の字軌道のプリセットです:G = 1、m1 = m2 = m3 = 1、ソフトニング ε = 0。位置 r1 = (−0.97000436, 0.24308753)、r2 = −r1、r3 = (0, 0);速度 v1 = v2 = (0.466203685, 0.43236573)、そして全運動量は0です。

    この問題をツールで開く

    1. 一見すると恣意的に見える長い小数ですが、自乗して足し合わせるとそうではないことが分かります。質点1と質点2はそれぞれ、質点3の在処である原点から正確に1単位の距離に位置しています。

      ∣r1∣=0.970004362+0.243087532=0.940908+0.059092=1.0000|\mathbf{r}_1|=\sqrt{0.97000436^{2}+0.24308753^{2}}=\sqrt{0.940908+0.059092}=1.0000
    2. 質点2は質点1を原点に関して点対称に移動した位置にあるため、この瞬間、3つの質点すべてが質点3を中点とする1本の直線上に並んでいます。すべての質点間距離はそこから導かれ、2つ目の平方根の計算は不要です。

      r2=−r1⟹r13=r23=1.0000,r12=2.0000\mathbf{r}_2=-\mathbf{r}_1\quad\Longrightarrow\quad r_{13}=r_{23}=1.0000,\qquad r_{12}=2.0000
    3. ニュートンの法則の適用に必要なのは、いま得られた離隔距離だけです。質点3は他の質点からそれぞれ1単位離れており、質点1と質点2の間の距離はその2倍であるため、2体間に働く引力は4分の1に弱まります。

      Fij=Gmimjrij2⟹F12=122=0.2500,F13=F23=1.0000F_{ij}=\frac{Gm_im_j}{r_{ij}^{2}}\quad\Longrightarrow\quad F_{12}=\frac{1}{2^{2}}=0.2500,\qquad F_{13}=F_{23}=1.0000
    4. この同直線上の配置を活用します。質点3は反対方向に等しく引かれるため、その合力は正確にゼロになります。質点1と質点2は、質点3がその間に位置しているため、双方の相手から同じ向きに引かれます。したがって、それら2つの力は成分分解を要せず、単なる数値として足し合わされます。

      ∣F3∣=1.0000−1.0000=0.0000,∣F1∣=1.0000+0.2500=1.2500|\mathbf{F}_3|=1.0000-1.0000=0.0000,\qquad |\mathbf{F}_1|=1.0000+0.2500=1.2500
    5. 全運動量がゼロであり、かつ質点1と質点2が同一の運動をしていることから、質点3の速度は −2v1 に固定されます。これにより質点3の運動エネルギーはそれらの4倍となり、3体全体の総和は単一の項にまとまります。

      ∑mivi=0, v1=v2 ⟹ v3=−2v1,  KE=12(1+1+4)v12=3v12{\textstyle\sum}m_i\mathbf{v}_i=0,\ \mathbf{v}_1=\mathbf{v}_2\ \Longrightarrow\ \mathbf{v}_3=-2\mathbf{v}_1,\ \ KE=\tfrac12(1+1+4)v_1^{2}=3v_1^{2}
    6. v1 の2つの成分を自乗して3倍することが、3体の運動エネルギー計算のすべてです。

      v12=0.4662036852+0.432365732=0.404286⟹KE=1.2129v_1^{2}=0.466203685^{2}+0.43236573^{2}=0.404286\quad\Longrightarrow\quad KE=1.2129
    7. 位置エネルギーの計算にもベクトルは不要です。3つの離隔距離の逆数をとるだけであり、ε = 0 であるため、それらの逆数は厳密な値となります。

      PE=−G(m1m22+m1m31+m2m31)=−2.5000PE=-G\left(\frac{m_1m_2}{2}+\frac{m_1m_3}{1}+\frac{m_2m_3}{1}\right)=-2.5000

    解答

    E = −1.2871、そしてパネルの t = 0 の表のすべての行 — 離隔距離、ペア間の力、合力、両エネルギー、そして全エネルギー — は、それら4つの小数と全運動量の消失のみから導かれました。各部分の変化に注目してください。t = 25.0 までに運動エネルギーは1.305まで増加し、位置エネルギーは−2.592まで低下しますが、全エネルギーは変化していません。この符号こそが決定的な意味を持ちます。束縛された3体系であってもメンバーを失うことはあり得ますが、タダではありません。もし1つの質点が余剰エネルギーなしに離脱した場合、残されたペアは−1.2871全体を担う必要があり、そのエネルギーを持つ2つの単位質量体の軌道半長軸は 1/(2 × 1.2871) = 0.3885 となります。これは、現在それらを隔てている 2.0000 の 1/5 未満です。放出とハードニングは同じ現象の表裏であり、その理由はエネルギー収支にあります。

    E=1.2129−2.5000=−1.2871⟹amax⁡=Gm22∣E∣=12(1.2871)=0.3885E=1.2129-2.5000=-1.2871\quad\Longrightarrow\quad a_{\max}=\frac{Gm^{2}}{2|E|}=\frac{1}{2(1.2871)}=0.3885
参考文献 (1)