本章介绍半经典光学理论，其中电磁场用经典方式描述，而物质系统用量子力学处理。这种方法成功解释了大多数激光-原子相互作用现象，包括吸收、受激发射、拉比振荡和饱和效应。我们将建立从经典渲染到量子光学的桥梁，为理解现代光学现象提供必要的理论基础。
完成本章后，您将能够：
在许多光学现象中，光场强度足够大，使得光子数不确定度相对较小。此时，可以将光场视为经典电磁波，而保持原子系统的量子描述。这种半经典近似在以下条件下有效：
| 相干态近似：激光场接近相干态 $ | \alpha\rangle$ |
定量分析：对于典型激光器（功率 P = 1 mW，波长 $\lambda = 632.8 \text{ nm}$），单模光子数：
$\langle n\rangle = P/(\hbar\omega\gamma) \approx 10^{12} \gg 1$
其中 $\gamma$ 是腔线宽。相对涨落：
$\Delta n/\langle n\rangle = 1/\sqrt{\langle n\rangle} \approx 10^{-6}$
因此半经典近似极其精确。
相干态的经典对应： 相干态 $|\alpha\rangle$ 是湮灭算符的本征态：$\hat{a}|\alpha\rangle = \alpha|\alpha\rangle$，其中 $\alpha = |\alpha|e^{i\varphi}$ 是复振幅。电场期望值：
| $\langle\alpha | \hat{E}(\mathbf{r},t) | \alpha\rangle = E_0\cos(\mathbf{k}\cdot\mathbf{r} - \omega t + \varphi)$ |
| 其中 $E_0 = 2 | \alpha | \sqrt{\hbar\omega/2\varepsilon_0V}$，V 是量子化体积。这正是经典电磁场的形式，证明了半经典处理的合理性。 |
半经典近似的适用边界：
从路径积分视角理解： 光场的路径积分在大光子数极限下由鞍点主导：
$Z = \int \mathcal{D}A \exp[iS[A]/\hbar] \approx \exp[iS[A_{cl}]/\hbar]$
其中 $A_{cl}$ 是经典场构型。量子修正项 $\sim O(\hbar/\langle n\rangle E)$，对于大光子数可忽略。
半经典近似的层次结构： 不同物理过程需要不同层次的量子描述：
与其他近似方法的关系：
| 绝热近似：时间缓变场，$ | \dot{\omega}/\omega^2 | \ll 1$ |
实验判据： 判断是否需要超越半经典的实验特征：
计算复杂度考虑：
考虑单色平面波：
$\mathbf{E}(\mathbf{r}, t) = \mathbf{E}_0 \cos(\mathbf{k}\cdot\mathbf{r} - \omega t + \varphi) = \frac{1}{2}\mathbf{E}_0[e^{i(\mathbf{k}\cdot\mathbf{r} - \omega t + \varphi)} + c.c.]$
其中：
在偶极近似下（原子尺度 $a_0 \ll$ 波长 $\lambda$），可以忽略空间变化：
$\mathbf{E}(t) = \mathbf{E}_0 \cos(\omega t - \varphi)$
偶极近似的有效性：对于可见光 $\lambda \approx 500 \text{ nm}$ 和 Bohr 半径 $a_0 \approx 0.05 \text{ nm}$：
$\mathbf{k}\cdot a_0 \approx 2\pi(a_0/\lambda) \approx 6 \times 10^{-4} \ll 1$
因此 $e^{i\mathbf{k}\cdot\mathbf{r}} \approx 1 + i\mathbf{k}\cdot\mathbf{r} \approx 1$ 在原子尺度上。
电磁场的完整数学描述： 从Maxwell方程组出发，在无源空间中：
$\nabla \times \mathbf{E} = -\partial\mathbf{B}/\partial t$ $\nabla \times \mathbf{B} = \mu_0\varepsilon_0\partial\mathbf{E}/\partial t$ $\nabla \cdot \mathbf{E} = 0$ $\nabla \cdot \mathbf{B} = 0$
波动方程：$\nabla^2\mathbf{E} - (1/c^2)\partial^2\mathbf{E}/\partial t^2 = 0$
平面波解的一般形式： $\mathbf{E}(\mathbf{r},t) = \text{Re}[\mathbf{E}_0 \exp(i\mathbf{k}\cdot\mathbf{r} - i\omega t)]$
其中色散关系：$\omega^2 = c^2k^2$
偏振态的完整描述： 任意偏振态可分解为两个正交分量：
$\mathbf{E} = E_x \hat{\mathbf{e}}_x + E_y \hat{\mathbf{e}}_y$
其中复振幅： $E_x = |E_x|\exp(i\varphi_x)$ $E_y = |E_y|\exp(i\varphi_y)$
偏振椭圆的参数：
| 椭圆度：$\tan(\chi) = \pm | E_y | / | E_x | $（当$\Delta\varphi = \pm\pi/2$） |
| 方位角：$\tan(2\psi) = 2 | E_x | E_y | \cos(\Delta\varphi)/( | E_x | ^2 - | E_y | ^2)$ |
特殊情况：
| 圆偏振：$ | E_x | = | E_y | $，$\Delta\varphi = \pm\pi/2$ |
超越偶极近似： 对于某些情况需要考虑高阶项：
$\mathbf{E}(\mathbf{r}, t) = \mathbf{E}(\mathbf{r}_0, t)[1 + i\mathbf{k}\cdot(\mathbf{r}-\mathbf{r}_0) - \frac{1}{2}(\mathbf{k}\cdot(\mathbf{r}-\mathbf{r}_0))^2 + …]$
光场的规范选择： 在Coulomb规范下（$\nabla\cdot\mathbf{A} = 0$），电磁场可表示为：
$\mathbf{E} = -\partial\mathbf{A}/\partial t$，$\mathbf{B} = \nabla \times \mathbf{A}$
矢势的量子化形式：
$\mathbf{A}(\mathbf{r},t) = \sum_{\mathbf{k},\lambda} \sqrt{\hbar/2\varepsilon_0\omega_{\mathbf{k}}V} \hat{\mathbf{e}}{\mathbf{k},\lambda} [\hat{a}{\mathbf{k},\lambda} e^{i(\mathbf{k}\cdot\mathbf{r}-\omega_{\mathbf{k}}t)} + h.c.]$
其中 $\hat{\mathbf{e}}{\mathbf{k},\lambda}$ 是偏振矢量，满足 $\mathbf{k}\cdot\hat{\mathbf{e}}{\mathbf{k},\lambda} = 0$。
脉冲光场的描述： 对于有限持续时间的激光脉冲：
$\mathbf{E}(\mathbf{r},t) = \mathbf{E}_0\varepsilon(t)\cos[\mathbf{k}\cdot\mathbf{r} - \int_0^t \omega(t’)dt’ + \varphi(t)]$
其中：
光强与Poynting矢量： 瞬时光强：$I(t) = \varepsilon_0c|\mathbf{E}(t)|^2$ 时间平均光强：$\langle I\rangle = \varepsilon_0c|\mathbf{E}_0|^2/2$ 能流密度：$\mathbf{S} = \mathbf{E} \times \mathbf{H} = \varepsilon_0c^2\mathbf{E} \times \mathbf{B}$
复杂光场的模式分解： 任意光场可分解为模式叠加：
$\mathbf{E}(\mathbf{r},t) = \sum_n c_n\mathbf{E}_n(\mathbf{r})e^{-i\omega_n t}$
常见模式基：
结构光场： 现代光学中的结构光场具有复杂的空间分布：
光场的相干性描述： 部分相干光需要用相关函数描述：
$\Gamma(\mathbf{r}_1,\mathbf{r}_2,\tau) = \langle E^*(\mathbf{r}_1,t)E(\mathbf{r}_2,t+\tau)\rangle$
相干度：$\gamma_{12}(\tau) = \Gamma(\mathbf{r}_1,\mathbf{r}_2,\tau)/\sqrt{[\Gamma(\mathbf{r}_1,\mathbf{r}_1,0)\Gamma(\mathbf{r}_2,\mathbf{r}_2,0)]}$
光场的统计性质： 不同光源的统计特性：
非线性光学效应的唯象描述： 强场下需考虑非线性极化：
$\mathbf{P} = \varepsilon_0[\chi^{(1)}\mathbf{E} + \chi^{(2)}\mathbf{E}\mathbf{E} + \chi^{(3)}\mathbf{E}\mathbf{E}\mathbf{E} + …]$
导致：
光场的时空耦合： 超短脉冲中时间和空间自由度耦合：
考虑具有离散能级的原子系统。系统状态由波函数描述：
| $ | \psi(t)\rangle = \sum_n c_n(t)e^{-iE_n t/\hbar} | n\rangle$ |
| 其中 $ | n\rangle$ 是能量本征态，$E_n$ 是对应能量。展开系数满足归一化： |
| $\sum_n | c_n(t) | ^2 = 1$ |
密度算符形式（纯态）：
| $\hat{\rho}(t) = | \psi(t)\rangle\langle\psi(t) | $ |
更一般的混合态：
| $\hat{\rho}(t) = \sum_i p_i | \psi_i(t)\rangle\langle\psi_i(t) | $ |
其中 $p_i$ 是系综权重，$\sum_i p_i = 1$。
密度矩阵元素：
| $\rho_{mn}(t) = \langle m | \hat{\rho}(t) | n\rangle = c_m(t)c_n^*(t)e^{-i(E_m-E_n)t/\hbar}$ |
物理意义：
| 对角元 $\rho_{nn}$：能级 $ | n\rangle$ 的布居数概率 |
量子态的几何表示： 在N维Hilbert空间中，纯态构成复射影空间$\text{CP}^{N-1}$。对于二能级系统：
| $ | \psi\rangle = \alpha | g\rangle + \beta | e\rangle$， $ | \alpha | ^2 + | \beta | ^2 = 1$ |
可参数化为： $\alpha = \cos(\theta/2)e^{i\varphi_1}$ $\beta = \sin(\theta/2)e^{i\varphi_2}$
整体相位不影响物理，取$\varphi_1 = 0$，得到Bloch球表示： $|\psi\rangle = \cos(\theta/2)|g\rangle + e^{i\varphi}\sin(\theta/2)|e\rangle$
密度矩阵的谱分解： 任意密度矩阵可对角化：
| $\hat{\rho} = \sum_i \lambda_i | \lambda_i\rangle\langle\lambda_i | $ |
| 其中$\lambda_i$是本征值（$0 \le \lambda_i \le 1$），$ | \lambda_i\rangle$是对应本征态。 |
纯度参数：$P = \text{Tr}(\hat{\rho}^2) = \sum_i \lambda_i^2$
量子态的纠缠测度： 对于复合系统$\hat{\rho}_{AB}$，部分迹给出约化密度矩阵：
$\hat{\rho}A = \text{Tr}_B(\hat{\rho}{AB})$
纠缠熵：$S(\hat{\rho}_A) = -\text{Tr}(\hat{\rho}_A \ln \hat{\rho}_A)$
密度矩阵的性质与约束：
| 正定性：$\langle\psi | \hat{\rho} | \psi\rangle \ge 0$ 对所有 $ | \psi\rangle$ |
von Neumann熵： 量子系统的熵定义为：
$S = -k_B \text{Tr}(\hat{\rho} \ln \hat{\rho}) = -k_B \sum_n \lambda_n \ln \lambda_n$
其中 $\lambda_n$ 是密度矩阵的本征值。
相干性的量化： $l_1$-范数相干性：$C_{l_1}(\rho) = \sum_{m\neq n} |\rho_{mn}|$ 相对熵相干性：$C_r(\rho) = S(\rho_{diag}) - S(\rho)$
其中 $\rho_{diag}$ 是去除非对角元后的密度矩阵。
开放系统动力学： 考虑系统与环境耦合，总密度矩阵演化：
$\hat{\rho}{total}(t) = \hat{U}(t)\hat{\rho}{total}(0)\hat{U}^\dagger(t)$
对环境求迹得约化密度矩阵：
$\hat{\rho}S(t) = \text{Tr}_E[\hat{\rho}{total}(t)]$
这导致非幺正演化和退相干。
主方程形式： 在马尔可夫近似下，系统演化遵循Lindblad主方程：
$d\hat{\rho}/dt = -i[\hat{H},\hat{\rho}]/\hbar + \sum_k \gamma_k(\hat{L}_k\hat{\rho}\hat{L}_k^\dagger - \frac{1}{2}{\hat{L}_k^\dagger\hat{L}_k,\hat{\rho}})$
其中 $\hat{L}_k$ 是Lindblad算符，描述不同的耗散通道。
原子能级结构的精细细节： 实际原子能级包含多重结构：
多电子原子的处理： 对于多电子系统，需要考虑：
中心场近似： $V_{eff}(\mathbf{r}) = V_{nuc}(\mathbf{r}) + V_{screen}(\mathbf{r})$
组态相互作用： $|\psi\rangle = \sum_i c_i|\Phi_i\rangle$（$\Phi_i$是Slater行列式）
LS耦合与jj耦合：
量子缺陷理论： 对于Rydberg态，能级可表示为：
$E_n = -R_\infty/(n - \delta_l)^2$
其中$\delta_l$是量子缺陷，反映核心电子的屏蔽效应。
时间依赖的量子系统： 对于含时哈密顿量$\hat{H}(t)$，演化算符：
$\hat{U}(t,t_0) = \mathcal{T} \exp[-i\int_{t_0}^t \hat{H}(t’)dt’/\hbar]$
其中$\mathcal{T}$是时间排序算符。
绝热定理与Berry相位： 缓慢变化的哈密顿量下，系统保持在瞬时本征态，但获得几何相位：
| $\gamma_n = i\oint\langle n(\mathbf{R}) | \nabla_{\mathbf{R}} | n(\mathbf{R})\rangle\cdot d\mathbf{R}$ |
这是Berry相位，纯几何起源，与动力学相位不同。
退相干机制： 主要退相干源包括：
在电偶极近似下，光-物质相互作用哈密顿量为：
$\hat{H}_{int} = -\hat{\mathbf{d}}\cdot\mathbf{E}(t)$
其中 $\hat{\mathbf{d}} = -e\hat{\mathbf{r}}$ 是电偶极矩算符（$e > 0$ 为基本电荷）。
总哈密顿量：
$\hat{H} = \hat{H}0 + \hat{H}{int} = \hat{H}_0 - \hat{\mathbf{d}}\cdot\mathbf{E}(t)$
偶极矩矩阵元素：
| $\mathbf{d}_{mn} = \langle m | \hat{\mathbf{d}} | n\rangle = -e\langle m | \hat{\mathbf{r}} | n\rangle$ |
选择定则：
相互作用能量尺度：
| $E_{int} \approx | \mathbf{d} | \mathbf{E} | \approx ea_0E_0$ |
对于 $E_0 = 10^6 \text{ V/m}$（典型激光场），$E_{int} \approx 10^{-23} \text{ J} \approx 10^{-4} \text{ eV}$。
精确的选择定则推导： 从宇称算符$\hat{P}$的性质出发：$\hat{P}|nlm\rangle = (-1)^l|nlm\rangle$
偶极算符具有奇宇称：$\hat{P}\hat{\mathbf{r}}\hat{P}^{-1} = -\hat{\mathbf{r}}$
因此跃迁矩阵元： $\langle n’l’m’|\hat{\mathbf{r}}|nlm\rangle \propto \langle l’m’|(-1)^{l’+l+1}|lm\rangle$
非零条件：$(-1)^{l’+l+1} = 1 \implies l’ + l = \text{奇数} \implies \Delta l = \pm 1, \pm 3, …$
但球谐函数的正交性进一步限制：$\Delta l = \pm 1$
相互作用能量的数量级估计： 不同尺度的比较：
强场判据：
偶极矩的对称性分析： 利用Wigner-Eckart定理，偶极矩矩阵元素可写为：
| $\langle n’l’m’ | \hat{d}_q | nlm\rangle = \langle n’l’ | \hat{d} | nl\rangle\langle l’m’ | lm;1q\rangle$ |
其中：
| $\langle n’l’ | \hat{d} | nl\rangle$ 是约化矩阵元素 |
| $\langle l’m’ | lm;1q\rangle$ 是Clebsch-Gordan系数 |
高阶多极矩展开： 完整的相互作用哈密顿量包含：
$\hat{H}_{int} = -\hat{\mathbf{d}}\cdot\mathbf{E} - \hat{\mathbf{Q}}:\nabla\mathbf{E} - \hat{\mathbf{m}}\cdot\mathbf{B} + …$
其中：
规范不变性： 在不同规范下，相互作用形式不同：
规范变换：$\hat{U} = \exp[ie \mathbf{A}(\mathbf{r},t)\cdot\mathbf{r}/\hbar]$
跃迁强度与振子强度： 振子强度定义：
| $f_{mn} = (2m/\hbar^2)\omega_{mn} | \langle m | \hat{\mathbf{r}} | n\rangle | ^2$ |
满足Thomas-Reich-Kuhn求和规则：
$\sum_n f_{n0} = N$（电子数）
相互作用图像： 在相互作用图像中，态矢和算符演化分离：
$|\psi_I(t)\rangle = e^{i\hat{H}_0 t/\hbar}|\psi_S(t)\rangle$ $\hat{O}_I(t) = e^{i\hat{H}_0 t/\hbar}\hat{O}_S e^{-i\hat{H}_0 t/\hbar}$
演化方程： $i\hbar\partial|\psi_I\rangle/\partial t = \hat{H}_{int,I}(t)|\psi_I\rangle$
| 在近共振条件下（$ | \omega - \omega_0 | \ll \omega_0$），可以采用旋转波近似（RWA），忽略快速振荡项。 |
详细推导：将电场写成正负频率分量：
$\mathbf{E}(t) = \mathbf{E}^+(t) + \mathbf{E}^-(t) = \frac{1}{2}\mathbf{E}_0e^{-i\omega t} + \frac{1}{2}\mathbf{E}_0^*e^{i\omega t}$
相互作用哈密顿量：
$\hat{H}{int} = -\hat{\mathbf{d}}\cdot\mathbf{E}(t) = -(\hat{d}+ + \hat{d}_-)\cdot(\mathbf{E}^+ + \mathbf{E}^-)$
| 其中 $\hat{d}_+ = | e\rangle\langle g | \mathbf{d}{eg}$，$\hat{d}- = | g\rangle\langle e | \mathbf{d}_{ge}$。 |
| 展开得四项，其中两项以频率 $\omega + \omega_0$ 快速振荡（反旋项），在时间尺度 $1/ | \omega - \omega_0 | $ 上平均为零。保留慢变项： |
$\hat{H}{int}^{RWA} = -\hbar\Omega/2(\hat{\sigma}+e^{-i\omega t} + \hat{\sigma}_-e^{i\omega t})$
其中：
| $\hat{\sigma}_+ = | e\rangle\langle g | $, $\hat{\sigma}_- = | g\rangle\langle e | $ 是升降算符 |
RWA有效条件：
| 近共振：$ | \omega - \omega_0 | \ll \omega_0$ |
| 考虑最简单的量子系统：二能级原子，基态 $ | g\rangle$ 和激发态 $ | e\rangle$： |
| 偶极矩：$\mathbf{d} = \langle e | \hat{\mathbf{d}} | g\rangle$ |
二能级模型的普适性：
完整哈密顿量（旋转坐标系中）：
$\hat{H} = \hbar\omega_0/2 \hat{\sigma}z - \hbar\Omega/2(\hat{\sigma}+e^{-i\omega t} + \hat{\sigma}_-e^{i\omega t})$
| 其中$\hat{\sigma}_z = | e\rangle\langle e | - | g\rangle\langle g | $ 是Pauli-z算符。 |
二能级系统的数学结构： Hilbert空间 $\mathcal{H} = \text{span}{|g\rangle, |e\rangle} \cong \mathbb{C}^2$
基矢的完备性和正交性： $|g\rangle\langle g| + |e\rangle\langle e| = \hat{I}$ $\langle g|e\rangle = \langle e|g\rangle = 0$ $\langle g|g\rangle = \langle e|e\rangle = 1$
任意算符的展开： $\hat{O} = O_{gg}|g\rangle\langle g| + O_{ee}|e\rangle\langle e| + O_{ge}|g\rangle\langle e| + O_{eg}|e\rangle\langle g|$
偶极矩的详细结构： 复偶极矩矢量：$\mathbf{d} = |\mathbf{d}|e^{i\varphi_d}\hat{\mathbf{e}}_d$
其中：
| $ | \mathbf{d} | = | \langle e | e\hat{\mathbf{r}} | g\rangle | $ 是偶极矩大小 |
实偶极矩的对称性： $\mathbf{d}{eg} = \langle e|\hat{\mathbf{d}}|g\rangle = \langle g|\hat{\mathbf{d}}|e\rangle^* = \mathbf{d}{ge}^*$
对角元为零（宇称禁戒）： $\langle g|\hat{\mathbf{d}}|g\rangle = \langle e|\hat{\mathbf{d}}|e\rangle = 0$
Pauli算符的完整代数： 定义Pauli算符：
| $\hat{\sigma}x = \hat{\sigma}+ + \hat{\sigma}_- = | e\rangle\langle g | + | g\rangle\langle e | $ |
| $\hat{\sigma}y = -i(\hat{\sigma}+ - \hat{\sigma}_-) = -i | e\rangle\langle g | + i | g\rangle\langle e | $ |
| $\hat{\sigma}_z = | e\rangle\langle e | - | g\rangle\langle g | $ |
满足对易关系： $[\hat{\sigma}i, \hat{\sigma}_j] = 2i\varepsilon{ijk}\hat{\sigma}_k$
反对易关系： ${\hat{\sigma}i, \hat{\sigma}_j} = 2\delta{ij}\hat{I}$
二能级系统的SU(2)对称性： 任意二能级态可参数化为：
| $ | \psi\rangle = \cos(\theta/2) | g\rangle + e^{i\varphi}\sin(\theta/2) | e\rangle$ |
这对应于Bloch球上的点 $(\theta, \varphi)$。幺正演化对应SO(3)旋转。
有效二能级系统的实现： 在多能级系统中，当满足以下条件时可约化为二能级：
| 大失谐条件：$ | \Delta_{ij} | \gg \Omega_{ij}$（其他跃迁失谐远大于拉比频率） |
绝热消除示例： 三能级系统 $|g\rangle, |e\rangle, |r\rangle$，当 $|\Delta_r| \gg \Omega_r$ 时：
| $\hat{H}_{eff} \approx \hbar\delta | e\rangle\langle e | - \hbar\Omega_{eff}/2( | e\rangle\langle g | + h.c.)$ |
| 其中 $\delta = \Delta + | \Omega_r | ^2/4\Delta_r$ 是AC Stark位移。 |
密度矩阵形式：
$\hat{\rho} = \begin{pmatrix} \rho_{ee} & \rho_{eg} \ \rho_{ge} & \rho_{gg} \end{pmatrix}$
满足约束：
物理意义的深化：
| 纯态：$\text{Tr}(\hat{\rho}^2) = 1$，如 $ | \psi\rangle = \alpha | g\rangle + \beta | e\rangle$ |
| 相干性：$ | \rho_{eg} | ^2 \le \rho_{ee}\rho_{gg}$（Cauchy-Schwarz不等式） |
可观测量的期望值：
$\langle\hat{O}\rangle = \text{Tr}(\hat{\rho}\hat{O}) = \sum_{ij} \rho_{ij}\hat{O}_{ji}$
从von Neumann方程出发：
$i\hbar\partial\hat{\rho}/\partial t = [\hat{H}, \hat{\rho}]$
详细推导过程：
计算对易子： $[\hat{H}, \hat{\rho}]{-e} = -\hbar\omega_0\rho{eg}/2 - \hbar\Omega e^{-i\omega t}(\rho_{gg} - \rho_{ee})/2$ $[\hat{H}, \hat{\rho}]{e-} = -\hbar\omega_0\rho{eg}/2 - \hbar\Omega e^{i\omega t}(\rho_{ee} - \rho_{gg})/2$ $[\hat{H}, \hat{\rho}]{–} = -[\hat{H}, \hat{\rho}]{ee} = -\hbar\Omega(\rho_{eg}e^{i\omega t} - \rho_{ge}e^{-i\omega t})/2$
引入旋转坐标系： $\tilde{\rho}{eg} = \rho{eg} e^{i\omega t}$ $\tilde{\rho}{ge} = \rho{ge} e^{-i\omega t}$
这消除了快速振荡的$e^{\pm i\omega_0 t}$因子。
加入弛豫项（唯象处理）：
最终得到光学Bloch方程：
$d\rho_{ee}/dt = i\Omega/2(\tilde{\rho}{eg} - \tilde{\rho}{ge}) - \Gamma\rho_{ee}$ $d\tilde{\rho}{eg}/dt = i\Delta\tilde{\rho}{eg} + i\Omega/2(\rho_{ee} - \rho_{gg}) - \gamma\tilde{\rho}_{eg}$
其中：
弛豫时间关系：
完整的矩阵形式推导： 将密度矩阵写成2×2形式：
$\hat{\rho} = \begin{pmatrix} \rho_{gg} & \rho_{ge} \ \rho_{eg} & \rho_{ee} \end{pmatrix}$
哈密顿量矩阵：
$\hat{H} = \hbar/2 \begin{pmatrix} -\omega_0 & -\Omega e^{-i\omega t} \ -\Omega e^{i\omega t} & \omega_0 \end{pmatrix}$
对易子$[\hat{H}, \hat{\rho}]$的矩阵元素：
$[\hat{H}, \hat{\rho}]{gg} = -\hbar\Omega/2(\rho{eg}e^{-i\omega t} - \rho_{ge}e^{i\omega t})$ $[\hat{H}, \hat{\rho}]{ee} = \hbar\Omega/2(\rho{eg}e^{-i\omega t} - \rho_{ge}e^{i\omega t})$ $[\hat{H}, \hat{\rho}]{ge} = \hbar\omega_0\rho{ge} + \hbar\Omega/2(\rho_{gg} - \rho_{ee})e^{i\omega t}$ $[\hat{H}, \hat{\rho}]{eg} = -\hbar\omega_0\rho{eg} - \hbar\Omega/2(\rho_{gg} - \rho_{ee})e^{-i\omega t}$
弛豫项的微观起源： 从系统-环境耦合的微观模型出发：
$\hat{H}{total} = \hat{H}_S + \hat{H}_E + \hat{H}{SE}$
其中$\hat{H}{SE} = \sum_k g_k(\hat{\sigma}+\hat{a}k + \hat{\sigma}-\hat{a}_k^\dagger)$描述系统与环境模式的耦合。
在Born-Markov近似下，导出主方程：
$d\hat{\rho}_S/dt = -i[\hat{H}_S, \hat{\rho}_S]/\hbar + \mathcal{L}[\hat{\rho}_S]$
其中Lindblad超算符：
$\mathcal{L}[\hat{\rho}] = \Gamma(\bar{n} + 1)(\hat{\sigma}-\hat{\rho}\hat{\sigma}+ - \frac{1}{2}{\hat{\sigma}+\hat{\sigma}-,\hat{\rho}})$ $+ \Gamma\bar{n}(\hat{\sigma}+\hat{\rho}\hat{\sigma}- - \frac{1}{2}{\hat{\sigma}-\hat{\sigma}+,\hat{\rho}})$
这里$\bar{n}$是热平均光子数，$\Gamma$是自发辐射率。
纯失相过程： 额外的失相机制可表示为：
$\mathcal{L}{deph}[\hat{\rho}] = \gamma\varphi/2(\hat{\sigma}_z\hat{\rho}\hat{\sigma}_z - \hat{\rho})$
| 导致：$d\rho_{eg}/dt | {deph} = -\gamma\varphi\rho_{eg}/2$ |
总的横向弛豫率：$\gamma = \Gamma/2 + \gamma_\varphi$
定义Bloch矢量分量：
$u = \tilde{\rho}{eg} + \tilde{\rho}{ge} = 2\text{Re}(\tilde{\rho}{eg})$ $v = i(\tilde{\rho}{eg} - \tilde{\rho}{ge}) = -2\text{Im}(\tilde{\rho}{eg})$ $w = \rho_{ee} - \rho_{gg}$
物理意义：
Bloch方程的矢量形式：
$du/dt = \Delta v - u/T_2$ $dv/dt = -\Delta u + \Omega w - v/T_2$ $dw/dt = -\Omega v - (w - w_0)/T_1$
其中 $w_0 = -1$ 是热平衡时的布居数差。
紧凑形式：
$d\mathbf{R}/dt = \mathbf{\Omega}_{eff} \times \mathbf{R} - \mathbf{\Gamma}\cdot(\mathbf{R} - \mathbf{R}_0)$
其中：
Bloch矢量 $\mathbf{R} = (u, v, w)$ 在单位球内运动：
| $ | \mathbf{R} | ^2 = u^2 + v^2 + w^2 \le 1$ |
| 纯态：$ | \mathbf{R} | = 1$（球面上） |
| 混合态：$ | \mathbf{R} | < 1$（球内部） |
特殊点与态的对应：
| 北极 (0,0,1)：激发态 $ | e\rangle$ |
| 南极 (0,0,-1)：基态 $ | g\rangle$ |
| 赤道 (cosφ,sinφ,0)：相干叠加态 $( | g\rangle+e^{i\varphi} | e\rangle)/\sqrt{2}$ |
动力学解释：
| 大失谐情况（$ | \Delta | \gg\Omega$）：绕z轴快速进动 |
几何计算示例： 任意纯态可写为： $|\psi\rangle = \cos(\theta/2)|g\rangle + e^{i\varphi}\sin(\theta/2)|e\rangle$
对应Bloch矢量： $\mathbf{R} = (\sin\theta\cos\varphi, \sin\theta\sin\varphi, \cos\theta)$
其中$\theta, \varphi$是球坐标角度。
在精确共振时，Bloch方程简化为：
$du/dt = -u/T_2$ $dv/dt = \Omega w - v/T_2$ $dw/dt = -\Omega v - (w - w_0)/T_1$
忽略弛豫（$T_1, T_2 \to \infty$），得到：
$d^2w/dt^2 + \Omega^2w = 0$
解为：
$w(t) = \cos(\Omega t)$ $v(t) = -\sin(\Omega t)$ $u(t) = 0$
这描述了布居数的拉比振荡。
考虑失谐但无弛豫的情况，定义广义拉比频率：
$\Omega’ = \sqrt{\Omega^2 + \Delta^2}$
布居数演化：
$w(t) = (\Delta^2/\Omega’^2) + (\Omega^2/\Omega’^2)\cos(\Omega’t)$
最大激发概率：
$P_{max} = \Omega^2/(\Omega^2 + \Delta^2)$
对于时变电场 $\mathbf{E}(t)$，定义脉冲面积：
$\theta = \int_{-\infty}^{\infty} \Omega(t)dt = (1/\hbar)\int_{-\infty}^{\infty} \mathbf{d}\cdot\mathbf{E}(t)dt$
特殊脉冲：
脉冲传播方程（McCall-Hahn）： 对于在二能级介质中传播的短脉冲：
$\partial\theta/\partial z + (1/v_g)\partial\theta/\partial t = \alpha_0\sin\theta$
其中：
| $\alpha_0 = 2\pi N | \mathbf{d} | ^2/(\hbar c)$ 是吸收系数 |
自感应透明（SIT）：
实验应用：
在旋转坐标系中，Bloch矢量绕有效磁场 $\mathbf{\Omega}_{eff} = (\Omega, 0, \Delta)$ 进动：
$d\mathbf{R}/dt = \mathbf{\Omega}_{eff} \times \mathbf{R} - \mathbf{\Gamma}\cdot\mathbf{R}$
其中 $\mathbf{\Gamma}$ 是弛豫张量。
几何理解：
| 进动角频率：$\Omega’ = | \mathbf{\Omega}_{eff} | = \sqrt{\Omega^2 + \Delta^2}$ |
章动频率测量： 通过测量偶极矩振荡的拍频，可以推断失谐：
| $\Delta\nu_{beat} = | \Omega’ - \omega_0 | = | \Delta | $ |
| 当场参数缓慢变化时（$ | d\Omega/dt | \ll \Omega^2$），系统绝热跟随瞬时本征态： |
$w_{ad}(t) \approx -\Delta(t)/\Omega’(t)$
这是绝热快速通道（ARP）和STIRAP技术的基础。
绝热条件（Landau-Zener）：
| 绝热参数：$Q = \Omega^2/( | d\Delta/dt | ) \gg 1$ |
非绝热跃迁概率：
| $P_{na} = \exp(-\pi\Omega^2/2 | d\Delta/dt | )$ |
STIRAP（受激拉曼绝热通道）： 三能级系统中的相干布居数转移：
应用领域：
在连续波激发下，设 $d/dt = 0$，得到稳态解：
$u_0 = -2\Omega\Delta T_2^2/(1 + \Delta^2 T_2^2 + \Omega^2 T_1 T_2)$ $v_0 = -2\Omega T_2/(1 + \Delta^2 T_2^2 + \Omega^2 T_1 T_2)$ $w_0 = -(1 + \Delta^2 T_2^2)/(1 + \Delta^2 T_2^2 + \Omega^2 T_1 T_2)$
定义饱和参数：
$s = \Omega^2 T_1 T_2/(1 + \Delta^2 T_2^2) = I/I_{sat}$
其中饱和强度：
| $I_{sat} = \hbar^2/(2 | \mathbf{d} | ^2 T_1 T_2) \times (1 + \Delta^2 T_2^2)$ |
稳态激发态布居数：
$\rho_{ee} = s/(2(1 + s))$
吸收线型：
$L(\omega) = (\gamma’/2\pi)/[(\omega - \omega_0)^2 + (\gamma’/2)^2]$
其中展宽的线宽：
$\gamma’ = \gamma\sqrt{1 + s} = \gamma\sqrt{1 + I/I_{sat}}$
这就是功率展宽效应。
均匀展宽：所有原子具有相同跃迁频率
非均匀展宽：原子跃迁频率分布
在非均匀展宽介质中，强激光选择性激发特定速度类原子，在吸收谱中”烧”出一个孔：
孔宽度：$\Delta\omega_{hole} \approx \gamma’$ 孔深度：$\propto s/(1 + s)$
半经典理论无法解释：
在半经典框架中，通过唯象地引入自发辐射率 $A_{21}$：
| $d\rho_{ee}/dt | {sp} = -A{21}\rho_{ee}$ |
Einstein A系数与偶极矩的关系：
| $A_{21} = \omega_0^3 | \mathbf{d} | ^2/(3\pi\varepsilon_0\hbar c^3)$ |
修正的Bloch方程：
$dw/dt = -\Omega v - (w + 1)/T_1 - A_{21}(w + 1)/2$
| 当失谐 $ | \Delta | \gg \gamma$ 或强场 $\Omega \gg \gamma$ 时，相干项快速衰减，可采用速率方程： |
$d\rho_{ee}/dt = R_{12}\rho_{gg} - R_{21}\rho_{ee} - A_{21}\rho_{ee}$ $d\rho_{gg}/dt = -R_{12}\rho_{gg} + R_{21}\rho_{ee} + A_{21}\rho_{ee}$
其中受激跃迁率：
$R_{12} = R_{21} = (\Omega^2/4)\cdot(\gamma/2\pi)/[(\omega - \omega_0)^2 + (\gamma/2)^2]$
将原子介质视为具有频率依赖吸收/发射的体积：
$dL(\omega)/ds = -\sigma_a(\omega)n[1 - \rho_{ee}(\omega)]L(\omega) + \sigma_e(\omega)n\rho_{ee}(\omega)L_e(\omega)$
其中：
这建立了微观原子物理与宏观渲染方程的联系。
完整量子理论预测的修正：
这些效应需要第27-28章的完整量子光学处理。
本章建立了半经典光-物质相互作用理论：
关键公式汇总：
| 饱和强度：$I_{sat} = \hbar^2/(2 | \mathbf{d} | ^2 T_1 T_2)$ |
21.1 推导二能级系统在共振$\pi$脉冲作用下的演化。证明初始处于基态的原子在脉冲后完全转移到激发态。
| 21.2 计算氢原子1S-2P跃迁的饱和强度。已知跃迁波长$\lambda = 121.6 \text{ nm}$，偶极矩$ | \mathbf{d} | = 2.5 \times 10^{-29} \text{ C}\cdot\text{m}$，$T_1 = T_2 = 1.6 \text{ ns}$。 |
| 21.3 证明在大失谐极限（$ | \Delta | \gg \Omega, \gamma$）下，二能级系统的AC Stark位移为 $\delta E = \hbar\Omega^2/(4\Delta)$。 |
21.4 分析双色激光场下的二能级系统动力学。考虑两个频率$\omega_1$和$\omega_2$的激光，推导有效拉比频率和共振条件。
21.5 推导包含自发辐射的光学Bloch方程的稳态荧光谱（Mollow谱）。说明为什么半经典理论无法完全解释三峰结构。
21.6 设计一个绝热快速通道（ARP）脉冲序列，实现99%的布居反转效率。给出脉冲形状和参数选择标准。
21.7 讨论如何将半经典光-物质相互作用理论应用于计算机图形学中的荧光材料渲染。考虑多能级系统和能量转移过程。
21.8 探讨超快激光脉冲（飞秒量级）与物质相互作用时，半经典理论的修正。这对渲染超快现象有何启示？
| 正确：检查条件 $\Omega, | \Delta | \ll \omega_0$ |
| 正确：验证 $ | \dot{\omega}/\omega^2 | \ll 1$ |
设计半经典光-物质相互作用系统时，确保：
继续到第22章：偏振光学基础 →