Beyond the Coulson–Fischer point: characterizing single excitation CI and TDDFT for excited states in single bond dissociations. Issue 39 (25th September 2019)

• Record Type: Journal Article
• Authors: Hait, Diptarka; Rettig, Adam; Head-Gordon, Martin
• Is Part Of: Physical chemistry chemical physics. Volume 21:Issue 39(2019)
• Page Start: 21761
• Page End: 21775
• Publication Date: 2019-09-25
• Subjects: Chemistry, Physical and theoretical -- Periodicals; 541.3
• DOI: 10.1039/c9cp04452c
• Languages: English
• ISSNs: 1463-9076

Abstract

Abstract : HF/DFT orbitals spin-polarize when single bonds are stretched past the Coulson–Fischer point. We report unphysical features in the excited state potential energy surfaces predicted by CIS/TDDFT in this regime, and characterize their origin. Abstract : Linear response time dependent density functional theory (TDDFT), which builds upon configuration interaction singles (CIS) and TD-Hartree–Fock (TDHF), is the most widely used class of excited state quantum chemistry methods and is often employed to study photochemical processes. This paper studies the behavior of the resulting excited state potential energy surfaces beyond the Coulson–Fischer (CF) point in single bond dissociations, when the optimal reference determinant is spin-polarized. Many excited states exhibit sharp kinks at the CF point, and connect to different dissociation limits via a zone of unphysical concave curvature. In particular, the unrestricted M S = 0 lowest triplet T1 state changes character, and does not dissociate into ground state fragments. The unrestricted M S = ±1 T1 CIS states better approximate the physical dissociation limit, but their degeneracy is broken beyond the CF point for most single bond dissociations. On the other hand, the M S = ±1 T1 TDHF states reach the asymptote too soon, by merging with the ground state from the CF point onwards. Use of local exchange–correlation functionals causes M S = ±1 T1 TDDFT states to resemble their unphysical M S = 0 counterpart. The 2 orbital, 2-electron model system of minimal basis H2 is analytically treated to understand the origin of these issues, revealing that the lack of double excitations is at the root of these remarkable observations. The behavior of excited state surfaces is also numerically examined for species like H2, NH3, C2 H6 and LiH in extended basis sets.