Learning objectives / 学习目标
- 用压缩因子 \(Z\) 定量判断实际气体偏离理想行为的方向与大小。
- 解释吸引作用常使 \(Z<1\)、短程排斥常使 \(Z>1\)，同时避免把这一经验趋势当成无条件定律。
- 使用 van der Waals 方程并从临界点拐点条件推导 \(V_{m,c},T_c,p_c,Z_c\)。
- 严格区分以摩尔密度或以压强展开的 virial 系数、单位和截断条件。
- 分析 Lennard–Jones 势的零点、极小值与力方向，连接束缚、凝结和临界现象。
- 通过量纲、极限、级数项大小和热力学稳定性审查模型。
Prerequisites / 先修知识
- Ideal-Gas Model 与 \(pV_m=RT\)；
- 导数、二阶导数、级数近似和无量纲化；
- 保守力 \(F_r=-dU/dr\)、势阱与稳定平衡。
Scope and conventions / 范围与约定
研究纯物质、单相、宏观平衡状态，使用 SI 和摩尔体积 \(V_m=V/n\)。压强均为绝对压强。除特别说明，virial 展开采用摩尔密度形式 \(Z=1+B_2/V_m+B_3/V_m^2+\cdots\)；若改用压强展开会另换符号。van der Waals 方程用于教学和定性结构，不替代高精度多参数状态方程。亚稳态、成核动力学、混合物相平衡与量子流体只作边界提示。
Notation and modeling contract / 符号与建模契约
| Symbol | Meaning | SI unit / convention |
|---|---|---|
| \(V_m\) | molar volume | \(\mathrm{m^3\,mol^{-1}}\) |
| \(Z=pV_m/(RT)\) | compressibility factor | dimensionless |
| \(a,b\) | van der Waals attraction and covolume parameters | \(\mathrm{Pa\,m^6\,mol^{-2}}\), \(\mathrm{m^3\,mol^{-1}}\) |
| \(B_2,B_3\) | density-form virial coefficients | \(\mathrm{m^3\,mol^{-1}}\), \(\mathrm{m^6\,mol^{-2}}\) |
| \(U(r),F_r\) | pair potential and radial force | J, N |
| \(\varepsilon,\sigma\) | Lennard–Jones well depth and zero-crossing length | J, m |
\(\sigma\) 在本课 K04 表示 Lennard–Jones 长度，不是 Stefan–Boltzmann 常数；\(b\) 是每摩尔的共体积参数，不等于单分子的几何体积。
Knowledge 1 — Deviations from ideal-gas behavior / 对理想气体行为的偏离
Knowledge ID: P03-U02-L03-K01
理想气体忽略分子体积和长程相互作用。真实分子在中等距离相吸、极短距离强烈排斥；因此高温、低密度极限通常最接近理想气体，而低温或高压时偏差增大。
定义压缩因子（compressibility factor）
\[Z=\frac{pV_m}{RT}.\]
\(Z=1\) 为理想值；在许多气体的中低压温区，吸引使给定 \(T,V_m\) 下的壁面压强降低，常见 \(Z<1\)；更高密度下排斥和排除体积占优，常见 \(Z>1\)。
该图只表达常见趋势：某些温度下曲线可从一开始就在 \(Z>1\)，因此必须依据数据或可靠状态方程，而不是凭“真实气体一定更易压缩”下结论。
在 \(p\to0\)（等价于 \(V_m\to\infty\)）的经典气体极限，分子相距很远，\(Z\to1\) 是重要一致性检查。
Knowledge 2 — van der Waals equation / 范德瓦耳斯方程
Knowledge ID: P03-U02-L03-K02
van der Waals 方程用两个参数对理想式作最简单修正：
\[\boxed{\left(p+\frac{a}{V_m^2}\right)(V_m-b)=RT,}\]
或
\[p=\frac{RT}{V_m-b}-\frac{a}{V_m^2}.\]
\(b>0\) 通过 \(V_m-b\) 描述不可用的排除体积；\(a>0\) 的项降低压强，代表平均场吸引修正。
方程自动恢复理想极限：当 \(V_m\gg b\) 且 \(a/V_m^2\ll p\) 时，\(pV_m\approx RT\)。但 \(V_m\le b\) 不在模型允许域内。为比较临界点两侧而不依赖任意的 \(a,b\) 数值，下面采用约化变量 \(p_r=p/p_c\)、\(V_r=V_m/V_{m,c}\)、\(T_r=T/T_c\)：
\[p_r=\frac{8T_r}{3V_r-1}-\frac{3}{V_r^2}, \qquad V_r>\frac13.\]
低于临界温度的 S 形“回环”含 \(\partial p/\partial V_m>0\) 的机械不稳定段，不能直接当作可观测均匀平衡态；真实系统在相共存区由液、气两相取代。Maxwell 等面积构造是该模型的一种热力学修正，本课只说明其角色。
临界点满足等温线水平拐点条件：
\[\left(\frac{\partial p}{\partial V_m}\right)_{T_c}=0, \qquad \left(\frac{\partial^2p}{\partial V_m^2}\right)_{T_c}=0.\]
对 van der Waals 方程求解得
\[V_{m,c}=3b, \qquad T_c=\frac{8a}{27Rb}, \qquad p_c=\frac{a}{27b^2}, \qquad Z_c=\frac38.\]
Knowledge 3 — Virial expansion and compressibility factor / 维里展开与压缩因子
Knowledge ID: P03-U02-L03-K03
virial 展开把低密度偏差按密度幂次系统排列。采用摩尔密度形式：
\[Z=\frac{pV_m}{RT} =1+\frac{B_2(T)}{V_m}+\frac{B_3(T)}{V_m^2}+\cdots.\]
\(B_2\) 单位为 \(\mathrm{m^3\,mol^{-1}}\)，\(B_3\) 为 \(\mathrm{m^6\,mol^{-2}}\)；每一项必须无量纲。\(B_2\) 首先编码两体相互作用，\(B_3\) 包含三体关联等更高阶效应。
另一常见写法以压强展开：
\[Z=1+B'_2(T)p+B'_3(T)p^2+\cdots.\]
带撇系数的单位与数值都不同，不能与密度形式的 \(B_2,B_3\) 混用。任何截断都必须检查后一项相对前一项是否足够小。
某温度若 \(B_2(T_B)=0\)，称第二 virial 意义下的 Boyle temperature；此时最低阶修正消失，但 \(B_3/V_m^2\) 等仍可非零，所以并不表示气体在所有压力下严格理想。
Knowledge 4 — Intermolecular potentials and condensation / 分子间势与凝结
Knowledge ID: P03-U02-L03-K04
一个常用的中性分子对势模型是 Lennard–Jones 12–6 势：
\[U(r)=4\varepsilon\left[\left(\frac{\sigma}{r}\right)^{12} -\left(\frac{\sigma}{r}\right)^6\right].\]
它在 \(r=\sigma\) 过零，在
\[r_{min}=2^{1/6}\sigma\]
取得极小值 \(U_{min}=-\varepsilon\)。径向力为
\[F_r=-\frac{dU}{dr} =\frac{24\varepsilon}{r}\left[2\left(\frac{\sigma}{r}\right)^{12} -\left(\frac{\sigma}{r}\right)^6\right].\]
当热扰动尺度 \(k_BT\) 相对势阱深度 \(\varepsilon\) 较小、密度又足够高时，吸引可使分子形成凝聚相；但“\(k_BT\approx\varepsilon\)”只给数量级直觉，临界温度和共存曲线是集体现象，不能由一对分子势的最低点直接等同得到。
临界点以上，液体与气体之间不再有一阶相界面，可连续改变密度；这不等于分子间作用突然消失，而是宏观相区分消失。
Worked Examples / 完整例题
Worked Example 1 — 用 \(Z\) 判断偏离方向
已知： \(p=5.00\ \mathrm{MPa}\)，\(T=300\ \mathrm K\)，测得 \(V_m=4.00\times10^{-4}\ \mathrm{m^3\,mol^{-1}}\)。
系统与参考系： 均匀单相气体，实验室参考系；状态量均为绝对量。
模型： 不预设状态方程，只用 \(Z\) 定义比较理想值。
计算：
\[Z=\frac{pV_m}{RT} =\frac{(5.00\times10^6)(4.00\times10^{-4})}{(8.314462618)(300)} =0.802.\]
检查与解释： \(Z\) 无量纲且低于 1，表示该状态的 \(pV_m\) 比理想值小 \(19.8\%\)；可说吸引效应在此状态显著，但不能仅凭单一点反推出唯一分子势。
Worked Example 2 — van der Waals 压强与理想值比较
已知： 取氮气教学参数 \(a=0.137\ \mathrm{Pa\,m^6\,mol^{-2}}\)，\(b=3.87\times10^{-5}\ \mathrm{m^3\,mol^{-1}}\)，\(T=300\ \mathrm K\)，\(V_m=1.00\times10^{-3}\ \mathrm{m^3\,mol^{-1}}\)。
系统与模型： 单相均匀气体，van der Waals 平均场；\(V_m>b\)。
计算：
\[p=\frac{RT}{V_m-b}-\frac{a}{V_m^2} =2.59\times10^6-1.37\times10^5 =2.46\times10^6\ \mathrm{Pa}.\]
理想值为 \(p_{\rm id}=RT/V_m=2.49\ \mathrm{MPa}\)，故该模型给 \(Z=p/p_{\rm id}=0.985\)。
检查： 两项单位均为 Pa；吸引项降低压强，而排除体积提高第一项。最终谁占优必须由数值决定。
Worked Example 3 — 从 \(a,b\) 推导临界常数
已知： 使用上例 \(a,b\)。
研究对象与模型： van der Waals 纯气体临界点，满足一、二阶体积导数为零。
计算：
\[V_{m,c}=3b=1.16\times10^{-4}\ \mathrm{m^3\,mol^{-1}},\]
\[T_c=\frac{8a}{27Rb}=126\ \mathrm K, \qquad p_c=\frac{a}{27b^2}=3.39\ \mathrm{MPa}.\]
进一步
\[Z_c=\frac{p_cV_{m,c}}{RT_c}=\frac38=0.375.\]
检查与解释： \(a/(Rb)\) 单位为 K，\(a/b^2\) 单位为 Pa。\(Z_c=3/8\) 是该两参数模型的普适预测，不是所有真实物质的精确常数。
Worked Example 4 — virial 截断与项大小
已知： \(T=300\ \mathrm K\)，\(V_m=1.00\times10^{-3}\ \mathrm{m^3\,mol^{-1}}\)，\(B_2=-1.20\times10^{-4}\ \mathrm{m^3\,mol^{-1}}\)，\(B_3=1.00\times10^{-8}\ \mathrm{m^6\,mol^{-2}}\)。
系统与模型： 单相低密度气体；截断到第三 virial 项。
计算：
\[Z=1+\frac{B_2}{V_m}+\frac{B_3}{V_m^2} =1-0.120+0.0100=0.890.\]
\[p=\frac{ZRT}{V_m}=2.22\times10^6\ \mathrm{Pa}.\]
截断检查： 第三项大小 \(0.0100\) 是第二项 \(0.120\) 的 \(8.3\%\)，并非完全可忽略；若未知第四项，就不能宣称高精度。所有修正项无量纲，符号显示两体吸引占主导。
Lesson synthesis / 全课整合
三层描述彼此互补：\(Z\) 是实验偏差的无量纲摘要；van der Waals 方程用 \(a,b\) 建立最小平均场模型并产生临界结构；virial 展开按低密度幂次系统组织实验或统计力学信息；分子对势解释这些系数为何随温度变化以及凝结为何可能发生。
可执行流程：先算 \(Z\) 判断偏差；若需快速定性用 van der Waals，但检查 \(V_m>b\) 与稳定性；若有 virial 数据，声明展开变量并比较逐项大小；若问微观原因，分析 \(U(r)\)、\(F_r\) 和 \(k_BT/\varepsilon\)，但不把两体直觉冒充完整相平衡理论。
Formula-and-decision table / 公式—决策表
| Goal | Relation | Use when | Main caution |
|---|---|---|---|
| quantify nonideality | \(Z=pV_m/(RT)\) | measured single-phase state | \(Z<1\) is a state-specific observation |
| two-parameter model | \(p=RT/(V_m-b)-a/V_m^2\) | qualitative real-gas modeling | not precision EOS; require \(V_m>b\) |
| vdW critical point | \(V_{m,c}=3b\), \(T_c=8a/(27Rb)\), \(p_c=a/(27b^2)\) | same vdW model | \(Z_c=3/8\) is model-specific |
| density virial series | \(Z=1+B_2/V_m+B_3/V_m^2+\cdots\) | low density with coefficients | check units and next-term size |
| pair interaction | \(U=4\varepsilon[(\sigma/r)^{12}-(\sigma/r)^6]\) | simple neutral-molecule model | exponents and parameters are empirical model choices |
| force direction | \(F_r=-dU/dr\) | central conservative pair potential | sign is radial, not scalar attraction label |
Model boundaries / 模型边界
- van der Waals 方程的回环含不稳定均匀态；共存区需相平衡条件或 Maxwell 构造。
- virial 级数是低密度展开，临界点附近通常收敛变差，有限项不可靠。
- \(Z<1\) 或 \(Z>1\) 只描述给定状态；不可把某物质永久分类为“吸引型”或“排斥型”。
- Lennard–Jones 12–6 势适合某些近球形中性粒子的定性或半定量模型；极性、氢键、多体效应需更丰富势能。
- 凝结与临界性是多体、集体和热力学极限现象；两分子形成束缚对并不等价于出现液相。
- 极高温下化学反应、电离改变组成；极低温或高密度下量子统计可能重要。
Misconception clinic / 易错概念诊所
| Misconception | Why it fails | Repair |
|---|---|---|
| 实际气体总有 \(Z<1\) | 高密度排斥可使 \(Z>1\) | 指定 \(p,T\) 后再判断 |
| 参数 \(b\) 就是分子总体积 | 它是排除体积的摩尔参数 | 保留模型定义与单位 |
| vdW 回环每一点都可稳定观测 | 正斜率段机械不稳定 | 用相共存替代 |
| Boyle temperature 使所有压力都理想 | 只消去第二 virial 项 | 高阶项仍存在 |
| 两种 virial 写法的 \(B_2\) 相同 | 展开变量不同导致单位与数值不同 | 明确密度式或压强式 |
| \(r=\sigma\) 是势能最低点 | 它只是过零点 | 最低点在 \(2^{1/6}\sigma\) |
| 临界点处分子间力消失 | 消失的是液气宏观区分 | 相互作用仍存在 |
Practice and retrieval / 练习与检索
- 写出 \(Z\) 的定义，说明 \(Z=0.92\) 与 \(Z=1.08\) 分别代表什么。
- 为什么 \(p\to0\) 时经典实际气体应趋向 \(Z=1\)？
- 在 van der Waals 方程中分别指出 \(a\) 与 \(b\) 修正的物理方向。
- 推导临界条件需要哪两个偏导数为零？为什么只令一阶导数为零不够？
- 若 \(a=0.200\ \mathrm{Pa\,m^6\,mol^{-2}}\)、\(b=5.00\times10^{-5}\ \mathrm{m^3\,mol^{-1}}\)，求 \(V_{m,c}\)。
- 密度式 virial 展开中，\(B_4\) 的 SI 单位是什么？
- \(B_2=-50\ \mathrm{cm^3\,mol^{-1}}\)、\(V_m=10.0\ \mathrm{L\,mol^{-1}}\)，只保留第二项，求 \(Z\)。
- 对 Lennard–Jones 势，写出过零点、极小点位置与极小值。
- 两分子间距小于 \(r_{\min}\) 时力方向怎样？用 \(F_r=-dU/dr\) 解释。
- 为什么不能用 \(k_BT=\varepsilon\) 直接宣称这就是精确临界温度？
Answer 1
\(Z=pV_m/(RT)\)。\(0.92\) 表示 \(pV_m\) 比理想值低 \(8\%\)；\(1.08\) 表示高 \(8\%\)，但原因仍需结合状态和模型判断。
Answer 2
密度趋零时平均粒子间距增大，有限程相互作用和排除体积对每个粒子的相对修正趋零，因此恢复理想极限。
Answer 3
\(b\) 使可用体积减小，从而提高 \(RT/(V_m-b)\)；\(a\) 项以 \(-a/V_m^2\) 降低压强。
Answer 4
需同时满足 \((\partial p/\partial V_m)_T=0\) 与 \((\partial^2p/\partial V_m^2)_T=0\)；第一条件只找水平切线，第二条件才选出临界拐点。
Answer 5
\(V_{m,c}=3b=1.50\times10^{-4}\ \mathrm{m^3\,mol^{-1}}\)。
Answer 6
\(B_4/V_m^3\) 必须无量纲，所以 \(B_4\) 单位为 \(\mathrm{m^9\,mol^{-3}}\)。
Answer 7
统一单位后 \(B_2/V_m=(-0.0500\ \mathrm{L\,mol^{-1}})/(10.0\ \mathrm{L\,mol^{-1}})=-0.00500\)，故 \(Z=0.995\)。
Answer 8
\(r=\sigma\) 时 \(U=0\)；\(r=2^{1/6}\sigma\) 时 \(U=-\varepsilon\)，这是极小点。
Answer 9
在 \(r<r_{\min}\) 的短程侧，\(F_r>0\)（径向向外），表现为排斥；负梯度把粒子推向更大间距。
Answer 10
\(k_BT/\varepsilon\) 只是两种能量尺度之比；临界性依密度、势的完整形状和多体关联，比例常数也由具体体系决定。
Exit ticket / 离堂检验
给你一组 \(p,V_m,T\) 数据和某个 \(B_2(T)\)。设计一条不超过六步的路线，判断是否可用二阶 virial、是否需要第三项，以及 van der Waals 或 Lennard–Jones 图景能解释什么、不能解释什么。
References and source notes / 参考资料与使用说明
Source-to-content map
| Source | Used for | Knowledge ID / Section | Model boundary or convention |
|---|---|---|---|
| OpenStax University Physics Vol. 2, §2.1 | ideal-model limits and van der Waals introduction | K01–K02 | qualitative undergraduate equation of state |
| MIT OCW 5.60, Lecture 1 notes | state variables, real-gas corrections, critical-point thermodynamic context | K01–K02 | equilibrium pure fluids |
| MIT OCW 5.62 — Virial Equation of State | density virial expansion and molecular interpretation | K03–K04 | low-density statistical mechanics |
| NIST — Virial Coefficients of Pure Gases | experimental virial-coefficient scope and data provenance | K03 | pure gases and temperature-dependent coefficients |
| NIST TDE — Gas-density equation descriptions | virial equation forms and coefficient conventions | K03 | equation-form convention |
Source notes
以上来源均于 2026-08-24 实际打开核对。MIT 5.62 与 NIST 用于防止混淆 density virial 与 pressure virial 约定；van der Waals 与 Lennard–Jones 图均按方程原创绘制。所有 \(a,b,B_i\) 示例值明确标为教学模型值，不冒充高精度 NIST 状态数据。
复习与拓展
- Unit 3.3 将把状态方程用于热力学路径、功和第一定律；
- 相平衡课程将引入化学势、Maxwell 构造和 Clapeyron 方程；
- 统计力学将从 Mayer 函数推导 \(B_2(T)\)；
- 数值物理可比较 van der Waals、virial 与高精度多参数 EOS 的残差。