Searching....

# 7.3.1 Brief Introduction to TDDFT

(December 20, 2021)

Excited states may be obtained from density functional theory by time-dependent density functional theory, 290 Dreuw A., Head-Gordon M.
Chem. Rev.
(2005), 105, pp. 4009.
which calculates poles in the response of the ground state density to a time-varying applied electric field. These poles are Bohr frequencies, or in other words the excitation energies. Operationally, this involves solution of an eigenvalue equation

 $\left(\begin{array}[]{cc}\mathbf{A}&\mathbf{B}\\ \mathbf{B}^{\dagger}&\mathbf{A}^{\dagger}\\ \end{array}\right)\left(\begin{array}[]{c}\mathbf{x}\\ \mathbf{y}\end{array}\right)=\omega\left(\begin{array}[]{cc}\mathbf{-1}&% \mathbf{0}\\ \mathbf{0}&\mathbf{1}\\ \end{array}\right)\left(\begin{array}[]{c}\mathbf{x}\\ \mathbf{y}\end{array}\right)$ (7.15)

where the elements of the matrix $\mathbf{A}$ similar to those used at the CIS level, Eq. (7.11), but with an exchange-correlation correction. 469 Hirata S., Head-Gordon M.
Chem. Phys. Lett.
(1999), 302, pp. 375.
Elements of $\mathbf{B}$ are similar. Equation (7.15) is solved iteratively for the lowest few excitation energies, $\omega$. Alternatively, one can make a CIS-like Tamm-Dancoff approximation (TDA), 470 Hirata S., Head-Gordon M.
Chem. Phys. Lett.
(1999), 314, pp. 291.
in which the “de-excitation” amplitudes $\mathbf{Y}$ are neglected, the $\mathbf{B}$ matrix is not required, and Eq. (7.15) reduces to $\mathbf{Ax}=\omega\mathbf{x}$.

TDDFT is popular because its computational cost is roughly similar to that of the simple CIS method, but a description of differential electron correlation effects is implicit in the method. It is advisable to only employ TDDFT for low-lying valence excited states that are below the first ionization potential of the molecule, or more conservatively, below the first Rydberg state, and in such cases the valence excitation energies are often remarkably improved relative to CIS, with an accuracy of $\sim$0.3 eV for many functionals. 638 Laurent A. D., Jacquemin D.
Int. J. Quantum Chem.
(2013), 113, pp. 2019.
The calculation of the nuclear gradients of full TDDFT and within the TDA is implemented. 698 Liu F. et al.
Mol. Phys.
(2010), 108, pp. 2791.

On the other hand, standard density functionals do not yield a potential with the correct long-range Coulomb tail, owing to the so-called self-interaction problem, and therefore excitation energies corresponding to states that sample this tail (e.g., diffuse Rydberg states and some charge transfer excited states) are not given accurately. 172 Casida M. E. et al.
J. Chem. Phys.
(1998), 108, pp. 4439.
, 1130 Tozer D. J., Handy N. C.
J. Chem. Phys.
(1998), 109, pp. 10180.
, 621 Lange A., Herbert J. M.
J. Chem. Theory Comput.
(2007), 3, pp. 1680.
The extent to which a particular excited state is characterized by charge transfer can be assessed using an a spatial overlap metric proposed by Peach, Benfield, Helgaker, and Tozer (PBHT). 865 Peach M. J. G. et al.
J. Chem. Phys.
(2008), 128, pp. 044118.
(However, see Ref.  966 Richard R. M., Herbert J. M.
J. Chem. Theory Comput.
(2011), 7, pp. 1296.
for a cautionary note regarding this metric.)

Standard TDDFT also does not yield a good description of static correlation effects (see Section 6.12), because it is based on a single reference configuration of Kohn-Sham orbitals. A new variation of TDDFT called spin-flip (SF) DFT was developed by Yihan Shao, Martin Head-Gordon and Anna Krylov to address this issue. 1022 Shao Y., Head-Gordon M., Krylov A. I.
J. Chem. Phys.
(2003), 118, pp. 4807.
SF-DFT is different from standard TDDFT in two ways:

• The reference is a high-spin triplet (quartet) for a system with an even (odd) number of electrons;

• One electron is spin-flipped from an alpha Kohn-Sham orbital to a beta orbital during the excitation.

SF-DFT can describe the ground state as well as a few low-lying excited states, and has been applied to bond-breaking processes, and di- and tri-radicals with degenerate or near-degenerate frontier orbitals. A SF-DFT method with a non-collinear exchange-correlation potential, originally developed by Ziegler and co-workers, 1180 Wang F., Ziegler T.
J. Chem. Phys.
(2004), 121, pp. 12191.
, 1019 Seth M., Mazur G., Ziegler T.
Theor. Chem. Acc.
(2011), 129, pp. 331.
has also been implemented. 93 Bernard Y. A., Shao Y., Krylov A. I.
J. Chem. Phys.
(2012), 136, pp. 204103.
This non-collinear version sometimes improves upon collinear SF-DFT for excitation energies but contains a factor of spin density ($\rho_{\alpha}-\rho_{\beta}$) in the denominator that sometimes causes stability problems. Early experience with SF-DFT suggested that best results are obtained using functionals with $\approx 50$% Hartree-Fock exchange, 1022 Shao Y., Head-Gordon M., Krylov A. I.
J. Chem. Phys.
(2003), 118, pp. 4807.
, 93 Bernard Y. A., Shao Y., Krylov A. I.
J. Chem. Phys.
(2012), 136, pp. 204103.
behavior that was later explained on theoretical grounds by Casida and co-workers. 493 Huix-Rotllant M. et al.
Phys. Chem. Chem. Phys.
(2010), 12, pp. 12811.
Becke’s half-and-half functional BH&HLYP has become something of a standard approach when using SF-DFT. A spin-adapted version of SF-DFT has been developed by Zhang and Herbert. 1273 Zhang X., Herbert J. M.
J. Chem. Phys.
(2015), 143, pp. 234107.
See also Section 7.10.5 for details on wave function-based spin-flip models.