MathLabs

幾何学

シンプレクティック幾何学

符号付き面積を測る形式に基づく幾何学で、古典力学の自然な舞台となる。

直観直感:位相空間における符号付き面積

振り子の位置と運動量を平面上の点(q,p)(q,p)として一緒にプロットしたところを想像してほしい。振り子が揺れると、この点は閉じたループを描き、そのループに囲まれた面積は、エネルギーが位置と運動量にどのように配分されても正確に一定のままである。シンプレクティック幾何学とは、この符号付き面積を測る規則であるシンプレクティック形式ω\omegaを備えた空間と、それを正確に保つ写像を研究する分野であり、古典力学の自然な数学的言語である。

位相空間内の可積分ハミルトン系の不変レベル集合を表すトーラスの3Dレンダリング。
位相空間内の不変トーラス:可積分ハミルトン系に対して、リウヴィル-アーノルドの定理は位相空間をこのようなトーラスの葉層に分解し、各トーラスは作用角変数(I,θ)(I,\theta)を持ち、その上で流れは直線的な回転となる。

大学定義:シンプレクティック多様体とハミルトンベクトル場

定義: シンプレクティック多様体

シンプレクティック多様体とは、組(M,ω)(M,\omega)であり、MMは偶数次元2n2nの滑らかな多様体で、ω\omegaは閉じた非退化な22次微分形式である:閉じているとはdω=0d\omega=0を意味し、非退化であるとはゼロでない任意の接ベクトルvvに対してω(v,w)≠0\omega(v,w)\neq 0となる接ベクトルwwが存在することを意味する。

dω=0d\omega = 0

閉性(dω=0d\omega=0)は、ω\omegaに局所的な源が存在しないことを意味し、任意の滑らかなエネルギー関数からハミルトン方程式を無矛盾に回復できるための条件である。非退化性は、ω\omegaが各点で接ベクトルと余接ベクトルの間の線形同型を与えることを意味し、エネルギー関数HHの微分dHdHを、以下の式を通じて実際のベクトル場XHX_Hに変える。

ιXHω=dH\iota_{X_H}\omega = dH

標準座標ω=∑i=1ndqi∧dpi\omega=\sum_{i=1}^n dq_i\wedge dp_iにおいて、定義方程式ιXHω=dH\iota_{X_H}\omega=dHはおなじみのハミルトン方程式q˙i=∂H/∂pi\dot q_i=\partial H/\partial p_i、p˙i=−∂H/∂qi\dot p_i=-\partial H/\partial q_iに展開される:抽象的なシンプレクティック形式論と古典力学の形式論は、表記が異なるだけで同じ主張である。

シンプレクティック幾何学とリーマン幾何学の比較
性質シンプレクティック (M,ω)(M,\omega)リーマン (M,g)(M,g)
双線形形式歪対称:ω(v,w)=−ω(w,v)\omega(v,w)=-\omega(w,v)対称:g(v,w)=g(w,v)g(v,w)=g(w,v)
局所モデル常に同じ:ω=∑i=1ndqi∧dpi\omega=\sum_{i=1}^n dq_i\wedge dp_i(ダルブー)曲率は点ごとに変化しうる
同型写像で保たれるもの面積・体積ωn\omega^nであり、長さや角度ではない長さ、角度、測地距離

大学重要な定理:剛性と保存則

次元2n2nのあらゆるシンプレクティック多様体(M,ω)(M,\omega)は局所的に標準モデルと同型である:各点の周りには座標(q1,…,qn,p1,…,pn)(q_1,\dots,q_n,p_1,\dots,p_n)が存在し、その中でω=∑i=1ndqi∧dpi\omega=\sum_{i=1}^n dq_i\wedge dp_iとなる。特に、シンプレクティック多様体は曲率のような局所不変量を一切持たない。

なぜ正しいのか?

リーマン多様体には曲率という局所不変量があり、ごく小さな近傍でも球面と平面を区別できる。ダルブーの定理は、シンプレクティック多様体にはそのような局所的な指紋が存在しないことを述べている:近くで見れば、あらゆるシンプレクティック多様体は標準形式を持つ平坦な位相空間R2n\mathbb{R}^{2n}とまったく同じに見える。したがって、シンプレクティック位相幾何学の興味深い内容はすべて大域的なものである。

証明

ステップ1(線形問題の設定)。線形座標変換によって、まず選んだ点p0p_0で形式ω\omegaが標準形式ω0\omega_0にちょうどp0p_0自身で一致するようにする。これは、R2n\mathbb{R}^{2n}上のあらゆる非退化な歪対称双線形形式が適切な基底により標準形に変換できるという線形代数の事実による。

ステップ2(補間する)。族ωt=(1−t)ω0+tω\omega_t=(1-t)\omega_0+t\omegaをt∈[0,1]t\in[0,1]に対して定義する。ω0\omega_0とω\omegaはp0p_0で一致し、そこで共に閉かつ非退化なので、ωt\omega_tも閉であり、必要なら近傍を縮小することで、すべてのttについて非退化となる。

ステップ3(モーザー方程式を解く)。ω\omega引くω0\omega_0は閉でp0p_0で消えるので、ポアンカレの補題によりdσ=ω−ω0d\sigma=\omega-\omega_0を満たす1次形式σ\sigmaが得られる。モーザーのトリックは、ιXtωt=−σ\iota_{X_t}\omega_t=-\sigmaを解く時間依存ベクトル場XtX_tを求めるものであり、ωt\omega_tが非退化であることからこれは各点で一意にXtX_tを定める。

ステップ4(流れを積分する)。XtX_tが生成する流れをφt\varphi_tとする。カルタンの公式による直接計算からLXtωt=dιXtωt+ιXtdωt=−dσ=−(ω−ω0)\mathcal{L}_{X_t}\omega_t=d\iota_{X_t}\omega_t+\iota_{X_t}d\omega_t=-d\sigma=-(\omega-\omega_0)が得られ、ω˙t=ω−ω0\dot\omega_t=\omega-\omega_0と合わせるとddt(φt∗ωt)=0\frac{d}{dt}(\varphi_t^*\omega_t)=0となる:引き戻しφt∗ωt\varphi_t^*\omega_tはttに関して一定であり、すべてのttについてω0\omega_0に等しい。

ステップ5(結論)。t=1t=1とおくと、縮小した近傍上でφ1∗ω=ω0\varphi_1^*\omega=\omega_0が得られ、φ1\varphi_1に沿って引き戻した座標が求めるダルブー座標である。

(M,ω)(M,\omega)を次元2n2nのシンプレクティック多様体とし、XHX_HをιXHω=dH\iota_{X_H}\omega=dHで定義されるHHのハミルトンベクトル場とし、ϕt\phi_tをそのフローとする。このときϕt\phi_tはシンプレクティック体積形式を保つ:すべてのttに対してϕt∗(ωn)=ωn\phi_t^{*}(\omega^n)=\omega^n。特に、ハミルトンフローは位相空間の体積を保存する。

なぜ正しいのか?

古典力学において、リウヴィルの定理は、ハミルトン方程式に従う粒子の気体が位相空間のより小さな領域へ自発的に圧縮できない数学的な理由である:流れはその領域を引き伸ばしたりねじったりして奇妙な形にすることはあっても、その体積が縮んだり増えたりすることは決してない。これは統計力学の基盤であり、この保存則を尊重しない素朴な数値積分法が長時間のシミュレーションで見かけ上のエネルギードリフトを生む理由を説明する。

証明

ステップ1(カルタンの魔法公式)。任意のベクトル場XXと形式ω\omegaに対して、カルタンの公式はLXω=d(ιXω)+ιX(dω)\mathcal{L}_X\omega=d(\iota_X\omega)+\iota_X(d\omega)を与える。XXをXHX_Hとして適用すると、ω\omegaは閉なのでdω=0d\omega=0であり、ハミルトンベクトル場の定義によりιXHω=dH\iota_{X_H}\omega=dHなので、LXHω=d(dH)+0\mathcal{L}_{X_H}\omega=d(dH)+0が得られる。

ステップ2(形式自体は変化しない)。外微分はすべての関数に対してd∘d=0d\circ d=0を満たすので、d(dH)=0d(dH)=0である。したがってLXHω=0\mathcal{L}_{X_H}\omega=0となる:ハミルトンフローは体積だけでなくシンプレクティック形式そのものを保つ。

ステップ3(最高次のべきへ移る)。体積形式はωn=ω∧⋯∧ω\omega^n=\omega\wedge\cdots\wedge\omega(nn個の因子)である。ウェッジ積上のリー微分に対するライプニッツ則からLXH(ωn)=n ωn−1∧LXHω\mathcal{L}_{X_H}(\omega^n)=n\,\omega^{n-1}\wedge\mathcal{L}_{X_H}\omegaが得られ、ステップ2によりLXHω=0\mathcal{L}_{X_H}\omega=0なので右辺はゼロになる:LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n)=0。

ステップ4(流れに沿って積分する)。ϕt\phi_tがXHX_Hの流れであれば、LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n)=0はまさにすべてのttについてddtϕt∗(ωn)=0\frac{d}{dt}\phi_t^{*}(\omega^n)=0を意味する。ϕ0=id\phi_0=\mathrm{id}からϕ0∗(ωn)=ωn\phi_0^{*}(\omega^n)=\omega^nが成り立つので、積分するとすべてのttについてϕt∗(ωn)=ωn\phi_t^{*}(\omega^n)=\omega^nが示され、これがリウヴィルの定理である。

発展発展:シンプレクティックな剛性とグロモフの非押し込み定理

グロモフは1985年に、シンプレクティック球B2n(r)={∣z∣≤r}B^{2n}(r)=\{|z|\le r\}が円柱B2(R)×R2n−2B^2(R)\times\mathbb{R}^{2n-2}へシンプレクティックに埋め込めるのはr≤Rr\le Rのとき、かつそのときに限ることを証明した——これはnnがどれほど大きくても、半径rrの22次元円板を半径RRの円板へ埋め込む条件とまったく同じである。この証明は擬正則曲線を用いており、シンプレクティック容量と呼ばれる不変量を生み出した。その最も単純なものがグロモフ幅c(B2n(r))=πr2c(B^{2n}(r))=\pi r^2であり、球を押し込める最小の円板因子の面積である。

c(B2n(r))=πr2c(B^{2n}(r)) = \pi r^2

大学実世界での応用と具体例

ハミルトン力学はシンプレクティックな言語で書かれているため、シンプレクティック幾何学はエネルギーと位相空間構造を保存する物理系のあらゆる場面に現れる:天体力学ではそれを用いて惑星軌道の長期安定性を研究し、粒子加速器の設計ではシンプレクティック写像を用いて、人為的な減衰や増大なしに数百万周にわたってビームを追跡し、分子動力学ソフトウェアではシンプレクティック積分法を用いて、シミュレートされた分子が数値誤差だけでエネルギーを得たり失ったりしないようにする。最適制御理論や数理物理学における幾何学的量子化もまた、シンプレクティック形式論の上に直接築かれている。

例: 回転がシンプレクティック写像であることの確認

(R2,ω=dq∧dp)(\mathbb{R}^2,\omega=dq\wedge dp)上で、固定角θ\thetaによる回転ϕθ(q,p)=(qcos⁡θ−psin⁡θ, qsin⁡θ+pcos⁡θ)\phi_\theta(q,p)=(q\cos\theta-p\sin\theta,\ q\sin\theta+p\cos\theta)を考える。ϕθ\phi_\thetaはシンプレクトモルフィズムか、すなわちϕθ∗ω=ω\phi_\theta^{*}\omega=\omegaを満たすか。

解答

ステップ1:ヤコビ行列を書く。この写像は線形であり、行列J=(cos⁡θ−sin⁡θsin⁡θcos⁡θ)J=\begin{pmatrix}\cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\end{pmatrix}が(q,p)(q,p)に作用する。

ステップ2:形式を引き戻す。R2\mathbb{R}^2上の線形写像に対してはϕθ∗(dq∧dp)=det⁡(J) dq∧dp\phi_\theta^{*}(dq\wedge dp)=\det(J)\,dq\wedge dpが成り立つので、JJの行列式を計算すればよい。

ステップ3:行列式を計算する。ピタゴラスの恒等式により、すべてのθ\thetaに対してdet⁡(J)=cos⁡2θ+sin⁡2θ=1\det(J)=\cos^2\theta+\sin^2\theta=1である。

ステップ4:結論。det⁡(J)=1\det(J)=1なので、すべてのθ\thetaに対してϕθ∗ω=ω\phi_\theta^{*}\omega=\omegaとなる:回転はシンプレクトモルフィズムであり、これは回転が面積を保つという幾何学的事実と正確に一致する。

例: 調和振動子に対するシンプレクティックオイラー法と陽的オイラー法の比較

調和振動子H(q,p)=12(p2+q2)H(q,p)=\tfrac12(p^2+q^2)について、シンプレクティックオイラー法の1ステップqn+1=qn+hpn, pn+1=pn−hqn+1q_{n+1}=q_n+hp_n,\ p_{n+1}=p_n-hq_{n+1}と、陽的(素朴な)オイラー法の1ステップqn+1=qn+hpn, pn+1=pn−hqnq_{n+1}=q_n+hp_n,\ p_{n+1}=p_n-hq_nを、小さいステップ幅h>0h>0で比較せよ。どちらがリウヴィルの定理と整合するか。

解答

ステップ1:シンプレクティックオイラー写像を(qn,pn)(q_n,p_n)の関数として書く。qn+1=qn+hpnq_{n+1}=q_n+hp_nを第二式に代入するとpn+1=pn−h(qn+hpn)=pn−hqn−h2pnp_{n+1}=p_n-h(q_n+hp_n)=p_n-hq_n-h^2p_nが得られる。

ステップ2:そのヤコビ行列を計算する。∂(qn+1,pn+1)∂(qn,pn)=(1h−h1−h2)\frac{\partial(q_{n+1},p_{n+1})}{\partial(q_n,p_n)}=\begin{pmatrix}1&h\\ -h&1-h^2\end{pmatrix}であり、行列式は1⋅(1−h2)−h⋅(−h)=1−h2+h2=11\cdot(1-h^2)-h\cdot(-h)=1-h^2+h^2=1となり、すべてのhhについて成り立つ。

ステップ3:陽的オイラー法のヤコビ行列を計算する。そこでは∂(qn+1,pn+1)∂(qn,pn)=(1h−h1)\frac{\partial(q_{n+1},p_{n+1})}{\partial(q_n,p_n)}=\begin{pmatrix}1&h\\ -h&1\end{pmatrix}であり、行列式は1+h21+h^2で、h≠0h\neq0である限り常に11より大きい。

ステップ4:解釈する。リウヴィルの定理により、厳密な流れは各ステップでヤコビ行列式がちょうど11である。シンプレクティックオイラー法はhhによらずこれと正確に一致するため、位相空間の体積を人為的に膨張させない。一方、陽的オイラー法は各ステップで面積を係数1+h21+h^2だけ拡大するため、多くのステップにわたってエネルギーが上方に漂う——これが長時間シミュレーションでシンプレクティック積分法が使われる理由である。

なぜシンプレクティック多様体は偶数次元2n2nでなければならないのか。

調和振動子H(q,p)=12(p2+q2)H(q,p)=\tfrac12(p^2+q^2)と標準形式ω=dq∧dp\omega=dq\wedge dpに対して、ιXHω=dH\iota_{X_H}\omega=dHで定義されるハミルトンベクトル場XH=(q˙,p˙)X_H=(\dot q,\dot p)は何か。

リウヴィルの定理によれば、シンプレクティック多様体上のハミルトンフローは何を保存するか。

太陽系の数百万年にわたる長期シミュレーションでは、標準的なルンゲ・クッタ法のような汎用手法ではなくシンプレクティック積分法が用いられる。なぜか。

参考文献

  1. Dusa McDuff, Dietmar Salamon (2017). Introduction to Symplectic Topology
  2. Mikhail Gromov (1985). Pseudo holomorphic curves in symplectic manifolds
  3. Pazit Haim-Kislev, Yaron Ostrover (2026). A Counterexample to Viterbo's Conjecture · arXiv:2405.16513