The mutation was properly modelized via double hydrogen bonding complexes of 7-azaindole (7AI) with protic solvents. AIM and NBO analyses and the aromaticity calculations of 7AI complexes confirm that the 7AI-formic acid complex has the strongest interaction to achieve the tautomerization. The results show partially electrostatic-partially covalent nature for N⋯H hydrogen bond in 7AI-formic acid complex. Localized molecular orbital energy decomposition reveals an competitive correlation between polarization and electrostatic components of attraction energy to determine the nature of hydrogen bonds.