气体分子运动论（Kinetic Theory of Gases） 通过统计力学方法，从分子运动的微观图像推导出宏观热力学性质，是连接微观粒子世界与宏观可观测现象的核心桥梁。
分子动理论（Kinetic Molecular Theory, KMT）是统计物理学的先驱，它假设物质由大量不断做无规则运动的粒子（分子或原子）组成，并利用统计平均的方法解释和预测宏观热力学量（如压强、温度、内能等）。该理论的核心思想是：宏观可观测量是微观粒子运动的统计平均结果。
| 时间 |
科学家 |
贡献 |
| 1738 |
丹尼尔·伯努利（Daniel Bernoulli） |
提出气体压强源于分子碰撞器壁的《流体动力学》思想 |
| 1857 |
鲁道夫·克劳修斯（Rudolf Clausius） |
引入"平均自由程"概念，建立气体动理论基本方程 |
| 1859 |
詹姆斯·克拉克·麦克斯韦（James Clerk Maxwell） |
导出分子速度分布律（麦克斯韦分布） |
| 1872 |
路德维希·玻尔兹曼（Ludwig Boltzmann） |
建立玻尔兹曼方程，给出 H 定理，奠定统计力学基础 |
| 1905 |
阿尔伯特·爱因斯坦（Albert Einstein） |
从分子动理论推导布朗运动理论，为分子存在提供定量证据 |
| 1908 |
让·佩兰（Jean Perrin） |
实验验证爱因斯坦的布朗运动理论，确定阿伏伽德罗常数 |
关键洞察：分子动理论不仅解释了理想气体状态方程 pV=nRT，更重要的是它将温度解释为分子平均动能的量度，将熵解释为系统微观状态数的对数，从而统一了力学与热学。
- 大量粒子假设：气体由数量极多的分子组成（1 mol 气体约有 6.022×1023 个分子），粒子的集体行为服从统计规律。
- 质点假设：分子本身的体积远小于气体所占体积，可忽略不计（分子直径 ∼0.3 nm，在标准状态下分子间平均距离 ∼3.3 nm，体积比约为 1:1000）。
- 无相互作用假设：分子间除弹性碰撞外无相互作用力（忽略分子间引力和斥力）。
- 完全弹性碰撞假设：分子之间、分子与器壁之间的碰撞是完全弹性的，动能守恒。
实际气体修正：真实气体不满足假设 2 和 3，范德瓦耳斯（van der Waals）通过引入分子体积修正项 b 和分子间引力项 a，将理想气体状态方程修正为 (p+a/Vm2)(Vm−b)=RT（见第六节）。
为了从微观运动推导宏观性质，我们需要两条基本公设：
- 各态历经假说（Ergodic Hypothesis）：经过足够长时间，系统会遍历所有可及的微观状态。
- 等概率原理：在平衡态下，系统处于每个可及微观状态的概率相等。
这两条公设使我们能够用系综平均替代时间平均，从而用统计方法处理大量粒子的运动。
考虑一个边长为 L 的立方体容器，内有 N 个质量为 m 的气体分子。每个分子的速度可以分解为三个分量：vx,vy,vz。
推导过程（一维简化）：
取分子 i 在 x 方向的分运动：
- 分子从右壁反弹时，动量变化为：Δpx=(−mvix)−(mvix)=−2mvix
- 器壁受到的动量变化为反作用：+2mvix
- 分子在 x 方向两次碰撞同一器壁的时间间隔：Δt=2L/vix
- 单位时间内分子 i 对器壁的平均作用力：Fi=Δt2mvix=Lmvix2
- N 个分子对器壁的总力：F=∑Lmvix2=Lm∑vix2
- 压强 p=L2F=L3m∑vix2=Vm∑vix2
注意 ∑vix2=31∑vi2（各向同性假设，三个方向的平方速度平均值相等），因此：
p=31Vm∑vi2=31VNmv2
其中 v2 是分子速度平方的平均值。定义方均根速度 vrms=v2，则有：
pV=31Nmvrms2
对比宏观理想气体状态方程 pV=NkBT，可得：
31Nmvrms2=NkBT
⇒21mvrms2=23kBT
这是分子动理论的核心结果： 气体分子的平均平动动能与热力学温度成正比。
Ek=21mv2=23kBT
已知：氮气 N2 在 T=300 K，分子质量 m=28×1.66×10−27=4.65×10−26 kg
计算：
方均根速度：vrms=m3kBT=4.65×10−263×1.38×10−23×300≈517 m/s
这个速度约 1860 km/h，接近声速的 1.5 倍！
| 气体 |
分子量 (g/mol) |
vrms 在 300K (m/s) |
与声速比 |
| H2 |
2 |
1930 |
5.6 倍 |
| He |
4 |
1365 |
4.0 倍 |
| N2 |
28 |
517 |
1.5 倍 |
| O2 |
32 |
483 |
1.4 倍 |
| CO2 |
44 |
412 |
1.2 倍 |
实际意义：分子热运动速度远大于日常生活中感受到的气体流速（如风速通常为 0-30 m/s），这说明气体中分子的定向流动速度远远小于其随机热运动速度。
麦克斯韦在 1859 年推导出，在热平衡状态下，气体分子的速度分布函数为：
f(v)dv=4πn(2πkBTm)3/2v2exp(−2kBTmv2)dv
其中 f(v)dv 表示速度在 v 到 v+dv 之间的分子数。
从分布函数可以定义三种具有物理意义的特征速度：
① 最概然速度 vp（分布函数的峰值位置）：
vp=m2kBT=M2RT
② 平均速度 vˉ：
vˉ=∫0∞vf(v)dv/∫0∞f(v)dv=πm8kBT=πM8RT
③ 方均根速度 vrms：
vrms=m3kBT=M3RT
三者之比为：vp:vˉ:vrms=1:1.128:1.225
以氮气为例，计算 T=300 K 时的三种特征速度：
| 特征速度 |
公式 |
数值 (m/s) |
| 最概然速度 vp |
M2RT |
422 |
| 平均速度 vˉ |
πM8RT |
476 |
| 方均根速度 vrms |
M3RT |
517 |
下图显示了不同温度下氧气分子速度分布的变化趋势（可用文字描述代替实际图像）：
分布曲线示意（氧气 $O_2$）：
f(v)
↑
| T=100K (v_p≈228 m/s) — 陡峭尖锐
| T=300K (v_p≈395 m/s) — 较宽
| T=600K (v_p≈559 m/s) — 宽而平缓
| .
| . . .
| . . . .
| . . . .
| . . . .
+----------------------------------→ v (m/s)
200 400 600 800 1000
关键观察：
- 温度升高，分布曲线变宽、峰值降低、向右移动
- 分布曲线下的总面积（总分子数）保持不变
- 温度越高，高能分子的比例越大——这是化学反应速率随温度升高的微观原因
行星能否保持大气取决于分子速度与逃逸速度的比值。以地球为例：
- 地球逃逸速度 vesc=11.2 km/s
- 氢气在 1000K 时 vrms≈3.5 km/s
- 虽然平均速度小于逃逸速度，但麦克斯韦分布尾部的分子速度可能超过 vesc
经验法则：当 vrms≥vesc/6 时，气体会在地质时间尺度上逃逸。地球的氢气 (vrms/vesc≈0.31) 因此会逐渐逃逸，而氧气和氮气则被保留。
| 行星 |
逃逸速度 (km/s) |
大气主要成分 |
保持能力 |
| 地球 |
11.2 |
N2,O2 |
保重重分子 |
| 月球 |
2.4 |
几乎无大气 |
全部逃逸 |
| 火星 |
5.0 |
稀薄 CO2 |
大部分逃逸 |
| 木星 |
59.5 |
H2,He |
保持轻分子 |
分子在两次连续碰撞之间自由飞行的平均距离称为平均自由程 λ，单位时间内一个分子与其他分子碰撞的平均次数称为碰撞频率 Z。
考虑一个半径为 r 的分子以相对速度 vˉrel 运动，其他分子视为静止的靶粒子：
有效碰撞截面：σ=πd2，其中 d=2r 是分子有效直径。
在一个时间 Δt 内，运动的分子扫过的碰撞体积为 Vcoll=σvˉrelΔt，其中的分子数为 nVcoll（n 是数密度）。
因此碰撞频率：
Z=nσvˉrel
考虑相对运动修正 vˉrel=2vˉ，可得：
Z=2nπd2vˉ
平均自由程：
λ=Zvˉ=2nπd21
代入 n=p/kBT：
λ=2πd2pkBT
条件：氮气在标准状态（T=273 K，p=1.013×105 Pa），分子直径 d≈0.37 nm
计算：
数密度：n=kBTp=1.38×10−23×2731.013×105≈2.69×1025 m−3
平均自由程：λ=2×π×(3.7×10−10)2×2.69×10251≈6.5×10−8 m=65 nm
碰撞频率：Z≈6.9×109 s−1
这意味着每个氮气分子每秒与其他分子碰撞约 70 亿次！
| 压强 |
λ (氮气, 273K) |
物理场景 |
| 105 Pa （常压） |
6.5×10−8 m |
日常环境 |
| 102 Pa （低真空） |
6.5×10−5 m |
真空管 |
| 1 Pa （高真空） |
6.5×10−3 m |
电子显微镜 |
| 10−4 Pa （超高真空） |
65 m |
粒子加速器 |
实际应用：真空技术的核心目标就是增大平均自由程，使粒子（电子、离子）在不受气体分子干扰的情况下长距离运动。电子显微镜所需的真空度正是基于这一原理。
能量均分定理（Equipartition Theorem）：在热平衡状态下，分子的每个自由度平均分配 21kBT 的能量。
| 运动形式 |
自由度数量 |
平均能量 |
| 平动 |
3 |
23kBT |
| 转动（双原子分子） |
2 |
kBT |
| 振动（双原子分子, 低T） |
2（动能+势能） |
kBT（高T时） |
根据自由度数量，可计算摩尔定容热容 CV,m：
| 分子类型 |
自由度 |
摩尔内能 Um |
CV,m |
| 单原子（He, Ne, Ar） |
3（平动） |
23RT |
23R |
| 双原子（N2,O2） |
5（3平动+2转动，室温） |
25RT |
25R |
| 双原子（高温） |
7（+2振动） |
27RT |
27R |
| 多原子（CO2,CH4） |
6（3平动+3转动） |
3RT |
3R |
理论值：C_V,m = \frac52R = \frac52 \times 8.314 = 20.785\ \textJ/(mol·K)
实验值（室温）：C_V,m \approx 20.74\ \textJ/(mol·K)
相对误差：0.2% —— 理论与实验高度吻合！
温度对氮气 $C_{V,m}$ 的影响：
C_{V,m} (J/mol·K)
↑
29.0 - · · (振动激活)
27.0 -
25.0 -
23.0 - · · · (转动激活完成)
21.0 - ·
19.0 - · (只有平动)
+-----------------------------------→ T (K)
50 200 500 1000 3000 5000
关键点：在低温下，某些自由度会被"冻结"——量子效应导致能量量子化，当 kBT 小于能级间隔时，该自由度无法被激发。这正是经典物理无法解释而需要量子统计的原因。
输运现象（Transport Phenomena）是当系统偏离平衡态时，分子通过碰撞和运动传递质量、动量和能量的过程。三种基本输运过程在分子层次上有统一的图像。
菲克第一定律：J=−Ddxdn
其中 J 是粒子流密度，D 是扩散系数。
从分子动理论推导的扩散系数：
D=31vˉλ
数值例：
CO2 在空气中的扩散（T=298 K, p=1 atm）：
vˉ≈380 m/s，λ≈40 nm
D=31×380×4×10−8≈1.6×10−5 m2/s
实验值约 1.6×10−5 m2/s，吻合良好。
| 扩散体系 |
扩散系数 (m2/s) |
物理含义 |
| 气体-气体（常压） |
10−5∼10−4 |
几秒内扩散数厘米 |
| 液体-液体 |
10−9∼10−8 |
比气体慢 10000 倍 |
| 固体-固体 |
10−14∼10−10 |
慢到几乎静止 |
实际应用：扩散是物质混合的根本机制。一杯咖啡不加搅拌时，糖分子的扩散约需数月才能分布到整杯——这就是为什么我们需要搅拌。
黏度 η 的分子动理论表达式：
η=31ρvˉλ
其中 ρ=nm 是密度。
关键预言：气体的黏度与压强无关！这由麦克斯韦首先预言，在当时极为反直觉。
验证实验：在 1 atm 和 0.1 atm 下测量气体黏度，结果几乎相同。麦克斯韦因此信服了分子动理论的正确性。
| 气体 |
η (\mu\textPa·s, 300K) |
温度依赖性 |
| 空气 |
18.5 |
∝T0.7 |
| N2 |
17.9 |
∝T0.7 |
| O2 |
20.8 |
∝T0.7 |
| He |
20.0 |
∝T0.6 |
| H2 |
8.9 |
∝T0.7 |
与液体的关键区别：气体黏度随温度升高而增大（分子运动加快，动量输运增强），而液体黏度随温度升高而减小（分子间作用力减弱）。
热导率 κ 的表达式：
κ=31ρvˉλcv
其中 cv 是单位质量的定容热容。
三种输运过程的统一形式：
| 输运过程 |
输运量 |
系数公式 |
宏观定律 |
| 扩散 |
粒子数 |
D=31vˉλ |
菲克定律 |
| 黏滞 |
动量 |
η=31ρvˉλ |
牛顿黏性定律 |
| 热传导 |
能量 |
κ=31ρvˉλcv |
傅里叶定律 |
所有三个系数都由 vˉλ 决定，体现了分子动理论的统一性。
1827 年，植物学家罗伯特·布朗（Robert Brown）在显微镜下观察到悬浮在水中的花粉颗粒做永不停息的无规则运动。当时他排除了活体生物、水流扰动、蒸发等各种可能原因，确定这是一种物理现象。
爱因斯坦在 1905 年发表了关于布朗运动的论文，从分子动理论角度给出了定量解释：
- 液体分子做热运动，不断撞击悬浮颗粒
- 由于分子运动的随机性，颗粒两侧受到的分子撞击力不平衡
- 这种净"推"力导致颗粒做随机运动
爱因斯坦的扩散公式：
颗粒的均方位移与时间成正比：
(Δx)2=2Dt
其中 D 是扩散系数。对于球形颗粒（半径 a），斯托克斯-爱因斯坦关系给出：
D=6πηakBT
将两个公式结合，得到可直接测量的关系：
(Δx)2=3πηakBTt
让·佩兰通过精心设计的实验，测量了布朗颗粒的位移分布：
实验数据（佩兰 1908，藤黄微粒，直径 0.53 μm）：
| 时间间隔 Δt (s) |
实测均方位移 (μm2) |
理论值 (μm2) |
| 30 |
112.5 |
113.0 |
| 60 |
223.8 |
226.0 |
| 90 |
340.1 |
339.0 |
| 120 |
452.9 |
452.0 |
结果：佩兰从这些数据计算出的阿伏伽德罗常数 NA=6.8×1023 mol−1，与现代值 6.022×1023 mol−1 在实验误差范围内吻合。这为分子的真实存在提供了无可辩驳的证据，爱因斯坦也因此获 1921 年诺贝尔物理学奖，佩兰获 1926 年诺贝尔物理学奖。
理想气体模型在某些条件下会失效：
- 高压：分子间距减小，分子体积不可忽略
- 低温：分子动能减小，分子间引力变得显著
- 相变点附近：分子间相互作用导致液化
范德瓦耳斯（van der Waals）在 1873 年对理想气体方程做了两项修正：
(p+Vm2a)(Vm−b)=RT
其中：
- b 是分子自身体积修正（b≈4NA⋅34πr3，约为分子本身体积的 4 倍）
- a 是分子间引力修正（a 越大，分子间引力越强）
体积修正 b：Vm−b 表示分子可自由运动的"有效体积"
压强修正 a/Vm2：器壁附近的气体分子受到内部分子的吸引力，减少了撞击器壁的动量，因此实际测得压强小于理想气体压强
| 气体 |
a (Pa·m6/mol2) |
b (×10−6 m3/mol) |
分子间引力强度 |
| He |
0.00341 |
23.7 |
极弱 |
| H2 |
0.0248 |
26.6 |
弱 |
| N2 |
0.137 |
38.7 |
中等 |
| O2 |
0.138 |
31.9 |
中等 |
| CO2 |
0.364 |
42.7 |
强 |
| H2O |
0.553 |
30.5 |
很强 |
范德瓦耳斯方程不仅修正了理想气体的偏差，还能预测气液相变。在临界温度 Tc 以下，方程给出的 p−V 曲线出现"凹凸"，对应气液共存区的相变。通过 ∂p/∂V=0 和 ∂2p/∂V2=0，可得临界点参数：
Tc=27Rb8a,pc=27b2a,Vm,c=3b
| 现象 |
分子动理论解释 |
| 闻到香水味 |
扩散——香味分子从高浓度区向低浓度区运动 |
| 气球受热膨胀 |
温度升高→分子速度增大→碰撞器壁更频繁→压强增大→体积膨胀 |
| 高压锅更快煮熟食物 |
压强增大→沸点升高→更高温度烹饪 |
| 摩擦生热 |
机械能通过分子碰撞转化为分子无规则运动（热能） |
| 乒乓球压扁后用热 |
温度升高→分子运动加剧→球内气体膨胀→将凹陷处推出水复原 |
- 气象学：大气压强随高度变化来源于空气分子受重力和热运动平衡
- 半导体：气体分子在真空室中的运动决定薄膜沉积质量
- 航空航天：高空中分子平均自由程增大，影响飞行器表面的气动加热
- 纳米技术：纳米尺度下，分子热运动成为主导因素（如 DNA 的分子动力学模拟）
┌──────────────────────────────────────────────────────┐
│ 热力学 │
│ (macroscopic, phenomenological) │
├──────────────────────────────────────────────────────┤
│ ↑ 统计平均 ↓ (桥梁) │
├──────────────────────────────────────────────────────┤
│ 分子动理论 │
│ (微观模型 + 统计方法 → 宏观可观测) │
├──────────────────────────────────────────────────────┤
│ ↑ 基础 ↓ │
├──────────────────────────────────────────────────────┤
│ 统计力学 / 量子统计 │
│ (一般理论框架, 涵盖更广泛系统) │
├──────────────────────────────────────────────────────┤
│ ↑ 微观基础 ↓ │
├──────────────────────────────────────────────────────┤
│ 经典力学 / 量子力学 (微观粒子的运动规律) │
└──────────────────────────────────────────────────────┘
| 物理量 |
公式 |
含义 |
| 压强微观表达式 |
p=31nmv2 |
压强 = 分子动量传递率 |
| 温度微观意义 |
21mv2=23kBT |
温度 = 平均动能 |
| 方均根速度 |
vrms=m3kBT |
分子热运动特征速度 |
| 平均自由程 |
λ=2πd2pkBT |
分子自由飞行距离 |
| 扩散系数 |
D=31vˉλ |
质量输运能力 |
| 黏度 |
η=31ρvˉλ |
动量输运能力 |
| 热导率 |
κ=31ρvˉλcv |
能量输运能力 |
| 能量均分 |
E=2fkBT |
每个自由度分配 21kBT |
| 范德瓦耳斯方程 |
(p+a/Vm2)(Vm−b)=RT |
真实气体修正 |
- 汪志诚，《热力学·统计物理》（第5版），高等教育出版社
- F. Reif, Fundamentals of Statistical and Thermal Physics, Waveland Press
- R. Baierlein, Thermal Physics, Cambridge University Press
- K. Huang, Statistical Mechanics, 2nd ed., Wiley
- Einstein, A. (1905). "On the Motion of Small Particles Suspended in a Stationary Liquid". Annalen der Physik
- Perrin, J. (1909). "Brownian Movement and Molecular Reality". Annales de Chimie et de Physique