Density functional theory, B3LYPr6-31G** and B3LYPr6-311 q G 2d,p , and ab initio MP2r6-31G** calculations have been carried out to investigate the conformers, transition states, and energy barriers of the conformational processes of oxalic acid and its anions. QCISDr6-31G** geometrical optimization is also performed in the stable forms. Its calculated energy differences between the two most stable conformers are very near to the related observed value at 7.0 kJrmol. It is found that the structures and relative energies of oxalic acid conformers predicted by these methods show similar Ž . results, and that the conformer L1 C with the double-interfunctional-groups hydrogen 2 h bonds is the most stable conformer. The magnitude of hydrogen bond energies depends on the energy differences of various optimized structures. The hydrogen bond energies will be about 32 kJrmol for interfunctional groups, 17 kJrmol for weak interfunctional Ž . groups, 24 kJrmol for intra-COOH in COOH , and 60 kJrmol for interfunctional 2 Ž . y1 Ž . groups in COOH COO ion if calculated using the B3LYPr6-311 q G 2d,p method.