自旋与全同粒子：电子的隐身特性与不可区分的魔法
电子除了轨道角动量还藏着一个「内禀角动量」自旋；多个相同的粒子一旦混在一起就完全无法区分。这两个概念塑造了原子结构和整个宇宙的统计规律。
作业 7.1–7.6；重点：自旋角动量大小与z分量本征值、泡利矩阵及其性质（厄米、平方为单位阵、反对易关系）、自旋函数α和β、自旋算符对易关系、角动量耦合与总角动量、全同粒子与泡利不相容原理、费米子反对称波函数、玻色子对称波函数、两电子自旋单态与三重态、仲氦与正氦、多粒子对称化波函数构造。课堂重点：N=3玻色子三单粒子态的对称化波函数（补充例题5）。
电子自旋是电子内禀的角动量（无经典对应），通过斯特恩-革拉赫实验被发现——银原子束在磁场中一分为二。自旋-轨道耦合导致原子光谱出现精细结构。全同粒子是质量电荷自旋完全相同的粒子，它们不可区分；交换两个粒子后，体系波函数要么对称（玻色子），要么反对称（费米子）。这导出泡利不相容原理（费米子不能占据同一量子态），是化学元素周期表的根源。
7.1 电子自旋与斯特恩-革拉赫实验
自旋是电子的内禀角动量，不是电子绕自身轴旋转（这样的速度会超光速，物理上不可能），而是电子作为基本粒子固有携带的角动量——一个量子力学专有的性质，无经典对应。
每个电子都具有自旋角动量，自旋量子数 \(s=\dfrac{1}{2}\)。自旋角动量的大小为：
自旋在任意方向（如z方向）的投影 \(S_z\) 只能取两个值：
即 \(S_z = +\dfrac{\hbar}{2}\)（自旋"向上"）或 \(S_z = -\dfrac{\hbar}{2}\)（自旋"向下"）。
📌斯特恩-革拉赫实验（Stern-Gerlach）▸
1922 年，一束银原子（基态电子外层只有一个 5s 电子，\(L=0\)）通过非均匀磁场。如果角动量是连续分布，应该看到一条模糊的细长束；实际上看到严格的两束分束——直接证明了自旋的空间量子化只有两个可能取向。这是自旋存在的直接证据。
推论：自旋角动量在z方向只能是 \(\pm\hbar/2\)，且外加磁场无法改变这一点——自旋取向是量子化的，无法连续变化。
7.2 自旋算符与泡利矩阵
自旋是一种角动量，必须满足角动量对易关系。记自旋算符为 \(\hat{S}_x, \hat{S}_y, \hat{S}_z\)，它们满足：
定义泡利算符 \(\hat{\sigma}_i\)，使得 \(\hat{S}_i = \dfrac{\hbar}{2}\hat{\sigma}_i\)。泡利矩阵为：
- 厄米性： \(\hat{\sigma}_i^\dagger = \hat{\sigma}_i\)（都是厄米矩阵，本征值必为实数）
- 平方为单位阵： \(\hat{\sigma}_i^2 = I\)（本征值为 \(\pm 1\)）
- 反对易关系： \(\{\hat{\sigma}_i, \hat{\sigma}_j\} = \hat{\sigma}_i\hat{\sigma}_j + \hat{\sigma}_j\hat{\sigma}_i = 2\delta_{ij}I\)
- 对易关系： \([\hat{\sigma}_i, \hat{\sigma}_j] = 2i\,\varepsilon_{ijk}\,\hat{\sigma}_k\)（例如 \([\hat{\sigma}_x, \hat{\sigma}_y] = 2i\hat{\sigma}_z\)）
自旋函数（自旋本征态）用二分量列向量表示：
它们是 \(\hat{S}_z\)（或 \(\hat{\sigma}_z\)）的本征态：\(\hat{S}_z\chi_\pm = \pm\frac{\hbar}{2}\chi_\pm\)。
7.3 简单塞曼效应
将类氢原子置于外部均匀强磁场中，原子的能级会劈裂。最简单的情况是只考虑自旋z分量与磁场的相互作用（轨道角动量已被外磁场充分作用）。设外磁场沿z方向，能量修正为：
其中 \(\mu_B = \dfrac{e\hbar}{2m_e}\) 是玻尔磁子（量子力学中磁矩的天然单位），\(g_e\approx 2\)（电子g因子），\(m_s = \pm 1/2\)。
对于一条原始的谱线（如 \(n_2 \to n_1\) 的跃迁），在磁场中劈裂成3条（一分三）：各电子自旋向上、向下、轨道角动量m值不同组合产生三条线。
7.4 角动量耦合与精细结构
原子中电子同时具有轨道角动量 \(\hat{L}\) 和自旋角动量 \(\hat{S}\)，两者会相互耦合。定义总角动量：
\(\hat{J}^2\) 和 \(\hat{J}_z\) 是守恒量。总角动量量子数 \(j\) 由以下规则确定：
对于 \(s = 1/2\)：
- 当 \(l \ne 0\) 时：\(j = l + \frac{1}{2}\) 或 \(j = l - \frac{1}{2}\)（两个j值）
- 当 \(l = 0\) 时：\(j = \frac{1}{2}\)（唯一的j值）
自旋-轨道耦合能将同一能级 \(n,l\) 的状态分裂成多个j值的子能级，这被称为光谱精细结构（如碱金属原子光谱的双线结构 D线）。
把轨道角动量和自旋角动量比作两个陀螺，它们不是独立的，而是互相"搂抱"转动（耦合），形成一个共同的旋转轴（总角动量）。这两个陀螺搂抱得越紧（耦合常数越大），分裂的能量差越大。
7.5 全同粒子与全同性原理
全同粒子（identical particles）：质量、电荷、自旋、磁矩等所有内禀属性完全相同的微观粒子。所有电子彼此全同；所有光子彼此全同；等等。
全同性原理：全同粒子完全不可区分。交换任意两个全同粒子，物理状态不变，这意味着体系的哈密顿量对粒子交换对称：\(\hat{H}(\mathbf{r}_1,\mathbf{r}_2,\ldots) = \hat{H}(\mathbf{r}_2,\mathbf{r}_1,\ldots)\)。
交换两粒子的数学操作用交换算符 \(\hat{P}_{12}\) 表示：\(\hat{P}_{12}\Psi(\mathbf{r}_1,\mathbf{r}_2) = \Psi(\mathbf{r}_2,\mathbf{r}_1)\)。由于 \([\hat{H}, \hat{P}_{12}] = 0\)，交换对称性是守恒量——一旦体系波函数具有某种交换对称性，在演化中始终保持。
全同粒子的波函数必须对称或反对称（不能是一般的线性组合），这是量子力学的基本要求。不是说我们可以选择对称或反对称——交换对称性由粒子的自旋决定，是自然界的硬性规律。
7.6 玻色子与费米子
| 类型 | 自旋量子数 | 波函数交换 | 统计 | 例子 |
|---|---|---|---|---|
| 玻色子(Boson) | 整数(0,1,2,...) | 对称 | 玻色-爱因斯坦 | 光子(s=1)、π介子(s=0)、基态氦原子(s=0) |
| 费米子(Fermion) | 半整数(1/2,3/2,...) | 反对称 | 费米-狄拉克 | 电子/质子/中子(s=1/2) |
这个分类是自旋统计定理（spin-statistics theorem）的结果——整数自旋→对称；半整数自旋→反对称。
7.7 泡利不相容原理
不能有两个（或多个）全同费米子处于同一个单粒子量子态。
换言之，对于费米子体系，总波函数对粒子交换必须反对称。这导致——原子中不能有两个电子具有相同的四个量子数 \((n,l,m_l,m_s)\)；每个单粒子轨道最多容纳2个电子（自旋向上和向下各一个）。
这个原理完全解释了原子的壳层结构和化学元素周期表：电子逐层填充各轨道，满足不相容原理，形成稳定的电子构型。
费米子像性格古怪的演员，每个演员（量子态）的房间里最多住两个人，而且这两个人必须面对面坐着（反对称波函数）。玻色子像友好的羊群，所有的羊都愿意挤在同一个房间里（对称波函数），甚至同一个位置也不在乎（这导致凝聚现象）。
7.8 两电子自旋函数与氦原子
两个电子体系的总自旋量子数 \(S\) 可以是 0 或 1：
关键性质：
- 单态 \(S=0\)：自旋反对称（两自旋相反）
- 三重态 \(S=1\)：自旋对称（两自旋平行或一对一反）
由于费米子总波函数必须反对称，而 \(\Psi_{\text{total}} = \Psi_{\text{spatial}} \times \Psi_{\text{spin}}\)，所以：
\[ \text{单态}(S=0)：\Psi_{\text{spatial}} \text{ 对称} \times \Psi_{\text{spin}} \text{ 反对称} = \text{反对称} \,\checkmark \\ \text{三重态}(S=1)：\Psi_{\text{spatial}} \text{ 反对称} \times \Psi_{\text{spin}} \text{ 对称} = \text{反对称} \,\checkmark \]
氦原子应用：
- 仲氦（orthohelium）： \(S=0\) 自旋单态，自旋反对称 → 空间必须对称 → 两电子可占据相同空间，能量低
- 正氦（parahelium）： \(S=1\) 自旋三重态，自旋对称 → 空间必须反对称 → 两电子互相排斥，能量高
这解释了氦的两个相对稳定的激发态能量差异。在实验中观测到正氦的三条谱线（S=1,m_s=+1,0,-1）和仲氦的复杂光谱。
7.9 N粒子对称化波函数的构造
对于N个全同粒子体系，若有k个可用的单粒子态 \(\psi_1, \psi_2, \ldots, \psi_k\)，体系需要用对称或反对称的组合表示。
玻色子（对称）的情况：每个单粒子态可被任意多个粒子占据。设第i个态被 \(n_i\) 个粒子占据（\(\sum n_i = N\)），对称化波函数为：
实际上，对称化就是对所有粒子标号的排列求和，然后除以重复因子以保证归一化。
费米子（反对称）的情况：每个单粒子态最多被1个粒子占据。N个粒子必须占据N个不同的单粒子态，波函数用Slater行列式表示：
行列式的性质自动保证：交换任意两行（粒子），行列式变号（反对称）；若两行（两个单粒子态）相同，行列式为零（泡利原理）。
先盖住解答自己动手，卡住了再点开看分步过程。👇
写出泡利矩阵
\[ \hat{\sigma}_x = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \quad \hat{\sigma}_y = \begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}, \quad \hat{\sigma}_z = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} \]
计算 \(\hat{\sigma}_x\hat{\sigma}_y\)
\[ \hat{\sigma}_x\hat{\sigma}_y = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}\begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix} = \begin{pmatrix} i & 0 \\ 0 & -i \end{pmatrix} \]
计算 \(\hat{\sigma}_y\hat{\sigma}_x\)
\[ \hat{\sigma}_y\hat{\sigma}_x = \begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} = \begin{pmatrix} -i & 0 \\ 0 & i \end{pmatrix} \]
求对易子
\[ [\hat{\sigma}_x, \hat{\sigma}_y] = \hat{\sigma}_x\hat{\sigma}_y - \hat{\sigma}_y\hat{\sigma}_x = \begin{pmatrix} i & 0 \\ 0 & -i \end{pmatrix} - \begin{pmatrix} -i & 0 \\ 0 & i \end{pmatrix} = \begin{pmatrix} 2i & 0 \\ 0 & -2i \end{pmatrix} = 2i\hat{\sigma}_z \,\checkmark \]
列本征方程与久期方程
设本征值为 \(\lambda\)，本征函数为 \(\begin{pmatrix} c_1 \\ c_2 \end{pmatrix}\)，则：
\[ \hat{\sigma}_x\begin{pmatrix} c_1 \\ c_2 \end{pmatrix} = \lambda\begin{pmatrix} c_1 \\ c_2 \end{pmatrix} \]
久期方程：\(\det(\hat{\sigma}_x - \lambda I) = 0\)，即 \(\begin{vmatrix} -\lambda & 1 \\ 1 & -\lambda \end{vmatrix} = \lambda^2 - 1 = 0\)。
求本征值
\[ \lambda^2 = 1 \quad \Rightarrow \quad \lambda = +1 \text{ 或 } -1 \]
求本征函数（\(\lambda = +1\)）
\[ (\hat{\sigma}_x - I)\begin{pmatrix} c_1 \\ c_2 \end{pmatrix} = 0 \quad \Rightarrow \quad \begin{pmatrix} -1 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} c_1 \\ c_2 \end{pmatrix} = 0 \]
得 \(c_1 = c_2\)，故本征函数为 \(\begin{pmatrix} 1 \\ 1 \end{pmatrix}\)。归一化：\(\begin{pmatrix} 1/\sqrt{2} \\ 1/\sqrt{2} \end{pmatrix}\)。
求本征函数（\(\lambda = -1\)）
\[ (\hat{\sigma}_x + I)\begin{pmatrix} c_1 \\ c_2 \end{pmatrix} = 0 \quad \Rightarrow \quad \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}\begin{pmatrix} c_1 \\ c_2 \end{pmatrix} = 0 \]
得 \(c_1 = -c_2\)，故本征函数为 \(\begin{pmatrix} 1 \\ -1 \end{pmatrix}\)。归一化：\(\begin{pmatrix} 1/\sqrt{2} \\ -1/\sqrt{2} \end{pmatrix}\)。
计算 \(\langle S_x \rangle\)
由 \(\hat{S}_x = \frac{\hbar}{2}\hat{\sigma}_x\)，
\[ \langle S_x \rangle = \frac{\hbar}{2}\langle\chi|\hat{\sigma}_x|\chi\rangle = \frac{\hbar}{2} \cdot \frac{1}{2}(\alpha + \beta)^\dagger \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}(\alpha + \beta) \]
\[ = \frac{\hbar}{2} \cdot \frac{1}{2}(\langle\alpha| + \langle\beta|)(\beta + \alpha) = \frac{\hbar}{2} \cdot \frac{1}{2} \cdot 2 = \frac{\hbar}{4} \]
实际上，该态是 \(\hat{\sigma}_x\) 的本征态，本征值 +1，所以 \(\langle S_x \rangle = \frac{\hbar}{2} \cdot 1 = \frac{\hbar}{2}\)。
计算 \(\langle S_y \rangle\)
\[ \langle S_y \rangle = \frac{\hbar}{2}\langle\chi|\hat{\sigma}_y|\chi\rangle = \frac{\hbar}{2} \cdot \frac{1}{2}(\langle\alpha| + \langle\beta|)\begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}(\alpha + \beta) \]
\[ = \frac{\hbar}{2} \cdot \frac{1}{2}(\langle\alpha|(-i\beta) + \langle\beta|(i\alpha)) = \frac{\hbar}{2} \cdot \frac{1}{2}(-i + i) = 0 \]
计算 \(\langle S_z \rangle\)
\[ \langle S_z \rangle = \frac{\hbar}{2}\langle\chi|\hat{\sigma}_z|\chi\rangle = \frac{\hbar}{2} \cdot \frac{1}{2}(\langle\alpha| + \langle\beta|)\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}(\alpha + \beta) \]
\[ = \frac{\hbar}{2} \cdot \frac{1}{2}(\langle\alpha|\alpha - \langle\beta|\beta) = \frac{\hbar}{2} \cdot \frac{1}{2}(1 - 1) = 0 \]
计算平方的平均值
由 \(\hat{\sigma}_x^2 = \hat{\sigma}_y^2 = I\)，得 \(\langle\hat{S}_x^2\rangle = \langle\hat{S}_y^2\rangle = \frac{\hbar^2}{4}\)（所有自旋态都满足）。
\[ \Delta S_x = \sqrt{\langle\hat{S}_x^2\rangle - \langle S_x\rangle^2} = \sqrt{\frac{\hbar^2}{4} - \frac{\hbar^2}{4}} = 0 \]
\[ \Delta S_y = \sqrt{\langle\hat{S}_y^2\rangle - \langle S_y\rangle^2} = \sqrt{\frac{\hbar^2}{4} - 0} = \frac{\hbar}{2} \]
识别量子数
第一项：\(n=2, l=1, m_l=1, m_s=+1/2\)。第二项：\(n=2, l=1, m_l=0, m_s=-1/2\)。
系数各为 \(1/\sqrt{2}\)，都归一化了。
求 \(\langle L_z \rangle\)
\[ \langle L_z \rangle = \frac{1}{2}\int d^3\mathbf{r}\,[\psi_{2,1,1}^*\alpha^\dagger + \psi_{2,1,0}^*\beta^\dagger] \hat{L}_z [\psi_{2,1,1}\alpha + \psi_{2,1,0}\beta] \]
由于 \(\hat{L}_z\) 只作用在空间部分，本征值分别为 \(m_l\hbar = 1\cdot\hbar\) 和 \(0\cdot\hbar\)，
\[ \langle L_z \rangle = \frac{1}{2}(1\cdot\hbar + 0\cdot\hbar) = \frac{\hbar}{2} \]
求 \(\langle S_z \rangle\)
第一项：\(\hat{S}_z\alpha = \frac{\hbar}{2}\alpha\)，概率 \(1/2\)；第二项：\(\hat{S}_z\beta = -\frac{\hbar}{2}\beta\)，概率 \(1/2\)。
\[ \langle S_z \rangle = \frac{1}{2}\cdot\frac{\hbar}{2} + \frac{1}{2}\cdot(-\frac{\hbar}{2}) = 0 \]
求总磁矩的z分量
原子的总磁矩为 \(\vec{\mu} = -\frac{e}{2m_e}(\vec{L} + 2\vec{S})\)（因子2来自电子g因子 \(\approx 2\)）。
\[ \langle\mu_z\rangle = -\frac{e}{2m_e}(\langle L_z\rangle + 2\langle S_z\rangle) = -\frac{e}{2m_e}\cdot\frac{\hbar}{2} = -\frac{1}{2}\cdot\frac{e\hbar}{2m_e} = -\frac{1}{2}\mu_B \]
（负号表示磁矩沿z负方向。用玻尔磁子 \(\mu_B = e\hbar/(2m_e)\)。）
列举粒子分布情况
三个粒子在两个态的分布用 \((n_1, n_2)\) 表示，其中 \(n_1 + n_2 = 3\)：
\((3,0), (2,1), (1,2), (0,3)\) —— 共 4 种。（答：体系有 4 个可能状态）
情况1：分布 (3,0) —— 三个粒子都在态1
\[ \Psi_1 = \psi_1(1)\psi_1(2)\psi_1(3) \]
（已是全对称）
情况2：分布 (2,1) —— 两个粒子在态1，一个在态2
对所有粒子位置的排列求和，除以重复因子 \(\sqrt{2!\cdot 1!}\)：
\[ \Psi_2 = \frac{1}{\sqrt{3}}[\psi_1(1)\psi_1(2)\psi_2(3) + \psi_1(1)\psi_2(2)\psi_1(3) + \psi_2(1)\psi_1(2)\psi_1(3)] \]
或简写为（对称化求和）：
\[ \Psi_2 = \frac{1}{\sqrt{3}}\sum_{\text{perm}} \psi_1(i_1)\psi_1(i_2)\psi_2(i_3) \]
情况3：分布 (1,2) —— 一个粒子在态1，两个在态2
\[ \Psi_3 = \frac{1}{\sqrt{3}}[\psi_1(1)\psi_2(2)\psi_2(3) + \psi_2(1)\psi_1(2)\psi_2(3) + \psi_2(1)\psi_2(2)\psi_1(3)] \]
情况4：分布 (0,3) —— 三个粒子都在态2
\[ \Psi_4 = \psi_2(1)\psi_2(2)\psi_2(3) \]
计算 \(\langle\chi_0|\chi_0\rangle\)（单态的归一化）
\[ \langle\chi_0|\chi_0\rangle = \frac{1}{2}(\langle\alpha_1\beta_2| - \langle\beta_1\alpha_2|)(\alpha_1\beta_2 - \beta_1\alpha_2) \]
\[ = \frac{1}{2}(1 + 1) = 1 \,\checkmark \]
计算 \(\langle\chi_1^0|\chi_1^0\rangle\)（三重态的归一化）
\[ \langle\chi_1^0|\chi_1^0\rangle = \frac{1}{2}(\langle\alpha_1\beta_2| + \langle\beta_1\alpha_2|)(\alpha_1\beta_2 + \beta_1\alpha_2) \]
\[ = \frac{1}{2}(1 + 1) = 1 \,\checkmark \]
计算正交性 \(\langle\chi_0|\chi_1^0\rangle\)
\[ \langle\chi_0|\chi_1^0\rangle = \frac{1}{2}(\langle\alpha_1\beta_2| - \langle\beta_1\alpha_2|)(\alpha_1\beta_2 + \beta_1\alpha_2) \]
\[ = \frac{1}{2}(\langle\alpha_1\beta_2|\alpha_1\beta_2\rangle + \langle\alpha_1\beta_2|\beta_1\alpha_2\rangle - \langle\beta_1\alpha_2|\alpha_1\beta_2\rangle - \langle\beta_1\alpha_2|\beta_1\alpha_2\rangle) \]
\[ = \frac{1}{2}(1 + 0 - 0 - 1) = 0 \,\checkmark \]
写出单粒子基态波函数（高斯波函数）
定义 \(\alpha_0 = \sqrt{m\omega/\hbar}\)，则
\[ \psi_{000}(\mathbf{r}) = \left(\frac{\alpha_0^3}{\pi^{3/4}}\right)^{1/2} \exp\left(-\frac{\alpha_0^2}{2}r^2\right) = \psi_0(x)\psi_0(y)\psi_0(z) \]
写出单粒子第一激发态（沿x方向）
\[ \psi_{100}(\mathbf{r}) = \psi_1(x)\psi_0(y)\psi_0(z) \]
其中 \(\psi_1(x) = \sqrt{2}\alpha_0\, x\cdot \psi_0(x)\)。
构造空间波函数（对称与反对称）
由于两电子不处于同一单粒子态，可以构造对称和反对称的空间部分：
\[ \Psi_{\text{sym}}(\mathbf{r}_1,\mathbf{r}_2) = \frac{1}{\sqrt{2}}[\psi_{000}(\mathbf{r}_1)\psi_{100}(\mathbf{r}_2) + \psi_{100}(\mathbf{r}_1)\psi_{000}(\mathbf{r}_2)] \]
\[ \Psi_{\text{asym}}(\mathbf{r}_1,\mathbf{r}_2) = \frac{1}{\sqrt{2}}[\psi_{000}(\mathbf{r}_1)\psi_{100}(\mathbf{r}_2) - \psi_{100}(\mathbf{r}_1)\psi_{000}(\mathbf{r}_2)] \]
匹配自旋对称性
由于费米子总波函数必须反对称，而 \(\Psi_{\text{total}} = \Psi_{\text{spatial}} \times \Psi_{\text{spin}}\)：
自旋单态 (S=0)：自旋部分反对称 \(\chi_0\)，空间部分必须对称 \(\Psi_{\text{sym}}\)。
自旋三重态 (S=1)：自旋部分对称 \(\chi_1\)，空间部分必须反对称 \(\Psi_{\text{asym}}\)。
写出完整波函数
单线态：
\[ \Psi_{\text{singlet}} = \Psi_{\text{sym}}(\mathbf{r}_1,\mathbf{r}_2) \times \frac{1}{\sqrt{2}}(\alpha_1\beta_2 - \beta_1\alpha_2) \]
三线态：
\[ \Psi_{\text{triplet}} = \Psi_{\text{asym}}(\mathbf{r}_1,\mathbf{r}_2) \times \begin{cases} \alpha_1\alpha_2 \\ \frac{1}{\sqrt{2}}(\alpha_1\beta_2 + \beta_1\alpha_2) \\ \beta_1\beta_2 \end{cases} \]
点选项即时判分，自动保存进度。多选题选完点「提交答案」。