Variational approach for the two-body problem in a multiband extended-Hubbard model

Considering a spin-up and a spin-down fermion in a generic tight-binding lattice with a multi-site basis, here we study the two-body problem using a multiband extended-Hubbard model with arbitrary but finite-ranged hopping and interaction parameters. We derive a linear eigenvalue problem for the entire two-body spectrum, alongside a nonlinear eigenvalue problem for the bound states in the form of a self-consistency equation. Our results are based on an exact variational approach and their versatility offers practical applications in a broad range of lattice geometries. As an illustration, we apply them to the linear-chain model and demonstrate that the resultant spin singlet and triplet bound states are in perfect agreement with the existing literature.