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

Hubbard model on the honeycomb lattice: From static and dynamical mean-field theories to lattice quantum Monte Carlo simulations

Marcin Raczkowski, Robert Peters, Phùng Thị Thu, Nayuta Takemori, Fakher F. Assaad, A. Honecker, Javad Vahedi

Published Mar 5, 2020• 57 citations
arXivHTMLPDF

Abstract

We study the one-band Hubbard model on the honeycomb lattice using a combination of quantum Monte Carlo (QMC) simulations and static as well as dynamical mean-field theory (DMFT). This model is known to show a quantum phase transition between a Dirac semimetal and the antiferromagnetic insulator. The aim of this paper is to provide a detailed comparison between these approaches by computing static properties, notably ground-state energy, single-particle gap, double occupancy, and staggered magnetization, as well as dynamical quantities such as the single-particle spectral function. At the static mean-field level, local moments cannot be generated without breaking the SU(2) spin symmetry. The DMFT approximation accounts for temporal fluctuations and thus captures both the evolution of the double occupancy and the resulting local moment formation in the paramagnetic phase. As a consequence, the DMFT approximation is found to be very accurate in the Dirac semimetallic phase where local moment formation is present and the spin correlation length small. However, in the vicinity of the fermion quantum critical point, the spin correlation length diverges and the spontaneous SU(2) symmetry breaking leads to low-lying Goldstone modes in the magnetically ordered phase. The impact of these spin fluctuations on the single-particle spectral function---waterfall features and narrow spin-polaron bands---is only visible in the lattice QMC approach.

DOI
10.1103/physrevb.101.125103

Related papers

Phase transitions of the ionic Hubbard model on the honeycomb lattice

Heng-Fu Lin, Haidi Liu, Hong-Shuai Tao +1 more

Many-body problem on the honeycomb lattice systems have been the subject of considerable experimental and theoretical interest. Here we investigate the phase transitions of the ionic Hubbard model on the honeycomb lattice with an alternate ionic potential for the half filling and hole doping cases by means of cellular dynamical mean field theory combining with continue time quantum Monte Carlo as an impurity solver. At half filling, as the increase of the interaction at a fixed ionic potential, we find the single particle gap decreases firstly, reaches a minimum at a critical interaction Uc, then increases upturn. At Uc, there is a band insulator to Mott insulator transition accompanying with the presence of the antiferromagnetic order. Away from half filing, the system shows three phases for the different values of hole density and interaction, paramagnetic metal, antiferromagnetic metal and ferromagnetic metal. Further, we present the staggered particle number, the double occupancy, the staggered magnetization, the uniform magnetization and the single particle spectral properties, which exhibit characteristic features for those phases.

Scientific Reports•May 11, 2015•16 citations

Quantum phase transitions in the honeycomb-lattice Hubbard model

Kazuhiro Seki, Y. Ohta

Quantum phase transitions in the Hubbard model on the honeycomb lattice are investigated in the variational cluster approximation. The critical interaction for the paramagnetic to antiferromagnetic phase transition is found to be in remarkable agreement with a recent large-scale quantum Monte Carlo simulation. Calculated staggered magnetization increases continuously with $U$ and thus we find the phase transition is of a second order. We also find that the semimetal-insulator transition occurs at infinitesimally small interaction and thus a paramagnetic insulating state appears in a wide interaction range. A crossover behavior of electrons from itinerant to localized character found in the calculated single-particle excitation spectra and short-range spin correlation functions indicates that an effective spin model for the paramagnetic insulating phase is far from a simple Heisenberg model with a nearest-neighbor exchange interaction.

arXiv:1209.2101•Sep 6, 2012•0 citations

Quantum Monte Carlo methods and strongly correlated electrons on honeycomb structures

Thomas C. Lang

In this thesis we apply recently developed, as well as sophisticated quantum Monte Carlo methods to numerically investigate models of strongly correlated electron systems on honeycomb structures. The latter are of particular interest owing to their unique properties when simulating electrons on them, like the relativistic dispersion, strong quantum fluctuations and their resistance against instabilities. This work covers several projects including the advancement of the weak-coupling continuous time quantum Monte Carlo and its application to zero temperature and phonons, quantum phase transitions of valence bond solids in spin-1/2 Heisenberg systems using projector quantum Monte Carlo in the valence bond basis, and the magnetic field induced transition to a canted antiferromagnet of the Hubbard model on the honeycomb lattice. The emphasis lies on two projects investigating the phase diagram of the SU(2) and the SU(N)-symmetric Hubbard model on the hexagonal lattice. At sufficiently low temperatures, condensed-matter systems tend to develop order. An exception are quantum spin-liquids, where fluctuations prevent a transition to an ordered state down to the lowest temperatures. Previously elusive in experimentally relevant microscopic two-dimensional models, we show by means of large-scale quantum Monte Carlo simulations of the SU(2) Hubbard model on the honeycomb lattice, 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. Inspired by the rich phase diagrams of SU(N) models we study the SU(N)-symmetric Hubbard Heisenberg quantum antiferromagnet on the honeycomb lattice to investigate the reliability of 1/N corrections to large-N results by means of numerically exact QMC simulations. We study the melting of phases as correlations increase with decreasing N and determine whether the quantum spin liquid found in the SU(2) Hubbard model at intermediate coupling is a specific feature, or also exists in the unconstrained t-J model and higher symmetries.

•Jan 1, 2010•2 citations

Pinning the Order : The Nature of Quantum Criticality in the Hubbard Model on Honeycomb Lattice

Fakher F. Assaad, Igor F. Herbut

In numerical simulations, spontaneously broken symmetry is often detected by computing two-point correlation functions of the appropriate local order parameter. This approach, however, computes the square of the local order parameter, and so when it is small, very large system sizes at high precisions are required to obtain reliable results. Alternatively, one can pin the order by introducing a local symmetrybreaking field and then measure the induced local order parameter infinitely far from the pinning center. The method is tested here at length for the Hubbard model on honeycomb lattice, within the realm of the projective auxiliary-field quantum Monte Carlo algorithm. With our enhanced resolution, we find a direct and continuous quantum phase transition between the semimetallic and the insulating antiferromagnetic states with increase of the interaction. The single-particle gap, measured in units of Hubbard U, tracks the staggered magnetization. An excellent data collapse is obtained by finite-size scaling, with the values of the critical exponents in accord with the Gross-Neveu universality class of the transition.

•Aug 26, 2013•277 citations

The Semimetal-Antiferromagnetic Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice

Johann Ostmeyer, Evan Berkowitz, Stefan Krieg +3 more

The Hubbard model on the honeycomb lattice undergoes a quantum phase transition from a semimetallic to a Mott insulating phase and from a disordered to an anti-ferromagnetically phase. We show that these transitions occur simultaneously and we calculate the critical coupling $U_c=3.835(14)$ as well as the critical exponents $\nu=1.181(43)$ and $\beta= 0.898(37)$ which are expected to fall into the $SU(2)$ Gross-Neveu universality class. For this we employ Hybrid Monte Carlo simulations, extrapolate the single particle gap and the spin structure factors to the thermodynamic and continuous time limits, and perform a data collapse fit. We also determine the zero temperature values of single particle gap and staggered magnetisation on both sides of the phase transition.

arXiv:2110.15441•International Symposium on Lattice Field Theory (LATTICE)•May 16, 2022•1 citation

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
We study the one-band Hubbard model on the honeycomb lattice using a combination of quantum Monte Carlo (QMC) simulations and static as well as dynamical mean-field theory (DMFT). This model is known to show a quantum phase transition between a Dirac semimetal and the antiferromagnetic insulator. The aim of this paper is to provide a detailed comparison between these approaches by computing static properties, notably ground-state energy, single-particle gap, double occupancy, and staggered magnetization, as well as dynamical quantities such as the single-particle spectral function. At the static mean-field level, local moments cannot be generated without breaking the SU(2) spin symmetry. The DMFT approximation accounts for temporal fluctuations and thus captures both the evolution of the double occupancy and the resulting local moment formation in the paramagnetic phase. As a consequence, the DMFT approximation is found to be very accurate in the Dirac semimetallic phase where local moment formation is present and the spin correlation length small. However, in the vicinity of the fermion quantum critical point, the spin correlation length diverges and the spontaneous SU(2) symmetry breaking leads to low-lying Goldstone modes in the magnetically ordered phase. The impact of these spin fluctuations on the single-particle spectral function---waterfall features and narrow spin-polaron bands---is only visible in the lattice QMC approach.
Many-body problem on the honeycomb lattice systems have been the subject of considerable experimental and theoretical interest. Here we investigate the phase transitions of the ionic Hubbard model on the honeycomb lattice with an alternate ionic potential for the half filling and hole doping cases by means of cellular dynamical mean field theory combining with continue time quantum Monte Carlo as an impurity solver. At half filling, as the increase of the interaction at a fixed ionic potential, we find the single particle gap decreases firstly, reaches a minimum at a critical interaction Uc, then increases upturn. At Uc, there is a band insulator to Mott insulator transition accompanying with the presence of the antiferromagnetic order. Away from half filing, the system shows three phases for the different values of hole density and interaction, paramagnetic metal, antiferromagnetic metal and ferromagnetic metal. Further, we present the staggered particle number, the double occupancy, the staggered magnetization, the uniform magnetization and the single particle spectral properties, which exhibit characteristic features for those phases.
Quantum phase transitions in the Hubbard model on the honeycomb lattice are investigated in the variational cluster approximation. The critical interaction for the paramagnetic to antiferromagnetic phase transition is found to be in remarkable agreement with a recent large-scale quantum Monte Carlo simulation. Calculated staggered magnetization increases continuously with UUU and thus we find the phase transition is of a second order. We also find that the semimetal-insulator transition occurs at infinitesimally small interaction and thus a paramagnetic insulating state appears in a wide interaction range. A crossover behavior of electrons from itinerant to localized character found in the calculated single-particle excitation spectra and short-range spin correlation functions indicates that an effective spin model for the paramagnetic insulating phase is far from a simple Heisenberg model with a nearest-neighbor exchange interaction.
In this thesis we apply recently developed, as well as sophisticated quantum Monte Carlo methods to numerically investigate models of strongly correlated electron systems on honeycomb structures. The latter are of particular interest owing to their unique properties when simulating electrons on them, like the relativistic dispersion, strong quantum fluctuations and their resistance against instabilities. This work covers several projects including the advancement of the weak-coupling continuous time quantum Monte Carlo and its application to zero temperature and phonons, quantum phase transitions of valence bond solids in spin-1/2 Heisenberg systems using projector quantum Monte Carlo in the valence bond basis, and the magnetic field induced transition to a canted antiferromagnet of the Hubbard model on the honeycomb lattice. The emphasis lies on two projects investigating the phase diagram of the SU(2) and the SU(N)-symmetric Hubbard model on the hexagonal lattice. At sufficiently low temperatures, condensed-matter systems tend to develop order. An exception are quantum spin-liquids, where fluctuations prevent a transition to an ordered state down to the lowest temperatures. Previously elusive in experimentally relevant microscopic two-dimensional models, we show by means of large-scale quantum Monte Carlo simulations of the SU(2) Hubbard model on the honeycomb lattice, 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. Inspired by the rich phase diagrams of SU(N) models we study the SU(N)-symmetric Hubbard Heisenberg quantum antiferromagnet on the honeycomb lattice to investigate the reliability of 1/N corrections to large-N results by means of numerically exact QMC simulations. We study the melting of phases as correlations increase with decreasing N and determine whether the quantum spin liquid found in the SU(2) Hubbard model at intermediate coupling is a specific feature, or also exists in the unconstrained t-J model and higher symmetries.
In numerical simulations, spontaneously broken symmetry is often detected by computing two-point correlation functions of the appropriate local order parameter. This approach, however, computes the square of the local order parameter, and so when it is small, very large system sizes at high precisions are required to obtain reliable results. Alternatively, one can pin the order by introducing a local symmetrybreaking field and then measure the induced local order parameter infinitely far from the pinning center. The method is tested here at length for the Hubbard model on honeycomb lattice, within the realm of the projective auxiliary-field quantum Monte Carlo algorithm. With our enhanced resolution, we find a direct and continuous quantum phase transition between the semimetallic and the insulating antiferromagnetic states with increase of the interaction. The single-particle gap, measured in units of Hubbard U, tracks the staggered magnetization. An excellent data collapse is obtained by finite-size scaling, with the values of the critical exponents in accord with the Gross-Neveu universality class of the transition.
The Hubbard model on the honeycomb lattice undergoes a quantum phase transition from a semimetallic to a Mott insulating phase and from a disordered to an anti-ferromagnetically phase. We show that these transitions occur simultaneously and we calculate the critical coupling Uc=3.835(14)U_c=3.835(14)Uc​=3.835(14) as well as the critical exponents ν=1.181(43)\nu=1.181(43)ν=1.181(43) and β=0.898(37)\beta= 0.898(37)β=0.898(37) which are expected to fall into the SU(2)SU(2)SU(2) Gross-Neveu universality class. For this we employ Hybrid Monte Carlo simulations, extrapolate the single particle gap and the spin structure factors to the thermodynamic and continuous time limits, and perform a data collapse fit. We also determine the zero temperature values of single particle gap and staggered magnetisation on both sides of the phase transition.
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.