第9章 时间依赖微扰理论:量子跃迁的发生

📑 目录

"世界不是静止的。当外部世界随时间变化,量子系统被迫响应——吸收光子、发射光子、在不同能级间跃迁。"

时间依赖微扰理论告诉我们:量子跃迁是如何发生的。


第9章 时间依赖微扰理论:量子跃迁的引擎

9.1 故事:量子灯光师的调光台

24世纪,"量子灯光师"不是控制舞台灯光,而是操控原子与光的相互作用。首席灯光师老陈正在调试一台新型量子计算机的光学接口。

"问题很简单,"他对实习生说,"我们要让原子从基态 a|a\rangle 跃迁到激发态 b|b\rangle,然后用这个态来存储量子比特。"

"直接用激光照射不就行了?"

"没那么简单。原子不会’立即’跳上去——它是一个渐进的积累过程。如果时间依赖微扰很小,原子在每个瞬间都只有很小的概率发生跃迁。只有在特定条件下——当激光频率匹配能级差时——跃迁概率才会随时间线性增长。"

"匹配能级差?"

"ωω0=(EbEa)/\omega \approx \omega_0 = (E_b - E_a)/\hbar。这叫共振。就像推秋千:只有在正确的节奏推,能量才能有效传递。"

实习生若有所思地看着屏幕上闪烁的能级图。老陈继续说:"但你要知道,这个理论框架诞生于一百多年前——1917年,爱因斯坦在一篇关于辐射与物质平衡的文章中,第一次系统地提出了吸收、受激发射和自发辐射三种过程。没有那篇文章,就没有今天的激光器。"

核心洞察:量子跃迁不是瞬发事件,而是微扰与系统持续相互作用的结果。频率匹配时,跃迁概率随时间积累。


9.2 前置知识:傅里叶变换基础

时间依赖微扰理论大量涉及频域分析,因此我们需要回顾傅里叶变换的核心概念。

9.2.1 傅里叶变换的定义

对于时间函数 f(t)f(t),其傅里叶变换为:

f~(ω)=f(t)eiωtdt\tilde{f}(\omega) = \int_{-\infty}^{\infty} f(t) e^{i\omega t} dt

逆变换:

f(t)=12πf~(ω)eiωtdωf(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} \tilde{f}(\omega) e^{-i\omega t} d\omega

物理意义f(t)f(t) 可以分解为不同频率 ω\omega 的简谐振动的叠加,f~(ω)\tilde{f}(\omega) 是频率 ω\omega 的"振幅"。

9.2.2 关键性质

  1. 线性F[af+bg]=af~+bg~\mathcal{F}[af + bg] = a\tilde{f} + b\tilde{g}
  2. 平移F[f(tt0)]=eiωt0f~(ω)\mathcal{F}[f(t-t_0)] = e^{i\omega t_0}\tilde{f}(\omega)
  3. 调制F[f(t)eiω0t]=f~(ωω0)\mathcal{F}[f(t)e^{-i\omega_0 t}] = \tilde{f}(\omega - \omega_0)
  4. 微分F[dfdt]=iωf~(ω)\mathcal{F}[\frac{df}{dt}] = -i\omega\tilde{f}(\omega)

9.2.3 delta函数的傅里叶表示

δ(ωω0)=12πei(ωω0)tdt\delta(\omega - \omega_0) = \frac{1}{2\pi} \int_{-\infty}^{\infty} e^{i(\omega - \omega_0)t} dt

这在含时微扰中至关重要——它揭示了能量守恒时间演化的深刻联系。

9.2.4 sinc函数与不确定性

矩形脉冲的傅里叶变换是sinc函数:

F[rect(t/T)]=Tsinc(ωT/2)=Tsin(ωT/2)ωT/2\mathcal{F}[\text{rect}(t/T)] = T \cdot \text{sinc}(\omega T/2) = T \cdot \frac{\sin(\omega T/2)}{\omega T/2}

关键洞察:脉冲越窄(TT 越小),频谱越宽——这正是能量-时间不确定性 ΔEΔt\Delta E \Delta t \sim \hbar 的数学体现。


9.3 含时微扰理论框架

9.3.1 问题设置与基本方程

设总哈密顿量可以写成:

H(t)=H0+H(t)H(t) = H_0 + H'(t)

其中 H0H_0 是未微扰哈密顿量(时间独立),我们已经完全求解了它的本征值问题:

H0n=EnnH_0 |n\rangle = E_n |n\rangle

H(t)H'(t) 是含时微扰,满足 HH0|H'| \ll |H_0|(微扰条件)。

含时薛定谔方程为:

itΨ(t)=H(t)Ψ(t)i\hbar \frac{\partial}{\partial t} |\Psi(t)\rangle = H(t) |\Psi(t)\rangle

我们的目标是:已知系统初始处于 H0H_0 的某个本征态,求它在微扰作用下跃迁到其他态的概率。

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
海森堡绘景态矢不演化idAdt=[A,H]i\hbar\frac{dA}{dt} = [A, H]
相互作用绘景仅微扰驱动演化H0H_0 驱动演化

定义相互作用绘景中的态矢:

ΨI(t)=eiH0t/Ψ(t)|\Psi_I(t)\rangle = e^{iH_0 t/\hbar} |\Psi(t)\rangle

它满足的方程是:

iddtΨI(t)=HI(t)ΨI(t)i\hbar \frac{d}{dt} |\Psi_I(t)\rangle = H'_I(t) |\Psi_I(t)\rangle

其中 HI(t)=eiH0t/H(t)eiH0t/H'_I(t) = e^{iH_0 t/\hbar} H'(t) e^{-iH_0 t/\hbar} 是相互作用绘景中的微扰。

物理直觉:相互作用绘景就像"站在未微扰系统的旋转参考系中"观察微扰效应。自由演化被"剔除"了,只剩下纯微扰驱动的变化。

ΨI(t)|\Psi_I(t)\rangle 按未微扰本征态展开:

ΨI(t)=ncn(t)n|\Psi_I(t)\rangle = \sum_n c_n(t) |n\rangle

代入方程,得到展开系数的耦合微分方程组:

c˙n(t)=imnHI(t)mcm(t)\dot{c}_n(t) = -\frac{i}{\hbar} \sum_m \langle n | H'_I(t) | m \rangle c_m(t)

=imnH(t)mei(EnEm)t/cm(t)= -\frac{i}{\hbar} \sum_m \langle n | H'(t) | m \rangle e^{i(E_n-E_m)t/\hbar} c_m(t)

定义 Bohr频率

ωnm=EnEm\omega_{nm} = \frac{E_n - E_m}{\hbar}

则方程变为:

\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 两态耦合模型

二能级系统是理解量子跃迁最清晰的场景。设未微扰系统只有两个本征态:

H0a=Eaa,H0b=EbbH_0 |a\rangle = E_a |a\rangle, \quad H_0 |b\rangle = E_b |b\rangle

定义能级间隔对应的频率:

ω0=EbEa\omega_0 = \frac{E_b - E_a}{\hbar}

系统的一般态可以展开为:

Ψ(t)=ca(t)aeiEat/+cb(t)beiEbt/|\Psi(t)\rangle = c_a(t) |a\rangle e^{-iE_a t/\hbar} + c_b(t) |b\rangle e^{-iE_b t/\hbar}

系数 ca(t)c_a(t)cb(t)c_b(t) 满足耦合方程:

c˙a=i[caHaa+cbHabeiω0t]\dot{c}_a = -\frac{i}{\hbar} \left[ c_a H'_{aa} + c_b H'_{ab} e^{-i\omega_0 t} \right]

c˙b=i[cbHbb+caHbaeiω0t]\dot{c}_b = -\frac{i}{\hbar} \left[ c_b H'_{bb} + c_a H'_{ba} e^{i\omega_0 t} \right]

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=0H'_{aa} = H'_{bb} = 0Hab=Hba=VH'_{ab} = H'_{ba} = V(常数,时间独立)。

从方程出发,对 c˙a\dot{c}_a 再求导:

c¨a+iω0c˙a+V22ca=0\ddot{c}_a + i\omega_0 \dot{c}_a + \frac{V^2}{\hbar^2} c_a = 0

特征方程:

λ2+iω0λ+V22=0\lambda^2 + i\omega_0 \lambda + \frac{V^2}{\hbar^2} = 0

解得:

λ=iω0±ω02+4V2/22\lambda = \frac{-i\omega_0 \pm \sqrt{-\omega_0^2 + 4V^2/\hbar^2}}{2}

Rabi频率

Ω=ω02+4V22\Omega = \sqrt{\omega_0^2 + \frac{4V^2}{\hbar^2}}

特别地,如果共振(ω0=0\omega_0 = 0,简并)

Pab(t)=sin2(Vt)P_{a \to b}(t) = \sin^2\left(\frac{Vt}{\hbar}\right)

物理图像:Rabi振荡——就像两个用弹簧连接的摆,能量在它们之间来回传递。

graph LR
    A["系统初始在|a⟩"] -->|"t=πℏ/2V"| B["完全跃迁到|b⟩"]
    B -->|"t=πℏ/V"| C["回到|a⟩"]
    C -->|周期性| A
    
    style B fill:#fff3e0

9.4.3 数值例子:Rabi振荡

问题:二能级系统,能级差 EbEa=1E_b - E_a = 1 eV,微扰矩阵元 V=0.1V = 0.1 eV。计算Rabi频率和振荡周期。

解答

步骤1:计算参数

ω0=1 eVω0=1=1.6×10191.055×10341.52×1015 rad/s\hbar \omega_0 = 1 \text{ eV} \Rightarrow \omega_0 = \frac{1}{\hbar} = \frac{1.6 \times 10^{-19}}{1.055 \times 10^{-34}} \approx 1.52 \times 10^{15} \text{ rad/s}

V=0.1 eV=0.1×1.6×1019 J=1.6×1020 JV = 0.1 \text{ eV} = 0.1 \times 1.6 \times 10^{-19} \text{ J} = 1.6 \times 10^{-20} \text{ J}

步骤2:计算Rabi频率

Ω=ω02+4V22=(1.52×1015)2+4×(1.52×1014)2\Omega = \sqrt{\omega_0^2 + \frac{4V^2}{\hbar^2}} = \sqrt{(1.52 \times 10^{15})^2 + 4 \times (1.52 \times 10^{14})^2}

=1.52×1015×1+0.041.55×1015 rad/s= 1.52 \times 10^{15} \times \sqrt{1 + 0.04} \approx 1.55 \times 10^{15} \text{ rad/s}

步骤3:振荡周期

TRabi=2πΩ2π1.55×10154.05×1015 s=4.05 fsT_{\text{Rabi}} = \frac{2\pi}{\Omega} \approx \frac{2\pi}{1.55 \times 10^{15}} \approx 4.05 \times 10^{-15} \text{ s} = 4.05 \text{ fs}

物理意义:在飞秒时间尺度上,系统在两个态之间周期性振荡。这正是量子计算中量子比特操控的基本原理。


9.5 正弦微扰与一阶近似

9.5.1 一阶含时微扰理论

最常见的情况是微扰随时间正弦振荡(比如电磁波照射原子):

H(t)=Vcos(ωt)=V2(eiωt+eiωt)H'(t) = V \cos(\omega t) = \frac{V}{2}\left(e^{i\omega t} + e^{-i\omega t}\right)

假设 VV 的对角元为零,非对角元 Vab=aVbV_{ab} = \langle a | V | b \rangle

一阶近似下:

cb(1)(t)=i0tbH(t)aeiω0tdtc_b^{(1)}(t) = -\frac{i}{\hbar} \int_0^t \langle b | H'(t') | a \rangle e^{i\omega_0 t'} dt'

代入 H(t)=Vcos(ωt)H'(t') = V \cos(\omega t')

cb(1)(t)=Vba2[ei(ω0+ω)t1ω0+ω+ei(ω0ω)t1ω0ω]c_b^{(1)}(t) = -\frac{V_{ba}}{2\hbar} \left[ \frac{e^{i(\omega_0 + \omega)t} - 1}{\omega_0 + \omega} + \frac{e^{i(\omega_0 - \omega)t} - 1}{\omega_0 - \omega} \right]

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\omega \approx \omega_0 时,第二项的分母接近零——共振

Pab(t)=cb(1)2(Vba2)2sin2[(ω0ω)t/2](ω0ω)2P_{a \to b}(t) = |c_b^{(1)}|^2 \approx \left(\frac{|V_{ba}|}{2\hbar}\right)^2 \frac{\sin^2[(\omega_0 - \omega)t/2]}{(\omega_0 - \omega)^2}

对于精确的共振 ω=ω0\omega = \omega_0

Pab=(Vba2)2t2P_{a \to b} = \left(\frac{|V_{ba}|}{2\hbar}\right)^2 t^2

问题:概率随时间二次增长,最终会超过1,违反概率守恒!这说明一阶近似只在短时间内有效。


9.6 费米黄金定则的完整推导

9.6.1 从分立态到连续谱

在实际物理系统中,末态往往不是单一态,而是一个连续谱(比如原子电离到自由电子态,或大量原子在不同环境中)。定义末态的态密度 ρ(Eb)\rho(E_b)

将跃迁概率对所有可能的末态积分:

Ptotal=(Vba2)2sin2[(ω0ω)t/2](ω0ω)2ρ(Eb)dEbP_{\text{total}} = \int_{-\infty}^{\infty} \left(\frac{|V_{ba}|}{2\hbar}\right)^2 \frac{\sin^2[(\omega_0 - \omega)t/2]}{(\omega_0 - \omega)^2} \rho(E_b) dE_b

9.6.2 sinc函数的积分

tt 足够大时,sin2(x)/x2\sin^2(x)/x^2 函数变得极其尖锐。利用:

sin2(αt/2)α2dα=πt2\int_{-\infty}^{\infty} \frac{\sin^2(\alpha t/2)}{\alpha^2} d\alpha = \frac{\pi t}{2}

9.6.3 费米黄金定则

于是得到单位时间的跃迁概率(跃迁速率):

\boxed{R_{a \to b} = \frac{2\pi}{\hbar} |V_{ba}|^2 \rho(E_b)}

这就是费米黄金定则(Fermi’s Golden Rule)。

关键特征

  1. 时间依赖:跃迁速率是常数
  2. 频率响应:只在能量守恒条件 EbEa+ωE_b \approx E_a + \hbar\omega 附近发生
  3. 矩阵元平方:速率正比于微扰矩阵元的平方
graph LR
    A["驱动频率 ω"] -->|"接近 ω₀"| B[共振跃迁]
    A -->|"远离 ω₀"| C["抑制
快速振荡抵消"] style B fill:#e8f5e9 style C fill:#ffebee

9.7 辐射的发射和吸收

9.7.1 电磁相互作用哈密顿量

原子与电磁场的相互作用来自最小耦合替换 ppqA\vec{p} \to \vec{p} - q\vec{A}。对于带电粒子(电荷 e-e),展开后得到:

H=emAp+e22mA2H' = \frac{e}{m}\vec{A} \cdot \vec{p} + \frac{e^2}{2m}\vec{A}^2

在弱场下忽略 A2A^2 项,得到电偶极近似

HdE(t)H' \approx -\vec{d} \cdot \vec{E}(t)

其中 d=er\vec{d} = -e\vec{r}电偶极矩算符

9.7.2 吸收过程

跃迁速率:

Rab=π3ε02dba2ρ(ω0)R_{a \to b} = \frac{\pi}{3\varepsilon_0\hbar^2} |\vec{d}_{ba}|^2 \rho(\omega_0)

其中 ρ(ω0)\rho(\omega_0) 是电磁场的能量密度。


9.8 氢原子的跃迁:详细计算

9.8.1 选择定则的推导

电偶极跃迁要求矩阵元 dba=ebra0\vec{d}_{ba} = -e\langle b | \vec{r} | a \rangle \neq 0

氢原子波函数:ψnlm(r,θ,ϕ)=Rnl(r)Ylm(θ,ϕ)\psi_{nlm}(r,\theta,\phi) = R_{nl}(r) Y_l^m(\theta,\phi)

位置算符 r\vec{r} 可以用球谐函数表示:

z=rcosθrY10z = r\cos\theta \propto r Y_1^0

x±iy=rsinθe±iϕrY1±1x \pm iy = r\sin\theta e^{\pm i\phi} \propto r Y_1^{\pm 1}

利用球谐函数的角向积分:

YlmY1qYlmdΩ\int Y_{l'}^{m'*} Y_1^{q} Y_l^m d\Omega

这个积分只有当满足以下条件时才不为零:

Δl=ll=±1\Delta l = l' - l = \pm 1

Δm=mm=0,±1\Delta m = m' - m = 0, \pm 1

这就是选择定则

9.8.2 数值例子:氢原子1s→2p跃迁

问题:计算氢原子从基态 1s|1s\rangle 到第一激发态 2p|2p\rangle 的电偶极矩阵元和跃迁速率。

解答

步骤1:波函数

ψ1s=1πa03er/a0\psi_{1s} = \frac{1}{\sqrt{\pi a_0^3}} e^{-r/a_0}

ψ2p,m=132πa05rer/2a0Y1m(θ,ϕ)\psi_{2p,m} = \frac{1}{\sqrt{32\pi a_0^5}} r e^{-r/2a_0} Y_1^m(\theta,\phi)

步骤2:计算矩阵元

zz 方向(Δm=0\Delta m = 0):

dz=e2p,0z1s=eψ2p,0(rcosθ)ψ1sd3rd_z = -e\langle 2p,0 | z | 1s \rangle = -e\int \psi_{2p,0}^* (r\cos\theta) \psi_{1s} d^3r

径向积分:

R=0r3R21(r)R10(r)dr=0r3r24a05er/2a02a03er/a0drR = \int_0^{\infty} r^3 R_{21}(r) R_{10}(r) dr = \int_0^{\infty} r^3 \cdot \frac{r}{\sqrt{24a_0^5}} e^{-r/2a_0} \cdot \frac{2}{\sqrt{a_0^3}} e^{-r/a_0} dr

=224a080r4e3r/2a0dr=224a044!(3/2a0)5=2×242432243a0= \frac{2}{\sqrt{24a_0^8}} \int_0^{\infty} r^4 e^{-3r/2a_0} dr = \frac{2}{\sqrt{24}a_0^4} \cdot \frac{4!}{(3/2a_0)^5} = \frac{2 \times 24}{\sqrt{24}} \cdot \frac{32}{243} a_0

=51224324a00.74a0= \frac{512}{243\sqrt{24}} a_0 \approx 0.74 a_0

角向积分:

Y10cosθY00dΩ=4π3Y10Y1014πdΩ=13\int Y_1^{0*} \cos\theta Y_0^0 d\Omega = \sqrt{\frac{4\pi}{3}} \int Y_1^{0*} Y_1^0 \frac{1}{\sqrt{4\pi}} d\Omega = \frac{1}{\sqrt{3}}

因此:

dz=e0.74a0130.43ea0d_z = -e \cdot 0.74 a_0 \cdot \frac{1}{\sqrt{3}} \approx -0.43 e a_0

步骤3:爱因斯坦A系数

A2p1s=ω03dba23πε0c3A_{2p \to 1s} = \frac{\omega_0^3 |d_{ba}|^2}{3\pi\varepsilon_0\hbar c^3}

其中:

  • ω0=E2E1=10.2\hbar\omega_0 = E_2 - E_1 = 10.2 eV =1.63×1018= 1.63 \times 10^{-18} J
  • ω0=2.46×1015\omega_0 = 2.46 \times 10^{15} rad/s
  • dba2=(0.43ea0)2=0.185e2a02=0.185×(1.6×1019)2×(5.29×1011)2|d_{ba}|^2 = (0.43 e a_0)^2 = 0.185 e^2 a_0^2 = 0.185 \times (1.6 \times 10^{-19})^2 \times (5.29 \times 10^{-11})^2

=0.185×2.56×1038×2.8×1021=1.33×1058 C2m2= 0.185 \times 2.56 \times 10^{-38} \times 2.8 \times 10^{-21} = 1.33 \times 10^{-58} \text{ C}^2\text{m}^2

代入:

A=(2.46×1015)3×1.33×10583π×8.85×1012×1.055×1034×(3×108)3A = \frac{(2.46 \times 10^{15})^3 \times 1.33 \times 10^{-58}}{3\pi \times 8.85 \times 10^{-12} \times 1.055 \times 10^{-34} \times (3 \times 10^8)^3}

=1.49×1047×1.33×10583π×8.85×1012×1.055×1034×2.7×1025= \frac{1.49 \times 10^{47} \times 1.33 \times 10^{-58}}{3\pi \times 8.85 \times 10^{-12} \times 1.055 \times 10^{-34} \times 2.7 \times 10^{25}}

6.3×108 s1\approx 6.3 \times 10^8 \text{ s}^{-1}

步骤4:激发态寿命

τ=1A16.3×1081.6×109 s=1.6 ns\tau = \frac{1}{A} \approx \frac{1}{6.3 \times 10^8} \approx 1.6 \times 10^{-9} \text{ s} = 1.6 \text{ ns}

实验值τ1.6\tau \approx 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]N_a B_{ab} \rho(\omega_0) = N_b \left[ B_{ba} \rho(\omega_0) + A_{ba} \right]

结合玻尔兹曼分布 Nb/Na=eω0/kBTN_b/N_a = e^{-\hbar\omega_0/k_B T},得到:

Bab=Bba,AbaBba=ω03π2c3B_{ab} = B_{ba}, \quad \frac{A_{ba}}{B_{ba}} = \frac{\hbar\omega_0^3}{\pi^2 c^3}

9.2 粒子数反转与光放大

激光需要粒子数反转Nb>NaN_b > N_a(非热平衡态)。

在粒子数反转介质中:

  • 一个入射光子诱发受激发射
  • 产生第二个同相位光子
  • 链式反应,形成相干光放大

增益系数

g=π3ε0cdba2ω0(NbNa)ρ(ω0)g = \frac{\pi}{3\varepsilon_0\hbar c} |\vec{d}_{ba}|^2 \omega_0 (N_b - N_a) \rho(\omega_0)

9.3 激光的五个要素

  1. 受激发射产生相干光子
  2. 粒子数反转维持放大条件
  3. 谐振腔提供反馈
  4. 泵浦维持能量输入
  5. 输出耦合提取激光
graph LR
    A[泵浦能量] -->|激发原子| B[粒子数反转]
    B -->|一个光子进入| C[受激发射链式反应]
    C -->|谐振腔反馈| D[相干激光输出]
    
    style B fill:#fff3e0
    style D fill:#e8f5e9

9.4 数值例子:氦氖激光器

问题:估算氦氖激光器中632.8 nm谱线的增益。

已知数据

  • 波长:λ=632.8\lambda = 632.8 nm
  • 频率:ν=c/λ=4.74×1014\nu = c/\lambda = 4.74 \times 10^{14} Hz
  • 跃迁:3s22p43s_2 \to 2p_4(Ne原子)
  • 线宽:Δν1.5\Delta\nu \approx 1.5 GHz(多普勒展宽)
  • 上能级寿命:τ2107\tau_2 \approx 10^{-7} s
  • 下能级寿命:τ1108\tau_1 \approx 10^{-8} s

解答

步骤1:计算爱因斯坦A系数

A21=1τ2107 s1A_{21} = \frac{1}{\tau_2} \approx 10^7 \text{ s}^{-1}

步骤2:计算B系数

A21B21=8πhν3c3\frac{A_{21}}{B_{21}} = \frac{8\pi h\nu^3}{c^3}

B21=A21c38πhν3=107×(3×108)38π×6.63×1034×(4.74×1014)3B_{21} = \frac{A_{21} c^3}{8\pi h\nu^3} = \frac{10^7 \times (3 \times 10^8)^3}{8\pi \times 6.63 \times 10^{-34} \times (4.74 \times 10^{14})^3}

=2.7×10328π×6.63×1034×1.07×10441.5×1020 m3J1s2= \frac{2.7 \times 10^{32}}{8\pi \times 6.63 \times 10^{-34} \times 1.07 \times 10^{44}} \approx 1.5 \times 10^{20} \text{ m}^3\text{J}^{-1}\text{s}^{-2}

步骤3:估算增益系数

假设粒子数密度差 ΔN=N2N11014\Delta N = N_2 - N_1 \approx 10^{14} cm3=1020^{-3} = 10^{20} m3^{-3}

g=hνcB21ΔN1Δνg = \frac{h\nu}{c} B_{21} \Delta N \cdot \frac{1}{\Delta\nu}

=6.63×1034×4.74×10143×108×1.5×1020×1020×11.5×109= \frac{6.63 \times 10^{-34} \times 4.74 \times 10^{14}}{3 \times 10^8} \times 1.5 \times 10^{20} \times 10^{20} \times \frac{1}{1.5 \times 10^9}

1.0×106×1.5×1020×1020×6.7×1010\approx 1.0 \times 10^{-6} \times 1.5 \times 10^{20} \times 10^{20} \times 6.7 \times 10^{-10}

0.01 m1\approx 0.01 \text{ 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

带走的三句话:

  1. 时间依赖微扰理论描述量子系统在时变外力下的响应:共振时跃迁概率最大。
  2. 光与物质的三种相互作用——吸收、受激发射、自发辐射——是激光和光谱学的物理基础。
  3. 自发辐射是量子电动力学的必然结果:即使没有外部电磁场,真空涨落也会"催促"激发态原子发射光子。

"自发辐射是量子电动力学最美丽的预言之一:真空不是空的,它充满了量子涨落,正是这些涨落’催促’激发态的原子发射光子。" —— Griffiths


9.11 练习与思考

  1. 共振宽度与能量-时间不确定性:二能级系统在正弦微扰下的跃迁概率作为频率的函数呈sinc²形状。证明这个分布的半高全宽(FWHM)约为 4π/t4\pi/t。这对应能量-时间不确定性关系 ΔEΔt\Delta E \Delta t \sim \hbar

  2. 选择定则与拉波特规则:证明对于电偶极跃迁,如果初态和末态的宇称相同,则 dba=0\vec{d}_{ba} = 0,跃迁被禁止。这就是拉波特选择定则(Laporte rule)。

  3. 氢原子寿命计算:氢原子 2p1s2p \to 1s 跃迁的自发辐射速率约为 A6×108A \approx 6 \times 10^8 s1^{-1}。对应的激发态寿命 τ=1/A\tau = 1/A 是多少?根据能量-时间不确定性 ΔE/τ\Delta E \approx \hbar/\tau,计算谱线的自然线宽(用频率和波长表示)。

  4. Rabi振荡:二能级系统,能级差 EbEa=2E_b - E_a = 2 eV,微扰矩阵元 V=0.2V = 0.2 eV。计算Rabi频率和第一次完全跃迁到 b|b\rangle 所需的时间。

  5. 费米黄金定则应用:电子在势阱中,基态能量 E1=1E_1 = 1 eV,受到频率 ω=(E2E1)/\omega = (E_2 - E_1)/\hbar 的正弦微扰。若微扰矩阵元 V21=0.01|V_{21}| = 0.01 eV,末态态密度 ρ(E2)=1020\rho(E_2) = 10^{20} J⁻¹,计算跃迁速率。