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

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

Yuzhi Liu, Wei Wang, Kai Sun, Zi Yang Meng

Published Feb 26, 2020• 41 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.101.064308

Related papers

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

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.

arXiv:2304.00034•Physical review. B./Physical review. B•Sep 22, 2023•13 citations

Fermion-enhanced first-order phase transition and chiral Gross-Neveu tricritical point

Yuzhi Liu, Zi Yang Meng, Shuai Yin

The fluctuations of massless Dirac fermion can not only turn a first-order bosonic phase transition (in the Landau sense) to a quantum critical point, but also work reversely to enhance the first-order transition itself, depending on the implementation of finite-size effects in the coupling corrections. Here, we report a case study of the latter by employing quantum Monte Carlo simulation upon a lattice model in which the bosonic part featuring the Landau-Devonshire first-order phase transition and Yukawa coupled to the Dirac fermions. We find that the parameter range for the first-order phase transition becomes larger as the Yukawa coupling increases, and the microscopic mechanism of this phenomena is revealed, at a quantitative level, as the interplay between the critical fluctuations and the finite-size effects. Moreover, the scaling behavior at the separation point between the first-order and the continuous phase transitions is found to belong to the chiral tricritical Gross-Neveu universality. Our results demonstrate that the interplay of massless Dirac fermions, critical fluctuations, and the finite-size effects could trigger a plethora of interesting phenomena, and therefore great care is called for when making generalizations.

arXiv:2012.00449•Physical review. B./Physical review. B•Feb 26, 2021•3 citations

Fermion Disorder Operator at Gross-Neveu and Deconfined Quantum Criticalities

Zi Hong Liu, Weilun Jiang, Bin-Bin Chen +5 more

The fermion disorder operator has been shown to reveal the entanglement information in 1D Luttinger liquids and 2D free and interacting Fermi and non-Fermi liquids emerging at quantum critical points (QCPs) [W. Jiang et al., arXiv:2209.07103]. Here we study, by means of large-scale quantum Monte Carlo simulation, the scaling behavior of the disorder operator in correlated Dirac systems. We first demonstrate the logarithmic scaling behavior of the disorder operator at the Gross-Neveu (GN) chiral Ising and Heisenberg QCPs, where consistent conformal field theory (CFT) content of the GN-QCP in its coefficient is found. Then we study a 2D monopole-free deconfined quantum critical point (DQCP) realized between a quantum-spin Hall insulator and a superconductor. Our data point to negative values of the logarithmic coefficients such that the DQCP does not correspond to a unitary CFT. Density matrix renormalization group calculations of the disorder operator on a 1D DQCP model also detect emergent continuous symmetries.

arXiv:2212.11821•Physical Review Letters•Jun 28, 2023•37 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

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

Random-mass disorder in the critical Gross-Neveu-Yukawa models

Hennadii Yerzhakov, Joseph Maciejko

An important yet largely unsolved problem in the statistical mechanics of disordered quantum systems is to understand how quenched disorder affects quantum phase transitions in systems of itinerant fermions. In the clean limit, continuous quantum phase transitions of the symmetry-breaking type in Dirac materials such as graphene and the surfaces of topological insulators are described by relativistic (2+1)-dimensional quantum field theories of the Gross-Neveu-Yukawa (GNY) type. We study the universal critical properties of the chiral Ising, XY, and Heisenberg GNY models perturbed by quenched random-mass disorder, both uncorrelated or with long-range power-law correlations. Using the replica method combined with a controlled triple epsilon expansion below four dimensions, we find a variety of new finite-randomness critical and multicritical points with nonzero Yukawa coupling between low-energy Dirac fields and bosonic order parameter fluctuations, and compute their universal critical exponents. Analyzing bifurcations of the renormalization-group flow, we find instances of the fixed-point annihilation scenario—continuously tuned by the power-law exponent of long-range disorder correlations and associated with an exponentially large crossover length—as well as the transcritical bifurcation and the supercritical Hopf bifurcation. The latter is accompanied by the birth of a stable limit cycle on the critical hypersurface, which represents the first instance of fermionic quantum criticality with emergent discrete scale invariance.

arXiv:2008.13663•Nuclear Physics B•Nov 17, 2020•12 citations
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×\else\texttimes\fi1616\ifmmode\times\else\texttimes\fi{}1616\ifmmode×\else\texttimes\fi16) 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.
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{\rtimes}{\mathbb{Z}}_{2}$ GNY N=8N=8N=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=8N=8N=8 island and could explore the chiral Gross-Neveu-Yukawa archipelago with ease.
The fluctuations of massless Dirac fermion can not only turn a first-order bosonic phase transition (in the Landau sense) to a quantum critical point, but also work reversely to enhance the first-order transition itself, depending on the implementation of finite-size effects in the coupling corrections. Here, we report a case study of the latter by employing quantum Monte Carlo simulation upon a lattice model in which the bosonic part featuring the Landau-Devonshire first-order phase transition and Yukawa coupled to the Dirac fermions. We find that the parameter range for the first-order phase transition becomes larger as the Yukawa coupling increases, and the microscopic mechanism of this phenomena is revealed, at a quantitative level, as the interplay between the critical fluctuations and the finite-size effects. Moreover, the scaling behavior at the separation point between the first-order and the continuous phase transitions is found to belong to the chiral tricritical Gross-Neveu universality. Our results demonstrate that the interplay of massless Dirac fermions, critical fluctuations, and the finite-size effects could trigger a plethora of interesting phenomena, and therefore great care is called for when making generalizations.
The fermion disorder operator has been shown to reveal the entanglement information in 1D Luttinger liquids and 2D free and interacting Fermi and non-Fermi liquids emerging at quantum critical points (QCPs) [W. Jiang et al., arXiv:2209.07103]. Here we study, by means of large-scale quantum Monte Carlo simulation, the scaling behavior of the disorder operator in correlated Dirac systems. We first demonstrate the logarithmic scaling behavior of the disorder operator at the Gross-Neveu (GN) chiral Ising and Heisenberg QCPs, where consistent conformal field theory (CFT) content of the GN-QCP in its coefficient is found. Then we study a 2D monopole-free deconfined quantum critical point (DQCP) realized between a quantum-spin Hall insulator and a superconductor. Our data point to negative values of the logarithmic coefficients such that the DQCP does not correspond to a unitary CFT. Density matrix renormalization group calculations of the disorder operator on a 1D DQCP model also detect emergent continuous symmetries.
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 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.
An important yet largely unsolved problem in the statistical mechanics of disordered quantum systems is to understand how quenched disorder affects quantum phase transitions in systems of itinerant fermions. In the clean limit, continuous quantum phase transitions of the symmetry-breaking type in Dirac materials such as graphene and the surfaces of topological insulators are described by relativistic (2+1)-dimensional quantum field theories of the Gross-Neveu-Yukawa (GNY) type. We study the universal critical properties of the chiral Ising, XY, and Heisenberg GNY models perturbed by quenched random-mass disorder, both uncorrelated or with long-range power-law correlations. Using the replica method combined with a controlled triple epsilon expansion below four dimensions, we find a variety of new finite-randomness critical and multicritical points with nonzero Yukawa coupling between low-energy Dirac fields and bosonic order parameter fluctuations, and compute their universal critical exponents. Analyzing bifurcations of the renormalization-group flow, we find instances of the fixed-point annihilation scenario—continuously tuned by the power-law exponent of long-range disorder correlations and associated with an exponentially large crossover length—as well as the transcritical bifurcation and the supercritical Hopf bifurcation. The latter is accompanied by the birth of a stable limit cycle on the critical hypersurface, which represents the first instance of fermionic quantum criticality with emergent discrete scale invariance.