[论文解读] Direct calculation of the solid-liquid Gibbs free energy difference in a single equilibrium simulation

Ulf R. Pedersen|arXiv (Cornell University)|Mar 21, 2013

Advanced Physical and Chemical Molecular Interactions参考文献 2被引用 34

一句话总结

本文提出界面钉扎（IP）方法，通过施加谐振势场以稳定两相共存构型，在单次平衡分子动力学模拟中直接计算固-液吉布斯自由能差。该方法通过计算势场对系统的平均作用力来确定自由能差，结合牛顿-拉夫森根求解法，可实现对Lennard-Jones模型熔点线的高精度测定，且滞后效应和有限尺寸效应极小。

ABSTRACT

Computing phase diagrams of model systems is an essential part of computational condensed matter physics. In this paper we discuss in detail the interface pinning (IP) method for calculation of the Gibbs free energy difference between a solid and a liquid. This is done in a single equilibrium simulation by applying a harmonic field that biases the system towards two-phase configurations. The Gibbs free energy difference between the phases is determined from the average force that the applied field exerts on the system. As a test system we study the Lennard-Jones model. It is shown that the coexistence line can be computed efficiently to a high precision when the IP method is combined with the Newton-Raphson method for finding roots. Statistical and systematic errors are investigated. Advantages and drawbacks of the IP method are discussed. The high pressure part of the temperature-density coexistence region is outlined by isomorphs.

研究动机与目标

• 开发一种稳健、高效的固-液吉布斯自由能差计算方法，无需依赖非平衡或多模拟方法。
• 克服传统方法（如Widom插入法、热力学积分法、反应路径积分法）在高密度或复杂体系中的局限性。
• 通过单次平衡模拟直接计算共存线，实现相图的精确构建，且滞后效应极小。
• 研究IP方法中的统计误差与系统误差，评估其在Lennard-Jones流体等模型体系中的可靠性。
• 展示该方法在高压区域的适用性，及其在研究液晶或准晶等复杂相态中的潜力。

提出的方法

• 界面钉扎（IP）方法通过施加与区分固相与液相的序参量耦合的谐振势场，在单次平衡模拟中稳定两相共存构型。
• 吉布斯自由能差通过谐振场对系统施加的平均力计算得出，依据热力学关系 ΔG = ⟨F⟩ × ΔQ。
• 序参量 Q 基于长程有序性定义，利用局域结构分析将粒子分类为固相、液相或界面相。
• 采用牛顿-拉夫森方法求解方程 ⟨Q⟩ = N/2，该方程对应共存条件，从而高效计算熔点线。
• 通过改变系统尺寸并观察 Δμ 对界面位置参数 a 的依赖关系，分析有限尺寸效应和界面迁移率带来的系统误差。
• 该方法应用于截断Lennard-Jones模型，对截断效应引起的势能偏差进行了修正，以提高精度。

实验结果

研究问题

• RQ1 是否可在单次平衡模拟中直接、准确地计算固-液吉布斯自由能差，而无需非平衡或多相模拟？
• RQ2 在高密度体系中，IP方法与Widom插入法或热力学积分法相比，在精度与效率方面表现如何？
• RQ3 IP方法中的主要统计误差与系统误差来源是什么？如何有效缓解？
• RQ4 IP方法在高压条件下是否仍保持高精度？是否可利用等构线（isomorphs）描述高压共存区域？
• RQ5 该方法能否推广至其他相变过程（如固-固或液-液相变）？在这些情况下存在哪些限制？

主要发现

• IP方法可从单次平衡模拟中直接计算固-液吉布斯自由能差，无需依赖非平衡或多模拟协议。
• 结合牛顿-拉夫森根求解算法，该方法在计算截断Lennard-Jones模型的熔点线时表现出高精度。
• 统计误差可控，有限尺寸效应引起的系统误差与界面位置参数 a 相关，尤其在小系统中表现明显。
• 温度-密度共存区域的高压部分可通过等构线良好描绘，表明该方法在探索复杂相行为方面的有效性。
• 该方法对滞后效应具有强鲁棒性，且可直接从同一模拟中计算界面性质（如表面能、晶体生长速率）。
• IP方法被证明是Widom插入法与热力学积分法的可行替代方案，尤其在高密度或复杂相行为体系中，后两者方法可能失效。