T times μ phase diagram from a fractal NJL model
Pith reviewed 2026-05-18 18:11 UTC · model grok-4.3
The pith
A μ-dependent coupling in the fractal NJL model produces a T-μ phase diagram compatible with lattice QCD and STAR data.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By fitting a μ-dependent coupling to lattice pseudo-critical temperatures at low μ and applying it to the fractal NJL model, the authors obtain a single parameter set that describes the T × μ phase diagram in good agreement with STAR data from heavy-ion collisions, for both extensive and non-extensive statistics.
What carries the argument
The μ-dependent coupling in the fractal NJL model, obtained by fitting lattice pseudo-critical temperatures at low μ, which accounts for gluon effects and allows direct computation of the dynamical mass, condensate, susceptibility, and phase boundary.
If this is right
- The model computes the dynamical quark mass and quark condensate at finite T and μ.
- It locates the phase transition via the thermal susceptibility in the T-μ plane.
- A single adjusted parameter set reproduces STAR data for both extensive and non-extensive statistics.
- The resulting phase diagram remains compatible with lattice QCD results at low μ.
Where Pith is reading between the lines
- Parameterizing gluon effects through a simple μ-dependent coupling may offer a practical way to extend other effective models into the high-density regime.
- The close match to STAR data suggests that fractal features in the quark interaction could capture relevant non-perturbative QCD physics.
- Predictions for the critical endpoint position could be compared against upcoming beam-energy scan results.
- The fitting procedure itself could be repeated for related models such as the Polyakov-NJL to test consistency across frameworks.
Load-bearing premise
The μ-dependence of the coupling fixed by fitting lattice data at low chemical potential can be transferred unchanged to compute the phase boundary at higher μ.
What would settle it
A lattice calculation or future heavy-ion measurement of the pseudo-critical temperature at a chemical potential well above the fitted range that lies substantially off the model's predicted curve.
Figures
read the original abstract
We propose a $\mu$-dependent coupling for a fractal effective model (FNJL) to make the results for the phase diagram compatible with the experimental data and lattice QCD calculations. The $\mu$-dependence of the coupling, which accounts for gluon effects, is obtained by fitting the lattice QCD results for the pseudo-critical temperature with the fractal model. We then use the new effective coupling in order to compute the dynamical mass, the quark condensate, the thermal susceptibility and, finally, the $T\times\mu$ phase diagram. We consider both extensive and non-extensive statistics, and with a slight variation in the $\mu$-dependent coupling parameters we provide a single result for our model which is able to describe incredibly well the data from STAR, considering the simplicity of the effective model.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a μ-dependent coupling in a fractal NJL (FNJL) model, obtained by fitting lattice QCD pseudo-critical temperatures at low μ, to compute the dynamical quark mass, condensate, thermal susceptibility, and the full T-μ phase diagram. Both extensive (Boltzmann-Gibbs) and non-extensive (Tsallis) statistics are considered. With a slight additional variation of the coupling parameters, the model is claimed to describe STAR freeze-out data well, offering a simple effective description compatible with lattice results.
Significance. If the extrapolation of the fitted μ-dependent coupling holds without further tuning, the construction supplies a compact effective model that interpolates lattice pseudo-critical temperatures to finite density and yields a phase boundary consistent with heavy-ion data. The fractal extension and non-extensive statistics constitute modest novelties, but the overall significance remains limited by the absence of independent validation for the coupling outside the fit window.
major comments (2)
- [Definition of μ-dependent coupling] § on definition of μ-dependent coupling: The functional form of G(μ) is fixed exclusively by fitting lattice pseudo-critical temperatures at small μ; the manuscript supplies neither a QCD-derived motivation for this functional dependence nor a cross-check against any other lattice observable (e.g., susceptibilities or condensates) at intermediate μ, rendering the subsequent phase-diagram computation at higher μ an extrapolation whose reliability is untested.
- [Abstract and STAR comparison] Abstract and results section on STAR comparison: The claim that 'with a slight variation in the μ-dependent coupling parameters' the model describes STAR data 'incredibly well' introduces an extra fitting step after the lattice calibration; this additional freedom weakens the assertion that the model predicts rather than reproduces the observed freeze-out points and increases the circularity burden on the central claim.
minor comments (2)
- [Figures] The legends and curve styles in the phase-diagram figures should explicitly label which lines correspond to the original lattice-fitted parameters versus the slightly varied set used for the STAR comparison.
- [Thermal susceptibility] Clarify whether the reported pseudo-critical temperatures are defined via the peak of the susceptibility or via a fixed value of the condensate; inconsistent definitions would affect the fitting procedure.
Simulated Author's Rebuttal
We thank the referee for the thorough review and constructive comments. We address each major point below and have revised the manuscript to improve clarity on the model's construction and claims.
read point-by-point responses
-
Referee: [Definition of μ-dependent coupling] § on definition of μ-dependent coupling: The functional form of G(μ) is fixed exclusively by fitting lattice pseudo-critical temperatures at small μ; the manuscript supplies neither a QCD-derived motivation for this functional dependence nor a cross-check against any other lattice observable (e.g., susceptibilities or condensates) at intermediate μ, rendering the subsequent phase-diagram computation at higher μ an extrapolation whose reliability is untested.
Authors: We agree that G(μ) is determined phenomenologically by fitting lattice pseudo-critical temperatures at low μ. As an effective model, the FNJL framework does not derive the functional form from first-principles QCD; instead, the μ-dependence is introduced to capture gluon screening effects not present in the standard NJL model. We have added a dedicated paragraph in the revised manuscript explaining this phenomenological motivation, the specific choice of functional form, and the limitations of extrapolating to higher μ. While the model computes the condensate and susceptibility, we have not performed direct comparisons to lattice results at intermediate μ, as the present scope focuses on the phase diagram; this is noted as a limitation for future work. revision: yes
-
Referee: [Abstract and STAR comparison] Abstract and results section on STAR comparison: The claim that 'with a slight variation in the μ-dependent coupling parameters' the model describes STAR data 'incredibly well' introduces an extra fitting step after the lattice calibration; this additional freedom weakens the assertion that the model predicts rather than reproduces the observed freeze-out points and increases the circularity burden on the central claim.
Authors: We acknowledge that the slight parameter variation for the STAR comparison constitutes an additional adjustment after the lattice calibration. This step was intended to illustrate the model's flexibility in describing experimental data. In the revised manuscript we have moderated the language in the abstract (removing 'incredibly well'), clarified that the primary calibration remains the lattice fit, and presented the STAR agreement as a consistency check within the uncertainties of the effective parameters rather than an independent prediction. The text now emphasizes that the phase diagram is obtained from the lattice-constrained coupling, with the data comparison serving to assess overall compatibility. revision: yes
Circularity Check
μ-dependent coupling fitted to lattice pseudo-critical temperatures then inserted into gap equation to obtain phase diagram and STAR agreement
specific steps
-
fitted input called prediction
[Abstract]
"The μ-dependence of the coupling, which accounts for gluon effects, is obtained by fitting the lattice QCD results for the pseudo-critical temperature with the fractal model. We then use the new effective coupling in order to compute the dynamical mass, the quark condensate, the thermal susceptibility and, finally, the T×μ phase diagram. ... with a slight variation in the μ-dependent coupling parameters we provide a single result for our model which is able to describe incredibly well the data from STAR"
G(μ) parameters are fixed by fitting lattice Tc(μ) at small μ; the same G(μ) is inserted into the gap equation to produce the critical line whose low-μ values must reproduce the fitted Tc by construction. The subsequent claim that this parametrization 'describes incredibly well' the STAR data therefore rests on an extrapolation whose functional form was already constrained by the lattice fit rather than derived from first-principles QCD.
full rationale
The derivation chain begins by fitting the functional form of the μ-dependent coupling to lattice QCD pseudo-critical temperatures at low μ. This fitted coupling is then substituted into the model's gap equation to generate the full T-μ critical line, dynamical masses, and condensates. The resulting phase boundary is presented as describing STAR freeze-out data well after minor parameter adjustment. Because the low-μ segment of the critical line is fixed by the same fit used for lattice Tc, agreement there is tautological; the finite-μ extrapolation inherits the parametrization without an independent QCD derivation or cross-validation against other lattice observables. This matches the fitted-input-called-prediction pattern with partial circularity on the central claim, but the model retains some independent structure outside the fit window, preventing a higher score.
Axiom & Free-Parameter Ledger
free parameters (1)
- μ-dependent coupling parameters
axioms (1)
- domain assumption The fractal NJL model captures the essential non-perturbative dynamics of QCD at finite temperature and density.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/FractalStructuresYangMillsthermofractal_q_derivation echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
q=1+3/(11Nc−2Nf). ... Geff(p)=Gq(1+(q−1)Ep/λ)^(−1/(q−1))
-
IndisputableMonolith/Cost/FunctionalEquationwashburn_uniqueness_aczel refines?
refinesRelation between the paper passage and the cited Recognition theorem.
The μ-dependence of the coupling... is obtained by fitting the lattice QCD results for the pseudo-critical temperature
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
The first second of the Universe,
D. Schwarz, “The first second of the Universe,”Annalen der Physik, vol. 515, p. 220–270, May 2003
work page 2003
-
[2]
The equation of state of hot, dense matter and neutron stars,
J. M. Lattimer and M. Prakash, “The equation of state of hot, dense matter and neutron stars,”Physics Reports, vol. 621, p. 127–164, Mar. 2016
work page 2016
-
[3]
Possible Generalization of Boltzmann-Gibbs Statistics,
C. Tsallis, “Possible Generalization of Boltzmann-Gibbs Statistics,”J. Statist. Phys., vol. 52, pp. 479–487, 1988
work page 1988
-
[4]
Tsallis-thermometer: a qgp indicator for large and small collisional systems,
G. Bíró, G. G. Barnaföldi, and T. S. Biró, “Tsallis-thermometer: a qgp indicator for large and small collisional systems,”Journal of Physics G: Nuclear and Particle Physics, vol. 47, p. 105002, Aug. 2020
work page 2020
-
[5]
G. Wilk and Z. Włodarczyk, “Interpretation of the nonextensivity paramete q in some applications of tsallis statistics and lévy distributions,”Physical Review Letters, vol. 84, p. 2770–2773, Mar. 2000
work page 2000
-
[6]
C. Beck and E. G. D. Cohen, “Superstatistics,”Physica A: Statistical Me- chanics and its Applications, vol. 322, pp. 267–275, 2003
work page 2003
-
[7]
A. Bialas and R. Peschanski, “Moments of rapidity distributions as a mea- sure of short-range fluctuations in high-energy collisions,”Nuclear Physics B, vol. 273, p. 703–718, Sept. 1986
work page 1986
-
[8]
K. Aamodt and et al., “Transverse momentum spectra of charged particles in proton–proton collisions at 900 gev with alice at the lhc,”Physics Letters B, vol. 693, p. 53–68, Sept. 2010
work page 2010
-
[9]
Asymptotic sum rules at infinite momentum,
J. D. Bjorken, “Asymptotic sum rules at infinite momentum,”Physical Review, vol. 179, p. 1547–1553, Mar. 1969
work page 1969
-
[10]
The renormalization group and theϵexpan- sion,
K. G. Wilson and J. Kogut, “The renormalization group and theϵexpan- sion,”Physics Reports, vol. 12, no. 2, pp. 75–199, 1974
work page 1974
-
[11]
Self-consistency in non-extensive thermodynamics of highly excited hadronic states,
A. Deppman, “Self-consistency in non-extensive thermodynamics of highly excited hadronic states,”Physica A - Statistical Mechanics and Its Appli- cations, vol. 391, no. 24, pp. 6380–6385, 2012. 13
work page 2012
-
[12]
A. Deppman, E. Megías, and D. P. Menezes, “Fractal Structures of Yang- Mills Fields and Non Extensive Statistics: Applications to High Energy Physics,”MDPI Physics, vol. 2, no. 3, pp. 455–480, 2020
work page 2020
-
[13]
Relativistic thermodynamics: Transverse mo- mentum distributions in high-energy physics,
J. Cleymans and D. Worku, “Relativistic thermodynamics: Transverse mo- mentum distributions in high-energy physics,”Journal of Physics G: Nu- clear and Particle Physics, vol. 40, no. 9, p. 095006, 2013
work page 2013
-
[14]
Tsallis statistics in the MIT bag model,
E. Megías, D. P. Menezes, and A. Deppman, “Tsallis statistics in the MIT bag model,”Nuclear Physics A, vol. 987, pp. 144–158, 2019
work page 2019
-
[15]
Nambu–Jona- Lasinio model with a fractal inspired coupling,
E. Megías, M. J. Teixeira, V. S. Timóteo, and A. Deppman, “Nambu–Jona- Lasinio model with a fractal inspired coupling,”Phys. Lett. B, vol. 860, p. 139192, 2025
work page 2025
-
[16]
A Nonlocal, covariant generalization of the NJL model,
R. D. Bowler and M. C. Birse, “A Nonlocal, covariant generalization of the NJL model,”Nucl. Phys. A, vol. 582, pp. 655–664, 1995
work page 1995
-
[17]
D. Gomez Dumm and N. N. Scoccola, “Chiral quark models with nonlocal separable interactions at finite temperature and chemical potential,”Phys. Rev. D, vol. 65, p. 074021, 2002
work page 2002
-
[18]
Dynamics and ther- modynamics of a non-local PNJL model with running coupling,
T. Hell, S. Roessner, M. Cristoforetti, and W. Weise, “Dynamics and ther- modynamics of a non-local PNJL model with running coupling,”Phys. Rev. D, vol. 79, p. 014022, 2009
work page 2009
-
[19]
Superstrong coupling NJL model in arbitrary spacetime dimensions,
A. Ahmad, A. Martínez, and A. Raya, “Superstrong coupling NJL model in arbitrary spacetime dimensions,”Phys. Rev. D, vol. 98, no. 5, p. 054027, 2018
work page 2018
-
[20]
The Nambu and Jona Lasinio model: Its implica- tions for hadrons and nuclei,
U. Vogl and W. Weise, “The Nambu and Jona Lasinio model: Its implica- tions for hadrons and nuclei,”Prog. Part. Nucl. Phys., vol. 27, pp. 195–272, 1991
work page 1991
-
[21]
The Nambu-Jona-Lasinio model of quantum chromody- namics,
S. P. Klevansky, “The Nambu-Jona-Lasinio model of quantum chromody- namics,”Rev. Mod. Phys., vol. 64, pp. 649–708, 1992
work page 1992
-
[22]
J. M. Conroy, H. G. Miller, and A. R. Plastino, “Thermodynamic Consis- tency of theq-Deformed Fermi-Dirac Distribution in Nonextensive Ther- mostatics,”Phys. Lett. A, vol. 374, pp. 4581–4584, 2010
work page 2010
-
[23]
Non extensive thermodynamics and neutron star properties,
D. P. Menezes, A. Deppman, E. Megías, and L. B. Castro, “Non extensive thermodynamics and neutron star properties,”Eur. Phys. J. A, vol. 51, no. 12, p. 155, 2015
work page 2015
-
[24]
Thermodynamics and critical behavior in the Nambu-Jona-Lasinio model of QCD,
P. Costa, M. C. Ruivo, and C. A. de Sousa, “Thermodynamics and critical behavior in the Nambu-Jona-Lasinio model of QCD,”Phys. Rev. D, vol. 77, p. 096001, 2008
work page 2008
-
[25]
Nonextensive effects in the Nambu-Jona-Lasinio model of QCD,
J. Rozynek and G. Wilk, “Nonextensive effects in the Nambu-Jona-Lasinio model of QCD,”J. Phys. G, vol. 36, p. 125108, 2009. 14
work page 2009
-
[26]
Thermo-magnetic effects in quark matter: Nambu–Jona-Lasinio model constrained by lattice QCD,
R. L. S. Farias, V. S. Timóteo, S. S. Avancini, M. B. Pinto, and G. Krein, “Thermo-magnetic effects in quark matter: Nambu–Jona-Lasinio model constrained by lattice QCD,”Eur. Phys. J. A, vol. 53, no. 5, p. 101, 2017
work page 2017
-
[27]
π0 pole mass calculation in a strong magnetic field and lattice constraints,
S. S. Avancini, R. L. S. Farias, M. Benghi Pinto, W. R. Tavares, and V. S. Timóteo, “π0 pole mass calculation in a strong magnetic field and lattice constraints,”Phys. Lett. B, vol. 767, pp. 247–252, 2017
work page 2017
-
[28]
W. R. Tavares, R. L. S. Farias, S. S. Avancini, V. S. Timóteo, M. B. Pinto, and G. Krein, “Nambu–Jona-Lasinio SU(3) model constrained by lattice QCD: thermomagnetic effects in the magnetization,”Eur. Phys. J. A, vol. 57, no. 9, p. 278, 2021
work page 2021
-
[29]
Chi- ral vortical catalysis constrained by LQCD simulations,
R. M. Nunes, R. L. S. Farias, W. R. Tavares, and V. S. Timóteo, “Chi- ral vortical catalysis constrained by LQCD simulations,”Phys. Rev. D, vol. 111, no. 5, p. 056026, 2025
work page 2025
-
[30]
Chiral effective model with the Polyakov loop,
K. Fukushima, “Chiral effective model with the Polyakov loop,”Phys. Lett. B, vol. 591, pp. 277–284, 2004
work page 2004
-
[31]
Polyakov loop in chiral quark models at finite temperature,
E. Megías, E. Ruiz Arriola, and L. L. Salcedo, “Polyakov loop in chiral quark models at finite temperature,”Phys. Rev. D, vol. 74, p. 065005, 2006
work page 2006
-
[32]
Phases of QCD: Lattice thermo- dynamics and a field theoretical model,
C. Ratti, M. A. Thaler, and W. Weise, “Phases of QCD: Lattice thermo- dynamics and a field theoretical model,”Phys. Rev. D, vol. 73, p. 014019, 2006
work page 2006
-
[33]
The Phase Structure of the Polyakov–Quark-Meson Model,
B.-J. Schaefer, J. M. Pawlowski, and J. Wambach, “The Phase Structure of the Polyakov–Quark-Meson Model,”Phys. Rev. D, vol. 76, p. 074023, 2007
work page 2007
-
[34]
Phase diagrams in the three-flavor Nambu-Jona-Lasinio model with the Polyakov loop,
K. Fukushima, “Phase diagrams in the three-flavor Nambu-Jona-Lasinio model with the Polyakov loop,”Phys. Rev. D, vol. 77, p. 114028, 2008. [Erratum: Phys.Rev.D 78, 039902 (2008)]
work page 2008
-
[35]
QCD Crossover at Finite Chem- ical Potential from Lattice Simulations,
S. Borsanyi, Z. Fodor, J. N. Guenther, R. Kara, S. D. Katz, P. Parotto, A. Pasztor, C. Ratti, and K. K. Szabo, “QCD Crossover at Finite Chem- ical Potential from Lattice Simulations,”Phys. Rev. Lett., vol. 125, no. 5, p. 052001, 2020
work page 2020
-
[36]
Hadron production in central nucleus-nucleus collisions at chemical freeze-out,
A. Andronic, P. Braun-Munzinger, and J. Stachel, “Hadron production in central nucleus-nucleus collisions at chemical freeze-out,”Nucl. Phys. A, vol. 772, pp. 167–199, 2006
work page 2006
-
[37]
Hadron Formation in Relativistic Nuclear Collisions and the QCD Phase Diagram,
F. Becattini, M. Bleicher, T. Kollegger, T. Schuster, J. Steinheimer, and R. Stock, “Hadron Formation in Relativistic Nuclear Collisions and the QCD Phase Diagram,”Phys. Rev. Lett., vol. 111, p. 082302, 2013. 15
work page 2013
-
[38]
Freeze-out conditions from net-proton and net-charge fluctuations at RHIC,
P. Alba, W. Alberico, R. Bellwied, M. Bluhm, V. Mantovani Sarti, M. Nahrgang, and C. Ratti, “Freeze-out conditions from net-proton and net-charge fluctuations at RHIC,”Phys. Lett. B, vol. 738, pp. 305–310, 2014
work page 2014
-
[39]
V. Vovchenko, V. V. Begun, and M. I. Gorenstein, “Hadron multiplicities and chemical freeze-out conditions in proton-proton and nucleus-nucleus collisions,”Phys. Rev. C, vol. 93, no. 6, p. 064906, 2016
work page 2016
-
[40]
L. Adamczyket al., “Bulk Properties of the Medium Produced in Rela- tivistic Heavy-Ion Collisions from the Beam Energy Scan Program,”Phys. Rev. C, vol. 96, no. 4, p. 044904, 2017
work page 2017
-
[41]
Non extensive thermodynam- ics for hadronic matter with finite chemical potentials,
E. Megías, D. P. Menezes, and A. Deppman, “Non extensive thermodynam- ics for hadronic matter with finite chemical potentials,”Physica A, vol. 421, pp. 15–24, 2015
work page 2015
-
[42]
A. Bazavovet al., “Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model,”Phys. Rev. D, vol. 86, p. 034509, 2012
work page 2012
- [43]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.