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:1901.06176 · Physical review. B./Physical review. B

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

Kazuhiro Seki, Yuichi Otsuka, Seiji Yunoki, Sandro Sorella

Published Mar 25, 2019• 12 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.99.125145

Related papers

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

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

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.

arXiv:1312.5728•Physical Review B•May 27, 2014•48 citations

Finite-temperature quantum Monte Carlo study of Dirac quantum criticality: Algorithmic advances and comprehensive analysis

Hou-Min Du, Qian Sun, Yuan-Yao He

We present a systematic study of the quantum phase transition from a Dirac semimetal to an antiferromagnetic Mott insulator in the honeycomb-lattice Hubbard model, using finite-temperature auxiliary-field quantum Monte Carlo (AFQMC) simulations. Given the dynamical critical exponent $z=1$ at criticality, we scale the inverse temperature $β$ with the linear system size $L$, i.e., $βt=L$ (with $t$ as the hopping amplitude), to approach the thermodynamic and zero-temperature limits simultaneously. On the algorithmic front, we generalize the fast Fourier transform in path propagation to multi-sublattice systems and employ the delayed update technique in configuration updates to accelerate the AFQMC simulations. These advances enable us to access unprecedented system sizes of up to $3528$ sites (with $L=42$ and $βt = 42$). We additionally use twist-averaged boundary conditions to reduce the finite-size effects in computed observables. On the physical side, we determine the Dirac quantum criticality via finite-size scaling analysis of the spin correlations and off-diagonal single-particle Green's function, and locate the critical point at $U_c/t=3.731(4)$ with the critical exponents $ν=1.097(7)$, $η_ϕ=0.72(2)$ and $η_ψ=0.155(6)$. We further corroborate the critical point by examining several relevant quantities, including the quasiparticle weight, Fermi-liquid parameter, charge compressibility, $U$-derivative of double occupancy, and fidelity susceptibility. Our work establishes an independent benchmark for the chiral Heisenberg Gross-Neveu-Yukawa criticality in the honeycomb Hubbard model, offering a complementary perspective to previous ground-state simulations.

arXiv:2609.14668•Sep 13, 2026•0 citations

Numerical evidence of Luttinger- and Fermi-liquid behavior in the two-dimensional Hubbard model

Sandro Sorella

The two-dimensional Hubbard model with a single spin-up electron interacting with a finite density of spin-down electrons is studied by use of the quantum Monte Carlo technique, the standard second-order perturbation theory and a very good variational wave function introduced recently by Edwards. We also present a new conjugate gradient scheme for the evaluation of this variational state. We performed simulations up to 242 sites at U/t=4 reaching the zero-temperature properties with no ``fermion-sign problem'' and found a surprisingly good accuracy of the Edwards state at low density or low doping. The conjugate gradient method was then applied to the system up to 3698 sites and infinite U. Fermi-liquid theory seems to remain stable in two dimensions for all cases studied with the exception of the half-filling case where a Luttinger-like behavior survives in the Hubbard model, yielding a vanishing quasiparticle weight in the thermodynamic limit.

arXiv:cond-mat/9308001•Physical review. B, Condensed matter•May 1, 1994•7 citations

Quantum criticality in the metal-superconductor transition of interacting Dirac fermions on a triangular lattice

Yuichi Otsuka, Kazuhiro Seki, Sandro Sorella +1 more

We investigate a semimetal-superconductor phase transition of two-dimensional Dirac electrons at zero temperature by large-scale and essentially unbiased quantum Monte Carlo simulations for the half-filled attractive Hubbard model on the triangular lattice, in the presence of alternating magnetic $\ensuremath{\pi}$ flux, that is introduced to construct two Dirac points in the one-particle bands at the Fermi level. This phase transition is expected to describe quantum criticality of the chiral XY class in the framework of the Gross-Neveu model, where, in the ordered phase, the $U(1)$ symmetry is spontaneously broken and a mass gap opens in the excitation spectrum. We compute the order parameter of the $s$-wave superconductivity and estimate the quasiparticle weight from the long-distance behavior of the single-particle Green's function. These calculations allow us to obtain the critical exponents of this transition in a reliable and accurate way. Our estimate for the critical exponents is in good agreement with those obtained for a transition to a Kekul\'e valence bond solid, where an emergent $U(1)$ symmetry is proposed [Z.-X. Li et al., Nat. Commun. 8, 314 (2017)].

arXiv:1803.02001•Physical review. B./Physical review. B•Jul 18, 2018•31 citations

Nonlocal Excitation Spectra in Two-Dimensional Doped Hubbard Model

Yoshiro Kakehashi, Peter Fulde

Single-particle excitation spectra of the two-dimensional Hubbard model on the square lattice near half filling and at zero temperature are investigated on the basis of the self-consistent projection operator method. The method guarantees a high accuracy of the spectra with high energy and high momentum resolutions. It takes into account long-range intersite correlations as well as the strong on-site correlations. Effects of nonlocal excitations are clarified by comparing the results with those of the single-site approximation. The calculated spectra verify the quantum Monte-Carlo results for finite temperatures. The spectra at the Fermi level yield a hole-like Fermi surface in the underdoped region and an electron-like Fermi surface in the overdoped region. From a numerical analysis of the momentum dependent effective mass and self-energy, it is concluded that a marginal Fermi-liquid like state persists even at finite doping concentrations in the strongly correlated region because a van Hove singularity is pinned to the Fermi surface. It is also found that a kink structure appears in the quasiparticle energy band in the same region. The kink is shown to be caused by a mixing between the quasiparticle band and an excitation band with strong short-range antiferromagnetic correlations. These results suggest an explanation for some of the unusual properties of the normal state in high-$T_{\rm c}$ cuprates.

arXiv:0705.0581•Journal of the Physical Society of Japan•Jun 25, 2007•4 citations

Quantum Monte Carlo study of the half-filled Hubbard model on the honeycomb lattice

Yuichi Otsuka, Seiji Yunoki, Sandro Sorella

The half-filled Hubbard model on the honeycomb lattice is investigated by numerically exact large-scale quantum Monte Carlo simulations for lattice sizes up to 2592 sites. By performing careful finite-size scaling for the spin-spin correlation functions, calculated with a high degree of accuracy, we find that the ground state is antiferromagnetically long-range ordered at U / t = 4, where U is the on-site Hubbard interaction and t is the nearest neighbor hopping. Our result is in sharp contrast to a recent report [Meng et al. , Nature 464, 847 (2010)], where instead strong evidence of spin liquid behavior is found.

Journal of Physics Conference Series•Aug 12, 2013•2 citations
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∼r\ensuremath−2\ensuremath{\sim}{r}^{\ensuremath{-}2}\ensuremath∼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 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 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 present a systematic study of the quantum phase transition from a Dirac semimetal to an antiferromagnetic Mott insulator in the honeycomb-lattice Hubbard model, using finite-temperature auxiliary-field quantum Monte Carlo (AFQMC) simulations. Given the dynamical critical exponent z=1z=1z=1 at criticality, we scale the inverse temperature βββ with the linear system size LLL, i.e., βt=Lβt=Lβt=L (with ttt as the hopping amplitude), to approach the thermodynamic and zero-temperature limits simultaneously. On the algorithmic front, we generalize the fast Fourier transform in path propagation to multi-sublattice systems and employ the delayed update technique in configuration updates to accelerate the AFQMC simulations. These advances enable us to access unprecedented system sizes of up to 352835283528 sites (with L=42L=42L=42 and βt=42βt = 42βt=42). We additionally use twist-averaged boundary conditions to reduce the finite-size effects in computed observables. On the physical side, we determine the Dirac quantum criticality via finite-size scaling analysis of the spin correlations and off-diagonal single-particle Green's function, and locate the critical point at Uc/t=3.731(4)U_c/t=3.731(4)Uc​/t=3.731(4) with the critical exponents ν=1.097(7)ν=1.097(7)ν=1.097(7), ηϕ=0.72(2)η_ϕ=0.72(2)ηϕ​=0.72(2) and ηψ=0.155(6)η_ψ=0.155(6)ηψ​=0.155(6). We further corroborate the critical point by examining several relevant quantities, including the quasiparticle weight, Fermi-liquid parameter, charge compressibility, UUU-derivative of double occupancy, and fidelity susceptibility. Our work establishes an independent benchmark for the chiral Heisenberg Gross-Neveu-Yukawa criticality in the honeycomb Hubbard model, offering a complementary perspective to previous ground-state simulations.
The two-dimensional Hubbard model with a single spin-up electron interacting with a finite density of spin-down electrons is studied by use of the quantum Monte Carlo technique, the standard second-order perturbation theory and a very good variational wave function introduced recently by Edwards. We also present a new conjugate gradient scheme for the evaluation of this variational state. We performed simulations up to 242 sites at U/t=4 reaching the zero-temperature properties with no ``fermion-sign problem'' and found a surprisingly good accuracy of the Edwards state at low density or low doping. The conjugate gradient method was then applied to the system up to 3698 sites and infinite U. Fermi-liquid theory seems to remain stable in two dimensions for all cases studied with the exception of the half-filling case where a Luttinger-like behavior survives in the Hubbard model, yielding a vanishing quasiparticle weight in the thermodynamic limit.
We investigate a semimetal-superconductor phase transition of two-dimensional Dirac electrons at zero temperature by large-scale and essentially unbiased quantum Monte Carlo simulations for the half-filled attractive Hubbard model on the triangular lattice, in the presence of alternating magnetic $\ensuremath{\pi}$ flux, that is introduced to construct two Dirac points in the one-particle bands at the Fermi level. This phase transition is expected to describe quantum criticality of the chiral XY class in the framework of the Gross-Neveu model, where, in the ordered phase, the U(1)U(1)U(1) symmetry is spontaneously broken and a mass gap opens in the excitation spectrum. We compute the order parameter of the sss-wave superconductivity and estimate the quasiparticle weight from the long-distance behavior of the single-particle Green's function. These calculations allow us to obtain the critical exponents of this transition in a reliable and accurate way. Our estimate for the critical exponents is in good agreement with those obtained for a transition to a Kekul\'e valence bond solid, where an emergent U(1)U(1)U(1) symmetry is proposed [Z.-X. Li et al., Nat. Commun. 8, 314 (2017)].
Single-particle excitation spectra of the two-dimensional Hubbard model on the square lattice near half filling and at zero temperature are investigated on the basis of the self-consistent projection operator method. The method guarantees a high accuracy of the spectra with high energy and high momentum resolutions. It takes into account long-range intersite correlations as well as the strong on-site correlations. Effects of nonlocal excitations are clarified by comparing the results with those of the single-site approximation. The calculated spectra verify the quantum Monte-Carlo results for finite temperatures. The spectra at the Fermi level yield a hole-like Fermi surface in the underdoped region and an electron-like Fermi surface in the overdoped region. From a numerical analysis of the momentum dependent effective mass and self-energy, it is concluded that a marginal Fermi-liquid like state persists even at finite doping concentrations in the strongly correlated region because a van Hove singularity is pinned to the Fermi surface. It is also found that a kink structure appears in the quasiparticle energy band in the same region. The kink is shown to be caused by a mixing between the quasiparticle band and an excitation band with strong short-range antiferromagnetic correlations. These results suggest an explanation for some of the unusual properties of the normal state in high-TcT_{\rm c}Tc​ cuprates.
The half-filled Hubbard model on the honeycomb lattice is investigated by numerically exact large-scale quantum Monte Carlo simulations for lattice sizes up to 2592 sites. By performing careful finite-size scaling for the spin-spin correlation functions, calculated with a high degree of accuracy, we find that the ground state is antiferromagnetically long-range ordered at U / t = 4, where U is the on-site Hubbard interaction and t is the nearest neighbor hopping. Our result is in sharp contrast to a recent report [Meng et al. , Nature 464, 847 (2010)], where instead strong evidence of spin liquid behavior is found.