「特殊函数」这个名字听起来吓人,其实它们只是几个常见微分方程的解,
被前人算好、编号、列成了表。薛定谔方程在几种规则势场里分离变量后,
剩下的常微分方程翻来覆去就那么几个——每个都对应一族特殊函数。
本节按族整理:解什么方程、怎么定义、递推关系、前几个长什么样、正文哪一章用。
先给一张总索引:
| 函数族 | 物理场景 | 出场章节 |
|---|
| Hermite Hn | 谐振子本征态 | 2.9 |
| Legendre Pℓ 与球谐 Yℓm | 角动量、任何中心力场的角向部分 | 5.1、5.2 |
| Laguerre Lnk | 氢原子径向波函数 | 6.3 |
| 球 Bessel jℓ,nℓ | 自由粒子的径向解、三维散射 | 6.1、6.2 |
Hermite 多项式
解决的问题:谐振子的定态方程。代入无量纲变量 ξ=mω/ℏx
并剥离高斯衰减因子 ψ=h(ξ)e−ξ2/2 后,h 满足 Hermite 方程
h′′−2ξh′+2nh=0(A.3.1)
只有 n 取非负整数时解才是多项式(否则级数不截断、ψ 发散)——
这正是谐振子能级 En=(n+21)ℏω 量子化的来源。
定义(Rodrigues 公式,一条式子生成全家):
Hn(ξ)=(−1)neξ2dξndne−ξ2(A.3.2)
递推与求导(数值和手算都靠这两条,不用背显式表达式):
Hn+1=2ξHn−2nHn−1,Hn′=2nHn−1(A.3.3)
前几个:
| n | Hn(ξ) |
|---|
| 0 | 1 |
| 1 | 2ξ |
| 2 | 4ξ2−2 |
| 3 | 8ξ3−12ξ |
| 4 | 16ξ4−48ξ2+12 |
正交性(权函数是高斯):
∫−∞∞HmHne−ξ2dξ=2nn!πδmn(A.3.4)
由此得到归一化的谐振子本征态
ψn(x)=(πℏmω)1/42nn!1Hn(ξ)e−ξ2/2。
Hn 是 n 次多项式、宇称为 (−1)n、有 n 个实零点——
对应 ψn 有 n 个节点。
Legendre 多项式与球谐函数
解决的问题:任何中心力场 V(r) 中分离变量后的角向方程——
它与 V 无关,所以角向解是通用的。θ 方向的方程(令 u=cosθ)是 Legendre 方程
dud[(1−u2)dudP]+ℓ(ℓ+1)P=0(A.3.5)
在两极 u=±1 处解要保持有限,逼出 ℓ=0,1,2,…——角动量量子化的来源。
定义与递推:
Pℓ(u)=2ℓℓ!1duℓdℓ(u2−1)ℓ,(ℓ+1)Pℓ+1=(2ℓ+1)uPℓ−ℓPℓ−1(A.3.6)
前几个:P0=1,P1=u,P2=21(3u2−1),P3=21(5u3−3u)。
带上 m 的连带 Legendre 函数 Pℓm(u)=(1−u2)∣m∣/2d∣m∣Pℓ/du∣m∣,
再配上 φ 方向的 eimφ 和归一化系数,就是球谐函数——
角动量算符 L^2,L^z 的公共本征函数:
L^2Yℓm=ℓ(ℓ+1)ℏ2Yℓm,L^zYℓm=mℏYℓm,m=−ℓ,…,ℓ(A.3.7)
最常用的几个(正文 5.2 节有三维图像):
| Yℓm | 表达式 |
|---|
| Y00 | 1/4π |
| Y10 | 3/4πcosθ |
| Y1±1 | ∓3/8πsinθe±iφ |
| Y20 | 5/16π(3cos2θ−1) |
正交归一(在整个球面上积分):
∫02π∫0πYℓ′m′∗Yℓmsinθdθdφ=δℓℓ′δmm′(A.3.8)
化学里的 px,py,pz 轨道就是 Y1±1,Y10 的实线性组合。
Laguerre 多项式
解决的问题:氢原子径向方程。剥离长程衰减 e−r/na 与短程行为 rℓ 之后,
剩下的多项式部分满足连带 Laguerre 方程。级数截断条件给出主量子数 n,
这就是氢原子能级 En=−13.6eV/n2 的来源
(完整推导见 6.3 节)。
定义(Rodrigues 公式;Lqk 是 q 次多项式):
Lq(x)=q!exdxqdq(xqe−x),Lqk(x)=(−1)kdxkdkLq+k(x)(A.3.9)
前几个:L0=1,L1=1−x,L2=1−2x+2x2;
L01=1,L11=2−x。
氢原子径向波函数的组装方式(a 为 Bohr 半径):
Rnℓ(r)∝(na2r)ℓLn−ℓ−12ℓ+1(na2r)e−r/na(A.3.10)
多项式部分有 n−ℓ−1 个零点,正是径向节点数。
球 Bessel 函数
解决的问题:势能为零(或常数)区域的径向方程。令 u=rR 之后径向方程变成
R′′+r2R′+[k2−r2ℓ(ℓ+1)]R=0,
解是球 Bessel 函数 jℓ(kr)(原点处规则)与球 Neumann 函数 nℓ(kr)(原点处发散)。
低阶显式表达式(都是三角函数除以幂次,可以手算):
| ℓ | jℓ(x) | nℓ(x) |
|---|
| 0 | xsinx | −xcosx |
| 1 | x2sinx−xcosx | −x2cosx−xsinx |
递推(fℓ 代表 jℓ 或 nℓ,同一条递推):
fℓ+1(x)=x2ℓ+1fℓ(x)−fℓ−1(x)(A.3.11)
两端行为(判断解的取舍全靠这两条):
jℓ(x)x→0(2ℓ+1)!!xℓ,jℓ(x)x→∞xsin(x−ℓπ/2)(A.3.12)
nℓ 在原点发散如 x−ℓ−1,所以包含原点的区域只能用 jℓ;
远处的正弦振荡(相位平移 ℓπ/2)是三维散射分波法里「相移」概念的参照零点。
∎本节关键公式
Hermite 递推
Hn+1=2ξHn−2nHn−1,H0=1, H1=2ξ 谐振子;ψₙ 有 n 个节点,宇称 (−1)ⁿ
Legendre 递推
(ℓ+1)Pℓ+1=(2ℓ+1)uPℓ−ℓPℓ−1 中心力场角向方程;两极有限 ⇒ ℓ 取整数
球谐本征方程
L^2Yℓm=ℓ(ℓ+1)ℏ2Yℓm,L^zYℓm=mℏYℓm 所有中心势共用的角向解
氢原子径向解
Rnℓ∝(2r/na)ℓLn−ℓ−12ℓ+1(2r/na)e−r/na 径向节点数 = n − ℓ − 1
球 Bessel 渐近
jℓx→0(2ℓ+1)!!xℓ,jℓx→∞xsin(x−ℓπ/2) 含原点区域舍弃 nℓ;远端相位是分波相移的零点
1.谐振子第 3 激发态的波函数 ψ₃ 有几个节点?宇称是?
2.为什么包含 r = 0 的区域里,径向解必须取 jℓ 而不能混入 nℓ?
3.氢原子 3p 态(n = 3, ℓ = 1)的径向波函数有几个径向节点?