热力学第一性原理：从 Frenkel & Smit 读懂分子模拟的基础

这是读 Frenkel & Smit《Understanding Molecular Simulation》第二章的笔记。这一章讲经典热力学——分子模拟的地基。

热和功不是"东西"

$$dE = \delta q + \delta w$$

Frenkel & Smit 用 δ 而不是 d 来写 q 和 w。这不是印刷偏好——热和功不是状态函数。

你不能说一个系统"含有多少功"。但你可以说它"做了多少功"。E 是状态函数（有 dE），q 和 w 是过程量（只有 δq, δw）。

对于最常见的 PV 系统：

$$dE = TdS - PdV + \mu dN$$

三种方式改变系统能量：

加热（TdS）
压缩（-PdV）
加粒子（μdN）

就这么三种。

熵不是神秘的"方向性"

$$dS = \frac{\delta q_{rev}}{T}$$

孤立系统自发趋向熵最大。这不是什么神秘的宇宙方向性——统计力学（§2.2）会告诉你，熵的本质是微观状态数的对数 Ω。

容易混淆的点：第二定律说的是孤立系统的熵最大。

如果系统不是孤立的（比如跟热浴接触），那就不是熵最大，而是其他热力学势最小。这就是 Legendre 变换出现的动机。

系统自己"知道"温度

给定 N, V, E 三个独立变量，S 就完全确定了：

$$S = S(N, V, E)$$

从 S 的微分可以直接读出 T, P, μ：

$$\frac{1}{T} = \left(\frac{\partial S}{\partial E}\right)_{N,V}$$

这是理解分子动力学的关键：在 NVE 系综（标准 MD）中，系统自己"知道"温度——它是 S 对 E 的导数，不是外部设定的。

你想要恒温？就要引入恒温器（thermostat），这发生在 §2.3。

广延量 vs 强度量

Extensive（广延量）：E, S, V, N — 随系统大小线性增长
Intensive（强度量）：T, P, μ — 跟系统大小无关

Gibbs-Duhem 关系告诉你一件事：T, P, μ 不是三个独立变量。

$$SdT - VdP + Nd\mu = 0$$

给定任意两个，第三个就被确定了。

水和蒸汽在相平衡时：T 和 P 固定 → μ 也被固定了（两相相等）。

对模拟的影响：μVT（巨正则系综）中，你固定 μ, V, T——但因为 Gibbs-Duhem，这等价于固定了 P。你不能同时固定 μ, V, T 和 P。

Legendre 变换：换变量，换势函数

很多时候我们在实验中控制的不是 N, V, E，而是 N, P, T——比如烧杯里的反应，大气压、室温、固定量。

Legendre 变换的做法：

原本：E = E(S, V, N)
我想把自变量从 S（不好控制）换成 T（好控制）

F = E - TS
→ F = F(T, V, N) ← 现在 T 是自变量了！

为什么减 TS？因为 ∂E/∂S = T，减掉 TS 就"抵消"了 S 的作用。

所以 dF = dE - TdS - SdT = -SdT - PdV + μdN

S 的微分项消失了。

四种系综对应四个势函数

系综	势函数	约束条件
NVE	E (内能)	N, V, E 固定
NVT	F (Helmholtz 自由能)	T, V, N 固定
NPT	G (Gibbs 自由能)	T, P, N 固定
μVT	Ω (巨势)	T, V, μ 固定

每种系综都对应一个 Legendre 变换后的势函数，平衡条件就是这个势取极小值（熵取极大值）。

两相平衡：T, P, μ 都相等

从两个子系统的推导可以得出一个统一原则：

在给定约束下，对应的热力学势取极小值（熵取极大值）

两相平衡时：

$$T_1 = T_2, \quad P_1 = P_2, \quad \mu_1 = \mu_2$$

温度相等、压力相等、化学势相等。三个条件缺一不可。

化学势 μ

$$\mu = \left(\frac{\partial E}{\partial N}\right)_{S,V}$$

加一个粒子的能量代价。

为什么模拟中 μ 很重要？

GCMC（巨正则 Monte Carlo）需要 μ 作为输入参数
相平衡（气-液、液-液）要求 μ 相等
酸在光刻胶中的扩散本质上是化学势梯度驱动的

最后一个跟我的工作直接相关。光刻胶里的酸扩散不是简单的浓度梯度驱动，而是化学势梯度。这意味着即使浓度均匀，如果局部环境（极性、自由体积）不同，酸仍然会迁移。

常见误区

"熵最大" vs "自由能最小" — 不是矛盾，是不同的约束条件
化学势的符号 — dE = +μdN 但 dS = -(μ/T)dN，注意正负号
Legendre 变换不是物理的 — 它只是数学换变量，系统的物理本质没变

下一步

§2.2 统计力学——熵的微观解释。§2.3 系综和恒温器。