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

Electronic ground state in bilayer graphene with realistic Coulomb interactions

Jia Ning Leaw, Ho Kin Tang, Pinaki Sengupta, Fakher F. Assaad, Igor F. Herbut, Shaffique Adam

Published Sep 6, 2019• 13 citations
arXivHTMLPDF

Abstract

Both insulating and conducting electronic behaviors have been experimentally seen in clean bilayer graphene samples at low temperature, and there is still no consensus on the nature of the interacting ground state at half filling and in the absence of a magnetic field. Theoretically, several possibilities for the insulating ground states have been predicted for weak interaction strength. However, a recent renormalization-group calculation on a Hubbard model for charge-neutral bilayer graphene with short-range interactions suggests the emergence of low-energy Dirac fermions that would stabilize the metallic phase for weak interactions. Using a nonperturbative projective quantum Monte Carlo, we calculate the ground state for bilayer graphene using a realistic model for the Coulomb interaction that includes both short-range and long-range contributions. We find that a finite critical on-site interaction is needed to gap bilayer graphene and the transition belongs to the Gross-Neveu universality class, thereby confirming the Hubbard model expectations even in the presence of a long-range Coulomb potential, in agreement with our theoretical analysis. In addition, we also find that the critical on-site interactions necessary to destabilize the metallic ground state decrease with increasing interlayer coupling.

DOI
10.1103/physrevb.100.125116

Related papers

Quantum Monte Carlo study of magnetic ordering and superconducting pairing symmetry in twisted bilayer graphene

Shi-Chao Fang, Guangkun Liu, Hai-Qing Lin +1 more

To understand correlated insulating and unconventional superconducting states in twisted bilayer graphene, we perform a systematic study of spin and pairing correlations in an effective two-orbital Hubbard model, by using the ground-state constrained-path quantum Monte Carlo method. Our numerical simulations reveal that when only the on-site Hubbard $U$ is considered, the long-range magnetic order develops for $U\ensuremath{\ge}3.2t$ at half filling. Upon doping away from half filling, an analysis based on the pairing correlations and corresponding vertex contributions identifies that pairing with the $d+id$ symmetry is the dominant pairing channel. As the on-site Coulomb interaction is increased, both spin correlation and effective pairing interaction are enhanced simultaneously, indicating the intimate relation between magnetism and superconductivity. An inclusion of nearest neighbor Coulomb interaction $V$ results in a suppression of the $d+id$ superconductivity, but it can survive in the small $V$ regime. Our findings are useful for clarifying the ongoing controversy on the superconducting pairing symmetry in twisted bilayer graphene.

Physical review. B./Physical review. B•Sep 16, 2019•16 citations

Interaction induced Dirac fermions from quadratic band touching in models of bilayer graphene

Sumiran Pujari, Thomas C. Lang, Ganpathy N. Murthy +1 more

We revisit the effect of local interactions on the quadratic band touching (QBT) of Bernal stacked bilayer graphene models using renormalization group (RG) arguments and quantum Monte Carlo simulations of the Hubbard model. We present an RG argument which predicts, contrary to previous studies, that weak interactions do not flow to strong coupling even if the free dispersion has a QBT. Instead they generate a linear term in the dispersion, which causes the interactions to flow back to weak coupling. Consistent with this RG scenario, in unbiased quantum Monte Carlo simulations of the Hubbard model we find compelling evidence that antiferromagnetism turns on at a finite $U/t$, despite the $U=0$ hopping problem having a QBT. The onset of antiferromagnetism takes place at a continuous transition which is consistent with a dynamical critical exponent $z=1$ as expected for 2+1 d Gross-Neveu criticality. We conclude that generically in models of bilayer graphene, even if the free dispersion has a QBT, small local interactions generate a Dirac phase with no symmetry breaking and that there is a finite-coupling transition out of this phase to a symmetry-broken state.

•Apr 13, 2016•0 citations

Valence Bond Orders at Charge Neutrality in a Possible Two-Orbital Extended Hubbard Model for Twisted Bilayer Graphene

Yuan Da Liao, Zi Yang Meng, Xiao Yan Xu

An extended Hubbard model on a honeycomb lattice with two orbitals per site at charge neutrality is investigated with unbiased large-scale quantum Monte Carlo simulations. The Fermi velocity of the Dirac fermions is renormalized as the cluster charge interaction increases, until a mass term emerges and a quantum phase transition from Dirac semimetal to valence bond solid (VBS) insulator is established. The quantum critical point is discovered to belong to the 3D N=4 Gross-Neveu chiral XY universality with the critical exponents obtained at high precision. Further enhancement of the interaction drives the system into two different VBS phases, the properties and transition between them are also revealed. Since the model is related to magic-angle twisted bilayer graphene, our results may have relevance towards the symmetry breaking order at the charge neutrality point of the material, and associate the wide range of universal strange metal behavior around it with quantum critical fluctuations.

arXiv:1901.11424•Physical Review Letters•Oct 8, 2019•89 citations

Fermionic Monte Carlo Study of a Realistic Model of Twisted Bilayer Graphene

Johannes Stephan Hofmann, Eslam Khalaf, Ashvin Vishwanath +2 more

The rich phenomenology of twisted bilayer graphene (TBG) near the magic angle is believed to arise from electron correlations in topological flat bands.An unbiased approach to this problem is highly desirable, but also particularly challenging, given the multiple electron flavors, the topological obstruction to defining tight-binding models, and the long-ranged Coulomb interactions.While numerical simulations of realistic models have thus far been confined to zero temperature, typically excluding some spin or valley species, analytic progress has relied on fixed point models away from the realistic limit.Here, we present unbiased Monte Carlo simulations of realistic models of magic-angle TBG at charge neutrality.We establish the absence of a sign problem for this model in a momentumspace approach and describe a computationally tractable formulation that applies even on breaking chiral symmetry and including band dispersion.Our results include (i) the emergence of an insulating Kramers intervalley coherent ground state in competition with a correlated semimetal phase, (ii) detailed temperature evolution of order parameters and electronic spectral functions that reveal a "pseudogap" regime, in which gap features are established at a higher temperature than the onset of order, and (iii) predictions for electronic tunneling spectra and their evolution with temperature.Our results pave the way towards uncovering the physics of magic-angle graphene through exact simulations of over a hundred electrons across a wide temperature range.

arXiv:2105.12112•Physical Review X•Mar 30, 2022•73 citations

Antiferromagnetism and competing charge instabilities of electrons in strained graphene from Coulomb interactions

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

We study the quantum many-body ground states of electrons on the half-filled honeycomb lattice with short- and long-ranged density-density interactions as a model for graphene. To this end, we employ the recently developed truncated-unity functional renormalization group (TU-fRG) approach which allows for a high resolution of the interaction vertex' wave vector dependence. We connect to previous lattice quantum Monte Carlo (QMC) results which predict a stabilization of the semimetallic phase for realistic ab initio interaction parameters and confirm that the application of a finite biaxial strain can induce a quantum phase transition towards an ordered ground state. In contrast to lattice QMC simulations, the TU-fRG is not limited in the choice of tight-binding and interaction parameters to avoid the occurrence of a sign problem. Therefore, we also investigate a range of parameters relevant to the realistic graphene material which are not accessible by numerically exact methods. Although a plethora of charge-density waves arises under medium-range interactions, we find the antiferromagnetic spin-density wave to be the prevailing instability for long-range interactions. We further explore the impact of an extended tight-binding Hamiltonian with second-nearest-neighbor hopping and a finite chemical potential for a more accurate description of the band structure of graphene's ${p}_{z}$ electrons.

arXiv:1709.02289•Physical review. B./Physical review. B•Nov 29, 2017•25 citations

Angle-tuned Gross-Neveu quantum criticality in twisted bilayer graphene

Cheng Huang, Nikolaos Parthenios, Maksim Ulybyshev +4 more

Abstract Twisted bilayber graphene (TBG) has revealed fascinating quantum many-body states, and unstrained chargeneutral TBG at 1.08° is understood, both experimentally and theoretically , as a correlated insulator due to the interplay of long-range Coulomb interactions and the quantum metrics of flat bands. However, the fate of the ground state as one continuously tunes the twisted angle is still largely unknown. Here, theoretically, by employing a newly developed momentum-space continuous-field quantum Monte Carlo method that can fully take into account the long-range Coulomb interactions and quantum metrics of flat bands with system sizes that were not possible before , we show that charge-neutral TBG realizes an angle-tuned quantum phase transition with Gross-Neveu criticality at Θc ~ 1.20(1)°, from the analyses of single-particle spectra, free energy and the order parameter of inter-valley coherence. This shows the exciting possibility that by tuning away from 1.08°, pristine TBG can mimic a fundamental phase transition of Dirac fermions.

arXiv:2412.11382•Nature Communications•Aug 4, 2025•15 citations
Both insulating and conducting electronic behaviors have been experimentally seen in clean bilayer graphene samples at low temperature, and there is still no consensus on the nature of the interacting ground state at half filling and in the absence of a magnetic field. Theoretically, several possibilities for the insulating ground states have been predicted for weak interaction strength. However, a recent renormalization-group calculation on a Hubbard model for charge-neutral bilayer graphene with short-range interactions suggests the emergence of low-energy Dirac fermions that would stabilize the metallic phase for weak interactions. Using a nonperturbative projective quantum Monte Carlo, we calculate the ground state for bilayer graphene using a realistic model for the Coulomb interaction that includes both short-range and long-range contributions. We find that a finite critical on-site interaction is needed to gap bilayer graphene and the transition belongs to the Gross-Neveu universality class, thereby confirming the Hubbard model expectations even in the presence of a long-range Coulomb potential, in agreement with our theoretical analysis. In addition, we also find that the critical on-site interactions necessary to destabilize the metallic ground state decrease with increasing interlayer coupling.
To understand correlated insulating and unconventional superconducting states in twisted bilayer graphene, we perform a systematic study of spin and pairing correlations in an effective two-orbital Hubbard model, by using the ground-state constrained-path quantum Monte Carlo method. Our numerical simulations reveal that when only the on-site Hubbard UUU is considered, the long-range magnetic order develops for $U\ensuremath{\ge}3.2t$ at half filling. Upon doping away from half filling, an analysis based on the pairing correlations and corresponding vertex contributions identifies that pairing with the d+idd+idd+id symmetry is the dominant pairing channel. As the on-site Coulomb interaction is increased, both spin correlation and effective pairing interaction are enhanced simultaneously, indicating the intimate relation between magnetism and superconductivity. An inclusion of nearest neighbor Coulomb interaction VVV results in a suppression of the d+idd+idd+id superconductivity, but it can survive in the small VVV regime. Our findings are useful for clarifying the ongoing controversy on the superconducting pairing symmetry in twisted bilayer graphene.
We revisit the effect of local interactions on the quadratic band touching (QBT) of Bernal stacked bilayer graphene models using renormalization group (RG) arguments and quantum Monte Carlo simulations of the Hubbard model. We present an RG argument which predicts, contrary to previous studies, that weak interactions do not flow to strong coupling even if the free dispersion has a QBT. Instead they generate a linear term in the dispersion, which causes the interactions to flow back to weak coupling. Consistent with this RG scenario, in unbiased quantum Monte Carlo simulations of the Hubbard model we find compelling evidence that antiferromagnetism turns on at a finite U/tU/tU/t, despite the U=0U=0U=0 hopping problem having a QBT. The onset of antiferromagnetism takes place at a continuous transition which is consistent with a dynamical critical exponent z=1z=1z=1 as expected for 2+1 d Gross-Neveu criticality. We conclude that generically in models of bilayer graphene, even if the free dispersion has a QBT, small local interactions generate a Dirac phase with no symmetry breaking and that there is a finite-coupling transition out of this phase to a symmetry-broken state.
An extended Hubbard model on a honeycomb lattice with two orbitals per site at charge neutrality is investigated with unbiased large-scale quantum Monte Carlo simulations. The Fermi velocity of the Dirac fermions is renormalized as the cluster charge interaction increases, until a mass term emerges and a quantum phase transition from Dirac semimetal to valence bond solid (VBS) insulator is established. The quantum critical point is discovered to belong to the 3D N=4 Gross-Neveu chiral XY universality with the critical exponents obtained at high precision. Further enhancement of the interaction drives the system into two different VBS phases, the properties and transition between them are also revealed. Since the model is related to magic-angle twisted bilayer graphene, our results may have relevance towards the symmetry breaking order at the charge neutrality point of the material, and associate the wide range of universal strange metal behavior around it with quantum critical fluctuations.
The rich phenomenology of twisted bilayer graphene (TBG) near the magic angle is believed to arise from electron correlations in topological flat bands.An unbiased approach to this problem is highly desirable, but also particularly challenging, given the multiple electron flavors, the topological obstruction to defining tight-binding models, and the long-ranged Coulomb interactions.While numerical simulations of realistic models have thus far been confined to zero temperature, typically excluding some spin or valley species, analytic progress has relied on fixed point models away from the realistic limit.Here, we present unbiased Monte Carlo simulations of realistic models of magic-angle TBG at charge neutrality.We establish the absence of a sign problem for this model in a momentumspace approach and describe a computationally tractable formulation that applies even on breaking chiral symmetry and including band dispersion.Our results include (i) the emergence of an insulating Kramers intervalley coherent ground state in competition with a correlated semimetal phase, (ii) detailed temperature evolution of order parameters and electronic spectral functions that reveal a "pseudogap" regime, in which gap features are established at a higher temperature than the onset of order, and (iii) predictions for electronic tunneling spectra and their evolution with temperature.Our results pave the way towards uncovering the physics of magic-angle graphene through exact simulations of over a hundred electrons across a wide temperature range.
We study the quantum many-body ground states of electrons on the half-filled honeycomb lattice with short- and long-ranged density-density interactions as a model for graphene. To this end, we employ the recently developed truncated-unity functional renormalization group (TU-fRG) approach which allows for a high resolution of the interaction vertex' wave vector dependence. We connect to previous lattice quantum Monte Carlo (QMC) results which predict a stabilization of the semimetallic phase for realistic ab initio interaction parameters and confirm that the application of a finite biaxial strain can induce a quantum phase transition towards an ordered ground state. In contrast to lattice QMC simulations, the TU-fRG is not limited in the choice of tight-binding and interaction parameters to avoid the occurrence of a sign problem. Therefore, we also investigate a range of parameters relevant to the realistic graphene material which are not accessible by numerically exact methods. Although a plethora of charge-density waves arises under medium-range interactions, we find the antiferromagnetic spin-density wave to be the prevailing instability for long-range interactions. We further explore the impact of an extended tight-binding Hamiltonian with second-nearest-neighbor hopping and a finite chemical potential for a more accurate description of the band structure of graphene's pz{p}_{z}pz​ electrons.
Abstract Twisted bilayber graphene (TBG) has revealed fascinating quantum many-body states, and unstrained chargeneutral TBG at 1.08° is understood, both experimentally and theoretically , as a correlated insulator due to the interplay of long-range Coulomb interactions and the quantum metrics of flat bands. However, the fate of the ground state as one continuously tunes the twisted angle is still largely unknown. Here, theoretically, by employing a newly developed momentum-space continuous-field quantum Monte Carlo method that can fully take into account the long-range Coulomb interactions and quantum metrics of flat bands with system sizes that were not possible before , we show that charge-neutral TBG realizes an angle-tuned quantum phase transition with Gross-Neveu criticality at Θc ~ 1.20(1)°, from the analyses of single-particle spectra, free energy and the order parameter of inter-valley coherence. This shows the exciting possibility that by tuning away from 1.08°, pristine TBG can mimic a fundamental phase transition of Dirac fermions.