We note that the splitting in the calculated spectrum in toluene is considerably overestimated compared to the experimental one (see Fig. ). There are several possible reasons that could explain this discrepancy: 1) an oversampling of contact ion configurations that cause large symmetry-breaking and S 1 -S 2 splitting; 2) errors in the CV-Cl -dynamics and computation of vertical excitation energies due to the treatment of the counterion as a classical point charge; 3) errors in the computed S 1 -S 2 splittings for the strongly symmetry-broken configuration due to the choice of DFT functional.