MathLabs

几何学

辛几何

基于测量带符号面积的形式所构建的几何学,是经典力学的自然舞台。

直观直觉:相空间中的带符号面积

设想将摆动摆的位置和动量一起画成平面上一点(q,p)(q,p)的运动。摆动时,这一点描出一条闭合曲线,而无论能量在位置与动量之间如何分配,该曲线所围的面积始终保持不变。辛几何研究配备了辛形式ω\omega——一种测量这种带符号面积的规则——的空间,以及精确保持它的映射,这正是经典力学的自然数学语言。

表示相空间中可积哈密顿系统不变水平集的环面三维渲染。
相空间中的不变环面:对于可积哈密顿系统,刘维尔-阿诺德定理将相空间叶状分解为这样的环面,每个环面都带有作用-角坐标(I,θ)(I,\theta),在其上流是直线旋转。

大学定义:辛流形与哈密顿向量场

定义: 辛流形

辛流形是一对(M,ω)(M,\omega),其中MM是偶数维2n2n的光滑流形,ω\omega是一个闭的、非退化的22次微分形式:闭指dω=0d\omega=0,非退化指对任意非零切向量vv,都存在切向量ww使得ω(v,w)≠0\omega(v,w)\neq 0。

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}完全相同。因此,辛拓扑中所有有趣的内容都是整体性的。

证明

第一步(建立线性问题)。通过线性坐标变换,首先安排在所选点p0p_0处,形式ω\omega恰好在p0p_0本身处等于标准形式ω0\omega_0,这是因为R2n\mathbb{R}^{2n}上任何非退化的斜对称双线性形式都可以通过合适的基化为标准形式。

第二步(插值)。对t∈[0,1]t\in[0,1]定义族ωt=(1−t)ω0+tω\omega_t=(1-t)\omega_0+t\omega。由于ω0\omega_0与ω\omega在p0p_0处相同,且两者在该处都是闭的、非退化的,故ωt\omega_t也是闭的,并且在必要时缩小邻域后,对每个tt都非退化。

第三步(求解莫泽方程)。由于ω\omega减去ω0\omega_0是闭的且在p0p_0处为零,由庞加莱引理可得一个1-形式σ\sigma满足dσ=ω−ω0d\sigma=\omega-\omega_0。莫泽技巧寻找一个随时间变化的向量场XtX_t求解ιXtωt=−σ\iota_{X_t}\omega_t=-\sigma;由于ωt\omega_t非退化,这在每一点唯一确定XtX_t。

第四步(对流积分)。设φt\varphi_t为XtX_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。

第五步(得出结论)。取t=1t=1可得在缩小后的邻域上φ1∗ω=ω0\varphi_1^*\omega=\omega_0,因此沿φ1\varphi_1拉回的坐标即为所求的达布坐标。

设XHX_H是辛流形(M,ω)(M,\omega)(维数为2n2n)上HH的哈密顿向量场,由ιXHω=dH\iota_{X_H}\omega=dH定义,设ϕt\phi_t为其流。则ϕt\phi_t保持辛体积形式:对任意tt都有ϕt∗(ωn)=ωn\phi_t^{*}(\omega^n)=\omega^n。特别地,哈密顿流保持相空间体积。

为什么成立?

在经典力学中,刘维尔定理是遵循哈密顿方程的粒子气体无法自发压缩到相空间更小区域的数学原因:流可以将该区域拉伸、扭曲成奇怪的形状,但其体积绝不会缩小或增大。这是统计力学的基础,也解释了为何不尊重此守恒律的朴素数值积分法会在长时间模拟中产生虚假的能量漂移。

证明

第一步(嘉当魔法公式)。对任意向量场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。

第二步(形式本身不变)。外微分对任意函数都满足d∘d=0d\circ d=0,故d(dH)=0d(dH)=0。因此LXHω=0\mathcal{L}_{X_H}\omega=0:哈密顿流不仅保持体积,还保持辛形式本身。

第三步(转到最高幂)。体积形式为ω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,由第二步LXHω=0\mathcal{L}_{X_H}\omega=0,右边为零:LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n)=0。

第四步(沿流积分)。若ϕ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——这与将半径为rr的22维圆盘嵌入半径为RR的圆盘的条件完全相同,无论nn有多大。该证明使用了伪全纯曲线,并催生了称为辛容量的不变量,其中最简单的是格罗莫夫宽度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?

解答

第一步:写出雅可比矩阵。该映射是线性的,矩阵为J=(cos⁡θ−sin⁡θsin⁡θcos⁡θ)J=\begin{pmatrix}\cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\end{pmatrix},作用于(q,p)(q,p)。

第二步:将形式拉回。对于R2\mathbb{R}^2上的线性映射,有ϕθ∗(dq∧dp)=det⁡(J) dq∧dp\phi_\theta^{*}(dq\wedge dp)=\det(J)\,dq\wedge dp,因此只需计算JJ的行列式。

第三步:计算行列式。由勾股恒等式,对任意θ\theta都有det⁡(J)=cos⁡2θ+sin⁡2θ=1\det(J)=\cos^2\theta+\sin^2\theta=1。

第四步:得出结论。由于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),比较辛欧拉法的一步qn+1=qn+hpn, pn+1=pn−hqn+1q_{n+1}=q_n+hp_n,\ p_{n+1}=p_n-hq_{n+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。哪一种方法与刘维尔定理相符?

解答

第一步:将辛欧拉映射写成(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。

第二步:计算其雅可比矩阵。∂(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都成立。

第三步:计算显式欧拉法的雅可比矩阵。此时∂(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。

第四步:解释。由刘维尔定理,精确的流在每一步的雅可比行列式恰好为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