The fully-polarizable KS-DFT/AMOEBA method with PWs can easily be extended to compute excited state properties of embedded systems with highly-accurate polarizable potentials, i.e. using a state-specific ∆SCF or a linear-response TD-DFT approach. The implementation of the KS-DFT/AMOEBA model is the first step towards the development of fully-polarizable QM/MM MD simulations for scalable HPC frameworks, such as MiMiC, aiming at feasible unprecedented accuracy multiscale simulations calculations of the complex biomolecular systems.