Performance of TD-DFT for Excited States of Open-Shell Transition Metal Compounds.

Time-dependent density functional theory (TD-DFT) has been very successful in accessing low-lying excited states of closed-shell systems. However, it is much less so for excited states of open-shell systems: unrestricted Kohn-Sham based TD-DFT (U-TD-DFT) often produces physically meaningless excited states due to heavy spin contaminations, whereas restricted Kohn-Sham based TD-DFT often misses those states of lower energies. A much better variant is the explicitly spin-adapted TD-DFT (X-TD-DFT) [J. Chem. Phys. 2011, 135, 194106] that can capture all the spin-adapted singly excited states yet without computational overhead over U-TD-DFT. While the superiority of X-TD-DFT over U-TD-DFT has been demonstrated for open-shell systems of main group elements, it remains to be seen if this is also the case for open-shell transition metal compounds. Taking as benchmark the results by MS-CASPT2 (multistate complete active space second-order perturbation theory) and ic-MRCISD (internally contracted multireference configuration interaction with singles and doubles), it is shown that X-TD-DFT is indeed superior to U-TD-DFT for the vertical excitation energies of ZnH, CdH, ScH2, YH2, YO, and NbO2. Admittedly, there exist a few cases where U-TD-DFT appears to be better than X-TD-DFT. However, this is due to a wrong reason: the underestimation (due to spin contamination) and the overestimation (due to either the exchange-correlation functional itself or the adiabatic approximation to the exchange-correlation kernel) happen to be compensated in the case of U-TD-DFT. As for [Cu(C6H6)2]2+, which goes beyond the capability of both MS-CASPT2 and ic-MRCISD, X-TD-DFT revises the U-TD-DFT assignment of the experimental spectrum.

[1]  Wenjian Liu,et al.  Erratum: “Spin-adapted open-shell time-dependent density functional theory. III. An even better and simpler formulation” [J. Chem. Phys. 135, 194106 (2011)] , 2013 .

[2]  Björn O. Roos,et al.  Second-order perturbation theory with a complete active space self-consistent field reference function , 1992 .

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

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

[5]  K. Evenson,et al.  Far-Infrared Rotational Spectra of ZnH and ZnD in theX2Σ+(v= 0) State , 1997 .

[6]  Wenjian Liu,et al.  First-order nonadiabatic coupling matrix elements between excited states: a Lagrangian formulation at the CIS, RPA, TD-HF, and TD-DFT levels. , 2014, The Journal of chemical physics.

[7]  Xuefeng Wang,et al.  Matrix Infrared Spectra and DFT Calculations of the Reactive MHx (x = 1, 2, and 3), (H2)MH2, MH2+, and MH4- (M = Sc, Y, and La) Species , 2002 .

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

[9]  P. Knowles,et al.  An efficient internally contracted multiconfiguration–reference configuration interaction method , 1988 .

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

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

[12]  Hans-Joachim Werner,et al.  Multireference perturbation theory for large restricted and selected active space reference wave functions , 2000 .

[13]  Celestino Angeli,et al.  A quasidegenerate formulation of the second order n-electron valence state perturbation theory approach. , 2004, The Journal of chemical physics.

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

[15]  Mark R. Hoffmann,et al.  SDS: the ‘static–dynamic–static’ framework for strongly correlated electrons , 2014, Theoretical Chemistry Accounts.

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

[17]  I. Ciofini,et al.  Assessment of Functionals for TD-DFT Calculations of Singlet-Triplet Transitions. , 2010, Journal of chemical theory and computation.

[18]  K. Balasubramanian Electronic states and potential energy surfaces of ScH2 , 1987 .

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

[20]  M. Hoffmann Canonical Van Vleck Quasidegenerate Perturbation Theory with Trigonometric Variables , 1996 .

[21]  Kirk A Peterson,et al.  Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc-Zn. , 2005, The Journal of chemical physics.

[22]  C. Jamorski,et al.  Theoretical determination of spectral lines for the Zn atom and the ZnH molecule , 1994 .

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

[24]  Zhengting Gan,et al.  Hole–particle correspondence in CI calculations , 2003 .

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

[26]  K. Hirao,et al.  A long-range correction scheme for generalized-gradient-approximation exchange functionals , 2001 .

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

[28]  M. Hoffmann,et al.  Explication and revision of generalized Van Vleck perturbation theory for molecular electronic structure , 2002 .

[29]  Celestino Angeli,et al.  Introduction of n-electron valence states for multireference perturbation theory , 2001 .

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

[31]  Mark R. Hoffmann,et al.  Further development of SDSPT2 for strongly correlated electrons , 2017 .

[32]  Angela K. Wilson,et al.  Gaussian basis sets for use in correlated molecular calculations. IX. The atoms gallium through krypton , 1993 .

[33]  Luis Serrano-Andrés,et al.  The multi-state CASPT2 method , 1998 .

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

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

[36]  C. David Sherrill,et al.  The Configuration Interaction Method: Advances in Highly Correlated Approaches , 1999 .

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

[38]  Evert Jan Baerends,et al.  Molecular calculations of excitation energies and (hyper)polarizabilities with a statistical average of orbital model exchange-correlation potentials , 2000 .

[39]  Wenjian Liu,et al.  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. , 2017, Physical chemistry chemical physics : PCCP.

[40]  Gaohong Zhai,et al.  Doubly contracted CI method , 2004 .

[41]  Yubin Wang,et al.  New schemes for internally contracted multi-reference configuration interaction. , 2014, The Journal of chemical physics.

[42]  B. Suo,et al.  High-resolution electronic spectra of yttrium oxide (YO): The D2Σ+-X2Σ+ transition. , 2017, The Journal of chemical physics.

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

[44]  Uttam Sinha Mahapatra,et al.  Development of a size-consistent state-specific multireference perturbation theory with relaxed model-space coefficients , 1999 .

[45]  C. Linton Photoluminescence of the A2Π-X2Σ+ system of the yttrium oxide molecule , 1978 .

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

[47]  Jochen Schirmer,et al.  Beyond the random-phase approximation: A new approximation scheme for the polarization propagator , 1982 .

[48]  S. Iwata,et al.  Theoretical studies of the ground and low-lying excited electronic states of the early transition metal dihydrides with state averaged MC SCF method , 1996 .

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

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

[51]  Benjamin T. Miller,et al.  A parallel implementation of the analytic nuclear gradient for time-dependent density functional theory within the Tamm–Dancoff approximation , 1999 .

[52]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

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

[54]  Rui Li,et al.  Spin–orbit all-electron configuration interaction study on the electronic structure and radiative lifetimes of low-lying excited states of CdH , 2014 .

[55]  K. Balasubramanian,et al.  Electronic states and potential energy surfaces of YH2 , 1988 .

[56]  Vincenzo Barone,et al.  Accurate excitation energies from time-dependent density functional theory: Assessing the PBE0 model , 1999 .

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

[58]  Wenjian Liu,et al.  A spin-adapted size-extensive state-specific multi-reference perturbation theory with various partitioning schemes. II. Molecular applications. , 2012, The Journal of chemical physics.

[59]  F. E. Jorge,et al.  Augmented Gaussian basis sets of triple and quadruple zeta valence quality for the atoms H and from Li to Ar: Applications in HF, MP2, and DFT calculations of molecular dipole moment and dipole (hyper)polarizability , 2007 .

[60]  W. T. Fernando,et al.  Fourier transform emission spectroscopy of the A2Π-X2Σ+ transition of ZnD , 1990 .

[61]  Wenjian Liu,et al.  A spin-adapted size-extensive state-specific multi-reference perturbation theory. I. Formal developments. , 2012, The Journal of chemical physics.

[62]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[63]  P. Bernath,et al.  Multi-isotopologue analyses of new vibration–rotation and pure rotation spectra of ZnH and CdH , 2006 .

[64]  Toru Shiozaki,et al.  Communication: extended multi-state complete active space second-order perturbation theory: energy and nuclear gradients. , 2011, The Journal of chemical physics.

[65]  Wenjian Liu,et al.  iCI: Iterative CI toward full CI. , 2016, Journal of chemical theory and computation.

[66]  T. Steimle,et al.  Fine and hyperfine structure in the X2Σ+ state of gas-phase yttrium monoxide , 1987 .

[67]  J. Gole,et al.  Single collision chemiluminescence studies of scandium and yttrium oxidation with O2, NO2, N2O and O3 , 1977 .

[68]  N. Handy,et al.  A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP) , 2004 .

[69]  Hans-Joachim Werner,et al.  A new internally contracted multi-reference configuration interaction method. , 2011, The Journal of chemical physics.

[70]  M. Head‐Gordon,et al.  Time-dependent density functional study on the electronic excitation energies of polycyclic aromatic hydrocarbon radical cations of naphthalene, anthracene, pyrene, and perylene , 1999 .

[71]  B. Roos,et al.  A modified definition of the zeroth-order Hamiltonian in multiconfigurational perturbation theory (CASPT2) , 2004 .

[72]  Bingbing Suo,et al.  First order nonadiabatic coupling matrix elements between excited states: implementation and application at the TD-DFT and pp-TDA levels. , 2014, The Journal of chemical physics.

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

[74]  H. Cox,et al.  Gas phase UV spectrum of a Cu(II)-bis(benzene) sandwich complex: experiment and theory. , 2014, The journal of physical chemistry. A.

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

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

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

[78]  Lemin Li,et al.  Parallelization of MRCI based on hole‐particle symmetry , 2005, J. Comput. Chem..

[79]  C. T. Campos,et al.  Augmented Gaussian basis sets for the elements K, Sc–Kr, Rb, and Y–Xe: Application in HF, MP2, and DFT calculations of molecular electric properties , 2013 .

[80]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[81]  Haruyuki Nakano,et al.  Quasidegenerate perturbation theory with multiconfigurational self‐consistent‐field reference functions , 1993 .

[82]  Andreas Dreuw,et al.  The algebraic diagrammatic construction scheme for the polarization propagator for the calculation of excited states , 2015 .

[83]  J. Gole,et al.  Chemiluminescence spectra of ScO and YO: Observation and analysis of the A′ 2Δ–X 2Σ+ band system , 1976 .

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

[85]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[86]  Hongbin Wu,et al.  Photoelectron Spectroscopy and Electronic Structure of ScOn- (n = 1−4) and YOn- (n = 1−5): Strong Electron Correlation Effects in ScO-and YO- , 1998 .

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

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

[89]  Carlo Adamo,et al.  Extensive TD-DFT Benchmark: Singlet-Excited States of Organic Molecules. , 2009, Journal of chemical theory and computation.