市面上有许多介绍Hamilton力学的书籍,无论它们来自于数学家还是物理学家。通常开篇就指出相空间是一个辛流形,接着讨论动力学过程如何用微分几何的语言描述,详细介绍如何从Lagrange力学恰当地变换到Hamilton力学,以及二者之间的联系的书籍却少有。本篇文章是关于如何从Lagrange力学下的位形流形过渡到Hamilton力学下的相流形的介绍,希望能架起这座桥桥梁。更细致深刻的讨论可以参看 Marsden 的 Introduction to Mechanics and Symmetry 及其中译版。
符号约定
由于物理学的惯用符号和微分几何学中的惯用符号有冲突,为避免混淆,我们先作符号约定
物理量
符号
说明
位形流形
M M M
力学系统的状态空间
位形流形上的点
q q q
q i : M → R q^i:M \to \mathbb{R} q i : M → R 为局部坐标
切向量/广义速度
v = v i ∂ ∂ q i v=v^i \frac{\partial}{\partial q^i} v = v i ∂ q i ∂
v i v^i v i 为切向量在局部坐标架下的分量
余切向量/广义动量
p = p i d q i p=p_i \mathrm{d}q^i p = p i d q i
p i p_i p i 为余切向量在局部坐标架下的分量
切丛中的元素
( q , v ) ∈ T M (q,v) \in TM ( q , v ) ∈ T M
将 ( q , v ) (q,v) ( q , v ) 的全体与 ⨆ p ∈ M T p M \bigsqcup_{p \in M} T_{p} M ⨆ p ∈ M T p M 自然认同
余切丛中的元素
( q , p ) ∈ T ∗ M (q,p) \in T^* M ( q , p ) ∈ T ∗ M
将 ( q , p ) (q,p) ( q , p ) 的全体与 ⨆ p ∈ M T p ∗ M \bigsqcup_{p \in M} T_p^* M ⨆ p ∈ M T p ∗ M 自然认同
Lagrange量
L = L ( q , v ) L = L(q,v) L = L ( q , v )
Lagrange量是 T M TM T M 到 R \mathbb{R} R 的映射
Hamilton量
H = H ( q , p ) H = H(q,p) H = H ( q , p )
Hamilton量是 T ∗ M T^* M T ∗ M 到 R \mathbb{R} R 的映射
拉格朗日力学系统
位形空间 一个力学系统可能处于的所有状态的空间称为位形空间。力学系统的状态可由一系列广义坐标 q i q^i q i 描述 ( i = 1 , ⋯ , n ) (i=1,\cdots,n) ( i = 1 , ⋯ , n ) ,其可视作流形 M M M 的局域坐标,因此位形空间也成为位形流形 (configuration manifold) 。位形空间 M M M 的维数 n = dim M n= \dim M n = dim M 称为力学系统的自由度。
拉格朗日力学系统 设流形 M M M 为位形空间,T M TM T M 为切丛,L : T M → R L:TM \to \mathbb{R} L : T M → R 是切丛上的函数。若力学系统的任一容许轨迹 γ : R → M \gamma: \mathbb{R} \to M γ : R → M 都使泛函
S [ γ ] = ∫ t 0 t 1 L ( γ , γ ˙ ) d t S[\gamma] = \int_{t_0}^{t_1} L(\gamma,\dot{\gamma}) \,\mathrm{d}t
S [ γ ] = ∫ t 0 t 1 L ( γ , γ ˙ ) d t
取极值,则称 L L L 为拉格朗日量 ,S S S 为作用量 ,γ \gamma γ 为拉格朗日力学系统下的运动轨迹。
/* 注:γ ˙ \dot{\gamma} γ ˙ 为曲线的切矢量,当轨迹确定时,其切矢也自然确定 */
定理:拉格朗日力学系统中的运动轨迹 γ \gamma γ 遵循拉格朗日方程
∂ L ∂ q i − d d t ∂ L ∂ q ˙ i = 0 \frac{\partial L}{\partial q^i}-\frac{\mathrm{d}}{\mathrm{d} t} \frac{\partial L}{\partial \dot{q}^i} = 0
∂ q i ∂ L − d t d ∂ q ˙ i ∂ L = 0
其中 q i q^i q i 为广义坐标/局域坐标,q ˙ i \dot{q}^i q ˙ i 为轨迹切矢 ∂ γ ∂ t \frac{\partial \gamma}{\partial t} ∂ t ∂ γ 在切空间局域坐标架 ∂ ∂ q i \frac{\partial}{\partial q^i} ∂ q i ∂ 上的分量。
需要说明的是,一般切丛 T M TM T M 上任意一点的局域坐标应当用 M M M 上一点的局域坐标 q i q^i q i 和该点切向量的分量 v i v^i v i 共同描述,拉格朗日函数作为 T M TM T M 到 R \mathbb{R} R 的映射也应当写成 L = L ( q 1 , ⋯ , q n , v 1 , ⋯ , v n ) L = L (q^1, \cdots,q^n,v^1,\cdots,v^n) L = L ( q 1 , ⋯ , q n , v 1 , ⋯ , v n ) 。但由于物理中解出的运动轨迹一定是 C ∞ C^\infty C ∞ 的,每一点的切向量都应当是真实光滑运动轨迹的切向量 q ˙ i ∂ ∂ q i \dot{q}^i \frac{\partial}{\partial q^i} q ˙ i ∂ q i ∂ ,而非可随意指定的 v i ∂ ∂ q i v^i \frac{\partial}{\partial q^i} v i ∂ q i ∂ 。因此物理中拉格朗日函数便干脆记作 L = L ( q 1 , ⋯ , q n , q ˙ 1 , ⋯ , q ˙ n ) L = L (q^1, \cdots,q^n,\dot{q}^1,\cdots,\dot{q}^n) L = L ( q 1 , ⋯ , q n , q ˙ 1 , ⋯ , q ˙ n ) ,但需谨记此式来源于轨迹的光滑性质,一般情况下每点切向量并非是自然选定的。
勒让德变换
勒让德变换的引入
理论力学中引入Hamilton量时,所作的勒让德变换通常写作
H ( q , p ) = ∑ p i q ˙ i − L ( q , q ˙ ) H(q,p) = \sum p_i \dot{q}^i - L(q,\dot{q})
H ( q , p ) = ∑ p i q ˙ i − L ( q , q ˙ )
也就是说,在某一个确定的位置 q q q 附近,给定关于 q ˙ \dot{q} q ˙ 的函数 L L L ,能够唯一地确定 p = ∂ L ∂ q ˙ p=\frac{\partial L}{\partial \dot{q}} p = ∂ q ˙ ∂ L ,从而完成变量替换,构成关于 p p p 的函数 H H H 。
纤维导数 给定 T M TM T M 上的函数 L ( q , v ) L(q,v) L ( q , v ) ,定义映射 F L : T M → T ∗ M \mathbb{F}L:TM \to T^*M F L : T M → T ∗ M ,且能将纤维 T q M T_q M T q M 映射到同一点的纤维 T q ∗ M T_q^*M T q ∗ M 之上,规则如下:对于任意 v ∈ T q M v \in T_q M v ∈ T q M ,规定 v v v 在映射 F L \mathbb{F}L F L 下的像 F L ( q , v , ⋅ ) ∈ T q ∗ M \mathbb{F}L(q,v,\cdot) \in T_q^*M F L ( q , v , ⋅ ) ∈ T q ∗ M 为
F L ( q , v , ⋅ ) : T q M → R w ↦ F L ( q , v , w ) = d d t ∣ t = 0 L ( v + t w ) \begin{aligned}
\mathbb{F}L(q,v,\cdot) : T_qM &\to \mathbb{R} \\
w &\mapsto \mathbb{F}L(q,v,w)=\frac{\mathrm{d}}{\mathrm{d}t}\bigg|_{t=0} L(v+tw)
\end{aligned}
F L ( q , v , ⋅ ) : T q M w → R ↦ F L ( q , v , w ) = d t d ∣ ∣ ∣ ∣ ∣ t = 0 L ( v + t w )
其中 w ∈ T q M w \in T_q M w ∈ T q M 为任意切矢量,称 F L ( q , v , w ) \mathbb{F}L(q,v,w) F L ( q , v , w ) 为 L L L 在纤维 T q M T_q M T q M 上 v v v 处沿着方向 w w w 的纤维导数 ,因此
F L : T M → T ∗ M ( q , v ) ↦ F L ( q , v , ⋅ ) = ( q , p ) \begin{aligned}
\mathbb{F}L:T M &\to T^*M \\
(q,v) &\mapsto \mathbb{F}L(q,v,\cdot)=(q,p)
\end{aligned}
F L : T M ( q , v ) → T ∗ M ↦ F L ( q , v , ⋅ ) = ( q , p )
/* 注:纤维导数类似于欧氏空间中的方向导数,取 M = R n M=\mathbb{R}^n M = R n ,则 F L ( q , v , w ) \mathbb{F}L(q,v,w) F L ( q , v , w ) 就是 ∂ L ∂ v ⋅ w \frac{\partial L}{\partial \boldsymbol{v}} \cdot \boldsymbol{w} ∂ v ∂ L ⋅ w */
在局部坐标表示下,纤维导数有表达式
F L ( q , v , ⋅ ) = ( q i , ∂ L ∂ v i d q i ) \mathbb{F}L(q,v,\cdot) = (q^i,\frac{\partial L}{\partial v^i} \mathrm{d}q^i)
F L ( q , v , ⋅ ) = ( q i , ∂ v i ∂ L d q i )
若记 p i = ∂ L ∂ v i p_i=\frac{\partial L}{\partial v^i} p i = ∂ v i ∂ L ,就得到熟知的广义动量。当 L L L 是非退化的拉格朗日函数时,映射 F L \mathbb{F}L F L 可逆,成为切丛到余切丛的局部微分同胚。
在能够定义勒让德变换之前,还有最后一步。由于 ( q , p ) (q,p) ( q , p ) 无法直接作用于 ( q , v ) (q,v) ( q , v ) 之上得到 p i v i p_i v^i p i v i ,还需要引入一个从 T ∗ M T^* M T ∗ M 到 T ∗ ( T M ) T^*(TM) T ∗ ( T M ) 的自然映射 ξ \xi ξ ,然后让 ξ ( q , p ) ∈ T ∗ ( T M ) \xi(q,p) \in T^*(TM) ξ ( q , p ) ∈ T ∗ ( T M ) 作用于 ( q , v ) (q,v) ( q , v ) 上:
ξ : T ∗ M → lift H ∗ ( T ∗ M ) → F L − 1 T ∗ ( T M ) ( q i , p i ) ↦ ( ( q i , p i ) , p i d q i ) ↦ ( ( q i , v i ) , p i d q i ) \begin{gathered}
\xi:T^* M \xrightarrow{\text{lift}} H^*(T^*M) \xrightarrow{\mathbb{F}L^{-1}} T^* (TM) \\
(q^i,p_i) \mapsto ((q^i,p_i),p_i \mathrm{d}q^i) \mapsto ((q^i,v^i),p_i \mathrm{d}q^i)
\end{gathered}
ξ : T ∗ M lift H ∗ ( T ∗ M ) F L − 1 T ∗ ( T M ) ( q i , p i ) ↦ ( ( q i , p i ) , p i d q i ) ↦ ( ( q i , v i ) , p i d q i )
其中 l i f t \mathrm{lift} l i f t 是余切丛上的水平提升,F L − 1 \mathbb{F}L^{-1} F L − 1 是 F L \mathbb{F}L F L 的逆映射。由于水平提升和 F L − 1 \mathbb{F}L^{-1} F L − 1 都是局部同构,因此 ξ \xi ξ 也可逆。
/* 注:忘记水平提升的回去复习纤维丛理论! */
勒让德变换的定义
勒让德变换 勒让德变换是从切丛上的函数 L L L 到余切丛上的函数 H H H 的一个映射,设 F L : ( q , v ) ↦ ( q , p ) \mathbb{F}L:(q,v) \mapsto (q,p) F L : ( q , v ) ↦ ( q , p ) 为局部微分同胚,则
H ( q , p ) = [ ξ ( q , p ) − L ] ∘ F L − 1 H(q,p) = \left[ \xi(q,p)-L \right] \circ\mathbb{F}L^{-1}
H ( q , p ) = [ ξ ( q , p ) − L ] ∘ F L − 1
函数 H ( q , p ) H(q,p) H ( q , p ) 称为Hamilton量 ,流形 T ∗ M T^*M T ∗ M 称为相空间 或相流形 ,力学系统称为Hamilton力学系统 。
哈密顿力学系统
辛空间
待更
相空间的辛结构
辛流形 设 N N N 是一个偶数维流形,dim N = 2 n \dim N=2n dim N = 2 n 。N N N 上一个辛结构/辛形式 是 N N N 上一个闭的非简并的 2-形式场 ω \omega ω ,即
d ω = 0 \mathrm{d}\omega=0 d ω = 0
∀ η ∈ Γ ( T N ) , ω ( ξ , η ) = 0 ⟺ ξ = 0 \forall \eta \in \varGamma(TN),\omega(\xi,\eta)=0 \iff \xi=0 ∀ η ∈ Γ ( T N ) , ω ( ξ , η ) = 0 ⟺ ξ = 0
赋予了辛结构的流形被称为辛流形。
正则辛形式 在 2 n 2n 2 n 维的一个辛流形 N N N 的每个点的一个邻域中,存在局域坐标 ( q i , p i ) (q^i,p_i) ( q i , p i ) ,其中 i = 1 , ⋯ , n i=1,\cdots,n i = 1 , ⋯ , n ,使得辛结构 ω \omega ω 可以被表述成正则形式
ω = d p i ∧ d q i \omega = \mathrm{d} p_i \wedge \mathrm{d} q^i
ω = d p i ∧ d q i