The new protocol was applied to the photoinduced dynamics of an Fe(II) complex. Here it should be noted that BSE@evGW gave a reasonable absorption spectrum if compared with experiment and TD-DFT using an optimized LC-BLYP functional. More importantly and in contrast to TD-DFT, however, BSE@evGW provided a reasonable assignment of the singlet and triplet transitions. Since the characters of transitions determine vibronic couplings, the BSE@evGW derived LVC parameters can be considered superior if compared to TD-DFT. ML-MCTDH dynamics after pulsed excitation yielded results for the MLCT populations in accord with experiment. The performance of different ML-trees was rather different in terms of convergence with respect to the number of SPFs and associated computational demands.