[論文レビュー] Topological Constraints and Rigidity of Network Glasses from Molecular Dynamics Simulations

Mathieu Bauchy|arXiv (Cornell University)|Jun 22, 2015

Complex Network Analysis Techniques被引用数 4

ひとこと要約

本論文は、原子の運動の統計的分析を通じて、原子スケールの結合伸長および結合角変形制約を定量化することにより、分子動力学（MD）シミュレーションにトポロジカル制約理論を厳密に適用する手法を提案する。角標準偏差を用いることで、完全な（低い角偏差）および破壊された（高い角偏差）制約を区別でき、正確な剛性予測が可能となり、理論が現実のガラス形成系への適用可能性が検証される。

ABSTRACT

Due to its non-crystalline nature, the glassy state has remained one the most exciting scientific challenges. To study such materials, Molecular Dynamics (MD) simulations have been extensively used because they provide a direct view into its microscopic structure. MD is therefore used not only to reproduce real system properties but also benefits from detailed atomic scale analysis. Unfortunately, MD shows inherent limitations because of the limited computational power. For instance, only the simulations of small systems are currently permitted, which prevents from studying small compositional changes, although it is well known that they can dramatically alter system properties. At this stage, it is tempting to follow topological constraint theory, which aims at describing macroscopic properties of the glass relying only on the connectivity of individual atoms, thus considering the complicated glass network as simple mechanical trusses. Thanks to only basic hand calculations, this theory has been successful in predicting complex composition and temperature behavior of glass properties, such as the glass forming ability, the viscosity of the corresponding liquid or the elasticity. The purpose of my PhD work is to connect these successful approaches based only on the topology of the underlying low-temperature network with properties that can be obtained from MD calculations. It should thus allow for an increased applicability of rigidity theory.

研究の動機と目的

• 制約数え上げを原子スケールの動的挙動に基づくものにすることで、トポロジカル制約理論と分子動力学シミュレーションを結びつける。
• 非四面体的または八電子則を満たさない原子に対して従来の制約数え上げに生じる不一致を、MDから得られる統計的指標を用いて解消する。
• 熱的または圧力誘発の応力下で、どの制約が機械的に関連している（完全な状態）か、破壊されたかを特定する。
• 温度および圧力関数として完全な結合角変形制約の割合を定量化することにより、ガラスにおける剛性転移の存在を検証する。
• 物理的根拠に基づいた、データ駆動型の方法を提供し、複雑な酸化物および chalcogenide ガラスへの剛性理論の適用を可能にする。

提案手法

• さまざまな温度および圧力条件下でのネットワークガラスの原子配置を分子動力学シミュレーションで生成する。
• 各結合角（BO）の時間的角標準偏差（σθ）を計算し、結合角変形（BB）制約の安定性を評価する。
• 二重ガウスフィッティングを角分布関数に適用し、完全な（σθが小さい）および破壊された（σθが大きい）制約を分離する。
• 二峰性分布の特徴に基づき、σθ ≈ 17° の閾値を、完全なと破壊されたBB制約を区別する基準として定義する。
• 二重ガウスモデルへのフィッティングを用いて、温度関数としての完全なBB制約の割合（q(T)）を追跡する。
• 第一回折りたたみピークの半値全幅（FWHM）といった構造的シグネチャと、剛性転移ウィンドウの存在を相関付ける。

実験結果

リサーチクエスチョン

• RQ1 分子動力学シミュレーションを用いて、ネットワークガラスにおいてどの結合角変形制約が機械的に完全か、破壊されたかを特定できるか？
• RQ2 完全な結合角変形制約の割合は温度に応じてどのように変化し、剛性転移と相関するか？
• RQ3 角偏差（σθ）の二峰性分布を用いて、完全な制約と破壊された制約の境界を定義できるか？
• RQ4 相関長（L）が最大となる密度ウィンドウが、剛性の増大または構造的秩序の向上を示す領域に対応するか？
• RQ5 圧力と組成が、制約統計から明らかになるネットワークガラスの剛性にどのように同時に影響するか？

主な発見

• 300 Kでは、結合角の角分布が低σθで鋭くピークを示しており、すべての結合角変形制約が完全であることが示唆される。
• 4500 Kでは、角分布が顕著に広がっており、熱的励起によってすべての結合角変形制約が破壊されていることが示唆される。
• ガラス転移温度（Tg）付近では、二峰性分布が出現し、一方のピークが完全な制約、もう一方が破壊された制約に対応する。
• 完全な結合角変形制約の割合（q(T)）は温度上昇に伴い減少し、マウロ＝ギプタ式によるフィッティングが良好に一致する。
• 二峰性分布の最小点に基づき、σθ ≈ 17° が完全な制約と破壊された制約の境界として特定された。
• 密度ρ = 2.1 から 3.5 g/cm³ の間のウィンドウでは、相関長（L = 2π / FWHM）が最大となり、粘度が最小で、ネットワーク形成原子の拡散度が最大となることと一致しており、構造的転移を示唆する。