This paper discusses a theorem concerning the variational description of the eigenfunctions and eigenvalues of the two complementary reduced density matrices for a many-particle system in a bound state.
For a system of N identical particles in a pure bound state, the density matrices of positive orders p and $N\ensuremath{-}p$ have the same nonzero eigenvalues with the same multiplicities. If the number of nonzero eigenvalues is finite, these density matrices are unitarily equivalent.
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.''
An explicit reduction in the number of terms of the expansion of the reduced two-particle density matrix for atomic systems is shown to occur when the expansion functions are eigenfunctions of Lz(1)+Lz(2) and Sz(1)+Sz(2).
We consider, for systems of indistinguishable fermions, approximate reconstruction of the three- and four-particle reduced density matrices (RDMs) from the one- and two-particle RDMs, $\mathit{\ensuremath{\gamma}}$ and $\mathit{\ensuremath{\Gamma}}.$ Our ansatz for reconstructing the four-particle RDM is the linear combination $a(\mathit{\ensuremath{\Gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\Gamma}})+b(\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\Gamma}})+c(\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\gamma}}),$ where ``\ensuremath{\wedge}'' denotes the antisymmetrized (Grassmann) product. This is a generalization of reconstruction functionals employed recently to perform direct RDM calculations without wave functions via the contracted Schr\"odinger equation. Here we consider relationships between the parameters a, b, and c that are required in order for the reconstruction functionals to respect the hierarchy of contraction relations between RDMs. To this end we establish several general theorems concerning contractions of antisymmetrized tensor products of $\mathit{\ensuremath{\gamma}},$ $\mathit{\ensuremath{\Gamma}},$ and various products thereof. The accuracy of proposed reconstruction functionals is evaluated using accurate density matrices for the ground state of Be.