6.2 热力学量的统计计算
6.2 热力学量的统计计算
本节摘要：上一节算出配分函数，本节把它"兑换成"热力学量。核心公式组：内能 U、自由能 A、熵 S、热容 C 都由配分函数 q 的各阶导数表示。正则配分函数连接 A = -kT ln q、S = -∂A/∂T、U = A + TS 等，从而把微观能级结构转化为宏观性质。玻尔兹曼解析为理想气体提供了从分子参数（转动常数、振动频率、质量）精确算摩尔热容的路径。本节用配分函数求数粒体系的能量、熵、热容，直观看到"微观→宏观"的翻译过程。
学习目标
阅读完本节，你应当能够：
- 用配分函数写出内能、熵、自由能的标准公式
- 说明自由度分解（平动/转动/振动/电子）如何叠加出分子配分函数
- 由平动配分函数推出理想气体的萨卡-泰特洛方程
- 计算简单体系的摩尔热容并解释经典均分极限
一、配分函数是"兑换表"
统计热力学答卷的固定动作是"用配分函数求宏观量"。最常用的一组关系：
A = -kT ln q S = k T (∂ln q/∂T)_V + k ln q U = k T² (∂ln q/∂T)_V C_V = (∂U/∂T)_V
从中可见：只要配分函数是温度的函数，四项热力学量就能全部解析出来。对多自由度的分子，配分函数按自由度分解相乘 q = q_trans·q_rot·q_vib·q_elec，因为能级可加→概率可乘→配分函数相乘。
二、平动配分函数与理想气体熵
理想气体的平动配分函数为 q_t = V(2πmkT/h²)^(3/2)。把它代入熵的公式并对 N 个粒子处理（考虑全同粒子不可分辨的统计因子 1/N!），就得到萨卡-泰特洛方程——可以根据分子量、体积、温度算出理想气体的绝对熵，无需任何实验拟合。能算出气体摩尔熵并与实验表的数值对得上，是统计方法最有力的自我证明。
三、转动与振动的贡献
- 转动配分函数 q_r = T/(Θ_rot·σ)，Θ_rot 是转动特征温度，σ 是对称数（对称分子如 CO₂ 的 σ>1 需除）。
- 振动配分函数 q_v = 1/(1 - e^(-Θ_vib/T))，Θ_vib 是振动特征温度。
在高温极限（T ≫ Θ），转动每自由度贡献 ½RT、振动每振子贡献 R 到内能——这就是经典的均分定理。低温下振动贡献快速“冻结”，热容低于均分预测，这解释了为什么异种气体的摩尔热容随温度先增后趋于定值。下面代码把配分函数的热容计算做出来，并注释经典极限。
import math import numpy as np kB = 1.380649e-23 h = 6.626e-34 NA = 6.022e23 def cv_vib(T, nu): """振动对摩尔热容的贡献: C(R) = (x² e^x/(e^x-1)²), x=Θ/T""" x = h*nu/(kB*T) return x*x*math.exp(x)/(math.exp(x)-1)**2 # 单位 R # 双原子分子振动特征温度 Θ ≈ hν/kB ; 取 ν~5e13 Hz (中等键) nu = 5.0e13 for T in [100, 300, 1000, 3000]: cv = cv_vib(T, nu) print(f"T={T:5d} K 振动贡献热容 ≈ {cv:5.2f} R") print("对比: 经典均分每振动自由度贡献 1 R；低温远低于此值")
低温（100 K）时振动热容远小于 1 R——振动自由度"冻结"；到 3000 K 附近逼近 1 R——均分极限恢复。这组数据解释了一个常见现象：分子摩尔热容并非恒定，而是随温度从“冻结到解冻”变化。
值得一提的是，化学平衡常数也能直接从配分函数之比算得。对气相反应 aA + bB ⇌ cC + dD，平衡常数可写成 K = Π(q_i/N_A)^ν · e^(-Δε₀/kT) 的乘积形式，其中 ν 为计量系数（产物为正、反应物为负）。如此就把反应的焓变与熵变统一进一个由配分函数构成的表达式里——我们无需翻阅平衡数据表，仅凭各分子的转动常数、振动频率与质量，就能从第一性原理预判反应偏向产物还是偏向反应物。这正是统计热力学把平衡常数从查表经验提升为逐分子计算的意义所在。
四、真实案例：解释氢气的"低温比热"
实验发现氢气在很低温度下摩尔等容热容只有平动+转动部分，振动完全"冻住"；升到高温后逐步加上振动贡献才接近经典值。统计计算把这条热容-温度曲线精确预测出来，无需经验参数，只需分子的转动/振动特征温度。航天推进、流体力学的燃气热力学计算，都靠这种"分子参数→热容→热力学"的链条，在数值代码里预置。
6.2 图：振动热容随温度的“冻结→均分”演化
五、与其他章节的关联
本节是 6.1 配分函数的具体输出，为 6.3 节理想气体/晶体的精确统计处理铺路。得到的摩尔热容等量反过来可查第 2 章热力学表，形成"统计算出→热力学表核对"的闭环。化学平衡常数也能从配分函数之比推导，是"计算化学"量化预测的基础。
本节要点回顾
- 兑换公式：A=-kT ln q、U=kT²(∂ln q/∂T)、S=kT(∂ln q/∂T)+k ln q
- 自由度分解：q=q_trans·q_rot·q_vib·q_elec，能级可加→概率可乘→q 相乘
- 萨卡-泰特洛方程：由平动配分函数算出理想气体绝对熵，无实验参数
- 均分与冻结：高温振动每振子约 R，低温热容骤降、低频自由度冻结
- 计算用途：分子参数→热容→热力学，是燃气/推进热力学计算的标准链条