[論文レビュー] Crystallization of Simple Fluids: Relative Stability of f.c.c. and b.c.c Structures

Swarn Lata Singh, Yashwant Singh|arXiv (Cornell University)|Feb 27, 2009

Material Dynamics and Properties被引用数 11

ひとこと要約

本稿では、オーロンシュタイン＝ツェルニケ方程式とロジャーズ＝ヤング閉じ込め方程式を用い、対称性保存部と対称性破れ部の両方を含む直接対相関関数（DPCF）の自由エネルギー関数を構築する。対称性破れ部は摂動展開で扱い、逆べき乗系においてn=12ではf.c.c.に、よりソフトなポテンシャル（n≤6）ではb.c.c.に凍結するという予測を正確に得られ、密度汎関数理論における長年の課題を解決する。

ABSTRACT

A free-energy functional for a crystal that contains both the symmetry conserved and symmetry broken parts of the direct pair correlation function is developed. The free-energy functional is used to investigate the crystallization of fluids interacting via the inverse power potential ; $u(r)=ε{(σ/r)}^n$. In agreement with simulation results we find that for $n=12$ the freezing is into close packed f.c.c structure while for soft repulsions $(n\leq 6)$ b.c.c phase is more stable.

研究の動機と目的

• 流体の結晶化におけるf.c.c.とb.c.c.相の相対的安定性を正しく捉える第一原理的統計力学理論を構築すること。
• 従来の密度汎関数理論が、直接対相関関数（DPCF）の対称性破れ寄与を無視しているため、流体-固体転移を正しく記述できないという問題を解決すること。
• 結晶相に現れる、性質的に新しい対称性破れ成分をDPCFに組み込むこと。この成分は融解時に消える。
• 逆べき乗ポテンシャル（n=12, 6, 4）のさまざまなソフトネスにおいて、モンテカルロシミュレーションデータと比較して、この改善された関数の予測能力を検証すること。

提案手法

• 自由エネルギー関数 A[ρ] = A_id[ρ] + A_ex[ρ] を定式化し、A_exにDPCFの対称性保存部 c^(0) と対称性破れ部 c^(b) の両方を含める。
• c^(0)(r) を、熱力学的整合性を保つために、ロジャーズ＝ヤング閉じ込め方程式を用いたオーロンシュタイン＝ツェルニケ方程式で計算する。
• c^(b)(r₁,r₂;[ρ]) を汎関数テイラー展開によりモデル化し、一次項は等方的流体の三体DPCFと密度ゆらぎ ρ_b(r) の積分として表現する。
• c^(b) を逆格子ベクトル G を用いたフーリエ空間で表現し、係数 c^(G)(r) を差分ベクトル r と結晶構造に依存する形で定式化する。
• 密度および秩序パラメータ場の関数的積分を、二重パラメータ化（λ と ξ）を用いて実行し、結晶の過剰自由エネルギー A_ex^(b) を計算する。
• 流体-固体共存点を特定するため、巨視的熱力学的関数 ΔW = W_l - W を用い、転移点では ΔW/N = 0 となる。

実験結果

リサーチクエスチョン

• RQ1 対称性保存部と対称性破れ部の両方を含む自由エネルギー関数が、逆べき乗系におけるf.c.c.とb.c.c.相の相対的安定性を正しく予測できるか？
• RQ2 特にソフトな反発ポテンシャルにおいて、対称性破れDPCFが自由エネルギー差に与える寄与はどの程度顕著か？
• RQ3 等方的流体の三体DPCFを用いたc^(b)の摂動的取り扱いが、固体相の安定性を信頼できる見積もりで与えられるか？
• RQ4 なぜ従来のDFT手法は、ソフト反発ポテンシャルの結晶化経路を正しく予測できないのか？
• RQ5 c^(b) を含めることで、n=12, 6, 4 におけるモンテカルロシミュレーション結果との一致がどの程度向上するか？

主な発見

• n=12 では、f.c.c. 構造がより安定であると予測され、モンテカルロシミュレーションと一致する。これは、対称性保存部と対称性破れ部の両方が有利に寄与しているためである。
• n=6 および n=4 では、b.c.c. 構造が f.c.c. よりも安定となり、シミュレーション結果を正しく再現する。対称性破れDPCFが安定化に重要な役割を果たしている。
• 自由エネルギー差 ΔW/N における対称性破れ寄与は無視できない。n=12 では対称性保存部の約25%、n=4 ではほぼ50%を占める。
• 対称性破れ項を含めない場合、特にソフトなポテンシャルでは流体の安定性が過剰に予測され、誤った予測につながる。
• 等方的流体の三体DPCFを用いたc^(b)の摂動的表現は、よりソフトなポテンシャル（n=6, 4）ではその大きさを過大評価し、シミュレーションよりも低いγ値で転移が予測される。
• f.c.c. と b.c.c. 構造における c^(G)(r) の関数形は、大きさとr依存性において顕著に異なるため、両者の安定性の違いが説明できる。