Perturbational treatment of spin-orbit coupling for generally applicable high-level multi-reference methods.

An efficient perturbational treatment of spin-orbit coupling within the framework of high-level multi-reference techniques has been implemented in the most recent version of the Columbus quantum chemistry package, extending the existing fully variational two-component (2c) multi-reference configuration interaction singles and doubles (MRCISD) method. The proposed scheme follows related implementations of quasi-degenerate perturbation theory (QDPT) model space techniques. Our model space is built either from uncontracted, large-scale scalar relativistic MRCISD wavefunctions or based on the scalar-relativistic solutions of the linear-response-theory-based multi-configurational averaged quadratic coupled cluster method (LRT-MRAQCC). The latter approach allows for a consistent, approximatively size-consistent and size-extensive treatment of spin-orbit coupling. The approach is described in detail and compared to a number of related techniques. The inherent accuracy of the QDPT approach is validated by comparing cuts of the potential energy surfaces of acrolein and its S, Se, and Te analoga with the corresponding data obtained from matching fully variational spin-orbit MRCISD calculations. The conceptual availability of approximate analytic gradients with respect to geometrical displacements is an attractive feature of the 2c-QDPT-MRCISD and 2c-QDPT-LRT-MRAQCC methods for structure optimization and ab inito molecular dynamics simulations.

[1]  Thomas Müller,et al.  Columbus—a program system for advanced multireference theory calculations , 2011 .

[2]  T. Fleig.  Invited review: Relativistic wave-function based electron correlation methods , 2012 .

[3]  B. Roos,et al.  Relativistic quantum chemistry: the multiconfigurational approach , 2004 .

[4]  E. Engel.  Relativistic Density Functional Theory: Foundations and Basic Formalism , 2002 .

[5]  H. Lischka,et al.  Excitation energies and transition moments by the multireference averaged quadratic coupled cluster (MR-AQCC) method , 2000 .

[6]  Evert Jan Baerends,et al.  Relativistic regular two-component Hamiltonians. , 1996 .

[7]  Rodney J. Bartlett,et al.  Multi-reference averaged quadratic coupled-cluster method: a size-extensive modification of multi-reference CI , 1993 .

[8]  Valérie Vallet,et al.  A two-step uncontracted determinantal effective Hamiltonian-based SO–CI method , 2000 .

[9]  G. Breit.  The Effect of Retardation on the Interaction of Two Electrons , 1929 .

[10]  Mostafa A. El-Sayed,et al.  Spin—Orbit Coupling and the Radiationless Processes in Nitrogen Heterocyclics , 1963 .

[11]  L. Seijo,et al.  The ab initio model potential method. Main group elements , 1987 .

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

[13]  G. Worth,et al.  Competing ultrafast intersystem crossing and internal conversion in the "channel 3" region of benzene. , 2010, Physical chemistry chemical physics : PCCP.

[14]  Frank Neese,et al.  The ORCA program system , 2012 .

[15]  W. Schwarz.  An Introduction to Relativistic Quantum Chemistry , 2010 .

[16]  P. L. de Boeij,et al.  Relativistic two-component formulation of time-dependent current-density functional theory: application to the linear response of solids. , 2007, The Journal of chemical physics.

[17]  Frank Neese,et al.  Advanced aspects of ab initio theoretical optical spectroscopy of transition metal complexes: Multiplets, spin-orbit coupling and resonance Raman intensities , 2007 .

[18]  J. Gauss,et al.  Spin-orbit coupling constants from coupled-cluster response theory , 2000 .

[19]  G. Worth,et al.  Quantum dynamics study of the competing ultrafast intersystem crossing and internal conversion in the "channel 3" region of benzene. , 2012, The Journal of chemical physics.

[20]  Trygve Helgaker,et al.  Principles of direct 4-component relativistic SCF: application to caesium auride , 1997 .

[21]  Thomas Müller,et al.  High-level multireference methods in the quantum-chemistry program system COLUMBUS: Analytic MR-CISD and MR-AQCC gradients and MR-AQCC-LRT for excited states, GUGA spin–orbit CI and parallel CI density , 2001 .

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

[23]  Russell M. Pitzer,et al.  Spin‐orbit (core) and core potential integrals , 1991 .

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

[25]  C. Crespo-Hernández,et al.  Environmental photochemistry of nitro-PAHs: direct observation of ultrafast intersystem crossing in 1-nitropyrene. , 2008, The journal of physical chemistry. A.

[26]  H. Lee,et al.  Spin-orbit effects calculated by two-component coupled-cluster methods: test calculations on AuH, Au2, TlH and Tl2 , 1998 .

[27]  B. A. Hess,et al.  Spin—orbit effects in the Br atom in the framework of the no-pair theory , 1991 .

[28]  C. Crespo-Hernández,et al.  Ultrafast excited-state dynamics in nucleic acids. , 2004, Chemical reviews.

[29]  K. Hirao,et al.  The Douglas-Kroll-Hess approach. , 2012, Chemical reviews.

[30]  Thomas Müller,et al.  Large-scale parallel uncontracted multireference-averaged quadratic coupled cluster: the ground state of the chromium dimer revisited. , 2009, The journal of physical chemistry. A.

[31]  D. Yarkony,et al.  On the effects of spin–orbit coupling on conical intersection seams in molecules with an odd number of electrons. II. Characterizing the local topography of the seam , 2001 .

[32]  C. Marian,et al.  Four-component relativistic coupled cluster and configuration interaction calculations on the ground and excited states of the RbYb molecule. , 2009, The journal of physical chemistry. A.

[33]  Michael Dolg,et al.  Relativistic pseudopotentials: their development and scope of applications. , 2012, Chemical reviews.

[34]  Per E. M. Siegbahn,et al.  The Direct Configuration Interaction Method from Molecular Integrals , 1977 .

[35]  F. Neese,et al.  A fully variational spin-orbit coupled complete active space self-consistent field approach: application to electron paramagnetic resonance g-tensors. , 2013, The Journal of chemical physics.

[36]  P. Marquetand,et al.  SHARC: ab Initio Molecular Dynamics with Surface Hopping in the Adiabatic Representation Including Arbitrary Couplings. , 2011, Journal of chemical theory and computation.

[37]  D. Yarkony,et al.  Spin-orbit coupling and conical intersections in molecules with an odd number of electrons. III. A perturbative determination of the electronic energies, derivative couplings and a rigorous diabatic representation near a conical intersection , 2002 .

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

[39]  P. Marquetand,et al.  Singlet and triplet excited-state dynamics study of the keto and enol tautomers of cytosine. , 2013, Chemphyschem : a European journal of chemical physics and physical chemistry.

[40]  A. B. Alekseyev,et al.  Comparison of spin-orbit configuration interaction methods employing relativistic effective core potentials for the calculation of zero-field splittings of heavy atoms with a 2Po ground state , 1998 .

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

[42]  H. Lischka,et al.  Multiconfiguration self-consistent field and multireference configuration interaction methods and applications. , 2012, Chemical reviews.

[43]  Trygve Helgaker,et al.  Four‐component relativistic Kohn–Sham theory , 2002, J. Comput. Chem..

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

[45]  P. A. Christiansen,et al.  IMPROVED Ab Initio EFFECTIVE CORE POTENTIALS FOR MOLECULAR CALCULATIONS , 1979 .

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

[47]  W. Jr. Pauli,et al.  Zur Quantenmechanik des magnetischen Elektrons , 1927 .

[48]  Bernd Schimmelpfennig,et al.  The restricted active space (RAS) state interaction approach with spin-orbit coupling , 2002 .

[49]  Isaiah Shavitt,et al.  The Graphical Unitary Group Approach and Its Application to Direct Configuration Interaction Calculations , 1981 .

[50]  K. Dyall,et al.  Introduction to Relativistic Quantum Chemistry , 2007 .

[51]  P. Marquetand,et al.  Femtosecond Intersystem Crossing in the DNA Nucleobase Cytosine. , 2012, The journal of physical chemistry letters.

[52]  H. Quiney,et al.  Chapter 3 - Relativistic Self-Consistent Fields , 2002 .

[53]  D. Yarkony,et al.  On the effects of spin-orbit coupling on conical intersection seams in molecules with an odd number of electrons. I. Locating the seam , 2001 .

[54]  Markus Reiher,et al.  The generalized Douglas–Kroll transformation , 2002 .

[55]  Markus Reiher,et al.  Douglas–Kroll–Hess Theory: a relativistic electrons-only theory for chemistry , 2006 .

[56]  G. Granucci,et al.  Competing ultrafast intersystem crossing and internal conversion: a time resolved picture for the deactivation of 6-thioguanine , 2014 .

[57]  C. Mead.  The ''noncrossing'' rule for electronic potential energy surfaces: The role of time-reversal invariance , 1979 .

[58]  P. Dirac.  The quantum theory of the electron , 1928 .

[59]  Fernando Bernardi,et al.  Excited-State Potential Surface Crossings in Acrolein - A Model for Understanding the Photochemistry and Photophysics of Alpha,Beta-Enones , 1994 .

[60]  Pekka Pyykkö,et al.  Relativistic effects in structural chemistry , 1988 .

[61]  Martin Klessinger,et al.  Excited states and photochemistry of organic molecules , 1995 .

[62]  W. C. Ermler,et al.  Spin-Orbit Coupling and Other Relativistic Effects in Atoms and Molecules , 1988 .

[63]  FRANCESCO AQUILANTE,et al.  MOLCAS 7: The Next Generation , 2010, J. Comput. Chem..

[64]  M. Klobukowski,et al.  Model Core Potentials: Theory and Applications , 1999 .

[65]  J. Olsen,et al.  The generalized active space concept for the relativistic treatment of electron correlation. I. Kramers-restricted two-component configuration interaction , 2001 .

[66]  A. Messiah.  Quantum Mechanics , 1961 .

[67]  Tom Ziegler,et al.  A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect. , 2005, The Journal of chemical physics.

[68]  Derivation of spin-orbit couplings in collinear linear-response TDDFT: a rigorous formulation. , 2014, The Journal of chemical physics.

[69]  Wim Klopper,et al.  Self-consistent treatment of spin-orbit interactions with efficient Hartree-Fock and density functional methods. , 2008, Physical chemistry chemical physics : PCCP.

[70]  E. Eliav,et al.  Relativistic Four-Component Multireference Coupled Cluster Methods: Towards A Covariant Approach , 2010 .

[71]  A. Walsh.  The absorption spectra of acrolein, crotonaldehyde and mesityl oxide in the vacuum ultra-violet , 1945 .

[72]  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.

[73]  P. A. Christiansen,et al.  Relativistic ab initio molecular structure calculations including configuration interaction with application to six states of TlH , 1982 .

[74]  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.

[75]  Philipp Marquetand,et al.  Non-adiabatic dynamics in SO2: II. The role of triplet states studied by surface-hopping simulations , 2013 .

[76]  C. Crespo-Hernández,et al.  Excited-state dynamics in nitro-naphthalene derivatives: intersystem crossing to the triplet manifold in hundreds of femtoseconds. , 2013, The journal of physical chemistry. A.