We present the multichannel Dyson equation that combines two or more many-body Green’s functions to describe the electronic structure of materials. In this work we use it to model photoemission spectra by coupling the one-body Green’s function with the three-body Green’s function. We demonstrate that, unlike methods using only the one-body Green’s function, our approach puts the description of quasiparticles and satellites on an equal footing. We propose a multichannel self-energy that is static and only contains the bare Coulomb interaction, making frequency convolutions and self-consistency unnecessary. Despite its simplicity, we demonstrate with a diagrammatic analysis that the physics it describes is extremely rich. Finally, we present a framework based on an effective Hamiltonian that can be solved for any many-body system using standard numerical tools. We illustrate our approach by applying it to the Hubbard dimer and show that it is exact both at $1/4$ and $1/2$ filling.

• Received 13 June 2023
• Revised 2 October 2023
• Accepted 20 October 2023

DOI:https://doi.org/10.1103/PhysRevLett.131.216401