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

Gross-Neveu-XY quantum criticality in moiré Dirac materials

Bilal Hawashin, Michael M. Scherer, Lukas Janssen

Published May 21, 2025• 7 citations
arXivHTMLPDF

Abstract

Two-dimensional van der Waals materials offer a highly tunable platform for engineering electronic band structures and interactions. By employing techniques such as twisting, gating, or applying pressure, these systems enable precise control over the electronic excitation spectrum. In moir\'e bilayer graphene, the tunability facilitates the transition from a symmetric Dirac semimetal phase through a quantum critical point into an interaction-induced long-range ordered phase with a finite band gap. At charge neutrality, the ordered state proposed to emerge from twist-angle tuning is the Kramers intervalley-coherent insulator. In this case, the transition falls into the quantum universality class of the relativistic Gross-Neveu-XY model in $2+1$ dimensions. Here, we refine estimates for the critical exponents characterizing this universality class using an expansion around the lower critical space-time dimension of two. We compute the order-parameter anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\varphi}}$ and the correlation-length exponent $\ensuremath{\nu}$ at one-loop order, and the fermion anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\psi}}$ at two-loop order. Combining these results with previous findings from the expansion around the upper critical dimension, we obtain improved estimates for the universal exponents in $2+1$ dimensions via Pad\'e interpolation. For ${N}_{\text{f}}=4$ four-component Dirac fermions, relevant to moir\'e bilayer graphene, we estimate $1/\ensuremath{\nu}=0.916(5), {\ensuremath{\eta}}_{\ensuremath{\varphi}}=0.926(13)$, and ${\ensuremath{\eta}}_{\ensuremath{\psi}}=0.0404(13)$. For ${N}_{\text{f}}=2$, potentially relevant to recent tetralayer ${\mathrm{WSe}}_{2}$ experiments, the Gross-Neveu-XY fixed point may be unstable owing to a fixed-point collision at ${N}_{\text{f,c}}$, with ${N}_{\text{f,c}}=1+\sqrt{2}+\mathcal{O}(\ensuremath{\epsilon})$ in the expansion around the lower critical dimension.

DOI
10.1103/physrevb.111.205129

Related papers

Gross-Neveu-Heisenberg criticality from 2 + ε expansion

Konstantinos Ladovrechis, Shouryya Ray, Tobias Meng +1 more

The Gross-Neveu-Heisenberg universality class describes a continuous quantum phase transition between a Dirac semimetal and an antiferromagnetic insulator. Such quantum critical points have originally been discussed in the context of Hubbard models on $\ensuremath{\pi}$-flux and honeycomb lattices, but more recently also in Bernal-stacked bilayer models, of potential relevance for bilayer graphene. Here, we demonstrate how the critical behavior of this fermionic universality class can be computed within an $\ensuremath{\epsilon}$ expansion around the lower critical space-time dimension of two. This approach is complementary to the previously studied expansion around the upper critical dimension of four. The crucial technical difference near the lower critical dimension is the presence of different four-fermion interaction channels at the critical point, which we take into account in a Fierz-complete way. By interpolating between the lower and upper critical dimensions, we obtain improved estimates for the critical exponents in $2+1$ space-time dimensions. For the situation relevant to single-layer graphene, we find an unusually small leading-correction-to-scaling exponent, arising from the competition between different interaction channels. This suggests that corrections to scaling may need to be taken into account when comparing analytical estimates with numerical data from finite-size extrapolations.

arXiv:2209.02734•Physical review. B./Physical review. B•Jan 30, 2023•20 citations

Fractionalized quantum criticality in spin-orbital liquids from field theory beyond the leading order

Shouryya Ray, Bernhard Ihrig, Daniel Kruti +3 more

Two-dimensional spin-orbital magnets with strong exchange frustration have recently been predicted to facilitate the realization of a quantum critical point in the Gross-Neveu-SO(3) universality class. In contrast to previously known Gross-Neveu-type universality classes, this quantum critical point separates a Dirac semimetal and a long-range-ordered phase, in which the fermion spectrum is only partially gapped out. Here, we characterize the quantum critical behavior of the Gross-Neveu-SO(3) universality class by employing three complementary field-theoretical techniques beyond their leading orders. We compute the correlation-length exponent $\ensuremath{\nu}$, the order-parameter anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\phi}}$, and the fermion anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\psi}}$ using a three-loop $\ensuremath{\epsilon}$ expansion around the upper critical space-time dimension of four, a second-order large-$N$ expansion (with the fermion anomalous dimension obtained even at the third order), as well as a functional renormalization group approach in the improved local potential approximation. For the physically relevant case of $N=3$ flavors of two-component Dirac fermions in $2+1$ space-time dimensions, we obtain the estimates $1/\ensuremath{\nu}=1.03(15), {\ensuremath{\eta}}_{\ensuremath{\phi}}=0.42(7)$, and ${\ensuremath{\eta}}_{\ensuremath{\psi}}=0.180(10)$ from averaging over the results of the different techniques, with the displayed uncertainty representing the degree of consistency among the three methods.

arXiv:2101.10335•Physical review. B./Physical review. B•Apr 30, 2021•20 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

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

Four-loop critical exponents for the Gross-Neveu-Yukawa models

Nikolai Zerf, Luminita N. Mihaila, Peter Marquard +2 more

We study the chiral Ising, the chiral XY, and the chiral Heisenberg models at four-loop order with the perturbative renormalization group in $4\ensuremath{-}\ensuremath{\epsilon}$ dimensions and compute critical exponents for the Gross-Neveu-Yukawa fixed points to order $\mathcal{O}({\ensuremath{\epsilon}}^{4})$. Further, we provide Pad\'e estimates for the correlation length exponent, the boson and fermion anomalous dimension, as well as the leading correction to scaling exponent in $2+1$ dimensions. We also confirm the emergence of supersymmetric field theories at four loops for the chiral Ising and the chiral XY models with $N=1/4$ and $N=1/2$ fermions, respectively. Furthermore, applications of our results relevant to various quantum transitions in the context of Dirac and Weyl semimetals are discussed, including interaction-induced transitions in graphene and surface states of topological insulators.

arXiv:1709.05057•Physical review. D/Physical review. D.•Nov 14, 2017•169 citations

Wilson-Fisher fixed points in the presence of Dirac fermions

Igor F. Herbut

Wilson–Fisher expansion near upper critical dimension has proven to be an invaluable conceptual and computational tool in our understanding of the universal critical behavior in the [Formula: see text] field theories that describe low-energy physics of the canonical models, such as Ising, XY, and Heisenberg. Here, I review its application to a class of the Gross–Neveu–Yukawa (GNY) field theories, which emerge as possible universal description of a number of quantum phase transitions in electronic two-dimensional systems such as graphene and d-wave superconductors. GNY field theories may be viewed as minimal modifications of the [Formula: see text] field theories in which the order parameter is coupled to relativistic Dirac fermions through Yukawa term and which still exhibit critical fixed points in the suitably formulated Wilson–Fisher [Formula: see text]-expansion. I discuss the unified GNY field theory for a set of different symmetry-breaking patterns, with focus on the semimetal-Néel-ordered-Mott insulator quantum phase transition in the half-filled Hubbard model on the honeycomb lattice, for which a comparison between the state-of-the-art [Formula: see text]-expansion, quantum Monte Carlo, large N, and functional renormalization-group calculations can be made.

arXiv:2304.07654•Modern Physics Letters B•Oct 8, 2024•16 citations
Two-dimensional van der Waals materials offer a highly tunable platform for engineering electronic band structures and interactions. By employing techniques such as twisting, gating, or applying pressure, these systems enable precise control over the electronic excitation spectrum. In moir'e bilayer graphene, the tunability facilitates the transition from a symmetric Dirac semimetal phase through a quantum critical point into an interaction-induced long-range ordered phase with a finite band gap. At charge neutrality, the ordered state proposed to emerge from twist-angle tuning is the Kramers intervalley-coherent insulator. In this case, the transition falls into the quantum universality class of the relativistic Gross-Neveu-XY model in 2+12+12+1 dimensions. Here, we refine estimates for the critical exponents characterizing this universality class using an expansion around the lower critical space-time dimension of two. We compute the order-parameter anomalous dimension \ensuremathη\ensuremathφ{\ensuremath{\eta}}_{\ensuremath{\varphi}}\ensuremathη\ensuremathφ​ and the correlation-length exponent \ensuremathν\ensuremath{\nu}\ensuremathν at one-loop order, and the fermion anomalous dimension \ensuremathη\ensuremathψ{\ensuremath{\eta}}_{\ensuremath{\psi}}\ensuremathη\ensuremathψ​ at two-loop order. Combining these results with previous findings from the expansion around the upper critical dimension, we obtain improved estimates for the universal exponents in 2+12+12+1 dimensions via Pad'e interpolation. For Nf=4{N}_{\text{f}}=4Nf​=4 four-component Dirac fermions, relevant to moir'e bilayer graphene, we estimate 1/\ensuremathν=0.916(5),\ensuremathη\ensuremathφ=0.926(13)1/\ensuremath{\nu}=0.916(5), {\ensuremath{\eta}}_{\ensuremath{\varphi}}=0.926(13)1/\ensuremathν=0.916(5),\ensuremathη\ensuremathφ​=0.926(13), and \ensuremathη\ensuremathψ=0.0404(13){\ensuremath{\eta}}_{\ensuremath{\psi}}=0.0404(13)\ensuremathη\ensuremathψ​=0.0404(13). For Nf=2{N}_{\text{f}}=2Nf​=2, potentially relevant to recent tetralayer WSe2{\mathrm{WSe}}_{2}WSe2​ experiments, the Gross-Neveu-XY fixed point may be unstable owing to a fixed-point collision at Nf,c{N}_{\text{f,c}}Nf,c​, with Nf,c=1+2+O(\ensuremathϵ){N}_{\text{f,c}}=1+\sqrt{2}+\mathcal{O}(\ensuremath{\epsilon})Nf,c​=1+2​+O(\ensuremathϵ) in the expansion around the lower critical dimension.
The Gross-Neveu-Heisenberg universality class describes a continuous quantum phase transition between a Dirac semimetal and an antiferromagnetic insulator. Such quantum critical points have originally been discussed in the context of Hubbard models on $\ensuremath{\pi}$-flux and honeycomb lattices, but more recently also in Bernal-stacked bilayer models, of potential relevance for bilayer graphene. Here, we demonstrate how the critical behavior of this fermionic universality class can be computed within an $\ensuremath{\epsilon}$ expansion around the lower critical space-time dimension of two. This approach is complementary to the previously studied expansion around the upper critical dimension of four. The crucial technical difference near the lower critical dimension is the presence of different four-fermion interaction channels at the critical point, which we take into account in a Fierz-complete way. By interpolating between the lower and upper critical dimensions, we obtain improved estimates for the critical exponents in 2+12+12+1 space-time dimensions. For the situation relevant to single-layer graphene, we find an unusually small leading-correction-to-scaling exponent, arising from the competition between different interaction channels. This suggests that corrections to scaling may need to be taken into account when comparing analytical estimates with numerical data from finite-size extrapolations.
Two-dimensional spin-orbital magnets with strong exchange frustration have recently been predicted to facilitate the realization of a quantum critical point in the Gross-Neveu-SO(3) universality class. In contrast to previously known Gross-Neveu-type universality classes, this quantum critical point separates a Dirac semimetal and a long-range-ordered phase, in which the fermion spectrum is only partially gapped out. Here, we characterize the quantum critical behavior of the Gross-Neveu-SO(3) universality class by employing three complementary field-theoretical techniques beyond their leading orders. We compute the correlation-length exponent $\ensuremath{\nu}$, the order-parameter anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\phi}}$, and the fermion anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\psi}}$ using a three-loop $\ensuremath{\epsilon}$ expansion around the upper critical space-time dimension of four, a second-order large-NNN expansion (with the fermion anomalous dimension obtained even at the third order), as well as a functional renormalization group approach in the improved local potential approximation. For the physically relevant case of N=3N=3N=3 flavors of two-component Dirac fermions in 2+12+12+1 space-time dimensions, we obtain the estimates $1/\ensuremath{\nu}=1.03(15), {\ensuremath{\eta}}_{\ensuremath{\phi}}=0.42(7)$, and ${\ensuremath{\eta}}_{\ensuremath{\psi}}=0.180(10)$ from averaging over the results of the different techniques, with the displayed uncertainty representing the degree of consistency among the three methods.
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.
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.
We study the chiral Ising, the chiral XY, and the chiral Heisenberg models at four-loop order with the perturbative renormalization group in $4\ensuremath{-}\ensuremath{\epsilon}$ dimensions and compute critical exponents for the Gross-Neveu-Yukawa fixed points to order $\mathcal{O}({\ensuremath{\epsilon}}^{4})$. Further, we provide Pad\'e estimates for the correlation length exponent, the boson and fermion anomalous dimension, as well as the leading correction to scaling exponent in 2+12+12+1 dimensions. We also confirm the emergence of supersymmetric field theories at four loops for the chiral Ising and the chiral XY models with N=1/4N=1/4N=1/4 and N=1/2N=1/2N=1/2 fermions, respectively. Furthermore, applications of our results relevant to various quantum transitions in the context of Dirac and Weyl semimetals are discussed, including interaction-induced transitions in graphene and surface states of topological insulators.
Wilson–Fisher expansion near upper critical dimension has proven to be an invaluable conceptual and computational tool in our understanding of the universal critical behavior in the [Formula: see text] field theories that describe low-energy physics of the canonical models, such as Ising, XY, and Heisenberg. Here, I review its application to a class of the Gross–Neveu–Yukawa (GNY) field theories, which emerge as possible universal description of a number of quantum phase transitions in electronic two-dimensional systems such as graphene and d-wave superconductors. GNY field theories may be viewed as minimal modifications of the [Formula: see text] field theories in which the order parameter is coupled to relativistic Dirac fermions through Yukawa term and which still exhibit critical fixed points in the suitably formulated Wilson–Fisher [Formula: see text]-expansion. I discuss the unified GNY field theory for a set of different symmetry-breaking patterns, with focus on the semimetal-Néel-ordered-Mott insulator quantum phase transition in the half-filled Hubbard model on the honeycomb lattice, for which a comparison between the state-of-the-art [Formula: see text]-expansion, quantum Monte Carlo, large N, and functional renormalization-group calculations can be made.