{"id":"4bf7ab25-c4cd-45bb-8746-baf8c63b2c6b","arxiv_id":"1908.08122","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"An extended Polyakov-Nambu-Jona-Lasinio model fitted to lattice thermodynamics predicts a first-order chiral transition ending at a critical point at T=110 MeV and quark chemical potential 320 MeV.","lead":"This paper builds a modified quark model that matches lattice QCD results at zero density and then predicts where quark matter changes phase at high density. The prediction matters because heavy-ion colliders and neutron star merger simulations need an equation of state across the whole temperature-density plane.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The CEP coordinates rest on a mu-independent T_phen(T) fitted to mu=0 lattice data, whose e/T term is admitted to compensate for missing mesons; no argument shows that compensation survives at finite mu.","rationale":"The reader's conditional verdict is appropriate. The strongest claim is quantitative CEP. I find the same weak point the reader identified, sharpened by the paper's own admission that T_phen is partly a stand-in for missing mesons. This makes the finite-mu extrapolation less secure than the abstract's 'does not require any new parameter' suggests. The mu=0 fit is a necessary but not sufficient validation; kappa agreement only probes small mu. The CEP at mu_q/T about 2.9 is in a regime where the fitted function and omitted channels are unconstrained. This is a correctness risk, not an accusation of error; the model could still be a useful phenomenological EOS. The proposed density-dependent generalization would quantify the sensitivity. No internal contradiction was found in the derivation itself, but the transfer of T_phen to finite mu is the place where the central claim could break. Verdict remains conditional; no adjustment.","tokens_in":14335,"tokens_out":9611,"duration_ms":98672,"concrete_test":"Recompute the CEP with Eq. (28) replaced by the minimal relativistic-density generalization T_phen(T,mu_q)=T_phen(sqrt(T^2+(mu_q/pi)^2)), leaving all other parameters and the solution procedure (Eqs. 13-14 and 34) unchanged. If the endpoint moves by more than about 30 MeV in either T or mu_q, the mu-independence of T_phen is load-bearing and the quoted coordinates are not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the quantitative CEP at T=110 MeV, mu_q=320 MeV, obtained by solving the gap equations and Eq. (34) with the quark-gluon coupling encoded in the fitted function T_phen(T)=a+bT+cT^2+dT^3+e/T (Eq. 28). The load-bearing assumption is not merely that this five-parameter curve is mu-independent; it is that the curve is a physical property of the gluon sector. The paper itself undercuts this: it states that the e/T term 'does not depend on the quark-gluon interaction' and 'compensates partially for the fact that only four types of mesons ... are taken into consideration' (Section III.A, after Eq. 28). Thus T_phen absorbs the missing hadron spectrum and mean-field defects at mu=0. There is no argument, and no test, that this compensation is the same at finite baryon density, where vector mesons, the di-quark sector, and baryonic degrees of freedom enter with different density dependence. The paper's own conclusion says vector mesons 'may change quantitatively the value of the critical temperatures' and that the di-quark sector 'will become important' at low T and high mu. The mu=0 EOS fit and the small-mu Taylor coefficient kappa=0.00989 do not constrain the endpoint: at the CEP, mu_q/T is about 2.9, far beyond the radius of convergence of a quadratic curvature fit (Eq. 30). Consequently the headline coordinates are a transfer of a mu=0 fudge function to an unvalidated regime; if T_phen has any density dependence, or if the missing channels shift the effective coupling, the CEP can move substantially.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the SU(3) Polyakov-Nambu-Jona-Lasinio (PNJL) model in two directions: it includes the next-to-leading order in the 1/Nc expansion, which adds mesonic contributions to the grand potential, and it replaces the constant Polyakov-loop scale T0 with a phenomenological temperature-dependent function T_phen(T)=a+bT+cT^2+dT^3+e/T whose coefficients are fitted to lattice data at zero chemical potential. With this input, the authors report agreement with lattice results for pressure, energy density, entropy density, interaction measure, and speed of sound at μ=0, and with the lattice Taylor-expansion coefficient κ toward finite baryon chemical potential. They then extend the model to finite chemical potential without new parameters, compare with pQCD at large μ, and predict a first-order chiral phase transition ending at a critical endpoint with T_CEP=110 MeV and μ_q=320 MeV. The paper also presents meson and quark masses as functions of T and μ and argues that the first-order transition may be reachable in heavy-ion experiments. The stated goal is to provide an equation of state in the whole (T, μ) plane for applications to heavy-ion collisions and neutron-star physics.","tokens_in":14723,"tokens_out":2860,"duration_ms":28951,"significance":"If the central claim is correct, the paper would deliver a complete, parameter-free-in-μ equation of state that reproduces the μ=0 lattice thermodynamics and the known large-μ pQCD limit, which is a useful input for hydrodynamic simulations and neutron-star modeling. The explicit inclusion of mesonic degrees of freedom at next-to-leading order in Nc and the temperature-dependent quark-gluon coupling are genuine technical improvements over the standard PNJL setup. The finite-μ predictions, including the meson Mott lines and the location of the critical endpoint, are falsifiable against future lattice and experimental data. However, the μ=0 agreement is partly a fit result because the five coefficients of T_phen are tuned to the very lattice curves shown in Figs. 3 and 6, and the paper itself identifies an artifact (the second minimum in the speed of sound) as a limitation of that tuning. The quantitative CEP coordinates rest on the unvalidated assumption that the fitted T_phen function remains valid at all chemical potentials, so the significance of the paper is conditional on addressing that assumption.","major_comments":[{"comment":"The five coefficients a through e in T_phen(T) are determined by comparison with the same lattice gauge calculations whose agreement is then exhibited in Fig. 3. The agreement at μ=0 is therefore a consistency check of the fit, not an independent prediction. The paper explicitly states that the e/T term does not depend on the quark-gluon interaction and 'compensates partially for the fact that only four types of mesons ... are taken into consideration'. That makes T_phen a μ=0 fudge that absorbs missing hadron species and mean-field defects. The central claim that T_phen is a μ-independent physical property of the gluon sector is asserted without any test. Because the CEP lies at μ_q/T≈2.9, far outside the radius of convergence of the quadratic curvature fitted in Eq. (30) (κ=0.00989), the small-μ lattice comparison does not constrain the endpoint. The authors should provide a concrete test of μ-independence, for example by comparing the model with lattice results at imaginary chemical potential or by computing higher-order Taylor coefficients, and should quantify the sensitivity of the CEP to variations in a-e within their fit uncertainties.","section":"§III.A, Eq. (28)"},{"comment":"The speed of sound shows a second minimum not present in the lattice data, and the text after Eq. (29) attributes this to the mean-field level of the gap equations and to the incomplete hadron spectrum. This explicitly confirms that T_phen at μ=0 compensates for model defects rather than representing a pure gluon back-reaction. Since the same fitted T_phen is used unchanged at all chemical potentials, the compensation is not guaranteed to survive at finite density, where vector mesons, baryons, and the di-quark sector enter with different density dependencies. The paper should show how sensitive the μ=0 equation of state is to the compensating term, for instance by repeating the calculation with e=0 or with additional hadronic channels, and should discuss how that sensitivity propagates to the finite-μ phase boundary.","section":"§IV.A, Fig. 6"},{"comment":"The critical endpoint is quoted as T_CEP=110 MeV and μ_q=320 MeV without any estimate of uncertainty. The solution of the six coupled equations in Eq. (34) depends on the specific parametrization of T_phen, on the restriction to pseudoscalar and scalar mesons, and on the mean-field treatment of the gap equations. The authors themselves state that vector mesons 'may change quantitatively the value of the critical temperatures' and that the di-quark sector 'will become important' at low temperature and high chemical potential. Without a parameter-sensitivity scan or an error estimate propagated from the fit of Eq. (28), the quantitative CEP prediction cannot be assessed against other PNJL-type models or against upcoming beam-energy-scan data. At minimum, the authors should show how T_CEP and μ_CEP vary when a-e are varied within their fit uncertainties and when the e/T compensation term is omitted.","section":"§IV.D, Eq. (34)"}],"minor_comments":[{"comment":"The phrase 'Suppress this' appears in the middle of a sentence before 'The temperature of the minimum of the speed of sound...' and should be removed; it is a leftover editorial instruction.","section":"§IV.D, text after Fig. 13"},{"comment":"The caption says 'pQMD calculations' and 'pQMD approach' in two places, but the text and reference [14] refer to pQCD; these are typographical errors that should be corrected to pQCD.","section":"§IV.C, Fig. 8 caption"},{"comment":"The caption begins with 'hase diagram' and should read 'Phase diagram'.","section":"Fig. 13 caption"},{"comment":"Several coefficient names (a0, a1, a2, a3, b3, b4, and the new a, b, c, d, e) are reused for different quantities in Tables I and II; the notation would be clearer if the new coefficients in Eq. (28) were given distinct symbols, since this reuse makes the text harder to follow.","section":"§II.A, Table I and §III.A, Table II"},{"comment":"Reference [24] is cited in the introduction for baryon masses, but the same reference is also listed as [41] later; the numbering should be made consistent.","section":"§I, Introduction"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a relevant problem and contains a useful model extension, but the headline CEP result needs substantially more support before it is suitable for publication. The main issue is not that the model disagrees with consensus, but that a central quantitative prediction is built on a μ=0 fit whose μ-independence is asserted rather than demonstrated. I would encourage the editor to send the paper back for a revision in which the authors provide (i) a direct sensitivity analysis of the CEP to the fitted coefficients, (ii) a test of the μ-independent T_phen assumption via imaginary-μ or higher-order Taylor comparisons, and (iii) a quantitative estimate of the effect of the compensating e/T term. The paper would then be a solid contribution; in its current form, the central claim is underdetermined."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the Polyakov NJL paper by Fuseau, Steinert, and Aichelin. The genuinely new piece is Eq. (28), a temperature-dependent T_phen(T) with five coefficients tuned to lattice data, and the resulting full T-mu phase diagram with a CEP at T=110 MeV, mu_q=320 MeV. The mesonic NLO machinery comes from refs. [16,17] and the back-reaction rescaling from Haas et al., but the combination plus the new parametrization is a real step for this model class.\n\nWhat it does well: with those five coefficients it reproduces the mu=0 pressure, energy density, entropy, interaction measure, and speed of sound over a wide temperature range. The Taylor coefficient kappa=0.00989 falls inside the lattice band. At very large mu the pressure matches pQCD. That is a useful consistency sandwich: lattice at one end, pQCD at the other. The paper is clearly written and the authors flag their own caveats.\n\nSoft spots, in order of size. First, the mu=0 agreement is a fit, not a prediction; the same lattice curves that fix a-e are shown as confirmation in Fig. 3. This is not fatal for a phenomenological model, but it should be said plainly. Second, and more load-bearing, the e/T term in T_phen is admitted to compensate for missing hadron species and mean-field defects. Since T_phen is assumed mu-independent, the CEP coordinates inherit that assumption. Nothing in the paper tests whether the compensation survives at mu_q/T ~ 2.9, far beyond the radius of the quadratic curvature fit. The stress-test note has this right, and the paper's own caveats about vector mesons and the di-quark sector undercut the quantitative endpoint. So I would call the CEP a model outcome, not a robust prediction. Third, the second minimum in the speed of sound is acknowledged but remains a defect in the EoS that matters for hydrodynamic applications. Minor: no code and no uncertainty estimate on the CEP.\n\nI do not think the central argument collapses; for a PNJL calculation this is honest work within the genre. But the headline number should be read as conditional on the mu-independence of an effective parameter fitted at mu=0.\n\nWho benefits: people building EoS tables for heavy-ion and neutron-star simulations, and PNJL practitioners calibrating models. It merits a real referee, with the request that the authors either justify or soften the mu-independence claim and quantify sensitivity of the CEP to the fit choices.\n\nRecommendation: send to peer review; expect revision.","headline":"A competent PNJL extension whose mu=0 floor is fitted, leaving the headline CEP a plausible but unvalidated extrapolation.","tokens_in":15261,"tokens_out":2011,"would_cite":false,"duration_ms":21488,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper predicts a critical end point at T = 110 MeV and $\\mu_q = 320$ MeV from a parameter-free extension of the PNJL model.","keywords":["Polyakov–Nambu–Jona-Lasinio model","QCD phase diagram","critical endpoint","chiral phase transition","quark-gluon coupling","equation of state","finite chemical potential","lattice QCD comparison"],"falsifier":"Compute the model's second-order baryon-number susceptibility and the next Taylor coefficient in $\\mu/T$ and compare with lattice results obtained by analytic continuation or higher-order Taylor expansion. If those coefficients disagree with the model at fixed $T_{\\mathrm{phen}}$, or if the lattice curvature requires an explicit density dependence in $T_{\\mathrm{phen}}$, then the zero-parameter extrapolation to finite density — and the predicted critical end point — would be ruled out. Alternatively, heavy-ion measurements in the 3–10 AGeV beam-energy range that show no non-monotonic behavior in net-proton kurtosis or pion multiplicity across the predicted first-order line would falsify the prediction.","tokens_in":14130,"feed_emoji":"⚛️","tokens_out":9288,"duration_ms":77399,"temperature":0.7,"pith_summary":"This paper argues that adding two ingredients to the SU(3) Polyakov–Nambu–Jona-Lasinio model – the next-to-leading-order term in the large-$N_c$ expansion, which brings in explicit meson degrees of freedom, and a temperature-dependent effective gluon–quark coupling encoded in a phenomenological critical temperature $T_{\\mathrm{phen}}(T)$ – lets the model reproduce the lattice equation of state at zero baryon chemical potential, including pressure, energy density, entropy density, interaction measure, and speed of sound. The same parameter set, extended to finite chemical potential without any new parameters, predicts a first-order chiral phase transition at low temperatures and high baryon density, ending at a critical end point with $T_{\\mathrm{CEP}}=110$ MeV and $\\mu_q=320$ MeV, while approaching perturbative QCD results at very large chemical potentials. If correct, the model supplies an equation of state over the whole $(T,\\mu)$ plane, which is a required input for dynamical simulations of heavy-ion collisions and for neutron star physics.","feed_headline":"A parameter-free extension finds the QCD critical endpoint at 110 MeV","feed_subtitle":"It reproduces lattice thermodynamics at zero density, then predicts a first-order chiral transition and a critical endpoint.","key_machinery":"The load-bearing object is the phenomenological temperature-dependent critical temperature $T_{\\mathrm{phen}}(T)=a+bT+cT^2+dT^3+e/T$ (Eq. 28), which enters the reduced temperature $t_{\\mathrm{phen}}=0.57\\,(T-T_{\\mathrm{phen}}(T))/T_{\\mathrm{phen}}(T)$ in the Polyakov-loop effective potential, thereby modifying the gluon sector beyond a simple constant rescaling of $T_0$. Its coefficients are determined by fitting the lattice equation of state at vanishing chemical potential and are then kept fixed at all $\\mu$; because the quark–gluon back-reaction is not otherwise computed, this function carries the entire coupling between the Polyakov loop and the quarks. The second piece is the next-to-leading-order mesonic grand potential (Eqs. (16)–(23)), which adds pion, kaon, $\\sigma$, and $a_0$ contributions to the pressure at low temperature, using the quark–antiquark phase shift to encode meson masses and widths.","core_discovery":"The central claim is that a modified PNJL model – with the partition function expanded to next-to-leading order in $1/N_c$ so that mesonic (scalar and pseudoscalar) contributions enter the grand potential, and with the Polyakov-loop potential rescaled by a temperature-dependent critical temperature $T_{\\mathrm{phen}}(T)$ fitted to lattice data at $\\mu=0$ – reproduces the lattice thermodynamics of strongly interacting matter at zero chemical potential. Carrying this same parameter set to finite chemical potential, with no additional input, the model predicts a crossover at small $\\mu$ that sharpens into a first-order chiral phase transition for low temperatures, with a critical end point at $T_{\\mathrm{CEP}}=110$ MeV and $\\mu_q=320$ MeV. The pressure computed at very large $\\mu$ agrees with perturbative QCD, and the predicted first-order line lies close to the chemical freeze-out curve extracted from statistical-model fits, suggesting the transition might be probed in heavy-ion collisions at beam energies of order 3–10 AGeV.","pith_inferences":["If the location of the critical end point is robust, the first-order line at low temperature implies a coexistence region in the phase diagram, which could seed density fluctuations and observable signatures in neutron star merger outflows.","The $T_{\\mathrm{phen}}$ parametrization is fitted at $\\mu=0$; computing the model's second-order baryon susceptibilities and comparing with lattice Taylor coefficients beyond $\\kappa$ would test whether the quark–gluon coupling genuinely remains $\\mu$-independent.","Including vector mesons or the di-quark/baryon sector, which the paper leaves for future work, could shift the critical end point; the values 110 MeV and 320 MeV should therefore be read as a baseline estimate rather than a firm prediction."],"forward_implications":["The model provides an equation of state for the full $(T,\\mu)$ plane with no new parameters, ready for use as input to hydrodynamic simulations of ultra-relativistic heavy-ion collisions.","The first-order transition line lies close to the chemical freeze-out curve, so the chiral transition may be observable in heavy-ion collisions at beam energies around 3–10 AGeV.","At very large chemical potentials the pressure agrees with perturbative QCD calculations, supporting the use of the equation of state in the density regime relevant for neutron star interiors and mergers.","Because the speed of sound develops a second minimum below the crossover temperature, the model's predictions for the softest point of the equation of state depend sensitively on the incomplete hadron spectrum and mean-field gap equations.","The dropping Mott temperatures of pions and kaons with increasing chemical potential trace the chiral/deconfinement boundary and can be compared with future lattice data from analytic continuation."],"supporting_citations":[{"why":"Provides the lattice equation of state at zero chemical potential that the modified PNJL model is fitted to reproduce.","marker":"[3]"},{"why":"Motivates the temperature-dependent rescaling of the Polyakov-loop potential and supplies the 0.57 factor.","marker":"[35]"},{"why":"Supplies the next-to-leading-order mesonic grand potential formalism used for the pressure at low temperature.","marker":"[16]"},{"why":"Provides the perturbative QCD pressure at large chemical potential used as a high-density consistency check.","marker":"[14]"},{"why":"Supplies the lattice value of the curvature coefficient $\\kappa$ for comparing the model's Taylor expansion at small chemical potential.","marker":"[40]"},{"why":"Provides the chemical freeze-out curve used to argue that the predicted first-order line may be reachable in heavy-ion experiments.","marker":"[6]"}],"fun_headline_variants":["PNJL model predicts QCD critical end point at 110 MeV","Improved PNJL matches lattice, predicts finite-mu phase transition","Next-to-leading-order PNJL yields critical endpoint and pQCD agreement","Lattice-consistent PNJL: first-order transition at high density","PNJL with mesonic terms pins down critical endpoint"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the temperature-dependent interaction strength fitted at zero density stays the same at all densities. If the true quark–gluon coupling depends on baryon density, the predicted phase boundary and critical end point would be off.","fun_headline_variants_meta":{"raw":{"variants":["PNJL model predicts QCD critical end point at 110 MeV","Improved PNJL matches lattice, predicts finite-mu phase transition","Next-to-leading-order PNJL yields critical endpoint and pQCD agreement","Lattice-consistent PNJL: first-order transition at high density","PNJL with mesonic terms pins down critical endpoint"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00064,"raw_usage":{"total_tokens":2982,"prompt_tokens":1016,"completion_tokens":1966,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":632,"completion_tokens_details":{"reasoning_tokens":1878}},"tokens_in":632,"tokens_out":1966,"duration_ms":13505,"temperature":1.0,"reasoning_tokens":1878,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:49:11.094732+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the model's second-order baryon-number susceptibility and the next Taylor coefficient in $\\mu/T$ and compare with lattice results obtained by analytic continuation or higher-order Taylor expansion. If those coefficients disagree with the model at fixed $T_{\\mathrm{phen}}$, or if the lattice curvature requires an explicit density dependence in $T_{\\mathrm{phen}}$, then the zero-parameter extrapolation to finite density — and the predicted critical end point — would be ruled out. Alternatively, heavy-ion measurements in the 3–10 AGeV beam-energy range that show no non-monotonic behavior in net-proton kurtosis or pion multiplicity across the predicted first-order line would falsify the prediction.","supporting_citations":[{"cited_title":"It tends to the pure Yang Mills reduced temperature for high temperatures, as it should, because asymptotically we expect a plasma of noninteracting quarks and gluons","cited_arxiv_id":null,"evidence_quote":"Motivates the temperature-dependent rescaling of the Polyakov-loop potential and supplies the 0.57 factor."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the next-to-leading-order mesonic grand potential formalism used for the pressure at low temperature."},{"cited_title":"Kurkela and A","cited_arxiv_id":null,"evidence_quote":"Provides the perturbative QCD pressure at large chemical potential used as a high-density consistency check."},{"cited_title":"Borsanyi, G","cited_arxiv_id":null,"evidence_quote":"Supplies the lattice value of the curvature coefficient $\\kappa$ for comparing the model's Taylor expansion at small chemical potential."},{"cited_title":"Cleymans, H","cited_arxiv_id":null,"evidence_quote":"Provides the chemical freeze-out curve used to argue that the predicted first-order line may be reachable in heavy-ion experiments."}],"review_version":1}