Deriving Force-Field Parameters from First Principles Using a Polarizable and Higher Order Dispersion Model

• Published in: Journal of chemical theory and computation Vol. 15; no. 3; pp. 1875 - 1883
• Main Authors: Visscher, Koen M, Geerke, Daan P
• Format: Journal Article
• Language: English
• Published: United States American Chemical Society 12.03.2019
• Subjects: Alkanes; Atomic radius; Computer simulation; Dipole moments; Dispersion; First principles; Mathematical models; Molecular dynamics; Parameterization; Parameters; Quantum mechanics
• ISSN: 1549-9618; 1549-9626
• DOI: 10.1021/acs.jctc.8b01105

In this work we propose a strategy based on quantum mechanical (QM) calculations to parametrize a polarizable force field for use in molecular dynamics (MD) simulations. We investigate the use of multiple atoms-in-molecules (AIM) strategies to partition QM determined molecular electron densities into atomic subregions. The partitioned atomic densities are subsequently used to compute atomic dispersion coefficients from effective exchange-hole-dipole moment (XDM) calculations. In order to derive values for the repulsive van der Waals parameters from first principles, we use a simple volume relation to scale effective atomic radii. Explicit inclusion of higher order dispersion coefficients was tested for a series of alkanes, and we show that combining C 6 and C 8 attractive terms together with a C 11 repulsive potential yields satisfying models when used in combination with our van der Waals parameters and electrostatic and bonded parameters as directly obtained from quantum calculations as well. This result highlights that explicit inclusion of higher order dispersion terms could be viable in simulation, and it suggests that currently available QM analysis methods allow for first-principles parametrization of molecular mechanics models.