跳转至

连带拉盖尔多项式

连带拉盖尔多项式与氢原子束缚态

物理问题引入:为什么氢原子的能级严格正比于 −1/n2-1/n^2?

1913 年尼尔斯·玻尔(Niels Bohr) 提出原子行星轨道模型,手工假设了角动量量子化条件 L=nℏL = nhbar,猜出了氢原子的里德伯能级公式:En∝−1/n2E_n propto -1/n^2. 然而,在严格的三维量子波动力学中,电子是一个遵从三维薛定谔方程演化的空间概率云:

−ℏ22μ∇2ψ(r)−e24πε0rψ(r)=Eψ(r)-hbar^22munabla^2psi(r) - e^24pivarepsilon_0 rpsi(r) = Epsi(r)

在球坐标系下分离变量后:

  • 角向部分的解正是前文探讨过的 球谐函数 Ylm(θ,ϕ)Y_l^m(theta, phi);
  • 剩下的最硬核部分是 径向波函数 R(r)R(r):它决定了电子在距离原子核不同半径处的发现概率;
  • 玻尔能级公式中的整数 nn(主量子数),到底是从什么样的数学边界中自然生长出来的?答案就是连带拉盖尔多项式(Associated Laguerre Polynomials).

1. 径向薛定谔方程与物理渐近剥离

对于束缚态(总能量 E<0E < 0),引入无量纲径向距离 ρ=2κrrho = 2kappa r(其中 κ=−2μEℏkappa = -2mu Ehbar),径向方程写为:

d2Rdρ2+2ρdRdρ+[λρ−14−l(l+1)ρ2]R=0d^2 Rdrho^2 + 2rhodRdrho + left[ rho - 14 - l(l+1)rho^2 right] R = 0

为了使解在全实空间保持良定,必须进行 两端渐近物理分析:

  1. 无穷远处渐近 (ρ→∞rho -> inf):方程退化为 R″−14R≈0R'' - 1/4R ~ 0,为保证波函数平方可积,必须保留指数衰减项 R(ρ)∝e−ρ/2R(rho) propto e^-rho/2;
  2. 原子核中心渐近 (ρ→0rho -> 0):离心势能垒项 −l(l+1)ρ2-l(l+1)/rho^2 占主导,代入试探解可得 R(ρ)∝ρlR(rho) propto rho^l.

因此,将波函数两端的渐近行为完全「剥离」出来,设全解形式为:

R(ρ)=ρle−ρ/2v(ρ)R(rho) = rho^l e^-rho/2 v(rho)

代入原方程,未知函数 v(ρ)v(rho) 所满足的微分方程化为标准的 连带拉盖尔方程.


2. 连带拉盖尔微分方程与级数截断

标准连带拉盖尔方程定义为:

xd2ydx2+(k+1−x)dydx+py=0x d^2 ydx^2 + (k + 1 - x)dydx + p , y = 0

若使用幂级数法求解 v(ρ)=∑amρmv(rho) = sum a_m rho^m,递推关系分析表明:若级数不截断,其远处发散行为如同 eρe^rho,与前面的因子 e−ρ/2e^-rho/2 相乘后将变成净增长的 e+ρ/2e^+rho/2,彻底破坏无穷远归一化!波函数物理归一化再次强制发出指令:级数必须在有限项截断! 截断条件要求参数 pp 必须是非负整数:

p=λ−l−1=nr∈{0,1,2,…}p = lambda - l - 1 = n_r in 0, 1, 2, dots

定义 主量子数(Principal Quantum Number) 为:

n=nr+l+1n = n_r + l + 1

由于径向节点数 nr≥0n_r >= 0,因此主量子数与轨道角动量量子数必须严格满足量子力学约束:

n≥l+1⟺l=0,1,2,…,n−1n >= l + 1 quad iff quad l = 0, 1, 2, dots, n-1

由此解出的有限项多项式就是 连带拉盖尔多项式 Lpk(x)L_p^k(x)(其中 k=2l+1,p=n−l−1k = 2l+1, p = n - l - 1).


3. 罗德里格斯公式与正交性质

3.1 罗德里格斯表示式

普通拉盖尔多项式由下式定义:

Lp(x)=exp!dpdxp(xpe−x)L_p(x) = e^xp! d^pdx^pleft( x^p e^-x right)

连带拉盖尔多项式为普通拉盖尔多项式的 kk 阶微商:

Lpk(x)=(−1)kdkdxkLp+k(x)=x−kexp!dpdxp(xp+ke−x)L_p^k(x) = (-1)^k d^kdx^k L_p+k(x) = x^-ke^xp! d^pdx^pleft( x^p+k e^-x right)

3.2 带权正交性

在半无穷区间 [0,∞)[0, inf) 上,以 w(x)=xke−xw(x) = x^k e^-x 为权函数:

∫0∞Lpk(x)Lqk(x)xke−xdx=(p+k)!p!δpqint_0^inf L_p^k(x) L_q^k(x) , x^k e^-x , dx = (p + k)!p! , delta_pq

4. 物理皇冠:玻尔能级与电子轨道云

将无量纲截断参数 λ=nlambda = n 还原回物理常数定义:

n=μe24πε0ℏ2κ⟹κ=μe24πε0ℏ21n=1na0n = mu e^24pivarepsilon_0 hbar^2 kappa implies kappa = mu e^24pivarepsilon_0 hbar^2 1n = 1n a_0

其中 a0=4πε0ℏ2μe2≈0.529 \AAa_0 = 4pivarepsilon_0 hbar^2mu e^2 ~ 0.529 AA 即为著名的 玻尔半径! 代入能量公式 E=−ℏ2κ22μE = -hbar^2 kappa^22mu,纯数学截断条件精准导出了宏伟的 里德伯 - 玻尔能级公式:

En=−μe432π2ε02ℏ21n2=−13.6 eVn2(n=1,2,3,…)E_n = -mu e^432pi^2 varepsilon_0^2 hbar^2 1n^2 = -13.6 eVn^2 quad (n = 1, 2, 3, dots)

氢原子的完整三维归一化定态波函数

ψnlm(r,θ,ϕ)=Rnl(r)Ylm(θ,ϕ)psi_nlm(r, theta, phi) = R_nl(r) Y_l^m(theta, phi)

其中径向波函数完全由连带拉盖尔多项式支配:

Rnl(r)=−(2na0)3(n−l−1)!2n[(n+l)!]3exp⁡(−rna0)(2rna0)lLn−l−12l+1(2rna0)R_nl(r) = - left(2n a_0right)^3 (n - l - 1)!2n [(n + l)!]^3 expleft(-rn a_0right) left(2rn a_0right)^l L_n-l-1^2l+1left(2rn a_0right)
  • 电子在距核 rr 到 r+drr+dr 球壳内的发现概率为 P(r)dr=r2|Rnl(r)|2drP(r)dr = r^2 |R_nl(r)|^2 dr;
  • 径向多项式的节点数 nr=n−l−1n_r = n - l - 1 精确决定了电子云同心球壳节面的数目!

至此,大一化学中背诵的电子轨道能量与其三维空间分布,在数学物理方法中完全显露出了优雅而必然的本原逻辑.


5. 学习衔接

  • 前驱复习:全系特殊函数本征正交性的统一母体,参见 Sturm–Liouville 理论;
  • 课程路线:至此,数理方法核心工具链(复分析 →-> 傅里叶变换 →-> PDE 分离变量 →-> Sturm-Liouville →-> 特殊函数)已全部闭环!


评论