Spin-adapted open-shell random phase approximation and time-dependent density functional theory. I. Theory.

The spin-adaptation of single-reference quantum chemical methods for excited states of open-shell systems has been nontrivial. The primary reason is that the configuration space, generated by a truncated rank of excitations from only one component of a reference multiplet, is spin-incomplete. Those "missing" configurations are of higher ranks and can, in principle, be recaptured by a particular class of excitation operators. However, the resulting formalisms are then quite involved and there are situations [e.g., time-dependent density functional theory (TD-DFT) under the adiabatic approximation] that prevent one from doing so. To solve this issue, we propose here a tensor-coupling scheme that invokes all the components of a reference multiplet (i.e., a tensor reference) rather than increases the excitation ranks. A minimal spin-adapted n-tuply excited configuration space can readily be constructed by tensor products between the n-tuple tensor excitation operators and the chosen tensor reference. Further combined with the tensor equation-of-motion formalism, very compact expressions for excitation energies can be obtained. As a first application of this general idea, a spin-adapted open-shell random phase approximation is first developed. The so-called "translation rule" is then adopted to formulate a spin-adapted, restricted open-shell Kohn-Sham (ROKS)-based TD-DFT (ROKS-TD-DFT). Here, a particular symmetry structure has to be imposed on the exchange-correlation kernel. While the standard ROKS-TD-DFT can access only excited states due to singlet-coupled single excitations, i.e., only some of the singly excited states of the same spin (S(i)) as the reference, the new scheme can capture all the excited states of spin S(i)-1, S(i), or S(i)+1 due to both singlet- and triplet-coupled single excitations. The actual implementation and computation are very much like the (spin-contaminated) unrestricted Kohn-Sham-based TD-DFT. It is also shown that spin-contaminated spin-flip configuration interaction approaches can easily be spin-adapted via the tensor-coupling scheme.

[1]  Debashis Mukherjee,et al.  A response-function approach to the direct calculation of the transition-energy in a multiple-cluster expansion formalism , 1979 .

[2]  D. Rowe,et al.  The tensor open-shell random phase approximation with application to the even nickel isotopes☆ , 1979 .

[3]  Ž. Rinkevičius,et al.  General excitations in time-dependent density functional theory. , 2007, The Journal of chemical physics.

[4]  Zilvinas Rinkevicius,et al.  Density functional theory for hyperfine coupling constants with the restricted-unrestricted approach. , 2004, The Journal of chemical physics.

[5]  Wenli Zou,et al.  Time-dependent quasirelativistic density-functional theory based on the zeroth-order regular approximation. , 2005, The Journal of chemical physics.

[6]  Società italiana di fisica,et al.  Struttura e reazioni nucleari nella teoria dei molti corpi = Many-body description of nuclear structure and reactions , 1966 .

[7]  Anna I Krylov,et al.  Equation-of-motion spin-flip coupled-cluster model with single and double substitutions: Theory and application to cyclobutadiene. , 2004, The Journal of chemical physics.

[8]  Anna I. Krylov,et al.  The spin–flip approach within time-dependent density functional theory: Theory and applications to diradicals , 2003 .

[9]  D R Yarkony,et al.  Modern electronic structure theory , 1995 .

[10]  E. Gross,et al.  Density-Functional Theory for Time-Dependent Systems , 1984 .

[11]  Anna I Krylov,et al.  Spin-flip equation-of-motion coupled-cluster electronic structure method for a description of excited states, bond breaking, diradicals, and triradicals. , 2006, Accounts of chemical research.

[12]  J. Gauss,et al.  Towards a spin-adapted coupled-cluster theory for high-spin open-shell states. , 2006, The Journal of chemical physics.

[13]  M. E. Casida,et al.  Propagator corrections to adiabatic time-dependent density-functional theory linear response theory. , 2005, The Journal of chemical physics.

[14]  So Hirata,et al.  Time-dependent density functional theory for radicals: An improved description of excited states with substantial double excitation character , 1999 .

[15]  N. Nakatsuji,et al.  Cluster expansion of the wavefunction. Excited states , 1978 .

[16]  N. Handy,et al.  Potential energy curves for PO, calculated using DFT and MRCI methodology , 1999 .

[17]  Tom Ziegler,et al.  Time-dependent density functional theory based on a noncollinear formulation of the exchange-correlation potential. , 2004, The Journal of chemical physics.

[18]  Jun Shen,et al.  Block correlated coupled cluster method with a complete-active-space self-consistent-field reference function: the implementation for low-lying excited states. , 2008, The Journal of chemical physics.

[19]  Henrik Koch,et al.  Coupled cluster response functions , 1990 .

[20]  D. Chong Recent Advances in Density Functional Methods Part III , 2002 .

[21]  Josef Paldus,et al.  Automation of the implementation of spin‐adapted open‐shell coupled‐cluster theories relying on the unitary group formalism , 1994 .

[22]  J. Nichols,et al.  Multiconfigurational electron propagator (MCEP) ionization potentials for general open shell systems , 1984 .

[23]  P. Wormer,et al.  Frequency‐dependent polarizabilities and van der Waals coefficients of half‐open‐shell systems in the time‐dependent coupled Hartree–Fock approximation , 1990 .

[24]  Wenli Zou,et al.  Time-dependent four-component relativistic density-functional theory for excitation energies. II. The exchange-correlation kernel. , 2005, The Journal of chemical physics.

[25]  D. Salahub,et al.  Time-dependent density-functional theory investigation of excitation spectra of open-shell molecules , 2000 .

[26]  E. Gross,et al.  Time-dependent density functional theory. , 2004, Annual review of physical chemistry.

[27]  D. Rowe,et al.  Tensor equations of motion for the excitations of rotationally invariant or charge-independent systems , 1975 .

[28]  P. Joergensen,et al.  Second Quantization-based Methods in Quantum Chemistry , 1981 .

[29]  Wenjian Liu,et al.  Comparison of Different Polarization Schemes in Open‐shell Relativistic Density Functional Calculations , 2003 .

[30]  H. Schaefer Methods of Electronic Structure Theory , 1977 .

[31]  H. Ågren,et al.  Spin-flip time dependent density functional theory for singlet–triplet splittings in σ,σ-biradicals , 2010 .

[32]  Tom Ziegler,et al.  The performance of time-dependent density functional theory based on a noncollinear exchange-correlation potential in the calculations of excitation energies. , 2005, The Journal of chemical physics.

[33]  V. McKoy,et al.  An equations of motion approach for open shell systems , 1975 .

[34]  M. Head‐Gordon,et al.  Configuration interaction with single substitutions for excited states of open-shell molecules , 1995 .

[35]  S. Pal,et al.  Use of Cluster Expansion Methods in the Open-Shell Correlation Problem , 1989 .

[36]  Rodney J. Bartlett,et al.  To Multireference or not to Multireference: That is the Question? , 2002 .

[37]  David Casanova,et al.  The spin-flip extended single excitation configuration interaction method. , 2008, The Journal of chemical physics.

[38]  F. Furche Developing the random phase approximation into a practical post-Kohn-Sham correlation model. , 2008, The Journal of chemical physics.

[39]  D. Yeager,et al.  The multiconfigurational spin‐tensor electron propagator method for determining vertical principal and shake‐up ionization potentials for open shell and highly correlated atoms and molecules , 1987 .

[40]  Chengbu Liu,et al.  Time-dependent four-component relativistic density functional theory for excitation energies. , 2004, The Journal of chemical physics.

[41]  Paweł Sałek,et al.  Restricted density functional theory of linear time-dependent properties in open-shell molecules , 2003 .

[42]  David J Rowe,et al.  EQUATIONS-OF-MOTION METHOD AND THE EXTENDED SHELL MODEL. , 1968 .

[43]  J. Olsen,et al.  Generalizations of the multiconfigurational time-dependent Hartree–Fock approach , 1984 .

[44]  Kimihiko Hirao,et al.  Multireference Møller—Plesset perturbation theory for high-spin open-shell systems , 1992 .

[45]  Georg Kresse,et al.  Hybrid functionals including random phase approximation correlation and second-order screened exchange. , 2010, The Journal of chemical physics.

[46]  B. Roos,et al.  A complete active space SCF method (CASSCF) using a density matrix formulated super-CI approach , 1980 .

[47]  Francesco A Evangelista,et al.  High-order excitations in state-universal and state-specific multireference coupled cluster theories: model systems. , 2006, The Journal of chemical physics.

[48]  Fan Zhang,et al.  A dressed TDDFT treatment of the 21Ag states of butadiene and hexatriene , 2004 .

[49]  K. Hirao,et al.  Cluster expansion of the wavefunction. Pseudo‐orbital theory based on the SAC expansion and its application to the spin density of open‐shell systems , 1978 .

[50]  Poul Jørgensen,et al.  The second-order approximate coupled cluster singles and doubles model CC2 , 1995 .

[51]  M. Seth,et al.  Calculation of excitation energies of open-shell molecules with spatially degenerate ground states. I. Transformed reference via an intermediate configuration Kohn-Sham density-functional theory and applications to d1 and d2 systems with octahedral and tetrahedral symmetries. , 2005, The Journal of chemical physics.

[52]  Mark E. Casida,et al.  Excited-state spin-contamination in time-dependent density-functional theory for molecules with open-shell ground states , 2009 .

[53]  Curtis L. Janssen,et al.  The automated solution of second quantization equations with applications to the coupled cluster approach , 1991 .

[54]  Anna I. Krylov,et al.  A spin-complete version of the spin-flip approach to bond breaking: What is the impact of obtaining spin eigenfunctions? , 2003 .

[55]  Poul Jo,et al.  Transition moments and dynamic polarizabilities in a second order polarization propagator approach , 1980 .

[56]  Grzegorz Mazur,et al.  Application of the dressed time‐dependent density functional theory for the excited states of linear polyenes , 2009, J. Comput. Chem..

[57]  John F. Stanton,et al.  The equation of motion coupled‐cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties , 1993 .

[58]  S. Chattopadhyay,et al.  State-specific multi-reference perturbation theories with relaxed coefficients: molecular applications , 2002 .

[59]  M. Head‐Gordon,et al.  On the Nature of Electronic Transitions in Radicals: An Extended Single Excitation Configuration Interaction Method , 1996 .

[60]  Anna I. Krylov,et al.  Size-consistent wave functions for bond-breaking: the equation-of-motion spin-flip model , 2001 .

[61]  V. McKoy,et al.  HIGHER RANDOM-PHASE APPROXIMATION AS AN APPROXIMATION TO THE EQUATIONS OF MOTION. , 1970 .

[62]  Zilvinas Rinkevicius,et al.  Time-dependent density functional theory with the generalized restricted-unrestricted approach. , 2006, The Journal of chemical physics.

[63]  Anna I Krylov,et al.  Double spin-flip approach within equation-of-motion coupled cluster and configuration interaction formalisms: Theory, implementation, and examples. , 2009, The Journal of chemical physics.

[64]  Michael W. Schmidt,et al.  Are atoms intrinsic to molecular electronic wavefunctions? I. The FORS model , 1982 .

[65]  V. Chernyak,et al.  Resonant nonlinear polarizabilities in the time-dependent density functional theory , 2003 .

[66]  J. Pople,et al.  Self‐Consistent Molecular Orbital Methods. X. Molecular Orbital Studies of Excited States with Minimal and Extended Basis Sets , 1971 .

[67]  Anna I. Krylov,et al.  Spin-flip configuration interaction: an electronic structure model that is both variational and size-consistent , 2001 .

[68]  P. Jørgensen,et al.  Polarization propagator methods in atomic and molecular calculations , 1984 .

[69]  Kimihiko Hirao,et al.  Multireference Møller-Plesset method , 1992 .

[70]  B. Jeziorski,et al.  Dispersion interaction of high-spin open-shell complexes in the random phase approximation , 2003 .

[71]  Kimihiko Hirao,et al.  Cluster expansion of the wavefunction. Symmetry-adapted-cluster expansion, its variational determination, and extension of open-shell orbital theory , 1978 .

[72]  Martin J. Packer,et al.  A new implementation of the second‐order polarization propagator approximation (SOPPA): The excitation spectra of benzene and naphthalene , 1996 .

[73]  M. E. Casida Time-Dependent Density Functional Response Theory for Molecules , 1995 .

[74]  Henrik Koch,et al.  Calculation of frequency-dependent polarizabilities using coupled-cluster response theory , 1994 .

[75]  M. Petersilka,et al.  Excitation energies from time-dependent density-functional theory. , 1996 .

[76]  Trygve Helgaker,et al.  Hartree-Fock and Kohn-Sham time-dependent response theory in a second-quantization atomic-orbital formalism suitable for linear scaling. , 2008, The Journal of chemical physics.

[77]  N. N. Breslavskaya,et al.  Koopmans' theorem in the ROHF method: canonical form for the Hartree-Fock Hamiltonian. , 2006, The Journal of chemical physics.

[78]  Donald C. Comeau,et al.  The equation-of-motion coupled-cluster method. Applications to open- and closed-shell reference states , 1993 .

[79]  Kieron Burke,et al.  Double excitations within time-dependent density functional theory linear response. , 2004, The Journal of chemical physics.

[80]  D Mukherjee,et al.  Molecular Applications of a Size-Consistent State-Specific Multireference Perturbation Theory with Relaxed Model-Space Coefficients. , 1999, The journal of physical chemistry. A.

[81]  Somnath Ghosh,et al.  A spin-adapted linear response theory in a coupled-cluster framework for direct calculation of spin-allowed and spin-forbidden transition energies , 1982 .