The density functional theory (DFT) provides us an effective way to calculate large cluster systems with moderate computational demands. We calculate potential energy surfaces (PES) with several different approaches of DFT. The PES in the ground electronic state are related to HCP's isomerization process. The calculated PES are compared with the “experimental” PES obtained by fitting from the experimental vibrational spectra and that given by the “accurate” quantum chemistry calculation with more expensive computations. The comparisons show that the potential surfaces calculated with DFT methods can reach the accuracy of less than 0.1 eV.