Direct variational determination of two-particle reduced density matrices for ground and excited states of N-boson systems using the dispersion operator approach

Abstract This work deals with the variational determination of two-particle reduced density matrices corresponding to eigenstates of N -boson systems. Stringent N -representability conditions have been imposed in the variational treatment, in which the resulting optimization problem is addressed within a standard semidefinite programming framework. A unified variational treatment, based on the dispersion operator technique, is proposed. This treatment allows to determine ground- and excited-state energies and their corresponding reduced density matrices, which have been evaluated in systems formulated by means of simple Hamiltonians. We report results from different N -representability levels of treatment. These results are compared with those arising from the full configuration interaction procedure, what allows to point out the relevance of the imposed N -representability conditions. We also highlight the differences found in applying this methodology to fermion and boson systems.