In this work, we use the fragment-based PE model to describe the effects of an environment on a central part of a molecular system. Here the central part is modeled using Kohn-Sham density functional theory (KS-DFT), where the environmental effects are included through an embedding potential. The potential is created by splitting the environment into small fragments, and then from individual KS-DFT calculations on each fragment, atomcentered multipoles and polarizabilities are obtained. The multipoles and polarizabilities are used to model interactions between the quantum subsystem and the permanent and induced charge distributions, respectively, of the fragments in the environment. The induced charge distributions are modeled by induced dipoles which are the result of the electric fields at the polarizable sites, i.e., atoms with a polarizability. The MIC can be used to calculate the shortest distance between particles defined within a simulation box with PBC. Here we use it to calculate interactions between the polarizable sites and single images of the other sites, i.e., using the MIC to find the minimum-image distance. This corresponds to centering a given polarizable site in the simulation box and folding back sites that fall outside the box as illustrated in Figure . Thus, only the first replicas of the simulation box are included. In the following, we will only describe the induced electrostatics because the permanent electrostatics remain the same as in the PE model without MIC. For a more detailed description of the PE model in general, we refer to refs 13-15. In practice, the PE model is combined with a QM method, in this case, KS-DFT. Eventually, an effective KS operator is obtained that consists of the usual KS operator ( f QM ) and a PE operator (v PE )