Efficient treatment of long-range electrostatics in charge equilibration approaches

A charge equilibration method based on real-space Gaussians as charge densities is presented. The implementation is part of the Electrode package available in the Large-scale Atomic/Molecular Massively Parallel Simulator and benefits from its efficient particle-mesh Ewald approach. A simple strategy required to switch from the previously used Slater-type orbital (STO) shielding to Gaussians is provided by fitting the Coulomb energy of two Gaussian charge distributions to the repulsion between two STOs. Their widths were optimized for O, Si, and Ti species, obtaining results consistent with previous studies using STOs in the case of SiO2 polymorphs. In the limit of sufficiently narrow Gaussians, it is shown that the implementation converges to electronegativity equalization method results for the case of Ti/TiOx interfaces. The method presented is implemented in a way that is potentially beneficial for the application of modern machine-learning force fields that include long-range electrostatic interactions.

DDC Class

541: Physical; Theoretical

530.41: Mechanics of Solids

004: Computer Sciences

620.11: Engineering Materials

Publication version

publishedVersion

Lizenz

https://creativecommons.org/licenses/by/4.0/

Name

174108_1_5.0255823.pdf

Size

4.95 MB

Format

Adobe PDF

Options

Efficient treatment of long-range electrostatics in charge equilibration approaches