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

Koen M. Visscher, Daan P. Geerke*

*Corresponding author for this work

Research output: Contribution to JournalArticleAcademicpeer-review

Abstract

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.

Original languageEnglish
Pages (from-to)1875–1883
Number of pages9
JournalJournal of Chemical Theory and Computation
Volume15
Issue number3
DOIs
Publication statusPublished - 14 Feb 2019

Fingerprint Dive into the research topics of 'Deriving Force-Field Parameters from First Principles Using a Polarizable and Higher Order Dispersion Model'. Together they form a unique fingerprint.

Cite this