SU(2) gauge theory of the Hubbard model and application to the honeycomb lattice

Abstract

Motivated by recent experiments on the triangular lattice Mott-Hubbard system $\ensuremath{\kappa}\text{\ensuremath{-}}{(\mathrm{BEDT}\text{\ensuremath{-}}\mathrm{TTF})}_{2}{\mathrm{Cu}}_{2}{(\mathrm{C}\mathrm{N})}_{3}$, we develop a general formalism to investigate quantum spin liquid insulators adjacent to the Mott transition in Hubbard models. This formalism, dubbed the SU(2) slave-rotor formulation, is an extension of the SU(2) gauge theory of the Heisenberg model to the case of the Hubbard model. Furthermore, we propose the honeycomb lattice Hubbard model (at half filling) as a candidate for a spin liquid ground state near the Mott transition; this is an appealing possibility, as this model can be studied via quantum Monte Carlo simulation without a sign problem. The pseudospin symmetry of Hubbard models on bipartite lattices turns out to play a crucial role in our analysis, and we develop our formalism primarily for the case of a bipartite lattice. We also sketch its development for a general Hubbard model. We develop a mean-field theory to describe spin liquids and some competing states, and apply it to the honeycomb lattice. On the insulating side of the Mott transition, we find an SU(2) algebraic spin liquid (ASL), described by gapless $S=1∕2$ Dirac fermions (spinons) coupled to a fluctuating SU(2) gauge field. This result contrasts with that obtained via a U(1) slave-rotor approach, which instead found a U(1) ASL. That formulation does not respect the pseudospin symmetry, and is therefore not correct on the honeycomb lattice. We construct a low-energy effective theory describing the ASL phase, the conducting semimetal phase and the Mott transition between them. This physics can be detected in numerical simulations via the simultaneous presence of substantial antiferromagnetic and valence-bond solid correlations. In the SU(2) ASL, these observables have slowly decaying fluctuations in space and time, described by power laws with the same critical exponent. We recover the results of the U(1) slave-rotor formulation in the presence of a strong breaking of pseudospin symmetry. Our analysis suggests that both a third-neighbor electron hopping, and/or pseudospin-breaking terms such as a nearest-neighbor density interaction, may help to stabilize a spin liquid phase.