For geometry optimization in explicit QM solvent and VEE/VDE calculations, the 2020 development version of TeraChem 110 was used. Frequency calculations in frozen QM solvent and vibronic spectra calculations were carried out using a developmental version of Gaussian16. The integral equation formalism for PCM as implemented in Gaussian16 was used for modeling solvent effects on implicit solvent FC absorption and emission spectra.