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

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 adenosinethymidine 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.