Abstract:Molecular dynamics simulations are performed to study spatially heterogeneous dynamics in a model of viscous silica above and below the critical temperature of the mode coupling theory, T(MCT) . Specifically, we follow the evolution of the dynamic heterogeneity as the temperature dependence of the transport coefficients shows a crossover from non-Arrhenius to Arrhenius behavior when the melt is cooled. It is demonstrated that, on intermediate time scales, a small fraction of oxygen and silicon atoms are more mobile than expected from a Gaussian approximation. These highly mobile particles form transient clusters larger than that resulting from random statistics, indicating that the dynamics are spatially heterogeneous. An analysis of the clusters reveals that the mean cluster size is maximum at times intermediate between ballistic and diffusive motion, and the maximum size increases with decreasing temperature. In particular, the growth of the clusters continues into the crossover to Arrhenius behavior for the transport coefficients. These findings imply that the structural relaxation in silica cannot be understood as a statistical bond breaking process. Although the mean cluster sizes for silica are at the lower end of the spectrum of values reported in the literature for other model liquids, we find that spatially heterogeneous dynamics in models of strong and fragile glass formers are similar on a qualitative level. However, unlike the results for fragile liquids, we show that correlated particle motion along quasi-one-dimensional, stringlike paths is of little importance for the structural relaxation in this model of silica, suggesting that stringlike motion is suppressed by the presence of a network structure. To study transient clusters of localized particles, we calculate a generalized susceptibility corresponding to the self-part of a four-point time dependent density correlation function. We find that this generalized susceptibility is maximum on the time scale of the structural relaxation, where a strong increase of the peak height indicates a growing length of spatial correlations between localized particles upon cooling. Characterizing the local structural environments of the most mobile and the most immobile particles, respectively, we show that high particle mobility is facilitated by, but not limited to, the vicinity of defects in the network structure.

Microscopic description of the low-temperature anomalies in silica and lithium silicate via computer simulations

Tunneling dynamics of side chains and defects in proteins, polymer glasses, and OH-doped network glasses

Characterization of the complex ion dynamics in lithium silicate glasses via computer simulations

Multilevel Tunnelling Systems and Fractal Clustering in the Low-Temperature Mixed Alkali-Silicate Glasses

Glassy dynamics and crystalline local order in two-dimensional amorphous silica

Decay of mean-values of multiplicative functions

Dynamics of Water in LiCl and CaCl2 Aqueous Solutions Confined in Silica Matrices: A Backscattering Neutron Spectroscopy Study

Molecular dynamics simulation of silica liquid and glass

Graph theory-based structural analysis on density anomaly of silica glass

The Nature of the Low-Temperature Crossover of Water in Hard Confinement

The Dynamics of Non-Crystalline Silica: Insight from Molecular Dynamics Computer Simulations

Molecular structural order and anomalies in liquid silica

Water-like anomalies as a function of tetrahedrality

Anomalous properties of the acoustic excitations in glasses on the mesoscopic length-scale

Temperature dependence of spatially heterogeneous dynamics in a model of viscous silica

Identification of structural motifs as tunneling two-level systems in amorphous alumina at low temperatures

Topological Data Analysis for Revealing the Structural Origin of Density Anomalies in Silica Glass

Depletion of two-level systems in ultrastable computer-generated glasses

Complex dynamics of partially freezable confined water revealed by combined experimental and computational studies

Influence of Random Internal Fields on the Tunneling of OH Defects in NaCl Crystals

Search for a liquid-liquid critical point in models of silica