[论文解读] 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导出的统计度量，解决传统约束计数在非四面体或不符合八隅体规则原子中的不一致问题。
• 识别在热应力或压力诱导下机械上相关的（完整）与断裂的约束。
• 通过量化随温度和压力变化的完整键角弯曲约束比例，验证玻璃中刚性转变的存在。
• 为复杂氧化物和硫属化物玻璃提供一种物理基础坚实、数据驱动的刚性理论应用方法。

提出的方法

• 使用分子动力学模拟，在不同温度和压力下生成网络玻璃的原子构型。
• 通过随时间计算每个键角（BO）的角度标准偏差（σθ），评估键角弯曲（BB）约束的稳定性。
• 对角度分布函数进行双高斯拟合，以分离完整（σθ较低）与断裂（σθ较高）的约束。
• 基于双峰分布特征，定义一个阈值（σθ ≈ 17°），以区分完整与断裂的BB约束。
• 利用双高斯模型拟合结果，追踪完整BB约束比例（q(T)）随温度的变化。
• 将结构特征（如第一衍射峰的半高宽（FWHM））与刚性转变窗口的存在相关联。

实验结果

研究问题

• RQ1 能否利用分子动力学模拟识别网络玻璃中哪些键角弯曲约束在机械上保持完整或已断裂？
• RQ2 完整键角弯曲约束的比例如何随温度演化，是否与刚性转变相关？
• RQ3 角度偏差（σθ）的双峰分布能否用于定义完整与断裂约束之间的阈值？
• RQ4 是否存在一个密度窗口，其相关长度（L）最大，是否对应于增强的刚性或结构有序性？
• RQ5 压力与成分如何共同影响网络玻璃的刚性，这一影响能否通过约束统计揭示？

主要发现

• 在300 K时，键角的角度分布集中在低σθ区域，表明所有键角弯曲约束均保持完整。
• 在4500 K时，角度分布显著展宽，表明所有键角弯曲约束均因热激活而断裂。
• 在玻璃化转变温度（Tg）附近，出现双峰分布，其中一个峰对应完整约束，另一个对应断裂约束。
• 完整键角弯曲约束比例（q(T)）随温度升高而降低，且数据与Mauro-Gupta公式拟合良好。
• 基于双峰分布的最小值，确定σθ ≈ 17°为完整与断裂约束的边界阈值。
• 在密度ρ = 2.1至3.5 g/cm³的窗口内，相关长度（L = 2π / FWHM）达到最大值，同时黏度最小，网络形成原子扩散最大，表明存在结构转变。