Improving the computational efficiency of the implementation is a major goal of future work on this method. The current implementation lacks an iterative diagonalization of the LASSI Hamiltonian, which is a necessary feature if systems with more than about 3 fragments are to be considered. Alternatively, we could consider layering on an approximation to the diagonalization of the LASSI Hamiltonian along the lines of, for example, block-correlated coupled cluster. We also anticipate integration with L-PDFT, 36 a more systematic exploration of the effects of orbital optimization, a possible generalization to a fully variational algorithm, and integration with the effects of spin-orbit coupling in future work.