{"id":"93cb6119-b3a1-4f13-b05f-e21aaeeda022","arxiv_id":"1908.04678","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Decay-corrected particle yield ratios are converted into light particle and quark chemical potentials, which show a maximum near 4 GeV that the yield ratios themselves do not exhibit.","lead":"The authors apply a standard thermal-model formula to measured particle yield ratios in heavy-ion and proton-proton collisions, after subtracting decay contributions, to extract chemical potentials of pions, kaons, protons, and light quarks. They report that the extracted potentials for protons and light quarks peak near 4 GeV collision energy, a feature not visible in the raw yield ratios themselves.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ~4 GeV maximum in µ_p, µ_u, µ_d, µ_s is inherited from the assumed T_ch(s) curve (Eq. 1) rather than established by the yield ratios; an alternative freeze-out temperature would shift or remove it.","rationale":"The reader's weakest assumption identifies exactly the load-bearing issue: the extracted chemical potentials, and especially the claimed maximum near 4 GeV, depend on the empirical T_ch(s) curve of Eq. (1). My read of the paper confirms this. The formalism in Eqs. (3)-(6) is algebraically transparent and internally consistent, and the authors honestly note the single-T_ch limitation and the absence of a measured pp freeze-out temperature. However, because the headline claim is the existence of a new energy-scale marker at about 4 GeV, and because that marker is generated by multiplying a monotonic T_ch by a monotonic log-ratio, the claim cannot be accepted as a physical feature without testing its sensitivity to the assumed T_ch(s). The paper provides no such test and no error bars. This does not require rejection; the paper's method is clear enough that a straightforward re-analysis with alternative T_ch curves would settle the question. The reader's CONDITIONAL verdict is therefore appropriate, and my stress-test does not change it.","tokens_in":12371,"tokens_out":9090,"duration_ms":96770,"concrete_test":"Recompute µ_p(s), µ_u(s), µ_d(s), and µ_s(s) from the same fitted k_π, k_K, k_p (Eqs. 15-21) using two alternative freeze-out temperature curves in place of Eq. (1): (i) the STAR BES chemical freeze-out temperatures reported in L. Adamczyk et al., Phys. Rev. C 96, 044904 (2017), and (ii) the parametrization from Andronic et al., Nature 561, 321 (2018). Also vary T_lim and the slope 0.45 in Eq. (1) within their published uncertainties. If the maximum remains at 4 ± 0.5 GeV in all cases, the claim is robust; if it moves by more than 1 GeV or disappears, the claimed energy-scale marker is an artifact of Eq. (1).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that µ_p, µ_u, µ_d, and µ_s peak near 4 GeV even though the measured yield ratios are monotonic. This is not a property of the data alone; it follows from Eq. (5), µ_p = -(1/2) T_ch(s) ln k_p(s), together with the externally imposed T_ch(s) of Eq. (1). T_ch(s) rises steeply below about 4 GeV and saturates above, while the fitted ln k_p is essentially a power law (Eq. 21 gives ln k_p ≈ -37.4 s^{-0.884} - 0.007). The product therefore has a maximum where the logarithmic slope of T_ch equals 0.884, which with Eq. (1) happens near s ≈ 4 GeV. The location and even the existence of the claimed maximum are thus controlled by the choice of T_ch(s), not by any extremum in the yield-ratio data. The paper explicitly acknowledges in Sec. 2 that a two- or multi-T_ch scenario is possible and that the single-T_ch treatment is adopted only because Eq. (1) is available; for pp collisions T_ch is not measured and is arbitrarily scaled by factors 0.9 and 0.8. Since no uncertainty is propagated from T_ch or from the fitted ratios into the extracted potentials, the claimed 4 GeV marker is currently an artifact of the assumed parametrization unless shown otherwise.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compiles published mid-rapidity yield ratios π-/π+, K-/K+, and anti-proton/proton from central Au-Au, central Pb-Pb, and pp collisions over √sNN from a few GeV to above 10 TeV, corrects them for strong and weak decay contributions using a recent statistical-thermal-model study [8], and then extracts particle chemical potentials μπ, μK, μp via Eq. (5) and quark chemical potentials μu, μd, μs via Eq. (6). The chemical freeze-out temperature Tch is taken from the empirical parametrization of Eq. (1). The central claim is that, unlike the monotonic yield ratios, the extracted μp, μu, μd, and μs exhibit a maximum at about 4 GeV, which the authors interpret as a possible limiting-fragmentation or phase-transition energy scale.","tokens_in":12790,"tokens_out":5153,"duration_ms":51634,"significance":"If the claimed ~4 GeV maximum in the extracted chemical potentials were robust, it would provide a potentially interesting energy-scale marker in heavy-ion freeze-out phenomenology, and the paper would offer a simple method for simultaneously determining light-particle and light-quark chemical potentials. The paper deserves credit for collecting a wide dataset, for explicitly separating decay-corrected and uncorrected ratios, and for clearly presenting the algebraic relations. However, the central claim is currently not established: the maximum is inherited from the assumed Tch parametrization rather than from the yield-ratio data, no uncertainties are propagated, and the fits that drive the conclusion have poor χ2. The paper's own treatment of pp collisions with Tch scaled by 0.9 and 0.8 shows that the extracted potentials are strongly Tch-dependent. With additional robustness checks, the result could become a useful observation, but as it stands the significance is limited.","major_comments":[{"comment":"The claimed maximum near 4 GeV in μp, μu, μd, and μs is not a property of the yield-ratio data alone, but is essentially imposed by the assumed Tch(s) of Eq. (1). Since μp = -(1/2)Tch ln kp and the fit of Eq. (21) gives ln kp ≈ -37.4 s^{-0.884} - 0.007, which is monotonic, the product Tch(s) s^{-0.884} has a maximum wherever d ln Tch/d ln s = 0.884. For the particular function in Eq. (1) that condition is satisfied near 4 GeV; a different, equally plausible Tch(s) would shift or remove the maximum. The paper itself acknowledges in Sec. 2 that the single-Tch scenario is adopted only because Eq. (1) is available. To support the abstract's central claim, the authors should demonstrate robustness by recomputing the potentials with alternative Tch parametrizations (e.g., from Refs. [3,4,35]) or by varying Tch by ±10%; if the maximum disappears or moves substantially, the claim should be reframed as a property of the assumed Tch model rather than of the data.","section":"Sec. 3, Eq. (5), Eqs. (1) and (21)"},{"comment":"The extracted chemical potentials are presented without any statistical uncertainties. The fit to kp has χ2/dof = 7.7 for Au-Au (Eq. (20)) and 5.3 for Pb-Pb (Eq. (21)), indicating that the fitted curves do not describe the data within the quoted point-to-point scatter; the kK fits also have χ2/dof above 2 (Eqs. (18)–(19)). Because the µp, µu, µd, and µs curves in Figs. 2 and 3 are derived from these fits through Eqs. (5)–(6), the 4 GeV maximum is not shown to be statistically significant relative to fit uncertainty or data scatter. The authors should propagate the fit parameter uncertainties (and, ideally, the experimental uncertainties of the yield ratios) into the extracted potentials and display error bands or error bars in Figs. 2 and 3.","section":"Sec. 3, Eqs. (15)–(21), Figs. 2–3"},{"comment":"The paper's treatment of pp collisions explicitly uses Tch, 0.9Tch, and 0.8Tch in Eq. (5) because the chemical freeze-out temperature is unavailable for pp, producing three sets of extracted potentials (Sec. 3, Fig. 2 caption). This is an explicit demonstration that the results depend sensitively on the assumed Tch value. The same sensitivity should be quantified for Au-Au and Pb-Pb, where Tch is also an external input from Eq. (1) rather than a measured quantity; the conclusion that the maximum at about 4 GeV is a common feature of central heavy-ion collisions cannot be evaluated without such a test.","section":"Sec. 2 and Sec. 3 (pp treatment)"}],"minor_comments":[{"comment":"The title contains a typo: 'negative ly' should be 'negatively'.","section":"Title"},{"comment":"The phrase 'yield rations' appears in Sec. 1 and Sec. 2 and should read 'yield ratios'.","section":"Sec. 1 and Sec. 2"},{"comment":"In Sec. 3, 'absorbtion' should be 'absorption'.","section":"Sec. 3"},{"comment":"The caption sentence describing the solid, dotted, and dashed curves for kπ is grammatically ambiguous: it says 'the solid (dotted) and dashed curves ... in central Au-Au (Pb-Pb) collisions without (with) the corrections of decays and in INEL or NSD pp collisions respectively', which can be read as three systems but only two curve styles. Please clarify which curve corresponds to which system and correction case.","section":"Fig. 1 caption"},{"comment":"The sentence 'According to the functions Eqs. (7)–(13), by using Eqs. (5) and (6)...' refers to the two- and multi-Tch formulas that are not used in the actual analysis; the text should refer to Eqs. (1), (3)–(6) instead to avoid confusion.","section":"Sec. 3, text near Eqs. (15)–(21)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this is an honest, transparent phenomenological paper, but its headline result—the ~4 GeV maximum in mu_p, mu_u, mu_d, mu_s—does not survive scrutiny. The maximum is a product of the empirical T_ch(s) of Eq. (1) and the fitted ln k_p(s); the yield ratios themselves are monotonic. Alternative T_ch curves would shift or remove the peak.\n\nWhat's new and good: they compile published particle yield ratios over a huge energy range (AGS to LHC) and apply the decay corrections from Yu & Luo to get \"primary\" ratios. The algebra relating ratios to particle and quark chemical potentials is simple, standard, and correctly laid out. They show their fits, give chi2 values, and explicitly acknowledge that a multi-T_ch scenario is possible—they use single T_ch simply because Eq. (1) is available. That is honest.\n\nSoft spots, in order of severity. First, the central claim. From Eq. (5), mu_p = -(1/2) T_ch ln k_p. With their fit, ln k_p ≈ -37.4 s^{-0.884} - 0.007, and T_ch from Eq. (1) rises steeply below ~4 GeV then saturates. The product necessarily peaks near where the logarithmic slopes match. The peak is not a feature of the data; it is a feature of the chosen T_ch curve. They do not test other T_ch parametrizations or propagate any uncertainty from T_ch. For pp collisions they simply scale T_ch by 0.9 and 0.8, which changes magnitudes but not the logic. Second, no error bars appear on any extracted chemical potential, even though the input ratios have uncertainties and the k_p fits have chi2/dof of 7.7 (Au-Au) and 5.3 (Pb-Pb). Third, the physical interpretation of the peak—limiting fragmentation, phase transition—is speculative and not tied to any independent observable. These are fixable, but as written the load-bearing claim is conditional.\n\nWho this is for: heavy-ion phenomenologists working on chemical freeze-out systematics. The compilation of decay-corrected ratios and the transparent extraction procedure could be useful. I would not cite the 4 GeV peak as a result, but I might cite the compilation if the authors add a table and uncertainties.\n\nRecommendation: this deserves a serious referee, not a desk reject. The method is sound and the data compilation is valuable; the authors need to add uncertainty propagation, test alternative T_ch(s), and either soften or support the 4 GeV claim. With those changes it could be a serviceable contribution.\n\nBest","headline":"Useful, transparent compilation, but the 4 GeV peak in extracted quark chemical potentials is inherited from an assumed T_ch(s) curve and is not established by the data.","tokens_in":13294,"tokens_out":2708,"would_cite":false,"duration_ms":27431,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.40.Aq","14.65.Bt","25.75.-q"],"model":"deepseek-v4-flash","headline":"The paper claims that after removing strong and weak decay backgrounds, the extracted chemical potentials of protons and up, down, and strange quarks peak near 4 GeV even though the raw yield ratios show no such extremum.","keywords":["light particle chemical potentials","light quark chemical potentials","yield ratios of negatively to positively charged particles","chemical freeze-out temperature","strong and weak decay corrections","heavy-ion collisions","proton-proton collisions","energy dependence"],"falsifier":"Recompute the proton and quark chemical potentials from the same decay-corrected ratios using a different accepted freeze-out temperature parameterization, such as separate temperatures for strange and non-strange particles; if the maximum near 4 GeV vanishes under that alternative, the peak is an artifact of the single-temperature curve rather than a physical energy scale.","tokens_in":12184,"feed_emoji":"⚛️","tokens_out":9663,"duration_ms":91916,"temperature":0.7,"pith_summary":"This paper tries to show that an energy scale sits near 4 GeV in the particle chemistry of high-energy collisions. From measured ratios of negatively to positively charged pions, kaons, and protons, with strong-decay and weak-decay backgrounds removed, the authors extract chemical potentials for those particles and for up, down, and strange quarks across collision energies from a few GeV to above 10 TeV. Their central result is that the proton, up-quark, down-quark, and strange-quark chemical potentials rise to a maximum around 4 GeV and then fall, while the raw yield ratios show no such extremum. If this is right, 4 GeV becomes a new marker for where the hot, dense collision zone changes character.","feed_headline":"Quark chemical potentials peak near 4 GeV after decay correction","feed_subtitle":"Removing decay backgrounds makes proton and quark chemical potentials peak, a possible new energy-scale marker.","key_machinery":"The machinery is a set of Boltzmann-statistics relations between yield ratios and chemical potentials. For any particle $j$, $k_j = \\exp(-2\\mu_j/T_{\\rm ch})$, so $\\mu_j = -\\frac12 T_{\\rm ch}\\ln k_j$; for light quarks, $\\mu_u = -\\frac16 T_{\\rm ch}(\\ln k_\\pi+\\ln k_p)$, $\\mu_d = -\\frac16 T_{\\rm ch}(-2\\ln k_\\pi+\\ln k_p)$, and $\\mu_s = -\\frac16 T_{\\rm ch}(\\ln k_\\pi-3\\ln k_K+\\ln k_p)$. The temperature $T_{\\rm ch}(\\sqrt{s_{NN}})$ is an empirical limiting-temperature curve from the literature, and the measured ratios are replaced by primary-production ratios taken from a published model calculation that removes strong and weak decays. Fitting smooth curves to the energy-dependent ratios and then applying these identities is what turns monotonic ratio trends into chemical potentials with a peak near 4 GeV.","core_discovery":"After correcting the measured yield ratios $k_\\pi=\\pi^-/\\pi^+$, $k_K=K^-/K^+$, and $k_p=\\bar p/p$ by removing strong-decay contributions from high-mass resonances and weak-decay contributions from heavy-flavor hadrons, the paper extracts chemical potentials from $k_j = \\exp(-2\\mu_j/T_{\\rm ch})$, giving $\\mu_j = -\\tfrac12 T_{\\rm ch}\\ln k_j$ for particles and linear combinations of the $\\ln k_j$ for quarks. Over $\\sqrt{s_{NN}}$ from a few GeV to above 10 TeV, the energy-dependent $\\mu_p$, $\\mu_u$, $\\mu_d$, and $\\mu_s$ show a maximum at about 4 GeV, whereas $k_\\pi$, $k_K$, $k_p$, $\\mu_\\pi$, and $\\mu_K$ do not show such an extremum. The paper ties this scale to the onset of hadronic limiting fragmentation and to a possible liquid-to-gas-like transition of nucleons and mesons in central nucleus-nucleus collisions.","pith_inferences":["Testing the same decay-corrected ratios with a different published freeze-out temperature curve would show whether the 4 GeV peak survives; if it does not, the peak is an artifact of the assumed $T_{\\rm ch}$.","A direct measurement of the freeze-out temperature in proton-proton collisions around 4 GeV, rather than an assumed fraction of $T_{\\rm ch}$, would settle whether the peak in pp data is physical.","The same formalism extends to D and B meson yield ratios, so charm and bottom quark chemical potentials could be extracted once those ratios are measured with decay corrections.","Because the decay corrections come from a single published calculation, an independent calculation of primary-production ratios would test whether the corrections, not the data, create the peak."],"forward_implications":["The 4 GeV maximum gives a specific energy to look for in beam-energy scans: particle and quark chemical potentials should change slope there even though the raw yield ratios do not.","The extraction supplies a way to estimate up, down, and strange quark chemical potentials directly from measured particle ratios, not just the baryon chemical potential.","At collision energies above the top RHIC energy, all extracted chemical potentials tend to zero and all ratios tend to one, so the method's predictions at the LHC are essentially vanishing potentials.","In proton-proton collisions, where the freeze-out temperature is not measured, using $T_{\\rm ch}$, $0.9T_{\\rm ch}$, or $0.8T_{\\rm ch}$ changes the size of the potentials but preserves the reported energy trends."],"supporting_citations":[{"why":"Supplies the model relations connecting yields, temperature, and chemical potentials.","marker":"[1–4]"},{"why":"Gives the empirical chemical freeze-out temperature $T_{\\rm ch}(\\sqrt{s_{NN}})$ used to convert every ratio into a potential.","marker":"[4, 5, 33–35]"},{"why":"Provides the exponential relation between antiproton-to-proton ratio, $T_{\\rm ch}$, and proton chemical potential.","marker":"[17, 36, 37]"},{"why":"Supplies the quark-composition expressions that let the paper extract up, down, and strange quark potentials from particle ratios.","marker":"[38, 39]"},{"why":"Provides the primary-production yield ratios with strong and weak decay contributions removed.","marker":"[8]"}],"fun_headline_variants":["Decay-corrected quark chemical potentials peak at 4 GeV","Proton and quark chemical potentials peak after decay removal","4 GeV peak in chemical potentials after removing decays","Yield ratios don't peak, but corrected chemical potentials do","Chemical potentials reveal 4 GeV scale after decay correction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole extraction rests on the empirical freeze-out temperature curve, which is assumed to hold for every collision system and energy in the study; if that curve is wrong or does not apply to proton-proton collisions, the 4 GeV maximum could move or disappear.","fun_headline_variants_meta":{"raw":{"variants":["Decay-corrected quark chemical potentials peak at 4 GeV","Proton and quark chemical potentials peak after decay removal","4 GeV peak in chemical potentials after removing decays","Yield ratios don't peak, but corrected chemical potentials do","Chemical potentials reveal 4 GeV scale after decay correction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000619,"raw_usage":{"total_tokens":2886,"prompt_tokens":974,"completion_tokens":1912,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":1833}},"tokens_in":590,"tokens_out":1912,"duration_ms":14782,"temperature":1.0,"reasoning_tokens":1833,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:59:11.215662+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the proton and quark chemical potentials from the same decay-corrected ratios using a different accepted freeze-out temperature parameterization, such as separate temperatures for strange and non-strange particles; if the maximum near 4 GeV vanishes under that alternative, the peak is an artifact of the single-temperature curve rather than a physical energy scale.","supporting_citations":[{"cited_title":"Yu and X","cited_arxiv_id":null,"evidence_quote":"Provides the primary-production yield ratios with strong and weak decay contributions removed."}],"review_version":1}