The simulation box used for MD simulations of solvated systems is sometimes larger than necessary for embedding calculations of molecular properties. In such cases, the system can be truncated to reduce the computational cost. For solute-solvent systems, the so-called droplet model can be a reasonable approach. In this model, all solvent molecules within a given distance from the core region are included while the remainder are removed. Here, we further partitioned the environment into two regions, the inner and outer regions, as depicted in Fig. , where the inner region is described through multipoles and polarizabilities, and the outer region through charges only. The main purpose of the outer region is to avoid artificial polarization at the outer boundary of the inner region, potentially allowing a smaller inner region and thus reducing the computational cost. However, as shown by Kvedaravičiūtė et al., this approach does not always work, and instead, the minimum-image convention could be used. However, in this work, due to a rather inefficient implementation, it is crucial to keep the size of the polarizable part of the environment as small as possible to further reduce the computational costs that arise predominantly during the calculation of the Hessian. We are working on a more efficient implementation, so this is not expected to be an issue in the future.