📖 核心概念 Core Concepts
1. 全同粒子与交换算符 Identical Particles & the Exchange Operator
两个粒子若所有内禀性质(质量、电荷、自旋、磁矩……)完全相同,则称为全同粒子 (identical particles):电子与电子全同,质子与质子全同,但电子与 \(\mu\) 子不全同。经典力学中即使两粒子全同,仍可「涂上颜色」沿轨道追踪「1 号粒子」「2 号粒子」——但在量子力学中没有轨迹,当两个全同粒子的波包重叠后再次分开,原则上无法通过任何测量区分「哪个是 1、哪个是 2」。这种不可分辨性 (indistinguishability) 给波函数加上了严苛的对称约束。
引入交换算符 (exchange operator) \(\hat{P}_{12}\):交换粒子 1 与粒子 2 的所有标签(坐标与自旋)
\[ \hat{P}_{12}\,\Psi(\vec{r}_1,s_1;\,\vec{r}_2,s_2)=\Psi(\vec{r}_2,s_2;\,\vec{r}_1,s_1). \]
再交换一次回到原态:\(\hat{P}_{12}^{2}=\mathbb{1}\),所以 \(\hat{P}_{12}\) 的本征值只有 \(\pm 1\)。又因两粒子全同,所有物理可观测量 \(\hat{A}\)(哈密顿量、动量、能量……)必须与 \(\hat{P}_{12}\) 对易:\([\hat{A},\hat{P}_{12}]=0\)——这意味着 \(\hat{P}_{12}\) 是运动常量,波函数的交换对称性不随时间演化改变。
identical particles (全同粒子), indistinguishability (不可分辨性), exchange operator \(\hat{P}_{12}\) (交换算符), eigenvalue \(\pm 1\) (本征值)
2. 对称化公设 The Symmetrization Postulate
量子力学在自旋-统计定理 (spin-statistics theorem,量子场论中由相对论 + 局域性 + 正能量推出) 中给出最终答案:自然界按自旋把全同粒子分成两类,对应 \(\hat{P}_{12}\) 的两个本征值:
\[ \boxed{\;\text{全同玻色子(整数自旋 }s=0,1,2,\dots\text{):}\hat{P}_{12}\Psi=+\Psi\;\;\text{(对称)}\;} \]
\[ \boxed{\;\text{全同费米子(半整数自旋 }s=\tfrac12,\tfrac32,\dots\text{):}\hat{P}_{12}\Psi=-\Psi\;\;\text{(反对称)}\;} \]
即对称化公设 (symmetrization postulate)。例子:电子、质子、中子、\(\mu\) 子、夸克均为费米子(\(s=1/2\)));\(\tau\)、\(\Omega^-\) 重子(\(s=3/2\))亦为费米子;光子(\(s=1\))、胶子(\(s=1\))、Higgs(\(s=0\))、\(^4\mathrm{He}\) 原子(2 质子+2 中子+2 电子=偶数个费米子合成)为玻色子;而 \({}^3\mathrm{He}\) 原子(2 质子+1 中子+2 电子=奇数个费米子)是费米子——同位素的统计性质截然不同,导致液氦 \({}^4\mathrm{He}\) 在 \(2.17\,\mathrm{K}\) 有超流转变而 \({}^3\mathrm{He}\) 要到 \(2.5\,\mathrm{mK}\) 才转变。
symmetrization postulate (对称化公设), boson (玻色子), fermion (费米子), symmetric (对称), antisymmetric (反对称), spin-statistics theorem (自旋-统计定理)
3. 构造两粒子态:空间 \(\times\) 自旋 Two-Particle States (Spatial \(\otimes\) Spin)
设两粒子可处于单粒子态 \(\psi_a\) 与 \(\psi_b\)。空间部分的对称与反对称组合:
\[ \Psi_{S}^{\text{sp}}(\vec{r}_1,\vec{r}_2)=\frac{1}{\sqrt{2}}\bigl[\psi_a(\vec{r}_1)\psi_b(\vec{r}_2)+\psi_a(\vec{r}_2)\psi_b(\vec{r}_1)\bigr]\;\;\text{(对称)}, \]
\[ \Psi_{A}^{\text{sp}}(\vec{r}_1,\vec{r}_2)=\frac{1}{\sqrt{2}}\bigl[\psi_a(\vec{r}_1)\psi_b(\vec{r}_2)-\psi_a(\vec{r}_2)\psi_b(\vec{r}_1)\bigr]\;\;\text{(反对称)}. \]
(若 \(a=b\),\(\Psi_A^{\text{sp}}\) 恒为零——空间上不可能让两费米子占据同一单粒子态。)两个电子(自旋 \(1/2\))的自旋部分按例 7.5 的结果:三重态 triplet(\(S=1\),对称)与单态 singlet(\(S=0\),反对称):
\[ \chi_{S}^{(1)}\in\text{triplet(对称,3 重)},\qquad \chi_{A}^{(0)}=\frac{1}{\sqrt{2}}\bigl(|\!\uparrow\downarrow\rangle-|\!\downarrow\uparrow\rangle\bigr)\;\text{(singlet,反对称,1 重)}. \]
费米子总波函数(空间 \(\otimes\) 自旋)必须反对称。两种组合都满足:
\[ \underbrace{\Psi_{A}^{\text{sp}}\otimes\chi_{S}^{(1)}}_{\text{反对称}\,\times\,\text{对称}=反对称}\;(\text{三重态正氦 ortho}),\qquad \underbrace{\Psi_{S}^{\text{sp}}\otimes\chi_{A}^{(0)}}_{\text{对称}\,\times\,\text{反对称}=反对称}\;(\text{单态仲氦 para}). \]
所以即便两电子处于同一空间轨道(如氦基态 \(1s^2\),\(\Psi^{\text{sp}}\) 对称),只要它们自旋构成 singlet(反向平行)总波函数仍反对称——这就是「电子配对填轨道」的量子起源。反之,两电子若自旋平行(triplet,对称),空间部分必须反对称,意味着它们倾向占据不同空间区域——这就是 Hund 第一规则的根源。
spatial part \(\Psi^{\text{sp}}\) (空间部分), spin part \(\chi\) (自旋部分), singlet (单态), triplet (三重态), para-helium (仲氦), ortho-helium (正氦), Hund's rule (洪德规则)
4. 斯莱特行列式与泡利不相容 Slater Determinant & Pauli Exclusion
把两费米子构造推广到 \(N\) 个费米子占据正交归一的单粒子态 \(\{\psi_1,\psi_2,\dots,\psi_N\}\):唯一反对称的归一化波函数是斯莱特行列式 (Slater determinant)
\[ \boxed{\;\Psi(\xi_1,\dots,\xi_N)=\frac{1}{\sqrt{N!}}\begin{vmatrix}\psi_1(\xi_1)&\psi_1(\xi_2)&\cdots&\psi_1(\xi_N)\\ \psi_2(\xi_1)&\psi_2(\xi_2)&\cdots&\psi_2(\xi_N)\\ \vdots& & &\vdots\\ \psi_N(\xi_1)&\psi_N(\xi_2)&\cdots&\psi_N(\xi_N)\end{vmatrix}\;} \]
其中 \(\xi_i\equiv(\vec{r}_i,s_i)\) 表示第 \(i\) 个粒子的坐标+自旋复合标签。三条性质自动由行列式给出:
- 反对称:交换两粒子 = 交换两列 = 行列式变号。
- 泡利不相容:若两行相同(即两个费米子占据同一单粒子态 \(\psi_i=\psi_j\)),行列式为零——不存在这种态。
- 归一化:行列式展开有 \(N!\) 项正交乘积,故归一化为 \(1/\sqrt{N!}\)。
泡利不相容原理 (Pauli exclusion principle, 1925) 的现代形式:每个单粒子量子态 \((n,l,m_l,m_s)\) 至多容纳一个电子。原子壳层容量为 \(2(2l+1)\):\(1s\) 容 2、\(2s\) 容 2、\(2p\) 容 6、\(3d\) 容 10……填满顺序对应周期表中元素的电子排布。化学的全部根基即是这一条量子规则。
Slater determinant (斯莱特行列式), Pauli exclusion principle (泡利不相容原理), shell capacity \(2(2l+1)\) (壳层容量), fermion not allowed in same state (不能同态)
5. 氦原子:基态、激发态、交换积分 Two-Electron Atoms — Helium
氦原子(\(Z=2\),核外两电子)是最简单的多电子原子,是检验全同粒子物理的试金石。哈密顿量(玻恩–奥本海默近似下,核固定):
\[ \hat{H}=\underbrace{\Bigl(-\frac{\hbar^{2}}{2m_e}\nabla_1^{2}-\frac{2e^{2}}{4\pi\epsilon_0 r_1}\Bigr)}_{\hat{h}_1\;(\text{氢原子状})}+\underbrace{\Bigl(-\frac{\hbar^{2}}{2m_e}\nabla_2^{2}-\frac{2e^{2}}{4\pi\epsilon_0 r_2}\Bigr)}_{\hat{h}_2}+\underbrace{\frac{e^{2}}{4\pi\epsilon_0 r_{12}}}_{\hat{V}_{ee}\;(\text{电子-电子排斥})}, \]
其中 \(r_{12}=|\vec{r}_1-\vec{r}_2|\)。前两项是两个独立的 \(\mathrm{He}^{+}\)(类氢,\(Z=2\)),可严格求解;第三项电子-电子排斥使问题不可分离变量,必须近似。
基态 \(1s^{2}\):忽略 \(\hat{V}_{ee}\) 时两电子都处于 \(\psi_{1s}\)(类氢 \(Z=2\),能量 \(-54.4\,\mathrm{eV}\) 各),总 \(-108.8\,\mathrm{eV}\)。考虑排斥后实际基态能量 \(-79.0\,\mathrm{eV}\)(实验:第一电离 \(24.6\,\mathrm{eV}\) + 第二电离 \(54.4\,\mathrm{eV}\))。空间部分对称(两电子同轨道),故自旋必为 singlet(仲氦,\(S=0\))。
激发态 \(1s\,2s\) 与 \(1s\,2p\):一个电子留在 \(1s\),另一个跃迁到 \(2s\) 或 \(2p\)。空间部分可对称(singlet,仲氦 para)或反对称(triplet,正氦 ortho)。两种态能量不同!原因:电子-电子排斥的期望依赖空间对称性。
定义直接积分 (direct integral) \(J\) 与交换积分 (exchange integral) \(K\):
\[ J=\int\!\!\int|\psi_{1s}(\vec{r}_1)|^{2}|\psi_{2s}(\vec{r}_2)|^{2}\,\frac{e^{2}}{4\pi\epsilon_0 r_{12}}\,\mathrm{d}^{3}r_1\,\mathrm{d}^{3}r_2\;\;(\text{经典型库仑自能}), \]
\[ K=\int\!\!\int\psi_{1s}^{*}(\vec{r}_1)\psi_{2s}^{*}(\vec{r}_2)\,\frac{e^{2}}{4\pi\epsilon_0 r_{12}}\,\psi_{1s}(\vec{r}_2)\psi_{2s}(\vec{r}_1)\,\mathrm{d}^{3}r_1\,\mathrm{d}^{3}r_2\;\;(\text{纯量子,无经典对应}). \]
则两态能量为
\[ E_{\text{singlet}}=E_{1s}+E_{2s}+J+K\;\;\text{(para, 空间对称)},\qquad E_{\text{triplet}}=E_{1s}+E_{2s}+J-K\;\;\text{(ortho, 空间反对称)}. \]
实验测得 \(1s\,2s\) 组态:\({}^{1}S\)(para)在基态之上 \(20.62\,\mathrm{eV}\),\({}^{3}S\)(ortho)在 \(19.82\,\mathrm{eV}\)。三重态更低,劈裂 \(\Delta E=E_s-E_t=2K=0.80\,\mathrm{eV}\),故 \(K\simeq 0.40\,\mathrm{eV}\)。物理上:空间反对称时电子「互相躲开」(\(\Psi(\vec{r},\vec{r})=0\),同一位置概率为零),平均距离更大,排斥能更小——这就是交换作用 (exchange interaction),是铁磁性、Hund 规则、化学键方向的根源,纯粹量子效应。
two-electron atom (双电子原子), direct integral \(J\) (直接积分), exchange integral \(K\) (交换积分), exchange interaction (交换作用), ortho/para-helium (正氦/仲氦)
6. 变分原理与 Hartree–Fock 概览 Variational Principle & Hartree–Fock
多电子原子无法严格求解,需要近似。变分原理 (variational principle) 给出最简单可靠的基态能量上界:
\[ \boxed{\;E[\psi_{\text{trial}}]\equiv\frac{\langle\psi_{\text{trial}}|\hat{H}|\psi_{\text{trial}}\rangle}{\langle\psi_{\text{trial}}|\psi_{\text{trial}}\rangle}\;\ge\;E_0\;}\qquad(\text{对任意归一化试探波函数}). \]
证明见推导 1.B:把 \(\psi_{\text{trial}}\) 按严格本征态 \(|n\rangle\) 展开 \(\psi_{\text{trial}}=\sum_n c_n|n\rangle\),则 \(E[\psi_{\text{trial}}]=\sum_n|c_n|^{2}E_n\ge\sum_n|c_n|^{2}E_0=E_0\)。变分法:选含可调参数 \(\{\alpha\}\) 的试探函数,最小化 \(E(\alpha)\) → \(\partial E/\partial\alpha=0\) → 最优 \(\alpha^{*}\) 给出最优上界。例 1.3 用此估氦基态能量。
Hartree–Fock 自洽场 (Hartree–Fock self-consistent field, HF) 是变分原理的系统化:用单个斯莱特行列式作为 \(N\) 电子试探态,对每个空间轨道 \(\phi_i(\vec{r})\) 变分得一组耦合本征方程(Hartree–Fock 方程),其中每个电子在其他 \(N-1\) 个电子的平均场中运动,并显式包含交换项(来自行列式的反对称性)。HF 是现代量子化学(密度泛函 DFT、后 Hartree–Fock 方法)的出发点。
variational principle (变分原理), trial wavefunction (试探波函数), upper bound (上界), Hartree–Fock (HF 自洽场), self-consistent field (自洽场)
🔭 推导 Derivation
推导 1.A — 从交换算符到泡利不相容 From the Exchange Operator to Pauli Exclusion
步骤 1(本征值约束):\(\hat{P}_{12}^{2}=\mathbb{1}\)(两次交换回到原态)。设 \(|\Psi\rangle\) 是 \(\hat{P}_{12}\) 的本征态,本征值 \(\lambda\):\(\hat{P}_{12}^{2}|\Psi\rangle=\lambda^{2}|\Psi\rangle=|\Psi\rangle\Rightarrow\lambda=\pm 1\)。
步骤 2(构造本征态):给定任意直积态 \(\psi_a(1)\psi_b(2)\)(用简写 \(1\equiv\xi_1\),\(2\equiv\xi_2\)),其交换为 \(\psi_a(2)\psi_b(1)\)。构造对称 / 反对称组合:
\[ |\Psi_{\pm}\rangle=\frac{1}{\sqrt{2}}\bigl(\psi_a(1)\psi_b(2)\pm\psi_a(2)\psi_b(1)\bigr). \]
步骤 3(验证本征值):\(\hat{P}_{12}|\Psi_{\pm}\rangle=\tfrac{1}{\sqrt{2}}\bigl(\psi_a(2)\psi_b(1)\pm\psi_a(1)\psi_b(2)\bigr)=\pm\tfrac{1}{\sqrt{2}}\bigl(\psi_a(1)\psi_b(2)\pm\psi_a(2)\psi_b(1)\bigr)=\pm|\Psi_{\pm}\rangle\) ✓(加号时直接等;减号时提出负号)。故 \(\Psi_+\) 对应 \(\lambda=+1\)(玻色子),\(\Psi_-\) 对应 \(\lambda=-1\)(费米子)。
步骤 4(泡利不相容):令 \(a=b\)(两费米子欲占同一单粒子态)。则 \(\Psi_{-}=\tfrac{1}{\sqrt{2}}(\psi_a(1)\psi_a(2)-\psi_a(2)\psi_a(1))=\mathbf{0}\)。波函数恒为零——不存在这种态。即「两个全同费米子不能处于完全相同的单粒子量子态」。但 \(\Psi_{+}=\sqrt{2}\,\psi_a(1)\psi_a(2)\ne 0\),玻色子不受此约束(Bose–Einstein 凝聚即大量玻色子凝聚到同一基态)。
步骤 5(坐标+自旋推广):当 \(\psi\) 含自旋指标 \(\xi=(\vec{r},s)\) 时,「同一单粒子态」指 \((n,l,m_l,m_s)\) 完全相同。两电子可同处 \(1s\) 轨道只要 \(m_s\) 不同(一 \(\uparrow\)、一 \(\downarrow\)),因为此时 \(\psi_a\ne\psi_b\),行列式不为零。这就是「电子成对填轨道」的量子起源。
步骤 6(与 \(N=2\) 斯莱特行列式一致):\(\Psi_{-}(1,2)=\tfrac{1}{\sqrt{2}}\begin{vmatrix}\psi_a(1)&\psi_a(2)\\\psi_b(1)&\psi_b(2)\end{vmatrix}\) 展开即 \(\tfrac{1}{\sqrt{2}}(\psi_a(1)\psi_b(2)-\psi_a(2)\psi_b(1))\) ✓。换列变号 = 反对称;两行相同(\(a=b\))= 行列式零 = 泡利。斯莱特行列式自动给出步骤 4–5 的所有结论。
推导 1.B — 变分原理证明与氦基态能量 Proof of Variational Principle & Helium Ground State
步骤 1(变分原理证明):设 \(\hat{H}|n\rangle=E_n|n\rangle\),\(E_0\le E_1\le\dots\) 为严格能谱。任意归一化 \(|\psi\rangle=\sum_n c_n|n\rangle\),\(\sum|c_n|^{2}=1\)。则
\[ \langle\psi|\hat{H}|\psi\rangle=\sum_n|c_n|^{2}E_n\;\ge\;\sum_n|c_n|^{2}E_0=E_0\;\;\surd. \]
等号当且仅当 \(c_0=1\)、其余为零(即 \(|\psi\rangle=|0\rangle\) 严格基态)。故任何试探态给出基态能量的上界。
步骤 2(氦的试探函数):忽略 \(\hat{V}_{ee}\) 时基态为两电子同处类氢 \(1s\)(核电荷 \(Z=2\))。考虑电子互相屏蔽,用有效核电荷 \(Z_{\text{eff}}\) 替代 \(Z\):
\[ \psi_{\text{trial}}(\vec{r}_1,\vec{r}_2)=\frac{Z_{\text{eff}}^{3}}{\pi a_0^{3}}\exp\!\Bigl[-\frac{Z_{\text{eff}}(r_1+r_2)}{a_0}\Bigr]\;\;\text{(已归一化)}, \]
把 \(Z_{\text{eff}}\) 当作变分参数。当 \(Z_{\text{eff}}=Z=2\) 时回归无屏蔽类氢近似。
步骤 3(动能期望):由维里定理,类氢 \(1s\)(参数 \(Z_{\text{eff}}\))满足 \(\langle T\rangle_{\text{单}}=-\tfrac12\langle V_{Z_{\text{eff}}}\rangle_{\text{单}}\),且总能量 \(-\tfrac12 Z_{\text{eff}}^{2}\alpha^{2}m_ec^{2}=-13.6\,Z_{\text{eff}}^{2}\,\mathrm{eV}\),故 \(\langle T\rangle_{\text{单}}=+13.6\,Z_{\text{eff}}^{2}\,\mathrm{eV}\)。两电子:
\[ \langle T\rangle=2\times 13.6\,Z_{\text{eff}}^{2}=27.2\,Z_{\text{eff}}^{2}\,\mathrm{eV}. \]
步骤 4(核-电子吸引期望):真实核电荷 \(Z=2\),但试探轨道用 \(Z_{\text{eff}}\),故 \(\langle 1/r\rangle_{Z_{\text{eff}}}=Z_{\text{eff}}/a_0\)。
\[ \langle V_{ne}\rangle_{\text{单}}=\Big\langle-\frac{Ze^{2}}{4\pi\epsilon_0 r}\Big\rangle=-Z\cdot\frac{e^{2}}{4\pi\epsilon_0}\cdot\frac{Z_{\text{eff}}}{a_0}=-27.2\,Z\,Z_{\text{eff}}\,\mathrm{eV}. \]
两电子(\(Z=2\)):\(\langle V_{ne}\rangle=-4\times 27.2\,Z_{\text{eff}}=-108.8\,Z_{\text{eff}}\,\mathrm{eV}\)。
步骤 5(电子-电子排斥期望):对两 \(1s\) 轨道(参数 \(Z_{\text{eff}}\)),标准结果
\[ \Big\langle\frac{e^{2}}{4\pi\epsilon_0 r_{12}}\Big\rangle=\frac{5}{8}\,Z_{\text{eff}}\cdot\frac{e^{2}}{4\pi\epsilon_0 a_0}=\frac{5}{8}Z_{\text{eff}}\times 27.2=17.0\,Z_{\text{eff}}\,\mathrm{eV}. \]
步骤 6(总能量与极小化):
\[ E(Z_{\text{eff}})=27.2\,Z_{\text{eff}}^{2}-108.8\,Z_{\text{eff}}+17.0\,Z_{\text{eff}}=27.2\,Z_{\text{eff}}^{2}-91.8\,Z_{\text{eff}}\;\;(\mathrm{eV}). \]
令 \(\mathrm{d}E/\mathrm{d}Z_{\text{eff}}=54.4\,Z_{\text{eff}}-91.8=0\),得
\[ \boxed{\;Z_{\text{eff}}^{*}=\frac{91.8}{54.4}=\frac{27}{16}=Z-\frac{5}{16}=1.6875\;}\;\;(\text{屏蔽掉了 }5/16\text{ 个核电荷}). \]
代入:\(E_{\min}=27.2\times(27/16)^{2}-91.8\times(27/16)=77.46-154.91=\mathbf{-77.5\,\mathrm{eV}}\)。实验值 \(-79.0\,\mathrm{eV}\),误差仅 \(1.5\,\mathrm{eV}\)(约 \(1.9\%\));而未变分(\(Z_{\text{eff}}=2\))给出 \(E=-74.8\,\mathrm{eV}\),误差 \(4.2\,\mathrm{eV}\)。简单的屏蔽变分把误差削减了一大半。