Linear scaling multireference singles and doubles configuration interaction

A linear scaling multireference singles and doubles configuration interaction (MRSDCI) method has been developed. By using localized bases to span the occupied and virtual subspace, local truncation schemes can be applied in tandem with integral screening to reduce the various bottlenecks in a MRSDCI calculation. Among these, the evaluation of electron repulsion integrals and their subsequent transformation, together with the diagonalization of the large CI Hamiltonian matrix, correspond to the most computationally intensive steps in a MRSDCI calculation. We show that linear scaling is possible within each step. The scaling of the method with system size is explored with a system of linear alkane chains and we proceed to demonstrate this method can produce smooth potential energy surfaces via calculating the dissociation of trans-6-dodecene $\left({\mathrm{C}}_{12}{\mathrm{H}}_{24}\right)$ along the central $\mathrm{C}\mathrm{C}$ bond.