lammps计算自由能

我第一次做 lammps计算自由能是把 A 原子慢慢”变”成 B 原子算固溶形成自由能，用 TI 积分 ⟨∂U/∂λ⟩，但忘了 LAMMPS 的 pair_style 在 λ 端点对长程库仑的处理，积分曲线在 λ→1 处翘尾巴，结果偏差 0.1 eV/atom。那次让我明白，lammps计算自由能的坑不在公式，在势函数怎么随 λ 平滑过渡。

一、为什么 lammps计算自由能先定方法

LAMMPS 里自由能主要三招：TI（对耦合参数求导积分）、FEP（指数平均，配合 fix adapt 调相互作用）、umbrella sampling + WHAM（沿反应坐标采）。材料/合金相变我默认 TI，吸附/相界面我用 umbrella。曾经我硬用 FEP 算一个二级相变，采样跨不过自由能垒，结果完全不可信。做 lammps计算自由能时，先判断”有没有势垒”比选软件更该优先——有垒必走伞采样。

二、核心原理：TI 的积分路径

TI 把两态用 λ 连接，ΔG = ∫₀¹ ⟨∂U/∂λ⟩_λ dλ。LAMMPS 用 fix adapt 或混排势参数实现 λ 耦合，每 λ 跑平衡后统计 ∂U/∂λ 的平均值。关键认知：积分曲线的端点最容易因截断/长程处理产生尖峰，需要加密端点窗口。我用梯形或 Simpson 积分时，默认端点附近 λ 间距减半，避免翘尾污染总面积。

三、关键技术要点：长程与混排

lammps计算自由能最致命的是长程库仑（PPPMC 或 Ewald）在 λ 端点因原子”消失/出现”导致的能量不连续。我做原子替换时会对被耦合原子用软核或逐步缩小电荷，平滑过渡。另一个坑是混排势（如 Lennard-Jones 的 λ 混排）在端点附近出现奇点，需要软核半径调大。曾经没调软核，λ=0.02 处能量发散，积分直接废掉。

四、实操流程：从势函数到 ΔG

我的标准动作：定义 λ 耦合方案（软核/混排）→ 规划 λ 序列（端点加密）→ 每 λ 平衡+生产 → 统计 ∂U/∂λ → 数值积分 + 误差估计 → 和实验或 DFT 对标。最容易被砍的是端点加密，但 lammps计算自由能 的积分质量全靠它，我默认保留。

补一个实算：Au-Cu 固溶形成自由能，均匀 λ 间距（11 点）算得 0.12 eV/atom，端点加密到 21 点后变成 0.09，和 DFT 热力学 0.08 接近。如果只交 11 点结果，你会高估固溶倾向，相图预测就偏。这之后我做 TI，窗口数先按”端点 ∂U/∂λ 是否平滑”倒推。

五、常见踩坑：截断与有限尺寸

lammps计算自由能常见错是势截断半径在 λ 变化时没保持一致，导致不同窗口能量基准漂移。另一个坑是盒子太小，受限体系的表面能混进体相自由能，需要做大盒或做表面修正。还有人忽略有限尺寸对熵的贡献，小体系算出来的自由能比热力学极限差一截。

六、复盘：材料自由能是相图的命门

回过头看，lammps计算自由能最值钱的是把相稳定性量化到可喂进相图计算，但前提是你把 λ 过渡、长程、端点三件事处理干净。我接项目默认先问有没有势垒、盒子多大，再定方法。把积分路径做扎实，ΔG 才敢拿去和 CALPHAD 对标。被证明有用的，是愿意为端点窗口多跑两档的耐心。下次有人给我一个”LAMMPS 算的形成能”，我会先问——你的 λ 端点和长程处理查了吗。