In order to obtain the polarization response of both the QM and MM regions, a selfconsistent procedure is employed. At step p of this procedure, the total energy of the coupled QM/MM system for the state S of interest (ground state S = 0, or excited states S > 0) is determined as