氢原子的束缚态波函数

                     

贡献者: addis

预备知识 球坐标系中的定态薛定谔方程,原子单位制

   本文使用原子单位制类氢原子(hydrogen-like atom)被定义为原子核有 $Z$ 个质子(核电荷为 $+Ze$)但只有一个核外电子的原子或离子,例如氢原子($Z=1$)、失去一个电子的氦原子 $\mathrm{He}^+$($Z=2$)和失去两个电子的锂离子 $\mathrm{Li}^{++}$($Z=3$)。

图
图 1:氢原子波函数的概率密度 $ \left\lvert \psi_{n,l,m}( \boldsymbol{\mathbf{r}} ) \right\rvert ^2$ 的 $x$-$z$ 截面,相对大小成比例。每个图中的三个数字分别是量子数 $n, l, m$,电子出现在白色圆圈内部的概率为 0.95。色条对应的数值是线性的,每个子图中的色条取值范围不相同。源码见 “氢原子波函数画图程序(Matlab)”。

   类氢原子的定态薛定谔方程为(式 1 ,原子单位制库仑势能 $V(r)=-Z/r$,电子质量 $m=1$)

\begin{equation} -\frac{1}{2m} \boldsymbol{\nabla}^2 \psi( \boldsymbol{\mathbf{r}} ) - \frac{Z}{r} \psi( \boldsymbol{\mathbf{r}} ) = E \psi( \boldsymbol{\mathbf{r}} )~, \end{equation}
类氢原子是唯一存在解析解的原子(离子)。该方程假设原子核不动,若要考虑原子核的运动,应该类似 “玻尔原子模型(约化质量)” 在式 1 使用约化质量 $\mu$ 代替电子质量 $m$。

   我们这里只讨论束缚态,即 $E < 0$ 的可归一化的解。1由于原子核的库仑势能具有球对称性,球坐标下的波函数具有最简单的形式。用分离变量法式 1 ,得波函数的表达式为

\begin{equation} \psi_{n,l,m}( \boldsymbol{\mathbf{r}} ) = \psi_{n,l,m}(r,\theta ,\phi) = R_{n,l}(r) Y_{l,m}(\theta, \phi)~. \end{equation}
其中 $n$ 是主量子数($n = 1, 2, \dots$),$l$ 是角量子数($l = 0, 1, \dots, n - 1$),$m$ 是磁量子数($m = -l, -l+1, \dots, l$)。$R_{n,l}(r)$ 是归一化的径向波函数(见下文),$Y_{l,m}(\theta, \phi)$ 是归一化的球谐函数(另见球谐函数表)。波函数概率密度的截面图如图 1 2

   波函数 $\psi_{n,l,m}$ 是能量本征态,对应的能级(本征值)为

\begin{equation} E_n = - \frac{m e^4}{32 \pi^2 \epsilon_0^2 \hbar ^2} \frac{Z^2}{n^2} \approx - 13.6 \,\mathrm{eV} \frac{Z^2}{n^2} \qquad (n = 1, 2, \dots)~, \end{equation}
考虑原子核运动时同样要把 $m$ 替换成 $\mu$。注意能级之和主量子数 $n$ 有关,多个不同的波函数具有相同能级的情况叫做简并,所以易得第 $n$ 能级有 $n^2$ 重简并。

   容易发现玻尔原子模型的能级和该式恰好相同,且都与氢原子光谱吻合,这也是当初人们接受玻尔模型的重要原因。但使用薛定谔方程才是真正符合量子力学的方法,因为波尔模型无法解释许多其他问题如多电子原子的能级。

1. 氢原子的径向波函数

   现在来解式 2 中的径向波函数 $R_{n,l}(r)$。球坐标系的薛定谔方程在分离变量后,可得(约化)径向方程式 3 ,其中 $V(r) = -Z/r$。对任意给定的 $E < 0$ 和 $l = 0,1,2,\dots$,令 $\rho = \kappa r$,$\eta = -Z/\kappa$,$\kappa = \sqrt{-2mE}$,径向方程化简为

\begin{equation} \frac{\mathrm{d}^{2}{u}}{\mathrm{d}{\rho}^{2}} + \left[-1 - \frac{2\eta}{\rho} - \frac{l(l+1)}{\rho^2} \right] u = 0~. \end{equation}
它的通解3
\begin{equation} \psi_{E,l}(x) = C_1 M_{-\eta, l+1/2}(2 \rho) + C_2 W_{-\eta, l+1/2}(2 \rho)~. \end{equation}
$M, W$ 分别是 Whittacker M 和 W 函数,$x = 0$ 时,$M$ 函数总是为零,而 $W$ 可能不为零,所以符合物理的解必须 $C_2 = 0$。另外可以证明只有对式 3 中离散的 $E_n$,波函数才能归一化得到束缚态。束缚态的径向波函数可以化简为(式 5 乘以 $1/r$)
\begin{equation} R_{n,l}(r) = \sqrt{ \left(\frac{2 Z}{n} \right) ^3 \frac{(n - l - 1)!}{2n (n + l)!}} \left(\frac{2Zr}{n} \right) ^l L_{n-l-1}^{2l+1} \left(\frac{2Zr}{n} \right) \mathrm{e} ^{-Zr/n}~, \end{equation}
其中 $L_n^l(x)$ 是连带拉盖尔多项式(associated Laguerre polynomial)。注意 $Z$ 的作用是把径向波函数关于原点缩小 $Z$ 倍,并保持波函数归一化。

   $Z = 1$ 时 $r R_{n,l}(r)$ 的函数图4图 2

图
图 2:径向波函数函数图(使用原子单位)

   以下给出前几个径向波函数,注意所有径向波函数的值都是实数。

\begin{equation} n = 1 \qquad R_{10}(r) = 2 Z^{3/2} \mathrm{e} ^{-Zr}~, \end{equation}
\begin{equation} n = 2 \qquad \left\{\begin{aligned} R_{20}(r) &= \frac{1}{\sqrt 2} Z^{3/2} \left(1 - \frac12 Zr \right) \mathrm{e} ^{-Zr/2}\\ R_{21}(r) &= \frac{1}{2\sqrt{6}} Z^{3/2} Zr \mathrm{e} ^{-Zr/2} \end{aligned}\right. ~,\end{equation}
\begin{equation} n = 3 \qquad \left\{\begin{aligned} R_{30}(r) &= \frac{2}{3\sqrt {3}} Z^{3/2} \left(1 - \frac23 Zr + \frac{2}{27} Z^2r^2 \right) \mathrm{e} ^{-Zr/3}\\ R_{31}(r) &= \frac{8}{27\sqrt 6} Z^{3/2} \left(1 - \frac16 Zr \right) Zr \mathrm{e} ^{-Zr/3}\\ R_{32}(r) &= \frac{4}{81\sqrt {30}} Z^{3/2} Z^2r^2 \mathrm{e} ^{-Zr/3} \end{aligned}\right. ~. \end{equation}
更多 $R_{n,l}$ 和 $\psi_{n,l,m}$ 表达式可以用 Mathematica 或者 Wolfram Alpha 生成,函数定义如
HydrogenR[Z_, n_, l_, r_] := Sqrt[(2Z/n)^3\
   Factorial[n-l-1]/(2n Factorial[n+l])] (2Z r/n)^l\
   LaguerreL[n-l-1, 2l+1, 2Z r/n] Exp[-Z r/n];
HydrogenPsi[Z_, n_, l_, m_, r_, θ_, ϕ_] := 
   HydrogenR[Z, n, l, r] SphericalHarmonicY[l, m, θ, ϕ];
或者用 Matlab 进行数值计算
代码 1:hydrogen_Rnl.m
function ret = hydrogen_Rnl(Z,n,l,r)
ret = sqrt((2*Z/n)^3*factorial(n-l-1)/(2*n*factorial(n+l))) *...
      (2*Z*r/n).^l .* laguerreL(n-l-1, 2*l+1, 2*Z*r/n) .* exp(-Z*r/n);
end
代码 2:hydrogen_Psi.m
function ret = hydrogen_Psi(Z,n,l,m,r,th,ph)
ret = hydrogen_Rnl(Z,n,l,r) .* SphHarm(l,m,th,ph);
end

2. 性质

   我们要求氢原子每个束缚态满足正交归一化条件

\begin{equation} \begin{aligned} &\int \Psi_{n',l',m'} ^* ( \boldsymbol{\mathbf{r}} )\Psi_{n,l,m}( \boldsymbol{\mathbf{r}} ) \,\mathrm{d}^{3}{r} \\ = &\int \int_0^\infty R_{n',l'} ^* (r) Y_{l',m'}( \hat{\boldsymbol{\mathbf{r}}} )R_{n,l}(r) Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} ) r^2 \,\mathrm{d}{r} \,\mathrm{d}{\Omega} = \delta_{n,n'}\delta_{l,l'}\delta_{m,m'}~. \end{aligned} \end{equation}
根据式 5 式 8 可得径向波函数的归一化条件
\begin{equation} \int_0^\infty R_{n', l}(r) R_{n,l}(r) r^2 \,\mathrm{d}{r} = \delta_{n,n'}~. \end{equation}
注意对不同的 $l$,无论 $n,n'$ 相同或不同,径向波函数都未必正交。

   Mathematica 代码如(式 10 中 $\Psi_{2,1,0}$ 和 $\Psi_{3,1,0}$ 的积分)

Z = 1; Integrate[
 Conjugate[HydrogenPsi[Z, 2, 1, 0, r, th, phi]] HydrogenPsi[
   Z, 3, 1, 0, r, th, phi] r^2 Sin[th], {r, 
  0, Infinity}, {th, 0, Pi}, {phi, 0, 2 Pi}]

3. 径向概率分布

   我们来求径向概率分布 $P(r)$。$P(r)$ 的定义为:发现粒子在 $r \in [a, b]$(厚球壳)内的概率为 $\int_a^b P(r) \,\mathrm{d}{r} $。由于波函数的模长平方就是三维的概率密度,有

\begin{equation} \int_a^b P(r) \,\mathrm{d}{r} = \int_a^b \int \left\lvert \frac1r \psi_{l,m}(r) Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} ) \right\rvert ^2 \,\mathrm{d}{\Omega} r^2 \,\mathrm{d}{r} = \int_a^b \left\lvert \psi_{l,m}(r) \right\rvert ^2 \,\mathrm{d}{r} ~. \end{equation}
对任意 $a, b > 0$ 都成立,所以有
\begin{equation} P(r) = \left\lvert \psi_{l,m}(r) \right\rvert ^2~. \end{equation}

   任意波函数可以表示为所有本征波函数的叠加

\begin{equation} \Psi( \boldsymbol{\mathbf{r}} ) = \frac1r \sum_{l,m} \psi_{l,m}(r) Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} )~. \end{equation}
其径向概率分布为
\begin{equation} \int_a^b P(r) \,\mathrm{d}{r} = \int_a^b \int \left\lvert \frac1r \sum_{l,m}\psi_{l,m}(r) Y_{l,m}( \hat{\boldsymbol{\mathbf{r}}} ) \right\rvert ^2 \,\mathrm{d}{\Omega} r^2 \,\mathrm{d}{r} = \sum_{l,m} \int_a^b \left\lvert \psi_{l,m}(r) \right\rvert ^2 \,\mathrm{d}{r} ~. \end{equation}
对任意 $a, b > 0$ 都成立,所以有
\begin{equation} P(r) = \sum_{l,m} \left\lvert \psi_{l,m}(r) \right\rvert ^2~. \end{equation}

4. 动量表象、动量分布

   要求动量表象下的波函数,我们需要将位置表象的波函数投影到归一化的动量的本征矢上,即三维傅里叶变换

\begin{equation} \psi_{n,l,m}( \boldsymbol{\mathbf{p}} ) = \left\langle \boldsymbol{\mathbf{p}} \middle| \psi \right\rangle = \frac{1}{\sqrt{2\pi}} \int \exp\left(- \mathrm{i} \boldsymbol{\mathbf{p}} \boldsymbol\cdot \boldsymbol{\mathbf{r}} \right) \psi( \boldsymbol{\mathbf{r}} ) \,\mathrm{d}^{3}{r} ~. \end{equation}
这个积分在球坐标中完成才是最方便的,具体方法我们将举例子说明(见例 1 )。

   正如位置表象下位置的分布函数是 $ \left\lvert \psi( \boldsymbol{\mathbf{r}} ) \right\rvert ^2$,动量表象下动量的分布函数是 $ \left\lvert \psi( \boldsymbol{\mathbf{p}} ) \right\rvert ^2$(也符合测量理论)。


1. ^ 如果不要求波函数可归一化,$E$ 取任意实数我们都能找到解。但只有当 $E$ 取特定离散值的时候这些波函数才能归一化,不能归一化的波函数没有物理意义。
2. ^ 为什么图中用 95% 而不是 100%?因为无论离原子核多远电子出现的概率都不会完全等于零,而是一个指数下降的尾巴(exponential tail)。
3. ^ 可以使用 Mathematica 求解
4. ^ 我们以后会看到 $r R_{l,m}(r)$ 比 $R_{l,m}(r)$ 更常用


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利