[論文レビュー] Reconsidering the structure of nucleation theories

Anja Kuhnhold, Hugues Meyer+3 more|arXiv (Cornell University)|Sep 26, 2018

nanoparticles nucleation surface interactions参考文献 2被引用数 5

ひとこと要約

本稿では、非平衡射影作用素形式を用いて、核化の正確な非マルコフ型運動方程式を導出しており、古典的核化理論（CNT）がその極限的状態であることを示している。Lennard-Jones液体の結晶化シミュレーションから、時間に依存する記憶核を持つ顕著な記憶効果が観測され、マルコフ型仮定を破る。これにより、CNTで用いられる標準的なフォッカー・プランク近似は無効であることが判明した。

ABSTRACT

We discuss the structure of the equation of motion that governs nucleation processes at first order phase transitions. From the underlying microscopic dynamics of a nucleating system, we derive by means of a non-equilibrium projection operator formalism the equation of motion for the size distribution of the nuclei. The equation is exact, i.e. the derivation does not contain approximations. To assess the impact of memory, we express the equation of motion in a form that allows for direct comparison to the Markovian limit. As a numerical test, we have simulated crystal nucleation from a supersaturated melt of particles interacting via a Lennard-Jones potential. The simulation data show effects of non-Markovian dynamics.

研究の動機と目的

• 微視的非平衡ダイナミクスから近似を一切用いずに、核のサイズ分布に対する正確な運動方程式を導出すること。
• 古典的核化理論（CNT）の中心的役割を果たすマルコフ型近似の妥当性を、ダイナミクスにおける記憶効果の分析を通じて評価すること。
• CNTで用いられる標準的なフォッカー・プランク方程式が、結晶核化におけるクラスターサイズ分布を正確に記述できるかを検証すること。
• Lennard-Jones液体の分子動力学的シミュレーションを用いて、非マルコフ型記憶の核化速度論への影響を定量化すること。
• 記憶効果の時間スケールが誘導時間と同等であるかどうかを特定し、物理的関連性を評価すること。

提案手法

• 微視的ダイナミクスから核のサイズ分布に対する正確な運動方程式を、非平衡射影作用素形式を用いて導出する。
• 時間に非局所的なKramers-Moyal展開形式で方程式を表現し、記憶効果を明示的に捉える。
• 記憶核の時間依存性を分析することで、方程式のマルコフ型極限と比較し、特にフォッカー・プランク方程式におけるデルタ関数近似の妥当性を検証する。
• 過飽和Lennard-Jones液体の分子動力学的シミュレーションを実施し、クラスターサイズ分布データを生成する。
• 記憶核および拡散関数のパラメータ化されたアンサンブルを用いて、シミュレーションデータをフォッカー・プランク方程式と完全な非マルコフ型方程式の両方にフィットさせる。
• シミュレーションデータとモデル予測の誤差を最小化するためにシミュレーテッド・アニーリングを用い、適合度の評価には正規化誤差指標Sを用いる。

実験結果

リサーチクエスチョン

• RQ1 記憶効果は、古典的核化理論におけるマルコフ型仮定をどの程度無効にするか？
• RQ2 標準的なフォッカー・プランク方程式は、結晶核化の分子動力学的シミュレーションで観測されたクラスターサイズ分布を正確に再現できるか？
• RQ3 核化における記憶効果の特徴的な時間スケールは何か？また、臨界核への最初の通過時間と比較するとどうなるか？
• RQ4 非マルコフ型記憶核は、フォッカー・プランク方程式で仮定されているデルタ関数形から著しく逸脱しているか？
• RQ5 非マルコフ型方程式を用いることで得られる適合度の向上は、統計的にも物理的にも意味のあるものか？

主な発見

• 導出された運動方程式は正確かつ非マルコフ型であり、記憶核が明示的に時間依存しており、デルタ関数分布ではないため、フォッカー・プランク方程式の仮定に反する。
• フォッカー・プランクモデルでは正規化誤差S = 0.02987となり、クラスターサイズ分布のシミュレーションデータとの一致が著しく悪いことが示された。
• 非マルコフ型モデルではS = 0.00975という著しく良好な適合度を達成し、記憶効果が正確な記述に不可欠であることが確認された。
• 主な駆動項の記憶核$ ilde{d}_0(n,t)$は時間に対してデルタ関数ではないため、マルコフ型極限が満たされていないことが確認された。
• 記憶効果の時間スケール（τ₀ ≈ 14–14.2、n ≈ 30）は、平均最初通過時間（t ≈ 80）と同程度であり、核化過程において物理的に関連性があることが示された。
• 非マルコフ型ダイナミクスは、核化障壁付近（n ≈ 33）で特に重要であり、モデル間の差異が最大となる。