Ab initio calculations of the pyrophosphate hydrolysis reaction

Ab initio quantum mechanical calculations were used to study the hydrolysis reaction H4P2O7 + H2O in equilibrium with 2H3PO4, as well as some molecular properties of the reactants and products. SCF calculations with several basis sets ranging from minimal to extended with polarization functions were used to look at the basis dependency of the reaction enthalpies and optimized geometries. Although the minimal basis sets yield erratic predictions of the enthalpy, when a more extended basis (3-21G*) was used for the geometry optimization, and the total energies of the reactants and products were computed with this and larger basis sets, we obtained more consistent predictions of the structural properties of the P-O-P bridge and of the heat of the hydrolysis reaction (delta E = -7.39 kcal/mol at the SCF/6-31G** level). A comparison is made with previous estimates performed with smaller basis sets and without taking into account the electron correlation effects, which are calculated in the present work. The inclusion of the zero point energy calculated using the harmonic approximation, and of the electronic correlation energy determined at the MBPT(2) level, raised the computed heat of the reaction to -3.83 kcal/mol, and when an estimate for the thermal energy was added, the value obtained was of -3.38 kcal/mol. In conclusion, we found that the hydrolysis of pyrophosphate should be exothermic in the gas phase. The implications of this result in relation to some recent theories about enzyme catalysis are discussed.