An additional source of inhomogeneous broadening is the full polarization of the chromophore-solvent environment. This polarization, along with potential charge transfer, can be achieved by using a QM description of the first solvent shell in the VEE calculations. For the same QM/MM trajectory snapshots, where there is no structural changes of the chromophore or solvent, we computed the vertical excitation energies for both the chromophore and all methanol molecules within 5 Å of the chromophore at the QM level, with the remaining methanol molecules at the MM level (≈ 230 QM atoms). This QM solvent shell treatment in the excitation energies leads to a shift to lower energies, see SI section S8, and some minor additional inhomogeneous broadening, see Figure . The corresponding spectral density, see SI section S9, shows that this QM treatment of the solvent in the VEE calculations leads to minimal change in the vibronic coupling of the high-frequency chromophore degrees of freedom. Despite the increased broadening with QM treatment of solvent for the VEEs, which is aligned with previous studies, the simulated spectrum based on the QM/MM trajectory produces an absorption spectrum more narrow than experiment.