Ab initio based polarizable force field generation and application to liquid silica and magnesia

We extend the program potfit, which generates effective atomic interaction potentials from ab initio data, to electrostatic interactions and induced dipoles. The potential parametrization algorithm uses the Wolf direct, pairwise summation method with spherical truncation. The polarizability of oxygen atoms is modeled with the Tangney-Scandolo interatomic force field approach. Due to the Wolf summation, the computational effort in simulation scales linearly in the number of particles, despite the presence of electrostatic interactions. Thus, this model allows to perform large-scale molecular dynamics simulations of metal oxides with realistic potentials. Details of the implementation are given, and the generation of potentials for SiO₂ and MgO is demonstrated. The approach is validated by simulations of microstructural, thermodynamic, and vibrational properties of liquid silica and magnesia.