The Free-Energy Barrier of Precritical Nuclei in Hard Spheres is Consistent with Predictions

技术摘要：硬球体系中前临界核的自由能垒

问题陈述
预测晶体成核速率仍然是凝聚态物理学中的一个重大挑战。在典型的硬球胶体体系中，实验测量值与计算机模拟之间存在巨大的差异。具体而言，在低过饱和度区域（即成核缓慢的区域），通过静态光散射推断出的实验成核速率与模拟预测值之间存在超过 10 个数量级的差异。这一差距至关重要，因为用于解决该问题的稀有事件采样技术是分子科学的基础，但这些技术却极少直接通过实验进行验证。以往试图利用粒子解析技术（共聚焦显微镜）来解决这一问题的尝试，受限于用于成像的大尺寸胶体极其缓慢的动力学特性，这使得在实验时间尺度内观察弱过饱和悬浮液中的稀有成核事件变得不可能。

方法论
作者通过使用密度匹配溶剂中的立体稳定 PMMA 球体（直径 2.0 µm，多分散性 4%）进行实空间共聚焦显微成像，开发了两种主要手段来应对这一挑战：

1. 高精度状态点映射： 为了确保实验与模拟之间的有效比较，作者开发了一种严谨的方法，将实验体积数 $\phi$ 映射到有效硬球体积数 $\phi_{eff}$。

   • 他们采用了一种深度学习粒子追踪算法（“Colloidoscope”），能够提取约 98% 的粒子坐标，且定位误差约为粒子直径的 5%，克服了传统算法在拥挤环境中的局限性。
   • 作者并未仅仅依赖于对体积数变化不敏感的径向分布函数 $g(r)$，而是利用了高阶相关函数。具体而言，他们分析了通过拓扑簇分类（TCC）识别出的特定几何簇的群体，TCC 能够最小化硬球体系中的局部自由能。
   • 通过对总对相关函数 $h(r)$ 和 TCC 簇群体进行加权比较，他们实现了前所未有的高精度状态点匹配（$\delta\phi < 0.0025$）。

2. 前临界核自由能的测量： 考虑到在差异区域内的临界核过于庞大且罕见，无法直接观察，作者将研究重点转向了前临界核。这些较小的簇形成频率足够高，可以在实验中被采样。

   • 他们利用观测到的核浓度，计算了作为函数形式的有效形成自由能 $\Delta G(n)/k_B T$，其中 $n$ 为核的大小：$\Delta G(n)/k_B T = -\ln(N_n/N)$。
   • 为了确保公平比较，模拟数据经过了后处理，以模拟实验条件，包括有限的成像体积、位置误差（高斯位移）以及“丢失”的粒子。
   • 核的大小通过 $q_6$ 键取向序参数（Ten Wolde 等人）进行确定，其连通性截断值为 0.7，且要求至少有 5 个类固体邻居。

关键结果

• 状态点一致性： 实验与模拟的径向分布函数及高阶簇群体表现出极佳的一致性，证实了实验体系已准确映射到硬球模型。带有残余静电相互作用（硬核尤卡瓦势）的模拟与实验簇群体不符，从而验证了硬球近似的有效性。
• 自由能一致性： 在多个体积数（$\phi = 0.495$ 至 $0.520$）下，测得的前临界核自由能剖面（$\Delta G_{obs}$）与无偏及有偏计算机模拟高度吻合。
  • 对于 $\phi \geq 0.530$ 的情况，观察到了约 $k_B T$ 的微小差异。作者指出，这对应于仅为 $e$ 倍（约 2.7）的成核速率差异，与之前光散射研究中观察到的 10 个数量级的差异相比，这种差异可以忽略不计。
• 核的性质： 对核几何形状的分析表明，实验与模拟之间的非球形度和回转半径（$R_g$）相匹配。较大的核趋向于更趋于球形，且 $R_g$ 随尺寸变化的比例关系遵循与模拟一致的幂律。

意义与主张
本文声称，这是首次在任何材料中成功实现核形成自由能的实验与模拟对比。其主要意义在于验证了广泛应用于模拟中的稀有事件采样技术；这种一致性表明这些计算方法是稳健且准确的。

作者认为，历史上的成核速率差异可能源于两个因素：

1. 状态点映射不准确： 先前的研究依赖于间接方法（如通过相图推导）来确定体积数，这可能受到成核、相分离和沉降等竞争时间尺度的影响。
2. 间接核检测： 静态光散射间接推断核的大小和数量，可能高估了核的计数。

通过使用直接实空间成像和精确的状态点映射，作者证明了硬球体系中前临界核的自由能垒与理论预测是一致的。尽管作者对向临界势垒高度（$\Delta G^*$）的外推保持谨慎，但前临界区域差异的消除表明，硬球成核的基础物理机制已被充分理解，而之前的巨大差异很可能是实验方法学和状态点确定过程中的人为误差。