Change search
ReferencesLink to record
Permanent link

Direct link
Paradynamics: an effective and reliable model for ab initio QM/MM free-energy calculations and related tasks
Stockholm University, Faculty of Science, Department of Organic Chemistry.
2011 (English)In: Journal of Physical Chemistry B, ISSN 1520-6106, E-ISSN 1520-5207, Vol. 115, no 24, 7950-7962 p.Article in journal (Refereed) Published
Abstract [en]

Recent years have seen tremendous effort in the development of approaches with which to obtain quantum mechanics/molecular mechanics (QM/MM) free energies for reactions in the condensed phase. Nevertheless, there remain significant challenges to address, particularly, the high computational cost involved in performing proper configurational sampling and, in particular, in obtaining ab initio QM/MM (QM(ai)/MM) free-energy surfaces. One increasingly popular approach that seems to offer an ideal way to progress in this direction is the elegant metadynamics (MTD) approach. However, in the current work, we point out the subtle efficiency problems associated with this approach and illustrate that we have at hand what is arguably a more powerful approach. More specifically, we demonstrate the effectiveness of an updated version of our original idea of using a classical reference potential for QM(ai)/MM calculations [J. Phys. Chem. 1995, 99, 17516)], which we refer to as paradynamics (PD). This approach is based on the use of an empirical valence bond (EVB) reference potential, which is already similar to the real ab initio potential. The reference potential is fitted to the ab initio potential by an iterative and, to a great degree, automated refinement procedure. The corresponding free-energy profile is then constructed using the refined EVB potential, and the linear response approximation (LRA) is used to evaluate the QM(ai)/MM activation free-energy barrier. The automated refinement of the EVB surface (and thus the reduction of the difference between the reference and ab initio potentials) is a key factor in accelerating the convergence of the LRA approach. We apply our PD approach to a test reaction, namely, the SN2 reaction between a chloride ion and methyl chloride, and demonstrate that, at present, this approach is far more powerful and cost-effective than the metadynamics approach (at least in its current implementation). We also discuss the general features of the PD approach in terms of its ability to explore complex systems and clarify that it is not a specialized approach limited to only accelerating QM(ai)/MM calculations with proper sampling, but rather can be used in a wide variety of applications. In fact, we point out that the use of a reference (CG) potential coupled with its PD refinement, as well as our renormalization approach, provides very general and powerful strategies that can be used very effectively to explore any property that has been studied by the MTD approach.

Place, publisher, year, edition, pages
2011. Vol. 115, no 24, 7950-7962 p.
National Category
Organic Chemistry
URN: urn:nbn:se:su:diva-62752DOI: 10.1021/jp201217bISI: 000291709500027OAI: diva2:444278
Swedish Research Council
Available from: 2011-09-28 Created: 2011-09-28 Last updated: 2012-01-24Bibliographically approved

Open Access in DiVA

No full text

Other links

Publisher's full text

Search in DiVA

By author/editor
Kamerlin, Shina C. L.
By organisation
Department of Organic Chemistry
In the same journal
Journal of Physical Chemistry B
Organic Chemistry

Search outside of DiVA

GoogleGoogle Scholar
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Altmetric score

Total: 32 hits
ReferencesLink to record
Permanent link

Direct link