Structural and dynamic properties of soda-lime-silica in the liquid phase

Abstract

Soda-lime-silica is a glassy system of strong industrial interest. In order to characterize its liquid state properties, we performed molecular dynamics simulations employing an aspherical ion model that includes atomic polarization and deformation effects. They allowed to study the structure and diffusion properties of the system at temperatures ranging from 1400 to 3000 K. We show that Na$^+$ and Ca$^{2+}$ ions adopt a different structural organization within the silica network, with Ca$^{2+}$ ions having a greater affinity for non-bridging oxygens than Na$^+$. We further link this structural behavior to their different diffusivities, suggesting that escaping from the first oxygen coordination shell is the limiting step for the diffusion. Na$^+$ diffuses faster than Ca$^{2+}$ because it is bonded to a smaller number of non-bridging oxygens. The formed ionic bonds are also less strong in the case of Na$^+$.