For the two solvated chromophores, one can perform QM/MM calculations with TDDFT/TDA applied to the QM region. Here we investigate two MM descriptions of water: (i) the non-polarizable TIP3P model, and (ii) the effective fragment potential (EFP) for water, which is a polarizable MM model. In the case of EFP, the excitation energies are collected without the polarization correction (denoted as TDA/EFP0) and with the correction (denoted as TDA/EFP). As shown in Fig. , TDA/MM models consistently overestimate the lowest vertical excitation energy with an RMSD of 0.1 eV or greater for both biochromophores in 100 representative solvent configurations. In terms of r 2 value, a better correlation between TDA/MM and full-TDA excitation energies is obtained for OLU -with both MM models, which might be related to the partici- and) pHBDI -anions in 30 water molecules compared against TDA excitation energies for the full chromophore-solvent complexes. TDA/EFP0 and TDDFT/EFP results were obtained without and with the polarization correction. For both biochromophores, 100 solvent configurations were collected from a 100 ps NVT trajectory. TDA calculations were performed using the ωB97X-D functional and 6-311++G(d,p) basis set. The red asterisks indicate the TDDFT/TDA excitation energy of each chromophore in the gas phase.