跃迁概率与斯塔克位移的量子力学计算
在原子物理与量子力学研究中，跃迁概率与斯塔克位移是两个核心概念，它们分别描述了原子在能级间相互作用的几率以及外电场对能级结构的扰动。深入理解并计算这两个参数，对于光谱分析、激光物理及等离子体诊断等领域至关重要。本文将系统阐述如何从量子力学基本原理出发，推导并计算这两项关键物理量。
费米黄金法则与跃迁概率计算
跃迁概率描述了原子从初始能态 $|i\rangle$ 跃迁至末态 $|f\rangle$ 的可能性。在微扰论框架下，这一过程通常由费米黄金法则（Fermi's Golden Rule）给出。该法则指出，单位时间内从初态跃迁到末态的概率 $W_{i \to f}$ 与微扰哈密顿量的矩阵元平方及末态态密度成正比。
对于电偶极跃迁，微扰哈密顿量 $\hat{H}'$ 通常取为 $-\hat{\mu} \cdot \mathbf{E}(t)$，其中 $\hat{\mu}$ 是电偶极矩算符，$\mathbf{E}(t)$ 是交变电场。若考虑单色光场，跃迁速率可表示为：
$$ W_{i \to f} = \frac{2\pi}{\hbar} |\langle f | \hat{H}' | i \rangle|^2 \rho(E_f) $$
在实际计算中，我们需要计算电偶极矩矩阵元 $\langle f | \hat{\mu} | i \rangle$。对于氢原子或类氢离子，利用解析解可以直接得出；而对于多电子原子，往往需要借助微扰论或数值方法。例如，在计算钠 D 线的共振跃迁概率时，需考虑自旋 - 轨道耦合对波函数的修正。若忽略自旋，主量子数 $n$、角量子数 $l$ 和磁量子数 $m$ 的选择定则为 $\Delta l = \pm 1$，$\Delta m = 0, \pm 1$。
此外，线型函数（Line Shape Function）也是计算总吸收或发射截面不可或缺的部分。在自然展宽占主导的情况下，线型通常为洛伦兹型，其半高全宽（FWHM）由自发辐射寿命 $\tau$ 决定，即 $\Gamma = \hbar / \tau$。通过积分线型函数，可以得到总的跃迁强度，进而推算出爱因斯坦 $A$ 系数。
斯塔克效应的量子力学处理
斯塔克效应是指原子能级在外加静电场作用下的分裂现象。对于非氢原子（如碱金属），由于存在量子亏损，其能级结构较为复杂，但在弱场近似下，仍可采用微扰论进行处理。
对于非简并能级，一级斯塔克位移为零，因为电偶极矩的期望值在空间反演对称性下为零。因此，必须考虑二级微扰，计算一级斯塔克位移 $\Delta E^{(1)}$：
$$ \Delta E^{(1)} = \sum_{n' \neq n} \frac{|\langle n' | \hat{\mu} \cdot \mathbf{E} | n \rangle|^2}{E_n - E_{n'}} $$
该公式表明，能级的移动取决于初态与所有其他可能态之间的电偶极耦合强度以及能级间距。对于氢原子，由于其能级简并度较高（$E_n$ 仅与主量子数 $n$ 有关），一级微扰论即可给出非零结果，导致能级线性分裂。分裂的大小与电场强度 $E$ 成正比，分裂间隔为 $3ea_0 n (k^2 - l^2/2)$，其中 $k$ 和 $l$ 为角动量量子数。
在强场情况下，微扰论失效，必须采用定态微扰论中的简并微扰法，或者使用数值对角化方法求解薛定谔方程。此时，新的本征态是原简并态的线性组合，称为斯塔克态。这些新态不再具有确定的 $m_l$ 量子数，但保留了总角动量在电场方向的分量特性。
数值计算实例与实验验证
为了直观理解上述理论，我们可以构建一个简单的数值计算示例。假设我们要计算氢原子 $n=2$ 能级在 $E=10^5$ V/m 电场下的斯塔克分裂。
首先，列出 $n=2$ 的简并态：$|2,0,0\rangle$, $|2,1,-1\rangle$, $|2,1,0\rangle$, $|2,1,1\rangle$。构建相应的哈密顿量矩阵，其中对角元为零（一级微扰），非对角元由电偶极矩矩阵元决定。例如，$\langle 2,0,0 | z | 2,1,-1 \rangle$ 不为零。
通过求解该 $4 \times 4$ 矩阵的本征值，可以得到分裂后的能级。计算结果显示，原本简并的能级分裂为三个能级，其中两个能量升高，一个能量降低，且分裂间距与电场强度成正比。这一结果与实验观测到的氢原子光谱在强电场下的分裂现象完全吻合。
综上所述，跃迁概率与斯塔克位移的计算不仅依赖于严格的量子力学推导，还需要结合具体的原子结构和实验条件进行数值求解。掌握这些计算方法，有助于我们更精准地解析复杂原子体系的光谱特征，推动相关领域的技术进步。