本文档报告了由Yijian Wu, Thomas Philippe, Aymane Graini以及Julien Lam共同完成的一项原创性研究。该研究发表于2026年1月5日的《物理评论快报》（*Physical Review Letters*）期刊上，标题为“nonclassical nucleation pathways in liquid condensation revealed by simulation and theory”。主要研究人员来自法国巴黎综合理工学院凝聚态物理实验室（CNRS, Ecole Polytechnique, Institut Polytechnique de Paris）以及法国里尔大学（Université Lille）的材料与转化单元（UMET）。
研究学术背景 此项研究聚焦于凝聚态物理和材料科学中的一个经典且关键的问题：成核（nucleation）过程，特别是液相从过饱和气相中凝结的均相成核。作为所有一级相变的触发机制，对成核的深入理解在材料工程、药物开发、气候控制等多个领域具有根本性重要性。长期以来，描述成核的主流理论是经典成核理论（Classical Nucleation Theory, CNT）。CNT以其简洁性著称，它假设成核过程仅由一个序参量——团簇尺寸——驱动，并采用毛细管近似（capillary approximation）模型来描述团簇的表面能。然而，尽管CNT提供了对成核的总体理解，但其在定量预测成核速率和临界团簇性质方面，无论是在实验还是数值模拟中，常常存在巨大偏差。 为了改进CNT，以往的研究工作主要集中在两个方面：一是修正毛细管模型，例如引入表面张力随团簇尺寸变化的托尔曼（Tolman）修正；二是发展弥散界面模型（diffuse-interface model），如基于方梯度近似（square-gradient approximation, SGA）的模型，它们能更好地与微观模拟结果定性吻合，但其高度多维的自由能面使得计算临界团簇性质和成核速率变得异常复杂和困难。在实验层面，杂质的存在（导致异相成核）和成核事件位置难以控制是主要挑战。在模拟层面，成核是一个稀有事件，在低过饱和度下通过常规的“暴力”（brute-force）分子动力学模拟捕获临界团簇极其困难且计算成本高昂。 因此，本研究旨在克服这些挑战，重新审视液相凝结这一经典问题。其核心目标是：1）利用先进的稀有事件采样模拟技术，在原子尺度上精确表征宽过饱和度范围内的临界成核团簇结构特征；2）基于模拟发现，构建一个既物理准确又易于计算的多维成核理论框架，以超越CNT的单变量描述局限，为成核过程提供更普适的理解和预测工具。
详细研究流程 本研究采用了“模拟-理论”紧密结合的研究范式，包含两个主要且相互关联的部分：高精度数值模拟和新型理论模型的构建与验证。
第一部分：稀有事件采样模拟 研究的首要任务是获取精确的临界成核团簇数据。研究对象是过饱和态下的伦纳德-琼斯（Lennard-Jones）气体系统，温度设定为T=0.8ε/kB。为了在宽泛的初始气相密度（ρ0）范围内（即不同过饱和度下）捕获临界团簇，研究者根据过饱和度的不同，创新性地组合并应用了两种互补的稀有事件采样技术，而非单一的“暴力”模拟。
低过饱和度区间（ρ0 < 0.015σ⁻³）：在此区域，预期的临界团簇尺寸较大。研究者采用了播种（Seeding）技术的一个特定实施方案。具体操作是：将一个预设半径和密度的球形液滴人工置入气相环境中。通过仔细选择系统的平均密度和模拟盒子尺寸，利用受限系统中的“超稳定化”（superstabilization）效应，液滴会自发地调整其半径和密度，并稳定在一个特定的状态。同时，气相密度也会通过有限尺寸效应和质量守恒调整到一个新的平衡值ρ0。在此条件下，所得到的液滴即为该ρ0对应的临界成核团簇。通过改变初始速度和原子位置，可以生成一组临界团簇的统计样本。这种方法的关键在于巧妙地利用系统约束来“固定”住临界团簇结构。
中等过饱和度区间（ρ0 ∈ [0.02σ⁻³; 0.04σ⁻³]）：此时临界团簇尺寸太小，播种技术不再适用。研究者转而采用一套结合能量偏置和轨迹偏置的增强采样技术组合。
- 第一步：导向分子动力学（Steering Molecular Dynamics）。使用配位数作为集体变量，对系统施加能量偏置，“推动”气体克服自由能势垒形成液滴。
- 第二步：承诺概率（Commitment Probability）计算。从上述能量偏置轨迹中提取中间结构，然后对每个结构进行多轮不同初始速度条件的“暴力”模拟，计算其分别演化为稳定液相或返回气相的“承诺概率”。
- 第三步：无目标射击（Aimless Shooting）。选取一条连接气相和液相的“暴力”模拟轨迹作为初始轨迹，启动无目标射击模拟。这种方法能高效地生成大量连接气相和液相的过渡路径，并从中精确地提取出位于自由能鞍点（saddle point）的临界核。
通过耦合这两种策略，研究团队针对每个研究的气相密度，成功获得了大约20个定义明确的临界团簇结构，从而能够精确地以密度和尺寸为指标来表征成核事件。所有分子动力学模拟均使用了LAMMPS等标准或改进的软件工具，具体的模拟和分析协议详见补充材料。
第二部分：二变量成核理论的构建 基于模拟获得的关键洞察——临界团簇的密度显著偏离宏观平衡液体的密度，研究团队构建了一个新颖的理论框架：二变量成核理论。该理论是对经典毛细管模型的扩展，其核心创新在于将团簇密度（ρ）与团簇尺寸（半径R）共同作为驱动成核过程的序参量。
模型构建：
- 形成功（Work of Formation）：团簇的形成功δω(R, ρ)被表达为体积自由能项和表面能项之和：δω = -(4πR³/3) * gₙ + 4πR²γ。其中，驱动势gₙ和表面张力γ都被处理为团簇密度ρ的函数。驱动势gₙ通过亥姆霍兹自由能密度f(ρ)计算得出。表面张力采用受方梯度近似启发的形式：γ = κ(ρ - ρ₀)²，其中κ为常数，由平衡态下的平面界面表面张力确定。此模型在初始阶段忽略了表面张力随尺寸的变化（即托尔曼修正）。
- 动力学方程：成核动力学通过描述团簇密度函数f(R, ρ, t)演化的福克-普朗克方程（Fokker-Planck equation）来研究。成核通量向量由热力学驱动力（自由能梯度）和动力学矩阵（扩散性矩阵D）共同决定。在鞍点（临界点Rc, ρc）附近，演化方程线性化为式（1），其中矩阵Z = D·H（H为海森矩阵）编码了热力学（H）与动力学（D）的复杂相互作用。
- 稳态成核率：基于朗格（Langer）和阿列克谢奇金（Alekseechkin）的理论框架，推导出稳态成核率I的表达式（式2），该速率取决于临界形成功δωc、矩阵Z的负特征值λ₁^z及其对应的不稳定方向（即成核通量方向）。
关键参数确定：
- 分布参数f₀：为解决高维相空间中f₀确定的难题，研究者提出了一个创新方法：利用亚稳态区域的质量守恒条件（式3）来计算f₀。这种方法避免了传统理论中的苛刻假设，使得计算在二维框架下变得可行且可归一化。
- 扩散性矩阵D：基于动态密度泛函理论（Dynamical Density-Functional Theory, DDFT）在鞍点附近的公式推导得出，其核心假设是粒子扩散的驱动力正比于实空间扩散常数D_diff和化学势梯度。D_diff的值通过分析播种模拟中临界后团簇早期生长阶段的动力学数据估计得出。
- 热力学输入：均相伦纳德-琼斯流体的亥姆霍兹自由能密度通过微扰理论计算获得，其预测的平衡密度与微观模拟结果高度一致。
- 表面能因子κ：由平衡温度下的表面张力数据确定。
对比模型：为了评估新模型的表现，研究结果与经典成核理论（CNT）以及方梯度近似弥散界面模型（SGA）的预测进行了系统比较。
主要研究结果 模拟和理论的结果相互印证，共同揭示了一种非经典的成核路径，并验证了新理论模型的强大预测能力。
临界团簇的结构特征：
- 模拟发现：图1(a-c)清晰地展示，临界团簇的密度并非恒定，而是随着过饱和度（初始气相密度ρ0）的增加而显著下降。这与经典成核理论假设临界团簇具有平衡液相密度的前提直接矛盾。同时，临界团簇的半径也随过饱和度增加而减小。
- 理论对比：
- 尺寸：所有模型（二变量理论、CNT、SGA）对临界半径的预测趋势相似，在低过饱和度下与模拟数据吻合良好，但在中等过饱和度下均略微低估了尺寸。
- 密度：CNT完全无法捕捉密度变化趋势，它预测的密度恒为平衡值。SGA模型虽然能定性地预测密度下降，但定量上高估了密度值，并且其预测曲线存在一个在模拟中未观察到的极大值。唯有本研究提出的二变量成核理论，能够同时在定性和定量上准确地重现模拟观察到的临界团簇密度和尺寸变化规律（图1）。
稳态成核率：
- 图2展示了理论预测的成核率与文献报道及本研究自身“暴力”模拟（采用Yasuoka-Matsumoto阈值法分析）结果的比较。
- 二变量理论在广泛的过饱和度范围内，尤其是CNT失效的高过饱和度区域，表现出卓越的预测精度，其预测与模拟数据高度吻合，误差在数个数量级以内。这证明了该理论在耦合了准确的热力学和动力学数据后，具备强大的、无拟合参数的预测能力。
- 此外，二变量理论预测的成核率在接近旋节线（spinodal）极限时会出现一个最大值，这一效应也在弥散界面模型中被观察到，但CNT无法预测。
- 研究还探讨了托尔曼修正（引入尺寸相关的表面张力γ® = γ∞(1 - 2δ/R)）的影响。结果显示，虽然它对临界团簇性质影响甚微，但由于成核率对形成功呈指数依赖，因此对成核率的预测有显著影响。分析表明，要保持与分子动力学结果的良好一致性，托尔曼长度δ必须保持较小的值（约0.05σ-0.1σ），这与先前研究在类似条件下的结果一致。
临界形成功：
- 图3对比了不同模型预测的成核自由能势垒δωc。在低过饱和度下，三个模型结果一致，因为此时临界密度接近平衡密度。在中等过饱和度区域，二变量理论和SGA模型预测的势垒远低于CNT的预测。这是因为形成一个密度较低的团簇虽然体积驱动力较小，但其表面能成本也大幅降低，总体上使得成核更容易发生。接近旋节线极限时，二变量理论预测临界团簇密度趋近于气相密度，团簇尺寸发散，形成功趋于零，平滑地过渡到旋节分解（spinodal decomposition）机制。CNT则完全无法描述这种向旋节分解的过渡。
成核路径与通量方向：
- 这是本研究最关键的发现之一。在CNT中，热力学和动力学都驱使团簇在恒定（平衡）密度下生长。然而，二变量理论预测，随着过饱和度增加，成核通量方向会偏离自由能最速下降方向。
- 图4展示了在ρ0 = 0.03σ⁻³条件下，理论计算的自由能面、预测的成核通量方向（红色虚线）以及从模拟中提取的四条生长轨迹。模拟轨迹清晰地显示，成核路径并非一条水平线，而是同时伴随着团簇的“生长”和“致密化”。更重要的是，模拟观察到的路径与二变量理论预测的通量方向高度一致，而非与能量最速下降方向一致。这直接证实了所提出的理论成功地捕获了非经典的、密度与尺寸协同演化的成核路径。
研究结论与意义 本研究通过整合创新的模拟策略和理论构建，重新审视并深化了对液相凝结成核过程的理解，取得了以下核心结论和价值：
揭示非经典成核路径：原子模拟首次以高精度证实，在液相凝结过程中，临界成核团簇的密度显著不同于宏观平衡液体，且随过饱和度变化。成核路径涉及团簇尺寸增长与密度增加的协同过程，这挑战了经典成核理论将密度视为常数的基本假设。
建立有效的二维理论框架：基于上述发现构建的二变量成核理论，成功地将团簇密度作为一个关键序参量纳入经典的毛细管近似模型中。该模型无需拟合参数，便能定量预测宽过饱和度范围内的成核速率、临界团簇尺寸和密度，并定性且定量地重现了模拟观察到的非经典成核路径。
指明经典理论的局限与改进方向：研究明确指出，忽略团簇密度作为一个集体变量是导致经典成核理论长期存在定量失败的主要原因之一。团簇密度是控制成核过程的一个决定性变量。
提供普适性研究框架：本研究建立的“数值模拟-理论构建”一体化框架具有广泛的适用性。它证明了通过结合先进的稀有事件采样技术和多维但易于处理的成核理论，可以克服传统CNT的局限，并为理解更广泛的相变成核问题（如结晶、气泡成核等）提供了坚实的理论基础和方法论借鉴。
研究亮点
- 方法创新：创新性地组合并应用了播种技术与导向分子动力学结合无目标射击的增强采样方案，成功攻克了在全过饱和度范围内（尤其是低过饱和度下）精确表征临界成核团簇这一长期存在的模拟难题。
- 理论突破：提出了一个物理图像清晰、计算上易处理的二变量成核理论。该理论的关键创新在于：a) 在毛细管近似中显式引入可变密度；b) 利用质量守恒条件巧妙地解决了高维相空间中分布参数f₀的确定问题；c) 通过福克-普朗克方程和DDFT自然地耦合了热力学与动力学。
- 重要发现：首次通过模拟直接观测并理论阐明了液相凝结中“生长与致密化同步”的非经典成核路径，定量揭示了临界团簇密度对过饱和度的强依赖性，并成功预测了成核率在接近旋节线时的最大值行为。
- 学科价值：这项研究不仅显著推进了对成核这一基础物理过程的理解，而且其发展的理论框架和模拟策略为材料科学、化学工程、大气科学等众多依赖成核过程的领域提供了更精准的理论预测工具和模拟分析范式。