Experimental and simulated MD-based cumulant spectra for coumarin-153 can be found in Fig. . We note that available experimental spectrum in Ref. 86 uses hexane as the non-polar solvent, rather than toluene as used in our calculations. Due to the fact that both are non-polar solvents with similar dielectric constants, we expect the spectra in the two solvents to match each other closely. Additionally, we note that the weakly interacting nature of the solvent means that this system is likely well-described by more commonly used and computationally affordable Franck-Condon type approaches. While strongly interacting solvents provide a more suitable application for the cumulant method, as direct solute-solvent interactions and slow collective chromophore-environment motions cannot be straightforwardly included in the Franck-Condon method, the well-defined vibronic finestructure of coumarin-153 in weakly interacting solvents allows for a more detailed analysis of the performance of the corrected cumulant approach and is therefore presented here as a test case.