贡献者： addis
预备知识 量子散射的微分截面
，跃迁概率（一阶微扰）
本文使用原子单位制。
1. 束缚态之间的跃迁截面
对于频率为 $\omega$ 的平面电磁波，截面可以想象成与电磁波传播方向垂直放置的一块面积为 $\sigma(\omega)$ 的面元，使得原子从电磁波中吸收的平均功率恰好等于电磁波经过该面元的功率。
所以自然地，可以用以下思路定义束缚态之间的跃迁截面：想象有许多相同的原子，在两个特定束缚态之间的跃迁需要的能量为 $\omega_{ij}$，所以平均每个原子在某个跃迁中吸收的能量为
\begin{equation}
E = \omega_{ij}P_{j\to i}~,
\end{equation}
$P_{j\to i}$ 是电子从束缚态 $ \left\lvert j \right\rangle $ 跃迁到另一个束缚态 $ \left\lvert i \right\rangle $ 的概率。
另一方面，电磁波在单位面积上单位频率的能量为 $s(\omega)$，那么定义每个原子的截面定义为
\begin{equation}
E = \sigma_{j\to i} s(\omega_{ij})~,
\end{equation}
事实上，这里的 $\sigma_{j\to i}$ 并不只是截面，而是截面乘以一个吸收谱线的宽度 $\Delta \omega$，否则量纲对不上。但为了方便一般还是简称为截面。
从式 15 来看，$P_{j\to i}$ 同样和 $s(\omega_{ij})$ 成正比，所以从两式中消去后截面 $\sigma_{j\to i}$ 和入射的波包无关：
\begin{equation}
\sigma_{j\to i} = \frac{4\pi^2 q^2}{c m^2 \omega_{ij}} \left\lvert \hat{\boldsymbol{\mathbf{e}}} \boldsymbol\cdot \left\langle i \middle| \boldsymbol{\mathbf{p}} \middle| j \right\rangle \right\rvert ^2~.
\end{equation}
这和
[1] 的式 4.46 的定义一样，除了偶极子近似。
2. 束缚态到连续态的跃迁截面
对于束缚态到连续态的跃迁，单位频率下原子吸收的能量为
\begin{equation}
\frac{\mathrm{d}{E}}{\mathrm{d}{\omega}} = \sigma(\omega) s(\omega)~,
\end{equation}
如果再对立体角微分得某频率的微分截面
\begin{equation}
\frac{\partial^2 E}{\partial \omega \partial \Omega} = \frac{\partial \sigma}{\partial \Omega} s(\omega)~,
\end{equation}
使用速度规范，当波包经过以后，$ \boldsymbol{\mathbf{A}} = 0$。令初始束缚态 $ \left\lvert j \right\rangle $ 的能量为 $-I_0$，把 $\omega$ 看作 $k$ 的函数
\begin{equation}
\omega(k) = \frac{k^2}{2m} + I_0~,
\end{equation}
积分得总能量
\begin{equation}
E = \iint \frac{\partial \sigma}{\partial \Omega} s(\omega) \,\mathrm{d}{\omega} \,\mathrm{d}{\Omega} = \iint \frac{\partial \sigma}{\partial \Omega} s(\omega(k)) \frac{k}{m} \,\mathrm{d}{k} \,\mathrm{d}{\Omega} ~.
\end{equation}
另一方面，每个原子从 $ \left\lvert j \right\rangle $ 电离平均吸收的能量为
\begin{equation}
E = \iint \omega(k) P_{j\to { \boldsymbol{\mathbf{k}} }} \cdot k^2 \,\mathrm{d}{k} \,\mathrm{d}{\Omega} ~.
\end{equation}
对比以上两式可以把微分截面表示为 $ \boldsymbol{\mathbf{k}} $ 的函数，也可以看成是单位矢量 $ \hat{\boldsymbol{\mathbf{k}}} $ 和 $\omega$ 的函数
\begin{equation}
\frac{\partial \sigma}{\partial \Omega} = \frac{km \omega}{s(\omega)} P_{j\to \boldsymbol{\mathbf{k}} }~,
\end{equation}
长度规范下有（式 8 ）
\begin{equation}
\frac{\partial \sigma}{\partial \Omega} = \frac{4\pi^2 m\omega k q^2}{c} \left\lvert \left\langle \boldsymbol{\mathbf{k}} \middle| \hat{\boldsymbol{\mathbf{e}}} \boldsymbol\cdot \boldsymbol{\mathbf{r}} \middle| j \right\rangle \right\rvert ^2~.
\end{equation}
速度规范下有（
式 15 ）
\begin{equation}
\frac{\partial \sigma}{\partial \Omega} = \frac{4\pi^2 k q^2}{c m \omega} \left\lvert \left\langle \boldsymbol{\mathbf{k}} \middle| \hat{\boldsymbol{\mathbf{e}}} \boldsymbol\cdot \boldsymbol{\mathbf{p}} \middle| j \right\rangle \right\rvert ^2~.
\end{equation}
这与
[2] 的式 19.86 完全相同。
代码 1：hydrogen_ioniz_dif_cross_section_z.m
% hydrogen single ionization differential cross section
% use <Coulomb plane wave|z|n,l,m> as transition matrix
% eq_SIcros_2
function [ret, mel] = hydrogen_ioniz_dif_cross_section_z(Z,n,l,m,k,r_max)
% fine structure constant = 1/(speed of light) in a.u.
alpha = 0.007297352569278;
mel = zeros(size(k));
parfor i = 1:numel(k)
mel(i) = hydrogen_trans_dipole_plane(Z, k(i), 0, 0, n, l, m, r_max);
end
ret = (4*pi^2*(k.^2/2+0.5*Z^2/n^2).*k)*alpha .* abs(mel).^2;
end
未完成：单位能量的总截面是多少？
[1] ^ Bransden, Physics of Atoms and Molecules, 2ed
[2] ^ Eugen Merzbacher. Quantum Mechanics 3ed
致读者： 小时百科一直以来坚持所有内容免费无广告，这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此，我们请求广大读者
热心打赏 ，使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元，我们一周就能脱离亏损， 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款， 他们的付出帮助了 99% 的读者免费获取知识， 我们在此表示感谢。