For all vibronic FC spectra, for both implicit and explicit solvent environments, the time-dependent implementation of one photon absorption and emission procedure was employed with the default adiabatic Hessian model as implemented in the development version of the Gaussian16 electronic structure program. A temperature of 300 K was used for the finite temperature FC spectra. Two advantages of the time-dependent FC implementation are that zero temperature spectra are obtained along with finite temperature spectra at no additional computational cost, and no approximations are made when summing over contributions for combination modes. The average FC lineshapes were generated using ten randomly selected snapshots and aligning their 0-0 transitions for avg 10 FTFC and avg 10 ZTFC methods. We tested the robustness of the lineshape generated with the randomly selected snapshots, see Figure , finding ten snapshots to be adequate.