Combining the spin-separated exact two-component relativistic Hamiltonian with the equation-of-motion coupled-cluster method for the treatment of spin-orbit splittings of light and heavy elements.

The spin-separated exact two-component (X2C) relativistic Hamiltonian [sf-X2C+so-DKHn, J. Chem. Phys., 2012, 137, 154114] is combined with the equation-of-motion coupled-cluster method with singles and doubles (EOM-CCSD) for the treatment of spin-orbit splittings of open-shell molecular systems. Scalar relativistic effects are treated to infinite order from the outset via the spin-free part of the X2C Hamiltonian (sf-X2C), whereas the spin-orbit couplings (SOC) are handled at the CC level via the first-order Douglas-Kroll-Hess (DKH) type of spin-orbit operator (so-DKH1). Since the exponential of single excitations, i.e., exp(T1), introduces sufficient spin orbital relaxations, the inclusion of SOC at the CC level is essentially the same in accuracy as the inclusion of SOC from the outset in terms of the two-component spinors determined variationally by the sf-X2C+so-DKH1 Hamiltonian, but is computationally more efficient. Therefore, such an approach (denoted as sf-X2C-EOM-CCSD(SOC)) can achieve uniform accuracy for the spin-orbit splittings of both light and heavy elements. For light elements, the treatment of SOC can even be postponed until the EOM step (denoted as sf-X2C-EOM(SOC)-CCSD), so as to further reduce the computational cost. To reveal the efficacy of sf-X2C-EOM-CCSD(SOC) and sf-X2C-EOM(SOC)-CCSD, the spin-orbit splittings of the 2Π states of monohydrides up to the sixth row of the periodic table are investigated. The results show that sf-X2C-EOM-CCSD(SOC) predicts very accurate results (within 5%) for elements up to the fifth row, whereas sf-X2C-EOM(SOC)-CCSD is useful only for light elements (up to the third row but with some exceptions). For comparison, the sf-X2C-S-TD-DFT-SOC approach [spin-adapted open-shell time-dependent density functional theory, Mol. Phys., 2013, 111, 3741] is applied to the same systems. The overall accuracy (1-10%) is satisfactory.

[1]  L. Visscher,et al.  Formulation and implementation of the relativistic Fock-space coupled cluster method for molecules , 2001 .

[2]  B. A. Hess,et al.  The relativistic Fock-space coupled-cluster method for molecules: CdH and its ions , 1998 .

[3]  J. Gauss,et al.  Communication: spin-orbit splittings in degenerate open-shell states via Mukherjee's multireference coupled-cluster theory: a measure for the coupling contribution. , 2012, The Journal of chemical physics.

[4]  H. Jensen,et al.  A direct relativistic four-component multiconfiguration self-consistent-field method for molecules. , 2008, The Journal of chemical physics.

[5]  Fan Wang,et al.  Equation of motion coupled cluster method for electron attached states with spin-orbit coupling , 2012 .

[6]  Wenjian Liu,et al.  Spin-adapted open-shell time-dependent density functional theory. III. An even better and simpler formulation. , 2011, The Journal of chemical physics.

[7]  Lucas Visscher,et al.  The molecular mean-field approach for correlated relativistic calculations. , 2009, The Journal of chemical physics.

[8]  Russell M. Pitzer,et al.  Electronic-structure methods for heavy-atom molecules , 1988 .

[9]  Hess,et al.  Relativistic electronic-structure calculations employing a two-component no-pair formalism with external-field projection operators. , 1986, Physical review. A, General physics.

[10]  Kirk A Peterson,et al.  Molecular core-valence correlation effects involving the post-d elements Ga-Rn: benchmarks and new pseudopotential-based correlation consistent basis sets. , 2010, The Journal of chemical physics.

[11]  Spin–orbit coupling and electron correlation in relativistic configuration interaction and coupled-cluster methods , 2012 .

[12]  Mark S. Gordon,et al.  MCSCF/6-31 G(d,p) Calculations of One-Electron Spin-Orbit Coupling Constants In Diatomic Molecules , 1992 .

[13]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[14]  Zhendong Li,et al.  Critical Assessment of TD-DFT for Excited States of Open-Shell Systems: I. Doublet-Doublet Transitions. , 2016, Journal of chemical theory and computation.

[15]  Fan Wang,et al.  The Beijing Density Functional (BDF) Program Package: Methodologies and Applications , 2003 .

[16]  P. Knowles,et al.  Spin-orbit matrix elements for internally contracted multireference configuration interaction wavefunctions , 2000 .

[17]  R. Bartlett,et al.  Coupled-cluster theory in quantum chemistry , 2007 .

[18]  C. Marian.  Spin–orbit coupling and intersystem crossing in molecules , 2012 .

[19]  Daoling Peng,et al.  Exact two-component Hamiltonians revisited. , 2009, The Journal of chemical physics.

[20]  J. Gauss,et al.  Spin-orbit couplings within the equation-of-motion coupled-cluster framework: Theory, implementation, and benchmark calculations. , 2015, The Journal of chemical physics.

[21]  J. Gauss,et al.  Closed-shell coupled-cluster theory with spin-orbit coupling. , 2008, The Journal of chemical physics.

[22]  J. Gauss,et al.  Perturbative calculation of spin-orbit splittings using the equation-of-motion ionization-potential coupled-cluster ansatz. , 2008, The Journal of chemical physics.

[23]  Russell M. Pitzer,et al.  Spin-Orbit Configuration Interaction Using the Graphical Unitary Group Approach and Relativistic Core Potential and Spin-Orbit Operators , 1999 .

[24]  Michael Dolg,et al.  The Beijing four-component density functional program package (BDF) and its application to EuO, EuS, YbO and YbS , 1997 .

[25]  Lee M J Huntington,et al.  Application of the multireference equation of motion coupled cluster method, including spin–orbit coupling, to the atomic spectra of Cr, Mn, Fe and Co , 2015 .

[26]  Roland Lindh,et al.  Main group atoms and dimers studied with a new relativistic ANO basis set , 2004 .

[27]  Wenjian Liu,et al.  Spin-adapted open-shell random phase approximation and time-dependent density functional theory. I. Theory. , 2010, The Journal of chemical physics.

[28]  Christoph van Wüllen,et al.  Accurate and efficient treatment of two-electron contributions in quasirelativistic high-order Douglas-Kroll density-functional calculations. , 2005, The Journal of chemical physics.

[29]  Werner Kutzelnigg,et al.  Quasirelativistic theory equivalent to fully relativistic theory. , 2005, The Journal of chemical physics.

[30]  Sandro Chiodo,et al.  One‐electron spin‐orbit contribution by effective nuclear charges , 2009, J. Comput. Chem..

[31]  Bernd A. Hess,et al.  Revision of the Douglas-Kroll transformation. , 1989, Physical review. A, General physics.

[32]  Wenjian Liu,et al.  On the spin separation of algebraic two-component relativistic Hamiltonians. , 2012, The Journal of chemical physics.

[33]  Christel M. Marian,et al.  A mean-field spin-orbit method applicable to correlated wavefunctions , 1996 .

[34]  Fan Wang,et al.  Symmetry exploitation in closed-shell coupled-cluster theory with spin-orbit coupling. , 2011, The Journal of chemical physics.

[35]  Wenjian Liu.  Advances in relativistic molecular quantum mechanics , 2014 .

[36]  W. C. Lineberger,et al.  Laser photoelectron spectrometry of PO - , PH - , and PH - 2 , 1976 .

[37]  Wenjian Liu,et al.  Theoretical and numerical assessments of spin-flip time-dependent density functional theory. , 2012, The Journal of chemical physics.

[38]  Marvin Douglas,et al.  Quantum electrodynamical corrections to the fine structure of helium , 1971 .

[39]  Wenjian Liu,et al.  Critical Assessment of Time-Dependent Density Functional Theory for Excited States of Open-Shell Systems: II. Doublet-Quartet Transitions. , 2016, Journal of chemical theory and computation.

[40]  U. Sassenberg,et al.  The Vacuum Ultraviolet Spectrum of Selenium Hydride I. Determination of the Ground State Spin-Orbit Coupling Constant , 1978 .

[41]  S. Knecht,et al.  Large-scale parallel configuration interaction. II. Two- and four-component double-group general active space implementation with application to BiH. , 2010, The Journal of chemical physics.

[42]  Hosung Sun,et al.  Spin–orbit splittings in the valence states of XH (X = K, Ca, Ga, Ge, As, Se, and Br) by the effective Hamiltonian approach , 2010 .

[43]  Yunlong Xiao,et al.  Combining spin-adapted open-shell TD-DFT with spin–orbit coupling , 2013 .

[44]  Wenjian Liu.  Big picture of relativistic molecular quantum mechanics , 2016 .

[45]  Lemin Li,et al.  RELATIVISTIC DENSITY FUNCTIONAL THEORY: THE BDF PROGRAM PACKAGE , 2004 .

[46]  Wenjian Liu,et al.  On the spin separation of algebraic two-component relativistic Hamiltonians: molecular properties. , 2014, The Journal of chemical physics.

[47]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[48]  Fan Wang,et al.  Equation-of-Motion Coupled-Cluster Theory for Excitation Energies of Closed-Shell Systems with Spin-Orbit Coupling. , 2014, Journal of chemical theory and computation.

[49]  W. C. Lineberger,et al.  Photoelectron spectroscopy of phosphorus hydride anions. , 2005, The Journal of chemical physics.

[50]  D. Fedorov,et al.  Performance of dynamically weighted multiconfiguration self-consistent field and spin-orbit coupling calculations of diatomic molecules of Group 14 elements. , 2011, The Journal of chemical physics.

[51]  Roland Lindh,et al.  New relativistic ANO basis sets for transition metal atoms. , 2005, The journal of physical chemistry. A.

[52]  Lan Cheng,et al.  Making four- and two-component relativistic density functional methods fully equivalent based on the idea of "from atoms to molecule". , 2007, The Journal of chemical physics.

[53]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[54]  R. Bartlett,et al.  Property evaluation and orbital relaxation in coupled cluster methods , 1987 .

[55]  D. Fedorov,et al.  Spin-orbit multireference multistate perturbation theory , 2001 .

[56]  M. Pélissier,et al.  Ab initio molecular calculations including spin-orbit coupling. I. Method and atomic tests , 1983 .

[57]  Anna I Krylov,et al.  Equation-of-motion coupled-cluster methods for open-shell and electronically excited species: the Hitchhiker's guide to Fock space. , 2008, Annual review of physical chemistry.

[58]  Michael Dolg,et al.  Relativistic energy‐consistent pseudopotentials—Recent developments , 2002, J. Comput. Chem..

[59]  Junji Seino,et al.  Expectation values in two-component relativistic theories. , 2010, The Journal of chemical physics.

[60]  Bingbing Suo,et al.  Spin-adapted open-shell time-dependent density functional theory. II. Theory and pilot application. , 2011, The Journal of chemical physics.

[61]  Fan Wang,et al.  Equation-of-motion coupled-cluster method for doubly ionized states with spin-orbit coupling. , 2012, The Journal of chemical physics.

[62]  Daoling Peng,et al.  Infinite-order quasirelativistic density functional method based on the exact matrix quasirelativistic theory. , 2006, The Journal of chemical physics.

[63]  Wenjian Liu.  Ideas of relativistic quantum chemistry , 2010 .

[64]  Yoon Sup Lee,et al.  Effects of the magnetic part of the Breit term on the 2Π states of diatomic hydrides , 1994 .

[65]  C. Marian,et al.  SPOCK.CI: a multireference spin-orbit configuration interaction method for large molecules. , 2006, The Journal of chemical physics.