Orbital optimization in nonorthogonal multiconfigurational self-consistent field applied to the study of conical intersections and avoided crossings

Nonorthogonal approaches to electronic structure methods have recently received renewed attention, with the hope that new forms of nonorthogonal wavefunction Ansätze may circumvent the computational bottleneck of orthogonal-based methods. The basis in which nonorthogonal configuration interaction is performed defines the compactness of the wavefunction description and hence the efficiency of the method. Within a molecular orbital approach, nonorthogonal configuration interaction is defined by a "different orbitals for different configurations" picture, with different methods being defined by their choice of determinant basis functions. However, identification of a suitable determinant basis is complicated, in practice, by (i) exponential scaling of the determinant space from which a suitable basis must be extracted, (ii) possible linear dependencies in the determinant basis, and (iii) inconsistent behavior in the determinant basis, such as disappearing or coalescing solutions, as a result of external perturbations, such as geometry change. An approach that avoids the aforementioned issues is to allow for basis determinant optimization starting from an arbitrarily constructed initial determinant set. In this work, we derive the equations required for performing such an optimization, extending previous work by accounting for changes in the orthogonality level (defined as the dimension of the orbital overlap kernel between two determinants) as a result of orbital perturbations. The performance of the resulting wavefunction for studying avoided crossings and conical intersections where strong correlation plays an important role is examined.