A general multireference configuration interaction gradient program

An efficient and general method for the computation of analytic energy gradients and energy response properties for general MRCI (multireference configuration interaction) and ACPF (averaged coupled pair functional) wave functions is presented. This methodology includes a general approach, based on successive orbital transformations, for the inclusion of the effects of various orbital resolution (canonicalization) constraints. Initial implementation in the columbus Program System demonstrates, particularly for large�scale multireference wave functions, that the additional computational effort required for the energy gradient is a small fraction of that required for the energy. For polyatomic molecules, the computational resources required for the energy gradient do not depend explicitly on the number of constituent atoms. This combination of features represents a major step forward in the computation and characterization of molecular potential energy surfaces.