Diatomic molecules , are common compositions in plasmas using the helium as the working gas, and play key roles in quenching O2 and N2. However, the potential energy curve (PEC) of the first state obtained at the CASSCF level of theory is not accuracy and data of higher excited states of are lacking. Here, we calculated PECs of four excited states of and nine lowest electronic states of by employing CASSCF+MRCI/Aug‐cc‐pV5Z and CCSD(T)/Aug‐cc‐pV5Z level of theory with sample internuclear distance interval of 0.005 Å. PECs were extrapolated to the complete basis sets. There are double local maxima in PECs of three excited states , , . The machine learning method based on the Gaussian process was used to interpolate and extrapolate the PECs, and then vibrational wave functions are solved numerically. This ab initio investigation provided rotational‐vibrational constants and new data for experimental research in the future.