Structural study of Na 2 O-B 2 O 3 -SiO 2 glasses from molecular simulations using a polarizable force field.

Sodium borosilicate glasses Na2 O-B2 O3 -SiO2 (NBS) are complex systems from a structural point of view. Three main building units are present: tetrahedral SiO4 and BO4 (BIV ) and triangular BO3 (BIII ). One of the salient features of these compounds is the change of the BIII /BIV 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 BIII /BIV with respect to the composition and the temperature.