Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules

Abstract

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems. © 2014 AIP Publishing LLC. [http://dx

Next‐generation quantum theory of atoms in molecules for the ground and excited state of the ring‐opening of cyclohexadiene

“…Regarding the choice of the density functional, see Refs. , where it is recommended to employ functional with sufficiently large fraction of the exact exchange to eliminate the undesired effect of the self‐interaction error of the conventional GGA functionals. The geometries of the S 1 and S 0 state minima, S 1 /S 0 conical intersections and the minimum energy paths (MEPs) were optimized using the SSR analytical energy derivatives formalism described in Ref.…”

Spin‐restricted ensemble‐referenced Kohn–Sham method: basic principles and application to strongly correlated ground and excited states of molecules

“…This is illustrated in Figure where, for a few organic molecules, the root mean square deviations (RMSDs) of the SI‐SA‐REKS geometries at the minimum energy CI points from the MR ab initio geometries are shown . The SI‐SA‐REKS method is also capable of accurately computing the non‐adiabatic coupling parameters between the ground and excited states in the form of the branching plane vectors of conical intersections . The branching plane vectors are obtained by differentiation of the respective parts of the SI‐SA‐REKS energy expression, see Refs , and for the detailed expressions.…”

“…To round up this section, ensemble DFT for excited states as implemented in the SI‐SA‐REKS method is a versatile and accurate approach to the calculation of various types of excitations in molecular systems. A wide range of excited states, which are otherwise not accessible with the use of TD‐DFT, can be studied, including the CT excitations, excitations in extended π ‐conjugated systems, such as cyanines and polyacenes, excitations in molecules undergoing bond breaking/bond formation, conical intersections between the ground and excited electronic states, etc. It is also noteworthy that the SI‐SA‐REKS results can be obtained at an essentially mean‐field cost, avoiding a steeper scaling of the linear response formalism of TD‐DFT.…”

Next‐generation quantum theory of atoms in molecules for the S1/S0 conical intersections in dynamics trajectories of a light‐driven rotary molecular motor

“…The SSR method enables one to describe electron correlation effects arising in due to multireference character of the ground and excited electronic states through the use of ensemble density functional theory . The SSR method has been rigorously tested in the past and proved to be capable of describing the PESs of the ground and excited electronic states with an accuracy matching the most advanced multireference wavefunction methods, including the crossings between the PESs, the so‐called CIs …”