Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1312.5728 · Physical Review B

Phase diagram and Fermi liquid properties of the extended Hubbard model on the honeycomb lattice

Wéi Wú, A.–M. S. Tremblay

Published May 27, 2014• 48 citations
arXivHTMLPDF

Abstract

The Hubbard and extended Hubbard models on the honeycomb lattice can be seen as prototype models of a single-layer graphene placed in a high dielectric constant environment that screens the Coulomb interaction. Taking advantage of the absence of a sign problem at half-filling, we study this problem with clusters up to 96 sites with the determinant quantum Monte Carlo method as an impurity solver for the dynamical cluster approximation at finite temperatures. After determining the stability of the semimetallic phase to interaction-induced spin-density wave (SDW), charge-density wave (CDW), and Mott insulating phases, we study the single-particle dynamics of the Dirac fermions. We show that when spontaneous symmetry breaking is avoided, the semimetallic phase is a stable Fermi liquid in the presence of repulsive interactions and that Kondo screening dominates the low temperature regime, even though there is a $\ensuremath{\rho}(\ensuremath{\omega})=|\ensuremath{\omega}|$ type local density of states. We also investigate the impact of the correlation effects on the renormalization of the Fermi velocity ${v}_{F}$. We find that ${v}_{F}$ is not renormalized when only on-site repulsion $U$ is present, but that near-neighbor repulsion $V$ does renormalize ${v}_{F}$. This may explain the variations between different measurements of ${v}_{F}$ in graphene.

DOI
10.1103/physrevb.89.205128

Related papers

Extended Hubbard model on a honeycomb lattice

Welberth Kennedy, Sebastião dos Anjos Sousa-Júnior, Natanael C. Costa +1 more

We investigate the phase diagram of the extended Hubbard model on a half-filled honeycomb lattice through determinant quantum Monte Carlo simulations. At zero temperature, the repulsive Hubbard model undergoes a transition to an antiferromagnetic state only for sufficiently strong on-site coupling $U$. However, the inclusion of a repulsive extended interaction $V$ between fermions in nearest-neighbor orbitals induces a charge-ordered phase. For attractive interactions, we also identify superconducting and phase-separated phases. We have determined a complete phase diagram with the aid of structure factors associated with various orderings, combined with measurements of double occupancy and the average sign of fermionic determinants. We have established that the semimetallic phase forms a zone from which ordered phases are excluded, preventing the stabilization of a superconducting phase with $d$-wave symmetry, i.e., only $s$-wave pairing is allowed.

arXiv:2411.07429•Physical review. B./Physical review. B•Apr 11, 2025•4 citations

Competing electronic instabilities of extended Hubbard models on the honeycomb lattice: A functional renormalization group calculation with high-wave-vector resolution

David Sánchez de la Peña, Julian Lichtenstein, Carsten Honerkamp A

We investigate the quantum many-body instabilities for electrons on the honeycomb lattice at half filling with extended interactions, motivated by a description of graphene and related materials. We employ a recently developed fermionic functional renormalization group scheme, which allows for highly resolved calculations of wave-vector dependences in the low-energy effective interactions. We encounter the expected anti-ferromagnetic spin density wave for a dominant on-site repulsion between electrons, and charge order with different modulations for dominant pure ${n}^{\text{th}}$ nearest-neighbor repulsive interactions. Novel instabilities towards incommensurate charge density waves take place when nonlocal density interactions among several bond distances are included simultaneously. Moreover, for more realistic Coulomb potentials in graphene including enough nonlocal terms there is a suppression of charge order due to competition effects between the different charge ordering tendencies, and if the on-site term fails to dominate, the semimetallic state is rendered stable. The possibility of a topological Mott insulator being the favored tendency for dominating second-nearest-neighbor interactions is not realized in our results with high momentum resolution.

arXiv:1606.01124•Physical review. B./Physical review. B•Feb 27, 2017•37 citations

Spectral functions of lattice fermions on the honeycomb lattice with Hubbard and long-range Coulomb interactions

Ho-Kin Tang, Indra Yudhistira, Udvas Chattopadhyay +4 more

The absence of screening of the nonlocal Coulomb interaction in Dirac systems at charge neutrality leads to the breakdown of the Fermi liquid and divergence of the Fermi velocity. On the other hand, Mott-Hubbard physics and the concomitant formation of local moments is dominated by the local effective Hubbard interaction. Using quantum Monte Carlo methods combined with stochastic analytical continuation, we compute the single particle spectral function of fermions on the honeycomb lattice for a realistic interaction that includes both the Hubbard interaction and long-ranged Coulomb repulsion. To a first approximation, we find that the generic high-energy features, such as the formation of the upper Hubbard band near the phase transition, are primarily determined by the local effective Hubbard interaction. In the weakly interacting regime, the long-range Coulomb interaction enhances the bandwidth of quasiparticles and suppresses their lifetime. Conversely, near the phase transition, the long-range Coulomb interaction suppresses the background antiferromagnetic fluctuation, which potentially promotes the propagation of spin polarons, leading to a slight enhancement of the quasiparticle spectral weight and lifetime.

arXiv:2311.12522•Physical review. B./Physical review. B•Oct 9, 2024•4 citations

Fermi-liquid ground state of interacting Dirac fermions in two dimensions

Kazuhiro Seki, Yuichi Otsuka, Seiji Yunoki +1 more

An unbiased zero-temperature auxiliary-field quantum Monte Carlo method is employed to analyze the nature of the semimetallic phase of the two-dimensional Hubbard model on the honeycomb lattice at half filling. It is shown that the quasiparticle weight $Z$ of the massless Dirac fermions at the Fermi level, which characterizes the coherence of zero-energy single-particle excitations, can be evaluated in terms of the long-distance equal-time single-particle Green's function. If this quantity remains finite in the thermodynamic limit, the low-energy single-particle excitations of the correlated semimetallic phase are described by a Fermi-liquid-type single-particle Green's function. Based on the unprecedentedly large-scale numerical simulations on finite-size clusters containing more than 10 000 sites, we show that the quasiparticle weight remains finite in the semimetallic phase below a critical interaction strength. This is also supported by the long-distance algebraic behavior ($\ensuremath{\sim}{r}^{\ensuremath{-}2}$, where $r$ is distance) of the equal-time single-particle Green's function that is expected for the Fermi liquid. Our result thus provides a numerical confirmation of Fermi-liquid theory in two-dimensional correlated metals.

arXiv:1901.06176•Physical review. B./Physical review. B•Mar 25, 2019•12 citations

Spin-Liquid Phase in the Hubbard Model on the Honeycomb Lattice.

Zi Yang Meng, Thomas C. Lang, Stefan Weßel +2 more

The Hubbard model encapsulates the physics of strongly correlated quantum systems in its most basic form. It has been studied intensively in the context of the high-temperature superconductivity. A number of novel phases were recently proposed for Hubbard-like models on the honeycomb lattice, the structure of graphene. We analyzed the Hubbard model of spin- \(\frac{1}{2}\) fermions on the honeycomb lattice at half-filling using large-scale quantum Monte Carlo simulations. We find that the weak coupling semimetal and the antiferromagnetic Mott insulator at strong interaction are separated by an extended gapped phase in an intermediate coupling regime. Exploring excitation gaps, various correlation functions as well as probing for flux quantization, we conclude that a quantum spin liquid, lacking any conventional order, emerges with local charge and spin correlations, best described by a resonating valence bonds state.

IEEE International Conference on High Performance Computing, Data, and Analytics•Jan 1, 2010•0 citations

Cluster dynamical mean field theory of quantum phases on a honeycomb lattice

Rong-Qiang He, Zhong-Yi Lu

We report the cluster dynamical mean field theory calculations performed for the ground state of the half filled Hubbard model on a honeycomb lattice with exact diagonalization on the cluster-impurity solver. Through using elaborate numerical analytic continuation, we identify the existence of a ``spin liquid'' from the on-site interaction $U=0$ to ${U}_{c}$ (between $4.6t$ and $4.85t$) with a smooth crossover correspondingly from the charge fluctuation dominating phase into the charge correlation dominating phase. The semimetallic state exits only at $U=0$. We further find that the magnetic phase transition at ${U}_{c}$ from the spin liquid to the N\'eel antiferromagnetic Mott insulating phase is a first-order quantum phase transition. We also show that the charge fluctuation plays a substantial role on keeping the spin liquid phase against the emergence of a magnetic order.

arXiv:1112.5124•Physical Review B•Jul 6, 2012•27 citations
The Hubbard and extended Hubbard models on the honeycomb lattice can be seen as prototype models of a single-layer graphene placed in a high dielectric constant environment that screens the Coulomb interaction. Taking advantage of the absence of a sign problem at half-filling, we study this problem with clusters up to 96 sites with the determinant quantum Monte Carlo method as an impurity solver for the dynamical cluster approximation at finite temperatures. After determining the stability of the semimetallic phase to interaction-induced spin-density wave (SDW), charge-density wave (CDW), and Mott insulating phases, we study the single-particle dynamics of the Dirac fermions. We show that when spontaneous symmetry breaking is avoided, the semimetallic phase is a stable Fermi liquid in the presence of repulsive interactions and that Kondo screening dominates the low temperature regime, even though there is a \ensuremathρ(\ensuremathω)=∣\ensuremathω∣\ensuremath{\rho}(\ensuremath{\omega})=|\ensuremath{\omega}|\ensuremathρ(\ensuremathω)=∣\ensuremathω∣ type local density of states. We also investigate the impact of the correlation effects on the renormalization of the Fermi velocity vF{v}_{F}vF​. We find that vF{v}_{F}vF​ is not renormalized when only on-site repulsion UUU is present, but that near-neighbor repulsion VVV does renormalize vF{v}_{F}vF​. This may explain the variations between different measurements of vF{v}_{F}vF​ in graphene.
We investigate the phase diagram of the extended Hubbard model on a half-filled honeycomb lattice through determinant quantum Monte Carlo simulations. At zero temperature, the repulsive Hubbard model undergoes a transition to an antiferromagnetic state only for sufficiently strong on-site coupling UUU. However, the inclusion of a repulsive extended interaction VVV between fermions in nearest-neighbor orbitals induces a charge-ordered phase. For attractive interactions, we also identify superconducting and phase-separated phases. We have determined a complete phase diagram with the aid of structure factors associated with various orderings, combined with measurements of double occupancy and the average sign of fermionic determinants. We have established that the semimetallic phase forms a zone from which ordered phases are excluded, preventing the stabilization of a superconducting phase with ddd-wave symmetry, i.e., only sss-wave pairing is allowed.
We investigate the quantum many-body instabilities for electrons on the honeycomb lattice at half filling with extended interactions, motivated by a description of graphene and related materials. We employ a recently developed fermionic functional renormalization group scheme, which allows for highly resolved calculations of wave-vector dependences in the low-energy effective interactions. We encounter the expected anti-ferromagnetic spin density wave for a dominant on-site repulsion between electrons, and charge order with different modulations for dominant pure nth{n}^{\text{th}}nth nearest-neighbor repulsive interactions. Novel instabilities towards incommensurate charge density waves take place when nonlocal density interactions among several bond distances are included simultaneously. Moreover, for more realistic Coulomb potentials in graphene including enough nonlocal terms there is a suppression of charge order due to competition effects between the different charge ordering tendencies, and if the on-site term fails to dominate, the semimetallic state is rendered stable. The possibility of a topological Mott insulator being the favored tendency for dominating second-nearest-neighbor interactions is not realized in our results with high momentum resolution.
The absence of screening of the nonlocal Coulomb interaction in Dirac systems at charge neutrality leads to the breakdown of the Fermi liquid and divergence of the Fermi velocity. On the other hand, Mott-Hubbard physics and the concomitant formation of local moments is dominated by the local effective Hubbard interaction. Using quantum Monte Carlo methods combined with stochastic analytical continuation, we compute the single particle spectral function of fermions on the honeycomb lattice for a realistic interaction that includes both the Hubbard interaction and long-ranged Coulomb repulsion. To a first approximation, we find that the generic high-energy features, such as the formation of the upper Hubbard band near the phase transition, are primarily determined by the local effective Hubbard interaction. In the weakly interacting regime, the long-range Coulomb interaction enhances the bandwidth of quasiparticles and suppresses their lifetime. Conversely, near the phase transition, the long-range Coulomb interaction suppresses the background antiferromagnetic fluctuation, which potentially promotes the propagation of spin polarons, leading to a slight enhancement of the quasiparticle spectral weight and lifetime.
An unbiased zero-temperature auxiliary-field quantum Monte Carlo method is employed to analyze the nature of the semimetallic phase of the two-dimensional Hubbard model on the honeycomb lattice at half filling. It is shown that the quasiparticle weight ZZZ of the massless Dirac fermions at the Fermi level, which characterizes the coherence of zero-energy single-particle excitations, can be evaluated in terms of the long-distance equal-time single-particle Green's function. If this quantity remains finite in the thermodynamic limit, the low-energy single-particle excitations of the correlated semimetallic phase are described by a Fermi-liquid-type single-particle Green's function. Based on the unprecedentedly large-scale numerical simulations on finite-size clusters containing more than 10 000 sites, we show that the quasiparticle weight remains finite in the semimetallic phase below a critical interaction strength. This is also supported by the long-distance algebraic behavior ($\ensuremath{\sim}{r}^{\ensuremath{-}2}$, where rrr is distance) of the equal-time single-particle Green's function that is expected for the Fermi liquid. Our result thus provides a numerical confirmation of Fermi-liquid theory in two-dimensional correlated metals.
The Hubbard model encapsulates the physics of strongly correlated quantum systems in its most basic form. It has been studied intensively in the context of the high-temperature superconductivity. A number of novel phases were recently proposed for Hubbard-like models on the honeycomb lattice, the structure of graphene. We analyzed the Hubbard model of spin- 12\frac{1}{2}21​ fermions on the honeycomb lattice at half-filling using large-scale quantum Monte Carlo simulations. We find that the weak coupling semimetal and the antiferromagnetic Mott insulator at strong interaction are separated by an extended gapped phase in an intermediate coupling regime. Exploring excitation gaps, various correlation functions as well as probing for flux quantization, we conclude that a quantum spin liquid, lacking any conventional order, emerges with local charge and spin correlations, best described by a resonating valence bonds state.
We report the cluster dynamical mean field theory calculations performed for the ground state of the half filled Hubbard model on a honeycomb lattice with exact diagonalization on the cluster-impurity solver. Through using elaborate numerical analytic continuation, we identify the existence of a ``spin liquid'' from the on-site interaction U=0U=0U=0 to Uc{U}_{c}Uc​ (between 4.6t4.6t4.6t and 4.85t4.85t4.85t) with a smooth crossover correspondingly from the charge fluctuation dominating phase into the charge correlation dominating phase. The semimetallic state exits only at U=0U=0U=0. We further find that the magnetic phase transition at Uc{U}_{c}Uc​ from the spin liquid to the N\'eel antiferromagnetic Mott insulating phase is a first-order quantum phase transition. We also show that the charge fluctuation plays a substantial role on keeping the spin liquid phase against the emergence of a magnetic order.