[論文レビュー] Supercooled liquids, the glass transition, and computer simulations

Walter Kob|arXiv (Cornell University)|Dec 14, 2002

Material Dynamics and Properties参考文献 13被引用数 27

ひとこと要約

本稿は、コンピュータシミュレーションを用いて過冷媒とガラス転移をレビューし、動的遅延を説明する枠組みとしてモード結合理論（MCT）に焦点を当てる。MCTは、三体相関（c₃）を含めることで、特に二酸化ケイ素において非エルゴード性パラメータとβ-relaxationダイナミクスを定量的に予測でき、シミュレーションデータと非常に良好な一致を示す。

ABSTRACT

Lecture notes for LES HOUCHES 2002 SUMMER SCHOOL - SESSION LXXVII on SLOW RELAXATIONS AND NONEQUILIBRIUM DYNAMICS IN CONDENSED MATTER; Les Houches: July 1-26, 2002

研究の動機と目的

• 分野に初めて入門する研究者向けに、過冷媒とガラス転移についての簡潔な紹介を提供すること。
• ガラス形成系を研究するにあたり、コンピュータシミュレーションの長所と短所を評価すること。
• モード結合理論（MCT）がガラス形成液体の静的および動的性質を記述するのに有効であるかどうかを検証すること。
• ネットワーク形成ガラス（例：二酸化ケイ素）に対して、MCTが高次相関が含まれる場合に定量的予測を可能にするかを調査すること。
• 特に非エルゴード性パラメータとβ-relaxationダイナミクスに関して、シミュレーション結果と理論的予測を比較すること。

提案手法

• 過冷媒（二酸化ケイ素およびBLJMモデルを含む）の構造的および動的性質を調査するため、分子動力学シミュレーションを用いる。
• 非一様中間散乱関数および非エルゴード性パラメータに注目して、MCTを用いて緩和ダイナミクスを分析する。
• 特にネットワーク形成系に対して理論的予測を改善するために、MCTの頂点に三体相関関数（c₃）を組み込む。
• 波数依存の非エルゴード性パラメータの分析を実施し、さまざまなq値におけるMCT予測の妥当性を検証する。
• β-relaxation関数のマスターカーブの収縮を用いて、MCTのスケーリング予測を検証し、β指数λを抽出する。
• c₃寄与を含む・含まない両方のMCT予測とシミュレーション結果を比較し、定量的精度を評価する。

実験結果

リサーチクエスチョン

• RQ1 モード結合理論（MCT）は、過冷媒における動的遅延をどの程度定量的に記述できるか？
• RQ2 三体相関関数（c₃）は、二酸化ケイ素のようなネットワーク形成ガラスにおけるMCT予測にどのように影響を与えるか？
• RQ3 c₃項を含めることで、非エルゴード性パラメータについてのMCTとシミュレーションデータの一致は改善されるか？
• RQ4 MCTは、二酸化ケイ素における非エルゴード性パラメータの波数依存性を正確に予測できるか？
• RQ5 β-relaxation領域は、ガラス形成液体の動的停止に果たす役割は何か？また、MCTはそのスケーリング行動をどの程度正確に捉えられるか？

主な発見

• MCTに三体相関関数（c₃）を含めることで、二酸化ケイ素における非エルゴード性パラメータについて、理論とシミュレーションの一致が顕著に向上し、特に中程度および大波数領域で顕著である。
• Si-Si相関関数に関しては、c₃項を含むMCTが非常に良好な定量的予測を示し、理論曲線がシミュレーションデータと非常に近接している。
• MCTが予測する臨界指数γ（2.35）は、シミュレーション値（約2.45）と非常に良好に一致しており、高次相関を含めた場合の理論の予測力が裏付けられる。
• 二酸化ケイ素におけるβ-relaxationダイナミクスはマスターカーブの性質を示し、温度が低下するにつれて収縮が顕著になる傾向を示し、MCTの予測と整合的である。
• マスターカーブから抽出したβ指数λは0.713であり、これにより理論的γ値が計算され、MCTスケーリング則との一貫性が確認された。
• シミュレーションでは、ガラス形成液体の粘度がわずかな温度変化の間で12〜14桁も増加しており、ガラス転移の特徴的である劇的な動的遅延が示された。