{"id":"12cf3d3a-93cd-4cb8-8fc1-66e1e6aedbe1","arxiv_id":"1908.05939","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In the Polyakov linear-sigma model, optimized perturbation theory agrees with lattice QCD better than the mean-field approximation for high-order cumulants of baryon, charge, and strangeness fluctuations.","lead":"This paper compares two approximation schemes for the Polyakov linear-sigma model, an effective theory of QCD matter, and checks them against lattice QCD results. It reports that the optimized perturbation theory (OPT) matches lattice data better than the simpler mean-field approximation, especially for higher-order fluctuations of conserved charges.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Polyakov potential parameters a and b in Eq. (5) are never specified, and U is not optimized; the claimed OPT improvement in chi6_B/chi2_B and chi8_B may depend on this unstated choice.","rationale":"The paper is a genuine attempt to compare two approximation schemes in the PLSM, and the qualitative direction of the claim is plausible: OPT modifies the quark sector, and higher-order cumulants are more sensitive to that sector. The comparison to lattice data is presented honestly, with several orders of moments shown. However, the most load-bearing condition for the central claim is the control of the Polyakov-loop potential. The manuscript leaves the parameters a and b in U(phi,phi*,T) unspecified, and explicitly does not optimize U in the OPT procedure. Since the higher-order cumulants that support the headline claim are obtained from the full free energy, the unstated parameters can affect the absolute magnitude and temperature dependence of these cumulants. Without specifying a and b, the results are not reproducible, and the claimed improvement of OPT over MFA for chi6_B/chi2_B and chi8_B is not shown to be robust against the allowed range of these parameters. The reader's conditional verdict is appropriate because the missing information and the lack of a parameter-robustness check are exactly the kind of conditions that should be satisfied before the claim is accepted. My read does not change the verdict from CONDITIONAL, so UNCHANGED is chosen. I partially agree with the reader's weakest_assumption: the reader also cited the reliability of the truncated delta-expansion, which is a separate concern not addressed here; my focus is the unoptimized, unspecified Polyakov potential, which is directly traceable in the text.","tokens_in":19026,"tokens_out":8062,"duration_ms":75962,"concrete_test":"Obtain the values of a and b used to generate Figs. 6-8 (from the authors, or from the PLSM parameter sets in Refs. [20] or [31]). Recompute chi6_B/chi2_B and chi8_B at mu_B=0 for MFA and OPT with a and b varied by ±20% around those values, or over the range spanned by previously published PLSM fits. If the OPT curve at T/Tchi > 1.2 stays closer to the lattice data than MFA for all tested (a,b), the claim is robust; if the ordering reverses for any plausible (a,b), the headline claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim rests on the higher-order cumulants chi6_B/chi2_B and chi8_B shown in Fig. 8. These are derivatives of the full free energy Eq. (9), which contains the Polyakov-loop potential U(phi,phi*,T) of Eq. (5). The parameters a and b in that potential are nowhere assigned numerical values in the manuscript; they are not listed in the parameter set determined at m_sigma = 800 MeV, and no reference gives them explicitly. Section II explicitly states that U is \"not directly impacted by the OPT approach\", so the same U is used in the MFA and OPT calculations. Consequently, every difference between OPT and MFA in the cumulant ratios is generated only by the quark/meson sector, and the absolute scale of the cumulants—and hence the visual agreement with lattice—depends on the unspecified a and b. If the authors' curves were produced with particular values of a and b taken from a previous fit, the central comparison is not reproducible from the manuscript, and the conclusion that OPT becomes \"more closer to QCD\" at higher order is not tested against the uncertainty in those parameters. The paper itself notes (Sec. III C) that the pseudo-critical temperature depends on \"the input parameters, and the Polyakov potential\", which is exactly the uncontrolled dependence that could affect the higher-order moments.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript compares two approximation schemes for the Polyakov linear-sigma model (PLSM): the mean-field approximation (MFA) and optimized perturbation theory (OPT), a δ-expansion with the variational mass parameter fixed by the principle of minimal sensitivity. The authors compute chiral condensates, Polyakov-loop expectation values, pseudo-critical temperatures, the thermodynamic pressure, and second- and higher-order cumulants of conserved charges, and compare these with lattice QCD results. The central claim is that, as one moves from lower- to higher-order moments, OPT agrees with lattice QCD better than MFA, with the improvement highlighted for χ6_B/χ2_B and χ8_B. The analysis is performed at finite temperature and finite quark chemical potential, and the paper also maps out the QCD phase diagram in both approximations.","tokens_in":19324,"tokens_out":2396,"duration_ms":25386,"significance":"If the central claim is supported, the paper would provide a useful demonstration that a truncated δ-expansion with PMS can outperform the mean-field treatment for higher-order fluctuations in a chiral effective model with Polyakov-loop dynamics. The comparison uses external lattice data that are not used to fit the OPT parameters, so the agreement is not circular in the strict sense. The paper also contains explicit expressions for the OPT free energy and the gap equations, which can be adapted by other practitioners. However, the significance is currently limited by the lack of a quantitative closeness metric, by the absence of specified values for the Polyakov-loop potential parameters a and b in Eq. (5), and by an internal tension between the pressure results (where OPT is farther from the Stefan-Boltzmann limit than MFA) and the claimed higher-order improvement. These issues are addressable, but they must be resolved before the headline conclusion can be accepted.","major_comments":[{"comment":"The parameters a and b of the Polyakov-loop potential U(φ,φ*,T) are nowhere given numerical values, and they are not listed among the parameters fixed at m_σ = 800 MeV. Since U enters the free energy in Eq. (9) and therefore all cumulant ratios, the absolute scale of χ6_B/χ2_B and χ8_B, and hence their visual agreement with lattice data, depends on these unstated inputs. The manuscript cannot be reproduced from the information provided, and no sensitivity analysis is given for the dependence of the higher-order moments on a and b. This is load-bearing for the central claim that OPT becomes closer to QCD at high order.","section":"Section II, Eq. (5)"},{"comment":"The headline conclusion that OPT agrees 'excellently' with lattice QCD for χ6_B/χ2_B and χ8_B is based entirely on visual inspection of Fig. 8. There is no quantitative measure of closeness, no error estimate or uncertainty band on the model curves, and no statement of the temperature range over which the agreement is claimed. Given that the lattice points carry sizable error bars and that the two model curves are close to each other in several regions, a numerical comparison (for example, a weighted squared deviation or a chi-square per degree of freedom over a stated T/Tχ interval) is needed to support the claim that OPT is actually better than MFA.","section":"Section III D 2, Fig. 8"},{"comment":"There is an internal consistency issue between the pressure results and the overall conclusion. In Section III C the authors state that at T ≤ 2.5Tχ the gap to the Stefan-Boltzmann limit is 31.8% for MFA and 34.9% for OPT, meaning OPT is farther from the ideal-gas limit for the pressure. Yet Section IV concludes that 'when moving to lower- to higher-order moments, the OPT approach becomes more and more reliable than the MFA.' The paper does not explain why a worse low-order quantity (pressure) is consistent with better high-order cumulants, or whether the pressure discrepancy is a signature of the truncation or of the fixed Polyakov-loop potential. This needs to be addressed explicitly, since the central narrative depends on the higher-order improvement being meaningful rather than accidental.","section":"Section III C and Section IV"}],"minor_comments":[{"comment":"The phrase 'OPT becomes more closer to QCD' is ungrammatical; it should read 'OPT becomes closer to QCD' or 'OPT moves closer to QCD'.","section":"Abstract and Section III D 2"},{"comment":"The label 'µf=0.150 MeV' appears to be a typo; it should presumably read 'µf=150.0 MeV' to match the other panels.","section":"Figure 2, panel (d)"},{"comment":"The caption repeats the label '(a)' and the lattice reference twice, and the embedded text in the figure image duplicates the MFA and OPT labels; the figure should be cleaned up.","section":"Figure 3, caption"},{"comment":"The sentence stating that the second term in Eq. (11) 'is approximately equal the derivative of the first term' is missing a word ('to') and would benefit from a more precise mathematical statement of the approximate relation.","section":"Section II, Eq. (11)"},{"comment":"Reference [21], 'A. N. Tawfik, A. M. Diab, M. T. Ghoneim, and H. Anwer, (2019)', lacks a journal or preprint identifier and cannot be located by the reader; a complete citation should be provided.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper's central comparison is not circular because the lattice data are not used to fit the OPT parameters, but the missing numerical values of a and b in Eq. (5) make the central comparison non-reproducible. The lack of a quantitative metric for 'closer to QCD' is a further obstacle. These issues are fixable within the manuscript's scope, so major revision is appropriate. I would also note that the paper relies heavily on visual comparison with lattice data from several different references with different normalizations, which weakens the force of the 'excellent agreement' statements."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the useful thing: this is the first direct comparison I know of between OPT and MFA within the Polyakov linear-sigma model focused on higher-order cumulants. The OPT setup (delta-expansion, PMS) is standard and clearly laid out. The paper also does an honest job of showing where OPT and MFA differ, e.g. Fig. 5d, and the phase diagram comparison with lattice is coherent.\n\nBut the central claim is not established. The headline says OPT becomes 'more closer' to QCD at higher moments, yet the paper's own Fig. 8 text says both approximations reproduce the lattice within its large error bars. For chi6_B/chi2_B and chi8_B, the MFA and OPT curves are close to each other and both are consistent with lattice at T/Tchi > 1.2. The conclusion that OPT is better seems to rest on eyeballing which dashed curve 'fits well'. There is no quantitative closeness metric, no error bars on the model curves, and no sensitivity analysis.\n\nThe more serious problem is reproducibility. The Polyakov-loop potential Eq. (5) has parameters a and b, and the Yukawa coupling g appears in the Lagrangian, but none of these are given a numerical value anywhere in the paper. The parameter list at m_sigma = 800 MeV covers the LSM couplings only. Since U is the same in both approximations, the relative OPT-vs-MFA comparison is not affected by a and b, but the absolute magnitude of the cumulants and the visual agreement with lattice certainly are. The paper even admits (Sec. III C) that the pseudo-critical temperature depends on the Polyakov potential. That is exactly the uncontrolled dependence that could affect the higher moments. If a and b come from a prior lattice fit, the paper needs to say so and test the sensitivity.\n\nOne more point: the pressure comparison works against the headline. MFA is closer to the SB limit (31.8% gap vs 34.9%), so the claim that OPT is globally closer to QCD is at best restricted to selected observables.\n\nOverall: a plausible but not demonstrated methodological recommendation. The comparison is legitimate, and the underlying calculation is probably correct, but the manuscript as written is not reproducible and the key conclusion overreaches the evidence. I would not desk-reject it; a serious referee could ask for the missing parameters, a quantitative agreement measure, and a more careful statement of what OPT actually improves. That could be a useful paper after revision.","headline":"A legitimate but under-supported comparison of OPT vs MFA in PLSM for high-order cumulants, undermined by unspecified Polyakov potential parameters and purely visual agreement claims.","tokens_in":19889,"tokens_out":3746,"would_cite":false,"duration_ms":35274,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.30.Rd","11.10.Wx","12.39.Fe","02.60.-x"],"model":"deepseek-v4-flash","headline":"In the Polyakov linear-sigma model, optimized perturbation theory follows lattice QCD more closely than the mean-field approximation does, and the gap grows with the order of the conserved-charge moments.","keywords":["chiral symmetries","chiral transition","chiral Lagrangian","numerical approximation and analysis","optimized perturbation theory","mean-field approximation","Polyakov linear-sigma model","conserved-charge cumulants"],"falsifier":"Tighter lattice QCD data on $\\chi_8^B$ in the range $T/T_\\chi \\approx 1.2$ to $2.0$—the current comparison is made within large error bars—would settle whether OPT really sits closer to the central values than MFA does. A second check would be to recompute the high-order cumulants with the gluon-sector potential also included in the optimization; if the OPT advantage disappears or reverses, the paper's conclusion would not survive.","tokens_in":18875,"feed_emoji":"⚛️","tokens_out":18644,"duration_ms":144692,"temperature":0.7,"pith_summary":"This paper asks whether optimized perturbation theory (OPT), a variationally tuned resummation scheme, improves on the standard mean-field approximation (MFA) when both are applied to the Polyakov linear-$\\sigma$ model, a low-energy effective description of the quark-hadron transition. It establishes that the two schemes give nearly identical chiral condensates, subtracted condensates, and pseudo-critical temperatures at small chemical potential, but that OPT is closer to lattice QCD for the thermodynamic pressure and for fluctuations and correlations of conserved charges. The gap widens with the order of the cumulant: the ratios $\\chi_6^B/\\chi_2^B$ and $\\chi_8^B$ computed in OPT agree with lattice QCD above about $1.2\\,T_\\chi$, where the mean-field curves drift away. Because these high-order moments are the proposed experimental signatures of the deconfinement crossover, the result matters for interpreting heavy-ion collision data with an affordable effective model.","feed_headline":"High-order moments put OPT ahead of mean field","feed_subtitle":"These cumulants are the proposed experimental signatures of the quark–gluon transition","key_machinery":"The engine is the $\\delta$-expansion of the PLSM Lagrangian. The full Lagrangian is rewritten as $L = L_0(\\eta) + \\delta \\,[L - L_0(\\eta)]$, where $L_0(\\eta)$ is a solvable free Lagrangian containing an arbitrary mass parameter $\\eta$; results are evaluated at $\\delta = 1$ after truncating at order $\\delta^2$. The variational parameter $\\eta$ is not chosen by hand but fixed by the principle of minimal sensitivity, $\\partial F_{\\mathrm{OPT}}/\\partial \\eta = 0$ at $\\delta = 1$, which makes the final answer approximately independent of the artificial mass. This $\\eta$ enters the quark dispersion relation $E_f = [|\\mathbf{P}|^2 + (m_f + \\eta)^2]^{1/2}$ and produces an extra term in the quark-antiquark thermodynamic potential that is absent in the mean-field limit, recovered by taking $\\eta = 0$ at $\\delta = 1$. The Polyakov-loop potential $U(\\varphi,\\varphi^*,T)$ is left outside the optimization and imported unchanged from a prior fit.","core_discovery":"The central claim is that, within the Polyakov linear-$\\sigma$ model, the optimized perturbation theory—a $\\delta$-expansion of the Lagrangian with an arbitrary mass parameter $\\eta$ fixed by the principle of minimal sensitivity—provides a systematically better account of the QCD equation of state and of conserved-charge fluctuations than the mean-field limit of the same model. The chiral order parameters and the crossover temperature are largely unaffected by the choice of approximation, but the story changes for derivatives of the free energy with respect to chemical potentials. For the quadratic susceptibilities of baryon number, electric charge, and strangeness, and for their mutual correlations, OPT sits closer to the lattice data than MFA does. The most distinctive result is at high order: $\\chi_6^B/\\chi_2^B$ and $\\chi_8^B$, computed in OPT, match the lattice QCD bands at $T/T_\\chi \\gtrsim 1.2$, a region where the MFA curves are visibly less accurate. The paper's conclusion is that moving from lower- to higher-order moments makes the OPT approach progressively more reliable than MFA.","pith_inferences":["Not pursued in the paper: applying the same $\\delta$-expansion to the gluon-sector Polyakov-loop potential rather than importing it from a prior fit; if the remaining gap between OPT and lattice QCD shrinks or grows, it would reveal how much of the agreement is carried by the unoptimized gluon part.","If the paper's ordering is correct—OPT improves with cumulant order—the method should be tested on even higher moments such as $\\chi_{10}^B$ and on strangeness-related high-order cumulants, where lattice data with small errors could separate the schemes more decisively.","An interpretation the paper does not make explicit is that the variational mass $\\eta$ resums quark thermal self-energy contributions that the mean-field approximation drops, which would explain why the improvements show up in chemical-potential derivatives of the free energy rather than in the order parameters themselves."],"forward_implications":["Within the Polyakov linear-sigma model at finite temperature and chemical potential, OPT is the approximation to use for the pressure and the cumulants of conserved charges, while MFA remains adequate for the order parameters and for the crossover line at small $\\mu_B$.","The ratios $\\chi_6^B/\\chi_2^B$ and $\\chi_8^B$ become practical discriminators between approximation schemes: OPT places them inside the lattice-QCD bands for $T/T_\\chi \\gtrsim 1.2$, where MFA drifts away.","At large baryon chemical potential the two schemes disagree on the pseudo-critical temperature, with OPT giving a slightly higher $T_\\chi$, so estimates of the phase boundary in that regime should carry the difference as a systematic uncertainty.","Because OPT reproduces the off-diagonal correlations $\\chi_{11}^{QS}$ and $-\\chi_{11}^{BS}$ better than MFA, these correlations become reliable tests of whether the model is in the right regime."],"supporting_citations":[{"why":"It introduces the optimized perturbation theory, the delta-expansion that the paper generalizes to the Polyakov linear-sigma model.","marker":"[23]"},{"why":"It formulates the principle of minimal sensitivity used to fix the variational mass parameter.","marker":"[28]"},{"why":"It supplies the OPT relation between the variational mass and the sigma-field expectation values.","marker":"[30]"},{"why":"It fixes the PLSM parameters used in both approximations.","marker":"[27]"},{"why":"It provides the lattice subtracted-condensate data used to validate the calculation in Fig. 3.","marker":"[43]"},{"why":"It provides the lattice pressure data that both approximations are compared against in Fig. 5.","marker":"[48]"},{"why":"It supplies the lattice second-order susceptibility and correlation data used in Fig. 6.","marker":"[52]"},{"why":"It provides the lattice higher-order susceptibility data, including the sixth- and eighth-order baryon results.","marker":"[53]"},{"why":"It offers the lattice results for the sixth- and eighth-order baryon susceptibilities that support the high-temperature agreement.","marker":"[61]"}],"fun_headline_variants":["OPT wins on high-order charge moments","Sixth and eighth cumulants favor OPT","OPT beats mean field on QCD moments","High-order moments side with OPT"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's conclusion rests on trusting the truncated optimization procedure, with its adjustable mass fixed by a minimal-sensitivity condition, to be a faithful approximation of the model's free energy, while the gluon-sector potential is taken from an earlier fit and left unoptimized; if either piece misrepresents the theory, the claimed improvement of OPT over MFA could be an artifact.","fun_headline_variants_meta":{"raw":{"variants":["OPT wins on high-order charge moments","Sixth and eighth cumulants favor OPT","OPT beats mean field on QCD moments","High-order moments side with OPT"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000742,"raw_usage":{"total_tokens":3268,"prompt_tokens":859,"completion_tokens":2409,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":475,"completion_tokens_details":{"reasoning_tokens":2357}},"tokens_in":475,"tokens_out":2409,"duration_ms":17574,"temperature":1.0,"reasoning_tokens":2357,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:59:58.139172+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Tighter lattice QCD data on $\\chi_8^B$ in the range $T/T_\\chi \\approx 1.2$ to $2.0$—the current comparison is made within large error bars—would settle whether OPT really sits closer to the central values than MFA does. A second check would be to recompute the high-order cumulants with the gluon-sector potential also included in the optimization; if the OPT advantage disappears or reverses, the paper's conclusion would not survive.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It introduces the optimized perturbation theory, the delta-expansion that the paper generalizes to the Polyakov linear-sigma model."},{"cited_title":"Kneur, M","cited_arxiv_id":null,"evidence_quote":"It supplies the OPT relation between the variational mass and the sigma-field expectation values."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It fixes the PLSM parameters used in both approximations."},{"cited_title":"Hansen et al., Phys","cited_arxiv_id":null,"evidence_quote":"It provides the lattice subtracted-condensate data used to validate the calculation in Fig. 3."},{"cited_title":"Borsanyi et al., Nature 539, 69 (2016)","cited_arxiv_id":null,"evidence_quote":"It provides the lattice pressure data that both approximations are compared against in Fig. 5."},{"cited_title":"The left panels of Fig","cited_arxiv_id":null,"evidence_quote":"It supplies the lattice second-order susceptibility and correlation data used in Fig. 6."},{"cited_title":"Andronic, P","cited_arxiv_id":null,"evidence_quote":"It provides the lattice higher-order susceptibility data, including the sixth- and eighth-order baryon results."}],"review_version":1}