Molecular dynamics simulations of sodium silicate glasses: Optimization and limits of the computational procedure

The performance of the molecular dynamics (MD) simulations to obtain the structure of silica glasses containing different concentrations of alkali oxides has been tested. An optimal MD simulation procedure (including cooling cycle, MD constants and ensemble used) has been developed by means of experimental design methodologies (DOE), firstly restricting the study to the 30% Na2O silica glass, for which experimental data are available to allow the comparison of the results. The optimization procedure led to simulations that well predict experimental density and short-range structure of glasses with different sodium content. On the contrary, the medium-range structure has been badly reproduced and it was not possible to determine a reliable correlation with the parameters of the simulation procedures used. Therefore, the correlation of medium-range properties with the structure has been studied by means of the PLS methods. The results showed to be useful to highlight the relationships among structural elements, such as radial distribution functions of specific bonds and angles, and Qn species, suggesting possible directions for force fields improvements.