Based on the S 0 optimised geometries, we compute the TDDFT vertical excitation energies of the first 7 excited states using two different solvation methods, LR-PCM and a perturbative statespecific solvation model (pt-SS). At the excited-state geometries we perform TDDFT computations using the LR-PCM and SS-PCM models. Specifically, the following workflow is used: All three PCM approaches start with a RKS/PCM ground-state calculation. LR-PCM, which is the default approach, proceeds with a TDDFT computation including a correction for non-equilibrium solvation added directly to the TDDFT response matrix. The ptSS approach, on the other hand, starts with a zeroth order TDDFT response matrix -not including any corrections for solvation. It proceeds by computing a state-specific correction term based on the zeroth order response vector and the relaxed ground state density. Oscillator strengths (f) were computed as f = 2/3 × ∆E × µ 2 by combining the original transition dipole moments µ with the perturbatively corrected excitation energies ∆E. Finally, SS-PCM is technically carried out via two subsequent TDDFT jobs: in the first step the solvent field is equilibrated to the state of interest and in the second step the DFT orbital optimisation as well as TDDFT are carried out in the equilibrated solvent field.