REVIEW 2 major objections 4 minor 1 cited by
The $(2+1)$-dimensional Gross-Neveu-Yukawa model at finite temperature, density, and magnetic field within the Functional Renormalization Group
T0 review · 2 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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.
Extended reading notes
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
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.
Editorial extensions
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.
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.
Assumptions & free parameters
free parameters (1)
- g^2 (four-fermion/Yukawa UV coupling) =
not quoted; tuned to give σ0=1
assumptions (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.
Cite this review
Pith. "Pith review of The $(2+1)$-dimensional Gross-Neveu-Yukawa model at finite temperature, density, and magnetic field within the Functional Renormalization Group." pith.science (2026). https://pith.science/paper/WCOOSU6C
@misc{pith2026260802280,
author = {Pith},
title = {Pith review of: The $(2+1)$-dimensional Gross-Neveu-Yukawa model at finite temperature, density, and magnetic field within the Functional Renormalization Group},
year = {2026},
howpublished = {\url{https://pith.science/paper/WCOOSU6C}},
note = {Machine review of arXiv:2608.02280}
}
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
Forward citations
Cited by 1 Pith paper
-
Magnetic catalysis and Hall conductivity of excitonic insulators in a planar four-Fermi model
In a planar Gross-Neveu model with an excitonic channel, a perpendicular magnetic field enhances the excitonic condensate and raises its critical temperature, and the Hall conductivity develops a characteristic slope ...
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]
-
[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]
arXiv 2013
Show all 44 references
-
[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]
2016 arXiv
-
[10]
Endrodi, Critical point in the QCD phase diagram for extremely strong background magnetic fields, JHEP07, 173, arXiv:1504.08280 [hep-lat]
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]
2015 arXiv
-
[12]
Ferreira, P
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]
2014 arXiv
-
[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]
2023 arXiv
-
[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]
2026 arXiv
-
[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]
2026
-
[17]
M. A. Stephanov, QCD Phase Diagram and the Criti- cal Point, Prog. Theor. Phys. Suppl.153, 139 (2004), arXiv:hep-ph/0402115
2004 arXiv
-
[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
1998 arXiv
-
[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]
Preis, A
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]
2016 arXiv
-
[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
1994 arXiv
-
[25]
Kneur, M
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]
2013 arXiv
-
[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]
2023 arXiv
-
[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]
Fukushima and J
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]
2012 arXiv
-
[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]
1993 arXiv
-
[31]
Grossi and N
E. Grossi and N. Wink, Resolving phase transitions with discontinuous Galerkin methods, SciPost Phys. Core6, 071 (2023), arXiv:1903.09503 [hep-th]
2023 arXiv
-
[32]
Stoll, N
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]
2021 arXiv
-
[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
2001 arXiv
-
[35]
I. A. Shovkovy, Magnetic Catalysis: A Review (2013) pp. 13–49, arXiv:1207.5081 [hep-ph]
2013 arXiv
-
[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]
Schaefer and J
B.-J. Schaefer and J. Wambach, The Phase diagram of the quark meson model, Nucl. Phys. A757, 479 (2005), arXiv:nucl-th/0403039
2005 arXiv
-
[41]
Tripolt, B.-J
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]
2018 arXiv
-
[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]
2012 arXiv
-
[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.040...
2022 arXiv
-
[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...
2020
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.