Finally, we emphasize that absorption spectra are only one way to characterize electronic transitions. For the performance, e.g. as photosensitizers, the dynamics after photoexcitation is important. Insight from simulations requires the availability of potential energy surfaces and their non-adiabatic and spin-orbit couplings. For example, in the present cases 2, 4, and 5 we have seen that the geometries of the singlet and triplet ground states differ, leading to a change in the ordering of the excited states according to their character. The dependence of the potential energy surface and dynamics for a transition metal complex on the chosen XC functional within TDDFT has recently been shown in Ref. 49. The present work has demonstrated that this particular limitation can be avoided by using BSE@GW, making BSE@GW an ideal candidate for generating potential energy surfaces for dynamics simulations. For complexes such as 2, 4, and 5, however, one has to resort to approximate methods such as the linear vibronic coupling model.