In any case, apart from the condition that the APCs minimize E 0 , we also desire that they mimic the perturbation that one fragment causes to other fragments in the full system. Hence, we force the APC of each atom to be equal to the APC-embedded ground-state ESP charge of the atom. Under this restriction, the variational optimization of E 0 can be solved by the following iterative procedure: