In practical terms, an important aspect of this so-called polarizable density embedding (PDE) is that the computationally rather demanding freeze-and-thaw iterations for mutual polarization of the subsystems is replaced by a polarizable MM treatment in terms of induced dipoles. A similar strategy was later on also used by Jensen and co-workers by combining a polarizable MM strategy with FDE and the level-shift operator technique to ensure external orthogonality.