Determination of reduced density matrices in the doubly occupied configuration interaction space: A Hellmann-Feynman theorem approach.

In this work, the Hellmann-Feynman theorem is extended within the doubly occupied configuration interaction space to enable practical calculations of reduced density matrices and expected values. This approach is straightforward, employing finite energy differences, yet remains reliable and accurate even with approximate energies from successive approximation methods. The method's validity is rigorously tested against the Richardson-Gaudin-Kitaev and reduced Bardeen-Cooper-Schrieffer models using approximate excitation energies procured from the Hermitian operator method within the same space, effectively proving the approach's reliability with median error rates for reduced density matrix calculations around 0.

Determination of electronic excitation energies within the doubly occupied configuration interaction space by means of the Hermitian operator method.

In this work, we formulate the equations of motion corresponding to the Hermitian operator method in the framework of the doubly occupied configuration interaction space. The resulting algorithms turn out to be considerably simpler than the equations provided by that method in more conventional spaces, enabling the determination of excitation energies in N-electron systems under an affordable polynomial computational cost. The implementation of this technique only requires to know the elements of low-order reduced density matrices of an N-electron reference state, which can be obtained from any approximate method.