Kapko Vitaliy, Matyushov Dmitry V, Angell C Austen
Center for Biological Physics, Arizona State University, PO Box 871504, Tempe, Arizona 85287-1504, USA.
J Chem Phys. 2008 Apr 14;128(14):144505. doi: 10.1063/1.2883693.
We report constant-volume and constant-pressure simulations of the thermodynamic and dynamic properties of the low-temperature liquid and crystalline phases of the modified Stillinger-Weber (SW) model. We have found an approximately linear temperature increase of the effective Gaussian width of the distribution of inherent structures. This effect comes from non-Gaussianity of the landscape and is consistent with the predictions of the Gaussian excitations model representing the thermodynamics of the configurational manifold as an ensemble of excitations, each carrying an excitation entropy. The SW model provides us with both the configurational and excess entropies, with the difference mostly attributed to vibrational anharmonicity. We therefore can address the distinction between the excess thermodynamic quantities, often used to interpret experiments, and configurational thermodynamics used to describe the dynamics in the Adam-Gibbs (AG) equation. However we are limited computationally to work at temperatures above the "crossover" temperature at which the breakdown in the Adam-Gibbs relation has been identified in laboratory studies. We find a new break in the slope of the constant pressure AG plot (in the same sense but at much higher temperature than with laboratory data) when the excess entropy is used in the AG equation. This break, which we associate with anharmonic vibrational effects, is not seen when the configurational entropy is used. The simulation diffusivity data are equally well fitted by the AG equation and by a new equation, derived within the Gaussian excitations model, that emphasizes enthalpy over entropy as the thermodynamic control variable for transport in viscous liquids. We show that the modified SW model has close links to the behavior observed for bulk metallic glasses, both in its diffusional and in its thermodynamic properties.
我们报告了对修正的斯蒂林格 - 韦伯(SW)模型的低温液相和晶相的热力学及动力学性质进行的恒容和恒压模拟。我们发现固有结构分布的有效高斯宽度随温度近似呈线性增加。这种效应源于势能面的非高斯性,并且与高斯激发模型的预测一致,该模型将构型流形的热力学表示为激发的集合,每个激发都携带一个激发熵。SW模型为我们提供了构型熵和过剩熵,其差异主要归因于振动非谐性。因此,我们能够探讨常用于解释实验的过剩热力学量与用于描述亚当 - 吉布斯(AG)方程中动力学的构型热力学之间的区别。然而,在计算上我们仅限于在高于“交叉”温度的温度下工作,在实验室研究中已确定在该温度下亚当 - 吉布斯关系会失效。当在AG方程中使用过剩熵时,我们发现恒压AG图的斜率出现了新的断点(与实验室数据的情况相同,但温度要高得多)。当使用构型熵时,未观察到这种与非谐振动效应相关的断点。模拟扩散率数据同样能很好地由AG方程以及在高斯激发模型内推导的一个新方程拟合,该新方程强调焓而非熵作为粘性液体中输运的热力学控制变量。我们表明,修正的SW模型在其扩散和热力学性质方面与大块金属玻璃所观察到的行为有密切联系。