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

Abstract

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-Sheffler scheme for dispersion interactions.

arXiv comment61 pages, 12 figures, 2 tables

References in corpus (6)

• Canonical sampling through velocity-rescaling
• Separable Dual Space Gaussian Pseudo-potentials
• A real-space grid implementation of the Projector Augmented Wave method
• $C_6$ coefficients and dipole polarizabilities for all atoms and many ions in rows 1-6 of the periodic table
• Extending Hirshfeld-I to bulk and periodic materials
• Atomic charges for modeling metal-organic frameworks: Why and how

Cited by in corpus (19)

• Non-covalent interactions across organic and biological subsets of chemical space: Physics-based potentials parametrized from machine learning
• Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles
• High rate nanofluidic energy absorption in porous zeolitic frameworks
• Force-Field-Enhanced Neural Network Interactions: from Local Equivariant Embedding to Atom-in-Molecule properties and long-range effects
• Modeling the structural and thermal properties of loaded metal-organic frameworks. An interplay of quantum and anharmonic fluctuations
• Accurate Deep Learning-aided Density-free Strategy for Many-Body Dispersion-corrected Density Functional Theory
• Hybrid Classical/Machine-Learning Force Fields for the Accurate Description of Molecular Condensed-Phase Systems
• FeNNol: an Efficient and Flexible Library for Building Force-field-enhanced Neural Network Potentials
• DASH: Dynamic Attention-Based Substructure Hierarchy for Partial Charge Assignment
• On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment
• Multi-center decomposition of molecular densities: a mathematical perspective
• Generalized Many-Body Dispersion Correction through Random-phase Approximation for Chemically Accurate Density Functional Theory
• A new framework for frequency-dependent polarizable force fields
• Embedding quantum statistical excitations in a classical force field
• Multi-center decomposition of molecular densities: A numerical perspective
• Approximations of the Iterative Stockholder Analysis scheme using exponential basis functions
• Minimally-corrected partial atomic charges for non-covalent electrostatic interactions
• Learning Atomic Multipoles: Prediction of the Electrostatic Potential with Equivariant Graph Neural Networks
• Accurate transferable polarization model derived from the monomer electron density