第 3 章 · 正则系综与配分函数：热力学的万能生成器

配分函数 $Z$ 是统计物理的"总开关"：算出一个积分（或求和），全部热力学量自动落袋。这一章把"从 $Z$ 到一切"的流水线标准化。

配分函数与自由能

自由能 $F = E - TS$ 为什么是正确的"生成元"：它在温度体积固定下取极小（平衡判据），且它的微分 $dF = -S\,dT - P\,dV$ 让熵与压强由求导获得：

能量 $E = F + TS = -T^2\partial(F/T)/\partial T$。一函定乾坤：算 $Z$ 的功换来全套热力学。

经典极限与理想气体

能量连续时 $Z = \int e^{-E(p,q)/kT}\frac{d^3p\,d^3q}{(2\pi\hbar)^{3N}N!}$。理想气体（$E = \sum p^2/2m$）：

状态方程从积分里掉出来；$\lambda = h/\sqrt{2\pi mkT}$ 是热波长——$N!$ 与 $\lambda^{3N}$ 是量子力学留给经典统计的两枚印章（吉布斯佯谬的解决）。

能量均分与比热

单粒子每个二次自由度贡献 $\frac{1}{2}kT$：单原子气体 $C_V = \frac{3}{2}Nk$、双原子再加转动与振动（常温只有转动激活、振动冻结——能级间隔 $> kT$ 时自由度"冻住"，经典困境的量子解答）。固体的杜隆-珀蒂值 $3Nk$ 与低温 $T^3$ 律（声子，第 4-5 章回收）在此埋线。

老师的提醒

• 算配分函数时的三件套检查：正确量子计数（$N!$、$(2\pi\hbar)^{-3N}$）、正确的能量函数、正确的边界条件（有限体积归一化）。
• 均分的"激活"判据 $kT$ vs 能级间隔是比热问题的万能钥匙：转动 $\theta_r$ 小、振动 $\theta_v$ 大——数据表自己会说话。

自检

• 由 $F = -kT\ln Z$ 推出 $S$ 与 $P$。
• $N!$ 因子解决什么佯谬？
• 双原子气体比热随温度的台阶形状？

更深一层：配分函数是"热力学的一键生成器"

$Z = \sum_n e^{-E_n/kT}$ 被称为"配分函数"（partition function），它的功能配得上这个程序化的名字：算出 $Z$ 之后，所有热力学量由求导一键生成——自由能 $F = -kT\ln Z$、熵 $S = -\partial F/\partial T$、压强 $P = -\partial F/\partial V$、内能 $E = -T^2\partial(F/T)/\partial T$、热容 $C_V = \partial E/\partial T$。"算 Z 的功"换"全套热力学"的账——这条流水线是统计物理的核心生产力。

配分函数的两个深层性质：广延性（$Z_{AB} = Z_A \cdot Z_B$，独立子系统相乘 ⟹ $F$ 可加）与单调性（$Z$ 随温度单调增）。理想气体的 $Z$ 由单粒子配分函数组装（$Z_N = Z_1^N/N!$——$N!$ 是全同粒子的量子计数修正），吉布斯佯谬的解决就藏在这个 $N!$ 里（不可分辨的粒子重复计数必须除掉）。

推导细读：能量均分的"激活阈值"读法

每个二次自由度贡献 $\frac12 kT$——这条均分定理的"激活阈值"解释是理解比热的关键：转动自由度在 $kT > \hbar^2/2I$（特征温度量级）时激活、振动自由度在 $kT > \hbar\omega$ 时激活。双原子气体比热的台阶（常温 $C_V = \frac52Nk$ 只含平动+转动、高温 $\frac72Nk$ 加振动）正是能级间隔穿过 $kT$ 时自由度"解冻"的记录。经典比热困境（均分预言错）与量子修正（能级间隔）的和解，全在 $kT$ 与能级间隔的对比里。

在真实世界里：配分函数的"工程接口"

配分函数是统计物理进入工程的接口：化工（反应平衡常数 $K = e^{-\Delta F/kT}$——配分函数比直接给出平衡移动方向）、材料（相图计算 = 两相 $F$ 的极小化竞争）、生物物理（蛋白质折叠的 $F$ 地形、配体结合的 $\Delta F$）、机器学习（玻尔兹曼分布与 softmax 同源——$p_i = e^{-E_i/kT}$ 就是 softmax）。看到 softmax 就想到配分函数，看到配分函数就想到自由能——两门学科共用一台发动机。

专题深入：正则系综的"热浴契约"完整版

正则系综的物理契约要讲透：系统与巨大热浴接触（温度固定、能量涨落），分布 $p_i \propto e^{-\beta E_i}$ 与系统自身微观细节无关——这就是为什么配分函数能"代换系统"：换一个物理系统只需换 $E_i$ 清单，统计机器不动。能量相对涨落 $\Delta E/E \sim 1/\sqrt{N}$ 在宏观极限消失——热力学量在 $N \to \infty$ 时锐化，涨落是 $1/\sqrt N$ 的残差。这条 $N^{-1/2}$ 律在相变处失效：临界点关联长度发散、涨落变得宏观（临界乳光——CO₂ 在临界点整管变浑浊），这是平均场理论失效、标度理论登场的信号。朗道自由能 $F = F_0 + a(T - T_c)\eta^2 + b\eta^4$ 的展开把对称性与序参量放进同一函数：对称性 ($\eta \to -\eta$) 限制展开项、系数变号驱动相变——"对称性 + 解析性"两假设的产出是全部平均场指数（$\beta = 1/2$ 等）。

推导细读：配分函数的"解剖术"完整版

配分函数计算的四板斧值得完整演练：两能级系统（$Z = 2\cosh(\beta\epsilon)$，热容在 $kT \sim \epsilon$ 处出现 Schottky 峰——顺磁盐、二能级缺陷、电子自旋的热容贡献全用它拟合）；谐振子（$Z = 1/(2\sinh(\beta\hbar\omega/2))$，高温回归经典均分 $\langle E\rangle = kT$、低温冻结在零点能——固体热容的爱因斯坦模型原型，解释低温热容指数下跌）；刚性转子（$Z_{\rm rot} = T/\sigma\theta_r$ 高温版，转动能级间隔决定同核双原子气体的热容"解冻温度"——氢气在 85 K 以下转动冻结、热容减半的实测案例）；理想气体平动（准连续能级给出 Sackur-Tetrode 熵公式，熵常数第一次被纯理论算出，吉布斯佯谬的分辨率——全同粒子可分辨性差一个 $N!$，熵是广延量的必要修正）。四板斧合起来："热容-温度曲线"是分子能谱的放大镜。

工程案例细节：热力学数据库与材料设计的工程版

第二章方法的工程出口：化学平衡常数的统计计算（$\Delta G^\circ = -RT\ln K$，$K$ 由配分函数比给出——化工流程模拟软件（Aspen Plus）的反应数据库底层就是光谱数据 + 统计力学；燃烧、合成氨、CVD 外延的平衡组成全部由此计算）；气相色谱与蒸馏设计（分配系数的温度依赖由配分函数决定——塔板数与回流比设计的物理输入）；热电材料优值（$ZT = S^2\sigma T/\kappa$，Seebeck 系数的玻尔兹曼输运计算——温差发电与固态制冷的材料筛选）；相图计算 CALPHAD（各相自由能的参数化模型拟合实验数据，多元合金相图外推——高温合金、锂电池电解质配方设计的工业标准工具）。"从谱到 Z、从 Z 到 F、从 F 到相图"的链条是计算材料学的主动脉。