and Ω p S = 0 for the ground state case. The whole procedure is repeated until the change of total energy is less than a preselected accuracy, typically 10 -5 Ha. The excitation energy Ω QM/MM S of a complex in the polarizable environment is then obtained as the difference