Abstract We investigate the fate of the non-supersymmetric Gross-Neveu-Yukawa fixed point found by Fei et al. in 4 − ϵ dimensions with a two-component Majorana fermion continued to two dimensions. Assuming that it is a fermionic minimal model which possesses a chiral ℤ2 symmetry (in addition to fermion number parity) and just two relevant singlet operators, we can zero in on four candidates. Assuming further that the least relevant deformation leads to the supersymmetric Gross-Neveu-Yukawa fixed point (i.e. fermionic tricritical Ising model), we can rule out two of them by matching the spin contents of the preserved topological defect lines. The final candidates are the fermionic (11, 4) minimal model if it is non-unitary, and the fermionic (E6, A10) minimal model if it is unitary. If we further use a constraint from the double braiding relation proposed by one of the authors, the former scenario is preferable.
Abstract We investigate the fate of the non-supersymmetric Gross-Neveu-Yukawa fixed point found by Fei et al. in 4 − ϵ dimensions with a two-component Majorana fermion continued to two dimensions. Assuming that it is a fermionic minimal model which possesses a chiral ℤ2 symmetry (in addition to fermion number parity) and just two relevant singlet operators, we can zero in on four candidates. Assuming further that the least relevant deformation leads to the supersymmetric Gross-Neveu-Yukawa fixed point (i.e. fermionic tricritical Ising model), we can rule out two of them by matching the spin contents of the preserved topological defect lines. The final candidates are the fermionic (11, 4) minimal model if it is non-unitary, and the fermionic (E6, A10) minimal model if it is unitary. If we further use a constraint from the double braiding relation proposed by one of the authors, the former scenario is preferable.
The canonical Gross-Neveu model for N two-component Dirac fermions in 2+1 dimensions suffers a continuous phase transition at a critical interaction gc1∼1/N at large N , at which its continuous symmetry SO(2N) is preserved and a discrete (Ising) symmetry becomes spontaneously broken. A recent mean-field calculation, however, points to an additional transition at a different critical gc2∼−Ngc1 , at which SO(2N)→SO(N)×SO(N) . To study the latter phase transition we rewrite the Gross-Neveu interaction g(ψ¯ψ)2 in terms of three different quartic terms for the single ( L=1 ) 4N -component real (Majorana) fermion, and then extend the theory to L>1 . This allows us to track the evolution of the fixed points of the renormalization group transformation starting from L≫1 , where one can discern three distinct critical points which correspond to continuous phase transitions into (1) SO(2N) -singlet mass-order-parameter, (2) SO(2N) -symmetric-tensor mass-order-parameters, and (3) SO(2N) -adjoint nematic-order-parameters, down to L=1 value that is relevant to the standard Gross-Neveu model. Below the critical value of Lc(N)≈0.35N for N≫1 only the Gross-Neveu critical point (1) still implies a diverging susceptibility for its corresponding ( SO(2N) -singlet) order parameter, whereas the two new critical points that existed at large L ultimately become equivalent to the Gaussian fixed point at L=1 . We interpret this metamorphosis of the SO(2N) -symmetric-tensor fixed point from critical to spurious as an indication that the transition at gc2 in the original Gross-Neveu model is turned first-order by fluctuations. Published by the American Physical Society 2024
It was recently established that the paradigmatic Gross-Neveu model with N copies of two-dimensional Dirac fermions features an SO ( 2 N ) symmetry if certain interactions are suppressed. This becomes evident when the theory is rewritten in terms of 2 N copies of two-dimensional Majorana fermions. Mean-field theory for the SO ( 2 N ) model predicts, besides the chiral Ising transition at g c 1 , a second critical point g c 2 where SO ( 2 N ) is broken down to SO ( N ) × SO ( N ) . A subsequent Wilsonian renormalization group analysis directly in d = 3 supports its existence in a generalized theory, where N f copies of the 4 N -component Majorana fermions are introduced. This allows us to track the evolution of a (i) quantum anomalous Hall Gross-Neveu-Ising, (ii) symmetric-tensor, and (iii) adjoint-nematic fixed point separately. However, it turns out that (ii) and (iii) lose their criticality when approaching N f = 1 , suggesting that the transition is first order. In this work, we approach the problem from the lower-critical dimension of two. We construct a Fierz-complete renormalizable Lagrangian, compute the leading order β functions, fermion anomalous dimension, as well as the order parameter anomalous dimensions, and resolve the three universality classes corresponding to (i)–(iii). Before becoming equal to the Gaussian fixed point at N f = 1 , (ii) remains critical for all values of N f > N f , c ST ( N ) ≈ 0.56 + 1.48 N + O ( ε ) , which compares well with the estimate of previous studies. We further find that (iii) becomes equal to (i) when approaching N f = 1 . An instability is, however, only present in the susceptibility corresponding to the Gross-Neveu-Ising order parameter.
We study the Gross-Neveu-Yukawa field theory for the SO(N) symmetric traceless rank-2 tensor order parameter coupled to Majorana fermions using the $\ensuremath{\epsilon}$-expansion around upper critical dimensions of 3+1 to two loops. Previously we established in the one-loop calculation that the theory does not exhibit a critical fixed point for $N\ensuremath{\ge}4$, but that nevertheless the stable fixed point inevitably emerges at a large number Nf of fermion flavors. For ${N}_{f}<{N}_{f,c1}\ensuremath{\approx}N/2$, no critical fixed point exists; for Nf,c1<Nf<Nf,c2, a real critical fixed point emerges from the complex plane but fails to satisfy the additional stability conditions necessary for a continuous phase transition; and finally only for ${N}_{f}>{N}_{f,c2}\ensuremath{\approx}N$ does the fixed point satisfy the stability conditions as well. In the present work we compute the $O(\ensuremath{\epsilon})$ (two-loop) corrections to the critical flavor numbers Nf,c1 and Nf,c2. Most importantly, we observe a sharp decrease in Nf,c2 from its one-loop value, which brings it closer to the point Nf=1 relevant to the standard Gross-Neveu model. Some three-loop results are also presented and discussed.
We study quantum multicritical behavior in a $(2+1)\text{\ensuremath{-}}\mathrm{dimensional}$ Gross-Neveu-Yukawa field theory with eight-component Dirac fermions coupled to two triplets of order parameters that act as Dirac masses and transform as (1,0)+(0,1) representation under the $\mathrm{S}\mathrm{O}(4)\ensuremath{\simeq} \mathrm{S}\mathrm{O}(3)\ifmmode\times\else\texttimes\fi{}\mathrm{SO}(3)$ symmetry group. This field theory is relevant to spin-1/2 fermions on honeycomb or $\ensuremath{\pi}$-flux lattices, for example, near the transition points between an s-wave superconductor and a charge-density wave, on one side, and N\'eel order, on the other. Two triplets of such order parameters always allow for a common pair of two other order parameters that would complete them to the maximal set of compatible (anticommuting) orders of five. We first derive a unitary transformation in the Nambu (particle-hole) space which maps any two such triplets, possibly containing some superconducting orders, onto purely insulating order parameters. This allows one to consider a universal SO(4) Gross-Neveu-Yukawa description of the multicriticality without any Nambu doubling. We then proceed to derive the renormalization-group flow of the coupling constants at one-loop order in $4\ensuremath{-}\ensuremath{\epsilon}$ space-time dimensions, allowing also a more general set of order parameters transforming under $\mathrm{S}\mathrm{O}({n}_{a})\ifmmode\times\else\texttimes\fi{}\mathrm{S}\mathrm{O}({n}_{b})$. While for na=nb>2 in the bosonic sector and with fermions decoupled there is a stable fixed point of the flow, the Yukawa coupling to fermions quickly leads to its elimination by a generic fixed-point collision in the relevant range of fermion flavor numbers Nf. This suggests the replacement of the critical behavior by a runaway flow in the physical case na=nb=3. The structure of the renormalization group flow at ${n}_{a}\ensuremath{\ne}{n}_{b}$ is also discussed, and some nonperturbative arguments in favor of the stability of the decoupled critical point when na=3 and nb=1 in D=2+1 are provided.
We establish the universal torus low-energy spectra at the free Dirac fixed point and at the strongly coupled chiral Ising fixed point and their subtle crossover behavior in the Gross-Neuveu-Yukawa field theory with nD=4 component Dirac spinors in D=(2+1) dimensions. These fixed points and the field theories are directly relevant for the long-wavelength physics of certain interacting Dirac systems, such as repulsive spinless fermions on the honeycomb lattice or $\ensuremath{\pi}$-flux square lattice. The torus energy spectrum has been shown previously to serve as a characteristic fingerprint of relativistic fixed points and is a powerful tool to discriminate quantum critical behavior in numerical simulations. Here, we use a combination of exact diagonalization and quantum Monte Carlo simulations of strongly interacting fermionic lattice models, to compute the critical torus energy spectrum on finite-size clusters with periodic boundaries and extrapolate them to the thermodynamic limit. Additionally, we compute the torus energy spectrum analytically using the perturbative expansion in $\ensuremath{\epsilon}=4\ensuremath{-}D$, which is in good agreement with the numerical results, thereby validating the presence of the chiral Ising fixed point in the lattice models at hand. We show that the strong interaction between the spinor field and the scalar order-parameter field strongly influences the critical torus energy spectrum and we observe prominent multiplicity features related to an emergent symmetry predicted from the quantum field theory. Building on these results we are able to address the subtle crossover physics of the low-energy spectrum flowing from the chiral Ising fixed point to the Dirac fixed point, and analyze earlier flawed attempts to extract Fermi velocity renormalizations from the low-energy spectrum.
We study the universal critical properties of the QED3-Gross-Neveu-Yukawa model with N flavors of four-component Dirac fermions coupled to a real scalar order parameter at four-loop order in the $\ensuremath{\epsilon}$ expansion below four dimensions. For N=1, the model is conjectured to be infrared dual to the SU(2)-symmetric noncompact CP1 model, which describes the deconfined quantum critical point of the N\'eel-valence-bond-solid transition of spin-1/2 quantum antiferromagnets on the two-dimensional square lattice. For N=2, the model describes a quantum phase transition between an algebraic spin liquid and a chiral spin liquid in the spin-1/2 kagome antiferromagnet. For general N we determine the order-parameter anomalous dimension, the correlation length exponent, the stability critical exponent, as well as the scaling dimensions of SU(N) singlet and adjoint fermion mass bilinears at the critical point. We use Pad\'e approximants to obtain estimates of critical properties in 2+1 dimensions.