Theoretical and Computational Investigations on Gas-Phase Thermal Tautomerization of Imidazole-Acetic Acid

Abstract

The gas-phase thermal tautomerization reaction between imidazole-4-acetic (I) and imidazole-5-acetic (II) acids was monitored using the traditional hybrid functional (B3LYP) and the long-range corrected functionals (CAM-B3LYP and ?B97XD) with 6-311++G** and aug-cc-pvdz basis sets. The roles of the long-range and dispersion corrections on their geometrical parameters, thermodynamic functions, kinetics, dipole moments, Highest Occupied Molecular Orbital–Lowest Unoccupied Molecular Orbital (HOMO–LUMO) energy gaps and total hyperpolarizability were investigated. All tested levels of theory predicted the preference of I over II by 0.750–0.877 kcal/mol. The origin of predilection of I is assigned to the H-bonding interaction ( ). This interaction stabilized I by 15.07 kcal/mol. The gas-phase interconversion between the two tautomers assumed a 1,2-proton shift mechanism, with two transition states (TS), TS1 and TS2, having energy barriers of 47.67–49.92 and 49.55–52.69 kcal/mol, respectively, and an sp³-type intermediate. A water-assisted 1,3-proton shift route brought the barrier height down to less than 20 kcal/mol in gas-phase and less than 12 kcal/mol in solution. The relatively high values of total hyperpolarizability of I compared to II were interpreted and discussed. The long-range correction is quite essential for total hyperpolarizability DFT evaluation as the traditional hybrid functional overestimate it. Finally, the computation of _tot values is basis set independent.