1. 麦克斯韦速度分布率
麦克斯韦速率分布律
在一个容器中装有大量气体分子,它们不断碰撞、改变速度.单个分子的速率是随机的、不可预测的,但当分子数目足够大时,它们的集体行为却服从一个确定的统计分布.这个分布由麦克斯韦于 1859 年首先推导出来,是统计物理中最基本的结果之一.本页将从概率论的基本概念出发,逐步建立速度空间和分布函数的概念,推导麦克斯韦速度分布和速率分布,并讨论其物理意义和应用.
本页面仍在建设中
本页面的内容仍在建设中,参考图片有待补充,部分推导和公式可能不完整或有待修正.请谨慎参考,并欢迎向我们反馈问题.
学习目标
读完本页后,你应该能够:
- 理解统计规律与分布函数的概念
- 理解速度空间和速度分布函数的定义
- 从统计独立和各向同性两个假设出发,推导麦克斯韦速度分布
- 写出麦克斯韦速率分布函数,理解其曲线特征
- 计算最概然速率、平均速率和方均根速率,并理解三者的物理意义
- 用速率分布解释大气逃逸、泻流等现象
§1 概率论预备知识
本页的推导需要用到概率论的基本概念.下面做一个简要回顾;如需详细学习,请参阅概率论相关页面.
前置知识:概率论速览
随机变量:将随机实验的结果映射为实数的函数.分为离散型和连续型.
概率密度函数 (PDF):对于连续随机变量 𝑋
,概率密度函数 𝑓(𝑥)
满足:
𝑃(𝑎≤𝑋≤𝑏)=∫𝑏𝑎𝑓(𝑥)d𝑥
归一化条件:∫+∞−∞𝑓(𝑥) d𝑥 =1
.概率密度函数在整个定义域上的积分必须等于 1.
数学期望:𝐸[𝑋] =∫+∞−∞𝑥 𝑓(𝑥) d𝑥
,度量分布的 "中心" 位置.
高斯积分:∫+∞0𝑥𝑛 𝑒−𝛼𝑥2 d𝑥
可通过递推公式计算,最基本的结果是 ∫+∞−∞𝑒−𝛼𝑥2 d𝑥 =√𝜋𝛼
.
详细内容请参阅:概率论的基本概念、一维随机变量及其分布、随机变量的数字特征.
§2 统计规律与分布函数
从伽尔顿板说起
想象一个伽尔顿板(Galton board):大量小球从顶部下落,每经过一层钉板时随机向左或向右偏转.单个小球最终落入哪个槽位是随机的——你无法预测下一个小球的位置.但是,当大量小球落下后,它们在各个槽位中堆积形成的 直方图 却呈现出一个确定的、光滑的钟形曲线.
这就是 统计规律 的核心:单个事件是随机的,但大量事件的集体行为服从确定的规律.
两个等价的视角:
- 一个粒子多次测量:对同一个分子反复测量其速率,记录各次结果的分布
- 大量粒子一次测量:在某一时刻,对容器中所有分子的速率进行统计
在平衡态下,这两种视角给出相同的结果.
从直方图到分布函数
设小球落入位置 𝑥
附近宽度为 Δ𝑥𝑖
的第 𝑖
个槽中的数目为 Δ𝑁𝑖
.当槽的宽度 Δ𝑥 →0
且小球总数 𝑁 →∞
时,离散的直方图过渡为连续的 分布函数:
d𝑁𝑁=𝑓(𝑥)d𝑥
其中 𝑓(𝑥)
称为 分布函数(更准确地说,是概率密度函数).它表示随机变量取值在 𝑥
附近单位区间内的 相对概率密度.
分布函数必须满足 归一化条件:
∫+∞−∞𝑓(𝑥)d𝑥=1
这对应 "所有小球最终都落入某个槽中" 的确定性事实.
任意物理量 𝐺(𝑥)
的统计平均值为:
――𝐺=∫+∞−∞𝐺(𝑥)𝑓(𝑥)d𝑥
§3 速度空间与速度分布函数
速度空间
在经典力学中,每个分子的速度可以用三个分量 (𝑣𝑥,𝑣𝑦,𝑣𝑧)
来描述.这三个分量构成一个三维空间,称为 速度空间.速度空间中的每个点对应一种可能的分子速度状态.
在平衡态下,容器中 𝑁
个分子的速度分布在速度空间中形成一个 "云"——某些区域分子密集(速度出现的概率大),某些区域分子稀疏.
速度分布函数
设 𝑓(𝑣𝑥,𝑣𝑦,𝑣𝑧)
为速度分布函数,其物理含义是:在速度空间中,速度分量落在 (𝑣𝑥,𝑣𝑦,𝑣𝑧)
附近体积元 d𝑣𝑥 d𝑣𝑦 d𝑣𝑧
内的分子数占总分子数的比例为:
d𝑁𝑁=𝑓(𝑣𝑥,𝑣𝑦,𝑣𝑧)d𝑣𝑥d𝑣𝑦d𝑣𝑧
归一化条件要求:
∭+∞−∞𝑓(𝑣𝑥,𝑣𝑦,𝑣𝑧)d𝑣𝑥d𝑣𝑦d𝑣𝑧=1
球坐标形式
由于在平衡态下气体是各向同性的(没有特殊方向),速度分布只依赖于速率 𝑣 =√𝑣2𝑥+𝑣2𝑦+𝑣2𝑧
,而不依赖于方向.将速度空间转化为球坐标 (𝑣,𝜃,𝜑)
,体积元为 𝑣2sin𝜃 d𝑣 d𝜃 d𝜑
.对角度积分后(∫2𝜋0d𝜑∫𝜋0sin𝜃 d𝜃 =4𝜋
),得到:
d𝑁𝑁=4𝜋𝑣2𝑓(𝑣)d𝑣
其中 𝑓(𝑣)
是只依赖于速率的分布函数.
约束条件
分布函数 𝑓(𝑣)
必须满足两个物理约束:
粒子数守恒(归一化):∫+∞04𝜋𝑣2𝑓(𝑣) d𝑣 =1
平均动能与温度的关系:由能均分定理,∫+∞012𝑚𝑣2 ⋅4𝜋𝑣2𝑓(𝑣) d𝑣 =32𝑘𝑇
这两个约束将用于确定分布函数中的待定常数.
§4 麦克斯韦速度分布的推导
两个基本假设
麦克斯韦在推导速度分布时,作了以下两个假设:
麦克斯韦的两个假设
假设 1:速度分量统计独立. 在热平衡态下,三个速度分量 𝑣𝑥
、𝑣𝑦
、𝑣𝑧
的分布彼此独立.即:
𝑓(𝑣𝑥,𝑣𝑦,𝑣𝑧)=𝑓(𝑣𝑥)𝑓(𝑣𝑦)𝑓(𝑣𝑧)
假设 2:各向同性. 对于宏观上静止的气体,速度分布不依赖于方向,只依赖于速率 𝑣
.即:
𝑓(𝑣𝑥,𝑣𝑦,𝑣𝑧)=Φ(𝑣2𝑥+𝑣2𝑦+𝑣2𝑧)
推导过程
将两个假设结合:𝑓(𝑣𝑥) 𝑓(𝑣𝑦) 𝑓(𝑣𝑧) =Φ(𝑣2𝑥 +𝑣2𝑦 +𝑣2𝑧)
.
两边取对数:
ln𝑓(𝑣𝑥)+ln𝑓(𝑣𝑦)+ln𝑓(𝑣𝑧)=lnΦ(𝑣2𝑥+𝑣2𝑦+𝑣2𝑧)
对 𝑣2𝑥
求偏导(注意左边只有 ln𝑓(𝑣𝑥)
项依赖于 𝑣𝑥
):
dln𝑓(𝑣𝑥)d(𝑣2𝑥)=𝜕lnΦ𝜕(𝑣2𝑥+𝑣2𝑦+𝑣2𝑧)≡−𝐵
由于左边只依赖于 𝑣𝑥
,右边只依赖于 𝑣2𝑥 +𝑣2𝑦 +𝑣2𝑧
,要使等式对任意 (𝑣𝑥,𝑣𝑦,𝑣𝑧)
成立,两边必须等于同一个常数 −𝐵
.同理对 𝑣𝑦
、𝑣𝑧
求导也得到相同的常数.
因此:
ln𝑓(𝑣𝑖)=−𝐵𝑣2𝑖+ln𝐶𝑖⟹𝑓(𝑣𝑖)=𝐶𝑖𝑒−𝐵𝑣2𝑖
其中 𝐵
和 𝐶𝑖
为待定常数.将三个分量相乘:
𝑓(𝑣𝑥,𝑣𝑦,𝑣𝑧)=𝐶𝑥𝐶𝑦𝐶𝑧𝑒−𝐵(𝑣2𝑥+𝑣2𝑦+𝑣2𝑧)=𝐶3𝑒−𝐵𝑣2
确定常数
常数 𝐵
由平均动能条件确定.考虑一个速度分量(如 𝑣𝑥
),其平均动能为:
――――12𝑚𝑣2𝑥=12𝑘𝑇
利用一维高斯分布的结果:―――𝑣2𝑥 =12𝐵
,因此:
12𝑚⋅12𝐵=12𝑘𝑇⟹𝐵=𝑚2𝑘𝑇
常数 𝐶
由归一化条件确定.对 𝑣𝑥
分量:
∫+∞−∞𝑓(𝑣𝑥)d𝑣𝑥=𝐶𝑥∫+∞−∞𝑒−𝐵𝑣2𝑥d𝑣𝑥=𝐶𝑥√𝜋𝐵=1
𝐶𝑥=√𝐵𝜋=√𝑚2𝜋𝑘𝑇
由各向同性,𝐶𝑥 =𝐶𝑦 =𝐶𝑧
,因此 𝐶3 =(𝑚2𝜋𝑘𝑇)3/2
.
最终结果
麦克斯韦速度分布函数
𝑓(𝒗)=(𝑚2𝜋𝑘𝑇)3/2exp(−𝑚(𝑣2𝑥+𝑣2𝑦+𝑣2𝑧)2𝑘𝑇)
其中 𝑚
为分子质量,𝑘
为玻尔兹曼常数,𝑇
为热力学温度.
各分量的分布是一维高斯分布:
𝑓(𝑣𝑖)=√𝑚2𝜋𝑘𝑇exp(−𝑚𝑣2𝑖2𝑘𝑇),𝑖=𝑥,𝑦,𝑧
§5 麦克斯韦速率分布
从速度分布到速率分布
在物理实验中,我们通常更关心分子的 速率 𝑣 =|𝒗|
(速度的大小),而不是速度的方向.将速度分布函数对所有方向积分,得到 速率分布函数:
麦克斯韦速率分布函数
𝐹(𝑣)=4𝜋(𝑚2𝜋𝑘𝑇)3/2𝑣2exp(−𝑚𝑣22𝑘𝑇),𝑣≥0
𝐹(𝑣) d𝑣
表示速率在 [𝑣,𝑣 +d𝑣]
区间内的分子数占总分子数的比例.
推导:在球坐标中,速率在 [𝑣,𝑣 +d𝑣]
内的分子占据速度空间中一个球壳,其体积为 4𝜋𝑣2 d𝑣
.由于 𝑓(𝒗)
只依赖于速率,可以直接写出:
𝐹(𝑣)d𝑣=𝑓(𝑣)⋅4𝜋𝑣2d𝑣=4𝜋(𝑚2𝜋𝑘𝑇)3/2𝑣2𝑒−𝑚𝑣2/2𝑘𝑇d𝑣
速率分布曲线的特征
𝐹(𝑣)
曲线具有以下特征:
- 起点:𝐹(0) =0
.速率为零的分子数为零(𝑣2
因子的贡献). - 终点:𝐹(∞) =0
.速率极大的分子数也为零(指数衰减占主导). - 存在极大值:在某一速率 𝑣𝑝
处,𝐹(𝑣)
取极大值.这个速率称为 最概然速率. - 归一化:∫+∞0𝐹(𝑣) d𝑣 =1
.曲线下的总面积为 1.
曲线随参数的变化
温度的影响(𝑚
固定):当温度 𝑇
升高时——
- 最概然速率 𝑣𝑝
增大(峰值右移) - 峰值高度 𝐹(𝑣𝑝)
降低(峰值变矮) - 分布变宽(高速分子比例增大)
- 由于总面积始终为 1,曲线变得更加 "矮胖"
质量的影响(𝑇
固定):当分子质量 𝑚
增大时——
- 最概然速率 𝑣𝑝
减小(峰值左移) - 峰值高度 𝐹(𝑣𝑝)
升高(峰值变高) - 分布变窄
- 曲线变得更加 "高瘦"
数值例:氧气在 273 K 下的速率分布
下表给出了氧气分子(𝑀 =32 g/mol
,𝑇 =273 K
)在各速率区间内的分子百分比:
| 速率区间 (m/s) | 分子百分比 (%) |
|---|
| 0–100 | 1.4 |
| 100–200 | 8.1 |
| 200–300 | 16.7 |
| 300–400 | 21.4 |
| 400–500 | 20.5 |
| 500–600 | 15.1 |
| 600–700 | 9.2 |
| 700–800 | 4.8 |
| 800–900 | 2.0 |
| 900–1000 | 0.7 |
| > 1000 | 0.1 |
可以看到,大多数分子的速率集中在 200–600 m/s 范围内,其中 300–400 m/s 区间占比最大(对应最概然速率附近).
§6 特征速率
速率分布函数 𝐹(𝑣)
描述了分子速率的完整统计信息.在实际问题中,我们常用几个有明确物理意义的 特征速率 来概括分布的主要特征.
最概然速率 𝑣𝑝
最概然速率是速率分布函数 𝐹(𝑣)
取极大值时对应的速率.通过求解 d𝐹d𝑣 =0
:
dd𝑣[𝑣2𝑒−𝑚𝑣2/2𝑘𝑇]=0⟹2𝑣−𝑚𝑘𝑇𝑣3=0![\dfrac{\mathrm{d}}{\mathrm{d}v}\!\left[v^2\,e^{-mv^2/2kT}\right] = 0 \quad \Longrightarrow \quad 2v - \dfrac{m}{kT}\,v^3 = 0](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
解得(舍去 𝑣 =0
):
𝑣𝑝=√2𝑘𝑇𝑚=√2𝑅𝑇𝑀mol≈1.41√𝑅𝑇𝑀mol
其中 𝑀mol
为摩尔质量,𝑅
为普适气体常量.
物理意义:𝑣𝑝
是概率密度最大的速率——在 𝑣𝑝
附近单位速率区间内的分子比例最高.
平均速率 ¯𝑣
平均速率为所有分子速率的统计平均值:
¯𝑣=⟨𝑣⟩=∫+∞0𝑣𝐹(𝑣)d𝑣=4𝜋(𝑚2𝜋𝑘𝑇)3/2∫+∞0𝑣3𝑒−𝑚𝑣2/2𝑘𝑇d𝑣
利用高斯积分 𝐼3(𝑚/2𝑘𝑇) =12(𝑚/2𝑘𝑇)2
,计算得:
¯𝑣=√8𝑘𝑇𝜋𝑚=√8𝑅𝑇𝜋𝑀mol≈1.60√𝑅𝑇𝑀mol
方均根速率 𝑣rms
方均根速率为速率平方平均值的平方根,与分子的平均动能直接相关:
――𝑣2=∫+∞0𝑣2𝐹(𝑣)d𝑣=3𝑘𝑇𝑚
𝑣rms=√――𝑣2=√3𝑘𝑇𝑚=√3𝑅𝑇𝑀mol≈1.73√𝑅𝑇𝑀mol
验证:一个分子的平均平动动能为 ―――𝐸𝑘 =12𝑚――𝑣2 =12𝑚 ⋅3𝑘𝑇𝑚 =32𝑘𝑇
,与能均分定理一致.
速度分量的方均根
对于单个速度分量(如 𝑣𝑥
),其分布是一维高斯分布,方均根速率为:
√―――𝑣2𝑥=√𝑘𝑇𝑚
特征速率汇总
特征速率公式
| 特征速率 | 公式 | 与 √𝑅𝑇/𝑀mol 的关系 |
|---|
最概然速率 𝑣𝑝 | √2𝑘𝑇𝑚 | ≈1.41√𝑅𝑇𝑀mol |
平均速率 ¯𝑣 | √8𝑘𝑇𝜋𝑚 | ≈1.60√𝑅𝑇𝑀mol |
方均根速率 𝑣rms | √3𝑘𝑇𝑚 | ≈1.73√𝑅𝑇𝑀mol |
三者满足不等式:𝑣𝑝 <¯𝑣 <𝑣rms
这个不等式的直观理解:𝑣2
对大速率赋予更大的权重,因此 𝑣rms >¯𝑣
;而 𝐹(𝑣)
的不对称性(右侧尾巴更长)使得 ¯𝑣 >𝑣𝑝
.
§7 应用
大气逃逸
行星大气中的气体分子是否能逃逸到太空,取决于分子的速率是否超过 逃逸速度 𝑣esc
.对于地球,𝑣esc ≈11.2 km/s
.
虽然平均速率远小于逃逸速度,但速率分布的高速 "尾巴" 中总有一些分子的速率超过 𝑣esc
.这些分子可以克服引力逃逸到太空,导致大气逐渐散失.
逃逸速率与最概然速率的比值 𝑣esc/𝑣𝑝
决定了大气逃逸的快慢:
- 若 𝑣esc ≫𝑣𝑝
(如地球上的氮气、氧气),逃逸极慢,大气可以长期保持 - 若 𝑣esc ≳𝑣𝑝
(如地球上的氢气、月球上的所有气体),逃逸较快,大气难以保持
这就是为什么地球大气中几乎没有氢气和氦气,而月球几乎没有大气.
泻流(Effusion)
当容器壁上有一个小孔(孔径远小于分子平均自由程)时,分子通过小孔逸出的现象称为 泻流.泻流速率与分子的平均速率成正比:
泻流速率∝𝑛¯𝑣=𝑛√8𝑘𝑇𝜋𝑚
其中 𝑛
为分子数密度.
泻流的一个重要应用是 同位素分离.在相同温度下,轻分子的平均速率大于重分子,因此轻分子通过小孔的泻流速率更快.利用这一原理,可以通过反复泻流来富集特定同位素.
例题:铀同位素分离
题目:UF6
气体中 235UF6
和 238UF6
的泻流速率之比是多少?
解答:泻流速率 ∝¯𝑣 ∝1/√𝑀mol
.因此:
¯𝑣235¯𝑣238=√𝑀238𝑀235=√349352≈1.0043
单次泻流的分离效果很微弱(仅约 0.43%),因此需要通过数千级的级联过程才能实现有效的同位素富集.
§8 例题
例题 1:氧气分子在 300 K 下的特征速率
题目:计算氧气分子(𝑀mol =32 g/mol =0.032 kg/mol
)在 𝑇 =300 K
下的最概然速率、平均速率和方均根速率.已知 𝑅 =8.31 J/(mol·K)
.
解答:
𝑣𝑝=√2𝑅𝑇𝑀mol=√2×8.31×3000.032=√155813≈395m/s
¯𝑣=√8𝑅𝑇𝜋𝑀mol=√8×8.31×300𝜋×0.032=√198493≈446m/s
𝑣rms=√3𝑅𝑇𝑀mol=√3×8.31×3000.032=√233719≈484m/s
讨论:𝑣𝑝 <¯𝑣 <𝑣rms
(395 <446 <484 m/s
),与理论预期一致.这些速率都远小于地球逃逸速度(11.2 km/s
),因此氧气分子在常温下不会逃逸地球大气.
例题 2:速率在某区间内的分子比例
题目:在 𝑇 =300 K
下,氮气分子(𝑀mol =28 g/mol
)中,速率在 400
到 500 m/s
之间的分子约占多少?
解答:最概然速率为:
𝑣𝑝=√2𝑅𝑇𝑀mol=√2×8.31×3000.028≈422m/s
区间 [400,500] m/s
大致在 𝑣𝑝
附近.可以用以下近似方法估计:由于 Δ𝑣 =100 m/s
较小,且区间在峰值附近,𝐹(𝑣)
变化不太剧烈,因此:
∫500400𝐹(𝑣)d𝑣≈𝐹(𝑣𝑝)×Δ𝑣×(修正因子)
更精确的计算需要数值积分或查速率分布表.根据氧气在 273 K 的数值表类推,此区间的分子比例约为 15–20%.
常见误区
常见误区
1. (✗) "最概然速率是分子出现最多的速率"
这个说法不够准确.𝐹(𝑣𝑝)
是概率 密度 的最大值,不是概率本身.𝐹(𝑣) d𝑣
才是概率.在 𝑣𝑝
附近的一个微小区间 [𝑣𝑝,𝑣𝑝 +d𝑣]
内的分子比例最高,但 " 速率为 𝑣𝑝
的分子 " 这种说法在连续分布下没有意义(单点概率为零).
2. (✗) "温度升高意味着所有分子都运动得更快"
温度升高使分布变宽,高速分子的比例增大,但同时也有一些分子的速率很低.温度是统计平均量,不能用来描述单个分子.
3. (✗) "𝐹(𝑣)
就是概率 "
𝐹(𝑣)
是概率 密度,不是概率.概率是 𝐹(𝑣) d𝑣
(密度乘以区间宽度).𝐹(𝑣)
本身的值可以大于 1(当 𝑣
接近 𝑣𝑝
且温度较低时).
4. (✗) "麦克斯韦分布只适用于理想气体"
麦克斯韦分布描述的是速度(或速率)的统计分布,它依赖的假设是分子间相互作用可以忽略(使得速度分量统计独立).对于实际气体,在密度不太高时,麦克斯韦分布仍然是很好的近似.
学习衔接
本页面最近更新:2026/7/3 22:14:58,更新历史
发现错误?想一起完善? 在 GitHub 上编辑此页!
本页面贡献者:Leafuke, Physics Learning Wiki, 匿名同学
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用