"世界不是静止的。当外部世界随时间变化,量子系统被迫响应——吸收光子、发射光子、在不同能级间跃迁。"
时间依赖微扰理论告诉我们:量子跃迁是如何发生的。
第9章 时间依赖微扰理论:量子跃迁的引擎
9.1 故事:量子灯光师的调光台
24世纪,"量子灯光师"不是控制舞台灯光,而是操控原子与光的相互作用。首席灯光师老陈正在调试一台新型量子计算机的光学接口。
"问题很简单,"他对实习生说,"我们要让原子从基态 ∣a⟩ 跃迁到激发态 ∣b⟩,然后用这个态来存储量子比特。"
"直接用激光照射不就行了?"
"没那么简单。原子不会’立即’跳上去——它是一个渐进的积累过程。如果时间依赖微扰很小,原子在每个瞬间都只有很小的概率发生跃迁。只有在特定条件下——当激光频率匹配能级差时——跃迁概率才会随时间线性增长。"
"匹配能级差?"
"ω≈ω0=(Eb−Ea)/ℏ。这叫共振。就像推秋千:只有在正确的节奏推,能量才能有效传递。"
实习生若有所思地看着屏幕上闪烁的能级图。老陈继续说:"但你要知道,这个理论框架诞生于一百多年前——1917年,爱因斯坦在一篇关于辐射与物质平衡的文章中,第一次系统地提出了吸收、受激发射和自发辐射三种过程。没有那篇文章,就没有今天的激光器。"
核心洞察:量子跃迁不是瞬发事件,而是微扰与系统持续相互作用的结果。频率匹配时,跃迁概率随时间积累。
9.2 前置知识:傅里叶变换基础
时间依赖微扰理论大量涉及频域分析,因此我们需要回顾傅里叶变换的核心概念。
9.2.1 傅里叶变换的定义
对于时间函数 f(t),其傅里叶变换为:
f~(ω)=∫−∞∞f(t)eiωtdt
逆变换:
f(t)=2π1∫−∞∞f~(ω)e−iωtdω
物理意义:f(t) 可以分解为不同频率 ω 的简谐振动的叠加,f~(ω) 是频率 ω 的"振幅"。
9.2.2 关键性质
- 线性:F[af+bg]=af~+bg~
- 平移:F[f(t−t0)]=eiωt0f~(ω)
- 调制:F[f(t)e−iω0t]=f~(ω−ω0)
- 微分:F[dtdf]=−iωf~(ω)
9.2.3 delta函数的傅里叶表示
δ(ω−ω0)=2π1∫−∞∞ei(ω−ω0)tdt
这在含时微扰中至关重要——它揭示了能量守恒与时间演化的深刻联系。
9.2.4 sinc函数与不确定性
矩形脉冲的傅里叶变换是sinc函数:
F[rect(t/T)]=T⋅sinc(ωT/2)=T⋅ωT/2sin(ωT/2)
关键洞察:脉冲越窄(T 越小),频谱越宽——这正是能量-时间不确定性 ΔEΔt∼ℏ 的数学体现。
9.3 含时微扰理论框架
9.3.1 问题设置与基本方程
设总哈密顿量可以写成:
H(t)=H0+H′(t)
其中 H0 是未微扰哈密顿量(时间独立),我们已经完全求解了它的本征值问题:
H0∣n⟩=En∣n⟩
H′(t) 是含时微扰,满足 ∣H′∣≪∣H0∣(微扰条件)。
含时薛定谔方程为:
iℏ∂t∂∣Ψ(t)⟩=H(t)∣Ψ(t)⟩
我们的目标是:已知系统初始处于 H0 的某个本征态,求它在微扰作用下跃迁到其他态的概率。
graph TD
A[含时微扰理论] --> B["已知: H₀的本征态|n⟩"]
A --> C["加入: H'(t)"]
B --> D[系统初始态]
C --> E["求解: 跃迁概率"]
D --> F["|Ψ(0)⟩ = |i⟩"]
E --> G["P(i→f,t) = |c_f(t)|²"]
style A fill:#e3f2fd
style G fill:#e8f5e9 9.3.2 相互作用绘景的详细解释
为了把微扰效应从自由演化中分离出来,我们引入相互作用绘景(Interaction Picture)。
三种绘景对比:
| 绘景 | 态矢演化 | 算符演化 |
|---|
| 薛定谔绘景 | $i\hbar\frac{\partial}{\partial t} | \Psi\rangle = H |
| 海森堡绘景 | 态矢不演化 | iℏdtdA=[A,H] |
| 相互作用绘景 | 仅微扰驱动演化 | 仅 H0 驱动演化 |
定义相互作用绘景中的态矢:
∣ΨI(t)⟩=eiH0t/ℏ∣Ψ(t)⟩
它满足的方程是:
iℏdtd∣ΨI(t)⟩=HI′(t)∣ΨI(t)⟩
其中 HI′(t)=eiH0t/ℏH′(t)e−iH0t/ℏ 是相互作用绘景中的微扰。
物理直觉:相互作用绘景就像"站在未微扰系统的旋转参考系中"观察微扰效应。自由演化被"剔除"了,只剩下纯微扰驱动的变化。
将 ∣ΨI(t)⟩ 按未微扰本征态展开:
∣ΨI(t)⟩=n∑cn(t)∣n⟩
代入方程,得到展开系数的耦合微分方程组:
c˙n(t)=−ℏim∑⟨n∣HI′(t)∣m⟩cm(t)
=−ℏim∑⟨n∣H′(t)∣m⟩ei(En−Em)t/ℏcm(t)
定义 Bohr频率:
ωnm=ℏEn−Em
则方程变为:
\boxed{\dot{c}_n(t) = -\frac{i}{\hbar} \sum_m H'_{nm}(t) e^{i\omega_{nm}t} c_m(t)}
这是含时微扰理论的基本方程。它是一组无穷多个耦合的一阶常微分方程。一般情况下无法精确求解,需要近似方法。
9.4 二能级系统
9.4.1 两态耦合模型
二能级系统是理解量子跃迁最清晰的场景。设未微扰系统只有两个本征态:
H0∣a⟩=Ea∣a⟩,H0∣b⟩=Eb∣b⟩
定义能级间隔对应的频率:
ω0=ℏEb−Ea
系统的一般态可以展开为:
∣Ψ(t)⟩=ca(t)∣a⟩e−iEat/ℏ+cb(t)∣b⟩e−iEbt/ℏ
系数 ca(t) 和 cb(t) 满足耦合方程:
c˙a=−ℏi[caHaa′+cbHab′e−iω0t]
c˙b=−ℏi[cbHbb′+caHba′eiω0t]
graph TD
A[二能级系统] -->|未微扰| B["|a⟩: E_a"]
A -->|未微扰| C["|b⟩: E_b"]
B -->|"微扰 H'"| C
C -->|"微扰 H'"| B
style B fill:#e3f2fd
style C fill:#fff3e0 9.4.2 简并二能级系统的精确解
考虑简化情形:Haa′=Hbb′=0,Hab′=Hba′=V(常数,时间独立)。
从方程出发,对 c˙a 再求导:
c¨a+iω0c˙a+ℏ2V2ca=0
特征方程:
λ2+iω0λ+ℏ2V2=0
解得:
λ=2−iω0±√−ω02+4V2/ℏ2
Rabi频率:
Ω=√ω02+ℏ24V2
特别地,如果共振(ω0=0,简并):
Pa→b(t)=sin2(ℏVt)
物理图像:Rabi振荡——就像两个用弹簧连接的摆,能量在它们之间来回传递。
graph LR
A["系统初始在|a⟩"] -->|"t=πℏ/2V"| B["完全跃迁到|b⟩"]
B -->|"t=πℏ/V"| C["回到|a⟩"]
C -->|周期性| A
style B fill:#fff3e0 9.4.3 数值例子:Rabi振荡
问题:二能级系统,能级差 Eb−Ea=1 eV,微扰矩阵元 V=0.1 eV。计算Rabi频率和振荡周期。
解答:
步骤1:计算参数
ℏω0=1 eV⇒ω0=ℏ1=1.055×10−341.6×10−19≈1.52×1015 rad/s
V=0.1 eV=0.1×1.6×10−19 J=1.6×10−20 J
步骤2:计算Rabi频率
Ω=√ω02+ℏ24V2=√(1.52×1015)2+4×(1.52×1014)2
=1.52×1015×√1+0.04≈1.55×1015 rad/s
步骤3:振荡周期
TRabi=Ω2π≈1.55×10152π≈4.05×10−15 s=4.05 fs
物理意义:在飞秒时间尺度上,系统在两个态之间周期性振荡。这正是量子计算中量子比特操控的基本原理。
9.5 正弦微扰与一阶近似
9.5.1 一阶含时微扰理论
最常见的情况是微扰随时间正弦振荡(比如电磁波照射原子):
H′(t)=Vcos(ωt)=2V(eiωt+e−iωt)
假设 V 的对角元为零,非对角元 Vab=⟨a∣V∣b⟩。
一阶近似下:
cb(1)(t)=−ℏi∫0t⟨b∣H′(t′)∣a⟩eiω0t′dt′
代入 H′(t′)=Vcos(ωt′):
cb(1)(t)=−2ℏVba[ω0+ωei(ω0+ω)t−1+ω0−ωei(ω0−ω)t−1]
graph TD
A["正弦微扰 H'(t)=Vcos(ωt)"] --> B[两项贡献]
B --> C["反共振项
(ω₀+ω)⁻¹"]
B --> D["共振项
(ω₀-ω)⁻¹"]
C --> E["快速振荡
平均贡献小"]
D --> F["分母趋零
主导跃迁"]
style D fill:#fff3e0
style F fill:#e8f5e9 9.5.2 共振条件
当驱动频率 ω≈ω0 时,第二项的分母接近零——共振!
Pa→b(t)=∣cb(1)∣2≈(2ℏ∣Vba∣)2(ω0−ω)2sin2[(ω0−ω)t/2]
对于精确的共振 ω=ω0:
Pa→b=(2ℏ∣Vba∣)2t2
问题:概率随时间二次增长,最终会超过1,违反概率守恒!这说明一阶近似只在短时间内有效。
9.6 费米黄金定则的完整推导
9.6.1 从分立态到连续谱
在实际物理系统中,末态往往不是单一态,而是一个连续谱(比如原子电离到自由电子态,或大量原子在不同环境中)。定义末态的态密度 ρ(Eb)。
将跃迁概率对所有可能的末态积分:
Ptotal=∫−∞∞(2ℏ∣Vba∣)2(ω0−ω)2sin2[(ω0−ω)t/2]ρ(Eb)dEb
9.6.2 sinc函数的积分
当 t 足够大时,sin2(x)/x2 函数变得极其尖锐。利用:
∫−∞∞α2sin2(αt/2)dα=2πt
9.6.3 费米黄金定则
于是得到单位时间的跃迁概率(跃迁速率):
\boxed{R_{a \to b} = \frac{2\pi}{\hbar} |V_{ba}|^2 \rho(E_b)}
这就是费米黄金定则(Fermi’s Golden Rule)。
关键特征:
- 时间依赖:跃迁速率是常数
- 频率响应:只在能量守恒条件 Eb≈Ea+ℏω 附近发生
- 矩阵元平方:速率正比于微扰矩阵元的平方
graph LR
A["驱动频率 ω"] -->|"接近 ω₀"| B[共振跃迁]
A -->|"远离 ω₀"| C["抑制
快速振荡抵消"]
style B fill:#e8f5e9
style C fill:#ffebee
9.7 辐射的发射和吸收
9.7.1 电磁相互作用哈密顿量
原子与电磁场的相互作用来自最小耦合替换 p⃗→p⃗−qA⃗。对于带电粒子(电荷 −e),展开后得到:
H′=meA⃗⋅p⃗+2me2A⃗2
在弱场下忽略 A2 项,得到电偶极近似:
H′≈−d⃗⋅E⃗(t)
其中 d⃗=−er⃗ 是电偶极矩算符。
9.7.2 吸收过程
跃迁速率:
Ra→b=3ε0ℏ2π∣d⃗ba∣2ρ(ω0)
其中 ρ(ω0) 是电磁场的能量密度。
9.8 氢原子的跃迁:详细计算
9.8.1 选择定则的推导
电偶极跃迁要求矩阵元 d⃗ba=−e⟨b∣r⃗∣a⟩≠0。
氢原子波函数:ψnlm(r,θ,ϕ)=Rnl(r)Ylm(θ,ϕ)
位置算符 r⃗ 可以用球谐函数表示:
z=rcosθ∝rY10
x±iy=rsinθe±iϕ∝rY1±1
利用球谐函数的角向积分:
∫Yl′m′∗Y1qYlmdΩ
这个积分只有当满足以下条件时才不为零:
Δl=l′−l=±1
Δm=m′−m=0,±1
这就是选择定则。
9.8.2 数值例子:氢原子1s→2p跃迁
问题:计算氢原子从基态 ∣1s⟩ 到第一激发态 ∣2p⟩ 的电偶极矩阵元和跃迁速率。
解答:
步骤1:波函数
ψ1s=√πa031e−r/a0
ψ2p,m=√32πa051re−r/2a0Y1m(θ,ϕ)
步骤2:计算矩阵元
取 z 方向(Δm=0):
dz=−e⟨2p,0∣z∣1s⟩=−e∫ψ2p,0∗(rcosθ)ψ1sd3r
径向积分:
R=∫0∞r3R21(r)R10(r)dr=∫0∞r3⋅√24a05re−r/2a0⋅√a032e−r/a0dr
=√24a082∫0∞r4e−3r/2a0dr=√24a042⋅(3/2a0)54!=√242×24⋅24332a0
=243√24512a0≈0.74a0
角向积分:
∫Y10∗cosθY00dΩ=√34π∫Y10∗Y10√4π1dΩ=√31
因此:
dz=−e⋅0.74a0⋅√31≈−0.43ea0
步骤3:爱因斯坦A系数
A2p→1s=3πε0ℏc3ω03∣dba∣2
其中:
- ℏω0=E2−E1=10.2 eV =1.63×10−18 J
- ω0=2.46×1015 rad/s
- ∣dba∣2=(0.43ea0)2=0.185e2a02=0.185×(1.6×10−19)2×(5.29×10−11)2
=0.185×2.56×10−38×2.8×10−21=1.33×10−58 C2m2
代入:
A=3π×8.85×10−12×1.055×10−34×(3×108)3(2.46×1015)3×1.33×10−58
=3π×8.85×10−12×1.055×10−34×2.7×10251.49×1047×1.33×10−58
≈6.3×108 s−1
步骤4:激发态寿命
τ=A1≈6.3×1081≈1.6×10−9 s=1.6 ns
实验值:τ≈1.6 ns,与计算高度一致!
9.8.3 三种跃迁过程
graph TD
subgraph 光与物质的三种相互作用
A[吸收] -->|"原子从低能态 → 高能态"| B["光子消失
能量存入原子"]
C[受激发射] -->|"原子从高能态 → 低能态"| D["产生光子
与入射光子相干"]
E[自发辐射] -->|"原子从高能态 → 低能态"| F["随机产生光子
无外部驱动"]
end
style B fill:#e8f5e9
style D fill:#fff3e0
style F fill:#fce4ec
9.9 激光的工作原理
9.1 爱因斯坦的洞察
1917年,爱因斯坦用纯粹的热力学推理发现了自发辐射的存在。
在热平衡中:
NaBabρ(ω0)=Nb[Bbaρ(ω0)+Aba]
结合玻尔兹曼分布 Nb/Na=e−ℏω0/kBT,得到:
Bab=Bba,BbaAba=π2c3ℏω03
9.2 粒子数反转与光放大
激光需要粒子数反转:Nb>Na(非热平衡态)。
在粒子数反转介质中:
- 一个入射光子诱发受激发射
- 产生第二个同相位光子
- 链式反应,形成相干光放大
增益系数:
g=3ε0ℏcπ∣d⃗ba∣2ω0(Nb−Na)ρ(ω0)
9.3 激光的五个要素
- 受激发射产生相干光子
- 粒子数反转维持放大条件
- 谐振腔提供反馈
- 泵浦维持能量输入
- 输出耦合提取激光
graph LR
A[泵浦能量] -->|激发原子| B[粒子数反转]
B -->|一个光子进入| C[受激发射链式反应]
C -->|谐振腔反馈| D[相干激光输出]
style B fill:#fff3e0
style D fill:#e8f5e9 9.4 数值例子:氦氖激光器
问题:估算氦氖激光器中632.8 nm谱线的增益。
已知数据:
- 波长:λ=632.8 nm
- 频率:ν=c/λ=4.74×1014 Hz
- 跃迁:3s2→2p4(Ne原子)
- 线宽:Δν≈1.5 GHz(多普勒展宽)
- 上能级寿命:τ2≈10−7 s
- 下能级寿命:τ1≈10−8 s
解答:
步骤1:计算爱因斯坦A系数
A21=τ21≈107 s−1
步骤2:计算B系数
B21A21=c38πhν3
B21=8πhν3A21c3=8π×6.63×10−34×(4.74×1014)3107×(3×108)3
=8π×6.63×10−34×1.07×10442.7×1032≈1.5×1020 m3J−1s−2
步骤3:估算增益系数
假设粒子数密度差 ΔN=N2−N1≈1014 cm−3=1020 m−3:
g=chνB21ΔN⋅Δν1
=3×1086.63×10−34×4.74×1014×1.5×1020×1020×1.5×1091
≈1.0×10−6×1.5×1020×1020×6.7×10−10
≈0.01 m−1
物理意义:每传播1米,光强增加约1%。在1米长的谐振腔中,单程增益约1%,配合高反射率镜面(99%反射),可以维持稳定的激光振荡。
9.10 本章总结
graph TD
A[第9章核心] --> B[含时微扰理论]
A --> C[共振条件]
A --> D[光与物质相互作用]
A --> E[自发辐射]
B -->|一阶近似| F["跃迁概率 ∝ |V_ba|²"]
C -->|"ω ≈ ω₀"| G["能量守恒
ΔE = ℏω"]
D -->|三种过程| H["吸收、受激发射、自发辐射"]
E -->|量子电动力学| I[爱因斯坦A系数]
F --> J["选择定则
决定哪些跃迁允许"]
G --> K[光谱线位置]
H --> L[激光原理]
I --> M["原子寿命
能级宽度"]
style J fill:#e3f2fd
style K fill:#e8f5e9
style L fill:#fff3e0
style M fill:#fce4ec带走的三句话:
- 时间依赖微扰理论描述量子系统在时变外力下的响应:共振时跃迁概率最大。
- 光与物质的三种相互作用——吸收、受激发射、自发辐射——是激光和光谱学的物理基础。
- 自发辐射是量子电动力学的必然结果:即使没有外部电磁场,真空涨落也会"催促"激发态原子发射光子。
"自发辐射是量子电动力学最美丽的预言之一:真空不是空的,它充满了量子涨落,正是这些涨落’催促’激发态的原子发射光子。" —— Griffiths
9.11 练习与思考
共振宽度与能量-时间不确定性:二能级系统在正弦微扰下的跃迁概率作为频率的函数呈sinc²形状。证明这个分布的半高全宽(FWHM)约为 4π/t。这对应能量-时间不确定性关系 ΔEΔt∼ℏ。
选择定则与拉波特规则:证明对于电偶极跃迁,如果初态和末态的宇称相同,则 d⃗ba=0,跃迁被禁止。这就是拉波特选择定则(Laporte rule)。
氢原子寿命计算:氢原子 2p→1s 跃迁的自发辐射速率约为 A≈6×108 s−1。对应的激发态寿命 τ=1/A 是多少?根据能量-时间不确定性 ΔE≈ℏ/τ,计算谱线的自然线宽(用频率和波长表示)。
Rabi振荡:二能级系统,能级差 Eb−Ea=2 eV,微扰矩阵元 V=0.2 eV。计算Rabi频率和第一次完全跃迁到 ∣b⟩ 所需的时间。
费米黄金定则应用:电子在势阱中,基态能量 E1=1 eV,受到频率 ω=(E2−E1)/ℏ 的正弦微扰。若微扰矩阵元 ∣V21∣=0.01 eV,末态态密度 ρ(E2)=1020 J⁻¹,计算跃迁速率。