The ground-state fermion second-order reduced density matrix (2-RDM) is determined variationally using itself as a basic variable. As necessary conditions of the N-representability, we used the positive semidefiniteness conditions, P, Q, and G conditions that are described in terms of the 2-RDM. The variational calculations are performed by using recently developed semidefinite programming algorithm (SDPA). The calculated energies of various closed- and open-shell atoms and molecules are excellent, overshooting only slightly the full-CI energies. There was no case where convergence was not achieved. The calculated properties also reproduce well the full-CI results.
With the use of the density equation [H. Nakatsuji, Phys. Rev. A 14, 41 (1976)], the second-order density matrices are directly determined without any use of the wave functions. The third- and fourth-order reduced density matrices (RDM's) are decoupled into lower-order ones using the Green's function technique. This method is applied to Be, Ne,
H2O,
H3O+, N
H3, C
H4, ${\mathrm{BH}}_{4}^{\phantom{\rule{0ex}{0ex}}\ensuremath{-}}$,
NH4+, and C
H3F, and the results are successfully compared with the full configuration interaction results. The convergence is fairly good, and the calculated second-order RDM's almost satisfy the necessary conditions of the N representability, the P, Q, and G conditions, and the first-order RDM's are exactly N representable. These results show that the present method is very promising.