“…[17][18][19] However, if accurate values of spectroscopic constants are wanted, one would better use iterative method instead of perturbation method to treat spin-orbit coupling. In fact, Kullie has already used a relativistic four-component Dirac-Coulomb Hamiltonian with spin-orbit coupling in the framework of time-dependent density functional theory and linear response approximation to calculate the electronic state of Zn2 20 and Cd2. 21 However, it is widely known that the results of spectroscopic constants are very dependent of the choice of functional.…”