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

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

---

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

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

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

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

$\dot{q}_i = \frac{\partial H}{\partial p_i}, \quad \dot{p}_i = -\frac{\partial H}{\partial q_i}$

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

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

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

这就是哈密顿-雅可比方程——一个关于 $S$ 的一阶偏微分方程。对于保守系统（$H$ 不含时），可以分离变量：

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

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

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

2. 作用量与经典路径

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

$S(q, t; q_0, t_0) = \int_{t_0}^{t} L(q_{cl}(\tau), \dot{q}_{cl}(\tau), \tau)\, d\tau$

其中 $q_{cl}(\tau)$ 是满足欧拉-拉格朗日方程的经典路径，边界条件为 $q(t_0) = q_0$，$q(t) = q$。

关键性质：

• $\frac{\partial S}{\partial q} = p$（末态动量）
• $\frac{\partial S}{\partial q_0} = -p_0$（初态动量）
• $\frac{\partial S}{\partial t} = -H$（末态哈密顿量）
• $\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\hbar\frac{\partial\psi}{\partial t} = \hat{H}\psi$

设 $\psi = e^{iS/\hbar}$，代入薛定谔方程，展开后取 $\hbar \to 0$ 极限，$S$ 恰好满足哈密顿-雅可比方程！

更精确地说，设 $\psi = A e^{iS/\hbar}$，其中 $A$ 和 $S$ 都是经典量（不依赖于 $\hbar$）。代入位置表象的薛定谔方程：

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

展开左边：$i\hbar(\dot{A} + iA\dot{S}/\hbar)e^{iS/\hbar} = (-A\dot{S} + i\hbar\dot{A})e^{iS/\hbar}$

展开右边动能项：

$-\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$ 的幂次整理：

$\hbar^0$ 阶：

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

这正是哈密顿-雅可比方程！（因为 $p = \nabla S$，所以 $H = p^2/2m + V = (\nabla S)^2/2m + V$）

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

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

4. 经典极限的直观图像

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

$S_{cl} = \int_{t_0}^{t} \frac{1}{2}mv^2\, d\tau = \frac{m(x-x_0)^2}{2(t-t_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 时间演化算符

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

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

对于保守系统（哈密顿量 $\hat{H}$ 不含时），时间演化算符为：

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

1.2 薛定谔方程

对 $|\psi(t)\rangle$ 求时间导数：

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

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

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

$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 概率守恒与幺正性

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

$\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 视角的转换

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

$\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}$

其中下标 $S$ 表示薛定谔绘景，$H$ 表示海森堡绘景。

2.2 海森堡运动方程

对 $\hat{A}_H(t)$ 求时间导数：

$\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 <-->|"期望等价
⟨ψ(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)$	$\hat{A}_H(t)$
运动方程	$\frac{dA}{dt} = \{A, H\}$	$\frac{d\hat{A}_H}{dt} = \frac{1}{i\hbar}[\hat{A}_H, \hat{H}]$
对易结构	泊松括号 $\{A,H\}$	对易子 $[\hat{A}_H, \hat{H}]/(i\hbar)$

这正是狄拉克的对应原理在动力学层面的体现：将经典泊松括号替换为量子对易子（除以 $i\hbar$），哈密顿方程就变成了海森堡方程。经典力学是量子力学在 $\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 期望值的不变性

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

• 薛定谔绘景：$\langle A\rangle_t = \langle\psi(t)|\hat{A}_S|\psi(t)\rangle$
• 海森堡绘景：$\langle A\rangle_t = \langle\psi_H|\hat{A}_H(t)|\psi_H\rangle$

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

$\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章中给出了比"期望值相等"更深刻的证明。两种绘景之间的联系是一个幺正变换：

$\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$

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

• 对易子：$[\hat{A}_H, \hat{B}_H] = \hat{U}^\dagger [\hat{A}_S, \hat{B}_S] \hat{U}$
• 本征值谱：如果 $\hat{A}_S |a\rangle = a|a\rangle$，则 $\hat{A}_H (\hat{U}^\dagger|a\rangle) = a(\hat{U}^\dagger|a\rangle)$
• 谱分解：$\hat{A}_H = \sum_a a |a\rangle_H\langle a|_H$，其中 $|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\rangle$ 时：

$\hat{H}|E\rangle = E|E\rangle$

薛定谔方程给出简单解：

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

4.2 定态的物理特征

定态有三重特殊之处：

1. 概率分布不变：$|\psi(q,t)|^2 = |\psi(q,0)|^2$，粒子在空间中的分布不随时间改变
2. 期望值不变：对任何不显含时的算符，$\langle A\rangle_t = \langle A\rangle_0$
3. 能量确定：能量测量总是给出精确值 $E$

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 定态与叠加态

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

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

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

---

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

5.1 自由粒子的哈密顿量

$\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}}

波包中心以群速度 $v_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 波包扩散的定量分析

考虑一个初始高斯波包：

$\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)$

其中 $\sigma_0$ 是初始位置不确定度，$p_0$ 是初始平均动量。

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

$\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)$

其中时变宽度为：

$\sigma_t = \sigma_0\sqrt{1 + \frac{\hbar^2 t^2}{4m^2\sigma_0^4}}$

物理意义：

• 波包中心以经典速度 $p_0/m$ 匀速运动：$\langle q \rangle_t = q_0 + \frac{p_0}{m}t$
• 波包宽度 $\sigma_t$ 随时间增长，且增长速度反比于初始宽度 $\sigma_0$
• 当 $t \gg \tau_{spread} = \frac{2m\sigma_0^2}{\hbar}$ 时，$\sigma_t \approx \frac{\hbar t}{2m\sigma_0}$，宽度与 $t$ 成正比

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

在海森堡绘景中，位置和动量算符为：

$\hat{q}_H(t) = \hat{q}_H(0) + \frac{\hat{p}_H(0)}{m}t$

$\hat{p}_H(t) = \hat{p}_H(0)$

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

---

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

6.1 定理的陈述

Ehrenfest定理（1927年）给出了期望值如何遵循经典运动方程：

$\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 证明

从海森堡方程出发，对 $\hat{q}$ 和 $\hat{p}$ 分别计算：

$\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}$

$\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 经典极限的条件

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

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

这只在以下情况成立：

1. 势场变化缓慢：在波包尺度上 $V(q)$ 近似线性（如宏观物体的重力场）
2. 波包很窄：$\Delta q \ll$ 势场特征变化尺度
3. 谐振子特殊情况：$\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 从算符到路径积分

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

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

$\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)$

其中 $L$ 是经典拉格朗日量。

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 经典极限

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

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

---

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

8.1 密度矩阵

对于混合态（经典统计与量子不确定性的结合），引入密度矩阵：

$\hat{\rho} = \sum_i p_i |\psi_i\rangle\langle\psi_i|$

其中 $p_i$ 是系统处于纯态 $|\psi_i\rangle$ 的经典概率。

8.2 量子刘维尔方程

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

$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 量子配分函数

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

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

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

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

---

数值例子

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

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

解答：

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

$\sigma_t = \sigma_0\sqrt{1 + \frac{\hbar^2 t^2}{4m^2\sigma_0^4}}$

令 $\sigma_t = 2\sigma_0$，即：

$\sqrt{1 + \frac{\hbar^2 t^2}{4m^2\sigma_0^4}} = 2$

$1 + \frac{\hbar^2 t^2}{4m^2\sigma_0^4} = 4$

$t = \frac{\sqrt{3} \cdot 2m\sigma_0^2}{\hbar} = \frac{2\sqrt{3}m\sigma_0^2}{\hbar}$

代入数值：

• $\hbar = 1.055 \times 10^{-34}$ J·s
• $m_e = 9.11 \times 10^{-31}$ kg
• $\sigma_0 = 1\text{ nm} = 10^{-9}$ m

$t = \frac{2 \times 1.732 \times 9.11 \times 10^{-31} \times (10^{-9})^2}{1.055 \times 10^{-34}}$

$t = \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 \times 10^{-14}$ s）。

物理讨论：

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

$t_{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}}$

$= \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 \times 10^{-17}$ s）是原子电子位置不确定度扩散的特征时间——这就是原子电子必须用量子力学描述的原因之一：在经典图像能形成之前，波包就已经弥散到整个原子尺度了。

---

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

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

解答：

谐振子参数：

• $m = 10^{-3}$ kg
• $\omega = 1$ rad/s
• $\hbar = 1.055 \times 10^{-34}$ J·s
• 特征长度：$a = \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频率）的特征长度是核子尺度的十亿分之一——说明宏观谐振子处于极其"经典"的量子态。

薛定谔绘景：

基态波函数：

$\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)$

时间演化：$\psi(q,t) = e^{-iE_0 t/\hbar}\psi_0(q)$，其中 $E_0 = \frac{1}{2}\hbar\omega = \frac{1}{2} \times 1.055 \times 10^{-34} \times 1 = 5.275 \times 10^{-35}$ J

位置期望值：

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

（因为基态波函数是偶函数，$q|\psi_0|^2$ 是奇函数，积分在对称区间上为零。）

位置方差：

$(\Delta q)^2 = \langle q^2 \rangle - \langle q \rangle^2 = \langle q^2 \rangle$

$\langle q^2 \rangle = \int q^2 |\psi_0(q)|^2 dq = \frac{a^2}{2} = \frac{\hbar}{2m\omega}$

$= \frac{1.055 \times 10^{-34}}{2 \times 10^{-3} \times 1} = 5.275 \times 10^{-32}\text{ m}^2$

$\Delta q = \sqrt{5.275 \times 10^{-32}} \approx 2.30 \times 10^{-16}\text{ m}$

海森堡绘景：

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

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

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

取期望值：

$\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)$

由于基态是宇称本征态（偶宇称），而 $\hat{q}$ 和 $\hat{p}$ 都是奇宇称算符：

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

因此：$\langle \hat{q}_H(t) \rangle = 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$

$\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)$

取基态期望值：

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

因此：

$\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}$

$\Delta q_H = \sqrt{\frac{\hbar}{2m\omega}} = \sqrt{5.275 \times 10^{-32}} \approx 2.30 \times 10^{-16}\text{ m}$

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

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

---

练习与思考

1. 绘景的显式转换

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

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

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

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

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

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

3. 路径积分的计算

利用路径积分公式计算自由粒子的传播子 $K(q_f,t_f;q_i,t_i) = \langle q_f,t_f|q_i,t_i\rangle$。

(a) 将时间间隔分成 $N$ 小段，每段 $\epsilon = (t_f-t_i)/N$。

(b) 对每段插入完备集 $|p\rangle\langle p|$ 和 $|q\rangle\langle q|$。

© 取 $N \to \infty$ 极限，证明结果为：

$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定理与维里定理

对于处于定态 $|\psi_n\rangle$ 的粒子，证明：

$\langle \psi_n | \hat{q} \frac{\partial V}{\partial \hat{q}} | \psi_n \rangle = 2\langle \psi_n | \hat{T} | \psi_n \rangle$

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

(a) 对谐振子 $V = \frac{1}{2}m\omega^2 q^2$，验证 $\langle T \rangle = \langle V \rangle = E_n/2$。

(b) 对库仑势 $V = -e^2/(4\pi\epsilon_0 r)$，验证 $2\langle T \rangle = -\langle V \rangle = -E_n$。

---

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