Phases of attractive spin-imbalanced fermions in square lattices

We determine the relative stability of different ground-state phases of spin-imbalanced populations of attractive fermions in square lattices. The ground state is determined within Hartree-Fock-Bogoliubov theory, with care taken to remove finite-size effects. The phases are systematically characterized by the symmetry of the order parameter and their real-and momentum-space structures. For quarter-to half-filled lattices, where the Fermi surfaces are most distorted from their spherical counterpart in the continuum, we predominantly find unidirectional Larkin-Ovchinikov-type phases. We discuss the effect of commensuration between the ordering wave vector and the density imbalance and describe the mechanism of Fermi surface reconstruction and pairing for various orders. A robust supersolid phase exists when the ordering wave vector is diagonally directed. Charge and pairing order coexist, rather than compete, and are responsible for the opening of the gap on different portions of the Fermi surface. A variational determination of the correct pair momentum of the Larkin-Ovchinikov phases shows that phase separation does not occur in the considered regime of density and magnetization.