REVIEW 2 major objections 4 minor 44 references
In the (2+1)-dimensional Gross-Neveu-Yukawa model, a beyond-mean-field renormalization group computation finds de Haas–van Alphen oscillations in the chiral phase boundary and a tricritical point that moves with magnetic field.
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 · deepseek-v4-flash
2026-08-04 10:03 UTC pith:WCOOSU6C
load-bearing objection First FRG computation of the GNY model with magnetic field, but the sharp-regulator floor could be inflating the dHvA features — a smooth-regulator check is needed before the central claim is fully solid. the 2 major comments →
The (2+1)-dimensional Gross-Neveu-Yukawa model at finite temperature, density, and magnetic field within the Functional Renormalization Group
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the full phase diagram of the (2+1)-dimensional GNY model, computed beyond mean field via the FRG in the local potential approximation, contains: (i) magnetic catalysis at large field and small chemical potential; (ii) inverse magnetic catalysis at intermediate chemical potential and small field; (iii) de Haas–van Alphen oscillations in the µ–|qB| plane at small field and large µ, realized as a sequence of first-order transitions into the broken phase followed by second-order restoring transitions; and (iv) a tricritical point that moves to higher T and lower µ with increasing |qB|. The bosonization of the Gross-Neveu model with a dynamical scalar and a Yukawa coupl
What carries the argument
The Landau-level truncated fermion loop in the flow equation for the effective potential, expressed by the floor function N_LL(k^2) = floor(k^2/(2|qB|)), which counts how many Landau levels contribute at a given renormalization-group scale k. This replacement converts the momentum integral into a sum over a finite number of levels that changes discretely as k flows, producing the oscillatory features. The boson loop is unaffected by the magnetic field, and the sharp momentum-space regulator makes the Landau-level sum finite.
Load-bearing premise
The sharp momentum-space regulator, which makes the number of Landau levels entering the flow a step function of k^2, is what produces the sharp oscillations; if a smooth regulator changes or washes out those transitions, the predicted de Haas–van Alphen structure could be an artifact of the regulator rather than a property of the model.
What would settle it
A direct check is to solve the same FRG flow equation with a smooth regulator shape function (e.g., an exponential or power-law cutoff) and see whether the de Haas–van Alphen oscillations in the µ–|qB| phase diagram and the associated first-order transitions survive. Alternatively, a lattice simulation of the (2+1)-dimensional GNY model at finite density and magnetic field could test the existence and location of the first-order transitions and the tricritical-point shift.
If this is right
- Magnetic catalysis: at zero chemical potential the critical temperature grows monotonically with |qB|.
- Inverse magnetic catalysis: at intermediate µ and small |qB|, the chemical potential at the second-order transition decreases before magnetic catalysis takes over at larger |qB|.
- De Haas–van Alphen oscillations: the chiral condensate shows repeated first-order transitions into the broken phase and second-order restoring transitions when |qB| is varied at fixed µ in the large-µ small-|qB| region.
- A tricritical point in the T–µ plane shifts to higher T and lower µ as |qB| increases, and a back-bending of the phase boundary appears at intermediate field strengths.
- These phenomena are already present in the Gross-Neveu model; the GNY model reproduces them beyond mean field with higher numerical precision due to the improved finite-volume numerical solver.
Where Pith is reading between the lines
- If the oscillatory features are regulator-induced (as the appendix acknowledges that a smooth regulator shape would make the Landau-level transitions more gradual), the de Haas–van Alphen oscillations may not be a genuine property of the model; a test with a smooth regulator would settle this.
- The finding that the lowest-Landau-level approximation becomes exact in the infrared (for k below the square root of 2|qB|) suggests LLL dominance is robust, but the first-order transitions at larger |qB| could depend on the sharp cutoff.
- The finite-volume numerical scheme used here could be applied to the (3+1)-dimensional quark-meson model with pions to resolve the deep-infrared region and potentially improve earlier functional renormalization group calculations of the QCD-type phase diagram.
- Since the chiral condensate is the order parameter, the predicted sequence of first-order transitions might be observable in condensed-matter analogues such as graphene in a magnetic field, if the GNY model captures the relevant universality class.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the (2+1)-dimensional Gross-Neveu-Yukawa model at finite temperature, chemical potential, and magnetic field using the Functional Renormalization Group in the local potential approximation. It derives the FRG flow equation for the effective potential, first without and then with a magnetic field, and solves it with a finite-volume Kurganov-Tadmor scheme, fixing the Yukawa coupling g^2 by requiring the vacuum condensate to sit at σ0=1. Results are presented as phase diagrams in the (T, μ), (T, |qB|), and (μ, |qB|) planes for N_f=2 and Λ=10^3. The main reported findings are magnetic catalysis at μ=0, inverse magnetic catalysis at intermediate μ and small |qB|, de Haas–van Alphen oscillations with sequences of first-order transitions, a TCP that shifts to higher T and lower μ with increasing |qB|, and a back-bending of the phase boundary.
Significance. If the reported phase structure is correct, this is a useful beyond-mean-field FRG confirmation that the GNY model reproduces features previously found in the Gross-Neveu model via mean-field, OPT, and lattice methods. The derivations in Appendix A are explicit, the numerical parameters are stated, and the coupling is fixed by a standard renormalization condition rather than by fitting the phase diagram, so there is no circularity. However, the central qualitative claims are not yet supported by an adequate regulator-sensitivity analysis: the non-analytic floor function in the Landau-level count enters directly into the fermionic loop, and the paper itself concedes that a smooth regulator would change the transition. For this reason the results are significant but currently conditional.
major comments (2)
- [§II.B, Eqs. (17)–(18), and footnote 2] The dHvA oscillations and the associated sequences of first-order transitions in Fig. 2(a)–(c) occur in the small-|qB|, large-μ region where N_LL(k^2)=floor(k^2/(2|qB|)) discontinuously truncates the Landau-level sum. Because this floor function enters directly into the fermionic loop (Eq. (18)), the non-analytic k-dependence of the flow is inherited from the sharp Litim regulator. Footnote 2 acknowledges that a smooth regulator would make the transition 'more gradual,' but no smooth-regulator run is provided. Since the conclusion explicitly claims these features 'also exist in the Gross-Neveu-Yukawa model,' the paper needs to show that the oscillations and first-order segments survive, with comparable amplitude and period, for a smooth regulator shape function (e.g., exponential or power-law). If they weaken or disappear, the phase-diagram claims must be rephrased as regulator-dependent
- [§III and Appendix B] The numerical setup is described in detail (Λ=10^3, k_IR=10^-2, Δσ=5×10^-3), but no convergence or truncation-error tests are reported. In particular, there is no variation of Δσ, k_IR, or Λ, and no test that the LPA truncation is adequate for the quantitatively claimed phase boundaries, such as the TCP location or the amplitude of the dHvA oscillations. The statement in §IV that the hydrodynamic scheme resolves features 'with unprecedented precision' is not supported by any quantitative comparison. Please add a convergence study and either provide numerical error estimates or soften the precision claim.
minor comments (4)
- [Abstract and Introduction] Language issues: 'allows for go further' and similar phrasing should be corrected. Also, 'allows to resolve' appears in §IV and should be 'allows us to resolve'.
- [Fig. 2(b)] The four curves are distinguished only by color in the text description. For print/accessibility, add explicit labels or markers for μ=0.92, 0.95, 1.01, 1.05.
- [§III, Fig. 3] The claimed TCP shift would be easier to verify if the TCP coordinates were tabulated or if the three phase boundaries were overlaid in one panel. As written, the shift is only apparent by visual comparison of separate panels.
- [Eq. (10)] The displayed equation '∂tU(t,σ) = −2' is incomplete in the text; presumably a diagrammatic term is missing in rendering. This should be fixed in the final version.
Circularity Check
No significant circularity: the phase diagram is computed from the FRG flow with only a standard scale-setting condition for g^2.
full rationale
The derivation chain is self-contained. The only fitted parameter is g^2, tuned via Newton-Raphson so that the IR effective potential has its minimum at sigma_0 = 1 (Appendix B, Eq. B4). This is a standard renormalization/scale-setting condition, not a fit to any target phase-diagram feature: the dHvA oscillations, TCP shift, and back-bending do not enter as inputs or constraints. All results are reported in units of h sigma_0, so this scale choice does not inject the predicted structure. The dHvA oscillations and first-order transitions are computed from the Landau-level sum in Eq. (18), with N_LL(k^2) = floor(k^2/(2|qB|)) in Eq. (17); the discrete Landau-level spectrum is part of the model input, and the floor function follows from the sharp Litim regulator rather than from a fitted or self-cited premise. The paper explicitly notes in footnote 2 that a smooth regulator would make the transition 'more gradual', which flags a genuine regulator-robustness concern for the sharpness of the oscillations, but that is a correctness/validity risk, not circularity: the output is not equivalent to the input by construction. The comparisons with mean-field, OPT, and lattice results in Refs. [20-22] are external benchmarks from independent groups, and the hydrodynamic solver is cited to Refs. [27,28], neither of which is a load-bearing self-citation by the authors. The LPA truncation and Litim regulator are stated approximations, not smuggled-in ans"atze. Overall, no step reduces to its own input, and no central claim is forced by a self-citation chain.
Axiom & Free-Parameter Ledger
free parameters (1)
- g^2 (four-fermion/Yukawa UV coupling) =
not quoted; tuned to give σ0=1
axioms (5)
- standard math The Wetterich flow equation (Eq. 3) is exact; LPA is the only truncation.
- domain assumption The Litim regulator shape functions (Eqs. 7-8) produce physical results; sharp-cutoff discontinuities do not generate spurious phase structure.
- domain assumption In a 2+1-dimensional model with a reducible 4-component Clifford representation, N_f=2 and dγ=4 define the physical flavor content.
- domain assumption Fermions move only in the transverse plane (p_z=0) with all flavors carrying the same charge q.
- domain assumption Adding a scalar kinetic term by hand in the Hubbard-Stratonovich bosonization yields the GNY model whose phase structure is representative of QCD-like chiral symmetry breaking.
read the original abstract
We investigate the phase diagram of the (2+1)-dimensional Gross-Neveu-Yukawa (GNY) model at finite temperature, density, and magnetic field beyond mean-field, using the Functional Renormalization Group (FRG) in the local potential approximation. Large magnetic fields result in magnetic catalysis, a dimensional reduction of the system, and enhancement of chiral symmetry breaking. We employ a hydrodynamical algorithm to solve the FRG flow equation for the effective potential, which allows for go further into the infrared region than with previously used methods. We find that the chiral condensate exhibits non-trivial behavior in various regions of the phase diagram: several first-order phase transitions and de Haas -- van Alphen oscillations at small magnetic field and large chemical potential, as well as a critical endpoint which shifts to higher temperature with increasing magnetic field.
Figures
Reference graph
Works this paper leans on
-
[1]
For vanishing magnetic field, Fourier-transforming the right-hand side of Eq
Zero Magnetic Field In LPA, only the effective potentialUdepends on the RG scale, such that the Wetterich equation (3) reads V2∂tU(t, σ) =∂tΓt[ ¯ψ, ψ, φ] φ=σ, ¯ψ=ψ=0 =STr 1 2 ∂tRt Γ(2) t [Φ] +R t −1 φ=σ, ¯ψ=0,ψ=0 ,(A1) whereV 2 =β(2π) 2δ(2)(0)is the two-dimensional volume element. For vanishing magnetic field, Fourier-transforming the right-hand side of E...
-
[2]
Therefore, for nonvanishing magnetic field, only the cal- culation of the fermion loop changes
Nonzero Magnetic Field In this model, only fermions carry electric charge. Therefore, for nonvanishing magnetic field, only the cal- culation of the fermion loop changes. To introduce the magnetic field we employ the replacements from Eqs. (15), (16), −2 =d γ Nf V2 |qB| 4π ∞X n=−∞ ∞X l=0 αl k2θ k2 −2|qB|l (νn +iµ) 2 +k 2 +m 2 f (t, σ) =d γ Nf βV2 |qB| 8π ...
-
[3]
This allows to use methods from numerical hydrodynamics to solve the flow equations
FRG Flow Equation in Conservative F orm Recently it has been found that the flow equation for the effective potential can be cast into conservative form [27]. This allows to use methods from numerical hydrodynamics to solve the flow equations. Here, we employ the Kurganov-Tadmor (KT) scheme [32], a high- resolution finite-volume method for solving hyperbo...
-
[4]
We fit the parameterg 2 using a Newton- Raphson root-finding algorithm so that the resulting ef- fective potential in the IR has a minimum atσ 0 = 1
Initial and Boundary Conditions The initial condition (UV limit) for the effective po- tential is given by its classical value in vacuum, UΛ(σ) = (hσ)2 2g2 .(B4) We seth= 1, and measure all dimensionful quantities in units ofσ 0. We fit the parameterg 2 using a Newton- Raphson root-finding algorithm so that the resulting ef- fective potential in the IR ha...
-
[5]
R. C. Duncan and C. Thompson, Formation of Very Strongly Magnetized Neutron Stars: Implications for Gamma-Ray Bursts, The Astrophysical Journal392, L9 (1992)
1992
-
[6]
D. Grasso and H. R. Rubinstein, Magnetic fields in the early universe, Phys. Rept.348, 163 (2001), arXiv:astro- ph/0009061
arXiv 2001
-
[7]
V. Skokov, A. Y. Illarionov, and V. Toneev, Estimate of the magnetic field strength in heavy-ion collisions, Int. J. Mod. Phys. A24, 5925 (2009), arXiv:0907.1396 [nucl-th]
Pith/arXiv arXiv 2009
-
[8]
K. Tuchin, Time and space dependence of the electro- magnetic field in relativistic heavy-ion collisions, Phys. Rev. C88, 024911 (2013), arXiv:1305.5806 [hep-ph]
Pith/arXiv arXiv 2013
-
[9]
J. O. Andersen, W. R. Naylor, and A. Tranberg, Phase diagram of QCD in a magnetic field: A review, Rev. Mod. Phys.88, 025001 (2016), arXiv:1411.7176 [hep-ph]
Pith/arXiv arXiv 2016
-
[10]
G. Endrodi, Critical point in the QCD phase diagram for extremely strong background magnetic fields, JHEP07, 173, arXiv:1504.08280 [hep-lat]
-
[11]
V. A. Miransky and I. A. Shovkovy, Quantum field the- ory in a magnetic field: From quantum chromodynamics to graphene and Dirac semimetals, Phys. Rept.576, 1 (2015), arXiv:1503.00732 [hep-ph]
Pith/arXiv arXiv 2015
-
[12]
M. Ferreira, P. Costa, O. Louren¸ co, T. Frederico, and C. Providˆ encia, Inverse magnetic catalysis in the (2+1)-flavor Nambu-Jona-Lasinio and Polyakov-Nambu- Jona-Lasinio models, Phys. Rev. D89, 116011 (2014), arXiv:1404.5577 [hep-ph]
Pith/arXiv arXiv 2014
-
[13]
R. Wen, S. Yin, W.-j. Fu, and M. Huang, Functional renormalization group study of neutral and charged pion under magnetic fields in the quark-meson model, Physical Review D108, 076020 (2023), arXiv:2306.04045 [hep- ph]
Pith/arXiv arXiv 2023
-
[14]
F. Gao, K. Huang, Y. Lu, and Y. Liu, From Magnetic to Inverse Magnetic Catalysis: The Interplay of Quark and Gluon Mass Generation in Magnetic Fields (2026), arXiv:2606.23736 [hep-ph]
Pith/arXiv arXiv 2026
-
[15]
R. D. Pisarski and F. Wilczek, Remarks on the Chiral Phase Transition in Chromodynamics, Phys. Rev. D29, 338 (1984)
1984
-
[16]
J. P. Klinger, R. Kaiser, O. Philipsen, and J. Schaible, On the phase structure of massless many-flavour QCD with staggered fermions, in42th International Symposium on Lattice Field Theory(2026) arXiv:2603.20099 [hep-lat]
arXiv 2026
-
[17]
M. A. Stephanov, QCD Phase Diagram and the Criti- cal Point, Prog. Theor. Phys. Suppl.153, 139 (2004), arXiv:hep-ph/0402115
Pith/arXiv arXiv 2004
-
[18]
M. A. Stephanov, K. Rajagopal, and E. V. Shuryak, Sig- natures of the tricritical point in QCD, Phys. Rev. Lett. 81, 4816 (1998), arXiv:hep-ph/9806219
Pith/arXiv arXiv 1998
-
[19]
V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Dy- namical flavor symmetry breaking by a magnetic field in 2+1 dimensions, Physical Review D52, 4718–4735 (1995)
1995
-
[20]
F. Preis, A. Rebhan, and A. Schmitt, Inverse magnetic catalysis in dense holographic matter, JHEP03, 033, arXiv:1012.4785 [hep-th]
-
[21]
G. S. Bali, F. Bruckmann, G. Endrodi, Z. Fodor, S. D. Katz, S. Krieg, A. Schafer, and K. K. Szabo, The QCD phase diagram for external magnetic fields, JHEP02, 044, arXiv:1111.4956 [hep-lat]
-
[23]
K.-I. Aoki, H. Uoi, and M. Yamada, Functional renormal- ization group study of the Nambu–Jona-Lasinio model at finite temperature and density in an external magnetic field, Phys. Lett. B753, 580 (2016), arXiv:1507.02527 [hep-ph]
Pith/arXiv arXiv 2016
-
[24]
V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Catalysis of dynamical flavor symmetry breaking by a magnetic field in (2+1)-dimensions, Phys. Rev. Lett.73, 3499 (1994), [Erratum: Phys.Rev.Lett. 76, 1005 (1996)], arXiv:hep-ph/9405262
Pith/arXiv arXiv 1994
-
[25]
J.-L. Kneur, M. B. Pinto, and R. O. Ramos, Phase diagram of the magnetized planar Gross-Neveu model beyond the large-N approximation, Phys. Rev. D88, 045005 (2013), arXiv:1306.2933 [hep-ph]
Pith/arXiv arXiv 2013
-
[26]
J. J. Lenz, M. Mandl, and A. Wipf, The magnetized (2+1)-dimensional Gross-Neveu model at finite density, Physical Review D108, 074508 (2023), arXiv:2304.14812 [hep-lat]
Pith/arXiv arXiv 2023
-
[27]
D. D. Scherer and H. Gies, Renormalization group study of magnetic catalysis in the 3dgross-neveu model, Phys. Rev. B85, 195417 (2012)
2012
-
[28]
K. Fukushima and J. M. Pawlowski, Magnetic catalysis in hot and dense quark matter and quantum fluctuations, Phys. Rev. D86, 076013 (2012), arXiv:1203.4330 [hep- ph]
Pith/arXiv arXiv 2012
-
[29]
D. J. Gross and A. Neveu, Dynamical symmetry breaking in asymptotically free field theories, Phys. Rev. D10, 3235 (1974)
1974
-
[30]
Wetterich, Exact evolution equation for the effective potential, Phys
C. Wetterich, Exact evolution equation for the effective potential, Phys. Lett. B301, 90 (1993), arXiv:1710.05815 [hep-th]
Pith/arXiv arXiv 1993
-
[31]
E. Grossi and N. Wink, Resolving phase transitions with discontinuous Galerkin methods, SciPost Phys. Core6, 071 (2023), arXiv:1903.09503 [hep-th]
Pith/arXiv arXiv 2023
-
[32]
J. Stoll, N. Zorbach, A. Koenigstein, M. J. Steil, and S. Rechenberger, Bosonic fluctuations in the (1 + 1)- dimensional Gross-Neveu(-Yukawa) model at varyingµ andTand finiteN, arXiv e-prints , arXiv:2108.10616 (2021), arXiv:2108.10616 [hep-ph]
Pith/arXiv arXiv 2021
-
[33]
Zinn-Justin, Four fermion interaction near four- dimensions, Nucl
J. Zinn-Justin, Four fermion interaction near four- dimensions, Nucl. Phys. B367, 105 (1991)
1991
-
[34]
D. F. Litim, Optimized renormalization group flows, Phys. Rev. D64, 105007 (2001), arXiv:hep-th/0103195
Pith/arXiv arXiv 2001
-
[35]
I. A. Shovkovy, Magnetic Catalysis: A Review (2013) pp. 13–49, arXiv:1207.5081 [hep-ph]
Pith/arXiv arXiv 2013
-
[36]
Kurganov and E
A. Kurganov and E. Tadmor, New high-resolution cen- tral schemes for nonlinear conservation laws and con- vection–diffusion equations, Journal of Computational Physics160, 241 (2000)
2000
-
[37]
Nijholt, J
B. Nijholt, J. Weston, J. Hoofwijk, and A. Akhmerov, Adaptive: parallel active learning of mathematical func- tions (2019)
2019
-
[38]
Bandyopadhyay and R
A. Bandyopadhyay and R. L. S. Farias, Inverse magnetic catalysis: how much do we know about?, The European Physical Journal Special Topics230, 719–728 (2021)
2021
-
[39]
Fraga, B
E. Fraga, B. Mintz, and J. Schaffner-Bielich, A search for inverse magnetic catalysis in thermal quark–meson models, Physics Letters B731, 154 (2014)
2014
-
[40]
B.-J. Schaefer and J. Wambach, The Phase diagram of the quark meson model, Nucl. Phys. A757, 479 (2005), arXiv:nucl-th/0403039
Pith/arXiv arXiv 2005
-
[41]
R.-A. Tripolt, B.-J. Schaefer, L. von Smekal, and J. Wambach, Low-temperature behavior of the quark- meson model, Phys. Rev. D97, 034022 (2018), arXiv:1709.05991 [hep-ph]
Pith/arXiv arXiv 2018
-
[42]
J. O. Andersen and A. Tranberg, The Chiral transi- tion in a magnetic background: Finite density effects and the functional renormalization group, JHEP08, 002, arXiv:1204.3360 [hep-ph]
-
[43]
Skokov, Phase diagram in an external magnetic field beyond a mean-field approximation, Phys
V. Skokov, Phase diagram in an external magnetic field beyond a mean-field approximation, Phys. Rev. D85, 034026 (2012), arXiv:1112.5137 [hep-ph]
Pith/arXiv arXiv 2012
-
[44]
M. J. Steil and A. Koenigstein, Numerical fluid dynam- ics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. III. Shock and rarefaction waves in RG flows reveal limitations of the N→∞limit in O(N)-type models, Phys. Rev. D106, 065014 (2022), arXiv:2108.04037 [cond-mat.stat-mech]
Pith/arXiv arXiv 2022
-
[45]
Virtanen, R
P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Pe- terson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, ˙I. Po- lat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henr...
2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.