Through Variational Monte Carlo simulation we show the d-wave RVB pairing in the Heisenberg model on triangular lattice can be better described in terms of a two component order parameter. The fully gapped chiral d-wave RVB state, which is predicted by the mean field theory to be the unique minimum of variational energy in the two dimensional representation space of d-wave pairing, is found to be actually a local maximum and the true minimum of energy is reached by the non-chiral
dxy state with line nodes. We also find that the usual Gutzwiller approximation, which enjoys great success for the square lattice system, fails badly on the triangular lattice as a result of the geometric frustration inherent of the system. An improved version of the Gutzwiiler approximation is proposed and is found to give a much better results than the usual one.
In this work, a long-range resonating valence bond state is proposed as a variational wave function for the ground state of the
S=1/2 antiferromagnetic Heisenberg model on the honeycomb lattice. Employing Variational Monte Carlo (VMC) method, we show that the ground state energy obtained from such RVB wave function, lies well below the energy of the Néel state and compares very well to the energies evaluated from spin-wave theory and series expansion method. We also obtain the spin-spin correlation function along zig-zag and armchair directions and find that the two correlations are different, which indicates the anisotropic nature of the system. We compare our results with the square lattice and we show that although the quantum fluctuations on honeycomb lattice are much stronger, but do not completely destroy the magnetic order.
We theoretically investigate, within Schwinger-boson mean-field theory, the transition from a gapped
Z2 quantum spin liquid, in a ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg spin-
21 system in a honeycomb lattice, to a chiral
Z2 spin-liquid phase under the presence of time-reversal symmetry-breaking scalar chiral interaction (with amplitude ${J}_{\ensuremath{\chi}}$). We numerically obtain a phase diagram of this ${J}_{1}\text{\ensuremath{-}}{J}_{2}\text{\ensuremath{-}}{J}_{\ensuremath{\chi}}$ system, where different ground states are distinguished based on the gap and the nature of excitation spectrum, topological invariant of the excitations, the nature of spin-spin correlation, and the symmetries of the mean-field parameters. The chiral
Z2 state is characterized by the nontrivial Chern number of the excitation bands and lack of long-range magnetic order, which leads to a large thermal Hall coefficient.
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.
We revisit the SU(6) Heisenberg model on the honeycomb lattice, which has been predicted to be a chiral spin liquid by mean-field theory [G. Szirmai et al., Phys. Rev. A 84, 011611(R) (2011)]. Using exact diagonalizations of finite clusters, infinite projected entangled pair state simulations, and variational Monte Carlo simulations based on Gutzwiller projected wave functions, we provide strong evidence that the model with one particle per site and nearest-neighbor exchange actually develops plaquette order. This is further confirmed by the investigation of the model with a ring-exchange term, which shows that there is a transition between the plaquette state and the chiral state at a finite value of the ring-exchange term.
We present a comprehensive computational study of the phase diagram of the frustrated
S=1/2 Heisenberg antiferromagnet on the honeycomb lattice, with second-nearest
(J2) and third-neighbor
(J3) couplings. Using a combination of exact diagonalizations (EDs) of the original spin model, of the Hamiltonian projected into the nearest-neighbor short-range valence-bond basis, and of an effective quantum dimer model, as well as a self-consistent cluster mean-field theory, we determine the boundaries of several magnetically ordered phases in the region ${J}_{2},{J}_{3}\ensuremath{\in}[0,1]$, and find a sizable magnetically disordered region in between. We characterize part of this magnetically disordered phase as a plaquette valence-bond crystal phase. At larger
J2, we locate a sizable region in which staggered valence-bond crystal correlations are found to be important, either due to genuine valence-bond crystal (VBC) ordering or as a consequence of magnetically ordered phases, which break lattice rotational symmetry. Furthermore, we find that a particular parameter-free Gutzwiller projected tight-binding wave function has remarkably accurate energies compared to finite-size extrapolated ED energies along the transition line from conventional N\'eel to plaquette VBC phases, a fact that points to possibly interesting critical behavior---such as a deconfined critical point---across this transition. We also comment on the relevance of this spin model to model the spin liquid region found in the half filled Hubbard model on the honeycomb lattice.