本文发表于 Physical Review Letters 期刊，出版于 2025年4月7日。该研究由Camilla Beneduce（第一作者）以及Diogo E. P. Pinto, Lorenzo Rovigatti, Flavio Romano, Petr Šulc, Francesco Sciortino, John Russo 等学者共同完成，研究机构主要来自意大利罗马大学、威尼斯大学以及美国亚利桑那州立大学等。
本研究的科学领域属于凝聚态物理、化学物理与材料科学的交叉领域，特别是相变成核理论与软物质自组装。研究的核心动机在于挑战并检验经典成核理论（Classical Nucleation Theory， CNT）的基石——毛细管近似（capillarity approximation）的普适性。
背景知识与研究动机：经典成核理论是描述从过冷液体或过饱和溶液中析出晶体等相变过程的奠基性理论。其核心假设（毛细管近似）认为，临界晶核的热力学性质（如界面自由能、密度）与宏观块体晶相相同，从而可以利用块体性质计算成核自由能垒和成核速率。尽管CNT因其简洁性和在某些情况下的成功应用而被广泛使用，但大量实验和模拟结果常常显示出与CNT预测的显著偏差。通常，这些偏差通过引入界面能随温度变化等唯象扩展来“修补”，这引发了一场关于CNT根本有效性的持续辩论。Oxtoby曾尖锐地指出：“成核理论是少数几个理论预测速率与测量速率在几个数量级内符合就被认为是重大成功的科学领域之一”，这凸显了精确检验CNT的困难所在。
以往检验CNT的方法多集中在直接比较预测与测量的成核速率上，但这种方法严重依赖于在成核条件下难以精确测量的量（如界面自由能γ），且结果可能因体系或状态点而异，难以得出普适性结论。
研究目的：本研究团队旨在提出一种对CNT乃至任何基于毛细管近似的成核理论的“可证伪性测试（falsifiability test）”。他们不追求精确预测成核速率，而是转换视角，聚焦于材料的多晶型（polymorphic） 组成。具体而言，他们设计了一种理想情况：在一个体系中，几种不同的晶体多晶型体具有完全相同的块体自由能和界面自由能。根据CNT的毛细管近似，这些多晶型体应该具有完全相同的成核性质。如果实验或模拟结果显示它们具有显著不同的成核行为，那么CNT的基本假设就被证伪了。本研究的核心目的，就是构建这样一个体系，并通过分子模拟演示这种与CNT预测的直接矛盾。
研究流程主要包括四个关键步骤：1）模型体系的设计与多晶型体的识别；2）多晶型体块体与界面性质的验证；3）直接成核模拟；4）自由能垒计算与熔体结构分析。
第一，模型体系的设计与多晶型枚举。 研究团队设计了一个二元四价补丁粒子（binary tetravalent patchy particles） 混合物模型，称为 N2C8（意为两种粒子、八种补丁类型）。如图1所示，两种粒子（青色和红色）表面有四个四面体排列的补丁，它们根据一个特异的相互作用矩阵（图1c）进行键合。这种特异性相互作用模拟了类似DNA碱基配对的可编程自组装，可通过DNA折纸技术实现。 该模型最初的设计目标是专门自组装成立方金刚石（cubic diamond）结构，避免形成六方金刚石。利用SAT组装算法（SAT-assembly algorithm），他们系统性地枚举了在48个粒子的立方金刚石晶格中，N2C8粒子所有可能的周期性排列方式。结果发现了三种不同的周期性排列，即三种多晶型体（图2）。这三种结构同属P1空间群（仅有平移对称性），但它们的单胞大小（即晶体中物种排列的最小重复单元）不同。研究团队据此将其命名为：dc-8（单胞8个粒子）、dc-16和dc-24。需要强调的是，在考虑的成核条件下，临界晶核尺寸远小于30个粒子，因此单胞大于48个粒子的更大周期结构未出现。
第二，验证多晶型体的热力学等同性。 这是检验的关键前提。研究者首先确认，这三种多晶型体都是等化学计量的，键合数量与类型完全相同，且晶格振动特性也一致，因此它们具有完全相同的块体自由能。接着，他们通过伞形采样（umbrella sampling） 和直接共存（direct coexistence） 模拟，精确计算了每种多晶型体与同一熔体相之间的固-液界面自由能。结果如图2(d)所示，三种多晶型体的界面自由能在约3%的误差范围内完全相同。至此，支撑CNT计算的所有宏观性质（化学势差δμ和界面能γ）对于三种多晶型体都是等同的。根据CNT，它们应当具有完全相同的成核自由能垒和速率。
第三，直接成核蒙特卡洛模拟。 为了观察实际的成核行为，研究者在两个选定的状态点（T=0.1，ρ=0.35 和 T=0.104，ρ=0.4）下，对包含500个等化学计量比粒子的体系进行了大量正则系综蒙特卡洛模拟。为了提高采样效率，部分模拟使用了聚集体积偏差（Aggregation-Volume Bias， AVB）移动算法。总共运行了数百条独立的模拟轨迹。为了识别和分析晶体，他们采用了基于球谐函数的总相干性（total coherence） 序参量来区分液相和晶相，并开发了一种新的序参量来精确区分三种不同的多晶型体。每条轨迹结束时，通过判断立方金刚石相粒子占比是否超过0.5来确认成核成功，并对最终晶体进行分类。
第四，自由能垒计算与熔体结构分析。 为了定量比较成核难度，研究者采用伞形采样技术，在固定温度和压力（T=0.104， P=0.018）下，通过施加偏向势，计算了形成不同尺寸晶核的自由能。他们小心地制备了不同尺寸（n0）的晶核作为种子，确保模拟过程中不发生自发成核或晶型转变，且相邻窗口有足够的重叠以保证采样的连续性。此外，为了探究成核差异的根源，他们分析了熔体中的结构涨落。他们定义了一个角度参数 α（见图4b插图），用于描述同一物种的次近邻粒子之间的补丁取向关系。在三种块体多晶型中，α值不同（dc-24约为109°，dc-16和dc-8约为9°）。他们从伞形采样模拟中提取了包含约50个粒子的晶核构型，计算了从晶核质心向外径向距离上的α值和总相干性序参量，以考察从晶核内部到界面再到熔体的结构变化。
直接成核模拟结果（图3a）与CNT预测截然相反。在T=0.1， ρ=0.35条件下，600条模拟轨迹中，dc-24（单胞最大）的成核事件最多（有AVB时43次，无AVB时30次），dc-16次之（15次和11次），而dc-8（单胞最小）几乎不形成（仅有1次无AVB的轨迹）。第二个状态点的结果类似。这表明，三种热力学性质完全等同的多晶型体，其成核速率存在数量级上的显著差异。图3b展示了两条典型轨迹的晶核尺寸随时间演化，可以看出dc-24和dc-16在临界晶核阶段就已呈现出结构差异。额外的生长分析排除了“交叉成核”（一种多晶型在另一种表面生长）的可能性，确认了成核过程的独立性。
自由能垒计算结果（图3c）为成核速率的差异提供了定量支撑。计算显示，dc-24的成核自由能垒最低，dc-16次之，dc-8的能垒最高。dc-24与dc-16之间的能垒差约为1 k_BT，这与它们成核轨迹数量的比值基本吻合。这一结果直接证实了成核行为的差异源于热力学能垒的不同，而非动力学偶然性。所有临界晶核尺寸nc均小于30个粒子，小于dc-24的单胞尺寸（48粒子）。
熔体结构分析结果（图4）揭示了成核差异的微观起源。图4a展示了总相干性的径向分布，清晰标出了晶核、界面和熔体区域。图4b的α角度径向分布图揭示了关键信息：在晶核核心区，α角度取各自块体多晶型的特征值（dc-24核为~109°，dc-16/dc-8核为~9°）。然而，在界面区域和熔体区域，无论所观察的晶核本身是哪种类型，α角度都稳定在dc-24的特征值（~109°）附近。这表明，熔体本身的局域取向序更接近于dc-24多晶型。也就是说，熔体的结构涨落并非各向同性，而是存在一种偏向性，这种“预组织”或“结构记忆”使得形成与熔体局部结构更匹配的dc-24晶核所需的熵代价更小，从而其成核能垒最低。补充材料中的进一步分析还排除了晶核表面存在特殊有利或不利键合位点的影响，进一步将原因指向了熔体本身的结构特性。
本研究的主要结论是：经典成核理论（CNT）的毛细管近似在解释多晶型选择问题上存在根本性缺陷。研究者通过设计一个精巧的“可证伪性测试”体系——其中不同多晶型体的块体和界面自由能完全相同——并利用分子模拟，无可辩驳地展示了这些多晶型体具有显著差异的成核速率和自由能垒。这一结果直接违反了CNT基于毛细管近似的预测，因为该理论无法解释这种差异。
研究的深层含义在于指出了CNT的一个主要局限：忽视了液相（熔体）内部的结构涨落。实验表明，正是熔体中存在的、与某种特定多晶型（在本研究中是dc-24）更接近的局域取向序，决定了哪种多晶型更容易成核。成核过程并非仅仅由宏观块体晶相的性质决定，而是强烈依赖于成核发生时的具体液相环境。这为超越CNT的理论框架提供了强有力的支持，这些新理论（如考虑液序参数、两步成核等）将液相的结构细节纳入了考量。
本研究的价值在于： 1. 理论价值：提出了一种新颖且强有力的理论检验范式（可证伪性测试），将焦点从难以精确定量的绝对成核速率预测，转移到更易于判断的相对多晶型排序上，避免了对精确界面能等参数的依赖。 2. 科学价值：清晰地揭示了CNT毛细管近似的失效场景，将问题的根源指向了液相结构，深化了对成核初期微观机制的理解，特别是对多晶型选择这一关键材料科学问题的理解。 3. 方法论价值：展示了结合“可编程”补丁粒子模型、SAT组装算法、先进增强采样模拟（伞形采样、AVB）和精细序参量分析，可以用于设计和研究高度复杂的自组装与成核问题。