Examining lineshapes computed with the E-sumFTFC approach, Figure , we note that both absorption and fluorescence spectra are generated using the same one hundred snapshots obtained from the ground and excited state AIMD as in the ensemble spectral simulations but with the ground and excited state geometry optimizations and frequency calculations performed with frozen QM solvent. The individual FTFC spectra (thin gray lines) show a significant degree of variability in the lineshape and position, confirming the importance of explicit environmental interactions on the resulting vibronic spectrum. The FWHM of individual explicit solvent FTFC spectra, see Figure , range between 0.260 and 0.780 eV, with an average FWHM value of 0.366 eV for absorption and 0.514 eV for emission. These average values are larger than the FWHM values for the implicit solvent FTFC spectra of 0.354 and 0.458 eV, suggesting that for NBD, the QM explicit solvent environment leads to stronger vibronic coupling than the implicit solvent environment, even more so for fluorescence.