Double Proton Transfer using Dissociable Force Fields
The construction, implementation, and use of dissociable classical force fields are discussed. Starting from zeroth-order interaction potentials for O2H5+ and N2H7+ calculated with MP2/6–311++G**, energy scaling of the potential energy surfaces allows adjustment of quantities such as the barrier heights to describe a range of physical situations observed in realistic systems. As an example, ‘potential morphing’ is used to investigate the dynamics of double proton transfer in 2-pyridone · 2-hydroxypyridine for which previous estimates of the barrier to tautomerization are available. Scaling factors to give barrier heights for double proton transfer between 3.6 and 17.6 kcal mol−1 are chosen to demonstrate the utility of the method to describe a range of different barrier heights and shapes. Considerable savings in computing time can be achieved compared to alternative methods such as mixed quantum/classical methods.