 The equations of motion are obtained by finding the stationary point of the Lagrangian under variations of "ψi" and R"I", with the orthogonality constraint.where Λ"ij" is a Lagrangian multiplier matrix to comply with the orthonormality constraint.In the formal limit where "μ" → 0, the equations of motion approach Born–Oppenheimer molecular dynamics.