For a proper treatment of L-edge spectra of transition metal compounds, multiconfigurational methods are mandatory. Therefore, the hemin monomer excitations were described using the RASSCF approach, but for the sake of demonstration, the coupling between different excitations on the separated monomers was included via the dipole approximation only. But also for valence excitations, there are countless examples of systems that require the use of a multiconfigurational wavefunction method. In recent work, Xie and Troisi 19 have applied multiconfigurational methods to a supermolecule followed by a diabatization procedure to parameterize the Frenkel model. As has been mentioned, such a scheme starting from the supermolecule is connected to significant computational costs. We are not aware that a multiconfigurational wavefunction approach to the direct parametrization of the FEH avoiding supermolecular calculation or dipole approximation exists beyond the example quoted above.