Pith. sign in

REVIEW 4 minor 46 references

One correction exponent in the chiral Heisenberg model diverges as d o3 at order 1/N^{2}; a resummation that mixes π∂^{2}π with four-fermion operators restores finite leading-order values that match a direct 3d calculation.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-13 20:31 UTC pith:HCU57GNE

load-bearing objection Solid 1/N^{2} correction exponents for the chiral Heisenberg model, with a clean d=3 pole analysis and resummation that actually matches a direct 3D calculation.

arxiv 2603.21989 v2 pith:HCU57GNE submitted 2026-03-23 hep-th cond-mat.str-el

Correction exponents in the chiral Heisenberg model at 1/N²: singular contributions and operator mixing

classification hep-th cond-mat.str-el
keywords chiral Heisenberg modelcorrection exponents1/N expansionoperator mixingfour-fermion operatorsresummationGross-Neveu-Yukawacritical exponents
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper computes the two correction exponents ω± (the slopes of the beta functions at the infrared fixed point) of the chiral Heisenberg model to order 1/N^{2} in general dimension d. These quantities agree with the known four-dimensional ε-expansion when expanded near d=4. One of them, however, develops a simple pole as d approaches 3. The authors show that the pole is produced by operator mixing: the composite operator π∂^{2}π can mix with a system of four-fermion operators that become degenerate with it precisely at d=3. By assembling the finite elementary symmetric polynomials of the anomalous-dimension matrix and solving the resulting cubic characteristic equation, they obtain a resummed set of eigenvalues that remain finite at d=3 and already differ from the naïve large-N result at leading order. An independent calculation performed directly in the three-dimensional theory reproduces exactly the same resummed values, confirming that the procedure correctly captures the infrared spectrum.

Core claim

At order 1/N^{2} the correction exponent ω– associated with the operator π∂^{2}π diverges as d o3. After the 3 imes3 mixing matrix of π∂^{2}π with the four-fermion operators (q̄σ_a γ_μ q)^{2} and (q̄γ_μν q)^{2} is constructed, the elementary symmetric coefficients of its characteristic polynomial remain finite; the roots of that polynomial at d=3 are γ_0=0 and γ_±=±(8 η_1/n)√(2/3). These values agree with a direct large-N calculation performed in three dimensions and already modify the leading-order exponents.

What carries the argument

The cubic characteristic equation P(λ)=-λ^{3}+A(d)λ^{2}-B(d)λ+C(d)=0 whose coefficients A,B,C are the elementary symmetric polynomials of the anomalous-dimension matrix of {π∂^{2}π,(q̄σ_a γ_μ q)^{2},(q̄γ_μν q)^{2}}; because the poles cancel in A,B,C the equation supplies a resummation that remains valid at d=3.

Load-bearing premise

That the three elementary symmetric polynomials built from the anomalous-dimension matrix stay free of poles at d=3 through the orders needed for a leading-order resummation, even after possible mixing with higher four-fermion blocks.

What would settle it

Compute the same correction exponents to order 1/N^{2} by an independent method that works directly at d=3 (for example a conformal bootstrap or a lattice large-N simulation) and check whether the numerical values of the two non-zero eigenvalues equal ±(8 η_1/n)√(2/3).

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The leading-order correction exponents that enter the approach to criticality in three-dimensional graphene-related models are already modified by four-fermion mixing and must be replaced by the resummed values.
  • Any higher-dimensional operator whose canonical dimension crosses that of a multi-fermion operator at an integer d will require an analogous resummation before its large-N series can be trusted at that dimension.
  • The same characteristic-polynomial construction supplies a practical algorithm for extracting finite critical exponents from divergent 1/N series in other Gross-Neveu-Yukawa-type models.
  • Agreement between the resummed general-d result and a pure three-dimensional calculation gives a non-trivial consistency check that can be repeated for related universality classes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same mixing-induced poles are likely to appear in the ε-expansion of the four-dimensional theory once it is continued past d=3, offering a diagnostic for operator-spectrum rearrangements already visible in the Ising model.
  • Because the degeneracy is lifted only at order 1/N^{2}, the leading large-N spectrum of all operators of dimension ≥4 in three-dimensional fermion models should be re-examined systematically for hidden four-fermion partners.
  • The method extends immediately to the full tower of higher-spin four-fermion operators once their 1/N^{2} poles are computed, potentially reorganising the entire subleading spectrum.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 4 minor

Summary. The paper computes the correction exponents ω± of the chiral Heisenberg (Gross–Neveu–Yukawa) model at order 1/N^{2} in arbitrary dimension d. These exponents are identified with the critical dimensions of the operators π∂^{2}π and π^{4}. The new 1/N^{2} expressions, when expanded in ε near d=4, reproduce the known four-loop series of Zerf et al. The authors also obtain the anomalous dimensions of the operators (π^{2})^m and the traceless tensors π^{a1}…π^{aℓ} at the same order. A pole appears in γ_- as d o3; it is traced to operator mixing of π∂^{2}π with the four-fermion operators (q̄σ_a γ_μ q)^{2} and (q̄γ_{μ u} q)^{2}. A resummation based on the characteristic polynomial of the 3 imes3 anomalous-dimension matrix, whose elementary symmetric coefficients remain finite, yields the leading-order d=3 eigenvalues γ_0=0 and γ_±=±(8 η_1/n)√(2/3). These values are confirmed by an independent calculation performed directly in three dimensions.

Significance. The work supplies the first complete 1/N^{2} results for the correction exponents of a model of direct relevance to graphene criticality, and it demonstrates a general mechanism by which poles at integer dimensions arise from changes in the mixing pattern of composite operators. The resummation procedure is concrete, multiplies checked against both the ε-expansion and a direct d=3 computation, and modifies the leading-order spectrum already at O(1/N). The diagram-by-diagram tabulation and the ancillary file of all critical indices make the results immediately usable for future large-N or conformal-bootstrap studies of fermionic CFTs.

minor comments (4)
  1. [Sect. 3, Eq. (9)] In the introduction and in Eq. (9) the relation n=N tr 1_l is stated, but the precise value of the spinor trace for non-integer d is left implicit; a short clarifying sentence would help readers who work only with the d=3 or d=4 conventions.
  2. [Fig. 3] The caption of Fig. 3 mentions that tr 1 is approximated as 2 for 2<d<3 and 2d-4 for 3<d<4; this choice should be stated once in the main text as well, so that the figure can be read independently.
  3. A few typographical inconsistencies remain (e.g., “ind=3”, missing spaces around “d=3”, and the occasional use of “1/n” versus “1/N”). A light copy-edit would remove them.
  4. [end of Sect. 1] The ancillary file is mentioned but not described; a one-sentence statement of its contents (list of all indices through 1/N^{2}) would be useful for archival purposes.

Circularity Check

0 steps flagged

No significant circularity: 1/N^{2} results and the d=3 resummation are independently computed and cross-checked against external ε-expansion and a direct d=3 mixing matrix.

full rationale

The derivation chain is self-contained. Correction exponents ω± are obtained from explicit 1/n and 1/n^{2} Feynman diagrams (Figs. 1, 4–8; Appendices A–C) for the operators π∂^{2}π and π⁴; the resulting series (36)–(37) match the independent four-loop ε-expansion of Zerf et al. [4] near d=4. The pole in γ− as d→3 is isolated to the Φ3 term in (35) and traced to double-box subgraphs that diverge only at d=3. The proposed resummation rewrites the characteristic polynomial of the 3×3 mixing matrix of {π∂^{2}π, (q̄σaγμq)², (q̄γμνq)²} so that the elementary symmetric coefficients A,B,C remain finite (verified with the singular 1/n^{2} pieces (40),(51)); solving at d=3 recovers the eigenvalues of the independently computed leading-order d=3 mixing matrix (47)–(48). Self-citations to Gracey [6] and Manashov–Strohmaier [21] supply only lower-order indices and diagram technology; they are not load-bearing for the new 1/n^{2} expressions or the resummation, which are checked against external and internal benchmarks with no fitted parameters. No step reduces a claimed prediction to its own input by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

Pure large-N field theory; no free parameters are fitted. All background technology (critical equivalence of the CH and GNY models, non-multiplicative renormalization up to 1/n^{2}, conformal bootstrap for dressed propagators) is taken from the established literature and used without modification. No new particles or forces are postulated.

axioms (3)
  • domain assumption Critical equivalence of the d-dimensional chiral Heisenberg model to the non-linear sigma-model-like action (8) for 2<d<4
    Invoked in Sect. 3; standard for Gross–Neveu–Yukawa models but not re-derived here.
  • domain assumption Anomalous dimensions extracted from the simple-pole residues of the Z-factors via γ=2u∂_u Z_1|u=1 up to O(1/n^{2})
    Eq. (18); justified by Derkachov–Manashov 1998 and used throughout.
  • ad hoc to paper Only the B_0 block of four-fermion operators produces poles at d=3 at O(1/n^{2}); higher blocks B_k (k≥2) remain finite
    Stated in Sect. 5 after examining the superficial degree of divergence of representative diagrams; not proven for all diagrams.

pith-pipeline@v1.1.0-grok45 · 23063 in / 2421 out tokens · 27215 ms · 2026-07-13T20:31:35.783624+00:00 · methodology

0 comments
read the original abstract

We calculate the correction exponents in the chiral Heisenberg model in the $1/N$ expansion. These exponents are related to the slopes of $\beta$ functions at the phase transition point. We present the results at order $1/N^2$ and check that they agree with the results of the $\epsilon$ expansion near $d = 4$. We find that one of the correction exponents diverges as $d \to 3$. We argue that the appearance of the pole is a rather general phenomenon and is associated with operator mixing involving the system of four-fermion operators. After analyzing the operator mixing structure, we propose a resummation procedure which modifies the exponents already at leading order. We also perform calculations directly in the three-dimensional model and find complete agreement with the resummed exponents.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

46 extracted references · 16 canonical work pages · 6 internal anchors

  1. [1]

    Kadanoff, Scaling laws for Ising models near T(c)

    L.P. Kadanoff, Scaling laws for Ising models near T(c). Physics Physique Fizika2, 263– 272 (1966). https://doi.org/10.1103/ PhysicsPhysiqueFizika.2.263

  2. [2]

    Wilson, J.B

    K.G. Wilson, J.B. Kogut, The Renormalization group and the epsilon expansion. Phys. Rept. 12, 75–199 (1974). https://doi.org/10.1016/ 0370-1573(74)90023-4

  3. [3]

    Mihaila, N

    L.N. Mihaila, N. Zerf, B. Ihrig, I.F. Herbut, M.M. Scherer, Gross-Neveu-Yukawa model at three loops and Ising critical behavior of Dirac systems. Phys. Rev. B96(16), 165133 (2017). https://doi.org/ 10.1103/PhysRevB.96.165133. arXiv:1703.08801 [cond-mat.str-el]

  4. [5]

    Gracey, Critical exponent ω in the Gross-Neveu- Yukawa model atO(1/N)

    J.A. Gracey, Critical exponent ω in the Gross-Neveu- Yukawa model atO(1/N). Phys. Rev. D96(6), 065015 (2017). https://doi.org/10.1103/PhysRevD.96. 065015. arXiv:1707.05275 [hep-th]

  5. [6]

    Gracey, Large N critical exponents for the chiral Heisenberg Gross-Neveu universality class

    J.A. Gracey, Large N critical exponents for the chiral Heisenberg Gross-Neveu universality class. Phys. Rev. D 97(10), 105009 (2018). https://doi.org/10.1103/ PhysRevD.97.105009. arXiv:1801.01320 [hep-th]

  6. [7]

    Gracey, A

    J.A. Gracey, A. Maier, P. Marquard, Y . Schr¨oder, Anoma- lous dimensions and critical exponents for the Gross- Neveu-Yukawa model at five loops. Phys. Rev. D 112(8), 085029 (2025). https://doi.org/10.1103/ lmpx-n3mj. arXiv:2507.22594 [hep-th]

  7. [8]

    Gross, A

    D.J. Gross, A. Neveu, Dynamical Symmetry Break- ing in Asymptotically Free Field Theories. Phys. Rev. D10, 3235 (1974). https://doi.org/10.1103/ PhysRevD.10.3235

  8. [9]

    Zamolodchikov, A.B

    A.B. Zamolodchikov, A.B. Zamolodchikov, Factorized s Matrices in Two-Dimensions as the Exact Solutions of Certain Relativistic Quantum Field Models. Annals Phys.120, 253–291 (1979). https://doi.org/10. 1016/0003-4916(79)90391-9

  9. [10]

    ZinnJustin, Four fermion interaction near four- dimensions

    J. ZinnJustin, Four fermion interaction near four- dimensions. Nucl. Phys. B367, 105–122 (1991). https: //doi.org/10.1016/0550-3213(91)90043-W

  10. [11]

    Gracey, Three loop calculations in the O(N) Gross- Neveu model

    J.A. Gracey, Three loop calculations in the O(N) Gross- Neveu model. Nucl. Phys. B341, 403–418 (1990). https://doi.org/10.1016/0550-3213(90) 90186-H

  11. [12]

    Gracey, Computation of the three loop Beta function of the O(N) Gross-Neveu model in minimal subtraction

    J.A. Gracey, Computation of the three loop Beta function of the O(N) Gross-Neveu model in minimal subtraction. Nucl. Phys. B367, 657–674 (1991). https://doi. org/10.1016/0550-3213(91)90012-M

  12. [13]

    Karkkainen, R

    L. Karkkainen, R. Lacaze, P. Lacock, B. Peters- son, Critical behavior of the three-dimensional Gross- Neveu and Higgs-Yukawa models. Nucl. Phys. B 415, 781–796 (1994). https://doi.org/10.1016/ 0550-3213(94)90309-3. [Erratum: Nucl.Phys.B 438, 650–650 (1995)]. arXiv:hep-lat/9310020

  13. [14]

    Gross-Neveu-Heisenberg criticality from $2+\boldsymbol{\epsilon}$ expansion

    K. Ladovrechis, S. Ray, T. Meng, L. Janssen, Gross- Neveu-Heisenberg criticality from 2+ε expansion. Phys. Rev. B107(3), 035151 (2023). https://doi.org/ 10.1103/PhysRevB.107.035151. arXiv:2209.02734 [cond-mat.str-el]

  14. [15]

    Huang, Q

    R. Huang, Q. Jin, Y . Li, On the seven-loop renor- malization of Gross-Neveu model. JHEP06, 134 (2025). https://doi.org/10.1007/JHEP06(2025)

  15. [16]

    arXiv:2504.00713 [hep-th]

  16. [17]

    Gracey, Calculation of exponent eta to O(1/N**2) in the O(N) Gross-Neveu model

    J.A. Gracey, Calculation of exponent eta to O(1/N**2) in the O(N) Gross-Neveu model. Int. J. Mod. Phys. A6, 395–408 (1991). https://doi.org/10.1142/ S0217751X91000241. [Erratum: Int.J.Mod.Phys.A 6, 2755 (1991)]

  17. [18]

    Gracey, Anomalous mass dimension at O(1/N**2) in the O(N) Gross-Neveu model

    J.A. Gracey, Anomalous mass dimension at O(1/N**2) in the O(N) Gross-Neveu model. Phys. Lett. B 297, 293–297 (1992). https://doi.org/10.1016/ 0370-2693(92)91265-B

  18. [19]

    Vasiliev, S.E

    A.N. Vasiliev, S.E. Derkachov, N.A. Kivel, A.S. Stepa- nenko, The 1/n expansion in the Gross-Neveu model: Conformal bootstrap calculation of the index eta in or- der 1/n**3. Theor. Math. Phys.94, 127–136 (1993). https://doi.org/10.1007/BF01019324

  19. [20]

    Gracey, Computation of critical exponent eta at O(1/N**3) in the four Fermi model in arbitrary dimen- sions

    J.A. Gracey, Computation of critical exponent eta at O(1/N**3) in the four Fermi model in arbitrary dimen- sions. Int. J. Mod. Phys. A9, 727–744 (1994). https:// doi.org/10.1142/S0217751X94000340. arXiv:hep- th/9306107

  20. [21]

    Vasiliev, A.S

    A.N. Vasiliev, A.S. Stepanenko, The 1/n expansion in the Gross-Neveu model: Conformal bootstrap calculation of the exponent 1/nu to the order 1/n**2. Theor. Math. Phys.97, 1349–1354 (1993). https://doi.org/10. 1007/BF01015764

  21. [22]

    Manashov, M

    A.N. Manashov, M. Strohmaier, Correction exponents in the Gross–Neveu–Yukawa model at1/N2. Eur. Phys. J. C78(6), 454 (2018). https://doi.org/10.1140/ epjc/s10052-018-5902-1 . arXiv:1711.02493 [hep- th]

  22. [24]

    Herbut, Interactions and phase transitions on graphene’s honeycomb lattice

    I.F. Herbut, Interactions and phase transitions on graphene’s honeycomb lattice. Phys. Rev. Lett. 97, 146401 (2006). https://doi.org/10.1103/ 13 PhysRevLett.97.146401. arXiv:cond-mat/0606195

  23. [26]

    Theory of interacting electrons on the honeycomb lattice

    I.F. Herbut, V . Juricic, B. Roy, Theory of interacting elec- trons on the honeycomb lattice. Phys. Rev. B79, 085116 (2009). https://doi.org/10.1103/PhysRevB.79. 085116. arXiv:0811.0610 [cond-mat.str-el]

  24. [27]

    Raghu, X.L

    S. Raghu, X.L. Qi, C. Honerkamp, S.C. Zhang, Topo- logical mott insulators. Physical Review Letters100(15) (2008). https://doi.org/10.1103/physrevlett. 100.156401. URL http://dx.doi.org/10.1103/ PhysRevLett.100.156401

  25. [28]

    Weyl semimetals with short-range interactions

    J. Maciejko, R. Nandkishore, Weyl semimetals with short-range interactions. Phys. Rev. B90(3), 035126 (2014). https://doi.org/10.1103/PhysRevB.90. 035126. arXiv:1311.7133 [cond-mat.str-el]

  26. [29]

    Kobayashi, T

    K. Kobayashi, T. Ohtsuki, K.I. Imura, I.F. Herbut, Den- sity of states scaling at the semimetal to metal transition in three dimensional topological insulators. Physical Review Letters112(1) (2014). https://doi.org/10. 1103/physrevlett.112.016402. URL http://dx. doi.org/10.1103/PhysRevLett.112.016402

  27. [30]

    New Quantum Transition in Weyl Semimetals with Correlated Disorder

    T. Louvet, D. Carpentier, A.A. Fedorenko, New quantum transition in Weyl semimetals with cor- related disorder. Phys. Rev. B95(1), 014204 (2017). https://doi.org/10.1103/PhysRevB.95. 014204. arXiv:1609.08368 [cond-mat.mes-hall]

  28. [31]

    Critical exponents at the unconventional disorder-driven transition in a Weyl semimetal

    S.V . Syzranov, P.M. Ostrovsky, V . Gurarie, L. Radz- ihovsky, Critical exponents at the unconventional disorder-driven transition in a Weyl semimetal. Phys. Rev. B93(15), 155113 (2016). https://doi.org/ 10.1103/PhysRevB.93.155113. arXiv:1512.04130 [cond-mat.mes-hall]

  29. [32]

    Vasilev,The field theoretic renormalization group in critical behavior theory and stochastic dynamics (2004)

    A.N. Vasilev,The field theoretic renormalization group in critical behavior theory and stochastic dynamics (2004)

  30. [33]

    Vasiliev, M.Y

    A.N. Vasiliev, M.Y . Nalimov, Analog of Dimensional Regularization for Calculation of the Renormalization Group Functions in the 1/n Expansion for Arbitrary Di- mension of Space. Theor. Math. Phys.55, 423–431 (1983).https://doi.org/10.1007/BF01015800

  31. [34]

    Derkachov, A.N

    S.E. Derkachov, A.N. Manashov, The Simple scheme for the calculation of the anomalous dimensions of com- posite operators in the 1/N expansion. Nucl. Phys. B 522, 301–320 (1998). https://doi.org/10.1016/ S0550-3213(98)00103-5. arXiv:hep-th/9710015

  32. [35]

    Korchemsky, On level crossing in conformal field theories

    G.P. Korchemsky, On level crossing in conformal field theories. JHEP03, 212 (2016). https://doi.org/10. 1007/JHEP03(2016)212. arXiv:1512.05362 [hep-th]

  33. [36]

    Derkachov, A.N

    S.E. Derkachov, A.N. Manashov, Critical dimensions of composite operators in the nonlinear sigma model. Theor. Math. Phys.116, 1034–1049 (1998). https: //doi.org/10.1007/BF02557145

  34. [37]

    Giombi, V

    S. Giombi, V . Kirilin, E. Skvortsov, Notes on Spin- ning Operators in Fermionic CFT. JHEP05, 041 (2017). https://doi.org/10.1007/JHEP05(2017)

  35. [38]

    arXiv:1701.06997 [hep-th]

  36. [39]

    Y . Ji, A.N. Manashov, Operator mixing in fermionic CFTs in noninteger dimensions. Phys. Rev. D 98(10), 105001 (2018). https://doi.org/10.1103/ PhysRevD.98.105001. arXiv:1809.00021 [hep-th]

  37. [40]

    Speer,Generalized Feynman Amplitudes

    E.R. Speer,Generalized Feynman Amplitudes. (AM-62), Volume 62(Princeton University Press, 1969)

  38. [41]

    Caron-Huot, M

    S. Caron-Huot, M. Kologlu, P. Kravchuk, D. Meltzer, D. Simmons-Duffin, Detectors in weakly-coupled field theories. JHEP04, 014 (2023). https://doi.org/10. 1007/JHEP04(2023)014. arXiv:2209.00008 [hep-th]

  39. [42]

    Double-logs, Gribov-Lipatov reciprocity and wrapping

    V .N. Velizhanin, Double-logs, Gribov-Lipatov reciprocity and wrapping. JHEP08, 092 (2011). https://doi.org/10.1007/JHEP08(2011)092. arXiv:1104.4100 [hep-th]

  40. [43]

    Velizhanin, Generalised double-logarithmic equa- tion in QCD

    V .N. Velizhanin, Generalised double-logarithmic equa- tion in QCD. Mod. Phys. Lett. A32(38), 1750213 (2017). https://doi.org/10.1142/ S0217732317502133. arXiv:1412.7143 [hep-ph]

  41. [44]

    Velizhanin, Exact result in N = 4 SYM theory: generalised double-logarithmic equation

    V .N. Velizhanin, Exact result in N = 4 SYM theory: generalised double-logarithmic equation. JHEP05, 176 (2022). https://doi.org/10.1007/JHEP05(2022)

  42. [45]

    arXiv:2201.04616 [hep-th]

  43. [46]

    Manashov, S

    A.N. Manashov, S. Moch, L.A. Shumilov, Anoma- lous dimensions at small spins. JHEP09, 106 (2025). https://doi.org/10.1007/JHEP09(2025)

  44. [47]

    arXiv:2506.05132 [hep-th]

  45. [48]

    Henriksson, S.R

    J. Henriksson, S.R. Kousvos, M. Reehorst, Spec- trum continuity and level repulsion: the Ising CFT from infinitesimal to finite ε. JHEP02, 218 (2023). https://doi.org/10.1007/JHEP02(2023)

  46. [49]

    arXiv:2207.10118 [hep-th]