Predicting Composition-Structure Relations in Alkali Borosilicate Glasses Using Statistical Mechanics

Predicting the atomic-scale structure of multicomponent glasses from their composition and thermal history would greatly accelerate the discovery of new engineering and functional glasses. A statistical mechanics-based approach has recently been applied to predict the composition-structure evolution in binary oxide glasses by determining the relative entropic and enthalpic contributions to the bonding preferences. In this work, we first establish the network modifier-former interaction parameters in sodium silicate and sodium borate glasses to predict the structural evolution in sodium borosilicate glasses. Due to the significant variations in the experimentally determined structural speciation in borosilicate glasses, we perform classical molecular dynamics (MD) simulations to establish and validate our structural model. We also show that the statistical mechanical model naturally accounts for the difference in structural speciation from MD simulations and NMR experiments, which in turn arises from the difference in cooling rate and thus thermal history of the glasses. Finally, we demonstrate the predictive capability of the model by accurately accounting for the structural evolution in potassium borosilicate glasses without using any adjustable model parameters. This is possible, because all the interaction parameters are already established in the potassium silicate, potassium borate, and sodium borosilicate glasses, respectively.