The multiscale strategy followed in this work is summarized in Figure and is further explained in each of the subsections of this section. It comprises five stages: after a system preparation, we perform MM-MD equilibration and production runs, followed by QM/MM-MD configurational sampling, and TD-DFT for the calculation of the MPA spectra combined with fragment-based polarizable embedding for describing the MM environment.