On the direct determination of constrained pure state one-electron density matrices: a new algorithm

An algorithm based on a "constrained variational principle" is suggested for the direct determination of constrained one-electron density matrices in a self-consistent manner. The algorithm uses both penalty-function and Lagrangian multiplier methods of incorporating equality constraints and can tackle any number of constraints. Results of preliminary calculations are presented