第V章:运动方程——两种时间演化

📑 目录

第V章:运动方程——两种时间演化

同样的物理,不同的舞台。薛定谔让波函数起舞,海森堡让算符歌唱——而狄拉克证明,它们是同一出戏的两个视角。


前置知识:经典力学中的哈密顿-雅可比方程

在进入量子力学的时间演化之前,我们必须先理解经典力学中一个极为优雅的理论框架——哈密顿-雅可比(Hamilton-Jacobi)理论。它不仅是分析力学的高峰,更是理解量子力学经典极限的钥匙。

1. 从哈密顿方程到雅可比方程

经典力学的哈密顿方程为:

q˙i=Hpi,p˙i=Hqi\dot{q}_i = \frac{\partial H}{\partial p_i}, \quad \dot{p}_i = -\frac{\partial H}{\partial q_i}

这些是一组 2n2n 个一阶微分方程。雅可比提出一个大胆的问题:能否找到一个单一函数 S(q,t)S(q, t),使得所有运动方程都可以从这个函数导出?

答案是肯定的。定义哈密顿主函数 S(q,t)S(q, t) 满足:

St+H(q,Sq,t)=0\frac{\partial S}{\partial t} + H\left(q, \frac{\partial S}{\partial q}, t\right) = 0

这就是哈密顿-雅可比方程——一个关于 SS 的一阶偏微分方程。对于保守系统(HH 不含时),可以分离变量:

S(q,t)=W(q)EtS(q, t) = W(q) - Et

其中 W(q)W(q) 满足哈密顿特征方程

H(q,Wq)=EH\left(q, \frac{\partial W}{\partial q}\right) = E

2. 作用量与经典路径

哈密顿主函数 SS 有深刻的物理意义:它等于沿经典路径计算的作用量积分

S(q,t;q0,t0)=t0tL(qcl(τ),q˙cl(τ),τ)dτS(q, t; q_0, t_0) = \int_{t_0}^{t} L(q_{cl}(\tau), \dot{q}_{cl}(\tau), \tau)\, d\tau

其中 qcl(τ)q_{cl}(\tau) 是满足欧拉-拉格朗日方程的经典路径,边界条件为 q(t0)=q0q(t_0) = q_0q(t)=qq(t) = q

关键性质:

  • Sq=p\frac{\partial S}{\partial q} = p(末态动量)
  • Sq0=p0\frac{\partial S}{\partial q_0} = -p_0(初态动量)
  • St=H\frac{\partial S}{\partial t} = -H(末态哈密顿量)
  • St0=H0\frac{\partial S}{\partial t_0} = H_0(初态哈密顿量)
graph TD
    A["经典路径 q_cl(τ)"] -->|"作用量积分"| B["S = ∫L dτ"]
    B --> C["∂S/∂q = p"]
    B --> D["∂S/∂t = -H"]
    C --> E["哈密顿-雅可比方程"]
    D --> E
    E --> F["∂S/∂t + H(q,∂S/∂q) = 0"]
    
    style A fill:#e3f2fd
    style B fill:#bbdefb
    style E fill:#90caf9
    style F fill:#64b5f6

3. 从经典到量子:Hamilton-Jacobi的量子类比

狄拉克和费曼发现,量子力学的经典极限与哈密顿-雅可比方程有深刻联系。薛定谔方程:

iψt=H^ψi\hbar\frac{\partial\psi}{\partial t} = \hat{H}\psi

ψ=eiS/\psi = e^{iS/\hbar},代入薛定谔方程,展开后取 0\hbar \to 0 极限,SS 恰好满足哈密顿-雅可比方程!

更精确地说,设 ψ=AeiS/\psi = A e^{iS/\hbar},其中 AASS 都是经典量(不依赖于 \hbar)。代入位置表象的薛定谔方程:

it(AeiS/)=[22m2+V]AeiS/i\hbar\frac{\partial}{\partial t}(A e^{iS/\hbar}) = \left[-\frac{\hbar^2}{2m}\nabla^2 + V\right]A e^{iS/\hbar}

展开左边:i(A˙+iAS˙/)eiS/=(AS˙+iA˙)eiS/i\hbar(\dot{A} + iA\dot{S}/\hbar)e^{iS/\hbar} = (-A\dot{S} + i\hbar\dot{A})e^{iS/\hbar}

展开右边动能项:

22m2(AeiS/)=22m[2A+2iAS+iA2SA2(S)2]eiS/-\frac{\hbar^2}{2m}\nabla^2(A e^{iS/\hbar}) = -\frac{\hbar^2}{2m}\left[\nabla^2 A + \frac{2i}{\hbar}\nabla A \cdot \nabla S + \frac{iA}{\hbar}\nabla^2 S - \frac{A}{\hbar^2}(\nabla S)^2\right]e^{iS/\hbar}

\hbar 的幂次整理:

0\hbar^0

St=12m(S)2+V-\frac{\partial S}{\partial t} = \frac{1}{2m}(\nabla S)^2 + V

这正是哈密顿-雅可比方程!(因为 p=Sp = \nabla S,所以 H=p2/2m+V=(S)2/2m+VH = p^2/2m + V = (\nabla S)^2/2m + V

1\hbar^1给出振幅 AA 的输运方程,描述了概率密度的经典演化。

核心洞察:量子力学中波函数的相位 SS 在经典极限下变成了经典作用量,而概率密度 A2|A|^2 在经典极限下收缩到经典轨迹附近。这就是WKB近似半经典量子力学的根基。

4. 经典极限的直观图像

想象一个自由粒子从 (x0,t0)(x_0, t_0) 传播到 (x,t)(x, t)。经典路径是直线匀速运动:xcl(τ)=x0+v(τt0)x_{cl}(\tau) = x_0 + v(\tau - t_0),其中 v=(xx0)/(tt0)v = (x-x_0)/(t-t_0)。经典作用量为:

Scl=t0t12mv2dτ=m(xx0)22(tt0)S_{cl} = \int_{t_0}^{t} \frac{1}{2}mv^2\, d\tau = \frac{m(x-x_0)^2}{2(t-t_0)}

量子力学的传播子(自由粒子)在 0\hbar \to 0 时:

K(x,t;x_0,t_0) = \sqrt{\frac{m}{2\pi i\hbar(t-t_0)}}\exp\left(\frac{im(x-x_0)^2}{2\hbar(t-t_0)}\right) \xrightarrow{\hbar\to 0} e^{iS_{cl}/\hbar}

指数因子精确还原了经典作用量!前置因子(平方根项)则来自路径积分中经典路径附近的高斯涨落——这是量子效应的经典前导修正。


时间银行的悖论

2154年,"永恒之城"是太阳系最大的时间交易中心。这里的商品不是金钱,而是时间状态。居民可以将自己的时间"存入"量子时间银行,选择以"个人视角"(薛定谔模式)或"事件视角"(海森堡模式)来体验时间的流逝。

年轻的交易员艾拉发现了一个悖论:在薛定谔模式下,她自身的"状态"不断变化——记忆累积,容颜渐改,而她所观测的物理量(如市场波动率)似乎是固定的算符。在海森堡模式下则相反——她的"状态"永远定格在开户那一刻,但周围的一切物理量都在随时间演化。

更令她困惑的是,无论选择哪种模式,她观测到的期望值——实际的交易结果——完全一致。"如果两种视角给出相同的现实,"艾拉在时间银行的咨询台前问道,"那么时间究竟是什么?是状态在流动,还是观测方式在流动?"

咨询台的AI沉默片刻,投影出一行字:"两种模式是同一幺正变换的镜像。时间不是流逝,而是关系的重组。"

这正是狄拉克在此章揭示的奥秘。


一、薛定谔绘景:态的舞蹈

1.1 时间演化算符

在薛定谔绘景中,量子态 ψ(t)|\psi(t)\rangle 随时间演化,算符保持不变。时间演化由时间演化算符 U^(t)\hat{U}(t) 描述:

ψ(t)=U^(t)ψ(0)|\psi(t)\rangle = \hat{U}(t)|\psi(0)\rangle

对于保守系统(哈密顿量 H^\hat{H} 不含时),时间演化算符为:

U^(t)=exp(iH^t)\hat{U}(t) = \exp\left(-\frac{i\hat{H}t}{\hbar}\right)

1.2 薛定谔方程

ψ(t)|\psi(t)\rangle 求时间导数:

itψ(t)=H^ψ(t)i\hbar\frac{\partial}{\partial t}|\psi(t)\rangle = \hat{H}|\psi(t)\rangle

这就是著名的薛定谔方程——量子力学的动力学核心。

在位置表象中,ψ(q,t)=qψ(t)\psi(q,t) = \langle q|\psi(t)\rangle,薛定谔方程变为:

iψt=22m2ψq2+V(q)ψi\hbar\frac{\partial\psi}{\partial t} = -\frac{\hbar^2}{2m}\frac{\partial^2\psi}{\partial q^2} + V(q)\psi

graph TD
    A["初始态 |ψ(0)⟩"] -->|"时间演化算符
Û(t) = e^(-iĤt/ℏ)"| B["t时刻态 |ψ(t)⟩"] B --> C["薛定谔方程
iℏ∂/∂t|ψ⟩ = Ĥ|ψ⟩"] C --> D["概率守恒
d/dt⟨ψ|ψ⟩ = 0"] D --> E["Û†Û = Î 保证幺正性"] style A fill:#e3f2fd style B fill:#bbdefb style C fill:#90caf9 style E fill:#64b5f6

1.3 概率守恒与幺正性

时间演化算符的幺正性(U^U^=I^\hat{U}^\dagger\hat{U} = \hat{I})直接保证了概率守恒:

ψ(t)ψ(t)=ψ(0)U^U^ψ(0)=ψ(0)ψ(0)=1\langle\psi(t)|\psi(t)\rangle = \langle\psi(0)|\hat{U}^\dagger\hat{U}|\psi(0)\rangle = \langle\psi(0)|\psi(0)\rangle = 1

概率不随时间改变——粒子不会凭空消失或出现。这是量子力学自洽性的基本保障。


二、海森堡绘景:算符的歌唱

2.1 视角的转换

在海森堡绘景中,态矢量保持固定:ψH=ψ(0)|\psi_H\rangle = |\psi(0)\rangle,而算符随时间演化:

A^H(t)=U^(t)A^SU^(t)=eiH^t/A^SeiH^t/\hat{A}_H(t) = \hat{U}^\dagger(t)\hat{A}_S\hat{U}(t) = e^{i\hat{H}t/\hbar}\hat{A}_S e^{-i\hat{H}t/\hbar}

其中下标 SS 表示薛定谔绘景,HH 表示海森堡绘景。

2.2 海森堡运动方程

A^H(t)\hat{A}_H(t) 求时间导数:

dA^Hdt=i[H^,A^H]+(A^t)H\frac{d\hat{A}_H}{dt} = \frac{i}{\hbar}[\hat{H}, \hat{A}_H] + \left(\frac{\partial \hat{A}}{\partial t}\right)_H

这就是海森堡运动方程。与薛定谔方程不同,它描述的是算符的演化,而非态的演化。

graph LR
    subgraph "薛定谔绘景"
        S1["态 |ψ(t)⟩ 随时间演化"]
        S2["算符 Â_S 固定"]
        S3["运动方程: iℏ∂|ψ⟩/∂t = Ĥ|ψ⟩"]
    end
    
    subgraph "海森堡绘景"
        H1["态 |ψ_H⟩ = |ψ(0)⟩ 固定"]
        H2["算符 Â_H(t) 随时间演化"]
        H3["运动方程: dÂ_H/dt = (i/ℏ)[Ĥ,Â_H]"]
    end
    
    S1 <-->|"#quot;期望等价
⟨ψ(t)"|Â_S|ψ(t)⟩ = ⟨ψ_H|Â_H(t)|ψ_H⟩"| H2 style S1 fill:#e8f5e9 style S2 fill:#c8e6c9 style H1 fill:#fff3e0 style H2 fill:#ffe0b2

2.3 与经典力学的惊人对应

海森堡运动方程与经典力学的哈密顿方程惊人地相似:

经典力学海森堡绘景(量子)
变量A(q,p)A(q,p)A^H(t)\hat{A}_H(t)
运动方程dAdt={A,H}\frac{dA}{dt} = \{A, H\}dA^Hdt=1i[A^H,H^]\frac{d\hat{A}_H}{dt} = \frac{1}{i\hbar}[\hat{A}_H, \hat{H}]
对易结构泊松括号 {A,H}\{A,H\}对易子 [A^H,H^]/(i)[\hat{A}_H, \hat{H}]/(i\hbar)

这正是狄拉克的对应原理在动力学层面的体现:将经典泊松括号替换为量子对易子(除以 ii\hbar),哈密顿方程就变成了海森堡方程。经典力学是量子力学在 0\hbar \to 0 极限下的近似。

graph TD
    A["经典力学"] -->|"哈密顿方程
dA/dt = {A,H}"| B["泊松括号代数"] B -->|"量子化:
{A,B} → [Â,B̂]/(iℏ)"| C["对易子代数"] C -->|"海森堡方程
dÂ/dt = (1/iℏ)[Â,Ĥ]"| D["量子海森堡绘景"] A -->|"ħ → 0 极限"| D style A fill:#c8e6c9 style B fill:#a5d6a7 style C fill:#90caf9 style D fill:#64b5f6

三、两种绘景的严格等价性

3.1 期望值的不变性

物理上可观测的只是算符的期望值。在两种绘景中:

  • 薛定谔绘景At=ψ(t)A^Sψ(t)\langle A\rangle_t = \langle\psi(t)|\hat{A}_S|\psi(t)\rangle
  • 海森堡绘景At=ψHA^H(t)ψH\langle A\rangle_t = \langle\psi_H|\hat{A}_H(t)|\psi_H\rangle

两者给出完全相同的数值:

ψ(t)A^Sψ(t)=ψ(0)U^A^SU^ψ(0)=ψHA^H(t)ψH\langle\psi(t)|\hat{A}_S|\psi(t)\rangle = \langle\psi(0)|\hat{U}^\dagger\hat{A}_S\hat{U}|\psi(0)\rangle = \langle\psi_H|\hat{A}_H(t)|\psi_H\rangle

3.2 等价性的深层证明

狄拉克在《量子力学原理》第V章中给出了比"期望值相等"更深刻的证明。两种绘景之间的联系是一个幺正变换

A^H=U^A^SU^,ψH=U^ψS(t)=ψS(0)\hat{A}_H = \hat{U}^\dagger \hat{A}_S \hat{U}, \quad |\psi_H\rangle = \hat{U}^\dagger |\psi_S(t)\rangle = |\psi_S(0)\rangle

所有代数结构在幺正变换下保持不变:

  • 对易子:[A^H,B^H]=U^[A^S,B^S]U^[\hat{A}_H, \hat{B}_H] = \hat{U}^\dagger [\hat{A}_S, \hat{B}_S] \hat{U}
  • 本征值谱:如果 A^Sa=aa\hat{A}_S |a\rangle = a|a\rangle,则 A^H(U^a)=a(U^a)\hat{A}_H (\hat{U}^\dagger|a\rangle) = a(\hat{U}^\dagger|a\rangle)
  • 谱分解:A^H=aaaHaH\hat{A}_H = \sum_a a |a\rangle_H\langle a|_H,其中 aH=U^aS|a\rangle_H = \hat{U}^\dagger |a\rangle_S

这意味着两种绘景不仅给出相同的期望值,而且在算符代数的所有结构上都完全等价。

3.3 表象无关的物理

两种绘景如同同一座山的两条登山路径。薛定谔路径让"观察者"移动,海森堡路径让"风景"变化,但到达任何观测点的"海拔高度"(期望值)完全一致。

graph TD
    A["物理现实"] --> B["薛定谔绘景"]
    A --> C["海森堡绘景"]
    A --> D["相互作用绘景"]
    
    B -->|"态演化
算符固定"| E["期望 ⟨ψ(t)|Â|ψ(t)⟩"] C -->|"算符演化
态固定"| E D -->|"自由部分算符演化
相互作用部分态演化"| E E --> F["同一物理结果"] style A fill:#f3e5f5 style E fill:#e1bee7 style F fill:#ce93d8

3.4 选择绘景的实用考量

虽然物理等价,不同绘景在不同问题中各有优势:

  • 薛定谔绘景:适合数值计算、散射问题、含时微扰
  • 海森堡绘景:适合与经典力学对比、格点场论、多体问题
  • 相互作用绘景:适合微扰论,自由部分用海森堡形式,相互作用部分用薛定谔形式

四、定态:能量本征的宁静

4.1 能量本征态的特殊角色

当系统处于哈密顿量的本征态 E|E\rangle 时:

H^E=EE\hat{H}|E\rangle = E|E\rangle

薛定谔方程给出简单解:

E(t)=eiEt/E|E(t)\rangle = e^{-iEt/\hbar}|E\rangle

4.2 定态的物理特征

定态有三重特殊之处:

  1. 概率分布不变ψ(q,t)2=ψ(q,0)2|\psi(q,t)|^2 = |\psi(q,0)|^2,粒子在空间中的分布不随时间改变
  2. 期望值不变:对任何不显含时的算符,At=A0\langle A\rangle_t = \langle A\rangle_0
  3. 能量确定:能量测量总是给出精确值 EE
graph TD
    A["能量本征态 |E⟩"] --> B["时间演化因子 e^(-iEt/ℏ)"]
    B --> C["概率密度不变 |ψ(q,t)|² = |ψ(q,0)|²"]
    B --> D["期望值不变 ⟨A⟩_t = ⟨A⟩_0"]
    B --> E["能量精确已知 ΔE = 0"]
    
    C --> F["定态: 粒子在空间中
概率分布静止"] D --> F E --> G["但 ΔE=0 意味着
Δt → ∞ (能量-时间不确定性)"] style A fill:#e8f5e9 style B fill:#c8e6c9 style F fill:#a5d6a7

4.3 定态与叠加态

一般的量子态是定态的叠加:

ψ(t)=ncneiEnt/En|\psi(t)\rangle = \sum_n c_n e^{-iE_nt/\hbar}|E_n\rangle

每个成分以各自的频率旋转,叠加后产生干涉,导致概率密度随时间变化——这就是量子动力学的丰富性来源。


五、自由粒子:最简单的舞蹈

5.1 自由粒子的哈密顿量

H^=p^22m\hat{H} = \frac{\hat{p}^2}{2m},没有势能项。能量本征态即动量本征态:

\hat{H}|p'\rangle = \frac{p'^2}{2m}|p'\rangle

5.2 波包的运动与扩散

一个局域化的波包是动量本征态的叠加:

\psi(q,t) = \int_{-\infty}^{\infty} \phi(p') e^{ip'q/\hbar} e^{-ip'^2t/(2m\hbar)} \frac{dp'}{\sqrt{2\pi\hbar}}

波包中心以群速度 vg=p0/mv_g = p_0/m 运动,同时扩散——初始位置越精确,扩散越快。

graph TD
    A["初始波包 ψ(q,0)"] -->|"动量空间展开"| B["φ(p') = ⟨p'|ψ(0)⟩"]
    B -->|"每个动量成分
以速度 p'/m 运动"| C["ψ(q,t) = ∫ φ(p')e^(i(p'q-Ept)/ℏ)dp'"] C --> D["波包中心运动"] C --> E["波包扩散"] D -->|"群速度 v_g = p₀/m"| F["经典轨迹 q = q₀ + (p₀/m)t"] E -->|"扩散速率∝(Δp)⁻¹"| G["初始位置越精确
扩散越快 (不确定性)"] style A fill:#e3f2fd style C fill:#bbdefb style F fill:#90caf9 style G fill:#ffccbc

5.3 波包扩散的定量分析

考虑一个初始高斯波包:

ψ(q,0)=(12πσ02)1/4exp((qq0)24σ02+ip0q)\psi(q,0) = \left(\frac{1}{2\pi\sigma_0^2}\right)^{1/4} \exp\left(-\frac{(q-q_0)^2}{4\sigma_0^2} + \frac{ip_0 q}{\hbar}\right)

其中 σ0\sigma_0 是初始位置不确定度,p0p_0 是初始平均动量。

自由粒子薛定谔方程的严格解给出:

ψ(q,t)=(12πσt2)1/4exp((qq0p0t/m)24σt2+ip0qip02t2m)\psi(q,t) = \left(\frac{1}{2\pi\sigma_t^2}\right)^{1/4} \exp\left(-\frac{(q-q_0-p_0t/m)^2}{4\sigma_t^2} + \frac{ip_0 q}{\hbar} - \frac{ip_0^2 t}{2m\hbar}\right)

其中时变宽度为:

σt=σ01+2t24m2σ04\sigma_t = \sigma_0\sqrt{1 + \frac{\hbar^2 t^2}{4m^2\sigma_0^4}}

物理意义

  • 波包中心以经典速度 p0/mp_0/m 匀速运动:qt=q0+p0mt\langle q \rangle_t = q_0 + \frac{p_0}{m}t
  • 波包宽度 σt\sigma_t 随时间增长,且增长速度反比于初始宽度 σ0\sigma_0
  • tτspread=2mσ02t \gg \tau_{spread} = \frac{2m\sigma_0^2}{\hbar} 时,σtt2mσ0\sigma_t \approx \frac{\hbar t}{2m\sigma_0},宽度与 tt 成正比

5.4 海森堡绘景中的自由粒子

在海森堡绘景中,位置和动量算符为:

q^H(t)=q^H(0)+p^H(0)mt\hat{q}_H(t) = \hat{q}_H(0) + \frac{\hat{p}_H(0)}{m}t

p^H(t)=p^H(0)\hat{p}_H(t) = \hat{p}_H(0)

这与经典力学的运动方程完全相同——自由粒子的动量守恒,位置匀速变化。


六、Ehrenfest定理:经典极限的桥梁

6.1 定理的陈述

Ehrenfest定理(1927年)给出了期望值如何遵循经典运动方程:

ddtq^=p^m,ddtp^=Vq^\frac{d}{dt}\langle \hat{q} \rangle = \frac{\langle \hat{p} \rangle}{m}, \quad \frac{d}{dt}\langle \hat{p} \rangle = -\left\langle \frac{\partial V}{\partial \hat{q}} \right\rangle

这组方程在结构上与经典哈密顿方程完全一致——只是把经典变量替换为期望值。

6.2 证明

从海森堡方程出发,对 q^\hat{q}p^\hat{p} 分别计算:

dq^Hdt=1i[q^,H^]=1i[q^,p^22m]=p^m\frac{d\hat{q}_H}{dt} = \frac{1}{i\hbar}[\hat{q}, \hat{H}] = \frac{1}{i\hbar}\left[\hat{q}, \frac{\hat{p}^2}{2m}\right] = \frac{\hat{p}}{m}

dp^Hdt=1i[p^,H^]=1i[p^,V(q^)]=Vq^\frac{d\hat{p}_H}{dt} = \frac{1}{i\hbar}[\hat{p}, \hat{H}] = \frac{1}{i\hbar}[\hat{p}, V(\hat{q})] = -\frac{\partial V}{\partial \hat{q}}

取期望值即得Ehrenfest定理。

6.3 经典极限的条件

然而,期望值遵循经典方程不等于量子系统表现为经典系统。关键条件是:

Vq^Vqq=q\left\langle \frac{\partial V}{\partial \hat{q}} \right\rangle \approx \frac{\partial V}{\partial q}\bigg|_{q=\langle q\rangle}

这只在以下情况成立:

  1. 势场变化缓慢:在波包尺度上 V(q)V(q) 近似线性(如宏观物体的重力场)
  2. 波包很窄Δq\Delta q \ll 势场特征变化尺度
  3. 谐振子特殊情况q3\langle q^3 \rangle 等奇次矩在相干态中恰好满足等号

对于强量子系统(如氢原子中的电子),势场在波包尺度上剧烈变化,上述近似完全失效——电子根本不沿任何经典轨道运动。

graph TD
    A["Ehrenfest定理"] --> B["d⟨q⟩/dt = ⟨p⟩/m"]
    A --> C["d⟨p⟩/dt = -⟨∂V/∂q⟩"]
    
    B --> D["经典极限条件:
⟨∂V/∂q⟩ ≈ ∂V/∂q|_{⟨q⟩}"] C --> D D -->|"满足"| E["经典行为
波包沿轨道运动"] D -->|"不满足"| F["量子行为
无经典对应"] E --> G["宏观物体
相干态"] F --> H["原子电子
势垒隧穿"] style A fill:#e3f2fd style D fill:#bbdefb style E fill:#c8e6c9 style F fill:#ffcdd2

七、作用量原理的路径积分形式

7.1 从算符到路径积分

狄拉克最早意识到时间演化算符的矩阵元与经典作用量的关系。费曼后来将其发展为路径积分形式。

跃迁振幅 qf,tfqi,ti\langle q_f, t_f|q_i, t_i\rangle 可以写成对所有路径的求和:

qf,tfqi,ti=D[q(t)]exp(ititfL(q,q˙)dt)\langle q_f, t_f|q_i, t_i\rangle = \int \mathcal{D}[q(t)] \exp\left(\frac{i}{\hbar}\int_{t_i}^{t_f} L(q,\dot{q})\,dt\right)

其中 LL 是经典拉格朗日量。

graph TD
    A["初始点 (qᵢ,tᵢ)"] -->|"所有可能路径"| B["最终点 (qƒ,tƒ)"]
    
    A --> C["路径1: e^(iS₁/ℏ)"]
    A --> D["路径2: e^(iS₂/ℏ)"]
    A --> E["路径N: e^(iSₙ/ℏ)"]
    
    C --> B
    D --> B
    E --> B
    
    B --> F["总振幅 = 对所有路径求和"]
    F --> G["经典极限 ℏ→0:
稳定相位 → 经典路径"] style A fill:#e8f5e9 style B fill:#c8e6c9 style F fill:#a5d6a7 style G fill:#ffd54f

7.2 经典极限

0\hbar \to 0 时,指数因子 eiS/e^{iS/\hbar} 振荡极快。相邻路径的贡献相消,只有使作用量 SS 取极值的经典路径附近的路径才有相长干涉。

路径积分揭示了经典力学的量子起源:经典路径是量子路径的干涉图样中最亮的条纹


八、吉布斯系综:统计力学的量子版本

8.1 密度矩阵

对于混合态(经典统计与量子不确定性的结合),引入密度矩阵

ρ^=ipiψiψi\hat{\rho} = \sum_i p_i |\psi_i\rangle\langle\psi_i|

其中 pip_i 是系统处于纯态 ψi|\psi_i\rangle 的经典概率。

8.2 量子刘维尔方程

密度矩阵的演化由量子版本的刘维尔方程描述:

iρ^t=[H^,ρ^]i\hbar\frac{\partial\hat{\rho}}{\partial t} = [\hat{H}, \hat{\rho}]

这与海森堡方程形式相似,但注意对易子的顺序:密度矩阵处于"中间"位置。

graph TD
    A["纯态 |ψ⟩"] -->|"量子不确定性"| B["投影 |ψ⟩⟨ψ|"]
    
    C["经典概率分布 ρ_classical"] -->|"经典统计不确定性"| D["ρ = Σ pᵢ|ψᵢ⟩⟨ψᵢ|"]
    
    B --> D
    
    D -->|"演化: iℏ∂ρ/∂t = [Ĥ,ρ]"| E["量子统计力学"]
    
    E -->|"热平衡: ρ ∝ e^(-Ĥ/kT)"| F["正则系综"]
    E -->|"巨正则系综"| G["量子巨配分函数"]
    
    style A fill:#e3f2fd
    style C fill:#fff8e1
    style D fill:#e8f5e9
    style E fill:#f3e5f5

8.3 量子配分函数

热平衡下的密度矩阵为 吉布斯分布

ρ^eq=1ZeH^/kBT\hat{\rho}_{eq} = \frac{1}{Z}e^{-\hat{H}/k_BT}

其中 Z=Tr(eH^/kBT)Z = \text{Tr}(e^{-\hat{H}/k_BT}) 是量子配分函数。迹(trace)运算给出了所有能量本征态的热权重求和:

Z=neEn/kBTZ = \sum_n e^{-E_n/k_BT}


数值例子

例1:自由电子波包的扩散时间尺度

问题:一个电子(质量 me=9.11×1031m_e = 9.11 \times 10^{-31} kg)被限制在初始位置不确定度为 σ0=1\sigma_0 = 1 nm 的波包中。计算波包宽度翻倍所需的时间。并讨论如果这个电子是原子中的价电子(被束缚在玻尔半径 a00.53a_0 \approx 0.53 Å),扩散意味着什么。

解答

波包宽度随时间的演化公式:

σt=σ01+2t24m2σ04\sigma_t = \sigma_0\sqrt{1 + \frac{\hbar^2 t^2}{4m^2\sigma_0^4}}

σt=2σ0\sigma_t = 2\sigma_0,即:

1+2t24m2σ04=2\sqrt{1 + \frac{\hbar^2 t^2}{4m^2\sigma_0^4}} = 2

1+2t24m2σ04=41 + \frac{\hbar^2 t^2}{4m^2\sigma_0^4} = 4

t=32mσ02=23mσ02t = \frac{\sqrt{3} \cdot 2m\sigma_0^2}{\hbar} = \frac{2\sqrt{3}m\sigma_0^2}{\hbar}

代入数值:

  • =1.055×1034\hbar = 1.055 \times 10^{-34} J·s
  • me=9.11×1031m_e = 9.11 \times 10^{-31} kg
  • σ0=1 nm=109\sigma_0 = 1\text{ nm} = 10^{-9} m

t=2×1.732×9.11×1031×(109)21.055×1034t = \frac{2 \times 1.732 \times 9.11 \times 10^{-31} \times (10^{-9})^2}{1.055 \times 10^{-34}}

t=3.464×9.11×1031×10181.055×1034=3.155×10481.055×10342.99×1014 st = \frac{3.464 \times 9.11 \times 10^{-31} \times 10^{-18}}{1.055 \times 10^{-34}} = \frac{3.155 \times 10^{-48}}{1.055 \times 10^{-34}} \approx 2.99 \times 10^{-14}\text{ s}

结果:波包宽度翻倍仅需约 30 飞秒3×10143 \times 10^{-14} s)。

物理讨论

  • 30飞秒是极短的时间尺度——原子振动周期约为10-100飞秒
  • 如果电子被束缚在玻尔半径 a0=0.53a_0 = 0.53 Å 内,初始 σ0a0\sigma_0 \approx a_0,则扩散时间:

ta0=23mea02=3.464×9.11×1031×(0.53×1010)21.055×1034t_{a_0} = \frac{2\sqrt{3}m_e a_0^2}{\hbar} = \frac{3.464 \times 9.11 \times 10^{-31} \times (0.53 \times 10^{-10})^2}{1.055 \times 10^{-34}}

=3.464×9.11×1031×2.81×10211.055×1034=8.86×10511.055×10348.4×1017 s= \frac{3.464 \times 9.11 \times 10^{-31} \times 2.81 \times 10^{-21}}{1.055 \times 10^{-34}} = \frac{8.86 \times 10^{-51}}{1.055 \times 10^{-34}} \approx 8.4 \times 10^{-17}\text{ s}

  • 84阿秒(8.4×10178.4 \times 10^{-17} s)是原子电子位置不确定度扩散的特征时间——这就是原子电子必须用量子力学描述的原因之一:在经典图像能形成之前,波包就已经弥散到整个原子尺度了。

例2:谐振子在两种绘景中的等价性数值验证

问题:一维谐振子,m=1m = 1 g,ω=1\omega = 1 rad/s,初始处于基态。在 t=1t = 1 s 时,计算位置算符的期望值 q^\langle \hat{q} \rangle 和方差 (Δq)2(\Delta q)^2。分别用薛定谔绘景和海森堡绘景计算,验证结果一致。

解答

谐振子参数

  • m=103m = 10^{-3} kg
  • ω=1\omega = 1 rad/s
  • =1.055×1034\hbar = 1.055 \times 10^{-34} J·s
  • 特征长度:a=mω=1.055×1034103×1=1.055×10313.25×1016a = \sqrt{\frac{\hbar}{m\omega}} = \sqrt{\frac{1.055 \times 10^{-34}}{10^{-3} \times 1}} = \sqrt{1.055 \times 10^{-31}} \approx 3.25 \times 10^{-16} m

这个宏观谐振子(1克质量,1 rad/s频率)的特征长度是核子尺度的十亿分之一——说明宏观谐振子处于极其"经典"的量子态。

薛定谔绘景

基态波函数:

ψ0(q)=(mωπ)1/4exp(mωq22)=1(πa2)1/4exp(q22a2)\psi_0(q) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\left(-\frac{m\omega q^2}{2\hbar}\right) = \frac{1}{(\pi a^2)^{1/4}}\exp\left(-\frac{q^2}{2a^2}\right)

时间演化:ψ(q,t)=eiE0t/ψ0(q)\psi(q,t) = e^{-iE_0 t/\hbar}\psi_0(q),其中 E0=12ω=12×1.055×1034×1=5.275×1035E_0 = \frac{1}{2}\hbar\omega = \frac{1}{2} \times 1.055 \times 10^{-34} \times 1 = 5.275 \times 10^{-35} J

位置期望值:

q^t=qψ(q,t)2dq=qψ0(q)2dq=0\langle \hat{q} \rangle_t = \int q |\psi(q,t)|^2 dq = \int q |\psi_0(q)|^2 dq = 0

(因为基态波函数是偶函数,qψ02q|\psi_0|^2 是奇函数,积分在对称区间上为零。)

位置方差:

(Δq)2=q2q2=q2(\Delta q)^2 = \langle q^2 \rangle - \langle q \rangle^2 = \langle q^2 \rangle

q2=q2ψ0(q)2dq=a22=2mω\langle q^2 \rangle = \int q^2 |\psi_0(q)|^2 dq = \frac{a^2}{2} = \frac{\hbar}{2m\omega}

=1.055×10342×103×1=5.275×1032 m2= \frac{1.055 \times 10^{-34}}{2 \times 10^{-3} \times 1} = 5.275 \times 10^{-32}\text{ m}^2

Δq=5.275×10322.30×1016 m\Delta q = \sqrt{5.275 \times 10^{-32}} \approx 2.30 \times 10^{-16}\text{ m}

海森堡绘景

海森堡绘景中的位置算符:

q^H(t)=q^H(0)cos(ωt)+p^H(0)mωsin(ωt)\hat{q}_H(t) = \hat{q}_H(0)\cos(\omega t) + \frac{\hat{p}_H(0)}{m\omega}\sin(\omega t)

对于谐振子,这是海森堡方程的精确解(因为 [H^,q^]=ip^/m[\hat{H}, \hat{q}] = -i\hbar\hat{p}/m[H^,p^]=imω2q^[\hat{H}, \hat{p}] = i\hbar m\omega^2 \hat{q},联立解得与经典谐振子相同形式的方程)。

取期望值:

q^H(t)=0q^H(0)0cos(ωt)+0p^H(0)0mωsin(ωt)\langle \hat{q}_H(t) \rangle = \langle 0|\hat{q}_H(0)|0\rangle\cos(\omega t) + \frac{\langle 0|\hat{p}_H(0)|0\rangle}{m\omega}\sin(\omega t)

由于基态是宇称本征态(偶宇称),而 q^\hat{q}p^\hat{p} 都是奇宇称算符:

0q^0=0,0p^0=0\langle 0|\hat{q}|0\rangle = 0, \quad \langle 0|\hat{p}|0\rangle = 0

因此:q^H(t)=0\langle \hat{q}_H(t) \rangle = 0,与薛定谔绘景结果一致。

位置方差:

(ΔqH)2=0q^H2(t)00q^H(t)02=0q^H2(t)0(\Delta q_H)^2 = \langle 0|\hat{q}_H^2(t)|0\rangle - \langle 0|\hat{q}_H(t)|0\rangle^2 = \langle 0|\hat{q}_H^2(t)|0\rangle

q^H2(t)=q^2(0)cos2(ωt)+p^2(0)(mω)2sin2(ωt)+q^(0)p^(0)+p^(0)q^(0)mωsin(ωt)cos(ωt)\hat{q}_H^2(t) = \hat{q}^2(0)\cos^2(\omega t) + \frac{\hat{p}^2(0)}{(m\omega)^2}\sin^2(\omega t) + \frac{\hat{q}(0)\hat{p}(0) + \hat{p}(0)\hat{q}(0)}{m\omega}\sin(\omega t)\cos(\omega t)

取基态期望值:

  • 0q^20=2mω=a22\langle 0|\hat{q}^2|0\rangle = \frac{\hbar}{2m\omega} = \frac{a^2}{2}
  • 0p^20=mω2\langle 0|\hat{p}^2|0\rangle = \frac{m\hbar\omega}{2}
  • 0(q^p^+p^q^)0=0\langle 0|(\hat{q}\hat{p} + \hat{p}\hat{q})|0\rangle = 0(对称化算符在定态中期望为零)

因此:

q^H2(t)=2mωcos2(ωt)+mω2(mω)2sin2(ωt)=2mω(cos2+sin2)=2mω\langle \hat{q}_H^2(t) \rangle = \frac{\hbar}{2m\omega}\cos^2(\omega t) + \frac{m\hbar\omega}{2(m\omega)^2}\sin^2(\omega t) = \frac{\hbar}{2m\omega}(\cos^2 + \sin^2) = \frac{\hbar}{2m\omega}

ΔqH=2mω=5.275×10322.30×1016 m\Delta q_H = \sqrt{\frac{\hbar}{2m\omega}} = \sqrt{5.275 \times 10^{-32}} \approx 2.30 \times 10^{-16}\text{ m}

验证结果:两种绘景给出完全相同的期望值和方差,数值均为 002.30×10162.30 \times 10^{-16} m,严格验证了绘景等价性。

物理讨论:对于这个宏观谐振子,量子零点涨落 Δq2.3×1016\Delta q \approx 2.3 \times 10^{-16} m 是原子核尺度(1015\sim 10^{-15} m)的约十分之一。这意味着宏观谐振子的量子效应在位置测量中几乎不可探测——这正是经典物理有效的原因。


本章总结

薛定谔绘景与海森堡绘景的等价性是量子力学最深刻的结果之一。

graph TD
    A["第V章核心结构"] --> B["薛定谔绘景"]
    A --> C["海森堡绘景"]
    A --> D["等价性证明"]
    
    B --> B1["态演化: |ψ(t)⟩ = e^(-iĤt/ℏ)|ψ(0)⟩"]
    B --> B2["薛定谔方程: iℏ∂|ψ⟩/∂t = Ĥ|ψ⟩"]
    B --> B3["算符固定"]
    
    C --> C1["态固定: |ψ_H⟩ = |ψ(0)⟩"]
    C --> C2["算符演化: Â_H(t) = e^(iĤt/ℏ)Âe^(-iĤt/ℏ)"]
    C --> C3["海森堡方程: dÂ/dt = (i/ℏ)[Ĥ,Â]"]
    
    D --> D1["期望相同: ⟨ψ(t)|Â|ψ(t)⟩ = ⟨ψ_H|Â_H(t)|ψ_H⟩"]
    D --> D2["经典对应: ħ→0 时海森堡→哈密顿"]
    
    A --> E["定态: |E(t)⟩ = e^(-iEt/ℏ)|E⟩"]
    A --> F["自由粒子: 波包运动与扩散"]
    A --> G["Ehrenfest定理: 经典极限"]
    A --> H["路径积分: 所有路径的干涉"]
    A --> I["吉布斯系综: 量子统计力学"]
    
    style B fill:#e3f2fd
    style C fill:#fff3e0
    style D fill:#e8f5e9
    style E fill:#f3e5f5
    style G fill:#ffecb3
    style H fill:#e1f5fe
    style I fill:#fce4ec

核心收获

  1. 哈密顿-雅可比方程是经典力学的高峰,其相位函数 SS 在量子力学的 0\hbar \to 0 极限下与波函数相位精确对应
  2. 薛定谔绘景中态随时间演化,海森堡绘景中算符随时间演化——两种描述通过幺正变换严格等价,不仅是期望值相等,整个算符代数结构都保持不变
  3. 海森堡方程 dA^/dt=(i/)[H^,A^]d\hat{A}/dt = (i/\hbar)[\hat{H},\hat{A}] 与经典哈密顿方程 dA/dt={A,H}dA/dt = \{A,H\} 在结构上一一对应,泊松括号 {A,B}\{A,B\} 被对易子 [A^,B^]/(i)[\hat{A},\hat{B}]/(i\hbar) 取代
  4. Ehrenfest定理表明期望值遵循经典方程,但经典行为要求势场在波包尺度上近似线性——这是宏观物体表现为经典的关键条件
  5. 定态是能量本征态,概率分布和期望值都不随时间改变;叠加态产生干涉导致动力学丰富性
  6. 自由粒子的波包以群速度运动并扩散,扩散速率与初始压缩度成反比——原子电子在10-100阿秒尺度上即弥散到整个原子尺度
  7. 路径积分将量子振幅表达为所有经典路径的贡献之和,经典极限对应稳定相位;传播子的前置因子来自经典路径附近的高斯涨落
  8. 密度矩阵将量子不确定性与经典统计不确定性统一,吉布斯系综推广了热力学,量子配分函数 Z=Tr(eH^/kBT)Z = \text{Tr}(e^{-\hat{H}/k_BT}) 将所有能级的热权重求和

练习与思考

1. 绘景的显式转换

考虑一维谐振子 H^=p^22m+12mω2q^2\hat{H} = \frac{\hat{p}^2}{2m} + \frac{1}{2}m\omega^2\hat{q}^2

(a) 在海森堡绘景中,证明 q^H(t)\hat{q}_H(t)p^H(t)\hat{p}_H(t) 满足的运动方程与经典谐振子的方程形式相同。

(b) 求解 q^H(t)\hat{q}_H(t)p^H(t)\hat{p}_H(t) 的显式表达式,验证它们与经典谐振子的解一致。

2. 不确定性关系的时间版本

能量-时间不确定性关系 ΔEΔt/2\Delta E \cdot \Delta t \geq \hbar/2 与位置-动量不确定性关系 ΔqΔp/2\Delta q \cdot \Delta p \geq \hbar/2 在数学结构上有何异同?

提示:时间 tt 在量子力学中不是算符,而是参数。"时间不确定性" Δt\Delta t 的含义是什么?考虑一个波包通过某一点所需的时间。

3. 路径积分的计算

利用路径积分公式计算自由粒子的传播子 K(qf,tf;qi,ti)=qf,tfqi,tiK(q_f,t_f;q_i,t_i) = \langle q_f,t_f|q_i,t_i\rangle

(a) 将时间间隔分成 NN 小段,每段 ϵ=(tfti)/N\epsilon = (t_f-t_i)/N

(b) 对每段插入完备集 pp|p\rangle\langle p|qq|q\rangle\langle q|

© 取 NN \to \infty 极限,证明结果为:

K=m2πi(tfti)exp(im(qfqi)22(tfti))K = \sqrt{\frac{m}{2\pi i\hbar(t_f-t_i)}}\exp\left(\frac{im(q_f-q_i)^2}{2\hbar(t_f-t_i)}\right)

并验证这个结果满足薛定谔方程。

4. Ehrenfest定理与维里定理

对于处于定态 ψn|\psi_n\rangle 的粒子,证明:

ψnq^Vq^ψn=2ψnT^ψn\langle \psi_n | \hat{q} \frac{\partial V}{\partial \hat{q}} | \psi_n \rangle = 2\langle \psi_n | \hat{T} | \psi_n \rangle

其中 T^=p^2/2m\hat{T} = \hat{p}^2/2m 是动能算符。这就是量子维里定理

(a) 对谐振子 V=12mω2q2V = \frac{1}{2}m\omega^2 q^2,验证 T=V=En/2\langle T \rangle = \langle V \rangle = E_n/2

(b) 对库仑势 V=e2/(4πϵ0r)V = -e^2/(4\pi\epsilon_0 r),验证 2T=V=En2\langle T \rangle = -\langle V \rangle = -E_n


"时间不是舞台,态也不是演员。时间是编织状态的织布机,而算符是丈量布匹的尺子。两种绘景不过是:一个人看织布机转动,另一个人看尺子移动——而布本身,从未改变。" —— 后记