外观
Hamilton 力学与 Poisson Bracket
第一次接触 Hamilton 力学时,常见的感受是:公式看起来很整齐,但不知道它到底比牛顿方程多说了什么。再往后看到 Poisson bracket,又容易觉得这只是一个“定义出来的代数运算”,和动力系统本身没有直接关系。其实这两件事是连在一起的。
Hamilton 形式做的事情可以概括成一句话:
不再直接写“受力等于质量乘加速度”,而是把动力系统写成一种由能量驱动的几何演化。
而 Poisson Bracket 则是在这套结构里回答另一个问题:
给定一个物理量,它随时间怎么变?
这篇文章只讲最基本的有限维情形。目标有三个:
- 解释 Hamilton 方程到底在描述什么;
- 说明 Poisson Bracket 为什么会自然出现;
- 顺手说明它为什么和守恒律、PDE 结构有关。
从相空间开始
在 Hamilton 力学里,一个系统的状态通常不用位置和速度写,而是用位置和动量写成
(q,p).
这里
- q 表示位置;
- p 表示动量。
所有可能的 (q,p) 组成的空间,叫做相空间。对一个 n 自由度系统来说,q∈Rn,p∈Rn,所以相空间是 2n 维的。
这一步看起来只是换变量,但含义其实很重要。牛顿方程通常写成二阶方程,而 Hamilton 形式把它改写成一阶系统。这样做的好处是:整个动力系统可以被统一看成相空间里的流动。
Hamiltonian
Hamiltonian 通常记作 H(q,p)。最常见的情形下,它就是系统总能量:
H(q,p)=T(p)+V(q),
其中 T 是动能,V 是势能。
例如一个经典粒子系统可以写成
H(q,p)=2m∣p∣2+V(q).
因为
p=mq˙,
所以
2m∣p∣2
正好就是动能。
这里最关键的点不是“它等于能量”这句话本身,而是:
一旦写出了 H(q,p),系统的演化规律就能直接从它读出来。
Hamilton 方程
Hamilton 力学把系统的时间演化写成
q˙=∂p∂H,p˙=−∂q∂H.
这就是 Hamilton 方程。
它的结构很值得注意:
- 位置的变化率由 H 对动量的导数决定;
- 动量的变化率由 H 对位置的导数决定,并且带一个负号。
这说明动力学不再是“分别给 q 和 p 写方程”,而是由一个统一的标量函数 H 同时生成。
用最经典的例子看一遍
若
H(q,p)=2mp2+V(q),
那么 Hamilton 方程给出
q˙=∂p∂H=mp,
以及
p˙=−∂q∂H=−∇V(q).
第一式就是
v=mp,
第二式则是
F=−∇V.
再结合 p=mq˙,就能恢复熟悉的牛顿方程。所以 Hamilton 形式并不是另一套物理学,而是对同一个动力系统的另一种写法。它真正的价值在于:这种写法把系统背后的几何结构暴露了出来。
Poisson Bracket 为什么会出现
到这里自然会冒出一个问题。
如果系统状态由 Hamilton 方程控制,那么一个一般的物理量
F(q,p)
会怎样随时间演化?
对 F(q,p) 直接用链式法则:
dtdF=i=1∑n(∂qi∂Fq˙i+∂pi∂Fp˙i).
再把 Hamilton 方程代进去:
dtdF=i=1∑n(∂qi∂F∂pi∂H−∂pi∂F∂qi∂H).
这个表达式太常见了,于是人们把它单独记成一个运算:
{F,G}=i=1∑n(∂qi∂F∂pi∂G−∂pi∂F∂qi∂G)
这就是标准的 Poisson Bracket。
因此 Hamilton 系统中任意物理量的演化都可以统一写成
dtdF={F,H}
这句话非常重要。它说明:
Hamiltonian H 不是只负责“能量解释”,它还是生成时间演化的那个函数。
Hamilton 方程可以从 Poisson Bracket 反推回来
上面的公式并不是抽象装饰,它真的把原来的动力学统一了。
取
F=qi,
则
{qi,H}=∂pi∂H,
所以
q˙i={qi,H}=∂pi∂H.
再取
F=pi,
则
{pi,H}=−∂qi∂H,
所以
p˙i={pi,H}=−∂qi∂H.
也就是说,Hamilton 方程本身就是
F˙={F,H}
在坐标函数 qi,pi 上的特例。
Poisson Bracket 的几个基本性质
Poisson Bracket 之所以重要,不只是因为它能把链式法则压缩成一个符号,还因为它满足一组和动力系统密切相关的结构性质。
双线性
Poisson Bracket 对两个输入都线性:
{aF+bG,H}=a{F,H}+b{G,H},
以及
{F,aG+bK}=a{F,G}+b{F,K}.
这保证了它确实像一个“可计算的代数运算”。
反对称性
{F,G}=−{G,F}
这条性质看起来简单,但它立刻带出一个核心结论。取 F=G=H,就有
{H,H}=0,
于是
H˙={H,H}=0.
也就是说,在没有显含时间的 Hamilton 系统里,Hamiltonian 自动守恒。能量守恒不是额外加上的,而是这套结构内生的结果。
Leibniz 规则
Poisson Bracket 对乘积满足类似导数的规则:
{FG,H}=F{G,H}+G{F,H}.
这说明它和“时间导数”之间的行为非常一致。
Jacobi 恒等式
{F,{G,H}}+{G,{H,F}}+{H,{F,G}}=0
这是 Poisson 结构最深的一条性质。它保证这种括号不是随便写出来的双线性算子,而是和几何上的相空间结构真正兼容。
从矩阵形式看,结构会更清楚
把相空间坐标记成
z=(q,p),
再引入标准辛矩阵
J=(0−II0),
那么 Poisson Bracket 可以写成
{F,G}=∇FTJ∇G
相应地,Hamilton 方程就变成
z˙=J∇H.
写成这个形式以后,Hamilton 系统的几何味道就更明显了:演化方向由梯度 ∇H 决定,但真正推动系统流动的是那个反对称矩阵 J。这就是标准辛几何结构在坐标里的体现。
为什么说 Casimir 不依赖具体 Hamiltonian
有些量的守恒并不是因为某个特定 Hamiltonian 的形式,而是因为 Poisson 结构本身。
如果一个量 C 满足
{C,G}=0∀G
那么 C 称为一个 Casimir。
它的含义是:无论系统选哪个 Hamiltonian,只要 Poisson 结构不变,这个量都会守恒。于是 Casimir 反映的就不是“某个能量函数的特殊性质”,而是几何结构本身的退化方向或不变量。
这类量在流体和等离子体问题里尤其常见。
为什么 PDE 里也会出现 Hamilton 结构
前面讲的是有限维系统,也就是状态由有限个坐标 (q,p) 描述的情形。到了 PDE 里,状态变量不再是有限维向量,而是函数,例如
u(x),h(x).
这时普通偏导数会换成泛函导数:
∂q∂⟶δuδ.
但核心思路没有变:仍然是
- 先写出一个 Hamiltonian;
- 再通过某种 Poisson 结构生成演化。
换句话说,PDE 里的 Hamilton 结构不是“完全不同的理论”,而是把有限维的思路推广到无穷维函数空间。
以浅水方程为例
浅水方程里,Hamiltonian 通常对应总能量:
H=∫(21h∣u∣2+21gh2)dA
其中
- 第一项是动能;
- 第二项是重力势能。
如果继续往 Hamilton 结构的方向写下去,就需要给出一个适合浅水系统的 Poisson Bracket。此时括号不再是有限维里那种对 (q,p) 的简单偏导,而会涉及速度、厚度、位涡等变量之间的耦合。
在离散浅水系统中,常见的形式会写成类似
{F,H}=⟨δuhδF,−qh(δuhδH)⊥⟩+⟨∇⋅δuhδF,δDhδH⟩−⟨δDhδF,∇⋅δuhδH⟩.
这类表达式看起来复杂,但本质仍然是同一件事:
Hamiltonian 负责告诉你系统储存了什么能量,Poisson 结构负责告诉你这些能量怎样在各个变量之间交换,并生成整个动力学。
也正因为如此,Hamilton 结构常常和下面这些性质绑在一起出现:
- 能量守恒;
- 位涡结构;
- 涡量或 enstrophy 之类的不变量;
- 数值格式中的结构保持性。
Almost-Poisson bracket 是什么
在很多 PDE 或离散系统里,人们会先写出一个“看起来像 Poisson Bracket”的东西,但还没有证明它满足 Jacobi 恒等式。这时更准确的叫法是 almost-Poisson bracket。
也就是说,它可能已经满足:
- 双线性;
- 反对称性;
- Leibniz 规则。
但如果 Jacobi 恒等式还没证出来,就不能直接把它叫成真正的 Poisson Bracket。这个区分在结构保持离散里很重要,因为很多守恒性质和几何解释最终都依赖 Jacobi 恒等式是否成立。
小结
把全文压缩成最核心的几句话,大概就是:
- Hamiltonian 描述系统的能量;
- Hamilton 方程把动力学写成相空间中的一阶流;
- Poisson Bracket 把“任意物理量如何演化”统一写成 F˙={F,H};
- 反对称性带来能量守恒,Jacobi 恒等式保证这套结构真正是几何性的;
- 到 PDE 里,变量从有限维坐标变成函数,但“能量 + Poisson 结构生成演化”的主线不变。
如果只记一句话,可以记成:
Hamiltonian 决定系统储存什么能量,Poisson Bracket 决定系统怎样演化。
版权所有
版权归属:Guisong Wu