In this paper, high-level ab initiocalculations by using multi-configuration self-consistent field method with atomic effective core potential, polarization potential, and uncontracted Gaussian basis function, are performed to compute the potential energy curves of a total of 36 low-lying ∧-S states with ∑ _g,u, Π _g,u, △ _g,usymmetries of Na ₂ ⁺cation associated with the lowest 9 dissociation limits Na (3s, 3p, 4s, 3d, 4p, 5s, 4d, 4f, 5p)+Na ⁺. On the basis of the potential energy curves, the spectroscopic constants ( T _e, R _e, ω _e, ω _e χ _e, B _e, α _e, D _e) of the bound states are determined, which are in good agreement with the existing available experimental and theoretical values. Our results indicate that 5 ²∑ _g ⁺-7 ²∑ _g ⁺, 3 ²∑ _u ⁺-7 ²∑ _u ⁺, 2 ²Π _g, 4 ²Π _g, 1 ²△ _uand 2 ²△ _ustates are repulsive, which supports Berriche's results, and we report 10 electron states for the first time, that is, 8 ²∑ _{g, u} ⁺-9 ²∑ _{g, u} ⁺, 5 ²Π _{g, u}-7 ²Π _{g, u}and 3 ²△ _{g, u}. The vibrational-rotational spectroscopic constants and lowest vibrational-rotational energy levels ( ν=0-20) of the bound states are also presented. Moreover, in order to illustrate the strong state interactions of adjacent states with same symmetry, the information about the avoided crossing points is shown in detail. Finally, the transition dipole moments from a few low-lying excited states (1 ²Π _u-3 ²Π _u) to the ground state X ²∑ _g ⁺are computed. Therefore, it is expected that our computational results in the present calculations are significant for the molecular spectroscopy, ion-atom interaction and molecular cold collision fields.