The seam space of a conical intersection is really a (hyper)surface of dimension , and while the two electronic states in question are degenerate at every point within this space, the electronic energy varies from one point to the next. To provide a simple picture of photochemical reaction pathways, it is often convenient to locate the minimum-energy crossing point (MECP) within this -dimensional seam. Two separate minimum-energy pathway searches, one on the excited state starting from the ground-state geometry and terminating at the MECP, and the other on the ground state starting from the MECP and terminating at the ground-state geometry, then affords a photochemical mechanism. (See Ref. 1236 for a simple example.) In some sense, then, the MECP is to photochemistry what the transition state is to reactions that occur on a single Born-Oppenheimer potential energy surface. One should be wary of pushing this analogy too far, because whereas a transition state reasonably be considered to be a bottleneck point on the reaction pathway, the path through a conical intersection may be downhill and perhaps therefore more likely to proceed from one surface to the other at a point “near" the intersection, and in addition there can be multiple conical intersections between the same pair of states so more than one photochemical mechanism may be at play. Such complexity could be explored, albeit at significantly increased cost, using non-adiabatic “surface hopping" ab initio molecular dynamics, as described in Section 9.9.6. Here we describe the computationally-simpler procedure of locating an MECP along a conical seam.
Recall that the branching space around a conical intersection between electronic states and is spanned by two vectors, [Eq. (9.6)] and [Eq. (9.7)]. While the former is readily available in analytic form for any electronic structure method that has analytic excited-state gradients, the non-adiabatic coupling vector is not available for most methods. For this reason, several algorithms have been developed to optimize MECPs without the need to evaluate , and three such algorithms are available in Q-Chem.
Martínez and coworkers
J. Phys. Chem. B
(2008), 112, pp. 405. developed a penalty-constrained MECP optimization algorithm that consists of minimizing the objective function
where is a fixed parameter to avoid singularities and is a Lagrange multiplier for a penalty function meant to drive the energy gap to zero. Minimization of is performed iteratively for increasingly large values .
A second MECP optimization algorithm is a simplification of the penalty-constrained approach that we call the “direct” method. Here, the gradient of the objective function is
is the mean energy gradient, with being the nuclear gradient for state , and
is the normalized difference gradient. Finally,
projects the gradient difference direction out of the mean energy gradient in Eq. (9.10). The algorithm then consists in minimizing along the gradient , with for the iterative cycle over a Lagrange multiplier, which can sometimes be slow to converge.
The third and final MECP optimization algorithm that is available in Q-Chem
is the branching-plane updating method developed by Morokuma and
J. Chem. Theory Comput.
(2010), 6, pp. 1538. and implemented in Q-Chem by Zhang and Herbert. 1236 J. Phys. Chem. B
(2014), 118, pp. 7806. This algorithm uses a gradient that is similar to that in Eq. (9.10) but projects out not just in Eq. (9.13) but also a second vector that is orthogonal to it, representing an iteratively-updated approximation to the branching space.
None of these three methods requires evaluation of non-adiabatic couplings, and
all three can be used to optimize MECPs at the CIS, SF-CIS, TDDFT, SF-TDDFT,
and SOS-CIS(D0) levels. The direct algorithm can also be used for EOM-XX-CCSD
methods (XX = EE, IP, or EA). It should be noted that since EOM-XX-CCSD is a
linear response method, it suffers from the same topology problem around
conical intersections involving the ground state that was described in regards
to TDDFT in Section 9.8.1. With spin-flip approaches,
correct topology is obtained.
J. Chem. Phys.
(2014), 141, pp. 064104.
Analytic derivative couplings are available for (SF-)CIS and (SF-)TDDFT, so for
these methods one can alternatively employ an optimization algorithm that makes
use of both and . Such an algorithm,
due to Schlegel and coworkers,
Chem. Phys. Lett.
(1994), 223, pp. 269. is available in Q-Chem and consists of optimization along the gradient in Eq. (9.10) but with a projector
in place of the projector in Eq. (9.13).
Equation (9.14) has the effect of projecting the span of
and (i.e., the branching space) out
of state-averaged gradient in Eq. (9.10). The tends to reduce
the number of iterations necessary to converge the MECP, and since calculation
of the (optional) vector represents only a slight amount
of overhead on top of the (required) vector, this last
algorithm tends to yield significant speed-ups relative to the other
J. Chem. Phys.
(2014), 141, pp. 064104. As such, it is the recommended choice for (SF-)CIS and (SF-)TDDFT.
It should be noted that while the spin-flip methods cure the topology problem
around conical intersections that involve the ground state, this method tends
to exacerbate spin contamination relative to the corresponding spin-conserving
J. Chem. Phys.
(2015), 143, pp. 234107. While spin contamination is certainly present in traditional, spin-conserving CIS and TDDFT, it presents the following unique challenge in spin-flip methods. Suppose, for definiteness, that one is interested in singlet excited states. Then the reference state for the spin-flip methods should be the high-spin triplet. A spin-flipping excitation will then generate S, S, S but will also generate the component of the triplet reference state, which therefore appears in what is ostensibly the singlet manifold. Q-Chem attempts to identify this automatically, based on a threshold for , but severe spin contamination can sometimes defeat this algorithm, 1236 J. Phys. Chem. B
(2014), 118, pp. 7806. hampering Q-Chem’s ability to distinguish singlets from triplets (in this particular example). An alternative might be the state-tracking procedure that is described in Section 9.8.5.