Variationally optimized density matrices in relation to chemical bonding

Variational optimization of the second-order density matrix is an alternative to wave function based ab initio methods. The main difficulty in this method is the 'N-representability problem', which is usually cannot be solved exactly. Does this method yield correct chemical properties? This issue is investigated using our own optimization program, applied to relevant chemical problems.