幾何学
シンプレクティック幾何学
符号付き面積を測る形式に基づく幾何学で、古典力学の自然な舞台となる。
直観直感:位相空間における符号付き面積
振り子の位置と運動量を平面上の点(q,p)として一緒にプロットしたところを想像してほしい。振り子が揺れると、この点は閉じたループを描き、そのループに囲まれた面積は、エネルギーが位置と運動量にどのように配分されても正確に一定のままである。シンプレクティック幾何学とは、この符号付き面積を測る規則であるシンプレクティック形式ωを備えた空間と、それを正確に保つ写像を研究する分野であり、古典力学の自然な数学的言語である。
位相空間内の不変トーラス:可積分ハミルトン系に対して、リウヴィル-アーノルドの定理は位相空間をこのようなトーラスの葉層に分解し、各トーラスは作用角変数(I,θ)を持ち、その上で流れは直線的な回転となる。大学定義:シンプレクティック多様体とハミルトンベクトル場
定義: シンプレクティック多様体
シンプレクティック多様体とは、組(M,ω)であり、Mは偶数次元2nの滑らかな多様体で、ωは閉じた非退化な2次微分形式である:閉じているとはdω=0を意味し、非退化であるとはゼロでない任意の接ベクトルvに対してω(v,w)=0となる接ベクトルwが存在することを意味する。
閉性(dω=0)は、ωに局所的な源が存在しないことを意味し、任意の滑らかなエネルギー関数からハミルトン方程式を無矛盾に回復できるための条件である。非退化性は、ωが各点で接ベクトルと余接ベクトルの間の線形同型を与えることを意味し、エネルギー関数Hの微分dHを、以下の式を通じて実際のベクトル場XHに変える。
ιXHω=dH 標準座標ω=∑i=1ndqi∧dpiにおいて、定義方程式ιXHω=dHはおなじみのハミルトン方程式q˙i=∂H/∂pi、p˙i=−∂H/∂qiに展開される:抽象的なシンプレクティック形式論と古典力学の形式論は、表記が異なるだけで同じ主張である。
シンプレクティック幾何学とリーマン幾何学の比較| 性質 | シンプレクティック (M,ω) | リーマン (M,g) |
|---|
| 双線形形式 | 歪対称:ω(v,w)=−ω(w,v) | 対称:g(v,w)=g(w,v) |
| 局所モデル | 常に同じ:ω=∑i=1ndqi∧dpi(ダルブー) | 曲率は点ごとに変化しうる |
| 同型写像で保たれるもの | 面積・体積ωnであり、長さや角度ではない | 長さ、角度、測地距離 |
大学重要な定理:剛性と保存則
次元2nのあらゆるシンプレクティック多様体(M,ω)は局所的に標準モデルと同型である:各点の周りには座標(q1,…,qn,p1,…,pn)が存在し、その中でω=∑i=1ndqi∧dpiとなる。特に、シンプレクティック多様体は曲率のような局所不変量を一切持たない。
なぜ正しいのか?
リーマン多様体には曲率という局所不変量があり、ごく小さな近傍でも球面と平面を区別できる。ダルブーの定理は、シンプレクティック多様体にはそのような局所的な指紋が存在しないことを述べている:近くで見れば、あらゆるシンプレクティック多様体は標準形式を持つ平坦な位相空間R2nとまったく同じに見える。したがって、シンプレクティック位相幾何学の興味深い内容はすべて大域的なものである。
証明
ステップ1(線形問題の設定)。線形座標変換によって、まず選んだ点p0で形式ωが標準形式ω0にちょうどp0自身で一致するようにする。これは、R2n上のあらゆる非退化な歪対称双線形形式が適切な基底により標準形に変換できるという線形代数の事実による。
ステップ2(補間する)。族ωt=(1−t)ω0+tωをt∈[0,1]に対して定義する。ω0とωはp0で一致し、そこで共に閉かつ非退化なので、ωtも閉であり、必要なら近傍を縮小することで、すべてのtについて非退化となる。
ステップ3(モーザー方程式を解く)。ω引くω0は閉でp0で消えるので、ポアンカレの補題によりdσ=ω−ω0を満たす1次形式σが得られる。モーザーのトリックは、ιXtωt=−σを解く時間依存ベクトル場Xtを求めるものであり、ωtが非退化であることからこれは各点で一意にXtを定める。
ステップ4(流れを積分する)。Xtが生成する流れをφtとする。カルタンの公式による直接計算からLXtωt=dιXtωt+ιXtdωt=−dσ=−(ω−ω0)が得られ、ω˙t=ω−ω0と合わせるとdtd(φt∗ωt)=0となる:引き戻しφt∗ωtはtに関して一定であり、すべてのtについてω0に等しい。
ステップ5(結論)。t=1とおくと、縮小した近傍上でφ1∗ω=ω0が得られ、φ1に沿って引き戻した座標が求めるダルブー座標である。
(M,ω)を次元2nのシンプレクティック多様体とし、XHをιXHω=dHで定義されるHのハミルトンベクトル場とし、ϕtをそのフローとする。このときϕtはシンプレクティック体積形式を保つ:すべてのtに対してϕt∗(ωn)=ωn。特に、ハミルトンフローは位相空間の体積を保存する。
なぜ正しいのか?
古典力学において、リウヴィルの定理は、ハミルトン方程式に従う粒子の気体が位相空間のより小さな領域へ自発的に圧縮できない数学的な理由である:流れはその領域を引き伸ばしたりねじったりして奇妙な形にすることはあっても、その体積が縮んだり増えたりすることは決してない。これは統計力学の基盤であり、この保存則を尊重しない素朴な数値積分法が長時間のシミュレーションで見かけ上のエネルギードリフトを生む理由を説明する。
証明
ステップ1(カルタンの魔法公式)。任意のベクトル場Xと形式ωに対して、カルタンの公式はLXω=d(ιXω)+ιX(dω)を与える。XをXHとして適用すると、ωは閉なのでdω=0であり、ハミルトンベクトル場の定義によりιXHω=dHなので、LXHω=d(dH)+0が得られる。
ステップ2(形式自体は変化しない)。外微分はすべての関数に対してd∘d=0を満たすので、d(dH)=0である。したがってLXHω=0となる:ハミルトンフローは体積だけでなくシンプレクティック形式そのものを保つ。
ステップ3(最高次のべきへ移る)。体積形式はωn=ω∧⋯∧ω(n個の因子)である。ウェッジ積上のリー微分に対するライプニッツ則からLXH(ωn)=nωn−1∧LXHωが得られ、ステップ2によりLXHω=0なので右辺はゼロになる:LXH(ωn)=0。
ステップ4(流れに沿って積分する)。ϕtがXHの流れであれば、LXH(ωn)=0はまさにすべてのtについてdtdϕt∗(ωn)=0を意味する。ϕ0=idからϕ0∗(ωn)=ωnが成り立つので、積分するとすべてのtについてϕt∗(ωn)=ωnが示され、これがリウヴィルの定理である。
発展発展:シンプレクティックな剛性とグロモフの非押し込み定理
グロモフは1985年に、シンプレクティック球B2n(r)={∣z∣≤r}が円柱B2(R)×R2n−2へシンプレクティックに埋め込めるのはr≤Rのとき、かつそのときに限ることを証明した——これはnがどれほど大きくても、半径rの2次元円板を半径Rの円板へ埋め込む条件とまったく同じである。この証明は擬正則曲線を用いており、シンプレクティック容量と呼ばれる不変量を生み出した。その最も単純なものがグロモフ幅c(B2n(r))=πr2であり、球を押し込める最小の円板因子の面積である。
c(B2n(r))=πr2 大学実世界での応用と具体例
ハミルトン力学はシンプレクティックな言語で書かれているため、シンプレクティック幾何学はエネルギーと位相空間構造を保存する物理系のあらゆる場面に現れる:天体力学ではそれを用いて惑星軌道の長期安定性を研究し、粒子加速器の設計ではシンプレクティック写像を用いて、人為的な減衰や増大なしに数百万周にわたってビームを追跡し、分子動力学ソフトウェアではシンプレクティック積分法を用いて、シミュレートされた分子が数値誤差だけでエネルギーを得たり失ったりしないようにする。最適制御理論や数理物理学における幾何学的量子化もまた、シンプレクティック形式論の上に直接築かれている。
例: 回転がシンプレクティック写像であることの確認
(R2,ω=dq∧dp)上で、固定角θによる回転ϕθ(q,p)=(qcosθ−psinθ, qsinθ+pcosθ)を考える。ϕθはシンプレクトモルフィズムか、すなわちϕθ∗ω=ωを満たすか。
解答
ステップ1:ヤコビ行列を書く。この写像は線形であり、行列J=(cosθsinθ−sinθcosθ)が(q,p)に作用する。
ステップ2:形式を引き戻す。R2上の線形写像に対してはϕθ∗(dq∧dp)=det(J)dq∧dpが成り立つので、Jの行列式を計算すればよい。
ステップ3:行列式を計算する。ピタゴラスの恒等式により、すべてのθに対してdet(J)=cos2θ+sin2θ=1である。
ステップ4:結論。det(J)=1なので、すべてのθに対してϕθ∗ω=ωとなる:回転はシンプレクトモルフィズムであり、これは回転が面積を保つという幾何学的事実と正確に一致する。
例: 調和振動子に対するシンプレクティックオイラー法と陽的オイラー法の比較
調和振動子H(q,p)=21(p2+q2)について、シンプレクティックオイラー法の1ステップqn+1=qn+hpn, pn+1=pn−hqn+1と、陽的(素朴な)オイラー法の1ステップqn+1=qn+hpn, pn+1=pn−hqnを、小さいステップ幅h>0で比較せよ。どちらがリウヴィルの定理と整合するか。
解答
ステップ1:シンプレクティックオイラー写像を(qn,pn)の関数として書く。qn+1=qn+hpnを第二式に代入するとpn+1=pn−h(qn+hpn)=pn−hqn−h2pnが得られる。
ステップ2:そのヤコビ行列を計算する。∂(qn,pn)∂(qn+1,pn+1)=(1−hh1−h2)であり、行列式は1⋅(1−h2)−h⋅(−h)=1−h2+h2=1となり、すべてのhについて成り立つ。
ステップ3:陽的オイラー法のヤコビ行列を計算する。そこでは∂(qn,pn)∂(qn+1,pn+1)=(1−hh1)であり、行列式は1+h2で、h=0である限り常に1より大きい。
ステップ4:解釈する。リウヴィルの定理により、厳密な流れは各ステップでヤコビ行列式がちょうど1である。シンプレクティックオイラー法はhによらずこれと正確に一致するため、位相空間の体積を人為的に膨張させない。一方、陽的オイラー法は各ステップで面積を係数1+h2だけ拡大するため、多くのステップにわたってエネルギーが上方に漂う——これが長時間シミュレーションでシンプレクティック積分法が使われる理由である。
なぜシンプレクティック多様体は偶数次元2nでなければならないのか。
調和振動子H(q,p)=21(p2+q2)と標準形式ω=dq∧dpに対して、ιXHω=dHで定義されるハミルトンベクトル場XH=(q˙,p˙)は何か。
リウヴィルの定理によれば、シンプレクティック多様体上のハミルトンフローは何を保存するか。
太陽系の数百万年にわたる長期シミュレーションでは、標準的なルンゲ・クッタ法のような汎用手法ではなくシンプレクティック積分法が用いられる。なぜか。