TY - JOUR
T1 - Quantitative Analysis of QM/MM Boundary Artifacts and Correction in Adaptive QM/MM Simulations
AU - Watanabe, Hiroshi C.
AU - Cui, Qiang
N1 - Funding Information:
We thank Yohich Suzuki and Yutaka Shikano for the fruitful discussions. This research was supported by JSPS KAKENHI Grant Number JP17K15101 and JST PRESTO Grant Number JPMJPR17GC.
Publisher Copyright:
© 2019 American Chemical Society.
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
C2 - 31095912
AN - SCOPUS:85069293895
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 -