# 经典力学中的哈密顿量

氢原子由带正电的原子核和带负电的电子组成。在经典力学中,其总哈密顿量可以表示为质心运动动能与质心系内电子、原子核动能以及它们之间的相互作用势能之和。

H=12MR2+12μr2+U(r)H = \frac12 M R^2 + \frac12 \mu r^2 + U(r)

其中,MM 是原子核和电子的总质量,RR 是质心位置;μ\mu 是约化质量,rr 是电子与原子核之间的相对位置。U(r)U(r) 是势能,通常表示为库仑势。

# 量子力学中的哈密顿量与波函数

在量子力学中,氢原子的哈密顿量可以写成算符形式:

H=22MR222μ2+U(r)H = -\frac{\hbar^2}{2M} \nabla^2_R - \frac{\hbar^2}{2\mu} \nabla^2 + U(r)

其中,第一项描述质心运动,第二项和第三项描述电子与原子核在质心系中的相对运动。通过分离变量法,其总波函数可以分解为质心波函数 ϕ(R)\phi(R)、相对波函数 ψ(r)\psi(r) 和时间演化部分。

Ψ(r,R,t)=ψ(r)ϕ(R)eiEt\Psi(r, R, t) = \psi(r) \phi(R) e^{-\frac{i}{\hbar} E_{总} t}

# 相对运动的动能算符

用于描述相对运动的动能算符可以进一步分解为径向部分和角动量部分:

22μ2=22μr2r(r2r)+12μr2L^2-\frac{\hbar^2}{2\mu} \nabla^2 = -\frac{\hbar^2}{2\mu r^2} \frac{\partial}{\partial r}\left( r^2 \frac{\partial}{\partial r} \right) + \frac{1}{2 \mu r^2} \hat L^2

其中,第一项为径向动能算符,第二项为角动量动能算符,L^2\hat L^2 是角动量算符的平方。


# 氢原子的本征函数和本征值

# 本征函数

氢原子在质心系中的本征函数 ψnlm(r,θ,ϕ)\psi_{nlm}(r, \theta, \phi) 可以用径向函数 Rnl(r)R_{nl}(r) 和球谐函数 Ylm(θ,ϕ)Y_{lm}(\theta, \phi) 的乘积来表示,即:

ψnlm(r,θ,ϕ)=Rnl(r)Ylm(θ,ϕ)\psi_{nlm}(r, \theta, \phi) = R_{nl}(r)Y_{lm}(\theta, \phi)

其中,量子数 n,l,mn, l, m 满足:nn 是主量子数,ll 是角量子数(l=0,1,,n1l = 0, 1, \cdots, n-1),mm 是磁量子数。

径向函数 Rnl(r)R_{nl}(r) 的表达式为:

Rnl(r)=(2na)3(nl1)!2n(n+l)!×(2rna)lernaLnl12l+1(2rna)R_{nl}(r) = \sqrt{\left( \frac{2}{na} \right)^3 \frac{(n-l-1)!}{2n(n+l)!}} \times \left( \frac{2r}{na} \right)^l e^{-\frac{r}{na}} L^{2l+1}_{n-l-1} \left( \frac{2r}{na} \right)

这里,aa 是玻尔半径,a=2μke20.53Aa = \frac{\hbar^2}{\mu k e^2} \approx 0.53 \overset{\circ}{A}Lmk(x)L_m^k(x) 是广义拉盖尔多项式,其定义为:

Lmk(x)=exm!xkdmdxm(xm+kex)L_m^k(x) = \frac{e^x}{m!x^k} \frac{d^m}{dx^m} \left( \frac{x^{m+k}}{e^x} \right)


# 本征值

氢原子的能量本征值 EnE_n 只与主量子数 nn 有关,其表达式为:

En=E1n2=μk2Z2e4221n2E_n = \frac{E_1}{n^2} = -\frac{\mu k^2Z^2e^4}{2\hbar^2} \frac{1}{n^2}

其中,k=14πε0k = \frac{1}{4\pi\varepsilon_0} 是库仑常数,ZZ 是原子序数(氢原子 Z=1Z=1)。基态能量 E113.6 eVE_1 \approx -13.6 \text{ eV}


# 氢原子本征函数举例

  • 基态 (n=1,l=0,m=0n=1, l=0, m=0)

    ψ100=4a3eraY00\psi_{100} = \sqrt{\frac{4}{a^3}} e^{-\frac{r}{a}} Y_{00}

  • 第一激发态 (n=2n=2)

    • l=0,m=0l=0, m=0

      ψ200=12a3(1r2a)er2aY00\psi_{200} = \sqrt{\frac{1}{2a^3}} \left( 1-\frac{r}{2a} \right) e^{-\frac{r}{2a}} Y_{00}

    • l=1l=1

      ψ21m=16a3(r2a)er2aY1m(m=1,0,1)\psi_{21m} = \sqrt{\frac{1}{6a^3}} \left( \frac{r}{2a} \right) e^{-\frac{r}{2a}} Y_{1m} \quad (m=-1, 0, 1)