物理 理论力学核心 连续系统与场的拉格朗日-哈密顿方法 连续系统与场的拉格朗日-哈密顿方法
前面各章所处理的系统,都有有限个自由度——无论是单个质点、刚体,还是由有限个粒子组成的力学系统。然而,弦的横向振动、声波在空气中的传播、电磁场在空间中的演化,这些物理现象本质上涉及无穷多个自由度。把拉格朗日力学和哈密顿力学推广到这类系统,既是经典力学体系的自然延伸,也是理解场论和现代物理的入口。从弦的振动出发,可以清晰地看到离散体系如何平滑过渡到连续场,以及对称性如何通过诺特定理与守恒律产生深层联系。
从离散系统到连续系统
以一根紧绷的弦为例。先把它离散化:将弦分成 N N N 段,每段长度 Δ x \Delta x Δ x ,每段端点处有一个质量为 m = μ Δ x m = \mu\Delta x m = μ Δ x 的质点(μ \mu μ 为线质量密度,单位 k g / m \mathrm{kg/m} kg/m ),相邻质点之间以劲度系数 κ \kappa κ 的弹簧连接。第 i i i 个质点的横向位移记为 y i ( t ) y_i(t) y i ( t ) ,其运动方程为:
μ Δ x y ¨ i = κ ( y i + 1 − y i ) − κ ( y i − y i − 1 ) = κ ( y i + 1 − 2 y i + y i − 1 ) \mu\Delta x\,\ddot{y}_i = \kappa\!\left(y_{i+1} - y_i\right) - \kappa\!\left(y_i - y_{i-1}\right) = \kappa\!\left(y_{i+1} - 2y_i + y_{i-1}\right) μ Δ x y ¨ i = κ ( y i + 1 − y i ) − κ ( y i − y i − 1 ) = κ ( y i + 1 − 2 y i + y i − 1 )
令 Δ x → 0 \Delta x \to 0 Δ x → 0 ,同时保持弦的张力 T = κ Δ x T = \kappa\Delta x T = κ Δ x 不变,右边的有限差分趋近于二阶偏导数:
y i + 1 − 2 y i + y i − 1 ( Δ x ) 2 ⟶ ∂ 2 y ∂ x 2 \frac{y_{i+1} - 2y_i + y_{i-1}}{(\Delta x)^2} \;\longrightarrow\; \frac{\partial^2 y}{\partial x^2} ( Δ x ) 2 y i + 1 − 2 y i + y i − 1 ⟶ ∂ x 2 ∂ 2 y
于是离散的运动方程平滑过渡为连续弦的波动方程 :
μ ∂ 2 y ∂ t 2 = T ∂ 2 y ∂ x 2 \mu\frac{\partial^2 y}{\partial t^2} = T\frac{\partial^2 y}{\partial x^2} μ ∂ t 2 ∂ 2 y = T ∂ x 2 ∂ 2 y
波速 v = T / μ v = \sqrt{T/\mu} v = T / μ (单位:m / s \mathrm{m/s} m/s ),这就是经典的一维波动方程。
与此同时,离散系统的拉格朗日量也经历同样的过渡。离散系统的拉格朗日量:
L = ∑ i [ 1 2 μ Δ x y ˙ i 2 − 1 2 κ ( y i + 1 − y i ) 2 ] L = \sum_i \left[\frac{1}{2}\mu\Delta x\,\dot{y}_i^2 - \frac{1}{2}\kappa\!\left(y_{i+1}-y_i\right)^2\right] L = i ∑ [ 2 1 μ Δ x y ˙ i 2 − 2 1 κ ( y i + 1 − y i ) 2 ]
连续极限下,对 x x x 的求和变成积分,得到:
L = ∫ L ( y , y ˙ , y ′ ) d x , L = 1 2 μ y ˙ 2 − 1 2 T ( ∂ y ∂ x ) 2 L = \int \mathcal{L}\!\left(y, \dot{y}, y'\right)dx, \quad \mathcal{L} = \frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T\!\left(\frac{\partial y}{\partial x}\right)^{\!2} L = ∫ L ( y , y ˙ , y ′ ) d x , L = 2 1 μ y ˙ 2 − 2 1 T ( ∂ x ∂ y ) 2
函数 L \mathcal{L} L 称为拉格朗日密度 (单位:J / m \mathrm{J/m} J/m ),它是描述连续系统动力学的基本量。
拉格朗日密度与场方程
对于一个广义场 ϕ ( x 1 , x 2 , x 3 , t ) = ϕ ( x μ ) \phi(x^1, x^2, x^3, t) = \phi(x^\mu) ϕ ( x 1 , x 2 , x 3 , t ) = ϕ ( x μ ) (μ = 0 , 1 , 2 , 3 \mu = 0,1,2,3 μ = 0 , 1 , 2 , 3 ),其中 x 0 = t x^0 = t x 0 = t 为时间坐标,系统的拉格朗日量写为:
L = ∫ L ( ϕ , ∂ μ ϕ ) d 3 x L = \int \mathcal{L}\!\left(\phi,\;\partial_\mu\phi\right) d^3x L = ∫ L ( ϕ , ∂ μ ϕ ) d 3 x
对应的作用量为:
S = ∫ L d t = ∫ L ( ϕ , ∂ μ ϕ ) d 4 x S = \int L\,dt = \int \mathcal{L}\!\left(\phi,\;\partial_\mu\phi\right) d^4x S = ∫ L d t = ∫ L ( ϕ , ∂ μ ϕ ) d 4 x
对作用量取极值 δ S = 0 \delta S = 0 δ S = 0 ,经过与离散情形完全平行的变分计算,得到欧拉-拉格朗日场方程 :
∂ L ∂ ϕ − ∂ μ ∂ L ∂ ( ∂ μ ϕ ) = 0 \frac{\partial \mathcal{L}}{\partial \phi} - \partial_\mu\frac{\partial \mathcal{L}}{\partial(\partial_\mu\phi)} = 0 ∂ ϕ ∂ L − ∂ μ ∂ ( ∂ μ ϕ ) ∂ L = 0
这是连续场的运动方程,地位与有限自由度系统中拉格朗日方程完全相同。
弦的拉格朗日密度为 L = 1 2 μ y ˙ 2 − 1 2 T ( y ′ ) 2 \mathcal{L} = \frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T(y')^2 L = 2 1 μ y ˙ 2 − 2 1 T ( y ′ ) 2 ,其中 y ˙ = ∂ y / ∂ t \dot{y} = \partial y/\partial t y ˙ = ∂ y / ∂ t ,y ′ = ∂ y / ∂ x y' = \partial y/\partial x y ′ = ∂ y / ∂ x 。
计算各偏导数:∂ L / ∂ y = 0 \partial\mathcal{L}/\partial y = 0 ∂ L / ∂ y = 0 ,∂ L / ∂ y ˙ = μ y ˙ \partial\mathcal{L}/\partial\dot{y} = \mu\dot{y} ∂ L / ∂ y ˙ = μ y ˙ ,∂ L / ∂ y ′ = − T y ′ \partial\mathcal{L}/\partial y' = -Ty' ∂ L / ∂ y ′ = − T y ′ 。代入欧拉-拉格朗日场方程:
0 − ∂ ∂ t ( μ y ˙ ) − ∂ ∂ x ( − T y ′ ) = 0 ⟹ μ y ¨ = T y ′ ′ 0 - \frac{\partial}{\partial t}(\mu\dot{y}) - \frac{\partial}{\partial x}(-Ty') = 0 \implies \mu\ddot{y} = Ty'' 0 − ∂ t ∂ ( μ y ˙ ) − ∂ x ∂ ( − T y ′ ) = 0 ⟹ μ y ¨ = T y ′′
正好就是前面得到的波动方程,验证了场方程的正确性。
在相对论场论中,质量为 m m m 的标量场 ϕ ( x μ ) \phi(x^\mu) ϕ ( x μ ) 的拉格朗日密度为:
L = 1 2 η μ ν ∂ μ ϕ ∂ ν ϕ − 1 2 ( m c ℏ ) 2 ϕ 2 \mathcal{L} = \frac{1}{2}\eta^{\mu\nu}\partial_\mu\phi\,\partial_\nu\phi - \frac{1}{2}\!\left(\frac{mc}{\hbar}\right)^{\!2}\phi^2 L = 2 1 η μν ∂ μ ϕ ∂ ν ϕ − 2 1 ( ℏ m c ) 2 ϕ 2
其中 η μ ν = d i a g ( + 1 , − 1 , − 1 , − 1 ) \eta^{\mu\nu} = \mathrm{diag}(+1,-1,-1,-1) η μν = diag ( + 1 , − 1 , − 1 , − 1 ) 为闵可夫斯基度规。代入场方程后得到:
□ ϕ + ( m c ℏ ) 2 ϕ = 0 , □ ≡ ∂ 2 ∂ t 2 − c 2 ∇ 2 \Box\phi + \left(\frac{mc}{\hbar}\right)^{\!2}\phi = 0, \quad \Box \equiv \frac{\partial^2}{\partial t^2} - c^2\nabla^2 □ ϕ + ( ℏ m c ) 2 ϕ = 0 , □ ≡ ∂ t 2 ∂ 2 − c 2 ∇ 2
这就是Klein-Gordon方程,描述相对论性自旋为零的粒子的传播,m = 0 m=0 m = 0 时退化为标准的波动方程。
欧拉-拉格朗日场方程 ∂ L / ∂ ϕ − ∂ μ ( ∂ L / ∂ ( ∂ μ ϕ ) ) = 0 \partial\mathcal{L}/\partial\phi - \partial_\mu(\partial\mathcal{L}/\partial(\partial_\mu\phi)) = 0 ∂ L / ∂ ϕ − ∂ μ ( ∂ L / ∂ ( ∂ μ ϕ )) = 0 与离散系统的拉格朗日方程在形式上完全对应:只需将广义坐标 q i q_i q i 替换为场 ϕ \phi ϕ ,对时间的导数替换为对全部时空坐标的偏导数,求和替换为对空间的积分。
连续系统的哈密顿形式
类比离散系统中从拉格朗日量到哈密顿量的过渡,连续系统中先定义正则动量密度 π \pi π ,其定义为:
π ( x , t ) = ∂ L ∂ ϕ ˙ \pi(x,t) = \frac{\partial \mathcal{L}}{\partial \dot{\phi}} π ( x , t ) = ∂ ϕ ˙ ∂ L
再通过对空间的勒让德变换得到哈密顿密度 (单位:J / m 3 \mathrm{J/m^3} J/ m 3 ):
H = π ϕ ˙ − L \mathcal{H} = \pi\dot{\phi} - \mathcal{L} H = π ϕ ˙ − L
系统总哈密顿量(即总能量)为:
H = ∫ H ( ϕ , π , ∇ ϕ ) d 3 x H = \int \mathcal{H}\!\left(\phi, \pi, \nabla\phi\right) d^3x H = ∫ H ( ϕ , π , ∇ ϕ ) d 3 x
连续系统的哈密顿方程为:
ϕ ˙ = δ H δ π , π ˙ = − δ H δ ϕ \dot{\phi} = \frac{\delta H}{\delta \pi}, \quad \dot{\pi} = -\frac{\delta H}{\delta \phi} ϕ ˙ = δ π δH , π ˙ = − δ ϕ δH
其中 δ H / δ ϕ \delta H/\delta\phi δH / δ ϕ 表示泛函导数(functional derivative),其定义与离散情形中的偏导数平行,只是需要同时考虑 H \mathcal{H} H 对 ϕ \phi ϕ 的直接依赖和通过 ∇ ϕ \nabla\phi ∇ ϕ 的间接依赖。
弦的正则动量密度 π = ∂ L / ∂ y ˙ = μ y ˙ \pi = \partial\mathcal{L}/\partial\dot{y} = \mu\dot{y} π = ∂ L / ∂ y ˙ = μ y ˙ ,单位为 k g / s \mathrm{kg/s} kg/s (即单位长度弦所携带的横向动量,[ k g / m ] × [ m / s ] = [ k g / s ] [\mathrm{kg/m}]\times[\mathrm{m/s}]=[\mathrm{kg/s}] [ kg/m ] × [ m/s ] = [ kg/s ] )。哈密顿密度为:
H = π y ˙ − L = μ y ˙ 2 − ( 1 2 μ y ˙ 2 − 1 2 T y ′ 2 ) = 1 2 μ y ˙ 2 + 1 2 T y ′ 2 \mathcal{H} = \pi\dot{y} - \mathcal{L} = \mu\dot{y}^2 - \left(\frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T y'^2\right) = \frac{1}{2}\mu\dot{y}^2 + \frac{1}{2}Ty'^2 H = π y ˙ − L = μ y ˙ 2 − ( 2 1 μ y ˙ 2 − 2 1 T y ′2 ) = 2 1 μ y ˙ 2 + 2 1 T y ′2
用 π \pi π 表示:H = π 2 / ( 2 μ ) + T y ′ 2 / 2 \mathcal{H} = \pi^2/(2\mu) + Ty'^2/2 H = π 2 / ( 2 μ ) + T y ′2 /2 。前一项是动能密度,后一项是弹性势能密度,总哈密顿量 H = ∫ H d x H = \int\mathcal{H}\,dx H = ∫ H d x 就是弦的总机械能,在无耗散情形下守恒。
应力-能量张量与守恒定理
对于时空均匀的系统,存在一个称为应力-能量张量 (或能量-动量张量)的量 T μ ν T^{\mu\nu} T μν ,定义为:
T μ ν = ∂ L ∂ ( ∂ μ ϕ ) ∂ ν ϕ − η μ ν L T^{\mu\nu} = \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}\partial^\nu\phi - \eta^{\mu\nu}\mathcal{L} T μν = ∂ ( ∂ μ ϕ ) ∂ L ∂ ν ϕ − η μν L
它是一个 4 × 4 4\times4 4 × 4 张量,包含了系统的能量密度、动量密度以及应力(动量流密度)等全部信息。
当场方程成立时,可以验证 T μ ν T^{\mu\nu} T μν 满足守恒方程:
∂ μ T μ ν = 0 \partial_\mu T^{\mu\nu} = 0 ∂ μ T μν = 0
这四个方程(ν = 0 , 1 , 2 , 3 \nu = 0, 1, 2, 3 ν = 0 , 1 , 2 , 3 )分别对应能量守恒和三个方向的动量守恒。
弦的能量密度 H = 1 2 μ y ˙ 2 + 1 2 T y ′ 2 \mathcal{H} = \frac{1}{2}\mu\dot{y}^2 + \frac{1}{2}Ty'^2 H = 2 1 μ y ˙ 2 + 2 1 T y ′2 。对时间求导:
∂ H ∂ t = μ y ˙ y ¨ + T y ′ y ˙ ′ = y ˙ ( μ y ¨ − T y ′ ′ ) + ∂ ∂ x ( T y ′ y ˙ ) \frac{\partial\mathcal{H}}{\partial t} = \mu\dot{y}\ddot{y} + Ty'\dot{y}' = \dot{y}(\mu\ddot{y} - Ty'') + \frac{\partial}{\partial x}(Ty'\dot{y}) ∂ t ∂ H = μ y ˙ y ¨ + T y ′ y ˙ ′ = y ˙ ( μ y ¨ − T y ′′ ) + ∂ x ∂ ( T y ′ y ˙ )
利用波动方程 μ y ¨ = T y ′ ′ \mu\ddot{y} = Ty'' μ y ¨ = T y ′′ ,括号内第一项为零,得到:
∂ H ∂ t + ∂ ∂ x ( − T y ′ y ˙ ) = 0 \frac{\partial\mathcal{H}}{\partial t} + \frac{\partial}{\partial x}\!\left(-Ty'\dot{y}\right) = 0 ∂ t ∂ H + ∂ x ∂ ( − T y ′ y ˙ ) = 0
这是能量守恒的局域形式:能量密度的时间变化率等于能量流 J E = − T y ′ y ˙ J_E = -Ty'\dot{y} J E = − T y ′ y ˙ (单位:W / m 2 \mathrm{W/m^2} W/ m 2 )在空间上的散度的负值。
应力-能量张量守恒方程 ∂ μ T μ ν = 0 \partial_\mu T^{\mu\nu} = 0 ∂ μ T μν = 0 是时空均匀性的直接后果:时间平移不变性给出能量守恒,空间平移不变性给出动量守恒。这一联系正是下一节诺特定理的具体体现。
诺特定理:对称性与守恒律
1915年,德国数学家埃米·诺特证明了一条极为深刻的定理:对于每一种连续对称性,系统都对应一个守恒量。 这条定理将对称性与守恒律的关系从经验现象提升为数学必然。
设作用量在如下无穷小变换下不变:
x μ → x μ + ϵ δ x μ , ϕ → ϕ + ϵ δ ϕ x^\mu \to x^\mu + \epsilon\,\delta x^\mu, \quad \phi \to \phi + \epsilon\,\delta\phi x μ → x μ + ϵ δ x μ , ϕ → ϕ + ϵ δ ϕ
其中 ϵ \epsilon ϵ 是无穷小参数,δ x μ \delta x^\mu δ x μ 和 δ ϕ \delta\phi δ ϕ 描述变换的方向。诺特定理指出,必然存在一个守恒流 j μ j^\mu j μ ,满足:
∂ μ j μ = 0 \partial_\mu j^\mu = 0 ∂ μ j μ = 0
对应的守恒荷(charge)为:
Q = ∫ j 0 d 3 x = 常数 Q = \int j^0\,d^3x = \text{常数} Q = ∫ j 0 d 3 x = 常数
若拉格朗日密度 L \mathcal{L} L 不显含时间 t t t ,则作用量在时间平移 t → t + ϵ t \to t + \epsilon t → t + ϵ 下不变,对应 δ x 0 = 1 \delta x^0 = 1 δ x 0 = 1 ,δ ϕ = 0 \delta\phi = 0 δ ϕ = 0 。诺特流的时间分量正好就是哈密顿密度:
j 0 = π ϕ ˙ − L = H j^0 = \pi\dot{\phi} - \mathcal{L} = \mathcal{H} j 0 = π ϕ ˙ − L = H
守恒条件 ∂ μ j μ = 0 \partial_\mu j^\mu = 0 ∂ μ j μ = 0 正是前面推导的能量守恒方程,守恒荷 Q = ∫ H d 3 x = H Q = \int\mathcal{H}\,d^3x = H Q = ∫ H d 3 x = H 就是总能量。
复标量场 ϕ \phi ϕ 的拉格朗日密度 L = ∣ ∂ μ ϕ ∣ 2 − m 2 ∣ ϕ ∣ 2 \mathcal{L} = |\partial_\mu\phi|^2 - m^2|\phi|^2 L = ∣ ∂ μ ϕ ∣ 2 − m 2 ∣ ϕ ∣ 2 在整体相位变换 ϕ → e i α ϕ \phi \to e^{i\alpha}\phi ϕ → e i α ϕ 下不变(α \alpha α 为实常数)。对应的诺特流为:
j μ = i ( ϕ ∗ ∂ μ ϕ − ϕ ∂ μ ϕ ∗ ) j^\mu = i\!\left(\phi^*\partial^\mu\phi - \phi\partial^\mu\phi^*\right) j μ = i ( ϕ ∗ ∂ μ ϕ − ϕ ∂ μ ϕ ∗ )
守恒荷 Q = ∫ j 0 d 3 x Q = \int j^0\,d^3x Q = ∫ j 0 d 3 x 在量子场论中对应粒子数之差(粒子数减去反粒子数),物理上就是电荷守恒定律的经典对应。
诺特定理的深刻之处在于:它不仅说明某个守恒量存在,还给出了守恒量的明确表达式,以及守恒的局域形式(守恒流满足连续性方程)。对称性越基本,对应的守恒律就越普遍,这一思想贯穿了整个现代理论物理。
相对论场论基础
将连续场的分析力学框架与狭义相对论结合,便进入相对论场论的领域。相对论场论的核心原则是:拉格朗日密度 L \mathcal{L} L 必须是洛伦兹标量 ,即在洛伦兹变换下保持不变,从而保证场方程在所有惯性参考系中具有相同的形式。
根据场的洛伦兹变换性质,场可以分类:
电磁场由四维矢量势 A μ = ( φ / c , A ) A^\mu = (\varphi/c, \mathbf{A}) A μ = ( φ / c , A ) 描述(单位:V ⋅ s / m \mathrm{V\cdot s/m} V ⋅ s/m ),场强张量定义为:
F μ ν = ∂ μ A ν − ∂ ν A μ F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu F μν = ∂ μ A ν − ∂ ν A μ
F μ ν F_{\mu\nu} F μν 的分量包含了电场 E \mathbf{E} E 和磁场 B \mathbf{B} B (单位:V / m \mathrm{V/m} V/m 和 T \mathrm{T} T )。电磁场的拉格朗日密度为:
L = − 1 4 μ 0 F μ ν F μ ν = ε 0 2 ( E 2 − c 2 B 2 ) \mathcal{L} = -\frac{1}{4\mu_0}F_{\mu\nu}F^{\mu\nu} = \frac{\varepsilon_0}{2}\!\left(E^2 - c^2B^2\right) L = − 4 μ 0 1 F μν F μν = 2 ε 0 ( E 2 − c 2 B 2 )
将这个 L \mathcal{L} L 代入欧拉-拉格朗日场方程,直接得到麦克斯韦方程组的前两组(有源方程):
∂ μ F μ ν = μ 0 j ν \partial_\mu F^{\mu\nu} = \mu_0 j^\nu ∂ μ F μν = μ 0 j ν
而后两组(无源方程,即 ∇ ⋅ B = 0 \nabla\cdot\mathbf{B}=0 ∇ ⋅ B = 0 和法拉第定律)则是 F μ ν F_{\mu\nu} F μν 反对称性的自动结果,无需额外引入。
相对论场论的内容仅作概念性介绍,完整的量子场论涉及场的量子化(正则量子化或路径积分量子化),超出经典力学的范围。这里的目的是展示拉格朗日-哈密顿方法的普适性:无论是一根振动的弦,还是充斥整个空间的电磁场,都可以纳入同一套变分原理的框架。
练习题
1. 关于从离散系统过渡到连续系统,下列说法正确的是:
A. 连续系统的拉格朗日量等于拉格朗日密度 L \mathcal{L} L 本身,不需要对空间积分
B. 对 N N N 个质点组成的弦离散模型取极限 Δ x → 0 \Delta x \to 0 Δ x → 0 时,差分方程变为偏微分波动方程
C. 拉格朗日密度 L \mathcal{L} L 的单位与拉格朗日量 L L L 的单位相同,均为 J \mathrm{J} J
D. 连续弦的波速 v = T / μ v = \sqrt{T/\mu} v = T / μ 中,T T T 是绝对温度
显示答案 答案:B
分析:连续系统的拉格朗日量是拉格朗日密度对空间的积分 L = ∫ L d x L = \int\mathcal{L}\,dx L = ∫ L d x ,A 错误。拉格朗日密度 L \mathcal{L} L 的单位是 J / m \mathrm{J/m} J/m (一维)或 J / m 3 \mathrm{J/m^3} J/ m 3 (三维),而 L L L 的单位是 J \mathrm{J} J ,C 错误。弦波速公式中 T T T 是弦的张力(单位 N \mathrm{N} N ),不是温度,D 错误。B 正确:对离散弦模型取连续极限,有限差分 ( y i + 1 − 2 y i + y i − 1 ) / ( Δ x ) 2 (y_{i+1}-2y_i+y_{i-1})/(\Delta x)^2 ( y i + 1 − 2 y i + y i − 1 ) / ( Δ x ) 2 趋近于 ∂ 2 y / ∂ x 2 \partial^2 y/\partial x^2 ∂ 2 y / ∂ x 2 ,运动方程变为波动方程 μ ∂ 2 y / ∂ t 2 = T ∂ 2 y / ∂ x 2 \mu\partial^2 y/\partial t^2 = T\partial^2 y/\partial x^2 μ ∂ 2 y / ∂ t 2 = T ∂ 2 y / ∂ x 2 。
2. 将欧拉-拉格朗日场方程应用于拉格朗日密度 L = 1 2 ( ∂ μ ϕ ) ( ∂ μ ϕ ) \mathcal{L} = \frac{1}{2}(\partial_\mu\phi)(\partial^\mu\phi) L = 2 1 ( ∂ μ ϕ ) ( ∂ μ ϕ ) 时,得到的场方程是:
A. ϕ = 0 \phi = 0 ϕ = 0 (ϕ \phi ϕ 等于零)
B. ∇ 2 ϕ = 0 \nabla^2\phi = 0 ∇ 2 ϕ = 0 (泊松方程,仅含空间导数)
C. □ ϕ = 0 \Box\phi = 0 □ ϕ = 0 ,即 ∂ 2 ϕ / ∂ t 2 = c 2 ∇ 2 ϕ \partial^2\phi/\partial t^2 = c^2\nabla^2\phi ∂ 2 ϕ / ∂ t 2 = c 2 ∇ 2 ϕ (波动方程)
D. ∂ ϕ / ∂ t = 0 \partial\phi/\partial t = 0 ∂ ϕ / ∂ t = 0 (ϕ \phi ϕ 不随时间变化)
显示答案 答案:C
分析:对 L = 1 2 ( ∂ μ ϕ ) ( ∂ μ ϕ ) \mathcal{L} = \frac{1}{2}(\partial_\mu\phi)(\partial^\mu\phi) L = 2 1 ( ∂ μ ϕ ) ( ∂ μ ϕ ) ,有 ∂ L / ∂ ϕ = 0 \partial\mathcal{L}/\partial\phi = 0 ∂ L / ∂ ϕ = 0 ,∂ L / ∂ ( ∂ μ ϕ ) = ∂ μ ϕ \partial\mathcal{L}/\partial(\partial_\mu\phi) = \partial^\mu\phi ∂ L / ∂ ( ∂ μ ϕ ) = ∂ μ ϕ 。代入场方程 ∂ μ ( ∂ μ ϕ ) = 0 \partial_\mu(\partial^\mu\phi) = 0 ∂ μ ( ∂ μ ϕ ) = 0 ,即 □ ϕ = 0 \Box\phi = 0 □ ϕ = 0 ,就是无质量标量场满足的波动方程,C 正确。这里 □ = ∂ 2 / ∂ t 2 − c 2 ∇ 2 \Box = \partial^2/\partial t^2 - c^2\nabla^2 □ = ∂ 2 / ∂ t 2 − c 2 ∇ 2 是达朗贝尔算符,包含时间和空间的二阶导数,B(只有空间导数)和 D(只有时间导数)均不完整。
A. 诺特定理仅适用于有限自由度系统,不适用于场论
B. 空间平移不变性对应能量守恒,时间平移不变性对应动量守恒
C. 每一种连续对称性对应一个守恒流 j μ j^\mu j μ ,满足 ∂ μ j μ = 0 \partial_\mu j^\mu = 0 ∂ μ j μ = 0 ,对应守恒荷 Q = ∫ j 0 d 3 x Q = \int j^0 d^3x Q = ∫ j 0 d 3 x
D. 规范对称性(相位变换不变性)对应角动量守恒
显示答案 答案:C
分析:诺特定理对有限自由度系统和场论都成立,A 错误。能量守恒来自时间平移不变性,动量守恒来自空间平移不变性,B 的说法正好相反,B 错误。规范对称性(如复场相位变换)对应电荷守恒,不是角动量守恒,D 错误。C 正确地描述了诺特定理的内容:连续对称性 ⇒ \Rightarrow ⇒ 守恒流 ∂ μ j μ = 0 \partial_\mu j^\mu = 0 ∂ μ j μ = 0 ⇒ \Rightarrow ⇒ 守恒荷 Q Q Q 为常数。
4. 弦的正则动量密度 π = ∂ L / ∂ y ˙ = μ y ˙ \pi = \partial\mathcal{L}/\partial\dot{y} = \mu\dot{y} π = ∂ L / ∂ y ˙ = μ y ˙ ,哈密顿密度 H = π 2 / ( 2 μ ) + T ( y ′ ) 2 / 2 \mathcal{H} = \pi^2/(2\mu) + T(y')^2/2 H = π 2 / ( 2 μ ) + T ( y ′ ) 2 /2 。下列关于弦的哈密顿量的说法,正确的是:
A. 总哈密顿量 H = ∫ H d x H = \int\mathcal{H}\,dx H = ∫ H d x 表示弦的总势能,不包含动能
B. H \mathcal{H} H 中的 π 2 / ( 2 μ ) \pi^2/(2\mu) π 2 / ( 2 μ ) 项对应弦的横向动能密度,T ( y ′ ) 2 / 2 T(y')^2/2 T ( y ′ ) 2 /2 对应弹性势能密度
C. 弦的哈密顿量 H H H 不守恒,因为弦是连续系统
D. 正则动量密度 π \pi π 的单位与广义动量(k g ⋅ m / s \mathrm{kg\cdot m/s} kg ⋅ m/s )相同
显示答案 答案:B
分析:哈密顿密度 H = π 2 / ( 2 μ ) + T ( y ′ ) 2 / 2 \mathcal{H} = \pi^2/(2\mu) + T(y')^2/2 H = π 2 / ( 2 μ ) + T ( y ′ ) 2 /2 中,第一项 = μ y ˙ 2 / 2 = \mu\dot{y}^2/2 = μ y ˙ 2 /2 是单位长度的动能(动能密度),第二项是单位长度的弹性势能(势能密度),总哈密顿量是两者之和,A 错误。在无外力和无耗散情形下,弦的总机械能守恒,C 错误。正则动量密度 π = μ y ˙ \pi = \mu\dot{y} π = μ y ˙ 的单位是 k g / ( m ⋅ s ) \mathrm{kg/(m\cdot s)} kg/ ( m ⋅ s ) (即动量除以长度),而不是 k g ⋅ m / s \mathrm{kg\cdot m/s} kg ⋅ m/s ,D 错误。B 正确。
5. 一根两端固定、长度为 L L L (单位 m \mathrm{m} m )的弦,线密度 μ \mu μ (k g / m \mathrm{kg/m} kg/m ),张力 T T T (N \mathrm{N} N ),自然振动的驻波解形如 y n ( x , t ) = A n sin ( n π x / L ) cos ( ω n t + φ n ) y_n(x,t) = A_n\sin(n\pi x/L)\cos(\omega_n t + \varphi_n) y n ( x , t ) = A n sin ( nπ x / L ) cos ( ω n t + φ n ) ,其中 n = 1 , 2 , 3 , … n = 1,2,3,\ldots n = 1 , 2 , 3 , … 。(一)将驻波解代入波动方程 μ y ¨ = T y ′ ′ \mu\ddot{y} = Ty'' μ y ¨ = T y ′′ ,求各阶角频率 ω n \omega_n ω n ;(二)写出第 n n n 阶模式的拉格朗日密度 L n \mathcal{L}_n L n ,并计算整弦的拉格朗日量 L n = ∫ 0 L L n d x L_n = \int_0^L \mathcal{L}_n\,dx L n = ∫ 0 L L n d x ;(三)说明为何不同阶模式之间的拉格朗日量相互独立(提示:利用正弦函数的正交性)。
显示答案 解:
(一)将 y n = A n sin ( n π x / L ) cos ( ω n t + φ n ) y_n = A_n\sin(n\pi x/L)\cos(\omega_n t+\varphi_n) y n = A n sin ( nπ x / L ) cos ( ω n t + φ n ) 代入 μ y ¨ = T y ′ ′ \mu\ddot{y} = Ty'' μ y ¨ = T y ′′ :
− μ ω n 2 A n sin n π x L cos ( ω n t + φ n ) = − T ( n π L ) 2 A n sin n π x L cos ( ω n t + φ n ) -\mu\omega_n^2 A_n\sin\frac{n\pi x}{L}\cos(\omega_n t+\varphi_n) = -T\!\left(\frac{n\pi}{L}\right)^{\!2}A_n\sin\frac{n\pi x}{L}\cos(\omega_n t+\varphi_n) − μ ω n 2 A n sin L nπ x cos ( ω n t + φ n ) = − T ( L nπ ) 2 A n sin L nπ x cos ( ω n t + φ n ) 两边约去公因子,得:
ω n 2 = T μ ( n π L ) 2 ⟹ ω n = n π L T μ = n π v L \omega_n^2 = \frac{T}{\mu}\!\left(\frac{n\pi}{L}\right)^{\!2} \implies \omega_n = \frac{n\pi}{L}\sqrt{\frac{T}{\mu}} = \frac{n\pi v}{L} ω n 2 = μ T ( L nπ ) 2 ⟹ ω n = L nπ μ T = L nπ v 其中 v = T / μ v = \sqrt{T/\mu} v = T / μ (m / s \mathrm{m/s} m/s )为波速,ω n \omega_n ω n 的单位为 r a d / s \mathrm{rad/s} rad/s 。基频(n = 1 n=1 n = 1 )为 ω 1 = π v / L \omega_1 = \pi v/L ω 1 = π v / L ,各阶频率均为基频的整数倍。
(二)第 n n n 阶模式的横向速度和斜率分别为:
y ˙ n = − A n ω n sin n π x L sin ( ω n t + φ n ) , y n ′ = A n n π L cos n π x L cos ( ω n t + φ n ) \dot{y}_n = -A_n\omega_n\sin\frac{n\pi x}{L}\sin(\omega_n t+\varphi_n), \quad y_n' = A_n\frac{n\pi}{L}\cos\frac{n\pi x}{L}\cos(\omega_n t+\varphi_n) y ˙ n = − A n ω n sin L nπ x sin ( ω n t + φ n ) , y n ′ = A n L nπ cos L nπ x cos ( ω n t + φ n ) 拉格朗日密度:
L n = 1 2 μ y ˙ n 2 − 1 2 T ( y n ′ ) 2 \mathcal{L}_n = \frac{1}{2}\mu\dot{y}_n^2 - \frac{1}{2}T(y_n')^2 L n = 2 1 μ y ˙ n 2 − 2 1 T ( y n ′ ) 2 = 1 2 μ A n 2 ω n 2 sin 2 n π x L sin 2 ( ω n t + φ n ) − 1 2 T A n 2 ( n π L ) 2 cos 2 n π x L cos 2 ( ω n t + φ n ) = \frac{1}{2}\mu A_n^2\omega_n^2\sin^2\!\frac{n\pi x}{L}\sin^2(\omega_n t+\varphi_n) - \frac{1}{2}T A_n^2\!\left(\frac{n\pi}{L}\right)^{\!2}\cos^2\!\frac{n\pi x}{L}\cos^2(\omega_n t+\varphi_n) = 2 1 μ A n 2 ω n 2 sin 2 L nπ x sin 2 ( ω n t + φ n ) − 2 1 T A n 2 ( L nπ ) 2 cos 2 L nπ x cos 2 ( ω n t + φ n ) 对 x x x 从 0 0 0 到 L L L 积分,利用 ∫ 0 L sin 2 ( n π x / L ) d x = ∫ 0 L cos 2 ( n π x / L ) d x = L / 2 \int_0^L\sin^2(n\pi x/L)\,dx = \int_0^L\cos^2(n\pi x/L)\,dx = L/2 ∫ 0 L sin 2 ( nπ x / L ) d x = ∫ 0 L cos 2 ( nπ x / L ) d x = L /2 :
L n = 1 4 μ A n 2 ω n 2 L sin 2 ( ω n t + φ n ) − 1 4 T A n 2 ( n π L ) 2 L cos 2 ( ω n t + φ n ) L_n = \frac{1}{4}\mu A_n^2\omega_n^2 L\sin^2(\omega_n t+\varphi_n) - \frac{1}{4}T A_n^2\!\left(\frac{n\pi}{L}\right)^{\!2} L\cos^2(\omega_n t+\varphi_n) L n = 4 1 μ A n 2 ω n 2 L sin 2 ( ω n t + φ n ) − 4 1 T A n 2 ( L nπ ) 2 L cos 2 ( ω n t + φ n ) 由 ω n 2 = T ( n π / L ) 2 / μ \omega_n^2 = T(n\pi/L)^2/\mu ω n 2 = T ( nπ / L ) 2 / μ ,两系数的大小相等,故:
L n = μ A n 2 ω n 2 L 4 [ sin 2 ( ω n t + φ n ) − cos 2 ( ω n t + φ n ) ] = − μ A n 2 ω n 2 L 4 cos ( 2 ω n t + 2 φ n ) L_n = \frac{\mu A_n^2\omega_n^2 L}{4}\!\left[\sin^2(\omega_n t+\varphi_n) - \cos^2(\omega_n t+\varphi_n)\right] = -\frac{\mu A_n^2\omega_n^2 L}{4}\cos(2\omega_n t+2\varphi_n) L n = 4 μ A n 2 ω n 2 L [ sin 2 ( ω n t + φ n ) − cos 2 ( ω n t + φ n ) ] = − 4 μ A n 2 ω n 2 L cos ( 2 ω n t + 2 φ n ) (三)对于两个不同阶数 m ≠ n m \neq n m = n ,混合项中含有 ∫ 0 L sin ( m π x / L ) sin ( n π x / L ) d x = 0 \int_0^L\sin(m\pi x/L)\sin(n\pi x/L)\,dx = 0 ∫ 0 L sin ( mπ x / L ) sin ( nπ x / L ) d x = 0 和 ∫ 0 L cos ( m π x / L ) cos ( n π x / L ) d x = 0 \int_0^L\cos(m\pi x/L)\cos(n\pi x/L)\,dx = 0 ∫ 0 L cos ( mπ x / L ) cos ( nπ x / L ) d x = 0 (正弦和余弦函数族在 [ 0 , L ] [0,L] [ 0 , L ] 上满足正交性)。因此,总拉格朗日量 L = ∫ 0 L L d x L = \int_0^L\mathcal{L}\,dx L = ∫ 0 L L d x 在展开时,不同模式之间的交叉项全部消失,各阶模式的拉格朗日量直接相加,彼此独立——连续弦等效于无穷多个相互独立的谐振子。
6. 对于拉格朗日密度 L = 1 2 μ y ˙ 2 − 1 2 T ( y ′ ) 2 \mathcal{L} = \frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T(y')^2 L = 2 1 μ y ˙ 2 − 2 1 T ( y ′ ) 2 ,用诺特定理验证空间平移不变性对应动量守恒。(一)写出对应空间平移 x → x + ϵ x \to x + \epsilon x → x + ϵ 的诺特流;(二)写出对应的守恒方程;(三)计算弦的总广义动量 P = ∫ j 0 d x P = \int j^0\,dx P = ∫ j 0 d x 并说明其物理意义。
显示答案 解:
(一)空间平移 x → x + ϵ x \to x + \epsilon x → x + ϵ 对应 δ x = 1 \delta x = 1 δ x = 1 ,δ t = 0 \delta t = 0 δ t = 0 ,场的变化 δ y = − y ′ \delta y = -y' δy = − y ′ (场随坐标平移的变化)。
对于空间平移,诺特流(动量流密度)的时间分量为:
j 0 = ∂ L ∂ y ˙ δ y = π ⋅ ( − y ′ ) = − μ y ˙ ⋅ y ′ j^0 = \frac{\partial\mathcal{L}}{\partial\dot{y}}\delta y = \pi\cdot(-y') = -\mu\dot{y}\cdot y' j 0 = ∂ y ˙ ∂ L δy = π ⋅ ( − y ′ ) = − μ y ˙ ⋅ y ′ 空间分量(即动量流密度,即应力)为:
j 1 = ∂ L ∂ y ′ δ y − L ⋅ δ x = ( − T y ′ ) ( − y ′ ) − L = T ( y ′ ) 2 − 1 2 μ y ˙ 2 + 1 2 T ( y ′ ) 2 = 3 2 T ( y ′ ) 2 − 1 2 μ y ˙ 2 j^1 = \frac{\partial\mathcal{L}}{\partial y'}\delta y - \mathcal{L}\cdot\delta x = (-Ty')(-y') - \mathcal{L} = T(y')^2 - \frac{1}{2}\mu\dot{y}^2 + \frac{1}{2}T(y')^2 = \frac{3}{2}T(y')^2 - \frac{1}{2}\mu\dot{y}^2 j 1 = ∂ y ′ ∂ L δy − L ⋅ δ x = ( − T y ′ ) ( − y ′ ) − L = T ( y ′ ) 2 − 2 1 μ y ˙ 2 + 2 1 T ( y ′ ) 2 = 2 3 T ( y ′ ) 2 − 2 1 μ y ˙ 2 更标准的推导给出应力-能量张量分量 T 10 = − T y ′ y ˙ T^{10} = -Ty'\dot{y} T 10 = − T y ′ y ˙ (动量流)和 T 00 = H T^{00} = \mathcal{H} T 00 = H (能量密度),这里动量守恒的诺特流为:
j 0 = − μ y ˙ y ′ , j 1 = T 11 = 1 2 μ y ˙ 2 − 1 2 T ( y ′ ) 2 = − L j^0 = -\mu\dot{y}\,y', \quad j^1 = T^{11} = \frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T(y')^2 = -\mathcal{L} j 0 = − μ y ˙ y ′ , j 1 = T 11 = 2 1 μ y ˙ 2 − 2 1 T ( y ′ ) 2 = − L (二)守恒方程为 ∂ t j 0 + ∂ x j 1 = 0 \partial_t j^0 + \partial_x j^1 = 0 ∂ t j 0 + ∂ x j 1 = 0 ,即:
∂ ∂ t ( − μ y ˙ y ′ ) + ∂ ∂ x ( − L ) = 0 \frac{\partial}{\partial t}(-\mu\dot{y}\,y') + \frac{\partial}{\partial x}\!\left(-\mathcal{L}\right) = 0 ∂ t ∂ ( − μ y ˙ y ′ ) + ∂ x ∂ ( − L ) = 0 验证:展开左边第一项:
∂ ∂ t ( − μ y ˙ y ′ ) = − μ y ¨ y ′ − μ y ˙ y ˙ ′ \frac{\partial}{\partial t}(-\mu\dot{y}\,y') = -\mu\ddot{y}\,y' - \mu\dot{y}\,\dot{y}' ∂ t ∂ ( − μ y ˙ y ′ ) = − μ y ¨ y ′ − μ y ˙ y ˙ ′ 第二项:
∂ ∂ x ( − L ) = − ∂ L ∂ x ∣ 显含 − ∂ L ∂ y y ′ − ∂ L ∂ y ˙ y ˙ ′ − ∂ L ∂ y ′ y ′ ′ = μ y ˙ y ˙ ′ + T y ′ y ′ ′ \frac{\partial}{\partial x}(-\mathcal{L}) = -\frac{\partial\mathcal{L}}{\partial x}\bigg|_{\text{显含}} - \frac{\partial\mathcal{L}}{\partial y}y' - \frac{\partial\mathcal{L}}{\partial\dot{y}}\dot{y}' - \frac{\partial\mathcal{L}}{\partial y'}y'' = \mu\dot{y}\dot{y}' + Ty'y'' ∂ x ∂ ( − L ) = − ∂ x ∂ L 显含 − ∂ y ∂ L y ′ − ∂ y ˙ ∂ L y ˙ ′ − ∂ y ′ ∂ L y ′′ = μ y ˙ y ˙ ′ + T y ′ y ′′ 两项之和为 − μ y ¨ y ′ + T y ′ y ′ ′ = y ′ ( T y ′ ′ − μ y ¨ ) = 0 -\mu\ddot{y}y' + Ty'y'' = y'(T y'' - \mu\ddot{y}) = 0 − μ y ¨ y ′ + T y ′ y ′′ = y ′ ( T y ′′ − μ y ¨ ) = 0 (利用波动方程),守恒方程成立。
(三)总广义动量:
P = ∫ 0 L j 0 d x = − ∫ 0 L μ y ˙ y ′ d x = − μ ∫ 0 L ∂ y ∂ t ∂ y ∂ x d x P = \int_0^L j^0\,dx = -\int_0^L \mu\dot{y}\,y'\,dx = -\mu\int_0^L\frac{\partial y}{\partial t}\frac{\partial y}{\partial x}\,dx P = ∫ 0 L j 0 d x = − ∫ 0 L μ y ˙ y ′ d x = − μ ∫ 0 L ∂ t ∂ y ∂ x ∂ y d x P P P 的单位为 k g ⋅ m / s \mathrm{kg\cdot m/s} kg ⋅ m/s ,称为弦的“准动量”(pseudo-momentum)或“场动量”,它描述弦上振动波型在空间中整体平移所对应的守恒量。对于两端固定的弦(边界条件 y ( 0 , t ) = y ( L , t ) = 0 y(0,t) = y(L,t) = 0 y ( 0 , t ) = y ( L , t ) = 0 ),驻波情形下左行波和右行波的动量大小相等、方向相反,总场动量为零,与空间平移对称性在有限边界下受到破坏的情况相吻合。