Structural study of NaO-BO-SiO glasses from molecular simulations using a polarizable force field.

Sodium borosilicate glasses NaO-BO-SiO (NBS) are complex systems from a structural point of view. Three main building units are present: tetrahedral SiO and BO (B) and triangular BO (B). One of the salient features of these compounds is the change of the B/B ratio with the alkali concentration, which is very difficult to capture in force fields-based molecular dynamics simulations. In this work, we develop a polarizable force field that is able to reproduce the boron coordination and more generally the structure of several NBS systems in the glass and in the melt. The parameters of the potential are fitted from density functional theory calculations only, in contrast with the existing empirical potentials for NBS systems. This ensures a strong improvement on the transferability of the parameters from one composition to another. Using this new force field, the structure of NBS systems is validated against neutron diffraction and nuclear magnetic resonance experiments. A special focus is given to the distribution of B/B with respect to the composition and the temperature.