第VI章：基本应用——精确解的凯旋

2026-05-16 科普 9.9k 字 49 分钟阅读
📑 目录
  1. 前置知识：特殊函数基础
    1. 1. Hermite多项式：谐振子的伙伴
    2. 2. Legendre多项式与球谐函数：角动量的几何
    3. 3. Laguerre多项式：氢原子的径向结构
  2. 钢琴师的量子调律
  3. 一、谐振子：代数解法的杰作
    1. 1.1 问题设定
    2. 1.2 升降算符的引入
    3. 1.3 Dirac的代数解法：与标准教材的区别
    4. 1.4 用升降算符求解
    5. 1.5 零点能与不确定性
  4. 二、角动量：本征值问题的几何
    1. 2.1 角动量算符的定义
    2. 2.2 角动量的代数结构
    3. 2.3 Dirac的升降算符方法
    4. 2.4 球谐函数
  5. 三、电子自旋：内禀的角动量
    1. 3.1 自旋的发现
    2. 3.2 自旋算符
    3. 3.3 泡利矩阵
    4. 3.4 自旋的物理效应
  6. 四、中心力场与氢原子
    1. 4.1 中心力场的对称性
    2. 4.2 分离变量
    3. 4.3 氢原子的能级
    4. 4.4 径向波函数与特殊函数
  7. 五、选择定则：跃迁的密码
    1. 5.1 光与物质的相互作用
    2. 5.2 选择定则的推导
    3. 5.3 自发辐射与受激辐射
  8. 六、塞曼效应：磁场中的分裂
    1. 6.1 磁矩与角动量
    2. 6.2 正常塞曼效应
    3. 6.3 反常塞曼效应
    4. 6.4 帕邢-巴克效应
  9. 数值例子
    1. 例1：谐振子能级与波函数的具体数值
    2. 例2：氢原子1s态和2p态的数值计算
  10. 本章总结
  11. 练习与思考

第VI章：基本应用——精确解的凯旋

当抽象的代数落在具体的物理系统上，算符的力量才真正显现。谐振子歌唱，角动量起舞，氢原子揭示了宇宙的谱线——这是精确解的凯旋。

前置知识：特殊函数基础

量子力学的精确解往往与经典数学中的特殊函数密切相关。在深入学习谐振子、角动量和氢原子之前，我们需要掌握三类核心特殊函数：Hermite多项式、Legendre多项式（连带）和Laguerre多项式（连带）。

1. Hermite多项式：谐振子的伙伴

Hermite多项式 Hn(x)H_n(x) 定义在区间 (−∞,∞)(-\infty, \infty) 上，权函数为 e−x2e^{-x^2}，满足正交归一条件：

∫−∞∞Hm(x)Hn(x)e−x2dx=π2nn!δmn\int_{-\infty}^{\infty} H_m(x) H_n(x) e^{-x^2} dx = \sqrt{\pi} 2^n n! \delta_{mn}

Rodrigues公式：

Hn(x)=(−1)nex2dndxne−x2H_n(x) = (-1)^n e^{x^2} \frac{d^n}{dx^n} e^{-x^2}

前几项：

  • H0(x)=1H_0(x) = 1
  • H1(x)=2xH_1(x) = 2x
  • H2(x)=4x2−2H_2(x) = 4x^2 - 2
  • H3(x)=8x3−12xH_3(x) = 8x^3 - 12x
  • H4(x)=16x4−48x2+12H_4(x) = 16x^4 - 48x^2 + 12

递推关系：

Hn+1(x)=2xHn(x)−2nHn−1(x)H_{n+1}(x) = 2x H_n(x) - 2n H_{n-1}(x)

Hn′(x)=2nHn−1(x)H_n'(x) = 2n H_{n-1}(x)

与谐振子的联系：一维谐振子的定态波函数为

ψn(q)=(mωπℏ)1/412nn!Hn(mωℏq)exp(−mωq22ℏ)\psi_n(q) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4} \frac{1}{\sqrt{2^n n!}} H_n\left(\sqrt{\frac{m\omega}{\hbar}} q\right) \exp\left(-\frac{m\omega q^2}{2\hbar}\right)

Hermite多项式的零点数目等于量子数 nn，而 e−x2e^{-x^2} 的高斯衰减确保了波函数可归一化。

2. Legendre多项式与球谐函数：角动量的几何

Legendre多项式 Pl(x)P_l(x) 定义在 [−1,1][-1, 1] 上，满足：

∫−11Pl(x)Pl′(x)dx=22l+1δll′\int_{-1}^{1} P_l(x) P_{l'}(x) dx = \frac{2}{2l+1} \delta_{ll'}

Rodrigues公式：

Pl(x)=12ll!dldxl(x2−1)lP_l(x) = \frac{1}{2^l l!} \frac{d^l}{dx^l}(x^2 - 1)^l

前几项：

  • P0(x)=1P_0(x) = 1
  • P1(x)=xP_1(x) = x
  • P2(x)=12(3x2−1)P_2(x) = \frac{1}{2}(3x^2 - 1)
  • P3(x)=12(5x3−3x)P_3(x) = \frac{1}{2}(5x^3 - 3x)

连带Legendre函数：

Plm(x)=(1−x2)m/2dmdxmPl(x),m=0,1,2,…,lP_l^m(x) = (1-x^2)^{m/2} \frac{d^m}{dx^m} P_l(x), \quad m = 0, 1, 2, \ldots, l

球谐函数是角动量本征态的完整表示：

Ylm(θ,ϕ)=(2l+1)4π(l−m)!(l+m)!Plm(cosθ)eimϕY_l^m(\theta, \phi) = \sqrt{\frac{(2l+1)}{4\pi} \frac{(l-m)!}{(l+m)!}} P_l^m(\cos\theta) e^{im\phi}

正交归一性：

∫02πdϕ∫0πdθsinθ(Ylm)∗Yl′m′=δll′δmm′\int_0^{2\pi} d\phi \int_0^{\pi} d\theta \sin\theta (Y_l^m)^* Y_{l'}^{m'} = \delta_{ll'} \delta_{mm'}

3. Laguerre多项式：氢原子的径向结构

Laguerre多项式 Ln(x)L_n(x) 和连带Laguerre多项式 Lnk(x)L_n^k(x) 定义在 [0,∞)[0, \infty) 上，权函数为 e−xe^{-x}。

Rodrigues公式：

Ln(x)=exdndxn(xne−x)L_n(x) = e^x \frac{d^n}{dx^n}(x^n e^{-x})

Lnk(x)=(−1)kdkdxkLn+k(x)L_n^k(x) = (-1)^k \frac{d^k}{dx^k} L_{n+k}(x)

前几项：

  • L0(x)=1L_0(x) = 1
  • L1(x)=−x+1L_1(x) = -x + 1
  • L2(x)=x2−4x+2L_2(x) = x^2 - 4x + 2

与氢原子的联系：氢原子径向波函数为

Rnl(r)=(2na0)3(n−l−1)!2n[(n+l)!](2rna0)lLn−l−12l+1(2rna0)e−r/(na0)R_{nl}(r) = \sqrt{\left(\frac{2}{na_0}\right)^3 \frac{(n-l-1)!}{2n[(n+l)!]}} \left(\frac{2r}{na_0}\right)^l L_{n-l-1}^{2l+1}\left(\frac{2r}{na_0}\right) e^{-r/(na_0)}

连带Laguerre多项式的阶数 n−l−1n-l-1 等于径向波函数的节点数（不包括原点）。

graph TD
    A["量子系统"] --> B["谐振子"]
    A --> C["角动量"]
    A --> D["氢原子"]
    
    B --> E["Hermite多项式 H_n(x)"]
    C --> F["Legendre多项式 P_l(x)"]
    C --> G["连带Legendre P_l^m(x)"]
    C --> H["球谐函数 Y_l^m"]
    D --> I["连带Laguerre L_{n-l-1}^{2l+1}(x)"]
    
    E --> J["权函数: e^{-x²}, 区间: (-∞,∞)"]
    F --> K["权函数: 1, 区间: [-1,1]"]
    I --> L["权函数: e^{-x}, 区间: [0,∞)"]
    
    style A fill:#e3f2fd
    style E fill:#bbdefb
    style H fill:#90caf9
    style I fill:#64b5f6

钢琴师的量子调律

2089年，"谐波之城"的中央音乐厅里，首席调律师陈默面对一架特殊的钢琴犯了难。这不是普通的乐器——它的每根弦都只有一个原子那么粗，振动遵循量子力学而非经典声学。

陈默发现，这架"量子钢琴"的音符不是连续的：你只能弹出 do、re、mi 的某些特定频率，中间不存在过渡音。更奇怪的是，每个音符内部还有"谐波"——当你弹出基音 do 时，会同时激发出 do 的高八度、高两个八度……但永远不会出现 fa 或 sol。

音乐厅的AI助手解释道："这是谐振子的量子化。能量不是连续的，而是 En=ℏω(n+12)E_n = \hbar\omega(n+\frac{1}{2})。基态有零点能 12ℏω\frac{1}{2}\hbar\omega，激发态是等间距的能级。"

"那角动量呢？"陈默指着旋转的舞台灯光系统问道。

"角动量也是量子化的，"AI回答，"而且更奇妙——你只能取到 \\sqrt{l(l+1)}\hbar 的大小，z分量只能是整数步长。电子自旋甚至不需要任何轨道运动，它自带 12ℏ\frac{1}{2}\hbar 的内禀角动量。"

陈默若有所思："所以氢原子的光谱……"

"正是谐振子、角动量和库仑势三者共同谱写的乐章。"AI投影出一组优美的谱线，"这就是精确解的凯旋——当数学遇见物理，宇宙开始歌唱。"

一、谐振子：代数解法的杰作

1.1 问题设定

一维量子谐振子的哈密顿量为：

H^=p^22m+12mω2q^2\hat{H} = \frac{\hat{p}^2}{2m} + \frac{1}{2}m\omega^2\hat{q}^2

经典谐振子是物理学的"Rosetta Stone"——几乎所有周期性系统在其平衡位置附近都可以近似为谐振子。

1.2 升降算符的引入

狄拉克的洞见：引入无量纲的升降算符（湮灭与产生算符）：

a^=mω2ℏ(q^+ip^mω)\hat{a} = \sqrt{\frac{m\omega}{2\hbar}}\left(\hat{q} + \frac{i\hat{p}}{m\omega}\right)

a^†=mω2ℏ(q^−ip^mω)\hat{a}^\dagger = \sqrt{\frac{m\omega}{2\hbar}}\left(\hat{q} - \frac{i\hat{p}}{m\omega}\right)

这两个算符不是厄米的（a^≠a^†\hat{a} \neq \hat{a}^\dagger），但它们的对易关系极其简洁：

[a^,a^†]=1[\hat{a}, \hat{a}^\dagger] = 1

graph TD
    A["位置 q̂ 和动量 p̂"] -->|"线性组合"| B["升降算符 â, â†"]
    B --> C["对易关系 [â,â†] = 1"]
    C --> D["哈密顿量简化:
Ĥ = ℏω(â†â + 1/2)"]
    D --> E["粒子数算符 n̂ = â†â"]
    E --> F["本征值 n = 0,1,2,..."]
    
    style A fill:#e3f2fd
    style B fill:#bbdefb
    style C fill:#90caf9
    style F fill:#64b5f6

1.3 Dirac的代数解法：与标准教材的区别

狄拉克在《量子力学原理》中对谐振子的处理有几个独特之处，与现代标准教材（如Griffiths或Shankar）有所区别：

Dirac的原始方法：

  1. 从哈密顿量直接构造：Dirac不先求解微分方程再回头定义升降算符，而是直接从哈密顿量的代数结构出发，识别出可以因式分解的形式
  2. 强调变换理论：Dirac将 a^\hat{a} 和 a^†\hat{a}^\dagger 视为从一个表象到另一个表象的变换生成元——这不是技巧，而是他整个量子力学框架的自然应用
  3. 不依赖具体表象：Dirac的推导完全不依赖位置表象或动量表象，是纯抽象的代数操作。本征值 En=ℏω(n+12)E_n = \hbar\omega(n+\frac{1}{2}) 的得出不需要写出任何微分方程
  4. 归一化作为构造的一部分：Dirac在定义升降算符的作用时，将归一化常数作为代数结构的自然产出，而非事后调整

标准教材的常见路径：

  1. 先写出位置表象的薛定谔方程：−ℏ22md2ψdq2+12mω2q2ψ=Eψ-\frac{\hbar^2}{2m}\frac{d^2\psi}{dq^2} + \frac{1}{2}m\omega^2 q^2\psi = E\psi
  2. 通过变量代换化为无量纲形式，发现方程的级数解必须在某处截断
  3. 截断条件给出 En=ℏω(n+12)E_n = \hbar\omega(n+\frac{1}{2})
  4. 然后引入升降算符作为"更优雅的方法"

Dirac的路径是从代数到分析，而标准教材是从分析到代数。Dirac的方法展示了量子力学的真正力量：如果你正确识别了算符的代数结构，连微分方程都不需要解。

1.4 用升降算符求解

哈密顿量用升降算符表示为：

H^=ℏω(a^†a^+12)=ℏω(n^+12)\hat{H} = \hbar\omega\left(\hat{a}^\dagger\hat{a} + \frac{1}{2}\right) = \hbar\omega\left(\hat{n} + \frac{1}{2}\right)

其中 n^=a^†a^\hat{n} = \hat{a}^\dagger\hat{a} 是粒子数算符。

关键性质：

  • 基态条件：若 ∣0⟩|0\rangle 是基态，则 a^∣0⟩=0\hat{a}|0\rangle = 0（否则能量可以更低）
  • 能级：En=ℏω(n+12)E_n = \hbar\omega(n + \frac{1}{2})
  • 等间距：ΔE=En+1−En=ℏω\Delta E = E_{n+1} - E_n = \hbar\omega
graph TD
    subgraph "能级图 ("ℏω为单位")"
        E0["n=0: E = 1/2 (零点能)"] <--> E1["n=1: E = 3/2"]
        E1 <--> E2["n=2: E = 5/2"]
        E2 <--> E3["n=3: E = 7/2"]
        E3 <--> E4["..."]
    end
    
    E0 -->|"#quot;â†"|0⟩ = |1⟩"| E1
    E1 -->|"#quot;â†"|1⟩ = √2|2⟩"| E2
    E2 -->|"#quot;â†"|2⟩ = √3|3⟩"| E3
    
    E1 -->|"#quot;â"|1⟩ = |0⟩"| E0
    E2 -->|"#quot;â"|2⟩ = √2|1⟩"| E1
    E3 -->|"#quot;â"|3⟩ = √3|2⟩"| E2
    
    style E0 fill:#e8f5e9
    style E1 fill:#c8e6c9
    style E2 fill:#a5d6a7
    style E3 fill:#81c784

1.5 零点能与不确定性

基态能量 E0=12ℏωE_0 = \frac{1}{2}\hbar\omega 不为零——这是零点能，量子涨落的最低限度。

基态波函数：

ψ0(q)=(mωπℏ)1/4exp(−mωq22ℏ)\psi_0(q) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\left(-\frac{m\omega q^2}{2\hbar}\right)

这是一个高斯波包，位置和动量的不确定性满足最小值：

Δq=ℏ2mω,Δp=mℏω2,Δq⋅Δp=ℏ2\Delta q = \sqrt{\frac{\hbar}{2m\omega}}, \quad \Delta p = \sqrt{\frac{m\hbar\omega}{2}}, \quad \Delta q \cdot \Delta p = \frac{\hbar}{2}

谐振子基态是最小不确定态，同时具有确定的零点能和不确定性的最小代价。

二、角动量：本征值问题的几何

2.1 角动量算符的定义

经典角动量 L⃗=r⃗×p⃗\vec{L} = \vec{r} \times \vec{p} 的量子对应是：

L^x=y^p^z−z^p^y\hat{L}_x = \hat{y}\hat{p}_z - \hat{z}\hat{p}_y

L^y=z^p^x−x^p^z\hat{L}_y = \hat{z}\hat{p}_x - \hat{x}\hat{p}_z

L^z=x^p^y−y^p^x\hat{L}_z = \hat{x}\hat{p}_y - \hat{y}\hat{p}_x

2.2 角动量的代数结构

角动量算符满足特殊的对易关系：

[L^x,L^y]=iℏL^z,[L^y,L^z]=iℏL^x,[L^z,L^x]=iℏL^y[\hat{L}_x, \hat{L}_y] = i\hbar\hat{L}_z, \quad [\hat{L}_y, \hat{L}_z] = i\hbar\hat{L}_x, \quad [\hat{L}_z, \hat{L}_x] = i\hbar\hat{L}_y

这个结构可以紧凑地写成：

[L^i,L^j]=iℏϵijkL^k[\hat{L}_i, \hat{L}_j] = i\hbar\epsilon_{ijk}\hat{L}_k

其中 ϵijk\epsilon_{ijk} 是Levi-Civita符号（完全反对称张量）。

graph TD
    A["角动量对易关系
[Lᵢ,Lⱼ] = iℏεᵢⱼₖLₖ"] --> B["总角动量平方 L̂² = L̂ₓ² + L̂ᵧ² + L̂ᵤ²"]
    B --> C["[L̂²,L̂ᵢ] = 0: 可以同时测量"]
    C --> D["选择 L̂² 和 L̂ᵤ 的共同本征态"]
    
    D --> E["本征值方程:
L̂²|l,m⟩ = ℏ²l(l+1)|l,m⟩"]
    D --> F["本征值方程:
L̂ᵤ|l,m⟩ = ℏm|l,m⟩"]
    
    E --> G["l = 0, 1/2, 1, 3/2, 2, ..."]
    F --> H["m = -l, -l+1, ..., l-1, l"]
    
    G --> I["每个 l 有 (2l+1) 个 m 值"]
    H --> I
    
    style A fill:#ffebee
    style B fill:#ffcdd2
    style E fill:#ef9a9a
    style F fill:#ef9a9a
    style I fill:#e57373

2.3 Dirac的升降算符方法

狄拉克对角动量本征值问题的处理同样展示了他偏爱的代数方法。定义角动量升降算符：

L^+=L^x+iL^y,L^−=L^x−iL^y\hat{L}_+ = \hat{L}_x + i\hat{L}_y, \quad \hat{L}_- = \hat{L}_x - i\hat{L}_y

它们满足：

[L^z,L^±]=±ℏL^±[\hat{L}_z, \hat{L}_\pm] = \pm\hbar\hat{L}_\pm

这意味着 L^±∣l,m⟩∝∣l,m±1⟩\hat{L}_\pm|l,m\rangle \propto |l, m\pm 1\rangle——升降算符将 mm 值增加或减少一个单位。

通过边界条件（不能超出 m=±lm = \pm l 的范围），得到：

  • ll 必须是整数或半整数
  • mm 取值范围：m=−l,−l+1,…,l−1,lm = -l, -l+1, \ldots, l-1, l（共 2l+12l+1 个值）

Dirac方法的独特之处：与谐振子一样，Dirac不依赖任何具体表象（如位置表象的微分方程），直接从对易关系的代数结构推导出全部本征值谱。ll 和 mm 的量子化是代数约束的自然结果，不是人为施加的边界条件。

2.4 球谐函数

在位置表象中，角动量本征态是球谐函数 Ylm(θ,ϕ)Y_l^m(\theta, \phi)：

L^2Ylm=ℏ2l(l+1)Ylm\hat{L}^2 Y_l^m = \hbar^2 l(l+1) Y_l^m

L^zYlm=ℏmYlm\hat{L}_z Y_l^m = \hbar m Y_l^m

前几例：

  • l=0l=0：Y00=14πY_0^0 = \frac{1}{\sqrt{4\pi}}（s波，球对称）
  • l=1l=1：Y10=34πcosθY_1^0 = \sqrt{\frac{3}{4\pi}}\cos\theta，Y1±1=∓38πsinθe±iϕY_1^{\pm 1} = \mp\sqrt{\frac{3}{8\pi}}\sin\theta e^{\pm i\phi}（p波）
  • l=2l=2：d波，5个分量
graph LR
    subgraph "s波 l=0"
        S0["Y₀⁰: 球对称"]
    end
    
    subgraph "p波 l=1"
        P0["Y₁⁰: 沿z轴"]
        P1["Y₁±¹: 在xy平面旋转"]
    end
    
    subgraph "d波 l=2"
        D0["Y₂⁰: 沿z轴椭圆"]
        D1["Y₂±¹: 复杂形状"]
        D2["Y₂±²: 四叶形"]
    end
    
    S0 -->|"激发"| P0
    P0 -->|"激发"| D0
    
    style S0 fill:#e8f5e9
    style P0 fill:#c8e6c9
    style D0 fill:#a5d6a7

三、电子自旋：内禀的角动量

3.1 自旋的发现

1922年施特恩-盖拉赫实验揭示了电子具有内禀角动量——不需要任何轨道运动，电子本身就携带 12ℏ\frac{1}{2}\hbar 的角动量。

自旋是量子力学最深刻的概念之一：它纯粹是量子效应，没有经典对应。一个"点粒子"怎么可能"自转"？答案是：自旋不是真实的转动，而是量子态的内禀标签。

3.2 自旋算符

自旋算符 S^x,S^y,S^z\hat{S}_x, \hat{S}_y, \hat{S}_z 满足与角动量相同的对易关系：

[S^i,S^j]=iℏϵijkS^k[\hat{S}_i, \hat{S}_j] = i\hbar\epsilon_{ijk}\hat{S}_k

对于电子，s=12s = \frac{1}{2}，因此：

  • ms=±12m_s = \pm\frac{1}{2}（自旋向上/向下）
  • S^2\hat{S}^2 的本征值：ℏ2s(s+1)=34ℏ2\hbar^2 s(s+1) = \frac{3}{4}\hbar^2
  • S^z\hat{S}_z 的本征值：±ℏ2\pm\frac{\hbar}{2}

3.3 泡利矩阵

自旋算符可以用 2×22\times 2 的泡利矩阵表示：

S^i=ℏ2σi\hat{S}_i = \frac{\hbar}{2}\sigma_i

其中：

σx=(0110),σy=(0−ii0),σz=(100−1)\sigma_x = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \quad \sigma_y = \begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}, \quad \sigma_z = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}

graph TD
    A["自旋 s = 1/2"] --> B["两个态: |↑⟩ = (1,0)ᵀ, |↓⟩ = (0,1)ᵀ"]
    B --> C["泡利矩阵 σₓ, σᵧ, σᵤ"]
    C --> D["Ŝᵢ = (ℏ/2)σᵢ"]
    
    D --> E["[σᵢ,σⱼ] = 2iεᵢⱼₖσₖ"]
    E --> F["反对易: {σᵢ,σⱼ} = 2δᵢⱼI"]
    
    F --> G["σₓσᵧ = iσᵤ (循环)"]
    G --> H["泡利矩阵的平方 = I"]
    
    C --> I["物理意义:"]
    I --> I1["σᵤ: 测量自旋z分量"]
    I --> I2["σₓ: 测量自旋x分量"]
    I --> I3["σᵧ: 测量自旋y分量"]
    
    style A fill:#fff3e0
    style B fill:#ffe0b2
    style C fill:#ffcc80
    style F fill:#ffb74d

3.4 自旋的物理效应

自旋带来一系列可观测效应：

  • 塞曼效应：自旋磁矩与外磁场的相互作用导致能级分裂
  • 泡利不相容原理：费米子（半整数自旋）不能占据相同的量子态
  • 自旋-轨道耦合：自旋磁矩与电子轨道运动产生的磁场相互作用

四、中心力场与氢原子

4.1 中心力场的对称性

中心力场 V(r)V(r) 只依赖于径向距离，与角度无关。哈密顿量：

H^=−ℏ22m∇2+V(r)\hat{H} = -\frac{\hbar^2}{2m}\nabla^2 + V(r)

由于势的球对称性，[H^,L^2]=0[\hat{H}, \hat{L}^2] = 0 和 [H^,L^z]=0[\hat{H}, \hat{L}_z] = 0——能量、总角动量和z分量角动量可以同时确定。

4.2 分离变量

在球坐标中，波函数分离为径向和角向部分：

ψ(r,θ,ϕ)=Rnl(r)Ylm(θ,ϕ)\psi(r,\theta,\phi) = R_{nl}(r)Y_l^m(\theta,\phi)

径向方程：

[−ℏ22m1r2ddr(r2ddr)+ℏ2l(l+1)2mr2+V(r)]Rnl=EnlRnl\left[-\frac{\hbar^2}{2m}\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{d}{dr}\right) + \frac{\hbar^2 l(l+1)}{2mr^2} + V(r)\right]R_{nl} = E_{nl}R_{nl}

其中 ℏ2l(l+1)2mr2\frac{\hbar^2 l(l+1)}{2mr^2} 是离心势能——角动量的运动学效应。

graph TD
    A["中心力场问题"] --> B["球对称 V(r)"]
    B --> C["守恒量: Ĥ, L̂², L̂ᵤ"]
    C --> D["分离变量: ψ = Rₙₗ(r)Yₗᵐ(θ,φ)"]
    
    D --> E["径向方程:
含离心势 ℏ²l(l+1)/(2mr²)"]
    D --> F["角向方程:
球谐函数 Yₗᵐ"]
    
    E --> G["量子数 n, l"]
    F --> H["量子数 l, m"]
    
    G --> I["主量子数 n = 1,2,3,..."]
    G --> J["角量子数 l = 0,1,...,n-1"]
    H --> K["磁量子数 m = -l,...,l"]
    
    style A fill:#e3f2fd
    style C fill:#bbdefb
    style E fill:#90caf9
    style F fill:#90caf9

4.3 氢原子的能级

对于库仑势 V(r)=−e24πϵ0rV(r) = -\frac{e^2}{4\pi\epsilon_0 r}，能级奇迹般地简并：

En=−mee42(4πϵ0)2ℏ21n2=−13.6 eVn2E_n = -\frac{m_e e^4}{2(4\pi\epsilon_0)^2\hbar^2}\frac{1}{n^2} = -\frac{13.6\text{ eV}}{n^2}

其中 n=1,2,3,…n = 1, 2, 3, \ldots 是主量子数。

意外之处：能量只依赖于 nn，而不依赖于 ll。这意味着 2s2s 和 2p2p（以及 3s,3p,3d3s, 3p, 3d 等）具有相同的能量——这是库仑势特有的"偶然简并"，源于 1/r1/r 势的额外对称性（拉普拉斯-龙格-楞次矢量守恒）。

graph TD
    subgraph "氢原子能级图"
        N1["n=1: -13.6 eV (基态)
        l=0 (1s)"] <--> N2["n=2: -3.4 eV
        l=0,1 ("2s,2p")"]
        N2 <--> N3["n=3: -1.51 eV
        l=0,1,2 ("3s,3p,3d")"]
        N3 <--> N4["n=4: -0.85 eV
        ..."]
        N4 <--> N5["n=∞: 0 eV (电离)"]
    end
    
    N2 -->|"莱曼系 (紫外)"| N1
    N3 -->|"巴尔末系 (可见光)"| N2
    N4 -->|"帕邢系 (红外)"| N3
    
    style N1 fill:#ffebee
    style N2 fill:#ffcdd2
    style N3 fill:#ef9a9a
    style N5 fill:#e8f5e9

4.4 径向波函数与特殊函数

氢原子的径向波函数 Rnl(r)R_{nl}(r) 与连带Laguerre多项式相关。前几例：

  • 1s1s (n=1,l=0n=1, l=0)：R10(r)=2(1a0)3/2e−r/a0R_{10}(r) = 2\left(\frac{1}{a_0}\right)^{3/2} e^{-r/a_0}
  • 2s2s (n=2,l=0n=2, l=0)：R20(r)=12(1a0)3/2(1−r2a0)e−r/(2a0)R_{20}(r) = \frac{1}{\sqrt{2}}\left(\frac{1}{a_0}\right)^{3/2} \left(1 - \frac{r}{2a_0}\right)e^{-r/(2a_0)}
  • 2p2p (n=2,l=1n=2, l=1)：R21(r)=126(1a0)3/2ra0e−r/(2a0)R_{21}(r) = \frac{1}{2\sqrt{6}}\left(\frac{1}{a_0}\right)^{3/2} \frac{r}{a_0} e^{-r/(2a_0)}

其中 a_0 = \frac{4\pi\epsilon_0\hbar^2}{m_e e^2} \approx 0.529\text{ Å} 是玻尔半径。

归一化条件：

∫0∞∣Rnl(r)∣2r2dr=1\int_0^\infty |R_{nl}(r)|^2 r^2 dr = 1

注意 r2r^2 因子来自球坐标体积元 dV=r2sinθdrdθdϕdV = r^2 \sin\theta dr d\theta d\phi。

五、选择定则：跃迁的密码

5.1 光与物质的相互作用

原子与电磁场的相互作用导致能级间的跃迁。跃迁概率正比于矩阵元：

⟨n′l′m′∣r⃗^∣nlm⟩\langle n'l'm'|\hat{\vec{r}}|nlm\rangle

5.2 选择定则的推导

利用角动量算符的性质和球谐函数的积分，可以得到选择定则：

Δl=±1,Δm=0,±1\Delta l = \pm 1, \quad \Delta m = 0, \pm 1

这意味着：

  • s波 (l=0l=0) 只能跃迁到 p波 (l=1l=1)
  • p波 (l=1l=1) 可以跃迁到 s波 (l=0l=0) 或 d波 (l=2l=2)
  • 不能从 s波直接跃迁到 d波（Δl=2\Delta l = 2 被禁止）
graph TD
    A["选择定则"] --> B["Δl = ±1"]
    A --> C["Δm = 0, ±1"]
    
    B --> D["电偶极跃迁
角动量守恒"]
    C --> E["光子携带
角动量 ±ℏ或0"]
    
    D --> F["禁戒跃迁:"]
    F --> F1["2s → 1s (Δl=0)"]
    F --> F2["3d → 2s (Δl=2)"]
    
    D --> G["允许跃迁:"]
    G --> G1["2p → 1s (Δl=-1)"]
    G --> G2["3p → 2s (Δl=-1)"]
    G --> G3["3d → 2p (Δl=-1)"]
    
    style A fill:#e8f5e9
    style D fill:#c8e6c9
    style F fill:#ffcdd2
    style G fill:#a5d6a7

5.3 自发辐射与受激辐射

  • 自发辐射：激发态原子自发跃迁到低能级，发射光子
  • 受激辐射：入射光子诱发激发态原子跃迁，发射相干光子——激光的原理
  • 吸收：原子吸收光子从低能级跃迁到高能级

爱因斯坦1917年的系数 A21A_{21}（自发辐射）、B21B_{21}（受激辐射）和 B12B_{12}（吸收）满足关系：

B21=B12,A21B21=8πhν3c3B_{21} = B_{12}, \quad \frac{A_{21}}{B_{21}} = \frac{8\pi h\nu^3}{c^3}

六、塞曼效应：磁场中的分裂

6.1 磁矩与角动量

带电粒子的轨道运动产生磁矩：

μ⃗^L=−e2meL⃗^\hat{\vec{\mu}}_L = -\frac{e}{2m_e}\hat{\vec{L}}

自旋也贡献磁矩：

μ⃗^S=−gse2meS⃗^\hat{\vec{\mu}}_S = -g_s\frac{e}{2m_e}\hat{\vec{S}}

其中 gs≈2g_s \approx 2 是电子的自旋g因子（量子电动力学给出 gs=2.002319...g_s = 2.002319...）。

6.2 正常塞曼效应

在弱磁场中，哈密顿量增加一项：

H^Z=−μ⃗^⋅B⃗=eB2me(L^z+2S^z)\hat{H}_Z = -\hat{\vec{\mu}} \cdot \vec{B} = \frac{eB}{2m_e}(\hat{L}_z + 2\hat{S}_z)

对于自旋为零的原子（单态，S=0S=0），能级分裂为：

ΔE=eB2meℏm=μBBm\Delta E = \frac{eB}{2m_e}\hbar m = \mu_B B m

其中 μB=eℏ2me\mu_B = \frac{e\hbar}{2m_e} 是玻尔磁子。每条能级分裂为 2l+12l+1 条。

graph TD
    subgraph "正常塞曼效应 ("S=0")"
        A0["B=0: 一条能级
        l=1, m=-1,0,1"]
        
        A1["B≠0: 分裂为三条
        m=1 ("ΔE=+μ_BB")"]
        A2["m=0 (ΔE=0)"]
        A3["m=-1 (ΔE=-μ_BB)"]
    end
    
    A0 -->|"磁场开启"| A1
    A0 -->|"磁场开启"| A2
    A0 -->|"磁场开启"| A3
    
    A1 -->|"跃迁"| B["光谱线分裂为三条"]
    A2 -->|"跃迁"| B
    A3 -->|"跃迁"| B
    
    style A0 fill:#e3f2fd
    style A1 fill:#c8e6c9
    style A2 fill:#a5d6a7
    style A3 fill:#c8e6c9

6.3 反常塞曼效应

当自旋不为零时，gs=2g_s = 2 导致能级分裂更复杂。对于一般的总角动量 J⃗=L⃗+S⃗\vec{J} = \vec{L} + \vec{S}，朗德g因子给出：

H^Z=μBBgJJ^z\hat{H}_Z = \mu_B B g_J \hat{J}_z

其中 gJg_J 依赖于 L,S,JL, S, J 的具体组合：

gJ=1+J(J+1)+S(S+1)−L(L+1)2J(J+1)g_J = 1 + \frac{J(J+1) + S(S+1) - L(L+1)}{2J(J+1)}

6.4 帕邢-巴克效应

在极强磁场中，自旋-轨道耦合被磁场压制，系统回到"正常塞曼效应"的行为——这是两个微扰竞争的经典例子。

数值例子

例1：谐振子能级与波函数的具体数值

问题：一个光学谐振腔中的电磁场模式可近似为一维谐振子，频率 ν=5×1014\nu = 5 \times 10^{14} Hz（可见光，绿光）。计算前四个能级的能量（以eV为单位）、相邻能级间距，以及基态波函数在 q=0q = 0 和 q=aq = a 处的比值，其中 a=ℏ/(mω)a = \sqrt{\hbar/(m\omega)} 是特征长度。

解答：

参数计算：

  • ω=2πν=2π×5×1014=3.142×1015\omega = 2\pi\nu = 2\pi \times 5 \times 10^{14} = 3.142 \times 10^{15} rad/s
  • 由于光子没有静止质量，我们考虑等效质量参数 mm 使得 mω2m\omega^2 对应谐振腔的等效弹性系数。实际上对电磁谐振子，ℏω\hbar\omega 直接给出能量量子：

ℏω=1.055×10−34×3.142×1015=3.31×10−19 J\hbar\omega = 1.055 \times 10^{-34} \times 3.142 \times 10^{15} = 3.31 \times 10^{-19}\text{ J}

转换为eV：

ℏω=3.31×10−191.602×10−19≈2.07 eV\hbar\omega = \frac{3.31 \times 10^{-19}}{1.602 \times 10^{-19}} \approx 2.07\text{ eV}

能级计算：

En=ℏω(n+12)E_n = \hbar\omega\left(n + \frac{1}{2}\right)

nnEnE_n (eV)与基态间距 (eV)
01.0350
13.1052.07
25.1754.14
37.2456.21

相邻能级间距：ΔE=En+1−En=ℏω=2.07\Delta E = E_{n+1} - E_n = \hbar\omega = 2.07 eV，对所有 nn 相同。

基态波函数：

ψ0(q)=(mωπℏ)1/4exp(−mωq22ℏ)=1(πa2)1/4exp(−q22a2)\psi_0(q) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4} \exp\left(-\frac{m\omega q^2}{2\hbar}\right) = \frac{1}{(\pi a^2)^{1/4}} \exp\left(-\frac{q^2}{2a^2}\right)

在 q=0q = 0：

ψ0(0)=(mωπℏ)1/4=1(πa2)1/4\psi_0(0) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4} = \frac{1}{(\pi a^2)^{1/4}}

在 q=aq = a：

ψ0(a)=1(πa2)1/4exp(−a22a2)=1(πa2)1/4e−1/2=ψ0(0)×e−0.5\psi_0(a) = \frac{1}{(\pi a^2)^{1/4}} \exp\left(-\frac{a^2}{2a^2}\right) = \frac{1}{(\pi a^2)^{1/4}} e^{-1/2} = \psi_0(0) \times e^{-0.5}

比值：

ψ0(a)ψ0(0)=e−0.5≈0.607\frac{\psi_0(a)}{\psi_0(0)} = e^{-0.5} \approx 0.607

即波函数在特征长度 aa 处衰减到峰值的约61%。概率密度比值为：

∣ψ0(a)∣2∣ψ0(0)∣2=e−1≈0.368\frac{|\psi_0(a)|^2}{|\psi_0(0)|^2} = e^{-1} \approx 0.368

物理讨论：

  • 光学谐振子的能级间距 2.07 eV 对应绿光光子能量——这解释了为什么激光的频率由谐振腔模式决定
  • 基态零点能 1.035 eV 是量子谐振子不可能"静止"的体现
  • 对于典型原子质量尺度（m∼10−25m \sim 10^{-25} kg）的机械谐振子，a∼10−12a \sim 10^{-12} m（皮米），波函数在原子尺度上高度局域

例2：氢原子1s态和2p态的数值计算

问题：计算氢原子基态（1s）电子的：
(a) 电离能（从基态到自由电子所需的能量）
(b) 最概然半径（径向概率密度最大处的 rr）
© 平均半径 ⟨r⟩\langle r \rangle 和 ⟨r2⟩\langle r^2 \rangle
(d) 2p态电子的角动量大小和z分量可能值

解答：

已知参数：

  • a0=0.529a_0 = 0.529 Å =0.529×10−10= 0.529 \times 10^{-10} m =5.29×10−11= 5.29 \times 10^{-11} m
  • En=−13.6/n2E_n = -13.6/n^2 eV

(a) 电离能：

基态能量 E1=−13.6E_1 = -13.6 eV

自由电子能量 E∞=0E_\infty = 0 eV

Eionization=E∞−E1=0−(−13.6)=13.6 eVE_{ionization} = E_\infty - E_1 = 0 - (-13.6) = 13.6\text{ eV}

这就是氢原子的里德伯能量，与实验值 13.598 eV 精确吻合。

(b) 最概然半径：

1s态径向波函数：R10(r)=2(1a0)3/2e−r/a0R_{10}(r) = 2\left(\frac{1}{a_0}\right)^{3/2} e^{-r/a_0}

径向概率密度：P(r)=r2∣R10(r)∣2=4r2a03e−2r/a0P(r) = r^2 |R_{10}(r)|^2 = \frac{4r^2}{a_0^3} e^{-2r/a_0}

求最大值：令 dPdr=0\frac{dP}{dr} = 0

ddr(r2e−2r/a0)=2re−2r/a0−2r2a0e−2r/a0=0\frac{d}{dr}\left(r^2 e^{-2r/a_0}\right) = 2r e^{-2r/a_0} - \frac{2r^2}{a_0}e^{-2r/a_0} = 0

2r(1−ra0)=02r\left(1 - \frac{r}{a_0}\right) = 0

解得 r=0r = 0（极小值）或 r=a0r = a_0（极大值）

最概然半径 = 玻尔半径 a0≈0.529a_0 \approx 0.529 Å

© 平均半径：

⟨r⟩1s=∫0∞rP(r)dr=∫0∞r⋅4r2a03e−2r/a0dr\langle r \rangle_{1s} = \int_0^\infty r P(r) dr = \int_0^\infty r \cdot \frac{4r^2}{a_0^3} e^{-2r/a_0} dr

利用积分公式 ∫0∞xne−axdx=n!an+1\int_0^\infty x^n e^{-ax} dx = \frac{n!}{a^{n+1}}：

⟨r⟩1s=4a03∫0∞r3e−2r/a0dr=4a03⋅3!(2/a0)4=4a03⋅6a0416=24a016=32a0\langle r \rangle_{1s} = \frac{4}{a_0^3} \int_0^\infty r^3 e^{-2r/a_0} dr = \frac{4}{a_0^3} \cdot \frac{3!}{(2/a_0)^4} = \frac{4}{a_0^3} \cdot \frac{6 a_0^4}{16} = \frac{24 a_0}{16} = \frac{3}{2}a_0

\langle r \rangle_{1s} = 1.5 \times 0.529\text{ Å} = 0.794\text{ Å}

⟨r2⟩1s=∫0∞r2P(r)dr=4a03∫0∞r4e−2r/a0dr=4a03⋅4!(2/a0)5=4a03⋅24a0532=3a02\langle r^2 \rangle_{1s} = \int_0^\infty r^2 P(r) dr = \frac{4}{a_0^3} \int_0^\infty r^4 e^{-2r/a_0} dr = \frac{4}{a_0^3} \cdot \frac{4!}{(2/a_0)^5} = \frac{4}{a_0^3} \cdot \frac{24 a_0^5}{32} = 3a_0^2

\langle r^2 \rangle_{1s} = 3 \times (0.529)^2\text{ Å}^2 = 3 \times 0.280\text{ Å}^2 = 0.840\text{ Å}^2 \Delta r = \sqrt{\langle r^2 \rangle - \langle r \rangle^2} = \sqrt{3a_0^2 - \frac{9}{4}a_0^2} = \sqrt{\frac{3}{4}}a_0 = \frac{\sqrt{3}}{2}a_0 \approx 0.458\text{ Å}

(d) 2p态角动量：

2p态：n=2,l=1n=2, l=1

总角动量大小：

∣L⃗∣=l(l+1)ℏ=1(1+1)ℏ=2ℏ≈1.414ℏ|\vec{L}| = \sqrt{l(l+1)}\hbar = \sqrt{1(1+1)}\hbar = \sqrt{2}\hbar \approx 1.414\hbar

z分量可能值：

Lz=mℏ,m=−1,0,+1L_z = m\hbar, \quad m = -1, 0, +1

即 Lz=−ℏ,0,+ℏL_z = -\hbar, 0, +\hbar

物理讨论：

  • 1s态最概然半径 0.529 Å 恰好等于玻尔半径——这是玻尔模型的伟大成功，也是其局限：玻尔模型预测的是"轨道半径"，而量子力学给出的是"最概然距离"
  • 平均半径 0.794 Å 大于最概然半径——概率密度分布右偏（尾巴拖向大r）
  • 2p态角动量大小 2ℏ\sqrt{2}\hbar 大于玻尔模型的 ℏ\hbar——量子力学修正了经典轨道图像
  • 注意 ∣Lz∣≤∣L⃗∣|L_z| \leq |\vec{L}|：最大z分量 ℏ\hbar 小于总大小 2ℏ\sqrt{2}\hbar，这意味着角动量矢量永远不能"完全指向"z轴——这是角动量量子化的几何体现

本章总结

精确解展示了量子力学代数方法的威力。

graph TD
    A["第VI章: 精确解的凯旋"] --> B["谐振子"]
    A --> C["角动量"]
    A --> D["自旋"]
    A --> E["氢原子"]
    A --> F["选择定则"]
    A --> G["塞曼效应"]
    A --> H["特殊函数基础"]
    
    B --> B1["升降算符 â, â†"]
    B --> B2["能级 Eₙ = ℏω(n+1/2)"]
    B --> B3["零点能 + 等间距"]
    B --> B4["Dirac代数法: 不求解微分方程"]
    
    C --> C1["对易关系 [L̂ᵢ,L̂ⱼ]=iℏεᵢⱼₖL̂ₖ"]
    C --> C2["本征值: L̂²→ℏ²l(l+1), L̂ᵤ→ℏm"]
    C --> C3["球谐函数 Yₗᵐ"]
    
    D --> D1["s = 1/2, 内禀角动量"]
    D --> D2["泡利矩阵 σₓ,σᵧ,σᵤ"]
    D --> D3["[Ŝᵢ,Ŝⱼ] = iℏεᵢⱼₖŜₖ"]
    
    E --> E1["Eₙ = -13.6 eV/n²"]
    E --> E2["库仑势的偶然简并"]
    E --> E3["径向波函数 Rₙₗ(r)"]
    E --> E4["连带Laguerre多项式"]
    
    F --> F1["Δl = ±1, Δm = 0,±1"]
    F --> F2["角动量守恒"]
    
    G --> G1["磁矩 μ̂ = -(e/2m)L̂"]
    G --> G2["能级分裂 ΔE = μ_BB m"]
    G --> G3["正常与反常塞曼效应"]
    G --> G4["朗德g因子"]
    
    H --> H1["Hermite: 谐振子"]
    H --> H2["Legendre/球谐: 角动量"]
    H --> H3["Laguerre: 氢原子径向"]
    
    style B fill:#e3f2fd
    style C fill:#e8f5e9
    style D fill:#fff3e0
    style E fill:#f3e5f5
    style F fill:#ffebee
    style G fill:#e1f5fe
    style H fill:#fce4ec

核心收获：

  1. 谐振子的代数解法（升降算符）展示了算符方法的优雅——不求解微分方程，直接从代数结构得到全部能级。Dirac的方法从哈密顿量的因式分解出发，完全不依赖具体表象，与标准教材"先解方程再引入算符"的路径截然不同
  2. 角动量的本征值问题揭示了量子化来自对易关系的代数约束，而非人为假设。ll 和 mm 的取值范围是升降算符边界条件的自然产出
  3. 电子自旋 s=12s=\frac{1}{2} 是纯量子现象，泡利矩阵提供了简洁的数学描述。自旋磁矩的g因子 gs≈2g_s \approx 2 在Dirac相对论量子力学中才能被完美解释
  4. 氢原子的能级公式 En=−13.6/n2E_n = -13.6/n^2 eV 是量子力学最辉煌的成就之一。库仑势的"偶然简并"（EE 不依赖 ll）源于 1/r1/r 势的隐藏对称性——拉普拉斯-龙格-楞次矢量守恒
  5. 选择定则来自角动量守恒，决定了原子光谱的"指纹"：Δl=±1,Δm=0,±1\Delta l = \pm 1, \Delta m = 0, \pm 1
  6. 塞曼效应展示了外场如何解除简并，光谱线的分裂揭示了磁矩与角动量的深层联系。反常塞曼效应中的朗德g因子 gJg_J 完美解释了复杂的分裂模式
  7. 三类特殊函数（Hermite、Legendre/球谐、连带Laguerre）是量子精确解的数学根基，它们的正交性和完备性保证了量子态空间的结构完整性

练习与思考

1. 谐振子的矩阵元

利用升降算符的性质，计算 ⟨n∣q^∣m⟩\langle n|\hat{q}|m\rangle 和 ⟨n∣p^∣m⟩\langle n|\hat{p}|m\rangle。

(a) 首先将 q^\hat{q} 和 p^\hat{p} 用 a^\hat{a} 和 a^†\hat{a}^\dagger 表示。

(b) 利用 a^∣n⟩=n∣n−1⟩\hat{a}|n\rangle = \sqrt{n}|n-1\rangle 和 a^†∣n⟩=n+1∣n+1⟩\hat{a}^\dagger|n\rangle = \sqrt{n+1}|n+1\rangle 求矩阵元。

© 验证这些矩阵元满足对易关系 [q^,p^]=iℏ[\hat{q}, \hat{p}] = i\hbar。

2. 角动量加法

考虑两个独立的角动量 J⃗1\vec{J}_1 和 J⃗2\vec{J}_2（例如轨道角动量和自旋）。总角动量 J⃗=J⃗1+J⃗2\vec{J} = \vec{J}_1 + \vec{J}_2。

(a) 证明 [J^i,J^j]=iℏϵijkJ^k[\hat{J}_i, \hat{J}_j] = i\hbar\epsilon_{ijk}\hat{J}_k，即总角动量也满足角动量对易关系。

(b) 若 j1=1j_1 = 1 且 j2=12j_2 = \frac{1}{2}，总角动量量子数 jj 可以取哪些值？验证态的总数为 (2j1+1)(2j2+1)(2j_1+1)(2j_2+1)。

© 用克莱布施-戈登系数将 ∣j1=1,m1⟩⊗∣j2=12,m2⟩|j_1=1, m_1\rangle \otimes |j_2=\frac{1}{2}, m_2\rangle 展开为 ∣j,m⟩|j, m\rangle 的线性组合。

3. 氢原子的维里定理

对于氢原子，动能和势能满足维里定理：⟨T⟩=−12⟨V⟩\langle T\rangle = -\frac{1}{2}\langle V\rangle，因此总能量 E=⟨T⟩+⟨V⟩=−⟨T⟩=12⟨V⟩E = \langle T\rangle + \langle V\rangle = -\langle T\rangle = \frac{1}{2}\langle V\rangle。

(a) 从量纲分析出发，解释为什么库仑势情况下维里定理给出 2⟨T⟩=−⟨V⟩2\langle T\rangle = -\langle V\rangle（而非谐振子情况的 ⟨T⟩=⟨V⟩\langle T\rangle = \langle V\rangle）。

(b) 利用氢原子能级公式 En=−13.6n2E_n = -\frac{13.6}{n^2} eV，计算 1s1s 态电子的平均动能 ⟨T⟩1s\langle T\rangle_{1s} 和平均势能 ⟨V⟩1s\langle V\rangle_{1s}。

© 计算 1s1s 态电子的平均距离 ⟨r⟩1s\langle r\rangle_{1s} 和方差 Δr=⟨r2⟩−⟨r⟩2\Delta r = \sqrt{\langle r^2\rangle - \langle r\rangle^2}。

4. 选择定则与球谐函数

利用球谐函数的具体形式，验证：

(a) ⟨Y10∣cosθ∣Y00⟩≠0\langle Y_1^0 | \cos\theta | Y_0^0 \rangle \neq 0（允许跃迁）

(b) ⟨Y20∣cosθ∣Y00⟩=0\langle Y_2^0 | \cos\theta | Y_0^0 \rangle = 0（禁戒跃迁，Δl=2\Delta l = 2）

© 计算 ⟨Y10∣z∣Y00⟩\langle Y_1^0 | z | Y_0^0 \rangle 的具体数值（用 a0a_0 表示）。

"谐振子教会我们升降算符的语言，角动量揭示了几何与代数的统一，氢原子则把这一切编织成光谱的诗篇。精确解之所以珍贵，不仅因为它们是’可解的’，更因为它们是所有近似方法的灯塔——当微扰失效、数值模糊，我们回到这些精确解，重新校准方向。" —— 后记