Analytic formulation of derivative coupling vectors for complete active space configuration interaction wavefunctions with floating occupation molecular orbitals.

The floating occupation molecular orbital complete active space configuration interaction (FOMO-CASCI) method is quite promising for the study of nonadiabatic processes. Use of this method directly in nonadiabatic dynamics simulations has been limited by the lack of available first-order nonadiabatic coupling vectors. Here, an analytic formulation of these derivative coupling vectors is presented for FOMO-CASCI wavefunctions using a simple Lagrangian-based approach. The derivative coupling vectors are applied in the optimization of minimum energy conical intersections of an aqueously solvated model compound for the chromophore of the green fluorescent protein (including 100 water molecules). The computational cost of the FOMO-CASCI derivative coupling vector is shown to scale quadratically, O(N2), with system size and is applied to systems with up to 1000 atoms.