Multipolar Ewald methods, 2: Applications using a quantum mechanical force field

Timothy J. Giese, Maria T. Panteva, Haoyuan Chen, Darrin York

Research output: Contribution to journalArticle

14 Citations (Scopus)

Abstract

A fully quantum mechanical force field (QMFF) based on a modified "divide-and-conquer" (mDC) framework is applied to a series of molecular simulation applications, using a generalized Particle Mesh Ewald method extended to multipolar charge densities. Simulation results are presented for three example applications: liquid water, p-nitrophenylphosphate reactivity in solution, and crystalline N,N-dimethylglycine. Simulations of liquid water using a parametrized mDC model are compared to TIP3P and TIP4P/Ew water models and experiment. The mDC model is shown to be superior for cluster binding energies and generally comparable for bulk properties. Examination of the dissociative pathway for dephosphorylation of p-nitrophenylphosphate shows that the mDC method evaluated with the DFTB3/3OB and DFTB3/OPhyd semiempirical models bracket the experimental barrier, whereas DFTB2 and AM1/d-PhoT QM/MM simulations exhibit deficiencies in the barriers, the latter for which is related, in part, to the anomalous underestimation of the p-nitrophenylate leaving group pKa. Simulations of crystalline N,N-dimethylglycine are performed and the overall structure and atomic fluctuations are compared with the experiment and the general AMBER force field (GAFF). The QMFF, which was not parametrized for this application, was shown to be in better agreement with crystallographic data than GAFF. Our simulations highlight some of the application areas that may benefit from using new QMFFs, and they demonstrate progress toward the development of accurate QMFFs using the recently developed mDC framework.

Original languageEnglish (US)
Pages (from-to)451-461
Number of pages11
JournalJournal of Chemical Theory and Computation
Volume11
Issue number2
DOIs
StatePublished - Feb 10 2015

Fingerprint

field theory (physics)
Water
simulation
Crystalline materials
Liquids
Charge density
Binding energy
water
Experiments
brackets
liquids
mesh
reactivity
examination
binding energy
nitrophenylphosphate
dimethylglycine

All Science Journal Classification (ASJC) codes

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Giese, Timothy J. ; Panteva, Maria T. ; Chen, Haoyuan ; York, Darrin. / Multipolar Ewald methods, 2 : Applications using a quantum mechanical force field. In: Journal of Chemical Theory and Computation. 2015 ; Vol. 11, No. 2. pp. 451-461.
@article{f98411cce7cf4ac59730ac6bb2aa081e,
title = "Multipolar Ewald methods, 2: Applications using a quantum mechanical force field",
abstract = "A fully quantum mechanical force field (QMFF) based on a modified {"}divide-and-conquer{"} (mDC) framework is applied to a series of molecular simulation applications, using a generalized Particle Mesh Ewald method extended to multipolar charge densities. Simulation results are presented for three example applications: liquid water, p-nitrophenylphosphate reactivity in solution, and crystalline N,N-dimethylglycine. Simulations of liquid water using a parametrized mDC model are compared to TIP3P and TIP4P/Ew water models and experiment. The mDC model is shown to be superior for cluster binding energies and generally comparable for bulk properties. Examination of the dissociative pathway for dephosphorylation of p-nitrophenylphosphate shows that the mDC method evaluated with the DFTB3/3OB and DFTB3/OPhyd semiempirical models bracket the experimental barrier, whereas DFTB2 and AM1/d-PhoT QM/MM simulations exhibit deficiencies in the barriers, the latter for which is related, in part, to the anomalous underestimation of the p-nitrophenylate leaving group pKa. Simulations of crystalline N,N-dimethylglycine are performed and the overall structure and atomic fluctuations are compared with the experiment and the general AMBER force field (GAFF). The QMFF, which was not parametrized for this application, was shown to be in better agreement with crystallographic data than GAFF. Our simulations highlight some of the application areas that may benefit from using new QMFFs, and they demonstrate progress toward the development of accurate QMFFs using the recently developed mDC framework.",
author = "Giese, {Timothy J.} and Panteva, {Maria T.} and Haoyuan Chen and Darrin York",
year = "2015",
month = "2",
day = "10",
doi = "10.1021/ct500799g",
language = "English (US)",
volume = "11",
pages = "451--461",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "2",

}

Multipolar Ewald methods, 2 : Applications using a quantum mechanical force field. / Giese, Timothy J.; Panteva, Maria T.; Chen, Haoyuan; York, Darrin.

In: Journal of Chemical Theory and Computation, Vol. 11, No. 2, 10.02.2015, p. 451-461.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Multipolar Ewald methods, 2

T2 - Applications using a quantum mechanical force field

AU - Giese, Timothy J.

AU - Panteva, Maria T.

AU - Chen, Haoyuan

AU - York, Darrin

PY - 2015/2/10

Y1 - 2015/2/10

N2 - A fully quantum mechanical force field (QMFF) based on a modified "divide-and-conquer" (mDC) framework is applied to a series of molecular simulation applications, using a generalized Particle Mesh Ewald method extended to multipolar charge densities. Simulation results are presented for three example applications: liquid water, p-nitrophenylphosphate reactivity in solution, and crystalline N,N-dimethylglycine. Simulations of liquid water using a parametrized mDC model are compared to TIP3P and TIP4P/Ew water models and experiment. The mDC model is shown to be superior for cluster binding energies and generally comparable for bulk properties. Examination of the dissociative pathway for dephosphorylation of p-nitrophenylphosphate shows that the mDC method evaluated with the DFTB3/3OB and DFTB3/OPhyd semiempirical models bracket the experimental barrier, whereas DFTB2 and AM1/d-PhoT QM/MM simulations exhibit deficiencies in the barriers, the latter for which is related, in part, to the anomalous underestimation of the p-nitrophenylate leaving group pKa. Simulations of crystalline N,N-dimethylglycine are performed and the overall structure and atomic fluctuations are compared with the experiment and the general AMBER force field (GAFF). The QMFF, which was not parametrized for this application, was shown to be in better agreement with crystallographic data than GAFF. Our simulations highlight some of the application areas that may benefit from using new QMFFs, and they demonstrate progress toward the development of accurate QMFFs using the recently developed mDC framework.

AB - A fully quantum mechanical force field (QMFF) based on a modified "divide-and-conquer" (mDC) framework is applied to a series of molecular simulation applications, using a generalized Particle Mesh Ewald method extended to multipolar charge densities. Simulation results are presented for three example applications: liquid water, p-nitrophenylphosphate reactivity in solution, and crystalline N,N-dimethylglycine. Simulations of liquid water using a parametrized mDC model are compared to TIP3P and TIP4P/Ew water models and experiment. The mDC model is shown to be superior for cluster binding energies and generally comparable for bulk properties. Examination of the dissociative pathway for dephosphorylation of p-nitrophenylphosphate shows that the mDC method evaluated with the DFTB3/3OB and DFTB3/OPhyd semiempirical models bracket the experimental barrier, whereas DFTB2 and AM1/d-PhoT QM/MM simulations exhibit deficiencies in the barriers, the latter for which is related, in part, to the anomalous underestimation of the p-nitrophenylate leaving group pKa. Simulations of crystalline N,N-dimethylglycine are performed and the overall structure and atomic fluctuations are compared with the experiment and the general AMBER force field (GAFF). The QMFF, which was not parametrized for this application, was shown to be in better agreement with crystallographic data than GAFF. Our simulations highlight some of the application areas that may benefit from using new QMFFs, and they demonstrate progress toward the development of accurate QMFFs using the recently developed mDC framework.

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

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

U2 - 10.1021/ct500799g

DO - 10.1021/ct500799g

M3 - Article

C2 - 25691830

AN - SCOPUS:84922639701

VL - 11

SP - 451

EP - 461

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 2

ER -