Fakher F. Assaad

Quantum spin liquid emerging in two-dimensional correlated Dirac fermions

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.

Correlation effects in two-dimensional topological insulators

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.

Dynamical signatures of edge-state magnetism on graphene nanoribbons.

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.

Correlation Effects in Quantum Spin-Hall Insulators: A Quantum Monte Carlo Study

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>U(c) 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.

Pinning the order: the nature of quantum criticality in the Hubbard model on honeycomb lattice

A new algorithm of quantum Monte Carlo simulations designed to detect very weak magnetic order allows high-resolution studies of the correlation between magnetic order and electrical insulation in Mott insulators. Specific investigations of the Hubbard model on honeycomb lattice---a paradigmatic model for Mott insulators---show that the correlation emerges through a continuous quantum phase transition belonging to a new universality class.

Spectral properties of the one-dimensional Hubbard model.

The spectral properties of the 1-D Hubbard model are obtained from quantum Monte Carlo simulations using the maximum entropy method. The one-particle excitations are characterized by dispersive cosine-like bands. Velocities for spin- and charge excitations are obtained that lead to a conformal charge c=0.98 +/- 0.05 for the largest system simulated (N=84). An exact sum-rule for the spin-excitations is fulfilled accurately with deviations of at most 10% only around 2 kF.

QUANTUM MONTE CARLO SIMULATIONS OF THE HALF-FILLED TWO-DIMENSIONAL KONDO LATTICE MODEL

The 2D half-filled Kondo lattice model with exchange J and nearest neighbor hopping t is considered. It is shown that this model belongs to a class of Hamiltonians for which zero-temperature auxiliary field Monte Carlo methods may be efficiently applied. We compute the staggered moment, spin and quasiparticle gaps on lattice sizes up to 12 X 12. The competition between the RKKY interaction and Kondo effect leads to a continuous quantum phase transition between antiferromagnetic and spin-gaped insulators. This transition occurs at J/t = 1.45 \pm 0.05.

Fermionic quantum criticality in honeycomb and π-flux Hubbard models: Finite-size scaling of renormalization-group-invariant observables from quantum Monte Carlo

We numerically investigate the critical behavior of the Hubbard model on the honeycomb and the

Single-hole dynamics in thet−Jmodel on a square lattice

\pi

Spin and charge dynamics of the ferromagnetic and antiferromagnetic two-dimensional half-filled Kondo lattice model

-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.