Spatially heterogeneous dynamics and dynamic facilitation in a model of viscous silica

We perform molecular dynamics simulations to study the structural relaxation dynamics of a model of viscous silica, the prototype of a strong glass former. We find that the melt dynamics are spatially heterogeneous regardless of whether the bulk relaxation is non-Arrhenius or Arrhenius, and cannot be understood as a statistical bondbreaking process. Further, we show that stringlike motion is suppressed by the covalent bondings, yet high particle mobility propagates continuously, supporting the concept of dynamic facilitation emphasized in recent theoretical work.