Beyond the Bethe-Salpeter equation in DFT based computational spectroscopy

We introduce a theoretical framework that accurately describes resonances beyond the reach of both multiplet-based and density-functional-theory (DFT)-based codes. When a resonance acquires strong continuum character, multiplet approaches struggle with the exponential growth of the Hilbert space driven by the increasing number of relevant orbitals, if they extend beyond a few atomic and ligand orbitals. Conversely, existing \emph{ab initio} codes, supplemented by diagrammatic techniques, remain largely confined to the Bethe--Salpeter equation, which tracks only two-particle excitations. However, many systems of interest, particularly those containing open $d$ or $f$ shells, require the propagation of an $N$-particle system, yielding a richer spectral landscape dominated by significant continuum effects. Here, we propose a theoretical framework, together with its numerical implementation, that bridges this gap and enables the rigorous exploration of such resonances.