D. E. Shaw Research , New York, New York 10036, United States.
J Phys Chem B. 2013 Oct 24;117(42):12898-907. doi: 10.1021/jp402102w. Epub 2013 Jul 10.
Understanding the nature of the glass transition--the dramatic slowing of dynamics and eventual emergence of a disordered solid from a cooling liquid--is a fundamental challenge in physical science. A central characteristic of glass-forming liquids is a non-exponential main relaxation process. The extent of deviation from exponential relaxation typically becomes more pronounced on cooling. Theories that predict a growth of spatially heterogeneous dynamics as temperature is lowered can explain these observations. In apparent contradiction to these theories, however, some experiments suggest that certain substances--notably including the intensely studied molecular glass-former ortho-terphenyl (OTP)--have a main relaxation process whose shape is essentially temperature independent, even though other observables predicted to be correlated with the degree of dynamical heterogeneity are temperature dependent. Here we report the first simulations based on an atomistic model of OTP that reach equilibrium at temperatures well into the supercooled regime. We first show that the results of these simulations are in reasonable quantitative agreement with experimental data for several basic properties over a wide range of temperatures. We then focus on rotational relaxation, finding nearly exponential behavior at high temperatures with clearly increasing deviations as temperature is lowered. The much weaker temperature dependence observed in light-scattering experiments also emerges from the same simulation data when we calculate correlation functions similar to those probed experimentally; this highlights the diversity of temperature dependencies that can be obtained with different probes. Further analysis suggests that the temperature insensitivity observed in the light-scattering experiments stems from the dependence of these measurements on internal as well as rotational molecular motion. Within the temperature range of our OTP simulations, our results strongly suggest that this archetypal glass-former behaves as anticipated by theories of the glass transition that predict increasing non-exponentiality with cooling, and our simulations thus strengthen the evidence supporting such theories.
理解玻璃化转变的本质——即液体冷却时动力学急剧减缓并最终出现无序固体——是物理科学的一个基本挑战。玻璃形成液体的一个核心特征是主要弛豫过程是非指数的。非指数弛豫的偏离程度通常随着冷却而变得更加明显。可以预测温度降低时空间异质动力学增长的理论可以解释这些观察结果。然而,与这些理论明显矛盾的是,一些实验表明,某些物质——特别是包括被深入研究的分子玻璃形成剂邻位三联苯(OTP)——的主要弛豫过程的形状基本上与温度无关,尽管其他被预测与动力学异质性程度相关的可观察量是与温度有关的。在这里,我们报告了第一个基于 OTP 原子模型的模拟,这些模拟在远低于过冷温度的温度下达到平衡。我们首先表明,这些模拟的结果在相当大的温度范围内与实验数据在几个基本性质上的定量一致性是合理的。然后,我们专注于旋转松弛,发现高温下几乎呈指数行为,随着温度降低,偏离度明显增加。在光散射实验中观察到的弱得多的温度依赖性也从相同的模拟数据中出现,当我们计算与实验中探测到的类似的相关函数时;这突出了不同探针可以获得的温度依赖性的多样性。进一步的分析表明,在光散射实验中观察到的温度不敏感性源于这些测量对内旋转和分子运动的依赖性。在我们的 OTP 模拟的温度范围内,我们的结果强烈表明,这种典型的玻璃形成剂的行为与预测冷却时非指数性增加的玻璃转变理论一致,我们的模拟因此加强了支持这些理论的证据。