第9-10章 不确定性关系与N自由度系统:从精密测量到多体世界

📑 目录
  1. 故事引入:星际导航员的困境
  2. 9.0 前置知识:不等式家族
    1. 9.0.1 Cauchy-Schwarz不等式:内积空间的基本约束
    2. 9.0.2 三角不等式:不确定度的几何直觉
    3. 9.0.3 从不等式到不确定性:一步之遥
  3. 9.1 不确定性关系的严格证明
    1. 9.1.1 从直觉到严格:Schrödinger不等式
    2. 9.1.2 代数证明的核心步骤
    3. 9.1.3 数值例子:验证高斯波包的不确定性关系
    4. 9.1.4 历史插曲:玻尔与爱因斯坦的测量之争
  4. 9.2 最小不确定态与高斯波包
    1. 9.2.1 何时取等号?
    2. 9.2.2 高斯波包的物理图像
    3. 9.2.3 高斯波包的演化
    4. 9.2.4 数值例子:电子波包的扩散时间尺度
  5. 9.3 能量-时间不确定性
    1. 9.3.1 特殊的"不确定性"
    2. 9.3.2 物理实例:π⁰介子的衰变
    3. 9.3.3 数值例子:激发态原子的寿命与谱线宽度
  6. 10.0 前置知识:多粒子系统的经典图像
    1. 10.0.1 经典多粒子系统:相空间的维度
    2. 10.0.2 质心系:经典起源
    3. 10.0.3 从经典到量子:什么变了?
  7. 10.1 N自由度系统的张量积结构
    1. 10.1.1 从单粒子到多粒子:Hilbert空间的直积
    2. 10.1.2 算符的张量积结构
  8. 10.2 两粒子系统与质心系
    1. 10.2.1 相对坐标与质心坐标
    2. 10.2.2 数值例子:氢原子的约化质量修正
    3. 10.2.3 物理直觉:为什么质心系如此有效?
  9. 10.3 全同粒子与交换对称性
    1. 10.3.1 量子力学中的"身份危机"
    2. 10.3.2 玻色子与费米子
    3. 10.3.3 Slater行列式与对称化
    4. 10.3.4 数值例子:两电子系统的Slater行列式
  10. 10.4 多粒子波函数的具体构造
    1. 10.4.1 可分离态与纠缠态
    2. 10.4.2 从单粒子到多体:占据数表象
  11. 本章总结
  12. 练习与思考

第9-10章 不确定性关系与N自由度系统:从精密测量到多体世界

故事引入:星际导航员的困境

2147年,"深空七号"科考船正以0.3倍光速穿越奥尔特云。导航员林薇盯着主控台上不断闪烁的红色警报——飞船的量子定位系统陷入了悖论。要精确测定飞船在星际介质中的位置,光子探测器的波长必须极短;但波长越短,光子动量越大,每次测量都会把飞船推离预定轨道。反之,用长波光子测量可以减小扰动,但位置精度又随之下降。林薇苦笑着在日志中写道:"宇宙似乎设置了一道不可逾越的墙——我们永远无法同时精确知道自己在哪,以及将要去向何方。"这道墙,正是海森堡不确定性原理的本质。而在飞船货舱里,三千个低温冷却的量子比特正在并行运算,它们构成了一个庞大的N自由度系统——每一个比特的状态都依赖于整体的纠缠结构,单独描述任何一个都是徒劳。从单粒子的测量极限到多体系统的张量宇宙,量子力学正在向我们展示一幅既受限又无限丰富的图景。


9.0 前置知识:不等式家族

在深入不确定性关系的严格证明之前,让我们先建立一个共同的数学基础。Shankar在证明中依赖的核心工具都来自一个"不等式家族",它们从纯代数结构出发,最终导向物理上可观测的极限。

9.0.1 Cauchy-Schwarz不等式:内积空间的基本约束

对于任意内积空间中的两个向量 u|u\ranglev|v\rangle,Cauchy-Schwarz不等式表述为:

uv2uuvv|\langle u | v \rangle|^2 \leq \langle u | u \rangle \langle v | v \rangle

等号成立当且仅当 u|u\ranglev|v\rangle 线性相关,即 u=λv|u\rangle = \lambda |v\rangleλ\lambda 为复数)。

这个不等式是纯代数的——它不依赖于任何物理解释。在三维欧几里得空间中,它退化为 abab\mathbf{a} \cdot \mathbf{b} \leq |\mathbf{a}| |\mathbf{b}|,即两个向量夹角余弦的绝对值不超过1。在量子力学的Hilbert空间中,它成为不确定性关系的数学根基。

证明思路:考虑 uλv|u\rangle - \lambda |v\rangle 的模平方非负:

uλvuλv0\langle u - \lambda v | u - \lambda v \rangle \geq 0

展开后选择 λ=vu/vv\lambda = \langle v | u \rangle / \langle v | v \rangle(即投影系数),即可得到Cauchy-Schwarz不等式。这个证明如此简洁,以至于几乎不可能出错——这也是Shankar选择它作为不确定性关系基础的原因。

9.0.2 三角不等式:不确定度的几何直觉

从内积空间的度量结构,我们还有三角不等式:

u+vu+vuu+vv\sqrt{\langle u + v | u + v \rangle} \leq \sqrt{\langle u | u \rangle} + \sqrt{\langle v | v \rangle}

或者用范数记号:u+vu+v\|u + v\| \leq \|u\| + \|v\|

在不确定性关系的语境下,当我们考虑两个可观测量偏差的叠加态时,三角不等式提供了不确定度增长的上界。如果一个量子态同时对两个可观测量都有较大偏差,那么对它们的线性组合也必然有相应的不确定度。

9.0.3 从不等式到不确定性:一步之遥

Cauchy-Schwarz不等式是数学,不确定性关系是物理。两者之间的桥梁是对易子 [A,B]=ABBA[A, B] = AB - BA

考虑两个厄米算符 AABB(可观测量),以及任意归一化态 ψ|\psi\rangle。定义偏移态:

α=(AA)ψ,β=(BB)ψ|\alpha\rangle = (A - \langle A \rangle)|\psi\rangle, \quad |\beta\rangle = (B - \langle B \rangle)|\psi\rangle

直接应用Cauchy-Schwarz不等式:

(ΔA)2(ΔB)2=ααββαβ2(\Delta A)^2 (\Delta B)^2 = \langle \alpha | \alpha \rangle \langle \beta | \beta \rangle \geq |\langle \alpha | \beta \rangle|^2

αβ\langle \alpha | \beta \rangle 分解为实部和虚部,实部关联反对易子,虚部关联对易子。保留对易子项(因为反对易子项非负,丢弃它只会使不等式更保守),就得到Robertson不等式:

ΔAΔB12[A,B]\Delta A \, \Delta B \geq \frac{1}{2}|\langle [A, B] \rangle|

这就是从纯数学不等式到物理不确定性关系的全部推导。Shankar强调:这个过程没有任何"测量扰动"的参与——不确定性是量子态本身的内禀属性。

graph TD
    A[不等式家族] --> B["Cauchy-Schwarz
纯代数约束"] A --> C["三角不等式
几何约束"] B --> D["应用于偏移态
|α⟩, |β⟩"] D --> E["分解实部/虚部"] E --> F["实部→反对易子
{A,B}"] E --> G["虚部→对易子
[A,B]"] G --> H["Robertson不等式
ΔAΔB ≥ ½|⟨[A,B]⟩|"] H --> I["位置-动量
ΔXΔP ≥ ℏ/2"] H --> J["能量-时间
ΔEΔt ≥ ℏ/2"] style A fill:#e3f2fd style H fill:#fff8e1 style I fill:#ffebee

9.1 不确定性关系的严格证明

9.1.1 从直觉到严格:Schrödinger不等式

海森堡在1927年提出不确定性原理时,其论证基于思想实验——γ射线显微镜观测电子位置时,光子动量转移带来的扰动。这是一个极具物理直觉的论证,但严格性不足。Shankar在这一章给出了完全基于量子力学公理的代数证明,不依赖任何具体的测量装置。

核心工具是Schrödinger不等式(也常称为Robertson-Schrödinger关系),它是不等式中最一般的形式。

AABB 是两个可观测量(厄米算符),对于任意量子态 ψ|\psi\rangle,定义:

ΔA=A2A2\Delta A = \sqrt{\langle A^2 \rangle - \langle A \rangle^2}

ΔB=B2B2\Delta B = \sqrt{\langle B^2 \rangle - \langle B \rangle^2}

其中 A=ψAψ\langle A \rangle = \langle \psi | A | \psi \rangle 是期望值。

Schrödinger证明的不等式为:

\boxed{(\Delta A)^2 (\Delta B)^2 \geq \left|\frac{1}{2i}\langle [A,B] \rangle\right|^2 + \left|\frac{1}{2}\langle \{A - \langle A \rangle, B - \langle B \rangle\} \rangle\right|^2}

这个公式看起来复杂,让我们拆解它:

  • [A,B]=ABBA[A,B] = AB - BA 是对易子,决定了两个量能否同时被精确测量
  • {A,B}=AB+BA\{A,B\} = AB + BA 是反对易子
  • 第一项是标准的Heisenberg不确定性,第二项是额外的"关联项"

如果我们只保留第一项,就得到更著名的Robertson不等式

ΔAΔB12[A,B]\Delta A \, \Delta B \geq \frac{1}{2} |\langle [A,B] \rangle|

对于位置 XX 和动量 PP,由于 [X,P]=i[X,P] = i\hbar,我们立刻得到:

ΔXΔP2\Delta X \, \Delta P \geq \frac{\hbar}{2}

9.1.2 代数证明的核心步骤

Shankar的证明极其优雅。定义态的偏移:

α=(AA)ψ,β=(BB)ψ|\alpha\rangle = (A - \langle A \rangle)|\psi\rangle, \quad |\beta\rangle = (B - \langle B \rangle)|\psi\rangle

这两个态分别代表"偏离期望值的波动"。由Cauchy-Schwarz不等式:

ααββαβ2\langle \alpha | \alpha \rangle \langle \beta | \beta \rangle \geq |\langle \alpha | \beta \rangle|^2

左边正是 (ΔA)2(ΔB)2(\Delta A)^2 (\Delta B)^2。右边展开后分为实部和虚部:

αβ=(AA)(BB)=12{AA,BB}+12i[A,B]\langle \alpha | \beta \rangle = \langle (A - \langle A \rangle)(B - \langle B \rangle) \rangle = \frac{1}{2}\langle \{A - \langle A \rangle, B - \langle B \rangle\} \rangle + \frac{1}{2i}\langle [A, B] \rangle

于是Schrödinger不等式自然浮现。整个证明仅用了内积的基本性质——没有测量、没有扰动、没有具体的物理装置。不确定性是量子态本身的内禀属性。

graph TD
    A[不确定性关系的来源] --> B{"是测量扰动吗?"}
    B -->|海森堡原始思想实验| C["γ射线显微镜"]
    B -->|Shankar严格证明| D[量子态内禀性质]
    C --> E[直觉但依赖具体装置]
    D --> F["基于Cauchy-Schwarz不等式"]
    F --> G[仅依赖算符对易关系]
    G --> H[与测量过程无关]
    H --> I[纯代数推导]
    I --> J["ΔA ΔB ≥ ½|⟨[A,B]⟩|"]
    style D fill:#e1f5fe
    style I fill:#fff3e0

9.1.3 数值例子:验证高斯波包的不确定性关系

让我们用一个具体的数值例子验证不确定性关系。考虑一个一维高斯波包:

ψ(x)=(12πσ2)1/4exp[(xx0)24σ2+ip0x]\psi(x) = \left(\frac{1}{2\pi \sigma^2}\right)^{1/4} \exp\left[-\frac{(x - x_0)^2}{4\sigma^2} + \frac{i p_0 x}{\hbar}\right]

计算步骤

1. 位置期望值

X=+xψ(x)2dx=x0\langle X \rangle = \int_{-\infty}^{+\infty} x |\psi(x)|^2 dx = x_0

2. 位置平方期望值

X2=+x212πσ2exp[(xx0)22σ2]dx=x02+σ2\langle X^2 \rangle = \int_{-\infty}^{+\infty} x^2 \sqrt{\frac{1}{2\pi \sigma^2}} \exp\left[-\frac{(x-x_0)^2}{2\sigma^2}\right] dx = x_0^2 + \sigma^2

(利用高斯分布的二阶矩公式)

3. 位置不确定度

ΔX=X2X2=σ\Delta X = \sqrt{\langle X^2 \rangle - \langle X \rangle^2} = \sigma

4. 动量空间波函数(Fourier变换):

ψ~(p)=2σ2π21/2exp[σ2(pp0)22+ix0(p0p)]\tilde{\psi}(p) = \sqrt{\frac{2\sigma^2}{\pi \hbar^2}}^{1/2} \exp\left[-\frac{\sigma^2 (p - p_0)^2}{\hbar^2} + \frac{i x_0 (p_0 - p)}{\hbar}\right]

这也是一个高斯型!

5. 动量不确定度

ΔP=2σ\Delta P = \frac{\hbar}{2\sigma}

6. 验证

ΔXΔP=σ2σ=2\Delta X \cdot \Delta P = \sigma \cdot \frac{\hbar}{2\sigma} = \frac{\hbar}{2}

恰好饱和不确定性下界!

具体数值:设 σ=1.0×1010\sigma = 1.0 \times 10^{-10} m(约1 Å,原子尺度),则:

  • ΔX=1.0×1010\Delta X = 1.0 \times 10^{-10} m
  • \Delta P = \frac{1.055 \times 10^{-34} \text{ J·s}}{2 \times 10^{-10} \text{ m}} = 5.27 \times 10^{-25} kg·m/s
  • 对应的动能不确定度:ΔE=(ΔP)22me(5.27×1025)22×9.11×10311.5×1019\Delta E = \frac{(\Delta P)^2}{2m_e} \approx \frac{(5.27 \times 10^{-25})^2}{2 \times 9.11 \times 10^{-31}} \approx 1.5 \times 10^{-19} J \approx 1 eV

这个数量级恰好是原子尺度电子运动的特征能量!

9.1.4 历史插曲:玻尔与爱因斯坦的测量之争

1927年索尔维会议上,爱因斯坦设计了一个精巧的思想实验:用弹簧悬挂的匣子释放光子,通过秤量匣子在释放前后的质量变化(利用 E=mc2E=mc^2)来确定光子能量,同时控制释放时间来挑战能量-时间不确定性。玻尔彻夜未眠,最终发现爱因斯坦自己的广义相对论帮了倒忙——引力红移效应恰好补偿了不确定性,使得 ΔEΔt/2\Delta E \, \Delta t \geq \hbar/2 依然成立。

这个故事的深层含义是:不确定性不是某种可以被"绕过"的技术限制,而是嵌入在物理定律的底层结构中。正如Shankar所说,尝试击败不确定性关系就像尝试找到一个比平面上的所有直线都短的曲线——这不是技术问题,而是数学不可能性。


9.2 最小不确定态与高斯波包

9.2.1 何时取等号?

Robertson不等式 ΔAΔB12[A,B]\Delta A \, \Delta B \geq \frac{1}{2}|\langle [A,B] \rangle| 何时取等号?答案藏在Cauchy-Schwarz不等式取等的条件中——当且仅当两个态成比例关系:

α=λβ|\alpha\rangle = \lambda |\beta\rangle

对于 XXPP,这意味着:

(XX)ψ=λ(PP)ψ(X - \langle X \rangle)|\psi\rangle = \lambda (P - \langle P \rangle)|\psi\rangle

取 \lambda = -i\hbar/(2\Delta X)^2 \cdot \text{(常数)},这个微分方程的解正是高斯波包

ψ(x)=(12π(ΔX)2)1/4exp[iPx(xX)24(ΔX)2]\psi(x) = \left(\frac{1}{2\pi (\Delta X)^2}\right)^{1/4} \exp\left[\frac{i\langle P \rangle x}{\hbar} - \frac{(x - \langle X \rangle)^2}{4(\Delta X)^2}\right]

让我们逐项解读这个公式:

  • 前面的常数是归一化因子,确保 ψ2dx=1\int |\psi|^2 dx = 1
  • exp[iPx]\exp\left[\frac{i\langle P \rangle x}{\hbar}\right] 是动量为 P\langle P \rangle 的平面波相位因子
  • exp[(xX)24(ΔX)2]\exp\left[-\frac{(x - \langle X \rangle)^2}{4(\Delta X)^2}\right] 是高斯型振幅,中心在 X\langle X \rangle,宽度由 ΔX\Delta X 决定

9.2.2 高斯波包的物理图像

高斯波包是量子力学中最接近经典"粒子"概念的态。它在位置空间中是一个钟形曲线:

ψ(x)2=12π(ΔX)2exp[(xX)22(ΔX)2]|\psi(x)|^2 = \sqrt{\frac{1}{2\pi (\Delta X)^2}} \exp\left[-\frac{(x - \langle X \rangle)^2}{2(\Delta X)^2}\right]

标准差 ΔX\Delta X 描述了位置的弥散程度。动量空间中,Fourier变换后仍然是高斯型:

ψ~(p)2=12π(ΔP)2exp[(pP)22(ΔP)2]|\tilde{\psi}(p)|^2 = \sqrt{\frac{1}{2\pi (\Delta P)^2}} \exp\left[-\frac{(p - \langle P \rangle)^2}{2(\Delta P)^2}\right]

ΔXΔP=/2\Delta X \, \Delta P = \hbar/2,恰好饱和不确定性下界!

graph LR
    subgraph 位置空间
        A1["高斯波包 ψ(x)"]
        B1["中心: ⟨X⟩"]
        C1["宽度: ΔX"]
        D1["|ψ(x)|² = 钟形曲线"]
    end
    subgraph Fourier变换
        E["⇅ 对偶变换"]
    end
    subgraph 动量空间
        A2["高斯波包 ψ̃(p)"]
        B2["中心: ⟨P⟩"]
        C2["宽度: ΔP"]
        D2["|ψ̃(p)|² = 钟形曲线"]
    end
    subgraph 不确定性
        F["ΔX · ΔP = ℏ/2"]
        G[最小不确定态]
    end
    A1 --> E
    E --> A2
    A1 --> F
    A2 --> F
    style F fill:#ffebee
    style E fill:#e8f5e9

9.2.3 高斯波包的演化

高斯波包不是定态,它会随时间演化。有趣的现象是:波包会扩散。初始宽度为 ΔX0\Delta X_0 的高斯波包,在自由空间中演化时间 tt 后宽度变为:

(ΔX(t))2=(ΔX0)2+(t2mΔX0)2(\Delta X(t))^2 = (\Delta X_0)^2 + \left(\frac{\hbar t}{2m \Delta X_0}\right)^2

这意味着:即使初始位置很确定,随着时间流逝,位置不确定性会增长。这是因为在动量空间中有弥散,不同动量的分量以不同速度传播,导致波包"散开"。

Shankar用了一个绝妙的比喻:想象一群跑步者从同一点出发,但速度不同。起初他们聚在一起,但随着时间推移,速度快的跑在前面,慢的落在后面,整个群体逐渐散开。量子波包的扩散正是同一现象——ΔP\Delta P 不为零意味着存在速度分布。

9.2.4 数值例子:电子波包的扩散时间尺度

设一个电子(me=9.11×1031m_e = 9.11 \times 10^{-31} kg)的初始波包宽度为 ΔX0=1.0×1010\Delta X_0 = 1.0 \times 10^{-10} m(原子尺度)。计算它扩散到宽度翻倍所需的时间:

ΔX(t)=2ΔX0(t2mΔX0)2=3(ΔX0)2\Delta X(t) = 2 \Delta X_0 \Rightarrow \left(\frac{\hbar t}{2m \Delta X_0}\right)^2 = 3(\Delta X_0)^2

t=23m(ΔX0)2=2×1.732×9.11×1031×10201.055×1034t = \frac{2\sqrt{3} m (\Delta X_0)^2}{\hbar} = \frac{2 \times 1.732 \times 9.11 \times 10^{-31} \times 10^{-20}}{1.055 \times 10^{-34}}

t3.0×1016 s=0.3 fst \approx 3.0 \times 10^{-16} \text{ s} = 0.3 \text{ fs}

这是一个飞秒量级的时间!说明原子尺度的电子波包扩散极快。但对于宏观尺度:

ΔX0=1.0×106\Delta X_0 = 1.0 \times 10^{-6} m(微米尺度):

t=23m(106)23.0×108 s=30 nst = \frac{2\sqrt{3} m (10^{-6})^2}{\hbar} \approx 3.0 \times 10^{-8} \text{ s} = 30 \text{ ns}

微米尺度的电子波包需要约30纳秒才扩散翻倍。这说明波包扩散与初始宽度的平方成正比——越"胖"的波包越稳定。


9.3 能量-时间不确定性

9.3.1 特殊的"不确定性"

能量-时间不确定性 ΔEΔt/2\Delta E \, \Delta t \geq \hbar/2 看起来和 ΔXΔP/2\Delta X \, \Delta P \geq \hbar/2 对称,但Shankar强调了一个微妙但关键的差异:时间不是可观测量。在量子力学中,位置 XX 是算符,有本征态和本征值;但时间只是一个参数,不是Hilbert空间中的算符。

那么 Δt\Delta t 到底是什么?Shankar给出了三种诠释:

  1. 外部参数诠释Δt\Delta t 是测量过程持续的时间。测量时间越长,能量测量越精确。
  2. 寿命诠释:若系统处于不稳定态,其寿命 τ\tau 与能量宽度 Γ\Gamma 满足 Γτ\Gamma \tau \sim \hbar
  3. 动力学诠释Δt=ΔA/dA/dt\Delta t = \Delta A / |d\langle A \rangle / dt|,即期望值变化一个标准差所需的时间。

第三种诠释最深刻。考虑一个可观测量的期望值变化率:

dAdt=i[H,A]\frac{d\langle A \rangle}{dt} = \frac{i}{\hbar}\langle [H, A] \rangle

结合Robertson不等式:

ΔHΔA12[H,A]=2dAdt\Delta H \, \Delta A \geq \frac{1}{2}|\langle [H, A] \rangle| = \frac{\hbar}{2}\left|\frac{d\langle A \rangle}{dt}\right|

整理得到:

ΔHΔAdA/dt2\Delta H \cdot \frac{\Delta A}{|d\langle A \rangle / dt|} \geq \frac{\hbar}{2}

定义 ΔtA=ΔA/dA/dt\Delta t_A = \Delta A / |d\langle A \rangle / dt| 为期望值变化一个"特征单位"所需的时间,则:

ΔEΔtA2\Delta E \, \Delta t_A \geq \frac{\hbar}{2}

这里的 ΔtA\Delta t_A 依赖于选择哪个观测量 AA 来"计时"。

graph TD
    A["能量-时间不确定性"] --> B{时间不是算符}
    B --> C["Δt的三种诠释"]
    C --> D[测量持续时间]
    C --> E[不稳定态寿命]
    C --> F[动力学诠释]
    F --> G["Δt_A = ΔA / |d⟨A⟩/dt|"]
    G --> H[系统演化特征时间]
    H --> I["能量越确定
演化越慢"] I --> J["定态: ΔE=0
Δt=∞"] D --> K["π⁰介子寿命
Γτ ≈ ℏ"] style F fill:#e3f2fd

9.3.2 物理实例:π⁰介子的衰变

中性π介子(π0\pi^0)通过电磁相互作用衰变为两个光子,寿命极短(约 8.4×10178.4 \times 10^{-17} 秒)。根据能量-时间不确定性,其质量有自然宽度:

\Gamma = \frac{\hbar}{\tau} = \frac{6.582 \times 10^{-16} \text{ eV·s}}{8.4 \times 10^{-17} \text{ s}} \approx 7.8 \text{ eV}

这意味着 π0\pi^0 的静止能量不是单一值 135.0 MeV135.0 \text{ MeV},而是一个以该值为中心、宽度约 7.8 eV7.8 \text{ eV} 的分布。实验上,确实观测到了这种展宽——它不是测量误差,而是量子力学的内禀属性。

9.3.3 数值例子:激发态原子的寿命与谱线宽度

氢原子的 2p1s2p \to 1s 跃迁发射Lyman-α\alpha光子(λ=121.6\lambda = 121.6 nm)。这个跃迁的寿命约为 τ1.6×109\tau \approx 1.6 \times 10^{-9} s(1.6 ns)。对应的能量宽度:

\Gamma = \frac{\hbar}{\tau} = \frac{6.582 \times 10^{-16} \text{ eV·s}}{1.6 \times 10^{-9} \text{ s}} \approx 4.1 \times 10^{-7} \text{ eV}

转换为频率宽度:

\Delta \nu = \frac{\Gamma}{h} = \frac{4.1 \times 10^{-7} \text{ eV}}{4.136 \times 10^{-15} \text{ eV·s}} \approx 100 \text{ MHz}

对应的波长展宽:

Δλ=λ2cΔν=(121.6×109)23×108×1085×1015 m=5 fm\Delta \lambda = \frac{\lambda^2}{c} \Delta \nu = \frac{(121.6 \times 10^{-9})^2}{3 \times 10^8} \times 10^8 \approx 5 \times 10^{-15} \text{ m} = 5 \text{ fm}

(注意:这里 10810^8 是将MHz转为Hz的近似)

这个展宽远小于光谱仪的分辨率,因此通常观测不到。但对于亚稳态(如 2s1s2s \to 1s 的双光子跃迁,寿命约0.12 s),宽度仅约 10810^{-8} eV,极端尖锐。


10.0 前置知识:多粒子系统的经典图像

在量子力学中,多自由度系统的数学结构是张量积。但在进入这个抽象构造之前,让我们先在经典力学中建立直觉。经典图像虽然不"正确",但它提供了我们理解多体问题的认知锚点。

10.0.1 经典多粒子系统:相空间的维度

在经典力学中,N个粒子在三维空间中运动,需要6N个变量描述:每个粒子有3个位置坐标和3个动量坐标。系统的状态是相空间中的一个点,相空间维度为6N。

相空间不是构型空间(仅位置)。构型空间是3N维的,但完整的动力学描述需要动量信息。这就是为什么Hamilton力学比Lagrange力学更适合向量子力学过渡——Hamilton力学直接在相空间中运作。

10.0.2 质心系:经典起源

对于两体问题,经典力学中引入质心坐标和相对坐标:

R=m1r1+m2r2m1+m2,r=r1r2\mathbf{R} = \frac{m_1 \mathbf{r}_1 + m_2 \mathbf{r}_2}{m_1 + m_2}, \quad \mathbf{r} = \mathbf{r}_1 - \mathbf{r}_2

总动能分解为:

T=12MR˙2+12μr˙2T = \frac{1}{2}M \dot{\mathbf{R}}^2 + \frac{1}{2}\mu \dot{\mathbf{r}}^2

其中 M=m1+m2M = m_1 + m_2μ=m1m2m1+m2\mu = \frac{m_1 m_2}{m_1 + m_2}。如果势能仅依赖于相对位置 V=V(r)V = V(r),则运动完全分离为质心的自由运动和相对运动在中心力场中的运动。

开普勒问题(行星绕太阳)正是这个结构的经典案例。太阳质量 MM_\odot \gg 行星质量 mm,所以 μm\mu \approx m,质心几乎与太阳中心重合。

10.0.3 从经典到量子:什么变了?

量子力学继承了经典多体问题的坐标结构,但关键差异在于:

  1. 状态描述:经典状态是相空间中的点;量子状态是Hilbert空间中的向量
  2. 可分离性:经典系统中,两个粒子的状态总是"可分离的"——每个粒子有独立的轨迹;量子系统中,粒子可以处于纠缠态
  3. 全同性:经典中相同粒子可以贴标签追踪;量子中全同粒子不可区分,波函数必须对称化或反对称化

这些差异不是技术性的,而是根本性的。量子多体系统的Hilbert空间是各粒子空间的张量积(不是直和),这为纠缠提供了数学舞台。


10.1 N自由度系统的张量积结构

10.1.1 从单粒子到多粒子:Hilbert空间的直积

当系统包含多个自由度(例如多个粒子,或一个粒子在三维空间中运动)时,其Hilbert空间不是简单相加,而是张量积(直积)。这是量子力学中最具构造性的数学结构之一。

对于两个独立子系统,其Hilbert空间为:

H=H1H2\mathcal{H} = \mathcal{H}_1 \otimes \mathcal{H}_2

{n1}\{|n\rangle_1\}H1\mathcal{H}_1 的基,{m2}\{|m\rangle_2\}H2\mathcal{H}_2 的基,则 H\mathcal{H} 的基为:

n1m2n,m|n\rangle_1 \otimes |m\rangle_2 \equiv |n, m\rangle

一般态可展开为:

ψ=n,mcnmn,m|\psi\rangle = \sum_{n,m} c_{nm} |n, m\rangle

对于N个粒子,每个在三维空间中运动:

H=H1(3D)H2(3D)HN(3D)\mathcal{H} = \mathcal{H}_1^{(3D)} \otimes \mathcal{H}_2^{(3D)} \otimes \cdots \otimes \mathcal{H}_N^{(3D)}

维度是指数增长的:dim(H)=()3N\dim(\mathcal{H}) = (\infty)^{3N}——这是量子多体问题困难性的根源。

10.1.2 算符的张量积结构

单粒子算符在张量积空间中的作用方式是局部的。例如,第i个粒子的动量算符:

P_i^{(x)} = I \otimes I \otimes \cdots \otimes \underbrace{(-i\hbar \frac{\partial}{\partial x_i})}_{\text{第i个位置}} \otimes \cdots \otimes I

总动量算符则是所有单粒子动量之和:

Ptotal=i=1NPiP_{\text{total}} = \sum_{i=1}^N P_i

总Hamiltonian通常包含动能项(单粒子)和势能项(可能涉及多粒子):

H=i=1NPi22mi+V(X1,X2,,XN)H = \sum_{i=1}^N \frac{P_i^2}{2m_i} + V(X_1, X_2, \ldots, X_N)

graph TD
    A[N自由度系统] --> B["单粒子空间 ℋ_i"]
    B --> C["总空间 = 张量积"]
    C --> D["ℋ = ℋ₁ ⊗ ℋ₂ ⊗ ... ⊗ ℋₙ"]
    D --> E[维度指数增长]
    E --> F[多体问题的困难根源]
    A --> G[算符的局部作用]
    G --> H[P_i 只作用于第i个粒子]
    H --> I["总动量 = ΣP_i"]
    I --> J[可分离 vs 纠缠态]
    J --> K["|ψ⟩ = |a⟩⊗|b⟩ 可分离"]
    J --> L["|ψ⟩ = Σcᵢⱼ|i⟩|j⟩ 纠缠"]
    style C fill:#f3e5f5
    style L fill:#ffebee

10.2 两粒子系统与质心系

10.2.1 相对坐标与质心坐标

对于两粒子系统,一个关键的简化是引入质心坐标和相对坐标:

R=m1X1+m2X2m1+m2,r=X1X2R = \frac{m_1 X_1 + m_2 X_2}{m_1 + m_2}, \quad r = X_1 - X_2

对应的共轭动量:

P=P1+P2,p=m2P1m1P2m1+m2P = P_1 + P_2, \quad p = \frac{m_2 P_1 - m_1 P_2}{m_1 + m_2}

总动能变成两项之和:

P122m1+P222m2=P22M+p22μ\frac{P_1^2}{2m_1} + \frac{P_2^2}{2m_2} = \frac{P^2}{2M} + \frac{p^2}{2\mu}

其中 M=m1+m2M = m_1 + m_2 是总质量,μ=m1m2m1+m2\mu = \frac{m_1 m_2}{m_1 + m_2}约化质量。如果势能只依赖于相对位置 V=V(r)V = V(r),则Hamiltonian完全分离:

H = \underbrace{\frac{P^2}{2M}}_{\text{质心运动}} + \underbrace{\left(\frac{p^2}{2\mu} + V(r)\right)}_{\text{相对运动}}

这是一个可分离变量的系统!总波函数可以写成:

Ψ(X1,X2)=Φ(R)ψ(r)\Psi(X_1, X_2) = \Phi(R) \cdot \psi(r)

质心部分 Φ(R)\Phi(R) 描述整体运动(通常可取为平面波),相对运动部分 ψ(r)\psi(r) 描述内部结构——氢原子、双原子分子等都采用这种分解。

10.2.2 数值例子:氢原子的约化质量修正

氢原子中,质子质量 mp=1.673×1027m_p = 1.673 \times 10^{-27} kg,电子质量 me=9.109×1031m_e = 9.109 \times 10^{-31} kg。

约化质量的精确计算

μ=mempme+mp=9.109×1031×1.673×10279.109×1031+1.673×1027\mu = \frac{m_e m_p}{m_e + m_p} = \frac{9.109 \times 10^{-31} \times 1.673 \times 10^{-27}}{9.109 \times 10^{-31} + 1.673 \times 10^{-27}}

μ=1.524×10571.6739×10279.104×1031 kg\mu = \frac{1.524 \times 10^{-57}}{1.6739 \times 10^{-27}} \approx 9.104 \times 10^{-31} \text{ kg}

与电子质量的比较

μme=9.1049.1090.999455\frac{\mu}{m_e} = \frac{9.104}{9.109} \approx 0.999455

μ0.999455me\mu \approx 0.999455 \, m_e

能级修正

Rydberg能量(用约化质量):

ER(μ)=μe422=μmemee4220.999455×13.6057 eV13.598 eVE_R(\mu) = \frac{\mu e^4}{2\hbar^2} = \frac{\mu}{m_e} \cdot \frac{m_e e^4}{2\hbar^2} \approx 0.999455 \times 13.6057 \text{ eV} \approx 13.598 \text{ eV}

与使用 mem_e 的结果相差约 0.054%0.054\%。这个差异在光谱学中是可以测量的,也是证明质子具有有限质量(而非无穷大)的量子力学效应之一。

验证展开式

μ=me(1memp)=me(19.109×10311.673×1027)=me(15.45×104)\mu = m_e \left(1 - \frac{m_e}{m_p}\right) = m_e \left(1 - \frac{9.109 \times 10^{-31}}{1.673 \times 10^{-27}}\right) = m_e (1 - 5.45 \times 10^{-4})

=0.999455me= 0.999455 \, m_e

与精确计算一致。

10.2.3 物理直觉:为什么质心系如此有效?

想象两个人在冰面上互相抛接一个球。如果你站在远处观察,会看到两个人一边整体滑动,一边彼此靠近又远离。但如果你坐在他们的质心上(假设可以),你会看到两个人在做纯粹的相对振动——整体平移被"扣除"了。

质心坐标变换正是这个思想的数学实现。在质心系中,两体问题完全等价于一个质量为 μ\mu 的假想粒子在势场 V(r)V(r) 中的运动。这是量子力学中最优美的简化之一。

graph LR
    subgraph 实验室系
        A1["粒子1: X₁, P₁"]
        B1["粒子2: X₂, P₂"]
        C1["H = P₁²/2m₁ + P₂²/2m₂ + V("\\"X₁-X₂\\"")"]
    end
    D[坐标变换]
    subgraph 质心系
        A2["质心: R, P"]
        B2["相对: r, p"]
        C2["H = P²/2M + p²/2μ + V(r)"]
        E2["可分离!"]
    end
    A1 --> D
    B1 --> D
    D --> A2
    D --> B2
    C2 --> E2
    style E2 fill:#e8f5e9

10.3 全同粒子与交换对称性

10.3.1 量子力学中的"身份危机"

在经典力学中,两个相同的台球(比如都是红球)虽然在所有物理属性上一样,但我们可以给它们贴标签——"球A"和"球B",追踪它们的轨迹。如果交换它们,我们得到了一个物理上不同但不可区分的状态。

但在量子力学中,粒子没有轨迹。两个电子在轨道重叠区域内,我们无法说"这是电子1,那是电子2"。交换两个全同粒子得到的不是新的物理状态,而是同一个状态最多差一个相位因子

设交换算符 P12P_{12} 定义为交换粒子1和2的所有坐标(包括自旋):

P12ψ(x1,x2)=ψ(x2,x1)P_{12} \psi(x_1, x_2) = \psi(x_2, x_1)

由于两次交换回到原态:

P122=IP_{12}^2 = I

所以 P12P_{12} 的本征值只能是 ±1\pm 1

P12ψ=±ψP_{12} |\psi\rangle = \pm |\psi\rangle

10.3.2 玻色子与费米子

实验上发现一个惊人规律:

  • 玻色子(自旋整数:s=0,1,2,s = 0, 1, 2, \ldots):交换对称 P12ψ=+ψP_{12}|\psi\rangle = +|\psi\rangle
  • 费米子(自旋半整数:s=1/2,3/2,s = 1/2, 3/2, \ldots):交换反对称 P12ψ=ψP_{12}|\psi\rangle = -|\psi\rangle

这个自旋-统计关系是量子场论(自旋-统计定理)的深刻结果,但在非相对论量子力学中作为实验事实引入。

对于费米子,反对称性直接导出Pauli不相容原理:两个费米子不能处于相同的单粒子态。因为如果 ψ=a1a2|\psi\rangle = |a\rangle_1 \otimes |a\rangle_2(两个粒子都在态 a|a\rangle),则:

P12ψ=+ψP_{12}|\psi\rangle = +|\psi\rangle

但费米子要求 ψ-|\psi\rangle,所以唯一可能是 ψ=0|\psi\rangle = 0——这样的态不存在!

10.3.3 Slater行列式与对称化

对于N个费米子,构造反对称波函数的标准方法是Slater行列式

ψ(x1,,xN)=1N!det(ϕ1(x1)ϕ1(x2)ϕ1(xN)ϕ2(x1)ϕ2(x2)ϕ2(xN)ϕN(x1)ϕN(x2)ϕN(xN))\psi(x_1, \ldots, x_N) = \frac{1}{\sqrt{N!}} \det\begin{pmatrix} \phi_1(x_1) & \phi_1(x_2) & \cdots & \phi_1(x_N) \\ \phi_2(x_1) & \phi_2(x_2) & \cdots & \phi_2(x_N) \\ \vdots & \vdots & \ddots & \vdots \\ \phi_N(x_1) & \phi_N(x_2) & \cdots & \phi_N(x_N) \end{pmatrix}

其中 ϕi\phi_i 是单粒子态。行列式的性质保证:

  • 交换两列(即两个粒子)引入负号——反对称
  • 两列相同则行列式为零——Pauli原理

对于玻色子,则使用永久式(Permanent,去掉行列式的符号权重),保证完全对称。

10.3.4 数值例子:两电子系统的Slater行列式

考虑氦原子的两个电子,假设它们分别占据1s和2s轨道(简化模型,忽略电子-电子相互作用)。单粒子波函数为:

ϕ1s(r)=1πa03er/a0,ϕ2s(r)=18πa03(1r2a0)er/(2a0)\phi_{1s}(r) = \frac{1}{\sqrt{\pi a_0^3}} e^{-r/a_0}, \quad \phi_{2s}(r) = \frac{1}{\sqrt{8\pi a_0^3}} \left(1 - \frac{r}{2a_0}\right) e^{-r/(2a_0)}

自旋部分:电子1自旋向上 |\uparrow\rangle,电子2自旋向下 |\downarrow\rangle。包含自旋的反对称波函数为:

Ψ(r1,s1;r2,s2)=12ϕ1s(r1)1ϕ1s(r2)2ϕ2s(r1)1ϕ2s(r2)2\Psi(\mathbf{r}_1, s_1; \mathbf{r}_2, s_2) = \frac{1}{\sqrt{2}} \begin{vmatrix} \phi_{1s}(\mathbf{r}_1)|\uparrow\rangle_1 & \phi_{1s}(\mathbf{r}_2)|\uparrow\rangle_2 \\ \phi_{2s}(\mathbf{r}_1)|\downarrow\rangle_1 & \phi_{2s}(\mathbf{r}_2)|\downarrow\rangle_2 \end{vmatrix}

=12[ϕ1s(r1)ϕ2s(r2)12ϕ1s(r2)ϕ2s(r1)21]= \frac{1}{\sqrt{2}} \left[\phi_{1s}(\mathbf{r}_1)\phi_{2s}(\mathbf{r}_2)|\uparrow\rangle_1|\downarrow\rangle_2 - \phi_{1s}(\mathbf{r}_2)\phi_{2s}(\mathbf{r}_1)|\uparrow\rangle_2|\downarrow\rangle_1\right]

验证交换反对称性:交换粒子1和2(r1r2\mathbf{r}_1 \leftrightarrow \mathbf{r}_2s1s2s_1 \leftrightarrow s_2):

P12Ψ=12[ϕ1s(r2)ϕ2s(r1)12ϕ1s(r1)ϕ2s(r2)21]P_{12}\Psi = \frac{1}{\sqrt{2}} \left[\phi_{1s}(\mathbf{r}_2)\phi_{2s}(\mathbf{r}_1)|\downarrow\rangle_1|\uparrow\rangle_2 - \phi_{1s}(\mathbf{r}_1)\phi_{2s}(\mathbf{r}_2)|\downarrow\rangle_2|\uparrow\rangle_1\right]

注意 12=21|\downarrow\rangle_1|\uparrow\rangle_2 = |\uparrow\rangle_2|\downarrow\rangle_1,第一项就是原来第二项的负号交换,所以:

P12Ψ=ΨP_{12}\Psi = -\Psi

验证成功!这个Slater行列式自动保证了反对称性。

graph TD
    A[全同粒子交换] --> B["两次交换 = 恒等变换"]
    B --> C["P₁₂² = I"]
    C --> D["本征值: ±1"]
    D --> E["玻色子: +1"]
    D --> F["费米子: -1"]
    E --> G["对称波函数
玻色-爱因斯坦统计"] F --> H["反对称波函数
费米-狄拉克统计"] H --> I[Pauli不相容原理] I --> J[电子壳层结构] J --> K["原子、分子、固体的
全部化学"] G --> L["玻色-爱因斯坦凝聚"] H --> M[Slater行列式] style I fill:#ffebee style K fill:#fff8e1

10.4 多粒子波函数的具体构造

10.4.1 可分离态与纠缠态

对于两粒子系统,一般态是:

ψ=i,jciji1j2|\psi\rangle = \sum_{i,j} c_{ij} |i\rangle_1 \otimes |j\rangle_2

如果 cijc_{ij} 可以分解为 cij=aibjc_{ij} = a_i b_j,则:

ψ=(iaii1)(jbjj2)=α1β2|\psi\rangle = \left(\sum_i a_i |i\rangle_1\right) \otimes \left(\sum_j b_j |j\rangle_2\right) = |\alpha\rangle_1 \otimes |\beta\rangle_2

这是可分离态(Product state)。两个子系统有独立的描述。

如果不能这样分解,则是纠缠态(Entangled state)。纠缠态最奇特的性质是:子系统没有独立的量子态——只有整个系统有纯态,每个子系统单独看都是混合态(由密度矩阵描述)。

Shankar著名的例子:两个自旋1/2粒子处于单态:

ψ=12()|\psi\rangle = \frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle - |\downarrow\uparrow\rangle)

测量粒子1的自旋:如果得到上,粒子2必定是下;如果得到下,粒子2必定是上。无论它们相隔多远,这种关联瞬时出现。这不是超光速信号(因为单独看每个粒子的测量结果都是随机的,只有比对数据才看到关联),但它挑战了我们对"局域实在论"的直觉。

10.4.2 从单粒子到多体:占据数表象

对于N个全同玻色子,由于对称性要求,可以用占据数 nin_i(态 ϕi\phi_i 中有多少个粒子)来标记基态:

n1,n2,,ini=N|n_1, n_2, \ldots\rangle, \quad \sum_i n_i = N

这引向了二次量子化(场论)的语言,但在非相对论量子力学中已经显现端倪。对于费米子,由于Pauli原理,ni{0,1}n_i \in \{0, 1\}

graph TD
    A[多粒子波函数] --> B{"是否可分离?"}
    B -->|是| C["直积态 |ψ⟩ = |α⟩⊗|β⟩"]
    C --> D[子系统独立描述]
    D --> E[经典关联]
    B -->|否| F[纠缠态]
    F --> G[子系统无纯态]
    G --> H[密度矩阵描述]
    H --> I[非经典关联]
    I --> J[EPR悖论]
    J --> K[Bell不等式违反]
    K --> L[量子信息革命]
    style F fill:#ffebee
    style L fill:#e8f5e9

本章总结

不确定性关系和N自由度系统是量子力学从单粒子向多体世界扩展的关键桥梁。不确定性关系告诉我们,精密测量的极限不是技术问题而是数学必然——它根植于Cauchy-Schwarz不等式和对易子的代数结构。高斯波包作为最小不确定态,是经典与量子世界的最佳接口:它同时具有确定的位置中心和确定的动量中心,只是两者的精度受到 /2\hbar/2 的约束。能量-时间不确定性的特殊性提醒我们,时间在非相对论量子力学中仍是一个外部参数,而非动力学变量。

当我们把目光从单个粒子转向N自由度系统,张量积结构揭示了多体量子世界的丰富性。两粒子问题在质心系中可分离为整体运动和相对运动,这是分析原子、分子散射的基石。全同粒子的交换对称性则引入了玻色子和费米子的根本区分,Pauli不相容原理主导了原子结构、化学键合和固体电子结构。从可分离态到纠缠态,量子关联超越了经典概率论的极限,为量子计算和量子通信埋下了伏笔。

graph TD
    A["第9-10章核心脉络"] --> B[不确定性关系]
    A --> C[N自由度系统]
    B --> D["Schrödinger不等式"]
    B --> E[高斯最小不确定态]
    B --> F["能量-时间诠释"]
    C --> G[张量积结构]
    C --> H[质心系简化]
    C --> I[全同粒子对称性]
    I --> J["玻色子 +1"]
    I --> K["费米子 -1"]
    K --> L[Pauli原理]
    L --> M[Slater行列式]
    G --> N[可分离态]
    G --> O[纠缠态]
    style B fill:#e3f2fd
    style C fill:#f3e5f5
    style O fill:#ffebee

练习与思考

  1. 最小不确定态的证明:证明高斯波包 ψ(x)=(2πσ2)1/4exp[(xx0)2/(4σ2)+ip0x/]\psi(x) = (2\pi\sigma^2)^{-1/4} \exp[-(x-x_0)^2/(4\sigma^2) + ip_0 x/\hbar] 确实满足 ΔXΔP=/2\Delta X \Delta P = \hbar/2。提示:直接计算 X\langle X \rangleX2\langle X^2 \rangleP\langle P \rangleP2\langle P^2 \rangle

  2. 交换算符的本征值:证明对于两个全同费米子,若它们处于相同的单粒子态(即 ψ=ϕ1ϕ2|\psi\rangle = |\phi\rangle_1 \otimes |\phi\rangle_2),则 P12ψ=ψP_{12}|\psi\rangle = |\psi\rangle,但这与费米子的反对称要求矛盾,因此这样的态不存在。这就是Pauli不相容原理。

  3. 质心系变换:对于氢原子(质子质量 mpm_p,电子质量 mem_e),计算约化质量 μ\mu 并与电子质量比较。证明 μme(1me/mp)\mu \approx m_e (1 - m_e/m_p),并讨论为什么化学中通常直接用 mem_e 代替 μ\mu 计算能级时误差很小。

  4. 波包扩散的数值估计:一个电子波包初始宽度为 ΔX0=109\Delta X_0 = 10^{-9} m(纳米尺度),估算它扩散到宽度 2ΔX02\Delta X_0 需要多长时间。并与原子尺度(101010^{-10} m)的结果比较,讨论扩散时间对初始宽度的依赖关系。