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

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

Sumiran Pujari, Thomas C. Lang, Ganpathy N. Murthy, Ribhu K. Kaul

Published Apr 13, 2016• 0 citations
OpenAlex

Abstract

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.

Related papers

Interaction-Induced Dirac Fermions from Quadratic Band Touching in 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 the Bernal honeycomb bilayer model using renormalization group (RG) arguments and quantum Monte Carlo (QMC) simulations. We present a 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 QMC 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 (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.

arXiv:1604.03876•Physical Review Letters•Aug 19, 2016•69 citations

Ubiquitous nematic Dirac semimetal emerging from interacting quadratic band touching systems

Hongyu Lu, Kai Sun, Zi Yang Meng +1 more

Quadratic band touching (QBT) points are widely observed in two- and three-dimensional (2D and 3D) materials, including bilayer graphene and Luttinger semimetals, and attract significant attention from theory to experiment. However, even in its simplest form, the 2D checkerboard lattice QBT model, the phase diagram characterized by temperature and interaction strength, still remains unknown beyond the weak-coupling regime. Intense debates persist regarding the existence of various interaction-driven insulating states in this system. To address these uncertainties, we employ thermal tensor network simulations, specifically exponential tensor renormalization group and tangent space tensor renormalization group, along with density matrix renormalization group calculations to provide a comprehensive finite-temperature phase diagram for this model and shed light on previous ambiguities. Notably, our findings reveal the emergence of a robust bond-nematic Dirac semimetal phase with distinct thermodynamic properties that set it part from the nematic insulating state and other symmetry-broken states. This previously overlooked feature is found to be ubiquitous in interacting QBT systems. We also discuss the implications of these results for experimental systems such as bilayer graphene and iridate compounds.

arXiv:2305.17189•Physical review. B./Physical review. B•Feb 12, 2024•5 citations

Interaction-induced nematic Dirac semimetal from quadratic band touching: A constrained-path quantum Monte Carlo study

Zi Hong Liu, Hongyu Lu, Zi Yang Meng +1 more

Electronic systems with quadratic band touchings, commonly found in two- and three-dimensional materials such as Bernal-stacked bilayer graphene, kagome metals, HgTe, and pyrochlore iridates, have attracted significant interest concerning the role of interactions in shaping their electronic properties. However, even in the simplest model of spinless fermions on a two-dimensional checkerboard lattice, the quantum phase diagram as a function of nearest-neighbor interaction remains under debate. We employ constrained-path quantum Monte Carlo simulations (CP-QMC) simulations to investigate the problem using a two-dimensional torus geometry. We cross-validate our results on small lattices by comparing them with density-matrix renormalization group calculations, finding quantitative agreement. In particular, we implement an improved optimization scheme within the CP-QMC simulations, enabling the identification of a bond-nematic Dirac semimetal phase that was found in tensor-network studies on cylindrical geometries, but remains inaccessible to Hartree-Fock mean-field methods. The CP-QMC approach makes it possible to establish the emergence of this phase in a geometry that preserves lattice rotational symmetry and permits extrapolation to the thermodynamic limit. Our results show that the quantum phase diagram of spinless fermions on the checkerboard lattice with nearest-neighbor repulsion features three interaction-induced phases at half filling: a quantum anomalous Hall insulator at weak coupling, a bond-nematic Dirac semimetal at intermediate coupling, and a site-nematic insulator at strong coupling.

arXiv:2507.15668•SciPost Physics•Jun 26, 2026•0 citations

Electronic ground state in bilayer graphene with realistic Coulomb interactions

Jia Ning Leaw, Ho Kin Tang, Pinaki Sengupta +3 more

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.

arXiv:1903.06177•Physical review. B./Physical review. B•Sep 6, 2019•13 citations

Quantum critical behavior of two-dimensional Fermi systems with quadratic band touching

Shouryya Ray, Matthias Vojta C, Lukas Janssen

At low temperatures, isotropic two-dimensional Fermi systems with quadratic band touching are unstable with respect to a symmetry-broken ground state when interactions are taken into account. Here, the authors demonstrate that in a setup with threefold symmetry only - as realized, for instance, in bilayer graphene - interactions induce a splitting of the quadratic band touching into four Dirac cones. This shifts the instability to a finite interaction strength and leads to a quantum critical point with emergent Lorentz invariance.

arXiv:1810.07695•Physical review. B./Physical review. B•Dec 19, 2018•27 citations

Emergent relativistic symmetry from interacting fermions on the honeycomb bilayer

Zi Hong Liu, Lukas Janssen

We study the phase diagram of interacting spinless fermions on the honeycomb bilayer at charge neutrality using large-scale quantum Monte Carlo simulations. In the noninteracting limit, the low-energy spectrum features quadratically dispersing bands that touch at the corners of the hexagonal Brillouin zone. Weak to intermediate interactions induce a splitting of each of the quadratic band touching points into four Dirac points, located along high-symmetry directions of the reciprocal lattice. Strong interactions lead to the formation of a layer-polarized charge density wave, which spontaneously breaks the $\mathbb Z_2$ layer inversion symmetry and opens an insulating gap in the spectrum. We show that the semimetal-to-insulator quantum phase transition as a function of interaction is continuous and characterized by emergent relativistic symmetry. Our results for the values of the correlation-length exponent $ν$, the order-parameter anomalous dimension $η_ϕ$, and the fermion anomalous dimension $η_ψ$ agree with those of the theoretically predicted 2+1D Gross-Neveu-Ising universality class with eight two-component Dirac fermions within less than 5\%\ deviation. We also determine the crossover scale as a function of interaction strength between the nonrelativistic semimetal state at high temperatures, characterized by dynamical critical exponent $z = 2$, and the Dirac semimetal state at intermediates temperatures, characterized by $z=1$. Further reducing the temperature below the crossover scale at a fixed value of the interaction strength above the quantum critical point results in a classical ordering transition in the 2D Ising universality class.

arXiv:2603.22259•Physical review. B./Physical review. B•Sep 22, 2026•0 citations
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.
We revisit the effect of local interactions on the quadratic band touching (QBT) of the Bernal honeycomb bilayer model using renormalization group (RG) arguments and quantum Monte Carlo (QMC) simulations. We present a 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 QMC 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 (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.
Quadratic band touching (QBT) points are widely observed in two- and three-dimensional (2D and 3D) materials, including bilayer graphene and Luttinger semimetals, and attract significant attention from theory to experiment. However, even in its simplest form, the 2D checkerboard lattice QBT model, the phase diagram characterized by temperature and interaction strength, still remains unknown beyond the weak-coupling regime. Intense debates persist regarding the existence of various interaction-driven insulating states in this system. To address these uncertainties, we employ thermal tensor network simulations, specifically exponential tensor renormalization group and tangent space tensor renormalization group, along with density matrix renormalization group calculations to provide a comprehensive finite-temperature phase diagram for this model and shed light on previous ambiguities. Notably, our findings reveal the emergence of a robust bond-nematic Dirac semimetal phase with distinct thermodynamic properties that set it part from the nematic insulating state and other symmetry-broken states. This previously overlooked feature is found to be ubiquitous in interacting QBT systems. We also discuss the implications of these results for experimental systems such as bilayer graphene and iridate compounds.
Electronic systems with quadratic band touchings, commonly found in two- and three-dimensional materials such as Bernal-stacked bilayer graphene, kagome metals, HgTe, and pyrochlore iridates, have attracted significant interest concerning the role of interactions in shaping their electronic properties. However, even in the simplest model of spinless fermions on a two-dimensional checkerboard lattice, the quantum phase diagram as a function of nearest-neighbor interaction remains under debate. We employ constrained-path quantum Monte Carlo simulations (CP-QMC) simulations to investigate the problem using a two-dimensional torus geometry. We cross-validate our results on small lattices by comparing them with density-matrix renormalization group calculations, finding quantitative agreement. In particular, we implement an improved optimization scheme within the CP-QMC simulations, enabling the identification of a bond-nematic Dirac semimetal phase that was found in tensor-network studies on cylindrical geometries, but remains inaccessible to Hartree-Fock mean-field methods. The CP-QMC approach makes it possible to establish the emergence of this phase in a geometry that preserves lattice rotational symmetry and permits extrapolation to the thermodynamic limit. Our results show that the quantum phase diagram of spinless fermions on the checkerboard lattice with nearest-neighbor repulsion features three interaction-induced phases at half filling: a quantum anomalous Hall insulator at weak coupling, a bond-nematic Dirac semimetal at intermediate coupling, and a site-nematic insulator at strong coupling.
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.
At low temperatures, isotropic two-dimensional Fermi systems with quadratic band touching are unstable with respect to a symmetry-broken ground state when interactions are taken into account. Here, the authors demonstrate that in a setup with threefold symmetry only - as realized, for instance, in bilayer graphene - interactions induce a splitting of the quadratic band touching into four Dirac cones. This shifts the instability to a finite interaction strength and leads to a quantum critical point with emergent Lorentz invariance.
We study the phase diagram of interacting spinless fermions on the honeycomb bilayer at charge neutrality using large-scale quantum Monte Carlo simulations. In the noninteracting limit, the low-energy spectrum features quadratically dispersing bands that touch at the corners of the hexagonal Brillouin zone. Weak to intermediate interactions induce a splitting of each of the quadratic band touching points into four Dirac points, located along high-symmetry directions of the reciprocal lattice. Strong interactions lead to the formation of a layer-polarized charge density wave, which spontaneously breaks the Z2\mathbb Z_2Z2​ layer inversion symmetry and opens an insulating gap in the spectrum. We show that the semimetal-to-insulator quantum phase transition as a function of interaction is continuous and characterized by emergent relativistic symmetry. Our results for the values of the correlation-length exponent ννν, the order-parameter anomalous dimension ηϕη_ϕηϕ​, and the fermion anomalous dimension ηψη_ψηψ​ agree with those of the theoretically predicted 2+1D Gross-Neveu-Ising universality class with eight two-component Dirac fermions within less than 5\%\ deviation. We also determine the crossover scale as a function of interaction strength between the nonrelativistic semimetal state at high temperatures, characterized by dynamical critical exponent z=2z = 2z=2, and the Dirac semimetal state at intermediates temperatures, characterized by z=1z=1z=1. Further reducing the temperature below the crossover scale at a fixed value of the interaction strength above the quantum critical point results in a classical ordering transition in the 2D Ising universality class.