In this work a large new class of extreme points of the set of N-representable 2-matrices is obtained. They correspond to the exact ground state of an operator B(g) which generalizes the BCS Hamiltonian. Further, three hitherto unpublished necessary conditions for N-representability are obtained.
Usually, variational calculations for the second-order reduced density matrix are performed subject to the constraint that the ``P'' and ``Q'' matrices are positive semidefinite, which only constrains the lowest eigenvalue of these matrices. We characterize the highest eigenvalue of these matrices and discuss how the associated constraint (which is related to the ground-state energy of the Hamiltonians $H=\ensuremath{-}P$ and $H=\ensuremath{-}Q$) can be implemented in practical calculations. This necessary condition for N-representability should help ensure that the second-order reduced density matrix is not ``overcorrelated.''
In this work, we analyze the effectiveness of different sets of well-known necessary N-representability conditions that are imposed on the variational determination of two-electron reduced density matrices corresponding to the ground and excited states of N-electron Hamiltonians within doubly occupied configuration interaction subspaces. The two-electron reduced density matrix elements are determined so that they are suitably related to those of higher-order reduced density matrices and other associated matrices, up to a given order, constituting positive semidefinite matrices. These conditions up to two-, three-, and four-order have been sequentially imposed on the variational calculation of the expectation values of the dispersion operator corresponding to several Hamiltonian models, obtaining the electronic spectra. These results have been contrasted with those arising from the full configuration interaction method, illuminating the sufficiency/insufficiency of the variational N-representability conditions applied.
We derive necessary conditions for any two-fermion reduced density matrix (2-RDM) to be representable by a pure N-fermion density matrix $\mathrm{\ensuremath{\Psi}}{\mathrm{\ensuremath{\Psi}}}^{*}$ where $\mathrm{\ensuremath{\Psi}}$ is the wave function. These pure N-representability conditions of the 2-RDM are important because they provide stringent constraints beyond those from the Pauli and the generalized Pauli constraints on the structure of many-fermion 2-RDMs and their wave functions. The pure 2-RDM conditions are derived as generalized Pauli constraints on effective one-fermion reduced density matrices (1-RDMs) generated by the removal or addition of a fermion from the wave function. Computationally, we show for four-electron molecules that the derived pure N-representability conditions are nontrivially active for exact ground-state 2-RDMs and that they provide significant restrictions beyond the D, Q, and G ensemble N-representability conditions. Constraints on higher-order p-RDMs where p>2 are derived in a similar fashion. The constraints have potentially significant applications to computing strongly correlated many-fermion states with enhanced accuracy and decreased computational complexity.
Given the set 𝒮N of all operator endomorphisms T with the property that whenever D is any 2-density operator, then T (D) is N-representable is considered. We specify a pair {Q2,PN} of extreme points of 𝒮N which are related to each other via the particle-hole duality (Theorem 1.2). The rest of the paper is devoted to a complete characterization of the set ŜN of all operator endomorphisms T with the property that whenever D is any 2-density operator then T (D) is quasi-N-representable, i.e., satisfies the necessary conditions Id, QN, BN, and CN.
The reduction of density matrices defined with respect to a finite basis set is considered. A basis set can be introduced into the vector space of Hermitian matrices acting on functions of the coordinates of p particles for each value of p. An analogy with the construction of spin eigenfunctions is used to obtain a particular basis set called the reducing basis. Reduction of one of these basis elements is either one to one or maps into the origin. Any element in the preimage of a density matrix can thus be resolved into two components, one uniquely determined by the density matrix and the other arbitrary within a certin subset. A restatement of the N-representability problem is given, and two sufficient conditions and one necessary condition for N representability are given in terms of distances and norms.