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 language | English |
|---|---|
| Pages (from-to) | 3917-3928 |
| Number of pages | 12 |
| Journal | Journal of chemical theory and computation |
| Volume | 15 |
| Issue number | 7 |
| DOIs | |
| Publication status | Published - 2019 Jul 9 |
Fingerprint
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 journal › Article
}
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 -