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:1302.6547 · Physical Review Letters

Comment on “Absence of Spin Liquid in Nonfrustrated Correlated Systems”

A. Liebsch

Published Jul 12, 2013• 6 citations
arXivHTMLPDF

Abstract

In a recent Letter, Hassan and Sénéchal [1] discussed the existence of a spin-liquid phase of the half-filled Hubbard model on the honeycomb lattice. Using schemes, such as the variational cluster approximation (VCA) and the cluster dynamical mean field theory (CDMFT) in combination with exact diagonalization (ED), they argued that a single bath orbital per site of the six-atom unit cell is insufficient and leads to the erroneous conclusion that the system is gapped for all nonzero values of the onsite Coulomb interaction $U$. In contrast, we point out here that, in the case of the honeycomb lattice, six bath levels per six-site unit cell are perfectly adequate for the description of short-range correlations. Instead, we demonstrate that it is the violation of long-range translation symmetry inherent in CDMFT-like schemes which opens a gap at Dirac points. The gap found at small $U$ therefore does not correspond to a Mott gap. As a result, present CDMFT schemes are not suitable for the identification of a spin-liquid phase on the honeycomb lattice. [1] S.R. Hassan and D. Sénéchal, Phys. Rev. Lett. 110, 096402 (2013).

DOI
10.1103/physrevlett.111.029701

Related papers

Absence of Spin Liquid in Nonfrustrated Correlated Systems

Syed R. Hassan, David Sénéchal

The question of the existence of a spin liquid state in the half-filled Hubbard model on the honeycomb (also known as graphene) lattice is revisited. The variational cluster approximation, the cluster dynamical mean field theory, and the cluster dynamical impurity approximation are applied to various cluster systems. Assuming that the spin liquid phase coincides with the Mott insulating phase in this nonfrustrated system, we find that the Mott transition is preempted by a magnetic transition occurring at a lower value of the interaction U, and therefore the spin liquid phase does not occur. This conclusion is obtained using clusters with two bath orbitals connected to each boundary cluster site. We argue that using a single bath orbital per boundary site is insufficient and leads to the erroneous conclusion that the system is gapped for all nonzero values of U.

arXiv:1301.7672•Physical Review Letters•Feb 26, 2013•53 citations

Coulomb correlations in the honeycomb lattice: Role of translation symmetry

A. Liebsch, Wéi Wú

The effect of Coulomb correlations in the half-filled Hubbard model of the honeycomb lattice is studied within the dynamical cluster approximation (DCA) combined with exact diagonalization (ED) and continuous-time quantum Monte Carlo (QMC), for unit cells consisting of six-site rings. The important difference between this approach and the previously employed cluster dynamical mean-field theory (CDMFT) is that DCA preserves the translation symmetry of the system, while CDMFT violates this symmetry. As the Dirac cones of the honeycomb lattice are the consequence of perfect long-range order, DCA yields semimetallic behavior at small on-site Coulomb interactions $U$, whereas CDMFT gives rise to a spurious excitation gap even for very small $U$. This basic difference between the two cluster approaches is found regardless of whether ED or QMC is used as the impurity solver. At larger values of $U$, the lack of translation symmetry becomes less important, so that the CDMFT reveals a Mott gap, in qualitative agreement with large-scale QMC calculations. In contrast, the semimetallic phase obtained in DCA persists even at $U$ values where CDMFT and large-scale QMC consistently show Mott-insulating behavior.

arXiv:1303.2214•Physical Review B•May 20, 2013•24 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

Effective Spin Couplings in the Mott Insulator of the Honeycomb Lattice Hubbard Model

Hong-Yu Yang, A. Fabricio Albuquerque, Andreas Martin Läuchli +1 more

Abstract. Motivated by the recent discovery of a spin liquid phase for the Hubbard model on the honeycomb lattice at half-filling [11], we apply both perturbative and non-perturbative techniques to derive effective spin Hamiltonians describing the low-energy physics of the Mott-insulating phase of the system. Exact diagonalizations of the so-derived models on small clusters are performed, in order to assess the quality of the effective low-energy theory in the spin-liquid regime. We show that six-spin interactions on the elementary loop of the honeycomb lattice are the dominant sub-leading effective couplings. A minimal spin model is shown to reproduce most of the energetic properties of the Hubbard model on the honeycomb lattice in its spin-liquid phase. Surprisingly, a more elaborate effective low-energy spin model obtained by a systematic graph expansion rather disagrees beyond a certain point with the numerical results for the Hubbard model at intermediate couplings.

•Aug 15, 2016•27 citations

Effective spin couplings in the Mott insulator of the honeycomb lattice Hubbard model

Hong-Yu Yang, A. Fabricio Albuquerque, Sylvain Capponi +2 more

Motivated by the recent discovery of a spin-liquid phase for the Hubbard model on the honeycomb lattice at half-filling (Meng et al 2010 Nature 88 487), we apply both perturbative and non-perturbative techniques to derive effective spin Hamiltonians describing the low-energy physics of the Mott-insulating phase of the system. Exact diagonalizations of the so-derived models on small clusters are performed, in order to assess the quality of the effective low-energy theory in the spin-liquid regime. We show that six-spin interactions on the elementary loop of the honeycomb lattice are the dominant sub-leading effective couplings. A minimal spin model is shown to reproduce most of the energetic properties of the Hubbard model on the honeycomb lattice in its spin-liquid phase. Surprisingly, a more elaborate effective low-energy spin model obtained by a systematic graph expansion rather disagrees beyond a certain point with the numerical results for the Hubbard model at intermediate couplings.

arXiv:1207.1072•New Journal of Physics•Nov 28, 2012•27 citations

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
In a recent Letter, Hassan and Sénéchal [1] discussed the existence of a spin-liquid phase of the half-filled Hubbard model on the honeycomb lattice. Using schemes, such as the variational cluster approximation (VCA) and the cluster dynamical mean field theory (CDMFT) in combination with exact diagonalization (ED), they argued that a single bath orbital per site of the six-atom unit cell is insufficient and leads to the erroneous conclusion that the system is gapped for all nonzero values of the onsite Coulomb interaction UUU. In contrast, we point out here that, in the case of the honeycomb lattice, six bath levels per six-site unit cell are perfectly adequate for the description of short-range correlations. Instead, we demonstrate that it is the violation of long-range translation symmetry inherent in CDMFT-like schemes which opens a gap at Dirac points. The gap found at small UUU therefore does not correspond to a Mott gap. As a result, present CDMFT schemes are not suitable for the identification of a spin-liquid phase on the honeycomb lattice. [1] S.R. Hassan and D. Sénéchal, Phys. Rev. Lett. 110, 096402 (2013).
The question of the existence of a spin liquid state in the half-filled Hubbard model on the honeycomb (also known as graphene) lattice is revisited. The variational cluster approximation, the cluster dynamical mean field theory, and the cluster dynamical impurity approximation are applied to various cluster systems. Assuming that the spin liquid phase coincides with the Mott insulating phase in this nonfrustrated system, we find that the Mott transition is preempted by a magnetic transition occurring at a lower value of the interaction U, and therefore the spin liquid phase does not occur. This conclusion is obtained using clusters with two bath orbitals connected to each boundary cluster site. We argue that using a single bath orbital per boundary site is insufficient and leads to the erroneous conclusion that the system is gapped for all nonzero values of U.
The effect of Coulomb correlations in the half-filled Hubbard model of the honeycomb lattice is studied within the dynamical cluster approximation (DCA) combined with exact diagonalization (ED) and continuous-time quantum Monte Carlo (QMC), for unit cells consisting of six-site rings. The important difference between this approach and the previously employed cluster dynamical mean-field theory (CDMFT) is that DCA preserves the translation symmetry of the system, while CDMFT violates this symmetry. As the Dirac cones of the honeycomb lattice are the consequence of perfect long-range order, DCA yields semimetallic behavior at small on-site Coulomb interactions UUU, whereas CDMFT gives rise to a spurious excitation gap even for very small UUU. This basic difference between the two cluster approaches is found regardless of whether ED or QMC is used as the impurity solver. At larger values of UUU, the lack of translation symmetry becomes less important, so that the CDMFT reveals a Mott gap, in qualitative agreement with large-scale QMC calculations. In contrast, the semimetallic phase obtained in DCA persists even at UUU values where CDMFT and large-scale QMC consistently show Mott-insulating behavior.
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.
Abstract. Motivated by the recent discovery of a spin liquid phase for the Hubbard model on the honeycomb lattice at half-filling [11], we apply both perturbative and non-perturbative techniques to derive effective spin Hamiltonians describing the low-energy physics of the Mott-insulating phase of the system. Exact diagonalizations of the so-derived models on small clusters are performed, in order to assess the quality of the effective low-energy theory in the spin-liquid regime. We show that six-spin interactions on the elementary loop of the honeycomb lattice are the dominant sub-leading effective couplings. A minimal spin model is shown to reproduce most of the energetic properties of the Hubbard model on the honeycomb lattice in its spin-liquid phase. Surprisingly, a more elaborate effective low-energy spin model obtained by a systematic graph expansion rather disagrees beyond a certain point with the numerical results for the Hubbard model at intermediate couplings.
Motivated by the recent discovery of a spin-liquid phase for the Hubbard model on the honeycomb lattice at half-filling (Meng et al 2010 Nature 88 487), we apply both perturbative and non-perturbative techniques to derive effective spin Hamiltonians describing the low-energy physics of the Mott-insulating phase of the system. Exact diagonalizations of the so-derived models on small clusters are performed, in order to assess the quality of the effective low-energy theory in the spin-liquid regime. We show that six-spin interactions on the elementary loop of the honeycomb lattice are the dominant sub-leading effective couplings. A minimal spin model is shown to reproduce most of the energetic properties of the Hubbard model on the honeycomb lattice in its spin-liquid phase. Surprisingly, a more elaborate effective low-energy spin model obtained by a systematic graph expansion rather disagrees beyond a certain point with the numerical results for the Hubbard model at intermediate couplings.
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.