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

技術要約：硬球における前臨界核の自由エネルギー障壁

問題提起
結晶核生成速度の予測は、凝縮系物理学における重要な課題である。典型的な系である硬球コロイドにおいて、実験による測定値とコンピュータシミュレーションの間には、極めて大きな乖離が存在する。具体的には、低過飽和領域（核生成が遅い領域）において、静的光散乱から推定される実験的な核生成速度は、シミュレーションの予測値と10桁以上の差がある。このギャップは、希少事象サンプリング手法（rare-event sampling techniques）が分子科学において極めて重要であるにもかかわらず、それらが実験に対して直接検証されたことがほとんどないという点において、非常に重要である。粒子分解能技術（共焦点顕微鏡）を用いたこれまでの試みは、画像化のために必要な大きなコロイドの運動性が極めて低いという制約により、弱過飽和懸濁液における希少な核生成イベントを実験的な時間スケールで観察することを不可能にしていた。

手法
著者らは、密度整合溶媒中のステリック安定化PMMA球（直径2.0 µm、多分散性4%）を用いた実空間共焦点顕微鏡法を用い、以下の2つの主要な開発を通じてこの課題に取り組んだ。

1. 高精度な状態点マッピング： 実験とシミュレーションの間の妥当な比較を保証するため、著者らは実験的な体積分率（$\phi$）を有効硬球体積分率（$\phi_{eff}$）へ厳密にマッピングする手法を開発した。

   • 著者らは、ディープラーニングによる粒子追跡アルゴリズム（「Colloidoscope」）を採用し、混雑した環境下での従来のアルゴリズムの限界を克服することで、粒子の局在誤差を直径の約5%に抑えつつ、約98%の粒子の座標を抽出した。
   • 径方向分布関数 $g(r)$ は、この領域における体積分率の変化に対して感度が低いため、これだけに頼るのではなく、高次相関関数を利用した。具体的には、硬球において局所自由エネルギーを最小化するトポロジカル・クラスター分類（TCC）によって特定された、特定の幾何学的クラスターの集団を分析した。
   • 全ペア相関関数 $h(r)$ とTCCクラスター集団の比較に重み付けを行うことで、前例のない精度での状態点マッチング（$\delta\phi < 0.0025$）を実現した。

2. 前臨界核の自由エネルギー測定： 乖離が生じている領域における臨界核は大きすぎて稀少であるため、直接観察することは困難である。そこで著者らは、前臨界核に着目した。これらのより小さなクラスターは、実験的にサンプリング可能な頻度で形成される。

   • 核のサイズ $n$ の関数としての有効形成自由エネルギー $\Delta G(n)/k_B T$ を、観察された核の濃度を用いて $\Delta G(n)/k_B T = -\ln(N_n/N)$ として算出した。
   • 公平な比較を確保するため、シミュレーションデータは、有限の撮像体積、位置誤差（ガウス型変位）、および「見落とされた」粒子を含む実験条件を模倣するように後処理された。
   • 核のサイズは、$q_6$ 配向秩序パラメータ（Ten Woldeらによる）を用い、連結性のためのカットオフ0.7および最小5個の固体的近接隣接数を基準として決定した。

主な結果

• 状態点の合致： 実験的な径方向分布関数および高次クラスター集団は、シミュレーションと極めて良好な一致を示し、実験系が硬球モデルに正確にマッピングされていることが確認された。残留静電相互作用（硬いコア・ユカワ・ポテンシャル）を持つシミュレーションは、実験的なクラスター集団とは不適合であり、硬球近似の妥当性を裏付けた。
• 自由エネルギーの一貫性： 前臨界核から測定された自由エネルギープロファイル（$\Delta G_{obs}$）は、複数の体積分率（$\phi = 0.495$ から $0.520$）において、アンバイアスおよびバイアスのかかったコンピュータシミュレーションと驚くほどよく一致した。
  • $\phi \geq 0.530$ では、約 $k_B T$ の小さな乖離が観察された。著者らは、これは核生成速度における $e$ （約2.7）倍程度の差に相当し、従来の静的光散乱研究で見られた10桁の乖離と比較すれば無視できる程度であると述べている。
• 核の特性： 核の幾何学的解析により、非球形性と回転半径（$R_g$）が実験とシミュレーション間で一致することが明らかになった。大きな核ほど球状に近づく傾向があり、$R_g$ のサイズに対するスケーリングはシミュレーションと一致する冪乗則に従った。

意義と主張
本論文は、核形成の自由エネルギーが、いかなる材料においても実験とシミュレーションの間でこれほど成功的に比較された最初の事例であると主張している。主な意義は、シミュレーションで広く用いられている希少事象サンプリング手法の検証にある。この一致は、これらの計算手法が堅牢かつ正確であることを示唆している。

著者らは、歴史的な乖離は以下の2つの要因に起因すると考えている：

1. 不正確な状態点マッピング： 過去の研究は、体積分率を決定するために間接的な方法（相図からの推論など）に依存しており、それが核生成、相分離、および沈降の競合するタイムスケールの影響を受けていた可能性がある。
2. 間接的な核検出： 静的光散乱は核のサイズと数を間接的に推測するため、核の数を過大評価する可能性がある。

直接的な実空間イメージングと精密な状態点マッピングを用いることで、著者らは、硬球における前臨界核の自由エネルギー障壁が理論的予測と一致していることを示した。著者らは、臨界障壁の高さ（$\Delta G^*$）への外挿については慎重な姿勢を示しているが、前臨界領域における乖離が解消されたことは、硬球核生成の基礎的な物理学が十分に理解されていること、そして以前の大きな乖離は、基礎理論の失敗ではなく、実験手法や状態点の決定におけるアーティファクトであった可能性が高いことを示唆している。