Regardless of whether a polarizable or non-polarizable description is adopted for the solvent molecules, a large number (N, ideally >10,000) of solvent configurations need to be sampled to properly account for the thermal fluctuation in the environment during the excitation and de-excitation processes. This makes QM/MM modeling of solvated chromophore properties roughly N times more expensive than continuum solvent calculations. In general, the higher computational cost of QM/MM models is justified by a more realistic description of the interaction of chromophore with the closest solvent molecules, especially those forming hydrogen bonds with the chromophore.