Pith. sign in

REVIEW 4 major objections 7 minor 28 references

Thermomagnetic Effects of Quark Matter in the NJL Model: Application of Regularization Schemes

T0 review · 4 major / 7 minor · reviewed 2026-08-11 · deepseek-v4-flash

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

Pith's one-line read The paper claims that in the SU(3) Nambu-Jona-Lasinio model, the choice and application of a regularization scheme changes the qualitative thermodynamics of hot magnetized quark matter, including a possible violation of causality when the…

desk verdict A plausible and useful warning about regularization scheme choice in magnetized NJL—undermined as written because the soft cut-off regulator multiplies already-integrated potentials, leaving the headline curves unreproducible. read the letter →

arxiv 2412.08674 v2 pith:QUAG3MZ7 submitted 2024-12-11 nucl-th

Xiang-Qiong Liu This is my paper · ORCID
classification nucl-th
keywords Nambu-Jona-Lasiniomodelthermomagneticquarkmatterregularizationschemesmagneticfield-independentsoftcutoffPauli-Villarssquaredspeedofsoundcausalityviolation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper asks whether the way ultraviolet divergences are tamed in the SU(3) Nambu-Jona-Lasinio model changes the predicted properties of quark matter at high temperature in a strong magnetic field. It compares three regularization schemes — magnetic field-independent cutoff, soft cutoff, and Pauli–Villars — and finds that they disagree not only in the size of corrections but in the qualitative shape of the thermodynamics: the pressure versus magnetic field rises then falls, falls then rises, or oscillates depending on the scheme, and the strange quark mass responds differently above $T\simeq 200$ MeV. The sharpest claim is that when the regularization is also applied to the temperature-dependent part of the soft cutoff and Pauli–Villars potentials, the pressure decreases at high temperature and the squared speed of sound $c_s^2$ grows sharply, which the paper says can lead to violations of causality. The upshot is that the choice of regularization is a first-order modeling decision in NJL studies of magnetized quark matter, not a minor technical detail.

What carries the argument

The central object is the mean-field thermodynamic potential of the SU(3) NJL model, split into vacuum, magnetic-field, and finite-temperature contributions, together with the gap equation that fixes the dynamical masses. The three regulators are the magnetic-field-independent cutoff expressed through Hurwitz zeta functions, the soft cutoff $f_\Lambda(p)=\left[1+\exp\left(\left(\sqrt{p^2+2|q_f|B\,n}-\Lambda\right)/(0.05\Lambda)\right)\right]^{-1}$, and the Pauli–Villars combination $f_\mathrm{PV}(M_f)=\sum_{j=0}^3 c_j\sqrt{M_f^2+j\Lambda^2}$ with $c_0=1$, $c_1=-3$, $c_2=3$, $c_3=-1$. The decisive mechanism is whether the regulator multiplies the whole potential, including the thermal-magnetic part, as in Eqs. (28) and (31), or only the vacuum part, as in Eqs. (29) and (32).

What would settle it

Recompute the soft cutoff thermodynamic potential and condensate with $f_\Lambda(p)$ placed inside each momentum integral and Landau sum before the integration and summation are performed, using the same parameters; if $P(T,eB)$ then rises with temperature at high $T$ and $c_s^2$ stays at or below $1/3$, the paper's central causality conclusion does not survive the corrected implementation.

Watch

Extended reading notes

Core claim

The paper's central claim is that the regularization scheme in the SU(3) NJL model changes the qualitative behavior of thermodynamic quantities in hot, magnetized quark matter, not just their numerical size. Starting from the mean-field thermodynamic potential $\Omega=\Omega_\text{vac}+\Omega_\text{mag}+\Omega_\text{med}$ and the gap equations for the constituent masses $M_u$, $M_d$, $M_s$, the paper computes pressure, energy density, specific heat, and squared speed of sound under three schemes. For the soft cutoff and Pauli–Villars schemes it distinguishes two variants: regulating only the vacuum part, and regulating the full temperature-dependent potential. The numerical comparison shows that the vacuum-only variants stay close to the magnetic-field-independent baseline, while the fully regulated variants make the pressure fall at high temperature, make energy density and specific heat drop instead of equilibrating, and drive $c_s^2$ sharply upward — behavior the paper reads as a possible violation of causality. The conclusion is that regularizing the temperature-dependent part is not a safe default, and the scheme choice must be made consciously for the physical question at hand.

Load-bearing premise

The paper's high-temperature conclusions rely on applying the soft cutoff regulator $f_\Lambda(p)$ as a multiplicative prefactor to already-integrated potentials in Eqs. (28) and (29), even though the regulator is only defined inside the momentum integral and Landau sum, so the regulated potential is not well defined unless the regulator is inserted at an earlier stage.

Editorial extensions

If this is right

  • If the claim holds, NJL calculations of magnetized quark matter should report which scheme was used and justify it, because the schemes give qualitatively different $P(eB)$, $c_s^2$, and strange quark mass behavior.
  • Applying regularization to the temperature-dependent part of the soft cutoff and Pauli–Villars potentials produces a falling high-temperature pressure and a sharply rising $c_s^2$; those variants should not be used for high-$T$ thermodynamics without checking causality.
  • Without temperature-part regularization, the soft cutoff and Pauli–Villars results are close to the MFIR results for energy density, specific heat, and speed of sound, with deviations concentrated in the strange quark sector above $T\simeq200$ MeV.
  • MFIR is the numerically simplest scheme and separates magnetic and non-magnetic contributions cleanly, so it is the natural starting point for exploratory NJL calculations.
  • There is no universal scheme: the appropriate choice depends on the quantity being studied, with the Pauli–Villars variant without temperature regularization favored for extensions such as the anomalous magnetic moment.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the scheme dependence is as strong as reported, then NJL-based conclusions about magnetic catalysis or inverse catalysis should come with a scheme-sensitivity test; otherwise the qualitative outcome may be an artifact of the cutoff.
  • A cleaner soft cutoff implementation would place $f_\Lambda(p)$ inside every momentum integrand and Landau sum before integrating; doing so would remove the ambiguity in Eqs. (28)–(29) and would test directly whether the high-temperature rise in $c_s^2$ survives.
  • The sharp rise of $c_s^2$ past $1/3$ when the temperature part is regulated may indicate that the fully regulated variant double-counts vacuum contributions at finite temperature rather than a physical instability; a flow-equation or renormalization-group regulator could separate those effects.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 7 minor

Summary. The manuscript studies three regularization schemes—magnetic field-independent regularization (MFIR), soft cut-off, and Pauli–Villars—in the SU(3) NJL model for quark matter at finite temperature and magnetic field. It writes the thermodynamic potential and gap equation in each scheme, solves the gap equation numerically, and compares quark masses, pressure, thermal susceptibility, energy density, specific heat, and squared sound speed as functions of T and eB. The central assertion is that the regularization choice is a first-order modeling decision: applying the soft cut-off or Pauli–Villars regulator also to the temperature-dependent part makes the pressure fall at high T and the sound speed rise sharply, with an associated remark about possible causality violations.

Significance. The question is legitimate and timely for effective-model studies of magnetized quark matter; different regularization prescriptions can indeed change thermodynamics, and the manuscript provides a useful catalog of parameter sets and standard MFIR formulas. The paper does not suffer from circularity in the fitted-input/output sense, since the parameters are taken from earlier fits and the conclusions are not obtained by fitting the outputs. However, the paper's quantitative value is undermined by two underdetermined regularization implementations and by comparing the schemes under different model parameters. If these issues were repaired, the comparison could be a useful reference; as it stands, the main quantitative message is not reproducible from the written equations, and no code or data tables are provided for independent checking.

major comments (4)
  1. [Section III.B, Eqs. (25)-(29)] The regulator f_Λ(p) is defined in Eq. (25) as a function of the integration variable p, the Landau level n, and the flavor through |q_f|. In Eqs. (28) and (29) it is used as a multiplicative prefactor multiplying Ω_vac, Ω_Tmag, φ_vac, and φ_Tmag, which have already been integrated over p and summed over n. Once the integrals and sums are performed, f_Λ(p) has no definite numerical value, so the 'Soft(With Regularization)' and 'Soft(Without Regularization)' results in Figs. 2 and 4-8 are not defined by the equations in the paper. The explicitly regulated integrands, with f_Λ inserted before ∫dp and Σ_n, are never shown; without them the soft-cut-off curves cannot be reproduced. This is load-bearing because the paper's high-temperature conclusions about P, ε, C_V, and c_s² depend directly on those curves.
  2. [Section III.C, Eqs. (30)-(32)] The Pauli–Villars regulator is written as f_P.V(M_f)=Σ_j c_j f(√(M_f²+jΛ²)) and then applied as a scalar f_P.V(M_f) multiplying already-integrated potentials. The standard Pauli–Villars construction is a linear combination of the full potential evaluated at shifted masses, Ω_PV(M)=Σ_j c_j Ω(√(M²+jΛ²)); the paper does not specify which operation is meant, and the notation c_j f(...) is not dimensionally transparent because the function f is not defined. Consequently the Pauli–Villars curves in Figs. 3-8 are unreproducible from the written equations, and the comparison of 'with' versus 'without' regularization in Eq. (31) versus Eq. (32) rests on an undefined prescription.
  3. [Section IV, parameter sets] The three schemes are compared with different model parameters: MFIR and soft cut-off use Λ=631.4 MeV, m_s=135.7 MeV, G=1.835/Λ², K=9.29/Λ², while Pauli–Villars uses Λ=781.2 MeV, m_s=236.9 MeV, G=4.90/Λ², K=129.8/Λ². The observed differences among the schemes in Figs. 1-8 are therefore not attributable solely to the regularization scheme; they could reflect the different parameter sets. Since the paper's central claim is that the choice of regularization scheme itself matters, this parameter/scheme confounding must be addressed, for example by repeating the comparison with a common parameter set or by demonstrating that the qualitative differences survive independent parameter variations.
  4. [Section IV, Fig. 8 and surrounding text] The statement that regularizing the temperature-dependent part 'can lead to violations of causality' is not supported by the displayed data: the vertical axis of Fig. 8 extends only to c_s²=0.5, far below the causal bound c_s²≤1. The figure may show a sharp increase in the sound speed, but not a causality violation. The claim should either be backed by results exceeding unity or removed and replaced by a more cautious statement about the scheme dependence of c_s².
minor comments (7)
  1. [Section III.A, Eq. (19)] The text refers to Γ(x_f) as 'Euler's totient function'; this should be the gamma function.
  2. [Section III.B, Eq. (29)] Eq. (29) is labeled the 'non-regularized form' but still multiplies the vacuum pieces by f_Λ(p); this labeling is internally inconsistent and should be clarified, since the comparison in Figs. 2 and 4-8 depends on what 'without regularization' means.
  3. [Section IV, figure captions] The figure captions do not state the fixed values of eB in Figs. 1-3 and Figs. 5-8, nor are the units of the right-hand panel of Fig. 4 explained beyond the label P×10^9; without these details the curves cannot be fully interpreted or reproduced.
  4. [References [9]-[12]] Several references in the regularization discussion (Refs. [9]-[12], and possibly [16]) appear to concern computer vision, fluid dynamics, or general gauge-theory techniques rather than NJL regularization; the reference list should be rechecked and replaced with appropriate QCD/NJL sources.
  5. [Throughout, Sec. IV] The manuscript repeatedly writes 'Pauli-Villas' instead of 'Pauli-Villars' (e.g., captions of Figs. 3 and 4 and several paragraphs); the spelling should be unified.
  6. [Section III.B, Eqs. (26)-(27)] Equations (26)-(27) repeat Eqs. (3) and (6) without explicitly defining the split into Ω_vac, Ω_Tmag, φ_vac, and φ_Tmag; defining these objects would remove part of the ambiguity in Eqs. (28)-(29).
  7. [Data Availability Statement] The statement 'No Data associated in the manuscript' is at odds with the many numerical figures; providing the numerical data underlying the figures, or at least tabulated values for representative curves, would improve reproducibility.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the regularization comparison is a defined modeling choice with outputs not fed back into the inputs.

full rationale

The paper compares three regularization schemes for the NJL thermodynamic potential and gap equation. Model parameters are quoted from earlier literature (Refs. [27,28]) and are not refit to the quantities being reported (P, epsilon, C_V, c_s^2). The 'regularized' versus 'non-regularized' versions in Eqs. (28)-(29) and (31)-(32) are defined by whether f_Lambda(p) or f_P.V.(M_f) multiplies the temperature-dependent piece; the resulting differences in P(T) and c_s^2(T) are numerical outputs of that defined choice, not quantities used to construct the choice. No equation is defined in terms of its own prediction, and no load-bearing claim is supported only by a self-citation. There is a genuine formal ambiguity: f_Lambda(p) in Eq. (25) depends on the integration variable p and the Landau level n, yet it is used in Eqs. (28)-(29) as a prefactor on already integrated potentials, so the soft-cut-off curves are not uniquely reproducible from the text. This is a correctness/reproducibility concern, not circularity, because the high-temperature behavior and causality remark in Sec. IV are consequences of the chosen regulator rather than inputs that determine it.

Assumptions & free parameters 9 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the NJL mean-field Lagrangian, the Landau-level representation of the magnetic field, the vacuum subtraction convention, and the specific form of each regulator. The model parameters Lambda, m_u, m_d, m_s, G, and K are imported from earlier fits rather than derived here, and they differ between the MFIR/soft and Pauli-Villars runs, which is a direct confound for the paper's scheme comparison. The soft cut-off width 0.05 Lambda is an unvaried hand-chosen number. No new particles, forces, dimensions, or conserved quantities are introduced.

free parameters (9)
  • Momentum cutoff Lambda (MFIR and soft cut-off) = 631.4 MeV
    Taken from Hatsuda and Kunihiro [27]; fitted to vacuum meson observables in previous work.
  • Momentum cutoff Lambda (Pauli-Villars) = 781.2 MeV
    Taken from Carignano and Buballa [28]; differs from the MFIR and soft cut-off value, confounding the scheme comparison.
  • Current quark masses m_u, m_d, m_s (MFIR and soft cut-off) = 5.5, 5.5, 135.7 MeV
    Inputs from Hatsuda and Kunihiro [27] used for the MFIR and soft cut-off runs.
  • Current quark masses m_u, m_d, m_s (Pauli-Villars) = 10.3, 10.3, 236.9 MeV
    Inputs from Carignano and Buballa [28] used for the Pauli-Villars runs.
  • Scalar coupling G (MFIR and soft cut-off) = 1.835 / Lambda^2
    Fitted to vacuum properties in prior literature; used without re-fitting.
  • Scalar coupling G (Pauli-Villars) = 4.90 / Lambda^2
    Fitted to vacuum properties in prior literature; differs substantially from the other scheme's G.
  • Determinant coupling K (MFIR and soft cut-off) = 9.29 / Lambda^2
    Fitted to vacuum properties in prior literature.
  • Determinant coupling K (Pauli-Villars) = 129.8 / Lambda^2
    Fitted to vacuum properties in prior literature; differs substantially from the other scheme's K.
  • Soft cut-off width prefactor 0.05 Lambda = 0.05 Lambda
    Chosen by hand in Eq. (25) to set the smoothness of the regulator; no sensitivity study is performed.
assumptions (4)
  • domain assumption Mean-field approximation for the NJL Lagrangian
    Eq. (2) replaces four-fermion interactions by condensate expectation values; the paper does not quantify corrections beyond mean field.
  • domain assumption Landau-level representation with gauge A_mu = delta_mu2 x_1 B
    Sec. II assumes a static uniform magnetic field in Landau gauge; the paper does not justify the adequacy of this representation for the temperature range studied.
  • domain assumption Vacuum normalization subtraction Omega_eff = Omega(T, mu, M, B) - Omega(0, 0, M, B)
    Eq. (4) sets the vacuum potential to zero by subtraction, a convention that affects the definition of pressure and all derived thermodynamic quantities.
  • ad hoc to paper The soft cut-off regulator f_Lambda(p) can be applied as a multiplicative prefactor to integrated potentials
    Eqs. (28) and (29) treat f_Lambda(p), which depends on momentum p, Landau level n, and charge q_f, as an overall factor on the already-integrated Omega_vac and Omega_Tmag. This is not a standard regulator insertion and is not derived in the paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Thermomagnetic Effects of Quark Matter in the NJL Model: Application of Regularization Schemes." pith.science (2026). https://pith.science/paper/QUAG3MZ7

@misc{pith2026241208674,
  author       = {Pith},
  title        = {Pith review of: Thermomagnetic Effects of Quark Matter in the NJL Model: Application of Regularization Schemes},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QUAG3MZ7}},
  note         = {Machine review of arXiv:2412.08674}
}
read the original abstract

In the SU(3) Nambu-Jona-Lasinio (NJL) model of a thermally magnetized medium, the regularization methods adopted for the thermodynamic potential and the mass gap equation are utilized to calculate the relevant thermodynamic quantities in thermomagnetic quark matter. When dealing with the thermodynamic quantities and the gap equation, three schemes can be chosen, namely magnetic field-independent regularization, soft cut-off regularization, and Pauli-Villars regularization. These three regularization schemes have different influences on the properties of thermomagnetic quark matter, and different choices of regularization schemes will also lead to different effects in the calculation of the properties of thermomagnetic quark matter.

Figures

Figures reproduced from arXiv: 2412.08674 by the authors.

Figure 1
Figure 1. FIG. 1. The relationship between M and T in magnetic-field-independent regularization [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. The relationship between M and T in the soft cutoff scheme: Left panel with regularization, [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. The relationship between M and T in the P.V cut-off scheme: Left panel with regularization, [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Left: The pressure varies with temperature. Right: The pressure varies with magnetic [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 6
Figure 6. Figure 6: FIG. 6. The variation of energy density with [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]
Figure 8
Figure 8. Figure 8: FIG. 8. The variation of sound velocity with [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

28 extracted references · 13 canonical work pages

  1. [9]

    Poggio, V

    T. Poggio, V. Torre, and C. Koch, Readings in computer vision , 638 (1987)

    work page 1987

  2. [11]

    Cortez, L

    R. Cortez, L. Fauci, and A. Medovikov, Physics of Fluids 17 (2005)

    work page 2005

  3. [1]

    S. P. Klevansky, Rev. Mod. Phys. 64, 649 (1992)

    1992

  4. [2]

    Buballa, Phys

    M. Buballa, Phys. Rept. 407, 205 (2005), arXiv:hep-ph/0402234

    arXiv 2005

  5. [3]

    Witten, Phys

    E. Witten, Phys. Rev. D 30, 272 (1984)

    1984

  6. [4]

    Nambu, Phys

    Y. Nambu, Phys. Rev. 117, 648 (1960)

    1960

  7. [5]

    Nambu and G

    Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961)

    1961

  8. [6]

    D. P. Menezes, M. Benghi Pinto, S. S. Avancini, A. Perez Martinez, and C. Providencia, Phys. Rev. C 79, 035807 (2009), arXiv:0811.3361 [nucl-th]

    arXiv 2009

Show all 28 references
  1. [7]

    D. P. Menezes, M. Benghi Pinto, S. S. Avancini, and C. Providencia, Phys. Rev. C 80, 065805 (2009), arXiv:0907.2607 [nucl-th]

    2009 arXiv

  2. [8]

    E. S. Fraga and A. J. Mizher, Phys. Rev. D 78, 025016 (2008), arXiv:0804.1452 [hep-ph]

    2008 arXiv

  3. [10]

    Leibbrandt, Rev

    G. Leibbrandt, Rev. Mod. Phys. 47, 849 (1975)

    1975

  4. [12]

    Girosi, M

    F. Girosi, M. Jones, and T. Poggio, Neural computation 7, 219 (1995)

    1995

  5. [13]

    Mao, Phys

    S. Mao, Phys. Lett. B 758, 195 (2016), arXiv:1602.06503 [hep-ph]

    2016 arXiv

  6. [14]

    Vogl and W

    U. Vogl and W. Weise, Prog. Part. Nucl. Phys. 27, 195 (1991)

    1991

  7. [15]

    Gandhi, A

    A. Gandhi, A. Sharma, R. Pachuau, B. Lalremruata, M. Mehta, P. N. Patil, S. V. Surya- narayana, L. S. Danu, B. K. Nayak, and A. Kumar, Eur. Phys. J. A 57, 1 (2021)

    2021

  8. [16]

    Kunihiro and T

    T. Kunihiro and T. Hatsuda, Phys. Lett. B 206, 385 (1988), [Erratum: Phys.Lett.B 210, 278–278 (1988)]

    1988

  9. [17]

    J. O. Andersen, W. R. Naylor, and A. Tranberg, Rev. Mod. Phys. 88, 025001 (2016), 13 arXiv:1411.7176 [hep-ph]

    2016 arXiv

  10. [18]

    V. A. Miransky and I. A. Shovkovy, Phys. Rept. 576, 1 (2015), arXiv:1503.00732 [hep-ph]

    2015 arXiv

  11. [19]

    F. W. Olver, NIST handbook of mathematical functions hardback and CD-ROM(Cambridge university press, 2010)

    2010

  12. [20]

    P. G. Allen, A. G. Grunfeld, and N. N. Scoccola, Phys. Rev. D 92, 074041 (2015), arXiv:1508.04724 [hep-ph]

    2015 arXiv

  13. [21]

    S. S. Avancini, R. L. S. Farias, M. B. Pinto, T. E. Restrepo, and W. R. Tavares, Phys. Rev. D 103, 056009 (2021), arXiv:2008.10720 [hep-ph]

    2021 arXiv

  14. [22]

    Adel and T

    A. Adel and T. Alharbi, Eur. Phys. J. A 53, 1 (2017)

    2017

  15. [23]

    Chaudhuri, S

    N. Chaudhuri, S. Ghosh, P. Roy, and S. Sarkar, Phys. Rev. D 106, 056020 (2022), arXiv:2209.02248 [hep-ph]

    2022 arXiv

  16. [24]

    Fayazbakhsh and N

    S. Fayazbakhsh and N. Sadooghi, Physical Review D—Particles, Fields, Gravitation, and Cosmology 83, 025026 (2011)

    2011

  17. [25]

    K. Xue, X. Yu, and X. Wang, Chin. Phys. C 46, 013103 (2022), arXiv:2105.14323 [hep-ph]

    2022 arXiv

  18. [26]

    Pauli and F

    W. Pauli and F. Villars, Rev. Mod. Phys. 21, 434 (1949)

    1949

  19. [27]

    Hatsuda and T

    T. Hatsuda and T. Kunihiro, Phys. Rept. 247, 221 (1994), arXiv:hep-ph/9401310

    1994 arXiv

  20. [28]

    Carignano and M

    S. Carignano and M. Buballa, Phys. Rev. D 101, 014026 (2020), arXiv:1910.03604 [hep-ph]. 14

    2020 arXiv

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.