Pith. sign in

REVIEW 4 major objections 4 minor 1 cited by

Four-loop renormalisation-group equations are here derived for every renormalisable three-dimensional scalar-fermion theory, along with a new totally attractive infrared fixed point in an SO(3)×U(N) model, controlled by ε=1/N.

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 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 19:05 UTC pith:CLYMZFGD

load-bearing objection The four-loop template coefficients look like a genuine, carefully cross-checked computation; the part I would be most careful about is not the RGEs themselves but the SO(3)×U(N) fixed point built on them. the 4 major comments →

arxiv 2512.02123 v2 pith:CLYMZFGD submitted 2025-12-01 hep-th cond-mat.str-elhep-ph

General Four-Loop Beta Function for Scalar-Fermion Theories in Three Dimensions

York Schr\"oder , Emmanuel Stamou , Tom Steudtner , Max Uetrecht This is my paper
classification hep-th cond-mat.str-elhep-ph PACS 11.10.Hi
keywords renormalisation grouptemplate RGEsthree-dimensional field theoryscalar-fermion theoryYukawa couplingsextic scalar couplinglarge-N expansioninfrared fixed point
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.

This paper derives the most general renormalisation-group equations available for three-dimensional physics: the four-loop beta functions and anomalous dimensions (field-rescaling factors) of any renormalisable theory built from scalars and fermions, packaged as a set of universal 'template' coefficients that need only be contracted with a given theory's Yukawa and scalar-sextic coupling tensors. With those templates in hand, the authors scan for critical points and find a new one: an SO(3)×U(N) model whose couplings are drawn, at long distances, to a totally attractive infrared fixed point at a Yukawa coupling of order 1/N, putting the fixed point under perturbative control in a way that has been rare in three dimensions. The calculation is cross-checked by demanding that N=1 and N=2 supersymmetry emerge for special couplings, which yields hundreds of relations among the coefficients that all hold. If the results are correct, any three-dimensional scalar-fermion critical phenomenon — condensed-matter critical points, high-temperature effective theories — can be followed to four loops by algebraic contraction alone, and the SO(3)×U(N) model supplies a concrete interacting critical theory controlled by ε=1/N.

Core claim

The paper's central claim is that the four-loop template RGE coefficients quoted for the MS scheme — the scalar and fermion anomalous dimensions (10)–(13), the 132 Yukawa vertex coefficients (26) and the 150 sextic vertex coefficients (37) — are correct for the general Lagrangian (1), which contains every renormalisable three-dimensional interaction of real scalars and spin-1/2 fermions. The supporting evidence built into the paper is threefold: the coefficients satisfy 214 independent N=1 and 46 independent N=2 supersymmetry relations that must hold in any supersymmetry-preserving scheme; the four-loop coefficients involving traces over an odd number of gamma matrices all vanish, so the res

What carries the argument

Template RGEs are the load-bearing object: every beta function and anomalous dimension is written as a sum of index contractions of the Yukawa tensor Y and the totally symmetric sextic tensor η with universal numerical coefficients, so one four-loop computation serves every model in the class. Computationally, infrared rearrangement with a common mass parameter turns all needed Feynman integrals into massive tadpoles, evaluated in d=3−2ε at four loops — numerically to 65 digits, analytically wherever an integer-relation search succeeds — and reduced to master integrals by integration-by-parts. The fixed-point search is carried by a trace identity: the two-loop rate of tr(Y²) is a sum of mani

Load-bearing premise

The load-bearing premise is the large-N power counting of Sec. 5.2: the alternating scalar-fermion bubble chains (67) are asserted to factorise into subgraphs without logarithmic divergences and therefore to drop out of the beta functions, leaving the retained two- and four-loop terms, each carrying one power of ε=1/N, as the entire leading contribution; if any omitted or subleading diagram contributed at the same order in ε, the fixed-point position (74) and its 'totally IR

What would settle it

Independently recompute one load-bearing coefficient in a fully independent setup — for instance the four-loop Yukawa coefficient β^(4)_{Y,23} = −π²/8 or the ε⁴ coefficient of the master integral I_952, the multiple-zeta combination quoted in (88); a single mismatch shifts the fixed point (74). Sharper still: compute the six-loop correction to β_{αY} in the SO(3)×U(N) model, since the paper's perturbative-control claim requires the first omitted terms to be suppressed by ε relative to the retained four-loop O(ε) term — an O(1) six-loop coefficient would break it — or test the 214 N=1 supersymm

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

If this is right

  • Any renormalisable three-dimensional scalar-fermion theory now has its RG flow at four loops available as a purely algebraic contraction of the two published coupling tensors, so specific models no longer require their own multi-loop computations.
  • The extracted supersymmetry relations (214 for N=1, 46 for N=2) hold in any renormalisation scheme and provide a non-trivial consistency test for any future recalculation or extension of these coefficients.
  • The SO(3)×U(N) model realises an infrared fixed point with two positive stability eigenvalues, meaning all nearby RG flows are drawn into it; at the fixed point the scalar potential is bounded from below since α*_η > 0.
  • The underlying mechanism — spinor representations of SO(n) produce negative two-loop Yukawa beta functions — gives a constructive criterion for searching out further perturbative fixed points in three dimensions, and the paper notes other critical phenomena may yet be found.
  • Continuing the model toward four dimensions by dimensional continuation, the UV fixed point turns marginal already around δ≈0.0175 and shows no sign of surviving at δ=1, so the theory does not extend to an asymptotically safe one in d=4.

Where Pith is reading between the lines

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

  • My inference: projecting the template coefficients onto specific condensed-matter universality classes (chiral Ising, Gross–Neveu-type semimetal transitions) would yield four-loop critical exponents directly in d=3, an order of precision not currently available for the most general Yukawa–sextic theories.
  • My inference: the paper notes antisymmetric leg corrections are 'still absent at four loops' but expected eventually, by analogy with four dimensions where they start at three loops; if they first appear at six loops, fixed-point analyses of theories with hidden flavour symmetries will need scheme-dependent shifts at that order, and the totally attractive character of (74) could be modified.
  • My inference: because all odd-gamma-trace four-loop coefficients vanish, this calculation is insensitive to the naive-dimensional-regularisation ambiguity tied to the 3D Levi–Civita tensor; if odd traces contribute at six loops, N=1 supersymmetry in the MS scheme may fail at that order, mirroring the four-dimensional N=1/2 story.
  • My inference: the two stability eigenvalues differ by a factor 1/ε (ϑ1∝ε² versus ϑ2∝ε), so the flow approaches the fixed point along a very slow direction; this predicts a long crossover regime near the critical point that lattice or functional-RG studies of the SO(3)×U(N) model could look for.

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

4 major / 4 minor

Summary. This paper computes, in the MS scheme, template anomalous dimensions and beta functions to four loops for the most general renormalisable three-dimensional scalar–fermion theory with couplings Y and η. The main results are the scalar and fermion anomalous dimensions (10)–(13), 132 four-loop Yukawa vertex coefficients (26), and 150 four-loop sextic vertex coefficients (37), together with an appendix of high-precision four-loop master integrals. The authors then impose N=1 and N=2 supersymmetry and report that all resulting relations (214 and 46, respectively) are satisfied by their direct results. Finally, they specialise to an SO(3)×U(N) model, quote its two- and four-loop beta functions (72)–(73), and find a totally IR-attractive fixed point (74) with eigenvalues (75) in a large-N expansion with ε=1/N.

Significance. If the template coefficients are correct, this is a valuable model-independent result for three-dimensional scalar–fermion theories, analogous to four-dimensional template RGEs, and the fixed-point application is potentially interesting. The paper has genuine strengths: the N=1 and N=2 SUSY compatibility checks provide strong internal consistency, agreement with [79] is cited, and the master-integral data are presented to 65 digits with PSLQ-based analytic identifications where available. The main weakness is that the bridge from the general template coefficients to the quoted SO(3)×U(N) beta functions is not shown, and one displayed template formula contains an internal inconsistency. These issues are fixable but currently block verification of the headline IR fixed point.

major comments (4)
  1. [Sec. 5.2, Eqs. (72)–(73)] The model-specific beta functions βαY and βαη are quoted without derivation from the template coefficients of Sec. 3. In particular, the coefficients 4(6+π²), 48, −288, and −288π² are nontrivial contractions of the 132 Yukawa and 150 sextic coefficients, and the large-N ordering that justifies retaining only these terms is only sketched via (67)–(69). An error in any of these contractions would shift the fixed point (74) and the stability eigenvalues (75). Please provide the surviving contraction list, the large-N power counting for each retained term, and, ideally, an auxiliary file with the reduction from (26) and (37) to (72)–(73). This is the load-bearing step for the paper's main phenomenological claim.
  2. [Sec. 3.3, Eq. (32) and coefficient list (37)] In the displayed expression for β̂η,ηY²Y², the coefficient βη,24 is printed twice and βη,23 is missing, while the coefficient list gives distinct values βη,23=−1/4 and βη,24=−π²/32. As printed, two different tensor contractions are assigned the same coefficient. This is an internal inconsistency in one of the central template results and should be corrected.
  3. [Sec. 3.1, p. 4] The claim that antisymmetric leg corrections 'are still absent at four loops' is stated without proof or a reference. Since the template anomalous dimensions (9) and (12) are written with only symmetric external-index contractions, any antisymmetric contribution would constitute missing coefficients. Because the completeness of the four-loop template is a central claim, this assertion needs to be substantiated—either by an explicit diagrammatic argument or by a direct computation of the antisymmetric part.
  4. [Sec. 5.2, footnote 2 and Eqs. (67)–(69)] The large-N argument relies on the assertion that the alternating scalar–fermion bubble chains factorise into subgraphs without logarithmic divergences and therefore drop out, and that the remaining leading large-N contributions are exactly those quoted in (72)–(73). This is plausible but not demonstrated in the text. A more explicit derivation of the N-counting for all classes of diagrams contributing at the orders kept, and for those discarded, is needed to make the claim 'under perturbative control in a large-N limit' verifiable.
minor comments (4)
  1. [Eq. (29)] The decomposition of β̂η⁽⁴⁾ lists two identical labels 'η,η3'; the second term should evidently be 'η,η2Y2'.
  2. [Sec. 4] The 214 N=1 and 46 N=2 SUSY relations are said to be listed in separate files, but these are not included in the visible submission. Please state clearly in the text how to obtain them or include them as ancillary material.
  3. [Conclusion, Sec. 6] The phrase 'renormalisation-scheme-agnostic template expressions' overstates the explicit results: the numerical coefficients are scheme-dependent (MS), although the tensor structure and the SUSY relations may be scheme-independent. Suggest rewording for precision.
  4. [Eq. (32) and surrounding formulas] There are several index-notational typos in the displayed four-loop contractions, e.g., repeated indices in some of the tr(...) arguments. A careful proofread of Sec. 3.3 would improve usability.

Circularity Check

0 steps flagged

No circularity: fixed-point and SUSY results are independent outputs of direct four-loop computations.

full rationale

Walking the claimed chain — general Lagrangian (1) → template RGEs (10)-(13), (26), (37) → SUSY relations in Sec. 4 → model beta functions (72)-(73) → fixed point (74) — I find no step where an output is identical by construction to an input. The 282 template coefficients are direct MS-scheme loop results from four-loop massive tadpole masters, not fitted to the fixed point; the SUSY section explicitly reports that 214 (N=1) and 46 (N=2) relations "are all compatible with the direct results" / "are all fulfilled", i.e. SUSY is used as an external cross-check, not as a constraint to determine unknown coefficients. The fixed point α*_Y = ε/(6+π²) is the algebraic root of the quoted (72) at α_η = O(ε²), with coefficients 4ε and 4(6+π²) presented as contractions of the template results; there is no indication these coefficients were tuned to force the root, and the stability eigenvalues (75) are likewise outputs. The main caveats are verifiability gaps, not circularity: Sec. 5.2 quotes (72)-(73) without displaying the surviving tensor contractions among the 132/150 coefficients, the supporting software FoRGEr is described as "unpublished" (Sec. 3), and the general large-N assertion that bubble chains factorise without logarithmic divergences (footnote 2) leans on the self-cited [84]. However, the latter is supported in-text by the directly computed vanishing coefficients β(2)_Y,1 and β(4)_Y,10, so the self-citation is not the sole load-bearing evidence. These omitted-derivation/reproducibility limitations do not make the prediction equivalent to its inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 6 axioms · 0 invented entities

The central claims rest on standard perturbative-QFT machinery plus a handful of scheme/model assumptions: dimensional regularisation in 3−2ε with MS (standard), the completeness of the interaction basis (1) (the paper's defined setting), the γ-algebra handling of odd traces in 3D (delicate, Sec. 3.2), the large-N counting where bubble chains are asserted to be divergence-free (footnote 2, Sec. 5.2), and the expected-but-unproved absence of conformality-breaking antisymmetric leg corrections at four loops (Sec. 3.1). No data are fitted anywhere: the fixed-point position is a prediction of the computed coefficients, with ε=1/N as an expansion bookkeeping parameter. No new entities (particles, forces, dimensions) are postulated.

free parameters (1)
  • ε = 1/N (large-N expansion parameter) = 1/N, N free integer
    Controls perturbativity of the IR fixed point (α*_Y = ε/(6+π²)). It is an expansion bookkeeping parameter, not fitted to data; however the fixed point's existence and position scale with it, so the fixed-point claim is parameterised by this choice.
axioms (6)
  • domain assumption Dimensional regularisation to d=3−2ε with MS and infrared rearrangement via a common mass parameter separates UV from IR divergences for the massive tadpole computation.
    Invoked in Sec. 2 to justify computing only tadpole integrals; standard within the community but a substantive scheme and technique choice.
  • domain assumption The Lagrangian (1) exhausts the marginal interactions of 3D renormalisable scalar–fermion theories (Yukawa φ²ψ² and sextic φ⁶; no gauge, no four-fermi, no Chern–Simons terms).
    Definition of the paper's setting; gauge sectors are explicitly deferred to future work (Sec. 6).
  • domain assumption tr(γμγνγρ) = −2i ε^{μνρ}, with the conclusion that all odd-γ-trace contractions vanish at four loops (β^(4)_{Y,36–43}=0).
    Sec. 3.2 below Eq. (23): 'After integration of the fermion loop with an odd number of γ matrices... only the contractions... do not vanish'; the quoted vanishing of all such coefficients is a delicate γ-algebra/NDR-sensitive point on which the scheme-insensitivity claim rests.
  • domain assumption Alternating scalar–fermion bubble chains (67) have no logarithmic divergences and hence contribute zero to the beta functions in the large-N limit.
    Sec. 5.2, footnote 2: 'each diagram of (67) factorises into subgraphs without logarithmic divergences, thus not contributing to the β-functions... the vanishing is renormalisation-scheme independent.' Load-bearing for the 1/N counting that makes the fixed point perturbative.
  • domain assumption Absence of antisymmetric leg corrections to field-strength renormalisation at four loops, so computed fixed points are conformal.
    Sec. 3.1: 'We eventually expect the same phenomena to also occur in three dimensions, though they are still absent at four loops.' The paper itself marks this as an expectation, not a proven fact.
  • standard math PSLQ identification of rational/π² RGE coefficients from 65-digit numerical expansions is exact.
    App. A.1; standard integer-relation practice, but formally a numerical conjecture for the six constants c_s that remain unidentified.

pith-pipeline@v1.3.0-alltime-deepseek · 43391 in / 31367 out tokens · 316505 ms · 2026-08-03T19:05:04.476612+00:00 · methodology

0 comments
read the original abstract

We present general four-loop template $\beta$-functions and anomalous field dimensions for renormalisable scalar-fermion theories in three dimensions. By imposing $\mathcal{N}=1$ and $\mathcal{N}=2$ supersymmetry, we obtain relations between the template RGE coefficients, valid in any renormalisation scheme. Directly in $d=3$, we identify a new theory with a non-trivial IR fixed point that is under perturbative control in a large-$N$ limit. We provide up-to-date numerical results for all required massive tadpole master integrals up to four loops and complement them with analytic expressions where available.

Figures

Figures reproduced from arXiv: 2512.02123 by Emmanuel Stamou, Max Uetrecht, Tom Steudtner, York Schr\"oder.

Figure 1
Figure 1. Figure 1: RG flow for the theory (70) with ε = 10−1 , showing the Gaussian fixed point (α ∗ Y = α ∗ η = 0), the infrared fixed point (74), and the purely scalar ultraviolet fixed point (α ∗ Y = 0) [80–82,84]. Arrows point from the UV to the IR. The UV fixed point is connected to the Gaussian only through the trajectory αY = 0 (violet). The RG flow along the red trajectory is slower–around the IR fixed point, it is s… view at source ↗
Figure 2
Figure 2. Figure 2: The master integrals up to four loops, labelled by their [PITH_FULL_IMAGE:figures/full_fig_p021_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Alternative choices for dotted masters. Note that the integral [PITH_FULL_IMAGE:figures/full_fig_p031_3.png] view at source ↗

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Local CFTs extremise $F$

    hep-th 2026-04 unverdicted novelty 7.0

    Local CFTs lie at the extrema of the sphere free energy tilde F for nonlocal CFT lines, and maximize it when unitary.

Reference graph

Works this paper leans on

16 extracted references · 7 linked inside Pith · cited by 1 Pith paper

  1. [1]

    Renormalization group and the Kadanoff scaling picture,Phys. Rev. B4(1971) 3174–3183. [2]K. G. Wilson,Renormalization group and critical phenomena

    1971

  2. [2]

    Phase space cell analysis of critical behavior,Phys. Rev. B4(1971) 3184–3205. [3]K. G. Wilson and M. E. Fisher,Critical exponents in 3.99 dimensions,Phys. Rev. Lett.28 (1972) 240–243. [4]K. G. Wilson and J. Kogut,The renormalization group and theεexpansion,Phys. Rept. 12(1974) 75–199. [5]J. C. Le Guillou and J. Zinn-Justin,Critical Exponents from Field Th...

    1971

  3. [6]

    divergent graphs

    in three-dimensions at large N,Phys. Rev. Lett.48(1982) 574–576. [83]J. S. Hager,Six-loop renormalization group functions of O(n)-symmetricφ 6-theory and epsilon-expansions of tricritical exponents up toε 3,J. Phys. A35(2002) 2703–2711. [84]S. Kvedarait ˙e, T . Steudtner and M. Uetrecht,Revisiting theφ 6 theory in three dimensions at large N,Phys. Rev. D1...

    arXiv 1982

  4. [52]

    Salam and J

    [118]A. Salam and J. Strathdee,On Superfields and Fermi-Bose Symmetry,Phys. Rev. D11 (1975) 1521–1535. [119]M. T . Grisaru, W . Siegel and M. Rocek,Improved Methods for Supergraphs,Nucl. Phys. B 159(1979)

    1975

  5. [60]

    Zoller,Top-Yukawa effects on theβ-function of the strong coupling in the SM at four-loop level,JHEP02(2016) 095,[1508.03624]

    35 SciPost Physics Submission [58]M. Zoller,Top-Yukawa effects on theβ-function of the strong coupling in the SM at four-loop level,JHEP02(2016) 095,[1508.03624]. [59]K. Chetyrkin and M. Zoller,Leading QCD-induced four-loop contributions to the β-function of the Higgs self-coupling in the SM and vacuum stability,JHEP06(2016) 175,[1604.00853]. [60]I. Schie...

    Pith/arXiv arXiv 2016

  6. [147]

    Chetyrkin, S

    32 SciPost Physics Submission [7]K. Chetyrkin, S. Gorishny , S. Larin and F . Tkachov,Five-loop renormalization group calculations in the gφ 4 theory,Phys. Lett. B132(1983)

    1983

  7. [351]

    Kleinert, J

    [8]H. Kleinert, J. Neu, N. Schulte-Frohlinde, K. Chetyrkin and S. Larin,Five-loop renormalization group functions of O(n)-symmetricϕ 4-theory andε-expansions of critical exponents up toε 5,Phys. Lett. B272(1991) 39–44,[hep-th/9503230]. [9]P . Calabrese and P . Parruccini,Five-loopεexpansion for O(n)×O(m)spin models,Nucl. Phys. B679(2004) 568–596,[cond-mat...

    Pith/arXiv arXiv 1991

  8. [429]

    [120]S. R. Coleman and D. J. Gross,Price of asymptotic freedom,Phys. Rev. Lett.31(1973) 851–854. [121]A. D. Bond and D. F . Litim,Theorems for Asymptotic Safety of Gauge Theories,Eur. Phys. J. C77(2017) 429,[1608.00519]. [122]A. D. Bond and D. F . Litim,Price of Asymptotic Safety,Phys. Rev. Lett.122(2019) 211601,[1801.08527]. [123]G. Veneziano,Some Aspect...

    Pith/arXiv arXiv 1973

  9. [1992]

    [102]H. R. P . Ferguson, D. H. Bailey and S. Arno,Analysis of PSLQ, an Integer Relation Finding Algorithm,Mathematics of Computation68(1999) 351–369. [103]Y. Schroder and A. Vuorinen,High precision evaluation of four loop vacuum bubbles in three-dimensions,hep-ph/0311323. [104]S. P . Martin and M. T . Vaughn,Two-loop renormalization group equations for so...

    Pith/arXiv arXiv 1999

  10. [2020]

    Di Pietro and E

    [127]L. Di Pietro and E. Stamou,Scaling dimensions in QED 3 from theε-expansion,JHEP12 (2017) 054,[1708.03740]. [128]L. Di Pietro and E. Stamou,Operator mixing in theε-expansion: Scheme and evanescent-operator independence,Phys. Rev. D97(2018) 065007,[1708.03739]. [129]L. Di Pietro, D. Gaiotto, E. Lauria and J. Wu,3d Abelian Gauge Theories at the Boundary...

    Pith/arXiv arXiv 2017

  11. [2305]

    Appelquist, M

    [32]T . Appelquist, M. J. Bowick, D. Karabali and L. C. R. Wijewardhana,Spontaneous Breaking of Parity in (2+1)-dimensional QED,Phys. Rev. D33(1986)

    1986

  12. [2575]

    Nash,Higher Order Corrections in (2+1)-Dimensional QED,Phys

    [35]D. Nash,Higher Order Corrections in (2+1)-Dimensional QED,Phys. Rev. Lett.62 (1989)

    1989

  13. [3024]

    Rosenstein, B

    [36]B. Rosenstein, B. Warr and S. H. Park,Dynamical symmetry breaking in four Fermi interaction models,Phys. Rept.205(1991) 59–108. [37]O. Aharony and D. Fleischer,IR Dualities in General 3d Supersymmetric SU(N) QCD Theories,JHEP02(2015) 162,[1411.5475]. [38]O. Aharony ,Baryons, monopoles and dualities in Chern-Simons-matter theories,JHEP 02(2016) 093,[15...

    Pith/arXiv arXiv 1991

  14. [3704]

    Appelquist, D

    [34]T . Appelquist, D. Nash and L. C. R. Wijewardhana,Critical Behavior in (2+1)-Dimensional QED,Phys. Rev. Lett.60(1988)

    1988

  15. [3774]

    [33]T . W . Appelquist, M. J. Bowick, D. Karabali and L. C. R. Wijewardhana,Spontaneous Chiral Symmetry Breaking in Three-Dimensional QED,Phys. Rev. D33(1986)

    1986

  16. [3998]

    [94]K. G. Chetyrkin, M. Misiak and M. Munz,Beta functions and anomalous dimensions up to three loops,Nucl. Phys. B518(1998) 473–494,[hep-ph/9711266]. [95]J. Brod, L. Hüdepohl, E. Stamou and T . Steudtner,MaRTIn – Manual for the ”Massive Recursive Tensor Integration”,Comput. Phys. Commun.306(2025) 109372, [2401.04033]. 37 SciPost Physics Submission [96]J. ...

    Pith/arXiv arXiv 1998