Combining an explicit description of the solvent environment with vibronic coupling is essential for an improved description of both absorption and fluorescence spectra. Next, we compare spectral lineshapes generated with the E-FC methods. The pictorial representation of the working of E-FC methods for NR in DMSO is shown in Figure . The individual FTFC spectra and resulting E-sum 100 FTFC lineshape are shown in Fig. . The individual FTFC spectra computed in the explicit solvent environment are quite similar to each other across configurations and they do not show the same intensity in the second vibronic peak as is seen with the implicit solvent FTFC spectra. The FWHM values for the individual spectra do not usually account for the contributions from the second vibronic peak, leading to very small values. Therefore, in Figure , we instead plot the fullwidth half-half max (FWHHM), taking the width at 25% of the lineshape maximum. At 25% of the maximum, we observe broadening primarily from two vibronic peaks in both absorption and emission lines. However, there are some single-peak outliers in both FWHHM distributions. For absorption, widths below 0.20 eV originate from a single peak, which in Figure we represent by light blue bars. For emission, widths below 0.20 eV are due to a combination of both single and double peaks; in Figure , these values are colored light blue with gray stripes. We computed the average value of this FWHHM using only the spectra for which there are two vibronic peaks contributing to the FWHHM value, which we determined to be 0.274 for absorption and 0.264 eV for fluorescence. These values are smaller than the corresponding implicit solvent FWHHM values of 0.452 and 0.470 eV taken at the 25% intensity value, showing decreased vibronic coupling in explicit solvent. The summation of the E-sum 100 FTFC method washes out the pronounced vibronic nature of the peaks, providing relatively smooth absorption and fluorescence lineshapes, as observed in the experimental lineshapes.