A novel two-component Equation-Of-Motion CCSD method for Ionization-Potential (EOMIP) based on a previously developed CC method with spin-orbital coupling (SOC)
[1-4] included in post-SCF part based on CCSD is implemented in the CFOUR program package
[5]. This approach could be use to calculate both ground and some excited states energies for open-shell systems with SOC. The EOMIP-CC method is computationally more efficient than methods that include SOC in SCF part
[6,7] and free from spin contaminant due to the closed-shell nature of reference states. SOC in the present approach is introduced either in the CC part or in the EOMIP part. The later implementation is particularly efficient since the computational effort for the EOMIP part is much smaller than that of the CCSD calculation without SOC. This approach has been used to calculate some open-shell systems. The splitting value of lowest
2P term for Kr
+, Xe
+, Rn
+, Br, I and At, as well as the splitting value of lowest
2D term for Zn
+, Cd
+, Hg
+, Cu, Ag, and Au. Meanwhile, the ionization energy contributed from s
2p
6 to s
2p
5 for Kr, Xe, Rn and the ionization energy contributed from d
10s
2 to d
9s
2 for Zn, Cd, Hg are calculated. Taking I
2+ as a typical molecule, the equilibrium bond length and harmonic frequency for five lowest state of
2Π
g,3/2,
2Π
g,1/2,
2Π
u,3/2,
2Π
u,1/2 and
2∑
g+ are also calculated. The computational results show that the method with SOC included in EOMIP can provide reliable results for atoms and molecules containing even sixth row elements.
[1] F. Wang, J. Gauss and C. van Wüllen, J. Chem. Phys., 129, 064113(2008)
[2] F. Wang and J. Gauss, J. Chem. Phys., 129, 174110(2008)
[3] F. Wang and J. Gauss, J. Chem. Phys., 131, 164113(2009)
[4] Z. Y. Tu, F. Wang and J. Gauss, unpublished
[5] http://slater.chemie.uni-mainz.de/cfour/
[6] L. Visscher, E. Eliav and U. Kaldor, J. Chem. Phys., 115, 9720(2001)
[7] S. Hirata, T. Yana, R. J. Harrison, M. Kamiya, and P. D. Fan, J. Chem. Phys., 126, 024104(2007)
