「一般解なし」を絶対的なカオスと混同してはならない
劉慈欣の小説『三体』の設定や世間の一般的なイメージにおいて、三体問題は往々にして「絶対的な予測不可能」の代名詞とされてきた。しかしこれは100年以上にわたる認識の誤解である。1880年代、アンリ・ポアンカレ（Henri Poincaré）は三体系において「一般的な解析解（代数的な初等関数による第一積分）が存在しない」ことを証明し、この結論が現代のカオス理論を生み出す契機となった。だが、普遍的な一般解の公式が存在しないからといって、カオス系の中に厳密な周期解が存在しないわけではない。
「Three Body Orbits」と名付けられたWebサイトは、13本の学術論文に散らばっていた3,915本もの厳密な周期解を1つのインタラクティブな図鑑として集約した。等しい質量を持つ3つの天体が、綿密に計算された初期条件のもとで互いに周回し、1周期を終えると寸分の狂いもなく元の位置と速度に戻ってくる。**「理論数学において一般解の公式が存在しないこと」と「工学的な数値計算において予測が不可能なこと」はまったく別の問題である。**十分な計算パワーさえあれば、カオスの大海原のなかにも秩序の孤島を描き出すことができるのだ。
図：Three Body Orbitsのトップページ全体像。各周期軌道が星図のように配置され、似た性質の軌道が島を形成している。出典：Three Body Orbitsサイト画面
32桁演算の極限に挑み、ブラウザ上のリアルタイムシミュレーションを断念
図鑑に収録された3,915本の軌道の背景には、極めて過酷な工学的再構築のプロセスが存在する。制作者は、論文に記載された初期条件をそのまま転記することはしなかった。学術論文は組版上の制約から、通常は有効数字を10〜15桁程度しか記載しない。しかし三体系のようなカオス系においては、たとえ10⁻¹⁵レベルの極小の丸め誤差であっても数周期のうちに指数関数的に増大し、天体の軌道は瞬く間に破綻して崩壊してしまう。
各軌道を始点と終点で継ぎ目なく閉じるため、制作者は高精度任意精度演算ライブラリとニュートン法（ニュートン・ラフソン法）を用い、微分方程式のポアンカレ断面写像（Poincaré section map）を再収束させた。32桁の十進高精度で初期位置と速度を反復補正した結果、軌道の閉合誤差の中央値を1.4×10⁻²⁴という驚異的な数値にまで抑え込むことに成功した。このスケールにおいて、天体が1周期を周回し終えた後の位置誤差は、原子1個の直径よりもはるかに小さい。収録された軌道のうち3,889本は10⁻⁹より優れた閉合精度を達成している。
この極限の精度によって得られた軌道は、フロントエンドの実装において一つのジレンマをもたらした。もし「ブラウザ上のリアルタイム数値シミュレーション」を選択し、JavaScriptのネイティブな64ビット倍精度浮動小数点数でクライアント側で実時間積分を行った場合、本来閉じるはずの軌道は浮動小数点誤差の蓄積によって数秒でドリフトを起こしてしまう。
そこでサイトは、「事前計算スプライン＋ブラウザ再生」というアーキテクチャを採用した。サーバー側で超高精度積分された真の座標データを、970万点を超える高次多項式スプラインのノット（制御点）として圧縮し、クライアントにストリーミング配信する。フロントエンドはタイムスタンプに基づいてこれらのノットを補間し、滑らかにレンダリングするだけで済む。事前計算によるストレージ容量をランタイムの正確性と引き換えにし、重い計算をバックエンドに隔離することで、ブラウザのリソースを毎秒60フレームの滑らかな描画へと全振りしたのである。
図：Figure-8（8の字形）軌道の表示画面。3つの天体が長時間露光の美しい軌跡を描いて交差する。出典：Three Body Orbitsサイト画面
三世代の研究者がバトンを繋いだ3,915本の厳密解
この3,915本の軌道には、過去数十年にわたる非線形力学の学術的な挑戦の歴史が凝縮されている。これは単一の研究室による成果ではなく、世界各国の研究チームがスーパーコンピュータを駆使してリレーのように繋いできた分類マップなのだ。
上海交通大学の廖世俊（Shijun Liao）教授率いるチーム（Li–Liaoグループ）が最大の主力であり、全体の過半数となる1,964本の軌道を提供している。彼らはスパコンと独自に開発した高精度数値解析手法「CNS（Clean Numerical Simulation）」を用い、未知のトポロジー構造を体系的に洗い出した。ブルガリアのソフィア大学（Sofia）チームは709本の複雑な軌道を報告し、高次高調波ファミリーを大幅に拡充した。数学界で名高い「コレオグラフィー（Choreographies、舞踏系）」は583本を占める。ここでは3つの天体がまるで同じワイヤーの上でジャグリングするように、単一の同一軌道曲線を共有しながら追随し合う。
さらに、自由落下の特徴を持つ311本の特殊解（Free-fall）、セルビアのベオグラード（Belgrade）研究グループ（Šuvakov and Dmitrašinović）が発見した幾何学的に美しい230本の解も含まれる。リストの中には、Roger Brouckeが1975年に早くも導出していた52本の古典的軌道も収録されている。これらに本サイトが分散コンピューティングで新たに発見した66本の軌道が加わり、包括的な力学カタログが完成した。座標の数値一つひとつが、スパコンクラスタで燃やされた数十万時間の計算リソースを物語っている。
スライダーの「キック」で試すFloquet安定限界
図鑑内の1,593本の軌道は「線形安定（Linearly Stable）」と明記されており、全体の40.7%を占める。力学系において周期軌道が外乱に耐えうるかを判断する核心的な指標が「フロケ乗数（Floquet multiplier）」である。これは系に微小な摂動（ノイズ）が加わった際、1周期を通過した後にその誤差がどのように成長するかを測定するものだ。
フロケ乗数の絶対値（大きさ）が1より大きい場合、誤差は指数関数的に累積し、系は瞬く間に崩壊する。一方で最大フロケ乗数が厳密に1.000である場合、系は安定境界にとどまる。この1.000は、摂動が周期的な振動を引き起こすにとどまり、軌道全体のトポロジーが維持されることを意味する。フロケ乗数が1.000であることこそ、カオスの荒波の中でこれらの軌道が自己の存在を繋ぎ止める数学的証拠なのだ。
安定軌道と不安定軌道の力学的な挙動の差を直感的に体験できるよう、サイトの右側には「Kick（キック）」スライダーが用意されている。ブラウザ上でいずれかの天体の速度ベクトルを特定方向にわずかにずらすと、フロントエンドのリアルタイム数値積分器が即座に介入し、安全な事前計算スプラインからリアルタイム演算へと切り替わる。
Broucke R12のような線形安定軌道であれば、小さなキックを与えても天体は元の軌道周辺で少し揺らめくだけで、自己の復元力によってバランスを取り戻し、三体の共舞を維持し続ける。しかし不安定な軌道にスライダーで軽くキックを加えると、三星系は数周期のうちに崩壊し、天体はランダムに宇宙の深淵へと放り出されてしまう。**このミニマルなインタラクティブ設計は、難解な微分方程式の固有値解析を、指先で触れられる直感的な物理体験へと昇華させた。**Hacker Newsの議論スレッドでは、物理学者やエンジニアが摂動の限界値について白热した議論を繰り広げている。
パラメータ探索を分散型の草の根クラウドソーシングへ
さらに最近、サイトはブラウザベースの計算能力提供機能「Huntモード」を公開した。ユーザーがブラウザでこのモードを有効にすると、指定されたパラメータ空間内で未知の軌道を探索する分散コンピューティングに参加できる。もしユーザーの端末が新たな線形安定周期解を収束・発見した場合、その発見者はアルゴリズム惑星系に自ら命名することができる。退屈な初期値の探索グリッドが、分散型のギーク向けクラウドソーシング・ゲームへと変貌を遂げたのだ。
世間一般では「三体問題には解がない」ことが絶対的な運命のように語られがちだ。しかしこの図鑑が示すのは、小数点以下24桁の精度で閉合解を検証できる現代の計算力があれば、カオスの中にも4,000本近い厳密な数学的秩序が存在するという事実である。三体問題はもはや哲学的な袋小路ではなく、具体的なデータと怒涛の計算力によって切り拓かれる、活気に満ちたエンジニアリングの実験場へと進化したのである。
参考リンク：
- Three Body Orbits 公式サイト
- Physical Review Letters 110, 114301
- Science China Physics, Mechanics & Astronomy 60, 129511
- Sofia University Research Group Publications
- Celestial Mechanics 12, 439-462
- Hacker Newsディスカッション (item?id=49670852)
