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.

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

No abstract

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.

We extend the weak-coupling diagrammatic determinantal algorithm to projective schemes as well as to the inclusion of phonon degrees of freedom. The projective approach provides a very efficient algorithm to access zero temperature properties. To implement phonons, we integrate them out in favor of a retarded density-density interaction and simulate the resulting purely electronic action with the weak-coupling diagrammatic determinantal algorithm. Both extensions are tested within the dynamical mean field approximation for the Hubbard and Hubbard-Holstein models.

We investigate the edge-state magnetism of graphene nanoribbons using projective quantum Monte Carlo simulations and a self-consistent mean-field approximation of the Hubbard model. The static magnetic correlations are found to be short ranged. Nevertheless, the correlation length increases with the width of the ribbon such that already for ribbons of moderate widths we observe a strong trend towards mean-field-type ferromagnetic correlations at a zigzag edge. These correlations are accompanied by a dominant low-energy peak in the local spectral function and we propose that this can be used to detect edge-state magnetism by scanning tunneling microscopy. The dynamic spin structure factor at the edge of a ribbon exhibits an approximately linearly dispersing collective magnonlike mode at low energies that decays into Stoner modes beyond the energy scale where it merges into the particle-hole continuum.

We employ quantum Monte Carlo techniques to calculate the Z2 topological invariant in a twodimensional model of interacting electrons that exhibits a quantum spin Hall topological insulator phase. In particular, we consider the parity invariant for inversion-symmetric systems, which can be obtained from the bulk's imaginary-time Green's function after an appropriate continuation to zero frequency. This topological invariant is used here in order to study the trivial-band to topological-insulator transitions in an interacting system with spin-orbit coupling and an explicit bond dimerization. We discuss the accessibility and behavior of this topological invariant within quantum Monte Carlo simulations.