pation of the solvent-delocalized LUMO+1 orbital in the S 1 state of pHBDI -. Additionally, we note that despite the systematic overestimation, TDA/EFP captures the span of the excitation energy better compared to TDA/TIP3P, as indicated by the fitted lines (red) that are more parallel to the diagonals of the parity plots in Fig. .