We generalize the recently introduced dual fermion (DF) formalism for disordered fermion systems by including the effect of interactions. For an interacting disordered system the contributions to the full vertex function have to be separated into ela
stic and inelastic scattering processes, and addressed differently when constructing the DF diagrams. By applying our approach to the Anderson-Falicov-Kimball model and systematically restoring the nonlocal correlations in the DF lattice calculation, we show a significant improvement over the Dynamical Mean-Field Theory and the Coherent Potential Approximation for both one-particle and two-particle quantities.
To reduce the rapidly growing computational cost of the dual fermion lattice calculation with increasing system size, we introduce two embedding schemes. One is the real fermion embedding, and the other is the dual fermion embedding. Our numerical te
sts show that the real fermion and dual fermion embedding approaches converge to essentially the same result. The application on the Anderson disorder and Hubbard models shows that these embedding algorithms converge more quickly with system size as compared to the conventional dual fermion method, for the calculation of both single-particle and two-particle quantities.
