Quantitative Analysis of QM/MM Boundary Artifacts and Correction in Adaptive QM/MM Simulations

Hiroshi Watanabe, Qiang Cui

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

A quantum chemical treatment of solvation effects using the standard quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations is challenging due to QM and MM solvent exchange near the QM solute. To this end, various adaptive QM/MM methods have been proposed; free solvent exchanges are allowed via flexible switching of their identities between QM and MM systems depending on their distances from the QM solute. However, temporal and spatial discontinuities remain in the standard implementations of adaptive QM/MM approaches and continue to hamper stable and accurate dynamics simulations. We previously demonstrated that the size-consistent multipartitioning (SCMP) method achieves temporal continuity while, to some extent, avoiding spatial discontinuities. In the present study, we demonstrate that residual spatial discontinuities may lead to severe artifacts under certain conditions. Through quantitative analyses, we show that all multiscale including QM/MM methods might experience these problems, which so far have not been investigated in depth. To alleviate these artifacts, we propose a correction scheme in the framework of the SCMP approach and demonstrate its effectiveness using bulk water simulations.

Original languageEnglish
Pages (from-to)3917-3928
Number of pages12
JournalJournal of chemical theory and computation
Volume15
Issue number7
DOIs
Publication statusPublished - 2019 Jul 9

Fingerprint

quantitative analysis
artifacts
discontinuity
Solvation
Computer simulation
solutes
Chemical analysis
Molecular dynamics
simulation
continuity
solvation
Water
molecular dynamics
water

ASJC Scopus subject areas

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Quantitative Analysis of QM/MM Boundary Artifacts and Correction in Adaptive QM/MM Simulations. / Watanabe, Hiroshi; Cui, Qiang.

In: Journal of chemical theory and computation, Vol. 15, No. 7, 09.07.2019, p. 3917-3928.

Research output: Contribution to journalArticle

@article{a1c0d21c5c7a45fab2dc235384d83f28,
title = "Quantitative Analysis of QM/MM Boundary Artifacts and Correction in Adaptive QM/MM Simulations",
abstract = "A quantum chemical treatment of solvation effects using the standard quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations is challenging due to QM and MM solvent exchange near the QM solute. To this end, various adaptive QM/MM methods have been proposed; free solvent exchanges are allowed via flexible switching of their identities between QM and MM systems depending on their distances from the QM solute. However, temporal and spatial discontinuities remain in the standard implementations of adaptive QM/MM approaches and continue to hamper stable and accurate dynamics simulations. We previously demonstrated that the size-consistent multipartitioning (SCMP) method achieves temporal continuity while, to some extent, avoiding spatial discontinuities. In the present study, we demonstrate that residual spatial discontinuities may lead to severe artifacts under certain conditions. Through quantitative analyses, we show that all multiscale including QM/MM methods might experience these problems, which so far have not been investigated in depth. To alleviate these artifacts, we propose a correction scheme in the framework of the SCMP approach and demonstrate its effectiveness using bulk water simulations.",
author = "Hiroshi Watanabe and Qiang Cui",
year = "2019",
month = "7",
day = "9",
doi = "10.1021/acs.jctc.9b00180",
language = "English",
volume = "15",
pages = "3917--3928",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "7",

}

TY - JOUR

T1 - Quantitative Analysis of QM/MM Boundary Artifacts and Correction in Adaptive QM/MM Simulations

AU - Watanabe, Hiroshi

AU - Cui, Qiang

PY - 2019/7/9

Y1 - 2019/7/9

N2 - A quantum chemical treatment of solvation effects using the standard quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations is challenging due to QM and MM solvent exchange near the QM solute. To this end, various adaptive QM/MM methods have been proposed; free solvent exchanges are allowed via flexible switching of their identities between QM and MM systems depending on their distances from the QM solute. However, temporal and spatial discontinuities remain in the standard implementations of adaptive QM/MM approaches and continue to hamper stable and accurate dynamics simulations. We previously demonstrated that the size-consistent multipartitioning (SCMP) method achieves temporal continuity while, to some extent, avoiding spatial discontinuities. In the present study, we demonstrate that residual spatial discontinuities may lead to severe artifacts under certain conditions. Through quantitative analyses, we show that all multiscale including QM/MM methods might experience these problems, which so far have not been investigated in depth. To alleviate these artifacts, we propose a correction scheme in the framework of the SCMP approach and demonstrate its effectiveness using bulk water simulations.

AB - A quantum chemical treatment of solvation effects using the standard quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations is challenging due to QM and MM solvent exchange near the QM solute. To this end, various adaptive QM/MM methods have been proposed; free solvent exchanges are allowed via flexible switching of their identities between QM and MM systems depending on their distances from the QM solute. However, temporal and spatial discontinuities remain in the standard implementations of adaptive QM/MM approaches and continue to hamper stable and accurate dynamics simulations. We previously demonstrated that the size-consistent multipartitioning (SCMP) method achieves temporal continuity while, to some extent, avoiding spatial discontinuities. In the present study, we demonstrate that residual spatial discontinuities may lead to severe artifacts under certain conditions. Through quantitative analyses, we show that all multiscale including QM/MM methods might experience these problems, which so far have not been investigated in depth. To alleviate these artifacts, we propose a correction scheme in the framework of the SCMP approach and demonstrate its effectiveness using bulk water simulations.

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

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

U2 - 10.1021/acs.jctc.9b00180

DO - 10.1021/acs.jctc.9b00180

M3 - Article

VL - 15

SP - 3917

EP - 3928

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 7

ER -