Bose-Hubbard model with occupation-parity couplings

Abstract

We study a Bose-Hubbard model having on-site repulsion, nearest-neighbor tunneling, and ferromagneticlike coupling between occupation parities of nearest-neighbor sites. For a uniform system in any dimension at zero tunneling, we obtain an exact phase diagram characterized by Mott-insulator (MI) and pair liquid phases and regions of phase separation of two MIs. For a general trapped system in one and two dimensions with finite tunneling, we perform quantum Monte Carlo and Gutzwiller mean-field calculations, both of which show the evolution of the system, as the parity coupling increases, from a superfluid to wedding-cake-structure MIs with their occupations jumping by 2. We also identify an exotic pair superfluid at relatively large tunneling strength. Our model ought to effectively describe recent findings in imbalanced Fermi gases in two-dimensional optical lattices and also potentially apply to an anisotropic version of bilinear-biquadratic spin systems.