LAMMPS计算自由能：固液界面TI-US双路径的λ策略与收敛判据

一个Cu(111)-水界面的自由能计算项目，第一次用热力学积分（TI）跑完11个等间距λ窗口（0至1，步长0.1），积分出来的固液界面自由能值是0.42 J/m²——和文献中实验值的0.29 J/m²差了40%以上。翻看每个λ窗口的∂H/∂λ曲线后找到了问题：λ接近0和1的端点窗口，∂H/∂λ的数值比相邻窗口陡峭三倍以上，11点梯形积分在两端各漏掉了约0.06 J/m²的积分面积。

问题在于等间距λ方案默认了∂H/∂λ在λ∈[0,1]上是平滑函数——但在固液相变这种一级相变附近的体系里，端点处的哈密顿量变化速率极端非线性。换用Gaussian-Quadrature型非均匀λ分布（两端密、中间疏）之后，同样的采样量下积分精度从0.13 J/m²的误差压到了0.015 J/m²——剩余的偏差已经在热噪声量级以下。这套λ方案在项目后续的Al(111)-水界面验证中也保持了一致的收敛表现。

每个λ窗口内的采样时间也需要按体系统一化。∂H/∂λ的自相关时间对含4000个原子左右的界面体系约为50-80 ps，每步窗口需要先跑0.5 ns平衡以衰减初始构型记忆，再跑至少1.5 ns的生产采样才能使running average的波动范围收窄到标准误差的三倍以内。相邻λ窗口之间的∂H/∂λ分布重叠度是另一项收敛判据——重叠面积低于30%时，无论running average看起来多平稳，中间都必须追加一个λ窗口。这个项目中在λ=0.85至0.95之间补插了λ=0.90窗口之后，积分数值的变化才降到0.01 J/m²以下。

伞形采样（US）在LAMMPS中通过fix colvars调用colvars模块实现。反应坐标选为表面法线方向的距离z，沿z方向布置24个窗口，每个窗口加谐振偏置势把水分子质心限制在预设的z值附近。力常数和窗口间距之间存在一个平衡：力常数k太小则窗口偏置兜不住水分子的热涨落导致重叠不足，k太大则窗口太窄需要更多窗口数。这个项目以k≈k_BT/(0.1 Å)²（约1.0 kcal/mol/Ų）为起点，窗口间距取谐振宽度σ≈sqrt(k_BT/k)的1.5倍（约0.15 Å），24个窗口覆盖了界面法线方向3.6 Å的范围。

WHAM重建PMF的收敛性通过”增样验证”来检验：连续两次新增50%采样帧数后PMF变化量小于1 kcal/mol。这个项目在采样量×1.5之后，PMF曲线在±0.5 kcal/mol范围内重合。WHAM本身的迭代收敛用10000次迭代做终检——1000次和10000次迭代的PMF结果差异小于0.02 kcal/mol，确认了迭代循环已经退出。

TI和US两条路径算出来的界面自由能在0.29至0.31 J/m²之间一致对上了——和实验值的偏差在统计误差范围内。之后这个项目又做了一条验证线：用FEP在相同体系统做自由能差分验证，结果在TI-US的置信区间内，三项交叉确认给了审稿人足够的信心。