We followed a protocol like the one described in our earlier work in order to extract the Q y excitation energies of the Chl molecules . The Hamiltonian R2018 and the "standard" parameter set were employed with a 1.0 E h threshold for the initial and the production runs with the requirement for 12 roots during the DFT/MRCI calculations. Furthermore, eight HOMO and six LUMO states were considered as active space including single and double excitations. In addition, all calculations were repeated in a QM/MM setup using the crystal structure coordinates for each pigment.