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

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

field theory (physics)
inclusions
Molecular mechanics
Alkanes
Dipole moment
coefficients
Paraffins
alkanes
Carrier concentration
Molecular dynamics
Electrostatics
partitions
dipole moments
simulation
electrostatics
molecular dynamics
Atoms
Molecules
radii
Computer simulation

Cite this

@article{92508efd0b8a4fe393d439e7b6417632,
title = "Deriving Force-Field Parameters from First Principles Using a Polarizable and Higher Order Dispersion Model",
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.",
author = "Visscher, {Koen M.} and Geerke, {Daan P.}",
year = "2019",
month = "2",
day = "14",
doi = "10.1021/acs.jctc.8b01105",
language = "English",
volume = "15",
pages = "1875–1883",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "3",

}

Deriving Force-Field Parameters from First Principles Using a Polarizable and Higher Order Dispersion Model. / Visscher, Koen M.; Geerke, Daan P.

In: Journal of Chemical Theory and Computation, Vol. 15, No. 3, 14.02.2019, p. 1875–1883.

Research output: Contribution to JournalArticleAcademicpeer-review

TY - JOUR

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

AU - Visscher, Koen M.

AU - Geerke, Daan P.

PY - 2019/2/14

Y1 - 2019/2/14

N2 - 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.

AB - 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.

UR - http://www.scopus.com/inward/record.url?scp=85061962516&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85061962516&partnerID=8YFLogxK

U2 - 10.1021/acs.jctc.8b01105

DO - 10.1021/acs.jctc.8b01105

M3 - Article

VL - 15

SP - 1875

EP - 1883

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 3

ER -