Firstorder derivative couplings between excited states from adiabatic TDDFT response theory
Abstract
We present a complete derivation of derivative couplings between excited states in the framework of adiabatic timedependent density functional response theory. Explicit working equations are given and the resulting derivative couplings are compared with derivative couplings from a pseudowavefunction ansatz. For degenerate excited states, i.e., close to a conical intersection (CI), the two approaches are identical apart from an antisymmetric overlap term. However, if the difference between two excitation energies equals another excitation energy, the couplings from response theory exhibit an unphysical divergence. This spurious behavior is a result of the adiabatic or static kernel approximation of timedependent density functional theory leading to an incorrect analytical structure of the quadratic response function. Numerical examples for couplings close to a CI and for wellseparated electronic states are given.
 Authors:
 Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104 (United States)
 Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 926972025 (United States)
 Publication Date:
 OSTI Identifier:
 22416099
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Chemical Physics; Journal Volume: 142; Journal Issue: 6; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; ADIABATIC APPROXIMATION; COMPARATIVE EVALUATIONS; DENSITY FUNCTIONAL METHOD; EQUATIONS; EXCITATION; EXCITED STATES; KERNELS; RESPONSE FUNCTIONS; TIME DEPENDENCE; WAVE FUNCTIONS
Citation Formats
Ou, Qi, Subotnik, Joseph E., Email: subotnik@sas.upenn.edu, Bellchambers, Gregory D., and Furche, Filipp, Email: filipp.furche@uci.edu. Firstorder derivative couplings between excited states from adiabatic TDDFT response theory. United States: N. p., 2015.
Web. doi:10.1063/1.4906941.
Ou, Qi, Subotnik, Joseph E., Email: subotnik@sas.upenn.edu, Bellchambers, Gregory D., & Furche, Filipp, Email: filipp.furche@uci.edu. Firstorder derivative couplings between excited states from adiabatic TDDFT response theory. United States. doi:10.1063/1.4906941.
Ou, Qi, Subotnik, Joseph E., Email: subotnik@sas.upenn.edu, Bellchambers, Gregory D., and Furche, Filipp, Email: filipp.furche@uci.edu. 2015.
"Firstorder derivative couplings between excited states from adiabatic TDDFT response theory". United States.
doi:10.1063/1.4906941.
@article{osti_22416099,
title = {Firstorder derivative couplings between excited states from adiabatic TDDFT response theory},
author = {Ou, Qi and Subotnik, Joseph E., Email: subotnik@sas.upenn.edu and Bellchambers, Gregory D. and Furche, Filipp, Email: filipp.furche@uci.edu},
abstractNote = {We present a complete derivation of derivative couplings between excited states in the framework of adiabatic timedependent density functional response theory. Explicit working equations are given and the resulting derivative couplings are compared with derivative couplings from a pseudowavefunction ansatz. For degenerate excited states, i.e., close to a conical intersection (CI), the two approaches are identical apart from an antisymmetric overlap term. However, if the difference between two excitation energies equals another excitation energy, the couplings from response theory exhibit an unphysical divergence. This spurious behavior is a result of the adiabatic or static kernel approximation of timedependent density functional theory leading to an incorrect analytical structure of the quadratic response function. Numerical examples for couplings close to a CI and for wellseparated electronic states are given.},
doi = {10.1063/1.4906941},
journal = {Journal of Chemical Physics},
number = 6,
volume = 142,
place = {United States},
year = 2015,
month = 2
}

Derivative couplings between TDDFT excited states obtained by direct differentiation in the TammDancoff approximation
Working within the TammDancoff approximation, we calculate the derivative couplings between timedependent densityfunctional theory excited states by assuming that the KohnSham superposition of singly excited determinants represents a true electronic wavefunction. All Pulay terms are included in our derivative coupling expression. The reasonability of our approach can be established by noting that, for closely separated electronic states in the infinite basis limit, our final expression agrees exactly with the ChernyakMukamel expression (with transition densities from response theory). Finally, we also validate our approach empirically by analyzing the behavior of the derivative couplings around the T{sub 1}/T{sub 2} conical intersection ofmore » 
Analytic derivative couplings in timedependent density functional theory: Quadratic response theory versus pseudowavefunction approach
We revisit the formalism for analytic derivative couplings between excited states in timedependent density functional theory (TDDFT). We derive and implement these couplings using quadratic response theory, then numerically compare this responsetheory formulation to couplings implemented previously based on a pseudowavefunction formalism and direct differentiation of the KohnSham determinant. Numerical results, including comparison to full configuration interaction calculations, suggest that the two approaches perform equally well for many molecular systems, provided that the underlying DFT method affords accurate potential energy surfaces. The response contributions are found to be important for certain systems with high symmetry, but can be calculated withmore » 
Firstorder nonadiabatic coupling matrix elements between excited states: A Lagrangian formulation at the CIS, RPA, TDHF, and TDDFT levels
Analytic expressions for the firstorder nonadiabatic coupling matrix elements between electronically excited states are first formulated exactly via both timeindependent equation of motion and timedependent response theory, and are then approximated at the configuration interaction singles, particlehole/particleparticle random phase approximation, and timedependent density functional theory/HartreeFock levels of theory. Note that, to get the Pulay terms arising from the derivatives of basis functions, the standard response theory designed for electronic perturbations has to be extended to nuclear derivatives. The results are further recast into a Lagrangian form that is similar to that for excitedstate energy gradients and allows to use atomicmore »