| Abstract: | A general formalism is presented for the determination of the optimum natural orbitals within the various strongly orthogonal geminals of a wavefunction describing a 2N-electron closed-shell molecule or atom. The relationship to Hartree–Fock–Roothaan theory is established; the algorithm that is developed is quadratically convergent to the desired result, and does not ignore off diagonal Lagrangian multipliers, or require an infinite series of 2 × 2 orthogonal transformations of the original basis. |
