欢迎光临散文网 会员登陆 & 注册

哈密顿力学概述||理论力学

2021-09-22 16:36 作者:湮灭的末影狐  | 我要投稿

// 笔者最近在物理学基础理论的学习上好像遇到了一个瓶颈。

// 越往后学,好像越容易接触到一些高度抽象的概念,光有高数线代数理方法搞不定,例如酉空间辛空间,群论,奇奇怪怪的流形...

// 其实,一个普普通通不搞纯理论的物理人,对这些东西好像也不需要掌握太多...但是我就是想搞清楚一些。

// 就比如这个理论力学,往简单了说会套公式就能应付了期末考,往难了说什么辛几何,相空间的余切丛,三维特殊正交群之类奇奇怪怪的概念都来了...

//总之,这些笔记笔者不想简单应付,而想加入自己的理解。最近各种专栏笔记都会更得比较慢,笔者需要等待一个突破...

0 Introduction

理论力学主要内容分为刚体力学、拉格朗日力学和哈密顿力学。

其中,拉格朗日力学最早引入了广义坐标的概念,用纯代数的方法处理了约束的影响,并给出拉格朗日方程:

%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20q%7D-%5Cfrac%7B%5Cmathrm%7Bd%7D%7D%7B%5Cmathrm%7Bd%7Dt%7D(%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20%5Cdot%20q%7D)%3D0

其中,L%3DT-V(动能与势能之差)是系统的拉格朗日量。

拉格朗日力学的优势在于它能完美解决约束带来的影响。笔者曾经借助拉格朗日方程完成混沌摆的建模和模拟,这也是笔者在计算机基础课程中的大作业内容。

而哈密顿力学中则将广义动量和广义坐标一样作为独立变量,构建出哈密顿正则方程,给出动量、坐标、哈密顿量的关系。哈密顿力学在本质上和拉格朗日力学和牛顿力学都是等价的,但其数学形式在其他领域得到广泛的应用,其意义是超越经典力学本身的。(当你看到在哈密顿力学中坐标与动量在数学形式上的对称关系,或许会想到量子力学...)

1 哈密顿正则方程

已知某个经典力学系统有拉格朗日量L,我们定义广义动量为:

p_i%3D%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20%5Cdot%20q_i%7D

然后,代入拉格朗日方程,不难得到

%5Cdot%20p_i%3D%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20q_i%7D

接下来,我们定义系统的哈密顿量:

H(p%2Cq%2Ct)%20%3D%20-L%20%2B%20%5Csum_%7Bi%7D%20p_i%20%5Cdot%20q_i

注意,虽然上式右边显含 %5Cdot%20q,但 H 的最终表达式是 p,q 的函数,也就是说最终 %5Cdot%20q 应当用 p,q 表示。接下来研究哈密顿量的微分:

%5Cmathrm%7Bd%7DH%20%3D%20-%5Cmathrm%7Bd%7DL%20%2B%20%5Csum_i%20(p_i%20%5Cmathrm%7Bd%7D%20%5Cdot%20q_i%2B%20%5Cdot%20q_i%20%5Cmathrm%7Bd%7D%20p_i)

而拉格朗日量的微分

%5Cbegin%7Baligned%7D%0A%20%20%20%20%5Cmathrm%7Bd%7DL%20%26%3D%20%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20t%7D%5Cmathrm%7Bd%7Dt%20%2B%20%5Csum_i%20(%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20q_i%7D%5Cmathrm%7Bd%7Dq_i%2B%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20%5Cdot%20q_i%7D%5Cmathrm%7Bd%7D%5Cdot%20q_i)%5C%5C%0A%20%20%20%20%26%3D%20%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20t%7D%5Cmathrm%7Bd%7Dt%20%2B%20%5Csum_i%20(%5Cdot%20p_i%5Cmathrm%7Bd%7Dq_i%2Bp_i%5Cmathrm%7Bd%7D%5Cdot%20q_i)%0A%5Cend%7Baligned%7D

联立两式,注意到 p_i%20%5Cmathrm%20d%20%5Cdot%20q_i项被消去,得

%5Cmathrm%20%20d%20H%3D%5Csum_i(%5Cdot%7Bq%7D_%7Bi%7D%5Cmathrm%20d%20p_%7Bi%7D-%5Cdot%7Bp%7D_%7Bi%7D%20%5Cmathrm%20d%20q_%7Bi%7D)-%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20t%7D%20%5Cmathrm%20d%20t

上式正是以 p,q,t 为独立变量得到的H的全微分,由此我们得到哈密顿正则方程:

%5Cdot%7Bq%7D_%7Bi%7D%3D%5Cfrac%7B%5Cpartial%20H%7D%7B%5Cpartial%20p_%7Bi%7D%7D%2C%20%5C%3B%5C%3B%20%5Cdot%7Bp%7D_%7Bi%7D%3D-%5Cfrac%7B%5Cpartial%20H%7D%7B%5Cpartial%20q_%7Bi%7D%7D

此外,如果考虑 H 对时间的微分,

%5Cfrac%7B%5Cmathrm%20d%20H%7D%7B%5Cmathrm%20d%20t%7D%3D%5Cfrac%7B%5Cpartial%20H%7D%7B%5Cpartial%20t%7D%2B%5Csum_i(%5Cfrac%7B%5Cpartial%20H%7D%7B%5Cpartial%20q_%7Bi%7D%7D%20%5Cdot%7Bq%7D_%7Bi%7D%2B%5Cfrac%7B%5Cpartial%20H%7D%7B%5Cpartial%20p_%7Bi%7D%7D%20%5Cdot%7Bp%7D_%7Bi%7D)

把前面的正则方程代入,恰好消去后两项,可得

%5Cfrac%7B%5Cmathrm%20d%20H%7D%7B%5Cmathrm%20d%20t%7D%3D%5Cfrac%7B%5Cpartial%20H%7D%7B%5Cpartial%20t%7D%3D-%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20t%7D

由此看出,如果哈密顿量不显含时间,则它是守恒量。事实上,哈密顿量就是系统的广义能量,它不显含时间意味着系统受到稳定约束,约束反力不做功。

哈密顿正则方程的数学结构很有意思。可以看到,仅从数学上看坐标与动量的地位是高度对称的,这也意味着从数学的角度我们并不能区分坐标与动量。

2 哈密顿原理

哈密顿原理,也被称作最小作用量原理。费曼教授在他的三部物理学讲义中鲜有提及分析力学的内容,但唯独用了一章专门介绍了这个最小作用量原理——这显然是一个让过去年轻的他感到相当印象深刻的理论。

作用量的定义是系统的拉格朗日量对某段时间的积分:

S%3D%5Cint_%7Bt_1%7D%5E%7Bt_2%7D%20L(q%2C%5Cdot%20q%2Ct)%5Cmathrm%20d%20t

哈密顿原理表示其变分为0:

%5Cdelta%20S%20%3D%200

通过展开 S 的变分可以推回拉格朗日方程,笔者在本文集前面的笔记提过了。

此外,作用量存在另一种包含广义动量的形式。前面的文章里面我们以及操作过一次变分了,

%5Cdelta%20S%20%3D%20%5Cdelta%5Cint%20L(q%2C%5Cdot%20q%20%2C%20t)%5Cmathrm%7Bd%7Dt%20%3D%20%5Csum_i%20%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20%5Cdot%20q_i%7D%5Cdelta%20q_i%7C_%7Bt_1%7D%5E%7Bt_2%7D%2B%5Cint_%7Bt_1%7D%5E%7Bt_2%7D%20(%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20q_i%7D-%5Cfrac%7B%5Cmathrm%7Bd%7D%7D%7B%5Cmathrm%7Bd%7Dt%7D%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20%5Cdot%20q_i%7D)%5Cdelta%20q_i%5Cmathrm%7Bd%7Dt

当时,我们令 %5Cdelta%20q(t_1)%3D%5Cdelta%20q(t_2)%3D0 而推出第二项为0. 而这一次我们选择把 t2 看作任意时间,看第一项易得

%5Cfrac%7B%5Cpartial%20S%7D%7B%5Cpartial%20q_i%7D%20%3D%20%5Cfrac%7B%5Cpartial%20L%7D%7B%5Cpartial%20%5Cdot%20q_i%7D%20%3D%20p_i

又,根据定义,%5Cfrac%7B%5Cmathrm%7Bd%7D%20S%7D%7B%5Cmathrm%7Bd%7D%20t%7D%20%3DL

由此联系又可以推出哈密顿-雅可比方程

%5Cfrac%7B%5Cpartial%20S%7D%7B%5Cpartial%20t%7D%20%2BH%3D0

此外,联系哈密顿量和作用量的定义,可以得到另一种形式的作用量

S%3D%5Cint%20(-H%5Cmathrm%20d%20t%2B%5Csum_i%20p_i%20%5Cmathrm%7Bd%7D%20%5Cdot%20q_i)

3 泊松括号

泊松括号似乎是某种数学结构。设存在广义坐标、动量、时间的函数f(p%2Cq%2Ct),它对时间的微分:

%5Cfrac%7B%5Cmathrm%20d%20f%7D%7B%5Cmathrm%20d%20t%7D%3D%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20t%7D%2B%5Csum_i(%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20q_%7Bi%7D%7D%20%5Cdot%7Bq%7D_%7Bi%7D%2B%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20p_%7Bi%7D%7D%20%5Cdot%7Bp%7D_%7Bi%7D)

把正则方程代入,并把求和里面的东西定义为泊松括号:

%5Bf%2Cg%5D%3D%5Csum_i(%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20q_i%7D%5Cfrac%7B%5Cpartial%20g%7D%7B%5Cpartial%20p_i%7D-%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20p_i%7D%20%5Cfrac%7B%5Cpartial%20g%7D%7B%5Cpartial%20q_i%7D%20%20%20)

则可以得到

%5Cfrac%7B%5Cmathrm%20d%20f%7D%7B%5Cmathrm%20d%20t%7D%3D%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20t%7D%2B%5Bf%2CH%5D

据说,量子力学里面常常遇到的对易子,起源和这个泊松括号有些关系。

泊松定理是说,如果已知系统的两个积分

%5Cphi(p%2Cq%2Ct)%3DC_1%2C%5C%3B%5C%3B%5Cpsi(p%2Cq%2Ct)%3DC_2

%5B%5Cphi%2C%5Cpsi%5D%3DC_3

也是系统的一个积分。(虽然,到目前为止笔者未能看到这有什么用)

4 正则变换

正则变换本质上是一种坐标变换。不过变换的不再是空间坐标,而是广义坐标和广义动量。

我们希望将本来的广义坐标与动量 q,p 变换为一组新的坐标 Q,P:

Q_i%20%3D%20Q_i(p%2Cq%2Ct)%2C%5C%3BP_i%20%3D%20P_i(p%2Cq%2Ct)

同时希望这个变换下正则方程保持不变:

%5Cdot%20Q_i%20%3D%20%5Cfrac%7B%5Cpartial%20H'%7D%7B%5Cpartial%20P_i%7D%2C%5C%3B%20%5Cdot%20P_i%20%3D%20-%5Cfrac%7B%5Cpartial%20H'%7D%7B%5Cpartial%20Q_i%7D

前文写过包含广义动量的作用量表达式,所以系统满足

%5Cdelta%20S%3D%5Cdelta%20%5Cint%5Cleft(p_%7Bi%7D%5Cmathrm%20%20d%20q_%7Bi%7D-H%20%5Cmathrm%20d%20t%5Cright)%3D0

同理对于变换后,

%5Cdelta%20%5Cint%5Cleft(P_%7Bi%7D%5Cmathrm%20%20d%20Q_%7Bi%7D-H'%20%5Cmathrm%20d%20t%5Cright)%3D0

两个变分相等,意味着它们的被积函数应当仅相差一个系统状态决定的任意函数的全微分dF. (因为这一项积分得到的是 F(末态)-F(初态),它的变分为0.)

p_%7Bi%7D%5Cmathrm%20%20d%20q_%7Bi%7D-H%20%5Cmathrm%20d%20t%3DP_%7Bi%7D%20%5Cmathrm%20d%20Q_%7Bi%7D-H%5E%7B%5Cprime%7D%5Cmathrm%20%20d%20t%2B%5Cmathrm%20d%20F

这里 F 称为生成函数,或母函数。移项不难得到 F 的全微分,从而得到:

p_%7Bi%7D%3D%5Cfrac%7B%5Cpartial%20F%7D%7B%5Cpartial%20q_%7Bi%7D%7D%2C%20%5Cquad%20P_%7Bi%7D%3D-%5Cfrac%7B%5Cpartial%20F%7D%7B%5Cpartial%20Q_%7Bi%7D%7D%2C%20H%5E%7B%5Cprime%7D%3DH%2B%5Cfrac%7B%5Cpartial%20F%7D%7B%5Cpartial%20t%7D

这里如果母函数恰好是 q,Q 的函数,上面就已经是正则变换的表达式了。

母函数还可能是:F(q, P, t), F(Q, p, t), F(p, P, t) (一定是一个新变量和一个旧变量的函数)

以 F(q,P,t) 为例,只需构造

%5CPhi(q%2CP%2Ct)%3DF%2B%5Csum_i%20P_iQ_i

接下来类似分部积分的原理,可以得到

p_%7Bi%7D%3D%5Cfrac%7B%5Cpartial%20%5CPhi%7D%7B%5Cpartial%20q_%7Bi%7D%7D%20.%20%5Cquad%20Q_%7Bi%7D%3D%5Cfrac%7B%5Cpartial%20%5CPhi%7D%7B%5Cpartial%20p_%7Bi%7D%7D%20.%20%5Cquad%20H%5E%7B%5Cprime%7D%3DH%2B%5Cfrac%7B%5Cpartial%20%5CPhi%7D%7B%5Cpartial%20t%7D%20.

正则变换后的坐标与动量已经变成抽象的概念,仅仅保留了原先的数学结构。

//这篇笔记暂时到此为止。理论力学到后面越发高深,当前有待更深入的理解。

//后面可能会讨论的内容:如何寻找母函数;正则变换的应用;哈密顿力学中的辛结构等

哈密顿力学概述||理论力学的评论 (共 条)

分享到微博请遵守国家法律