Coexisting Néel and charge density wave orders in attractive three-color fermions

In optical lattices attractive ultracold fermions with three hyperfine-spin components (colors) can form three fermionic configurations depending on the interactions---unbound fermions, on-site trions, and off-site trions---leading to the coexistence of multiple fermionic species in the ordered phase, which demonstrates that the attractive three-color fermions are unique from other correlated fermion systems and may host intriguing phases and phase transitions. At temperatures below the superexchange energy scale, we employ the determinant quantum Monte Carlo (QMC) method to investigate the phases and phase transitions in the half-filled attractive three-color Hubbard model on a honeycomb lattice where the Hubbard interactions are color dependent (anisotropic interactions) and the coupling between color 3 and colors (1, 2) serves as a control parameter. In the coupling regime where on-site and off-site trions coexist, our QMC simulations demonstrate coexisting N\'eel and charge density wave (CDW) orders which are common in condensed matter but rare in ultracold atoms. At weak coupling where the color superfluid (CSF) order is scattered by color-3 fermions, we find that very small coupling of color 3 with colors (1, 2) can destroy the CSF order and that the vanishing of the CSF order is not immediately accompanied by the emergence of the on-site trionic phase, which strikingly disagrees with the prevalent results of dynamical mean-field theory. The underlying mechanisms of the coexisting N\'eel/CDW orders and of the CSF order breaking are also presented based on intuitive physical pictures.