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

Quantum Monte Carlo calculation of critical exponents of the Gross-Neveu-Yukawa on a two-dimensional fermion lattice model

Ting-Tung Wang, Zi Yang Meng

Published Sep 22, 2023• 13 citations
arXivHTMLPDF

Abstract

It is expected that the Gross-Neveu-Yukawa (GNY) chiral Ising transition of Dirac fermions coupled with a scalar field in $(2+1)$ dimensions will be the first fermionic quantum critical point that various methods, such as conformal bootstrap, perturbative renormalization group, and quantum Monte Carlo (QMC) simulations, would yield converged critical exponents---serving the same role as the Ising and $O(N)$ models in the textbooks of statistical and quantum physics. However, such an expectation has not been fully realized from the lattice QMC simulations due to the obstacles introduced by the UV finite-size effect. In this Letter, by means of the elective-momentum ultrasize (EMUS)-QMC method, we compute the critical exponents of the $O{(N/2)}^{2}\ensuremath{\rtimes}{\mathbb{Z}}_{2}$ GNY $N=8$ chiral Ising transition on a two-dimensional $\ensuremath{\pi}$-flux fermion lattice model between Dirac semimetal and quantum spin Hall insulator phases. With the matching of fermionic and bosonic momentum transfer and collective update in momentum space, our QMC results provide fully consistent exponents with those obtained from the bootstrap and perturbative approaches. In this way, the EMUS now live happily on the $N=8$ island and could explore the chiral Gross-Neveu-Yukawa archipelago with ease.

DOI
10.1103/physrevb.108.l121112

Related papers

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

Chiral Ising Gross-Neveu Criticality of a Single Dirac Cone: A Quantum Monte Carlo Study

S. Mojtaba Tabatabaei, Amir-Reza Negari, Joseph Maciejko +1 more

We perform large-scale quantum Monte Carlo simulations of SLAC fermions on a two-dimensional square lattice at half filling with a single Dirac cone with N=2 spinor components and repulsive on-site interactions. Despite the presence of a sign problem, we accurately identify the critical interaction strength U_{c}=7.28±0.02 in units of the hopping amplitude, for a continuous quantum phase transition between a paramagnetic Dirac semimetal and a ferromagnetic insulator. Using finite-size scaling, we extract the critical exponents for the corresponding N=2 chiral Ising Gross-Neveu universality class: the inverse correlation length exponent ν^{-1}=1.19±0.03, the order parameter anomalous dimension η_{ϕ}=0.31±0.01, and the fermion anomalous dimension η_{ψ}=0.136±0.005.

arXiv:2112.09209•Physical Review Letters•Jun 3, 2022•25 citations

Designer Monte Carlo simulation for the Gross-Neveu-Yukawa transition

Yuzhi Liu, Wei Wang, Kai Sun +1 more

In this paper, we study the quantum criticality of Dirac fermions via large-scale numerical simulations, focusing on the Gross-Neveu-Yukawa chiral-Ising quantum critical point (QCP) with critical bosonic modes coupled with Dirac fermions. We show that finite-size effects at this QCP can be efficiently minimized via model design, which maximizes the ultraviolet cutoff and at the same time places the bare control parameters closer to the nontrivial fixed point to better expose the critical region. Combined with the efficient self-learning quantum Monte Carlo algorithm, which enables a nonlocal update of the bosonic field, we find that moderately large system size (up to $16\ifmmode\times\else\texttimes\fi{}16$) is already sufficient to produce robust scaling behavior and critical exponents. The conductance of free Dirac fermions is also calculated, and its frequency dependence is found to be consistent with the scaling behavior predicted by the conformal field theory. The methods and model-design principles developed for this study can be generalized to other fermionic QCPs, and thus provide a promising direction for controlled studies of strongly correlated itinerant systems.

arXiv:1910.07430•Physical review. B./Physical review. B•Feb 26, 2020•41 citations

Quantum Monte Carlo Simulation of the Chiral Heisenberg Gross-Neveu-Yukawa Phase Transition with a Single Dirac Cone

Thomas C. Lang, Andreas Martin Läuchli

We present quantum Monte Carlo simulations for the chiral Heisenberg Gross-Neveu-Yukawa quantum phase transition of relativistic fermions with N=4 Dirac spinor components subject to a repulsive, local four fermion interaction in (2+1)D. Here we employ a two-dimensional lattice Hamiltonian with a single, spin-degenerate Dirac cone, which exactly reproduces a linear energy-momentum relation for all finite size lattice momenta in the absence of interactions. This allows us to significantly reduce finite size corrections compared to the widely studied honeycomb and π-flux lattices. A Hubbard term dynamically generates a mass beyond a critical coupling of U_{c}=6.76(1) as the system acquires antiferromagnetic order and SU(2) spin rotational symmetry is spontaneously broken. At the quantum phase transition, we extract a self-consistent set of critical exponents ν=0.98(1), η_{ϕ}=0.53(1), η_{ψ}=0.18(1), and β=0.75(1). We provide evidence for the continuous degradation of the quasiparticle weight of the fermionic excitations as the critical point is approached from the semimetallic phase. Finally, we study the effective "speed of light" of the low-energy relativistic description, which depends on the interaction U, but is expected to be regular across the quantum phase transition. We illustrate that the strongly coupled bosonic and fermionic excitations share a common velocity at the critical point.

arXiv:1808.01230•Physical Review Letters•Sep 25, 2019•73 citations

Boundary Criticality for the Gross-Neveu-Yukawa Models

Huan Jiang, Ge Yang, Shao-Kai Jian

We study the boundary criticality for the Gross-Neveu-Yukawa (GNY) models. Employing interacting Dirac fermions on a honeycomb lattice with armchair boundaries, we use determinant quantum Monte Carlo simulation to uncover rich boundary criticalities at the quantum phase transition to a charge density wave (CDW) insulator, including the ordinary, special, and extraordinary transitions. The Dirac fermions satisfy a Dirichlet boundary condition, while the boson field, representing the CDW order, obeys Dirichlet and Neumann conditions at the ordinary and special transitions, respectively, thereby enriching the critical GNY model. We develop a perturbative 4-ε renormalization group approach to compute the boundary critical exponents. Our framework generalizes to other GNY universality class variants and provides theoretical predictions for experiments.

arXiv:2503.13247•Physical Review Letters•Sep 10, 2025•3 citations

Fermion-sign-free Majarana-quantum-Monte-Carlo studies of quantum critical phenomena of Dirac fermions in two dimensions

Zixiang Li, Yi‐Fan Jiang, Hong Yao

Quantum critical phenomena may be qualitatively different when massless Dirac fermions are present at criticality. Using our recently-discovered fermion-sign-free Majorana quantum Monte Carlo method introduced by us in (Li et al 2015 Phys. Rev. B 91 241117 ), we investigate the quantum critical phenomena of spinless Dirac fermions at their charge-density-wave phase transitions on the honeycomb lattice having sites with largest L = 24. By finite-size scaling, we accurately obtain critical exponents of this so-called Gross–Neveu chiral-Ising universality class of two (two-component) Dirac fermions in 2+1D: , , and , which are qualitatively different from the mean-field results but are reasonably close to the ones obtained from renormalization group calculations.

arXiv:1411.7383•New Journal of Physics•Aug 6, 2015•104 citations
It is expected that the Gross-Neveu-Yukawa (GNY) chiral Ising transition of Dirac fermions coupled with a scalar field in (2+1)(2+1)(2+1) dimensions will be the first fermionic quantum critical point that various methods, such as conformal bootstrap, perturbative renormalization group, and quantum Monte Carlo (QMC) simulations, would yield converged critical exponents---serving the same role as the Ising and O(N)O(N)O(N) models in the textbooks of statistical and quantum physics. However, such an expectation has not been fully realized from the lattice QMC simulations due to the obstacles introduced by the UV finite-size effect. In this Letter, by means of the elective-momentum ultrasize (EMUS)-QMC method, we compute the critical exponents of the O(N/2)2\ensuremath⋊Z2O{(N/2)}^{2}\ensuremath{\rtimes}{\mathbb{Z}}_{2}O(N/2)2\ensuremath⋊Z2​ GNY N=8N=8N=8 chiral Ising transition on a two-dimensional \ensuremathπ\ensuremath{\pi}\ensuremathπ-flux fermion lattice model between Dirac semimetal and quantum spin Hall insulator phases. With the matching of fermionic and bosonic momentum transfer and collective update in momentum space, our QMC results provide fully consistent exponents with those obtained from the bootstrap and perturbative approaches. In this way, the EMUS now live happily on the N=8N=8N=8 island and could explore the chiral Gross-Neveu-Yukawa archipelago with ease.
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.
We perform large-scale quantum Monte Carlo simulations of SLAC fermions on a two-dimensional square lattice at half filling with a single Dirac cone with N=2 spinor components and repulsive on-site interactions. Despite the presence of a sign problem, we accurately identify the critical interaction strength U_{c}=7.28±0.02 in units of the hopping amplitude, for a continuous quantum phase transition between a paramagnetic Dirac semimetal and a ferromagnetic insulator. Using finite-size scaling, we extract the critical exponents for the corresponding N=2 chiral Ising Gross-Neveu universality class: the inverse correlation length exponent ν^{-1}=1.19±0.03, the order parameter anomalous dimension η_{ϕ}=0.31±0.01, and the fermion anomalous dimension η_{ψ}=0.136±0.005.
In this paper, we study the quantum criticality of Dirac fermions via large-scale numerical simulations, focusing on the Gross-Neveu-Yukawa chiral-Ising quantum critical point (QCP) with critical bosonic modes coupled with Dirac fermions. We show that finite-size effects at this QCP can be efficiently minimized via model design, which maximizes the ultraviolet cutoff and at the same time places the bare control parameters closer to the nontrivial fixed point to better expose the critical region. Combined with the efficient self-learning quantum Monte Carlo algorithm, which enables a nonlocal update of the bosonic field, we find that moderately large system size (up to $16\ifmmode\times\else\texttimes\fi{}16$) is already sufficient to produce robust scaling behavior and critical exponents. The conductance of free Dirac fermions is also calculated, and its frequency dependence is found to be consistent with the scaling behavior predicted by the conformal field theory. The methods and model-design principles developed for this study can be generalized to other fermionic QCPs, and thus provide a promising direction for controlled studies of strongly correlated itinerant systems.
We present quantum Monte Carlo simulations for the chiral Heisenberg Gross-Neveu-Yukawa quantum phase transition of relativistic fermions with N=4 Dirac spinor components subject to a repulsive, local four fermion interaction in (2+1)D. Here we employ a two-dimensional lattice Hamiltonian with a single, spin-degenerate Dirac cone, which exactly reproduces a linear energy-momentum relation for all finite size lattice momenta in the absence of interactions. This allows us to significantly reduce finite size corrections compared to the widely studied honeycomb and π-flux lattices. A Hubbard term dynamically generates a mass beyond a critical coupling of U_{c}=6.76(1) as the system acquires antiferromagnetic order and SU(2) spin rotational symmetry is spontaneously broken. At the quantum phase transition, we extract a self-consistent set of critical exponents ν=0.98(1), η_{ϕ}=0.53(1), η_{ψ}=0.18(1), and β=0.75(1). We provide evidence for the continuous degradation of the quasiparticle weight of the fermionic excitations as the critical point is approached from the semimetallic phase. Finally, we study the effective "speed of light" of the low-energy relativistic description, which depends on the interaction U, but is expected to be regular across the quantum phase transition. We illustrate that the strongly coupled bosonic and fermionic excitations share a common velocity at the critical point.
We study the boundary criticality for the Gross-Neveu-Yukawa (GNY) models. Employing interacting Dirac fermions on a honeycomb lattice with armchair boundaries, we use determinant quantum Monte Carlo simulation to uncover rich boundary criticalities at the quantum phase transition to a charge density wave (CDW) insulator, including the ordinary, special, and extraordinary transitions. The Dirac fermions satisfy a Dirichlet boundary condition, while the boson field, representing the CDW order, obeys Dirichlet and Neumann conditions at the ordinary and special transitions, respectively, thereby enriching the critical GNY model. We develop a perturbative 4-ε renormalization group approach to compute the boundary critical exponents. Our framework generalizes to other GNY universality class variants and provides theoretical predictions for experiments.
Quantum critical phenomena may be qualitatively different when massless Dirac fermions are present at criticality. Using our recently-discovered fermion-sign-free Majorana quantum Monte Carlo method introduced by us in (Li et al 2015 Phys. Rev. B 91 241117 ), we investigate the quantum critical phenomena of spinless Dirac fermions at their charge-density-wave phase transitions on the honeycomb lattice having sites with largest L = 24. By finite-size scaling, we accurately obtain critical exponents of this so-called Gross–Neveu chiral-Ising universality class of two (two-component) Dirac fermions in 2+1D: , , and , which are qualitatively different from the mean-field results but are reasonably close to the ones obtained from renormalization group calculations.