古代原子思想与现代物理萌芽
在人类探索物质本质的漫长征途中，对“不可分割单元”的思辨早已超越朴素哲学，演变为精密科学。本章将系统追溯从古希腊的原子猜想，到19世纪化学与物理学的实证突破，最终聚焦于20世纪初那场颠覆性的实验革命——电子的发现与原子核式结构的确立。这不仅是物理学史的关键转折，更是整个粒子物理学大厦的奠基石。我们将深入核心实验的原始设计与数据分析，严格推导卢瑟福散射微分截面公式，并通过Python数值模拟，直观展现原子核电荷数对散射角分布的决定性影响。同时，剖析黑体辐射与光电效应如何撕开经典物理的坚固防线，迫使量子概念横空出世，为后续标准模型的孕育铺设必然之路。
---
**一、哲学猜想与科学实证的千年回响**
**1.1 德谟克利特的“原子”：不可分的哲学实在**
约公元前5世纪，古希腊哲学家留基伯与德谟克利特提出，宇宙万物由一种不可再分、永恒不变的微小粒子——“原子”（ἄτομος，意为“不可分割”）构成。这些原子在形状、大小和质量上存在差异，在虚空中永恒运动，通过碰撞与结合形成千变万化的物质。这一思想的核心在于机械决定论与实在论：物质世界是原子的机械组合，其性质由原子本身的物理属性决定。然而，受限于时代，这一构想缺乏任何实验依据，仅是对世界统一性与多样性的哲学解释，与同时期的四元素说（水、火、土、气）处于同等的思辨层面。它如同一颗深埋的种子，等待科学方法的春雨唤醒。
**1.2 道尔顿的复兴：化学原子论的实证基石**
1803年，英国化学家约翰·道尔顿基于质量守恒定律、定比定律和倍比定律等化学实验事实，提出了首个具有现代科学意义的原子理论。其核心要点包括： 1. 元素由特定原子构成：同种元素的原子在质量和性质上完全相同，不同元素原子则不同。 2. 化合物是原子的机械组合：化合物由不同元素原子以简单整数比结合而成。 3. 化学反应是原子的重新排列：反应中原子种类、数目不变，仅改变结合方式。
道尔顿理论最革命性的贡献在于将原子从一个定性哲学概念转变为一个定量科学实体，其“原子量”的提出首次赋予了原子可测量、可比较的物理属性。例如，通过对水（H₂O）和氨（NH₃）等化合物组成比例的精确分析，他推导出氢、氧、氮等元素的相对原子质量。尽管他错误地假设了简单化合物的最简式（如认为水是HO），但其框架为后续化学的发展提供了坚实的语言和逻辑基础。化学原子论的成功，强烈暗示了原子作为物质结构基本单元的真实性，但原子的内部结构、是否可分，依然是巨大的未知。
---
**二、19世纪末的物理学晴空与两朵“乌云”**
**2.1 经典物理学的宏伟统一**
19世纪下半叶，物理学呈现出近乎完美的理论图景：
- **电磁学的大一统**：詹姆斯·克拉克·麦克斯韦于1860年代将电、磁、光现象统一于一组优雅的偏微分方程——麦克斯韦方程组。他预言了电磁波的存在，并证明光就是一种电磁波，其传播速度 $c$ 由真空介电常数 $\epsilon_0$ 和磁导率 $\mu_0$ 决定：$c = 1/\sqrt{\epsilon_0 \mu_0}$。这一定律的预言与赫兹1887年的实验验证，标志着经典物理在描述宏观电磁现象上的终极成功。
- **热力学的成熟**：热力学第一定律（能量守恒）与第二定律（熵增原理）确立了能量转换与过程方向性的普遍法则。统计力学（玻尔兹曼、吉布斯）则从微观粒子运动的角度成功推导出宏观热力学定律，将概率概念引入物理学。
- **光学的波动说胜利**：杨氏双缝实验与菲涅尔的衍射理论彻底击败了牛顿的微粒说，光的波动性成为共识。
此时，物理学家普遍相信，所有自然现象均可归结为连续、确定、可预测的经典力学与电磁学方程。正如开尔文勋爵1900年所言，物理学天空中“只有两朵小小的乌云”：一是迈克尔逊-莫雷实验对“以太”风的零结果，二是黑体辐射的“紫外灾难”。正是这两朵乌云，最终引发了席卷整个物理学根基的量子风暴。
**2.2 “以太”的困境与实验物理学的勃兴**
“以太”作为光波传播的绝对介质，其存在本身即隐含了一个可被探测的“绝对参考系”。阿尔伯特·迈克尔逊与爱德华·莫雷于1887年设计的干涉仪实验，旨在测量地球相对于以太的运动速度。其结果是革命性的：未观测到任何显著的干涉条纹移动。这一“零结果”直接否定了以太漂移的存在，为爱因斯坦1905年提出狭义相对论（抛弃以太，确立光速不变与相对性原理）埋下了伏笔。更重要的是，该实验代表了精密测量在物理学中的核心地位日益凸显，为后续一系列关键实验（如汤姆逊、卢瑟福实验）树立了方法论典范：用可重复、可量化的实验数据，裁决理论之争。
---
**三、电子的发现：阴极射线与J.J.汤姆逊的实验**
**3.1 阴极射线的争议：粒子还是波？**
19世纪中叶，研究者发现，在低压气体放电管中，阴极会发出一种神秘的射线（阴极射线），能使管对面玻璃发出荧光，并可被磁场偏转。其本质引发激烈争论：
- **“粒子说”**：认为阴极射线是带负电的粒子流（“微粒”）。
- **“波说”**：认为其本质是某种在以太中传播的电磁波（类似光）。
关键分歧在于荷质比（e/m）。若为粒子，其 $e/m$ 应远大于氢离子的 $e/m$（因为 $e$ 相同，$m$ 应极小）；若为波，则无需讨论此量。
**3.2 汤姆逊的突破性实验（1897年）**
J.J.汤姆逊设计了一系列精妙的实验，使用改进的克鲁克斯管， definitive 地证明了阴极射线的粒子性并测定了其 $e/m$。
实验装置与原理： 1. 真空度：将放电管抽至高真空（$<10^{-6}$ mmHg），消除气体分子干扰，使阴极射线成为在真空中传播的准直束。 2. 电场偏转：在射线路径上施加已知强度的横向静电场（由平行板电容器产生）。射线受洛伦兹力 $F_e = eE$ 发生偏转。通过测量偏转距离，可计算 $e/m$。 3. 磁场偏转：单独施加已知强度的横向静磁场（由电磁铁产生）。射线受洛伦兹力 $F_m = evB$ 做圆周运动，偏转半径 $r = mv/(eB)$。通过测量 $r$ 和已知的 $v$（或结合电场测量），同样可计算 $e/m$。 4. 电场与磁场平衡法：同时施加 $E$ 和 $B$，调节至射线不发生偏转（即 $eE = evB$），此时 $v = E/B$。将此速度代入磁场偏转公式 $r = mv/(eB)$，即可直接解出 $e/m = E/(B^2 r)$，无需预先知道 $v$ 或 $e$、$m$ 的单独值。
核心数据与结论： 汤姆逊测得阴极射线粒子的 $e/m \approx 1.7 \times 10^{11}$ C/kg，比当时已知的最小离子（氢离子）的 $e/m$ 大约1800倍。由于电荷 $e$ 被认为基本单位（与电解定律相符），这意味着该粒子的质量 $m$ 远小于氢原子质量（约1/1800）。他据此宣布发现了一种普遍存在的、比原子更基本的带负电粒子，并将其命名为“微粒”（corpuscle），后称电子（electron）。
历史意义：
- **首次揭示亚原子结构**：电子是人类发现的第一个基本粒子，证明原子是可分的。
- **确立“普适性”**：无论阴极材料（金属、气体）或射线来源（光电效应、热辐射）如何，所产生粒子的 $e/m$ 恒定，表明电子是所有物质的共同组成部分。
- **方法论典范**：通过**交叉验证**（电场、磁场、平衡法）和**精密测量**，解决了根本性质疑，为后续实验树立标杆。
---
**四、卢瑟福散射与原子核式模型：从“葡萄干布丁”到“太阳系”**
**4.1 汤姆逊模型与它的困境**
发现电子后，汤姆逊于1904年提出原子的“葡萄干布丁”模型：原子是一个带正电的、均匀分布的球体（“布丁”），电子（“葡萄干”）镶嵌其中，处于静电平衡。该模型能解释原子的电中性，但无法解释以下关键现象：
- **α粒子的大角度散射**：盖革和马斯登在卢瑟福指导下，用放射性源（如镭）发射的α粒子（氦原子核，$^{4}\text{He}^{2+}$，质量约7300倍电子）轰击极薄金箔（厚度约1000原子层）。按汤姆逊模型，由于正电荷均匀分布，α粒子穿过时只会经历微弱的、渐进的偏转（多级小角度散射），**大角度（>90°）散射概率应几乎为零**。然而实验观测到**约1/8000的α粒子发生超过90°的大角度偏转，甚至有的被反弹回来**。这是汤姆逊模型完全无法容纳的。
**4.2 卢瑟福的洞察与核式模型（1911年）**
欧内斯特·卢瑟福在分析盖革-马斯登实验数据后，做出了革命性推断：只有原子中绝大部分质量和正电荷都集中在一个极小的核心（原子核）内，才能产生如此强大的库仑排斥力，使重α粒子发生大角度偏转。他于1911年提出行星模型（核式模型）：
- 原子中心有一个**体积极小（直径约 $10^{-14}$ m，比原子 $10^{-10}$ m 小万倍）、质量极大（集中99.95%以上原子质量）、带正电**的原子核。
- 电子在核外广阔空间（$10^{-10}$ m 量级）中绕核运动，类似行星绕太阳。
- 原子内部绝大部分是**空旷空间**，这解释了为何绝大多数α粒子能直接穿过金箔（未靠近原子核）。
模型的成功： 1. 自然解释大角度散射：当α粒子非常接近原子核时，巨大的库仑斥力使其发生剧烈偏转。 2. 定量预言：卢瑟福基于单次碰撞的经典力学（忽略电子影响，因电子质量太小），推导出α粒子被原子核散射的微分截面公式，完美符合实验数据。
---
**五、卢瑟福散射微分截面公式的严格推导**
**5.1 物理图像与基本假设**
考虑一个质量为 $m_\alpha$、电荷为 $Z_1 e$（对α粒子，$Z_1=2$）的入射粒子，以初速度 $v$、瞄准参数（碰撞参数）$b$ 飞向一个质量为 $M$、电荷为 $Z_2 e$（对金核，$Z_2=79$）的静止靶核。假设： 1. 库仑力主导：相互作用仅为静电库仑力 $\vec{F} = \frac{Z_1 Z_2 e^2}{4\pi\epsilon_0 r^2} \hat{r}$，且为排斥力（$Z_1, Z_2 > 0$）。 2. 核静止：由于 $M \gg m_\alpha$（金核质量约α粒子的50倍），近似认为靶核固定不动。 3. 单次碰撞：忽略多次散射。 4. 无能量损失：纯库仑散射是弹性碰撞，动能守恒。
**5.2 经典轨道与散射角关系**
在中心力场中，粒子的轨迹是双曲线的一支，原子核位于双曲线的焦点。瞄准参数 $b$ 与散射角 $\theta$（出射方向与入射方向夹角）存在一一对应关系。
推导步骤： 1. 角动量守恒：$L = m_\alpha v b = \text{常数}$。 2. 能量守恒：$E = \frac{1}{2} m_\alpha v^2 = \text{常数}$。 3. 轨道方程：在中心力场中，轨道极坐标方程为： $$ u(\phi) = \frac{1}{r} = \frac{Z_1 Z_2 e^2}{4\pi\epsilon_0 L^2} [1 + e \cos(\phi - \phi_0)] $$ 其中 $u=1/r$，$\phi$ 是极角，$e$ 是轨道偏心率，$\phi_0$ 是相位常数。由初始条件（$r \to \infty$ 时，$\phi \to \pm\phi_\infty$，且 $dr/d\phi = \mp b$），可确定： $$ e = \sqrt{1 + \frac{2 E L^2}{m_\alpha (Z_1 Z_2 e^2 / 4\pi\epsilon_0)^2}} $$ 代入 $E$ 和 $L$，得： $$ e = \sqrt{1 + \left( \frac{2 b E}{Z_1 Z_2 e^2 / 4\pi\epsilon_0} \right)^2} = \sqrt{1 + \left( \frac{b v^2 m_\alpha}{Z_1 Z_2 e^2 / 4\pi\epsilon_0} \right)^2} $$
4. 散射角与瞄准参数关系：双曲线的渐近线方向与焦点连线的夹角 $\phi_\infty$ 满足 $\cos\phi_\infty = -1/e$。散射角 $\theta = \pi - 2\phi_\infty$，故： $$ \cot\frac{\theta}{2} = \frac{2 b E}{Z_1 Z_2 e^2 / 4\pi\epsilon_0} = \frac{b v^2 m_\alpha}{Z_1 Z_2 e^2 / 4\pi\epsilon_0} $$ 这是核心几何关系：瞄准参数 $b$ 越大，散射角 $\theta$ 越小。
**5.3 微分截面的定义与推导**
微分截面 $d\sigma/d\Omega$ 是散射物理学核心概念：它表示单位立体角内被散射的粒子数占入射粒子总数的比例，物理意义是有效靶面积。若瞄准参数在 $b$ 到 $b+db$ 之间的粒子，被散射到极角 $\theta$ 到 $\theta+d\theta$、方位角 $\phi$ 到 $\phi+d\phi$ 的立体角 $d\Omega = \sin\theta d\theta d\phi$ 内，则： $$ \frac{d\sigma}{d\Omega} db \cdot 2\pi = \frac{N(\theta)}{N_{\text{总}}} $$ （利用轴对称性，方位角积分得 $2\pi$）。因此： $$ \frac{d\sigma}{d\Omega} = \frac{b}{\sin\theta} \left| \frac{db}{d\theta} \right| $$
由 $\cot(\theta/2) = K b$，其中 $K = \frac{m_\alpha v^2}{Z_1 Z_2 e^2 / 4\pi\epsilon_0}$ 为常数，可得： $$ b = \frac{1}{K} \cot\frac{\theta}{2} $$ 求导： $$ \frac{db}{d\theta} = -\frac{1}{2K} \csc^2\frac{\theta}{2} $$ 取绝对值并代入： $$ \frac{d\sigma}{d\Omega} = \frac{ \frac{1}{K} \cot\frac{\theta}{2} }{ \sin\theta } \cdot \frac{1}{2K} \csc^2\frac{\theta}{2} $$ 利用三角恒等式 $\sin\theta = 2\sin(\theta/2)\cos(\theta/2)$，化简： $$ \frac{d\sigma}{d\Omega} = \frac{1}{4K^2} \cdot \frac{ \cot\frac{\theta}{2} \csc^2\frac{\theta}{2} }{ 2\sin\frac{\theta}{2}\cos\frac{\theta}{2} } = \frac{1}{4K^2} \cdot \frac{ \cos\frac{\theta}{2} / \sin\frac{\theta}{2} \cdot 1/\sin^2\frac{\theta}{2} }{ 2\sin\frac{\theta}{2}\cos\frac{\theta}{2} } = \frac{1}{4K^2} \cdot \frac{1}{4 \sin^4(\theta/2)} $$ 最终得到卢瑟福散射公式： $$ \boxed{ \frac{d\sigma}{d\Omega} = \left( \frac{Z_1 Z_2 e^2}{8\pi\epsilon_0 m_\alpha v^2} \right)^2 \frac{1}{\sin^4(\theta/2)} } $$
公式解读：
- 微分截面与 $\sin^4(\theta/2)$ 成反比，意味着**小角度散射占绝对主导**，大角度散射极其稀少（符合实验：1/8000）。
- 与 $Z_1^2 Z_2^2$ 成正比：核电荷数 $Z_2$ 越大，大角度散射概率急剧增加（金箔 $Z_2=79$ 产生显著大角度散射，若用铝箔 $Z_2=13$ 则几乎观测不到）。
- 与入射动能 $(\propto m_\alpha v^2)$ 成反比：能量越高，越不易偏转。
- **与电子完全无关**：公式中不包含电子质量或电荷，因为电子质量太小，对α粒子运动轨迹无显著影响。这从理论上解释了为何汤姆逊模型失败——它试图用电子云来解释大角度散射，但电子无法提供足够的反冲力。
---
**六、Python模拟：可视化卢瑟福散射角分布**
为直观展示卢瑟福公式的预言，我们编写Python代码，计算不同原子核电荷数 $Z_2$ 下，α粒子（$Z_1=2$）的微分截面随散射角 $\theta$ 的变化，并绘制分布图。
**6.1 模拟原理与代码**
模拟固定入射能量（即固定 $m_\alpha v^2$），计算 $d\sigma/d\Omega$ 的相对值（因常数因子不影响角分布形状）。关键步骤： 1. 定义物理常数（$e, \epsilon_0$）和粒子参数（$m_\alpha, Z_1$）。 2. 定义卢瑟福公式函数，输入 $\theta$（弧度）、$Z_2$，输出 $d\sigma/d\Omega$。 3. 生成 $\theta$ 从 $1^\circ$ 到 $180^\circ$ 的数组（避开 $0^\circ$ 奇点）。 4. 对几个典型 $Z_2$（如13-铝、29-铜、79-金）计算并绘图。 5. 使用对数坐标（因数值跨越多个量级）。
```python import numpy as np import matplotlib.pyplot as plt
物理常数 (SI单位)
e = 1.602176634e-19 # 基本电荷, C epsilon_0 = 8.8541878128e-12 # 真空介电常数, F/m m_alpha = 6.6446573357e-27 # α粒子质量 (≈4*u), kg Z1 = 2 # α粒子电荷数