Accelerated Computation of Free Energy Profile at Ab Initio Quantum Mechanical/Molecular Mechanics Accuracy via a Semiempirical Reference-Potential. 4. Adaptive QM/MM

03 November 2020, Version 2
This content is a preprint and has not undergone peer review at the time of posting.


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.


reference-potential method

Supplementary materials



Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.