[論文レビュー] Local vs. Cooperative: Unraveling Glass Transition Mechanisms with SEER

Massimo Pica Ciamarra, Wencheng Ji+1 more|arXiv (Cornell University)|Feb 10, 2023

Material Dynamics and Properties被引用数 4

ひとこと要約

本稿では、過冷却液体の構造から局所的励起状態とそのエネルギー障壁を体系的に抽出する新規アルゴリズムSEERを提案する。冷却に伴いエネルギー障壁が上昇する局所的エネルギー障壁が、ミリ秒スケールまで緩和ダイナミクスを定量的に制御していることが判明し、ガラス転移における協同的理論に挑戦するものであり、局所的障壁が脆さの支配的要因である可能性を示唆する。

ABSTRACT

Which phenomenon slows down the dynamics in super-cooled liquids and turns them into glasses is a long-standing question of condensed-matter. Most popular theories posit that as the temperature decreases, many events must occur in a coordinated fashion on a growing length scale for relaxation to occur. Instead, other approaches consider that local barriers associated with the elementary rearrangement of a few particles or `excitations' govern the dynamics. To resolve this conundrum, our central result is to introduce an algorithm, SEER, which can systematically extract hundreds of excitations and their energy from any given configuration. We also provide a novel measurement of the activation energy, characterizing the liquid dynamics, based on fast quenching and reheating. We use these two methods in a popular liquid model of polydisperse particles. Such polydisperse models are known to capture the hallmarks of the glass transition and can be equilibrated efficiently up to millisecond time scales. The analysis reveals that cooperative effects do not control the fragility of such liquids: the change of energy of local barriers determines the change of activation energy. More generally, these methods can now be used to measure the degree of cooperativity of any liquid model.

研究の動機と目的

• ガラス転移付近の遅れを支配するのは協同的ダイナミクスか、局所的エネルギー障壁かという長年の議論を解決すること。
• 平衡状態からの過冷却液体における局所的励起状態の同定と定量的評価の体系的手法を開発すること。
• 局所的励起状態のエネルギー障壁の温度依存性のみで、ミリ秒スケールまで観測された緩和ダイナミクスを定量的に説明できるかを検証すること。
• 協同的理論（例：Adam-Gibbs理論、KCMs）を、局所的障壁モデルの予測とシミュレーションデータを定量的に比較することで挑戦すること。

提案手法

• SEERは、与えられた構造から最小エネルギー経路を探索するためのナーヴド・エラスティックバンド法を用いて、局所的励起状態を体系的に同定する。
• 各同定された励起状態について、そのエネルギー障壁を、励起状態を抑制する最小の剛性パラメータ $ k_e $ を特定することで計算する。これにより、安定性と収束性が保証される。
• アルゴリズムは、異なる待ち時間 $ t_w $ を用いた熱的サイクルと、50回の独立ランのアンサンブル平均を用いて、励起状態カタログの統計的安定性と収束性を確保する。
• 重複する励起状態は距離閾値 $ d_{ij} $ を適用することで除去され、最終的なカタログには、明確に分離された独立した励起状態のみが保持される。
• 固有構造のエネルギーシフト $ riangle E $ は、平衡状態を低温にクエンチし、初期状態と最終状態のエネルギー差を解析することで推定される。
• 時間-温度スーパー・ポジションを適用して相関関数を重ね合わせ、$ eta acksimeq 0.52 $ のストレッチド指数的減衰が確認され、低温における緩和時間の推定が可能になった。

Figure 1: a Temperature dependence of the relaxation time $\tau$ estimated from the decay of the overlap function in our poly-disperse (see Sec. III ). It is customary to define the onset temperature $T_{\rm on}$ as the temperature below which the energy of inherent structures starts to decrease. $T

実験結果

リサーチクエスチョン

• RQ1 ガラス転移付近のダイナミクスの遅れは、増大する協同的長尺スケールではなく、局所的エネルギー障壁によって支配されているか？
• RQ2 局所的励起状態のエネルギー障壁の温度依存性のみで、平衡状態の過冷却液体におけるミリ秒スケールまでの緩和時間の予測が定量的に可能か？
• RQ3 異なる熱的歴史や系サイズにおいて、SEERアルゴリズムが一貫して異なる励起状態を抽出できるか、その頑健性はいかがなものか？
• RQ4 局所的障壁がダイナミクスにおいて支配的であるにもかかわらず、協同的効果が観察される脆さにどの程度寄与しているか？

主な発見

• 冷却に伴い励起状態の密度状態が高エネルギー側にシフトし、エネルギー障壁の上昇が緩和ダイナミクスの遅れと直接的に関連していることが明らかになった。
• 局所的障壁エネルギーのシフトに基づく予測された緩和時間は、ミリ秒スケールまでシミュレーションデータと定量的に一致した。
• 発見された励起状態の数は $ N_E $ の増加とともに収束し、待ち時間 $ t_w $ に依存しないため、アルゴリズムの安定性と頑健性が確認された。
• 高温で異常な高さの $ riangle E $ を示す外れ値の構造は、即時のクエンチを避けることで除去され、サンプル間の一貫性が向上した。
• 単一の最大値を持つ基本的励起状態を用いることで、偶然的な組み合わせに起因する誤差の多い複雑な励起状態よりも信頼性の高いギャップエネルギー推定が可能になった。
• 時間-温度スーパー・ポジションにより、重ね合わせ関数のストレッチド指数的減衰（$ eta acksimeq 0.52 $）が確認され、低温における緩和時間の信頼性ある推定が可能になった。

Figure 2: Schematic illustration of SEER’s working principle. a The algorithm suppresses a discovered excitation by transforming it into a saddle bifurcation; b in real space, it corresponds to adding restoring forces opposite to displacements along the suppressed excitations.