Motivated by the recent numerical evidence [Z. Meng, T. Lang, S. Wessel, F. Assaad, and A. Muramatsu, Nature (London) 464, 847 (2010)] of a short-range resonating valence bond state in the honeycomb lattice Hubbard model, we consider Schwinger boson mean field theories of possible spin liquid states on honeycomb lattice. From general stability considerations the possible spin liquids will have gapped spinons coupled to
Z2 gauge field. We apply the projective symmetry group method to classify possible
Z2 spin liquid states within this formalism on honeycomb lattice. It is found that there are only two relevant
Z2 states, differed by the value of gauge flux, zero or $\ensuremath{\pi}$, in the elementary hexagon. The zero-flux state is a promising candidate for the observed spin liquid and continuous phase transition into commensurate N\'eel order. We also derive the critical field theory for this transition, which is the well-studied O(4) invariant theory [A. V. Chubukov, T. Senthil, and S. Sachdev, Phys. Rev. Lett. 72, 2089 (1994); A. V. Chubukov, S. Sachdev, and T. Senthil, Nucl. Phys. B 426, 601 (1994); S. V. Isakov, T. Senthil, and Y. B. Kim, Phys. Rev. B 72, 174417 (2005)], and has an irrelevant coupling between Higgs and boson fields with cubic power of spatial derivatives as required by lattice symmetry. This is in sharp contrast to the conventional theory [S. Sachdev and N. Read, Int. J. Mod. Phys. B 5, 219 (1991)], where such transition generically leads to incommensurate magnetic order. In this scenario the
Z2 spin liquid could be close to a tricritical point. Soft boson modes will exist at seven different wave vectors. This will show up as low-frequency dynamical spin susceptibility peaks not only at the $\ensuremath{\Gamma}$ point (the N\'eel order wave vector) but also at Brillouin-zone-edge center
M points and twelve other points. Some simple properties of the $\ensuremath{\pi}$-flux state are studied as well. Symmetry allowed further neighbor mean field ansatz are derived in appendices which can be used in future theoretical works along this direction.
Spin liquids are novel states of matter with fractionalized excitations. A recent numerical study of Hubbard model on a honeycomb lattice\cite{Meng2010} indicates that a gapped spin liquid phase exists close to the Mott transition. Using Projective Symmetry Group, we classify all the possible spin liquid states by Schwinger fermion mean-field approach. We find there is only one fully gapped spin liquid candidate state: "Sublattice Pairing State" that can be realized up to the 3rd neighbor mean-field amplitudes, and is in the neighborhood of the Mott transition. We propose this state as the spin liquid phase discovered in the numerical work. To understand whether SPS can be realized in the Hubbard model, we study the mean-field phase diagram in the
J1−J2 spin-1/2 model and find an s-wave pairing state. We argue that s-wave pairing state is not a stable phase and the true ground state may be SPS. A scenario of a continuous phase transition from SPS to the semimetal phase is proposed. This work also provides guideline for future variational studies of Gutzwiller projected wavefunctions.
We study incompressible ground states of bosons in a two-dimensional rotating square optical lattice. The system can be described by the Bose-Hubbard model in an effective uniform magnetic field present due to the lattice rotation. To study ground states of the system, we map it to a frustrated spin model, followed by Schwinger boson mean-field theory and projective symmetry group analysis. Using symmetry analysis, we identify bosonic fractional quantum Hall states, predicted for bosonic atoms in rotating optical lattices, with possible stable gapped spin-liquid states within the Schwinger boson formalism. In particular, our results indicate that previously obtained fractional quantum Hall states induced by the lattice potential, and with no counterpart in the continuum [M\"oller and Cooper, Phys. Rev. Lett. 103, 105303 (2009)], correspond to ``$\ensuremath{\pi}$-flux'' spin-liquid states of the frustrated spin model.