The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry.

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g.

Analytic non-adiabatic derivative coupling terms for spin-orbit MRCI wavefunctions. II. Derivative coupling terms and coupling angle for KHeAΠ⇔KHeBΣ.

A method for calculating the analytic nonadiabatic derivative coupling terms (DCTs) for spin-orbit multi-reference configuration interaction wavefunctions is reviewed. The results of a sample calculation using a Stuttgart basis for KHe are presented. Additionally, the DCTs are compared with a simple calculation based on the Nikitin's 3 × 3 description of the coupling between the Σ and Π surfaces, as well as a method based on Werner's analysis of configuration interaction coefficients.

Analytic non-adiabatic derivative coupling terms for spin-orbit MRCI wavefunctions. I. Formalism.

Analytic gradients of electronic eigenvalues require one calculation per nuclear geometry, compared to at least 3n + 1 calculations for finite difference methods, where n is the number of nuclei. Analytic nonadiabatic derivative coupling terms (DCTs), which are calculated in a similar fashion, are used to remove nondiagonal contributions to the kinetic energy operator, leading to more accurate nuclear dynamics calculations than those that employ the Born-Oppenheimer approximation, i.e.