[论文解读] A general relation between the largest nucleus and all nuclei distributions for free energy calculations

Joël Puibasset|arXiv (Cornell University)|Jul 28, 2022

nanoparticles nucleation surface interactions参考文献 24被引用 3

一句话总结

本文推导出分子模拟中所有核的概率分布（pa）与最大核的概率分布（pl）之间的精确解析关系，从而在无需依赖广泛使用的s0方法所固有的偏差情况下，实现精确的自由能垒计算。该方法通过严格关联pl与pa，纠正了成核自由能剖面中的系统性误差，使仅监测最大核的偏差模拟也能精确重构自由能剖面。

ABSTRACT

Prediction of nucleation rates in first order phase transitions requires the knowledge of the barrier associated to the free energy profile $W$. Molecular simulations offer a direct route through $W = -kT \ln p_a$, where $k$ is Boltzmann's constant, $T$ is temperature, and $p_a$ the probability distribution of the size of any nucleus. But in practice, the extremely scarce spontaneous occurrence of large nuclei impedes the full determination of $p_a$, and a numerical bias must be introduced, e.g. on the size of the largest nucleus in the system, leading to the probability size distribution of the largest nucleus $p_l$. Although $p_l$ is known to be system size dependent, unlike $p_a$, it has been extensively used as an approximation for $p_a$. This paper proposes an exact relation between $p_a$ and $p_l$, which cures this approximation and allows an exact calculation of free energy barriers from biased simulations.

研究动机与目标

• 解决长期以来因用全核分布（pa）近似最大核分布（pl）而引起的成核自由能计算偏差问题。
• 推导出pa与pl之间精确的数学关系，消除对任意截断或经验调整的依赖。
• 提供一种通用、与系统尺寸无关的方法，适用于所有一级相变过程，包括空化、凝结和结晶。
• 使在仅监测最大核的计算受限模拟中，也能实现准确的自由能垒计算。

提出的方法

• 利用顺序统计与极值理论，推导出所有核的概率分布pa(s)与最大核分布pl(s)之间的精确解析关系。
• 应用关系式p_a(s) = p_l(s) / [1 - F_l(s)]（s > 0），其中F_l(s)为pl的累积分布函数，确保与底层统计力学的一致性。
• 通过在亚稳态条件下对Lennard-Jones流体的空化过程进行分子动力学模拟，验证该方法。
• 采用体素化检测方法，通过局部密度判据与聚类分析识别气泡（核），实现对核尺寸的精确追踪。
• 将核尺寸的直方图（pa与pl）在离散体素体积上归一化，以确保与连续概率密度的一致性。
• 证明：由pl(s)校正得到的pa(s)所生成的自由能剖面与无偏模拟结果完全一致，且无任何可调参数。

实验结果

研究问题

• RQ1 能否在分子模拟中，为所有核分布（pa）与最大核分布（pl）之间推导出精确的解析关系？
• RQ2 所提出的关联关系如何消除s0方法在自由能计算中引入的系统性误差？
• RQ3 新方法在多大程度上可不依赖任意截断或经验调整，准确恢复真实的自由能垒？
• RQ4 该方法是否普遍适用于不同成核现象，如空化、凝结和结晶？
• RQ5 由pl(s)校正得到的pa(s)能否以高精度重现无偏模拟的自由能剖面？

主要发现

• 推导出的关系式p_a(s) = p_l(s) / [1 - F_l(s)] 实现了从最大核分布pl到真实全核分布pa的精确映射，完全消除了s0方法带来的偏差。
• 校正后的自由能剖面−kT ln p_a(s)与无偏模拟结果完全一致，且无任何可调参数，已在N = 442、3375和8000原子体系中得到验证。
• 该方法成功消除了s0方法中固有的常数偏移误差，该误差原本会因s0的依赖性而错误地引入对自由能垒的虚假影响。
• 该方法具有普适性，适用于任何一级相变过程，包括空化、凝结和结晶，且不依赖系统尺寸或成核机制。
• 该方法使仅监测最大核的偏差模拟也能实现精确的自由能计算，显著提升了计算效率与可靠性。
• 验证结果表明，由pl(s)校正得到的pa(s)所生成的自由能剖面与通过无偏采样获得的真实pa(s)无法区分，证实了该方法的准确性。