Deep machine learning interatomic potential for liquid silica

The use of machine learning to develop neural network potentials (NNP) representing the interatomic potential energy surface allows us to achieve an optimal balance between accuracy and efficiency in computer simulation of materials. A key point in developing such potentials is the preparation of a training dataset of ab initio trajectories. Here we apply a deep potential molecular dynamics (DeePMD) approach to develop NNP for silica, which is the representative glassformer widely used as a model system for simulating network-forming liquids and glasses. We show that the use of a relatively small training dataset of high-temperature ab initio configurations is enough to fabricate NNP, which describes well both structural and dynamical properties of liquid silica. In particular, we calculate the pair correlation functions, angular distribution function, velocity autocorrelation functions, vibrational density of states, and mean-square displacement and reveal a close agreement with ab initio data. We show that NNP allows us to expand significantly the time-space scales achievable in simulations and thus calculating dynamical and transport properties with more accuracy than that for ab initio methods. We find that developed NNP allows us to describe the structure of the glassy silica with satisfactory accuracy even though no low-temperature configurations were included in the training procedure. The results obtained open up prospects for simulating structural and dynamical properties of liquids and glasses via NNP.