Multi-orbital two-particle self-consistent approach -- strengths and limitations

Extending many-body numerical techniques which are powerful in the context of simple model calculations to the realm of realistic material simulations can be a challenging task. Realistic systems often involve multiple active orbitals, which increases the complexity and numerical cost because of the large local Hilbert space and the large number of interaction terms or sign-changing off-diagonal Green's functions. The two-particle self-consistent approach (TPSC) is one such many-body numerical technique, for which multi-orbital extensions have proven to be involved due to the substantially more complex structure of the local interaction tensor. In this paper we extend earlier multi-orbital generalizations of TPSC by setting up two different variants of a fully self-consistent theory for TPSC in multi-orbital systems. We first investigate the strengths and limitations of the approach analytically and then benchmark both variants against dynamical mean-field theory (DMFT) and D-TRILEX results. We find that the exact behavior of the system can be faithfully reproduced in the weak-coupling regime, while at stronger couplings the performance of the two TPSC variants strongly depends on details of the system.