in which M q is the effective charge of the M-th MM atom. Note that instead of point charges, Gaussian charges can also be used, which have been shown to be able to avoid overpolarization of the QM electron density. Because the QM/MM electrostatic interaction energy depends on the electron density, the solution to the electron density must come from the variational optimization of