These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
Accelerated Computation of Free Energy Profile at Ab Initio Quantum Mechanical/Molecular Mechanics Accuracy via a Semiempirical Reference-Potential. 4. Adaptive QM/MM
preprintrevised on 01.11.2020, 13:55 and posted on 03.11.2020, 05:38 by Jia-Ning Wang, Wei Liu, Pengfei Li, Yan Mo, Wenxin Hu, Jun Zheng, Xiaoliang Pan, Yihan Shao, Ye Mei
Although Quantum Mechanical/Molecular Mechanics (QM/MM) methods are now routinely applied to the studies of chemical reactions in condensed phases and enzymatic reactions, they may confront technical difficulties when the reactive region is varying over time. For instance, when the solvent molecules are participating in the reaction, the exchange of water molecules between the QM and MM regions may occur on a time scale that is comparable to that of the reaction. Several adaptive QM/MM schemes have been proposed to cope with this situation. However, these methods either significantly increase the computational cost or introducing unrealistic restraints to the system. In this work, we developed a novel adaptive QM/MM scheme and applied it to a study of the nucleophilic addition reaction. In this approach, the simulation was performed with a small QM region (without solvent molecules), and the thermodynamic properties under other potential energy functions with larger QM regions (with a different number of solvent molecules and/or different level of QM theory) are computed via extrapolation using the reference-potential method. The results show that this reweighting process is numerically stable, at least for the case studied in this work. Furthermore, this method also offers an inexpensive way to examine the convergence of the QM/MM calculation with respect to the size of the QM region.