where E HF is the Hartree-Fock energy and F pq and S pq refer to the Fock and overlap matrices, respectively, in the (canonical) MO basis. The Lagrange multipliers Z pq and W pq are determined by using stationary conditions as described later. General MOs are denoted by p and q (without the subscript of "s"). For single-configuration methods, the independent orbital rotations ("indep") consist of (p, q) = (secondary, inactive). Importantly, the MO coefficient matrix C obtained as a result of the self-consistent field calculation satisfies the stationary condition of L HF :