Minimal Basis Iterative Stockholder: Atoms in Molecules for Force-Field Development

Atomic partial charges appear in the Coulomb term of many force-field models and can be derived from electronic structure calculations with a myriad of atoms-in-molecules (AIM) methods. More advanced models have also been proposed, using the distributed nature of the electron cloud and atomic multip...

Full description

Saved in:
Bibliographic Details
Published in:Journal of chemical theory and computation 2016-08, Vol.12 (8), p.3894-3912
Main Authors: Verstraelen, Toon, Vandenbrande, Steven, Heidar-Zadeh, Farnaz, Vanduyfhuys, Louis, Van Speybroeck, Veronique, Waroquier, Michel, Ayers, Paul W
Format: Article
Language:English
Subjects:
Approximation
Coulomb friction
Electron density
Electronic structure
Electrostatics
Iterative methods
Mathematical analysis
Mathematical models
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Atomic partial charges appear in the Coulomb term of many force-field models and can be derived from electronic structure calculations with a myriad of atoms-in-molecules (AIM) methods. More advanced models have also been proposed, using the distributed nature of the electron cloud and atomic multipoles. In this work, an electrostatic force field is defined through a concise approximation of the electron density, for which the Coulomb interaction is trivially evaluated. This approximate “pro-density” is expanded in a minimal basis of atom-centered s-type Slater density functions, whose parameters are optimized by minimizing the Kullback–Leibler divergence of the pro-density from a reference electron density, e.g., obtained from an electronic structure calculation. The proposed method, Minimal Basis Iterative Stockholder (MBIS), is a variant of the Hirshfeld AIM method, but it can also be used as a density-fitting technique. An iterative algorithm to refine the pro-density is easily implemented with a linear-scaling computational cost, enabling applications to supramolecular systems. The benefits of the MBIS method are demonstrated with systematic applications to molecular databases and extended models of condensed phases. A comparison to 14 other AIM methods shows its effectiveness when modeling electrostatic interactions. MBIS is also suitable for rescaling atomic polarizabilities in the Tkatchenko–Scheffler scheme for dispersion interactions.
ISSN:1549-9618
1549-9626
DOI:10.1021/acs.jctc.6b00456