To account for effects of electronic excitations in complex molecular environment, a quantum (QM) region with the excited state complex is embedded in a classical, polarizable atomistic (MM) model for the environment. The QM/MM scheme in VOTCA-XTP makes use of a distributed atomic multipole representation for molecules in the MM region, which allows treatment of both the effects of static electric fields and the polarization response as a selfconsistent reaction field. Specifically, this classical MM energy for the system is evaluated as