At sufficiently low temperatures, condensed-matter systems tend to develop order. A notable exception to this behaviour is the case of quantum spin liquids, in which quantum fluctuations prevent a transition to an ordered state down to the lowest temperatures. There have now been tentative observations of such states in some two-dimensional organic compounds, yet quantum spin liquids remain elusive in microscopic two-dimensional models that are relevant to experiments. Here we show, by means of large-scale quantum Monte Carlo simulations of correlated fermions on a honeycomb lattice (a structure realized in, for example, graphene), that a quantum spin liquid emerges between the state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. This unexpected quantum-disordered state is found to be a short-range resonating valence-bond liquid, akin to the one proposed for high-temperature superconductors: the possibility of unconventional superconductivity through doping therefore arises in our system. We foresee the experimental realization of this model system using ultra-cold atoms, or group IV elements arranged in honeycomb lattices.

In numerical simulations, spontaneously broken symmetry is often detected by computing twopoint correlation functions of the appropriate local order parameter. This approach, however, computes the square of the local order parameter, and so when it is small, very large system sizes at high precisions are required to obtain reliable results. Alternatively, one can pin the order by introducing a local symmetry breaking field, and then measure the induced local order parameter infinitely far from the pinning center. The method is tested here at length for the Hubbard model on honeycomb lattice, within the realm of the projective auxiliary field quantum Monte Carlo algorithm. With our enhanced resolution we find a direct and continuous quantum phase transition between the semi-metallic and the insulating antiferromagnetic states with increase of the interaction. The single particle gap in units of the Hubbard U tracks the staggered magnetization. An excellent data collapse is obtained by finite size scaling, with the values of the critical exponents in accord with the Gross-Neveu universality class of the transition.

No abstract

We study the two-dimensional Kane-Mele-Hubbard model at half filling by means of quantum Monte Carlo simulations. We present a refined phase boundary for the quantum spin liquid. The topological insulator at finite Hubbard interaction strength is adiabatically connected to the groundstate of the Kane-Mele model. In the presence of spin-orbit coupling, magnetic order at large Hubbard U is restricted to the transverse direction. The transition from the topological band insulator to the antiferromagnetic Mott insulator is in the universality class of the three-dimensional XY model. The numerical data suggest that the spin liquid to topological insulator and spin liquid to Mott insulator transitions are both continuous.Comment: 13 pages, 10 figures; final version; new Figs. 4(b) and 8(b

Topological insulators have become one of the most active research areas in condensed matter physics. This article reviews progress on the topic of electronic correlation effects in the two-dimensional case, with a focus on systems with intrinsic spin-orbit coupling and numerical results. Topics addressed include an introduction to the noninteracting case, an overview of theoretical models, correlated topological band insulators, interaction-driven phase transitions, topological Mott insulators and fractional topological states, correlation effects on helical edge states, and topological invariants of interacting systems.

We consider the Kane-Mele model supplemented by a Hubbard U term. The phase diagram is mapped out using projective auxiliary field quantum Monte Carlo simulations. The quantum spin liquid of the Hubbard model is robust against weak spin-orbit interaction, and is not adiabatically connected to the spin-Hall insulating state. Beyond a critical value of U > Uc both states are unstable toward magnetic ordering. In the quantum spin-Hall state we study the spin, charge and single-particle dynamics of the helical Luttinger liquid by retaining the Hubbard interaction only on a ribbon edge. The Hubbard interaction greatly suppresses charge currents along the edge and promotes edge magnetism, but leaves the single-particle signatures of the helical liquid intact.

No abstract

We numerically investigate the critical behavior of the Hubbard model on the honeycomb and the π-flux lattice, which exhibits a direct transition from a Dirac semimetal to an antiferromagnetically ordered Mott insulator. We use projective auxiliary-field quantum Monte Carlo simulations and a careful finite-size scaling analysis that exploits approximately improved renormalization-group-invariant observables. This approach, which is successfully verified for the three-dimensional XY transition of the Kane-Mele-Hubbard model, allows us to extract estimates for the critical couplings and the critical exponents. The results confirm that the critical behavior for the semimetal to Mott insulator transition in the Hubbard model belongs to the Gross-Neveu-Heisenberg universality class on both lattices.