This eigenvalue equation can be solved iteratively together with the expectation constraint as well as the normalization constraint, leading to solutions of f 0 , |A 0 ⟩, and Ẽ0 . When the constraints are satisfied, the total energy, which is a function of the nuclear expectation positions, will serve as the effective PES, or the CMES: