分析:
-
为什么角动量算符能得到角动量的量子化
-
为什么径向方程能得到波尔的能量量子化一样的结果
氢原子薛定谔方程解的详细推导
氢原子是量子力学中少数几个可以精确求解的体系之一。下面将详细介绍如何通过分离变量法求解氢原子的定态薛定谔方程,得到其波函数和能级。
1. 氢原子的薛定谔方程
氢原子的哈密顿量为: \[ \hat{H} = -\frac{\hbar^2}{2\mu}\nabla^2 - \frac{e^2}{4\pi\epsilon_0 r} \] 其中 \(\mu = \frac{m_em_p}{m_e+m_p}\) 是约化质量,\(r\) 是电子与质子间的距离。
定态薛定谔方程为: \[ \left[ -\frac{\hbar^2}{2\mu}\nabla^2 - \frac{e^2}{4\pi\epsilon_0 r} \right]\psi(\mathbf{r}) = E\psi(\mathbf{r}) \]
2. 球坐标系下的分离变量
在球坐标系 \((r,\theta,\phi)\) 中,拉普拉斯算符为: \[ \nabla^2 = \frac{1}{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial}{\partial r}\right) + \frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{r^2\sin^2\theta}\frac{\partial^2}{\partial\phi^2} \]
设波函数为: \[ \psi(r,\theta,\phi) = R(r)Y(\theta,\phi) \]
分离变量后得到:
-
径向方程: \[ -\frac{\hbar^2}{2\mu}\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{dR}{dr}\right) + \left[ \frac{\hbar^2\ell(\ell+1)}{2\mu r^2} - \frac{e^2}{4\pi\epsilon_0 r} \right]R = ER \]
-
角向方程(球谐函数方程): \[ \frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial Y}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2 Y}{\partial\phi^2} = -\ell(\ell+1)Y \]
3. 角向部分的解:球谐函数
角向方程的解是球谐函数 \(Y_\ell^m(\theta,\phi)\): \[ Y_\ell^m(\theta,\phi) = \sqrt{\frac{(2\ell+1)(\ell-|m|)!}{4\pi(\ell+|m|)!}} P_\ell^{|m|}(\cos\theta) e^{im\phi} \] 其中:
- \(\ell = 0,1,2,\cdots\) 是角量子数
- \(m = -\ell,-\ell+1,\cdots,\ell\) 是磁量子数
- \(P_\ell^m\) 是连带勒让德多项式
4. 径向方程的求解
引入无量纲变量: \[ \rho = \frac{r}{a_0}, \quad a_0 = \frac{4\pi\epsilon_0\hbar^2}{\mu e^2} \text{(玻尔半径)}, \quad \lambda^2 = -\frac{2\mu E a_0^2}{\hbar^2} \]
径向方程变为: \[ \frac{d^2 u}{d\rho^2} + \left[ \frac{2}{\rho} - \frac{\ell(\ell+1)}{\rho^2} - \lambda^2 \right]u = 0 \] 其中 \(u(\rho) = \rho R(\rho)\)。
4.1 渐近行为分析
- 当 \(\rho \to \infty\):方程近似为 \(\frac{d^2 u}{d\rho^2} - \lambda^2 u = 0\),解为 \(u \sim e^{-\lambda\rho}\)
- 当 \(\rho \to 0\):方程近似为 \(\frac{d^2 u}{d\rho^2} - \frac{\ell(\ell+1)}{\rho^2}u = 0\),解为 \(u \sim \rho^{\ell+1}\)
因此设通解形式: \[ u(\rho) = \rho^{\ell+1} e^{-\lambda\rho} F(\rho) \]
4.2 级数解法
代入后得到 \(F(\rho)\) 的方程: \[ \rho \frac{d^2 F}{d\rho^2} + [2(\ell+1) - 2\lambda\rho] \frac{dF}{d\rho} + [2 - 2\lambda(\ell+1)] F = 0 \]
设级数解: \[ F(\rho) = \sum_{k=0}^\infty c_k \rho^k \]
得到递推关系: \[ c_{k+1} = \frac{2\lambda(k+\ell+1)-2}{(k+1)(k+2\ell+2)} c_k \]
4.3 量子化条件
为使级数截断为多项式,要求: \[ 2\lambda(n_r + \ell + 1) - 2 = 0 \] 其中 \(n_r\) 是径向量子数。定义主量子数 \(n = n_r + \ell + 1\),则: \[ \lambda = \frac{1}{n} \quad \Rightarrow \quad E_n = -\frac{\mu e^4}{2(4\pi\epsilon_0)^2\hbar^2 n^2} \]
4.4 径向波函数
最终径向波函数为: \[ R_{n\ell}(r) = \sqrt{\left(\frac{2}{na_0}\right)^3 \frac{(n-\ell-1)!}{2n[(n+\ell)!]^3}} \left(\frac{2r}{na_0}\right)^\ell e^{-r/na_0} L_{n-\ell-1}^{2\ell+1}\left(\frac{2r}{na_0}\right) \] 其中 \(L_p^q\) 是连带拉盖尔多项式。
5. 完整的氢原子波函数
结合径向和角向部分,氢原子的定态波函数为: \[ \psi_{n\ell m}(r,\theta,\phi) = R_{n\ell}(r) Y_\ell^m(\theta,\phi) \]
6. 能级和量子数
- 能级仅与主量子数 \(n\) 有关: \[ E_n = -\frac{13.6,\text{eV}}{n^2} \]
- 量子数的取值范围:
- \(n = 1,2,3,\cdots\)
- \(\ell = 0,1,\cdots,n-1\)
- \(m = -\ell,-\ell+1,\cdots,\ell\)
7. 前几个波函数的显式形式
-
基态 (\(n=1,\ell=0,m=0\)): \[ \psi_{100} = \frac{1}{\sqrt{\pi a_0^3}} e^{-r/a_0} \]
-
第一激发态 (\(n=2\)):
- \(\ell=0\): \[ \psi_{200} = \frac{1}{4\sqrt{2\pi a_0^3}} \left(2 - \frac{r}{a_0}\right) e^{-r/2a_0} \]
- \(\ell=1\): \[ \psi_{21m} = \frac{1}{4\sqrt{2\pi a_0^3}} \frac{r}{a_0} e^{-r/2a_0} Y_1^m(\theta,\phi) \]
8. 解的物理意义
- 量子化能级:能量只能取离散值,解释了氢原子光谱的分立性。
- 角动量量子化:轨道角动量 \(L^2 = \hbar^2\ell(\ell+1)\)。
- 空间量子化:角动量在z方向的投影 \(L_z = \hbar m\)。
- 电子云分布:\(|\psi_{n\ell m}|^2\) 给出电子在不同位置出现的概率密度。
=====================qwen:
以下是薛定谔方程求解氢原子的完整、详细推导过程,从建立方程到最终得到能级和波函数,每一步都清晰解释其物理意义与数学逻辑。这是量子力学中最经典、最重要的精确解之一。
🧩 薛定谔方程求解氢原子全过程
一、物理模型与哈密顿量
氢原子由一个质子(电荷 \(+e\))和一个电子(电荷 \(-e\))组成。忽略相对论效应和自旋,电子在库仑势中运动。
1. 势能函数(库仑势):
\[ V(r) = -\frac{1}{4\pi\varepsilon_0} \frac{e^2}{r} \]
2. 哈密顿算符(动能 + 势能):
\[ \hat{H} = -\frac{\hbar^2}{2\mu} \nabla^2 + V(r) \] 其中 \(\mu\) 是电子与质子的约化质量: \[ \mu = \frac{m_e m_p}{m_e + m_p} \approx m_e \quad (\text{因为 } m_p \gg m_e) \]
二、不含时薛定谔方程
我们求定态解: \[ \hat{H} \psi(\vec{r}) = E \psi(\vec{r}) \] 即: \[ \left( -\frac{\hbar^2}{2\mu} \nabla^2 - \frac{e^2}{4\pi\varepsilon_0 r} \right) \psi = E \psi \]
三、使用球坐标分离变量
由于势能 \(V(r)\) 是球对称的,使用球坐标 \((r, \theta, \varphi)\),拉普拉斯算符为:
\[ \nabla^2 = \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial}{\partial r} \right) + \frac{1}{r^2 \sin\theta} \frac{\partial}{\partial \theta} \left( \sin\theta \frac{\partial}{\partial \theta} \right) + \frac{1}{r^2 \sin^2\theta} \frac{\partial^2}{\partial \varphi^2} \]
设波函数可分离变量: \[ \psi(r, \theta, \varphi) = R(r) Y(\theta, \varphi) \]
更进一步,角向部分也可分离: \[ Y(\theta, \varphi) = \Theta(\theta) \Phi(\varphi) \]
但更标准的做法是直接设: \[ \psi(r, \theta, \varphi) = R(r) Y_{\ell m}(\theta, \varphi) \] 其中 \(Y_{\ell m}\) 是球谐函数。
四、代入方程并分离变量
将 \(\psi = R(r) Y_{\ell m}(\theta, \varphi)\) 代入薛定谔方程,并利用角动量算符性质:
\[ \hat{L}^2 Y_{\ell m} = \ell(\ell+1)\hbar^2 Y_{\ell m} \]
可得径向方程(推导略,结果如下):
\[ -\frac{\hbar^2}{2\mu} \frac{d^2 u}{dr^2} + \left[ -\frac{e^2}{4\pi\varepsilon_0 r} + \frac{\ell(\ell+1)\hbar^2}{2\mu r^2} \right] u(r) = E u(r) \]
其中定义了径向波函数: \[ u(r) = r R(r) \]
这个方程称为径向薛定谔方程,\(u(r)\) 的物理意义是:在半径 \(r\) 处球壳内找到电子的概率正比于 \(|u(r)|^2\)。
有效势为: \[ V_{\text{eff}}(r) = -\frac{k}{r} + \frac{\ell(\ell+1)\hbar^2}{2\mu r^2}, \quad k = \frac{e^2}{4\pi\varepsilon_0} \]
- 第一项:库仑吸引
- 第二项:离心势(阻止电子撞向核)
五、分析边界条件
我们要求波函数物理合理:
- \(u(r) \to 0\) 当 \(r \to 0\)
- \(u(r) \to 0\) 当 \(r \to \infty\)(束缚态)
- 连续、单值、有限、可归一化
1. 当 \(r \to 0\):
- 若 \(\ell > 0\),离心势主导,解得 \(u(r) \sim r^{\ell+1}\)
- 若 \(\ell = 0\),\(u(r) \sim r\)
所以一般有: \[ u(r) \sim r^{\ell+1} \quad \text{as } r \to 0 \]
2. 当 \(r \to \infty\):
- \(V(r) \to 0\),方程近似为: \[ \frac{d^2 u}{dr^2} \approx -\frac{2\mu E}{\hbar^2} u \]
- 对束缚态 \(E < 0\),令 \(\kappa = \sqrt{-2\mu E}/\hbar\),则: \[ u(r) \sim e^{-\kappa r} \quad \text{as } r \to \infty \]
六、变量替换与无量纲化
为简化方程,引入无量纲变量。
定义: \[ \rho = \kappa r = \frac{\sqrt{-2\mu E}}{\hbar} r \]
再引入另一个尺度参数。定义: \[ \lambda = \frac{\mu k}{\hbar^2 \kappa} = \frac{\mu}{\hbar^2} \cdot \frac{e^2}{4\pi\varepsilon_0} \cdot \frac{\hbar}{\sqrt{-2\mu E}} = \frac{\mu e^2}{4\pi\varepsilon_0 \hbar \sqrt{-2\mu E}} \]
这个 \(\lambda\) 将成为关键参数。
现在设: \[ u(r) = e^{-\rho} \rho^{\ell+1} s(\rho) \] 这个形式满足:
- \(e^{-\rho}\):保证 \(r \to \infty\) 时衰减
- \(\rho^{\ell+1}\):保证 \(r \to 0\) 时行为正确
- \(s(\rho)\):待定函数
七、代入并化为合流超几何方程
将 \(u = e^{-\rho} \rho^{\ell+1} s(\rho)\) 代入径向方程,经过繁琐但标准的微分运算(略),得到:
\[ \rho \frac{d^2 s}{d\rho^2} + (2\ell + 2 - \rho) \frac{ds}{d\rho} + (\lambda - \ell - 1) s = 0 \]
这正是合流超几何方程(Confluent Hypergeometric Equation)的标准形式,其通解为合流超几何函数: \[ s(\rho) = {}_1F_1(\lambda - \ell - 1, 2\ell + 2; \rho) \]
八、物理条件:要求解可归一化 → 量子化!
合流超几何函数 \({}_1F_1(a,c;\rho)\) 在 \(\rho \to \infty\) 时行为为: \[ {}_1F_1(a,c;\rho) \sim e^{\rho} \rho^{a-c} \]
但我们已经有一个 \(e^{-\rho}\) 因子,所以总行为: \[ u(r) \sim e^{-\rho} \cdot e^{\rho} \cdot \rho^{a} = \rho^{a} \] 除非 \(a\) 是非正整数,否则 \(u(r)\) 在无穷远处不衰减!
要使解物理,必须让级数在有限项后截断,变成一个多项式。这要求: \[ \lambda - \ell - 1 = n_r \in \mathbb{Z}_{\geq 0} \] 即: \[ \lambda = n_r + \ell + 1 \]
定义主量子数: \[ n = n_r + \ell + 1, \quad n = 1,2,3,\dots \] 其中 \(n_r = 0,1,2,\dots\) 是径向量子数(节点数)。
九、导出能量量子化(关键一步!)
由: \[ \lambda = \frac{\mu e^2}{4\pi\varepsilon_0 \hbar \sqrt{-2\mu E}} = n \]
解出 \(E\):
先平方两边: \[ \lambda^2 = \frac{\mu^2 e^4}{(4\pi\varepsilon_0)^2 \hbar^2 \cdot (-2\mu E)} = \frac{\mu e^4}{(4\pi\varepsilon_0)^2 \hbar^2 \cdot (-2E)} \]
所以: \[ E = -\frac{\mu e^4}{2(4\pi\varepsilon_0)^2 \hbar^2} \cdot \frac{1}{\lambda^2} = -\frac{\mu e^4}{2(4\pi\varepsilon_0)^2 \hbar^2} \cdot \frac{1}{n^2} \]
定义里德伯能量: \[ E_{\text{Ryd}} = \frac{\mu e^4}{2(4\pi\varepsilon_0)^2 \hbar^2} \approx 13.6\ \text{eV} \]
最终得到氢原子能级: \[ \boxed{E_n = -\frac{13.6\ \text{eV}}{n^2}} \]
十、波函数的完整形式
1. 径向部分 \(R_{n\ell}(r)\)
由 \(u(r) = r R(r) = e^{-\rho} \rho^{\ell+1} s(\rho)\),且 \(s(\rho)\) 是多项式(广义拉盖尔多项式),可得:
\[ R_{n\ell}(r) = N_{n\ell} , e^{-\rho/2} \rho^{\ell} L_{n-\ell-1}^{2\ell+1}(\rho) \]
其中:
- \(\rho = \frac{2r}{n a_0}\),\(a_0 = \frac{4\pi\varepsilon_0 \hbar^2}{\mu e^2}\) 是玻尔半径
- \(L_{k}^{\alpha}(\rho)\) 是广义拉盖尔多项式
- \(N_{n\ell}\) 是归一化常数
2. 角向部分:球谐函数
\[ Y_{\ell m}(\theta, \varphi) = \text{球谐函数} \]
3. 完整波函数:
\[ \psi_{n\ell m}(r,\theta,\varphi) = R_{n\ell}(r) Y_{\ell m}(\theta,\varphi) \]
十一、量子数总结
量子数 | 符号 | 取值范围 | 物理意义 |
---|---|---|---|
主量子数 | \(n\) | \(1,2,3,\dots\) | 决定能量 \(E_n \propto -1/n^2\) |
角量子数 | \(\ell\) | \(0,1,2,\dots,n-1\) | 决定角动量大小 \(L = \sqrt{\ell(\ell+1)}\hbar\) |
磁量子数 | \(m\) | \(-\ell,\dots,+\ell\) | 决定角动量方向 \(L_z = m\hbar\) |
十二、物理结果验证
1. 能级与实验一致
- 计算出的 \(E_n\) 与氢原子光谱(如巴尔末系)完全吻合
- 里德伯常数理论值与实验值高度一致
2. 波函数可归一化
- 所有 \(\psi_{n\ell m}\) 满足 \(\int |\psi|^2 d^3r = 1\)
3. 概率分布合理
- 如 \(1s\) 态:最概然半径为 \(a_0\)
- \(2p\) 态:哑铃形电子云,符合化学图像
✅ 总结:求解流程图
薛定谔方程
↓
球坐标分离变量 → ψ = R(r) Y(θ,φ)
↓
角向解:球谐函数 Y_{ℓm} → ℓ, m 量子化
↓
径向方程 + 边界条件
↓
渐近分析:u ∼ r^{ℓ+1} (r→0), e^{-κr} (r→∞)
↓
设 u = e^{-ρ} ρ^{ℓ+1} s(ρ)
↓
得到合流超几何方程
↓
要求 s(ρ) 为多项式 → λ - ℓ - 1 = n_r
↓
定义 n = n_r + ℓ + 1
↓
λ = n → 代入能量表达式
↓
E_n = -常数 / n²
↓
得到完整波函数 ψ_{nℓm}(r,θ,φ)
📚 参考资料
- Griffiths, D.J. Introduction to Quantum Mechanics(第4章)
- Landau & Lifshitz, Vol. 3 Quantum Mechanics
- Schrödinger, E. (1926). Quantisierung als Eigenwertproblem
- Cohen-Tannoudji, Diu & Laloë Quantum Mechanics, Vol. 1
🌟 结语
氢原子的薛定谔方程求解,是人类第一次从第一性原理出发,精确解释原子光谱的壮举。它不仅验证了波动力学的正确性,更揭示了量子化、对称性、概率解释的深刻本质。这个解至今仍是量子物理的“基石”与“灯塔”。
\[ \frac{2m_e E a_0^2}{\hbar^2} = -\frac{1}{n^2} \] (注意:由于 \(E < 0\),左边应为负值)
我们一步步来验证这个表达式是否成立,以及它是否隐含在玻尔的原始推导中。
一、先回顾玻尔模型的基本结果
1. 玻尔假设:
- 角动量量子化:\( L = m_e v r = n \hbar \)
- 库仑力提供向心力:\( \dfrac{e^2}{4\pi\varepsilon_0 r^2} = \dfrac{m_e v^2}{r} \)
2. 推导出轨道半径:
\[ r_n = \frac{4\pi\varepsilon_0 \hbar^2}{m_e e^2} n^2 = a_0 n^2 \] 其中玻尔半径: \[ a_0 = \frac{4\pi\varepsilon_0 \hbar^2}{m_e e^2} \]
3. 推导出能量:
总能量 = 动能 + 势能: \[ E = \frac{1}{2} m_e v^2 - \frac{1}{4\pi\varepsilon_0} \frac{e^2}{r} \]
由向心力关系:\( \frac{1}{2} m_e v^2 = \frac{1}{2} \cdot \frac{1}{4\pi\varepsilon_0} \frac{e^2}{r} \)
所以: \[ E = \frac{1}{2} \cdot \frac{e^2}{4\pi\varepsilon_0 r} - \frac{e^2}{4\pi\varepsilon_0 r} = -\frac{1}{2} \cdot \frac{e^2}{4\pi\varepsilon_0 r} \]
代入 \( r = a_0 n^2 \),得: \[ E_n = -\frac{1}{2} \cdot \frac{e^2}{4\pi\varepsilon_0} \cdot \frac{1}{a_0 n^2} \]
现在我们来处理这个表达式。
二、将 \( E_n \) 表达成 \( \dfrac{2m_e E a_0^2}{\hbar^2} \) 的形式
我们计算: \[ \frac{2m_e E_n a_0^2}{\hbar^2} \]
先代入 \( E_n = -\dfrac{1}{2} \cdot \dfrac{e^2}{4\pi\varepsilon_0} \cdot \dfrac{1}{a_0 n^2} \)
所以: \[ \frac{2m_e}{\hbar^2} \cdot \left( -\frac{1}{2} \cdot \frac{e^2}{4\pi\varepsilon_0} \cdot \frac{1}{a_0 n^2} \right) \cdot a_0^2 = -\frac{m_e}{\hbar^2} \cdot \frac{e^2}{4\pi\varepsilon_0} \cdot \frac{a_0}{n^2} \]
现在代入 \( a_0 = \dfrac{4\pi\varepsilon_0 \hbar^2}{m_e e^2} \),得:
\[ \frac{a_0}{\hbar^2} = \frac{4\pi\varepsilon_0}{m_e e^2} \quad \Rightarrow \quad \frac{m_e}{\hbar^2} \cdot \frac{e^2}{4\pi\varepsilon_0} \cdot a_0 = \frac{m_e}{\hbar^2} \cdot \frac{e^2}{4\pi\varepsilon_0} \cdot \frac{4\pi\varepsilon_0 \hbar^2}{m_e e^2} = 1 \]
所以: \[ \frac{2m_e E_n a_0^2}{\hbar^2} = -\frac{1}{n^2} \]
✅ 完全正确!
三、结论:玻尔模型中确实有
\[ \boxed{ \frac{2m_e E_n a_0^2}{\hbar^2} = -\frac{1}{n^2} } \]
这个关系可以从玻尔的原始推导中得出,虽然他本人可能没有这样写,但它隐含在他对 \(E_n\) 和 \(a_0\) 的表达式中。
四、物理意义解读
这个公式可以理解为:
- 左边是一个无量纲能量表达式:\(\dfrac{2m E a_0^2}{\hbar^2}\)
- 它等于 \(-1/n^2\),说明能量量子化表现为一个简单的倒平方律
这与薛定谔方程解中的形式惊人地一致!
在量子力学中,我们也会遇到类似结构,例如在径向方程中定义: \[ \kappa = \frac{\sqrt{-2mE}}{\hbar}, \quad \rho = \kappa r \] 而 \(a_0\) 是自然长度尺度,所以 \(\kappa a_0 \propto 1/n\),平方后就是 \(1/n^2\)。
五、与量子力学的联系
在薛定谔方程求解中,我们定义: \[ \lambda = \frac{\mu k}{\hbar^2 \kappa}, \quad k = \frac{e^2}{4\pi\varepsilon_0}, \quad \kappa = \frac{\sqrt{-2\mu E}}{\hbar} \] 最终得到 \(\lambda = n\),并导出 \(E_n \propto -1/n^2\)
而如果你代入 \(a_0\) 的定义,你会发现: \[ \frac{2m_e |E| a_0^2}{\hbar^2} = \frac{1}{n^2} \] 正是你提出的公式的绝对值形式。
👉 所以:这个关系在玻尔模型和量子力学中都成立,是氢原子能量结构的一个“不变量”表达。
六、量纲验证(快速检查)
我们验证左边是否无量纲:
- \(m_e\): \([M]\)
- \(E\): \([M L^2 / T^2]\)
- \(a_0^2\): \([L^2]\)
- \(\hbar^2\): \([M^2 L^4 / T^2]\)
所以: \[ \frac{m E a^2}{\hbar^2} \sim \frac{[M] [M L^2 / T^2] [L^2]}{[M^2 L^4 / T^2]} = \frac{M^2 L^4 / T^2}{M^2 L^4 / T^2} = 1 \]
✅ 无量纲,合理。
✅ 最终结论
是的,玻尔在推导氢原子光谱时,虽然没有明确写出 \[ \frac{2m_e E a_0^2}{\hbar^2} = -\frac{1}{n^2} \] 但这个关系完全可以从他的结果中推导出来,是其理论的自然推论。
它体现了:
- 玻尔模型的数学自洽性
- \(a_0\) 作为自然长度尺度的核心地位
- 能量量子化 \(E_n \propto -1/n^2\) 的简洁表达
这个公式也在量子力学中继续成立,说明玻尔模型虽然被取代,但其核心洞察——用基本常数构造出量子化的能量结构——是深刻而持久的。
📌 延伸思考
这个表达式也可以写成: \[ \sqrt{-\frac{2m_e E}{\hbar^2}} \cdot a_0 = \frac{1}{n} \] 左边是“衰减波数”乘以玻尔半径,等于 \(1/n\),这在量子力学中对应归一化条件。
👉 它暗示了:即使是旧量子论,也已经触碰到波动力学的边缘。
1. 角向方程(球谐函数方程)
方程形式: \[ \frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial Y}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2 Y}{\partial\phi^2} = -\ell(\ell+1)Y \]
物理意义
-
轨道角动量量子化
角向方程的解是 球谐函数 \(Y_\ell^m(\theta,\phi)\),其本征值 \(\ell(\ell+1)\hbar^2\) 对应电子轨道角动量平方 \(\hat{L}^2\) 的量子化:
\[ \hat{L}^2 Y_\ell^m = \ell(\ell+1)\hbar^2 Y_\ell^m
\]
其中:- \(\ell\) 为 角量子数(\(\ell = 0,1,2,\dots\)),决定轨道形状(s, p, d, …轨道)。
- \(m\) 为 磁量子数(\(m = -\ell, \dots, \ell\)),表示角动量在z轴的分量 \(L_z = m\hbar\)。
-
空间取向的量子化
球谐函数 \(Y_\ell^m(\theta,\phi)\) 的复数相位因子 \(e^{im\phi}\) 反映了电子波函数的方位角依赖性,说明轨道角动量在空间中的取向是量子化的(空间量子化)。 -
对称性与电子云分布
球谐函数的模平方 \(|Y_\ell^m|^2\) 描述电子在角度 \((\theta,\phi)\) 方向的概率分布。例如:- \(\ell=0\)(s轨道):球对称分布。
- \(\ell=1\)(p轨道):哑铃形分布。
2. 径向方程
方程形式: \[ -\frac{\hbar^2}{2\mu}\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{dR}{dr}\right) + \left[ \frac{\hbar^2\ell(\ell+1)}{2\mu r^2} - \frac{e^2}{4\pi\epsilon_0 r} \right]R = ER \]
物理意义
-
电子在径向的势能运动
径向方程描述电子在 有效势能场 中的运动:
\[ V_\text{eff}(r) = \underbrace{-\frac{e^2}{4\pi\epsilon_0 r}}_{\text{库仑势}} + \underbrace{\frac{\hbar^2\ell(\ell+1)}{2\mu r^2}}_{\text{离心势垒}}
\]- 库仑势:质子与电子的吸引势能(负值)。
- 离心势垒:角动量引起的排斥势能(\(\ell\) 越大,电子越难靠近原子核)。
-
能量量子化与主量子数 \(n\)
径向方程的解 \(R_{n\ell}(r)\) 要求能量 \(E\) 必须满足量子化条件:
\[ E_n = -\frac{\mu e^4}{2(4\pi\epsilon_0)^2\hbar^2 n^2} \quad (n=1,2,3,\dots)
\]- \(n\) 为 主量子数,决定能级(如 \(n=1\) 对应基态)。
- \(n \geq \ell +1\),即径向波函数的节点数为 \(n-\ell-1\)。
-
电子径向概率分布
径向概率密度 \(P(r) = r^2 |R_{n\ell}(r)|^2\) 表示电子在距离核 \(r\) 处出现的概率。例如:- 基态(\(n=1,\ell=0\)):概率峰值在玻尔半径 \(a_0\) 处。
- 激发态(如 \(n=2,\ell=0\)):存在径向节点。
3. 径向与角向的协同作用
-
完整波函数的构建
总波函数 \(\psi_{n\ell m}(r,\theta,\phi) = R_{n\ell}(r) Y_\ell^m(\theta,\phi)\) 结合了径向和角向信息,完整描述电子的量子态。 -
简并度
对于给定的 \(n\),不同的 \((\ell,m)\) 组合(共 \(n^2\) 种)对应同一能级 \(E_n\),称为 简并。例如:- \(n=2\) 时,\(\ell=0\)(1种)和 \(\ell=1\)(3种,对应 \(m=-1,0,1\))共4个简并态。
-
物理可观测量的分离
- 角向部分决定角动量及其取向。
- 径向部分决定能量和电子与核的平均距离。
4. 直观对比
特征 | 角向方程 | 径向方程 |
---|---|---|
描述对象 | 电子运动的角度分布 \((\theta,\phi)\) | 电子运动的径向分布 \(r\) |
量子数 | \(\ell, m\) | \(n, \ell\) |
核心物理量 | 轨道角动量 \(\hat{L}^2, L_z\) | 能量 \(E\)、电子-核距离 |
势能项 | 无显式势能 | 库仑势 + 离心势垒 |
解的函数 | 球谐函数 \(Y_\ell^m(\theta,\phi)\) | 拉盖尔多项式 \(R_{n\ell}(r)\) |
对称性 | 旋转对称性(SO(3)群) | 径向对称性(依赖 \(n,\ell\)) |
5. 为什么需要分离变量?
- 数学简化:将三维偏微分方程分解为独立的常微分方程。
- 物理清晰性:角向部分揭示角动量特性,径向部分揭示能级和束缚态性质。
- 实验对应:角向解解释原子光谱的偏振特性(如塞曼效应),径向解解释能级间隔(如巴尔末系)。
通过这种分离,氢原子的量子力学解完美复现了实验观测到的光谱规律(如里德伯公式),并奠定了多电子原子结构理论的基础。