Transferable potential for molecular dynamics simulations of borosilicate glasses and structural comparison of machine learning optimized parameters

The simulation of borosilicate glasses is challenging due to the composition and temperature dependent coordination state of boron atoms. Here, we present a new machine-learning optimized classical two-body Buckingham potential for molecular dynamics simulations that achieves transferability across diverse borosilicate glass compositions. Our potential accurately predicts the glass structural variations in short- and medium-range order in different glass compositions, including validating our potential against experimental X-ray structure factor data. Notably, these data are not included in the optimization framework, which focuses exclusively on density and four-fold coordinated boron fraction. We further investigate the impact of empirical parameters in the force field formulation on the microscopic bond lengths, bond angles and the macroscopic densities, providing new insights into the relationship between interatomic potentials and bulk glass behaviors.