X

Search Results

Searching....

9.9 Nonadiabatic Couplings and Optimization of Minimum-Energy Crossing Points

9.9.1 Nonadiabatic Couplings

(April 13, 2024)

Conical intersections are degeneracies between Born-Oppenheimer potential energy surfaces that facilitate nonadiabatic transitions between excited states, i.e., internal conversion and intersystem crossing processes, both of which represent a breakdown of the Born-Oppenheimer approximation. 823 Matsika S., Krause P.
Annu. Rev. Phys. Chem.
(2011), 62, pp. 621.
Link
, 502 Herbert J. M. et al.
Acc. Chem. Res.
(2016), 49, pp. 931.
Link
Although simultaneous intersections between more than two electronic states are possible, 823 Matsika S., Krause P.
Annu. Rev. Phys. Chem.
(2011), 62, pp. 621.
Link
consider for convenience the two-state case, and let

𝐇=(HJJ(𝐑)HJK(𝐑)HJK(𝐑)HKK(𝐑)). (9.46)

denote the matrix representation of the vibronic (vibrational + electronic) Hamiltonian [Eq. (4.2)] in a basis of two electronic states, J and K. (Electronic degrees of freedom have been integrated out of this expression, and 𝐑 represents the remaining, nuclear coordinates.) By definition, the Born-Oppenheimer states are the ones that diagonalize 𝐇 at a particular molecular geometry 𝐑, and thus two conditions must be satisfied in order to obtain degeneracy in the Born-Oppenheimer representation: HJJ=HKK and HJK=0. As such, degeneracies between two Born-Oppenheimer potential energy surfaces exist in subspaces of dimension Nint-2, where Nint=3Natoms-6 is the number of internal (vibrational) degrees of freedom (assuming the molecule is non-linear). This (Nint-2)-dimensional subspace is known as the seam space because the two states are degenerate everywhere within this space. In the remaining two degrees of freedom, known as the branching space, the degeneracy between Born-Oppenheimer surfaces is lifted by an infinitesimal displacement, which in a three-dimensional plot resembles a double cone about the point of intersection, hence the name conical intersection.

The branching space is defined by the span of a pair of vectors 𝐠JK and 𝐡JK. The former is simply the difference in the gradient vectors of the two states in question,

𝐠JK=EJ𝐑-EK𝐑, (9.47)

and is readily evaluated at any level of theory for which analytic energy gradients are available (or less-readily, via finite difference, if they are not!). The definition of the nonadiabatic coupling vector 𝐡JK, on the other hand, is more involved and not directly amenable to finite-difference calculations:

𝐡JK=ΨJ|(H^/𝐑)|ΨK. (9.48)

This is closely related to the derivative coupling vector

𝐝JK=ΨJ|(/𝐑)|ΨK=𝐡JKEJ-EK. (9.49)

The latter expression for 𝐝JK demonstrates that the coupling between states becomes large in regions of the potential surface where the two states are nearly degenerate. The relative orientation and magnitudes of the vectors 𝐠JK and 𝐡JK determined the topography around the intersection, i.e., whether the intersection is “peaked” or “sloped”; 50 Atchity G. J., Xantheas S. S., Ruedenberg K.
J. Chem. Phys.
(1991), 95, pp. 1862.
Link
, 1363 Yarkony D. R.
J. Chem. Phys.
(2001), 114, pp. 2601.
Link
see Ref.  502 Herbert J. M. et al.
Acc. Chem. Res.
(2016), 49, pp. 931.
Link
for a pedagogical overview.

Algorithms to compute the nonadiabatic couplings 𝐝JK are not widely available in quantum chemistry codes, but thanks to the efforts of the Herbert and Subotnik groups, they are available in Q-Chem when the wave functions ΨJ and ΨK, and corresponding electronic energies EJ and EK, are computed at the CIS or TDDFT level, 338 Fatehi S. et al.
J. Chem. Phys.
(2011), 135, pp. 234105.
Link
, 1378 Zhang X., Herbert J. M.
J. Chem. Phys.
(2014), 141, pp. 064104.
Link
, 1379 Zhang X., Herbert J. M.
J. Chem. Phys.
(2015), 142, pp. 064109.
Link
, 921 Ou Q. et al.
J. Chem. Phys.
(2015), 142, pp. 064114.
Link
or at the corresponding spin-flip (SF) levels of theory (SF-CIS or SF-TDDFT). The spin-flip implementation 1378 Zhang X., Herbert J. M.
J. Chem. Phys.
(2014), 141, pp. 064104.
Link
is particularly significant, because only that approach—and not traditional spin-conserving CIS or TDDFT—affords correct topology around conical intersections that involve the ground state.

To understand why, suppose that J in Eq. (9.46) represents the ground state; call it J=0 for definiteness. In linear response theory (TDDFT) or in CIS (by virtue of Brillouin’s theorem), the coupling matrix elements between the reference (ground) state and all of the excited states vanish identically, hence H0K(𝐑)0. This means that there is only one condition to satisfy in order to obtain degeneracy, hence the branching space is one- rather than two-dimensional, for any conical intersection that involves the ground state. 730 Levine B. G. et al.
Mol. Phys.
(2006), 104, pp. 1039.
Link
(For intersections between two excited states, the topology should be correct.) In the spin-flip approach, however, the reference state has a different spin multiplicity than the target states; if the latter have spin quantum number S, then the reference state has spin S+1. This has the effect that the ground state of interest (spin S) is treated as an excitation, and thus on a more equal footing with excited states of the same spin, and it rigorously fixes the topology problem around conical intersections. 1378 Zhang X., Herbert J. M.
J. Chem. Phys.
(2014), 141, pp. 064104.
Link
, 502 Herbert J. M. et al.
Acc. Chem. Res.
(2016), 49, pp. 931.
Link
This can have consequences for simulation of internal conversion to the ground state, where for example S1 lifetimes may be significantly affected by warping of the potential energy surfaces around the S1/S0 intersection in conventional TDDFT. 1381 Zhang X., Herbert J. M.
J. Chem. Phys.
(2021), 155, pp. 124111.
Link

Nonadiabatic (derivative) couplings are available for both CIS and TDDFT. The CIS nonadiabatic couplings can be obtained from direct differentiation of the wave functions with respect to nuclear positions. 338 Fatehi S. et al.
J. Chem. Phys.
(2011), 135, pp. 234105.
Link
, 1378 Zhang X., Herbert J. M.
J. Chem. Phys.
(2014), 141, pp. 064104.
Link
For TDDFT, the same procedure can be carried out to calculate the approximate nonadiabatic couplings, in what has been termed the “pseudo-wave function” approach. 922 Ou Q. et al.
J. Chem. Phys.
(2014), 141, pp. 024114.
Link
, 1378 Zhang X., Herbert J. M.
J. Chem. Phys.
(2014), 141, pp. 064104.
Link
Formally more rigorous TDDFT nonadiabatic couplings derived from quadratic response theory are also available, although they are subject to certain undesirable, accidental singularities if for the two states J and K in Eq. (9.48), the energy difference |EJ-EK| is quasi-degenerate with the excitation energy ωI=EI-E0 for some third state, I. 1379 Zhang X., Herbert J. M.
J. Chem. Phys.
(2015), 142, pp. 064109.
Link
, 921 Ou Q. et al.
J. Chem. Phys.
(2015), 142, pp. 064114.
Link
As such, the pseudo-wave function method is the recommended approach for computing nonadiabatic couplings with TDDFT, although in the spin flip case the pseudo-wave function approach is rigorously equivalent to the pseudo-wave function approach, and is free of singularities. 1379 Zhang X., Herbert J. M.
J. Chem. Phys.
(2015), 142, pp. 064109.
Link

Finally, we note that there is mounting evidence that SF-TDDFT calculations are most accurate when used with functionals containing 50% Hartree-Fock exchange, 1117 Shao Y., Head-Gordon M., Krylov A. I.
J. Chem. Phys.
(2003), 118, pp. 4807.
Link
, 539 Huix-Rotllant M. et al.
Phys. Chem. Chem. Phys.
(2010), 12, pp. 12811.
Link
and many studies with this method have used the BH&HLYP functional; see Refs.  502 Herbert J. M. et al.
Acc. Chem. Res.
(2016), 49, pp. 931.
Link
and for reviews. The BH&HLYP functional combines LYP correlation is combined with Becke’s “half and half” (BH&H) exchange functional, consisting of 50% Hartree-Fock exchange and 50% Becke88 exchange (EXCHANGE = BHHLYP in Q-Chem.)