[论文解读] Structure, diffusion and orientational freezing in lithium metasilicate

Cristian Balbuena, Carolina Brito+1 more|arXiv (Cornell University)|Aug 7, 2013

Glass properties and applications参考文献 21被引用 3

一句话总结

本研究利用分子动力学模拟研究了锂 metasilicate (Li₂SiO₃) 中的动态与结构转变，识别出两个关键温度：Tc ≈ 1500 K，与模式耦合理论相关，且相对四面体运动停止；Tg ≈ 1000 K，标志着玻璃化转变，此时四面体内部原子位置停止松弛。关键发现是取向序参数——用于衡量相对四面体取向与局部键序——在 Tc 和 Tg 处表现出显著变化，将动态冻结与超出标准模式耦合理论的结构重排联系起来。

ABSTRACT

We report on the dynamic and structural characterization of lithium metasilicate $Li_2SiO_3$, a network forming ionic glass, by means of molecular dynamics simulations. The system is characterized by a network of $SiO_4$ tetrahedra disrupted by $Li$ ions which diffuse through the network. Measures of mean square displacement of $Si$ and $O$ atoms allow us to identify a temperature at which tetrahedra stop moving relative to each other. This temperature $T_c\approx 1500\,K$ can be characterized within the framework of mode coupling theory. At a much lower temperature $T_g\approx 1000\,K$, a change in the slope of the volume versus temperature data allows to single out the glass transition. We find signatures of both transitions in structural order parameters, related to the orientation of tetrahedra. Going down in temperature we find that, around the mode coupling transition temperature, a set of order parameters which measure the relative orientation of tetrahedra cease to increase and stay constant below $T_c$. Another well known measure of orientational order, the bond orientational order parameter, which in the studied system measures local order within single tetrahedrons, is found to continue growing below $T_c$ until $T_g$, below which it remains constant. Our results allow to relate two characteristic dynamic transitions with corresponding structural transitions, as observed in two different orientational order parameters. Furthermore, the results indicate that the network of thetrahedra continue to relax well below the point where neighboring tetrahedra cannot rearrange relative to each other, and the glass is reached only upon a process of relaxation of atoms which form the thetrahedron, as quantified by the change in the bond orientational order parameters.

研究动机与目标

• 理解网络形成离子玻璃 Li₂SiO₃ 中的动态与结构转变。
• 研究在玻璃化转变附近动态减缓时是否存在结构特征。
• 考察四面体网络有序性与取向相关性在锂 metasilicate 动力学中的作用。
• 检验模式耦合理论（MCT）在具有强短程有序性的体系（如硅酸盐网络）中的适用性。
• 利用超越局部键序的取向序参数，探索动态冻结与结构重排之间的关联。

提出的方法

• 对 Li₂SiO₃ 进行分子动力学模拟，以研究温度变化过程中的原子运动与结构演化。
• 通过分析均方位移与自中间散射函数，提取扩散系数与弛豫时间，从而识别模式耦合临界温度 Tc。
• 通过体积随温度曲线的斜率变化确定玻璃化转变温度 Tg。
• 采用两个取向序参数：用于评估局部四面体有序性的 Steinhardt 键取向序参数，以及用于量化相邻四面体之间相对取向的 Rey 参数。
• 跟踪这些序参数随温度的变化，以检测与 Tc 和 Tg 对应的结构转变。
• 在模式耦合理论框架内解释结果，并与热力学与动力学可观测量进行比较。

实验结果

研究问题

• RQ1 锂 metasilicate 中模式耦合转变 Tc ≈ 1500 K 是否对应四面体网络中可检测到的结构变化？
• RQ2 在 Tg ≈ 1000 K 的玻璃化转变处，除了原子运动的动态冻结外，是否存在其他结构特征？
• RQ3 取向序参数——特别是衡量相对四面体取向与局部键序的参数——在 Tc 和 Tg 附近温度变化时如何响应？
• RQ4 在 Tc 以下，四面体转动与内部弛豫的动力学在多大程度上发生解耦？
• RQ5 系统中观察到的两步弛豫过程是否可通过取向序参数与不同的结构转变相联系？

主要发现

• 从扩散与弛豫时间的幂律标度中识别出模式耦合临界温度 Tc ≈ 1500 K，与 MCT 预测一致。
• 在 Tc 处，通过 Rey 参数测量的相邻 SiO₄ 四面体之间的相对取向不再增加并保持恒定，表明四面体重取向发生结构冻结。
• 玻璃化转变温度 Tg ≈ 1000 K 由体积随温度曲线的斜率变化所标记，对应于四面体内部原子松弛的停止。
• 键取向序参数在 Tc 以下继续增长，仅在 Tg 处趋于饱和，表明四面体原子的内部松弛在四面体间相对运动停止后仍持续进行。
• Tc 与 Tg 之间的温度窗口较大，表明在四面体旋转动态冻结开始后，系统仍经历了显著的结构弛豫。
• 本研究建立了动态转变（Tc 与 Tg）与取向序中明显结构变化之间的直接联系，证明两种转变在四面体网络中均有清晰的结构特征。