The implementation in the SHARC package was used. Modes with imaginary frequencies or frequencies lower than 100 cm -1 were discarded from the sampling procedure. An ensemble of 700 distorted geometries was generated for further single-point calculations. The ensembles are shown in Figure . Excited states at each geometry in the ensemble were calculated with the time-dependent density functional theory (TD-DFT) utilizing the TURBOMOLE package. The range-separated density functional ωB97X-V was employed together with the def2-TZVP basis set. The calculations employed the resolution of identity approximation, an integration grid of size m3, and the convergence criterion of the self-consistent field method was set to 10 -8 E h . Seven excited states were calculated for each geometry of the ensemble without the Tamm-Dancoff approximation.