Viscosity, breakdown of Stokes-Einstein relation, and dynamical heterogeneity in supercooled liquid Ge2Sb2Te5 from simulations with a neural network potential

Phase change materials are exploited in non-volatile electronic memories and photonic devices that rely on a fast and reversible transformation between the amorphous and crystalline phases upon heating. Recrystallization of the amorphous phase under the operation conditions of the memories occurs in the supercooled liquid phase above the glass transition temperature Tg. The dynamics of the supercooled liquid is thus of great relevance for the operation of the devices and, close to Tg, also for the structural relaxations of the glass that affect the performances of the memories. Information on the atomic dynamics is provided by the diffusion coefficient (D) and by the viscosity (η), which are, however, both difficult to be measured experimentally under the operation conditions of the devices due to fast crystallization. In this work, we leverage a machine learning interatomic potential for the flagship phase change compound Ge2Sb2Te5 to compute η, D, and the α-relaxation time in a wide temperature range from 1200 K to about 100 K above Tg. Large-scale molecular dynamics simulations allowed quantifying the fragility of the liquid and the occurrence of a breakdown of the Stokes-Einstein relation between η and D in the supercooled phase. Isoconfigurational analysis provided a visualization of the emergence of dynamical heterogeneities responsible for the breakdown of the Stokes-Einstein relation. The analysis revealed that the regions of most mobile atoms are related to the presence of Ge atoms with particular local environments.