“…Thus, if one increases the basis set size (for example, including polarization and even diffuse functions in the trial function), then the potential curve of the ground-state H þ 2 should be improved significantly in all the regions of R. In Figure 1a, we demonstrate the effect of basis set size by using Hartree-Fock (HF) method implemented in the GAUSSIAN 09 electronic structure package [16]. It is seen that the numerical potential curve for the ground-state H þ 2 is greatly improved and approaching to the most accurate data [13,14] as the basis size is increased from GAUSSIANtype basis set 3-21G (2 sets of GAUSSIAN functions in the valence region), to 6-31G (d,p) (including polarization p and d functions), and then to 6-311+G (3df,2pd) (including 2p and 1d polarization functions and diffuse functions) [16].…”