Structural organisation in oxide glasses from molecular dynamics modelling

Classical molecular dynamics modelling has been used to obtain new models of 50CaO·50P₂O₅ and 50MgO·50SiO₂ glasses and, together with previously published models of 63CaO·37Al₂O₃, and 50CaO·50SiO₂ glasses, these have been inspected to evaluate structural features. For the first time, models of glasses near the eutectic in three systems, aluminate, silicate, and phosphate, with the same modifier, Ca, have been compared. All have short range order which is similar to that in crystals of the same composition, 5CaO·3Al₂O₃, CaSiO₃ and Ca(PO₃)₂. There is a clear trend in bonding of bridging oxygen to Ca, which is dominant in aluminate glass, common in silicate glass, and absent in phosphate glass. Preliminary results for 50MgO·50SiO₂ glass show unusual behaviour because ~ 5% of oxygen is present as “non-network” oxygen, i.e. bonded only to Mg. The models show broader Qⁿ distributions than seen in NMR experiments, and this remains an area for improvement of MD modelling of glasses. The distributions of Ca in the models have been studied using the pair distribution function T_CaCa(r) which is found to be similar in the three glasses, and also similar to the previous experimental measurement for 50CaO·50SiO₂ glass. The distributions of Ca are markedly different in the glasses compared to the crystals, being isotropic in the former and anisotropic in the latter, which should be a factor in glass forming ability.