First-principles study of a sodium borosilicate glass-former. I. The liquid state

We use ab initio simulations to study the static and dynamic properties of a sodium borosilicate liquid with composition 3Na(2)O-B2O3-6SiO(2), i.e., a system that is the basis of many glass-forming materials. In particular, we focus on the question how boron is embedded into the local structure of the silicate network liquid. From the partial structure factors we conclude that there is a weak nanoscale phase separation between silicon and boron and that the sodium atoms form channel-like structures as they have been found in previous studies of sodosilicate glass-formers. Our results for the x-ray and neutron structure factor show that this feature is basically not detectable in the former but should be visible in the latter as a small peak at small wave vectors. At high temperatures we find a high concentration of threefold coordinated boron atoms which decreases rapidly with decreasing T, whereas the number of fourfold coordinated boron atoms increases. Therefore, we conclude that at the experimental glass transition temperature most boron atoms will be fourfold coordinated. We show that the transformation of ([3]) B into ([4]) B with decreasing T is not just related to the diminution of nonbridging oxygen atoms as claimed in previous studies, but to a restructuring of the silicate matrix. The diffusion constants of the various elements show an Arrhenius behavior and we find that the one for boron has the same value as the one of oxygen and is significantly larger than the one of silicon. This shows that these two network-formers have rather different dynamical properties, a result that is also confirmed from the time dependence of the van Hove functions. Finally, we show that the coherent intermediate scattering function for the sodium atoms is very different from the incoherent one and that it tracks the one of the matrix atoms.