Cooling rate and size effects on the medium-range structure of multicomponent oxide glasses simulated by molecular dynamics

Abstract

A set of molecular dynamics simulations were performed to investigate the effect of cooling rate and system size on the medium-range structure of melt-derived multicomponent silicate glasses, represented by the quaternary 45S5 Bioglass composition. Given the significant impact of the glass degradation on applications of these materials in biomedicine and nuclear waste disposal, bulk structural features which directly affect the glass dissolution process are of particular interest. Connectivity of the silicate matrix, ion clustering and nanosegregation, distribution of ring and chain structural patterns represent critical features in this context, which can be directly extracted from the models. A key issue is represented by the effect of the computational approach on the corresponding glass models, especially in light of recent indications questioning the suitability of conventional MD approaches (that is, involving melt-and-quench of systems containing ~10(3) atoms at cooling rates of 5-10 K/ps) when applied to model these glasses. The analysis presented here compares MD models obtained with conventional and nonconventional cooling rates and system sizes, highlighting the trend and range of convergence of specific structural features in the medium range. The present results show that time-consuming computational approaches involving much lower cooling rates and/or significantly larger system sizes are in most cases not necessary in order to obtain a reliable description of the medium-range structure of multicomponent glasses. We identify the convergence range for specific properties and use them to discuss models of several glass compositions for which a possible influence of cooling-rate or size effects had been previously hypothesized. The trends highlighted here represent an important reference to obtain reliable models of multicomponent glasses and extract converged medium-range structural features which affect the glass degradation and thus their application in different fields. In addition, as a first application of the present findings, the fully converged structure of the 45S5 glass was further analyzed to shed new light on several dissolution-related features whose interpretation has been rather controversial in the past.