Direct Validation of the Single Step Classical to Quantum Free Energy Perturbation

详细信息    查看全文

• 作者：Christopher Cave-Ayland ; Chris-Kriton Skylaris ; Jonathan W. Essex
• 刊名：Journal of Physical Chemistry B
• 出版年：2015
• 出版时间：January 22, 2015
• 年：2015
• 卷：119
• 期：3
• 页码：1017-1025
• 全文大小：408K
• ISSN：1520-5207

文摘

The use of the Zwanzig equation in the calculation of single-step perturbations to provide first-principles (ab initio) quantum mechanics (QM) correction terms to molecular mechanics (MM) free energy cycles is well established. A rigorous test of the ability to converge such calculations would be very useful in this context. In this work, we perform a direct assessment of the convergence of the MM to QM perturbation, by attempting the reverse QM to MM perturbation. This required the generation of extensive QM molecular dynamics trajectories, using density functional theory (DFT), within the representative biological system of a DNA adenosine鈥搕hymidine dimer. Over 100 ps of dynamics with the PBE functional and 6.25 ps with the LDA functional were generated. We demonstrate that calculations with total potential energies are very poorly convergent due to a lack of overlap of phase space distributions between ensembles. While not theoretically rigorous, the use of interaction energies provides far superior convergence, despite the presence of nonclassical charge transfer effects within the DFT trajectories. The source of poor phase space overlap for total energies is diagnosed, the approximate quantification of overlaps suggesting that even for the comparatively simple system considered here convergence of total energy calculations within a reasonable simulation time is unfeasible.