The present paper is organized as follows. In Section 2, the methods for the determination of the kinetic energy operator, the electronic structure, the potential energy surface and the computational details are presented. The results are discussed and concluding remarks are given in Section 3. In the Appendix, the fitting procedure of PES is given.