{"id":"fa80e20f-6eaa-4d3d-b679-d7328c577044","arxiv_id":"hep-ph/9903292","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A set of thermodynamic observables and analytic tools that turn event-by-event pion fluctuations into a search strategy for the QCD critical endpoint, with order-of-magnitude estimates that the critical signal should be observable.","lead":"This paper develops the theoretical toolkit for using event-by-event fluctuations of pion multiplicity and mean transverse momentum to locate the QCD critical endpoint in heavy ion collisions. It is the long companion to the authors' 1998 letter and became the reference framework for the experimental Beam Energy Scan program.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The 'easily observable' conclusion rests on the adopted xi_max ~ 6 fm, but the paper's own finite-time scaling (Eq. 112) gives ~2.5 fm; since the signal scales as xi^2, this is the load-bearing quantitative uncertainty, and G is a co-equal factor.","rationale":"The reader's weakest_assumption correctly identifies the finite-time estimate of xi_max as the load-bearing quantitative uncertainty for the 'easily observable' claim, and the paper itself calls this estimate 'suspect for several reasons.' I agree with the verdict of ACCEPT: the analytic core (Eqs. 50 and 88, the finite-multiplicity estimator of Eq. 145) is sound and checked against known limits, the NA49 comparison is framed as broad agreement rather than a precision test, and the qualitative nonmonotonic prediction survives conservative inputs. My sharpening is that the paper's own Eq. (112) yields xi <~ 2.5 fm, not the adopted 6 fm, and the jump is justified only by an unquantified compensating effect; combined with the quadratic dependence on G (itself a tricritical extrapolation), the raw signal at conservative inputs is comparable to the systematic background rather than 50x the statistical error. This does not overturn the central claim, because the nonmonotonic variation with control parameters is the actual detection strategy and systematics are roughly energy-independent, but the abstract's 'easily observable' phrasing is stronger than the body supports in the most conservative case. The concern is acknowledged by the authors and does not rise to a correctness error, so the verdict stands unchanged. The preliminary, private NA49 data (G. Roland, private communication) limit independent reproduction of the data comparison, but that is a data-availability caveat, not a flaw in the argument.","tokens_in":49431,"tokens_out":6201,"duration_ms":111850,"concrete_test":"Recompute the predicted sqrt(F_sigma) - 1 from Eq. (115) using xi = 2.5 fm (the paper's own Eq. 112 finite-time bound) and G = 300 MeV with mu_pi = 0, including the factor-of-two dilution from non-direct pions, and compare it to the combined statistical plus systematic uncertainty on sqrt(F) at a single collision energy. If the single-energy signal-to-background ratio is below unity, the 'easily observable' claim in abstract:s9 rests entirely on the nonmonotonic control-parameter scan, and the paper should say so explicitly rather than quoting the 50x statistical-error figure. A dynamical (hydro + Langevin) estimate of the actual xi reached at freeze-out would settle the same question more physically.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is Eq. (115): F_sigma - 1 = 0.35 (G/300 MeV)^2 (xi/6 fm)^2 (mu_pi = 60 MeV), with the multiplicity analogue in Eq. (117). The qualitative nonmonotonic prediction is robust, but the 'easily observable' language of abstract:s9 depends on the product G*xi, and both factors are uncertain by factors of a few. The reader correctly flags xi. I would sharpen it: there is an internal tension the paper does not resolve. Eq. (112), the finite-time scaling bound, gives xi/l < (tau_H/l)^{1/Delta_t} ~ 2.5, i.e. xi <~ 2.5 fm. Eq. (110), the finite-size bound, gives xi <~ 12 fm. The binding constraint is therefore 2.5 fm, yet the paper adopts 6 fm, justified only qualitatively by an unquantified 'compensating good effect' (longer dwell time near T_c from the large specific heat) and an unknown prefactor. Taking Eq. (112) at face value, F_sigma - 1 falls from 0.35 to ~0.06 (mu_pi = 60) or from 0.14 to ~0.024 (mu_pi = 0). After the factor-of-two dilution from non-direct pions, sqrt(F_sigma) - 1 is ~0.6-1.5%, comparable to the systematic background the paper itself estimates (Bose ~2%, energy-conservation/thermal-contact ~1%), not the 50x statistical error quoted. The nonmonotonic scan can still separate signal from a roughly constant background, but that argument is not made quantitatively. Separately, G enters quadratically and its value at freeze-out comes from extrapolating tricritical scaling (G ~ m^{3/5}) from m = 0 to physical quark mass, which the authors themselves call uncertain; a G of 150 MeV instead of 300 MeV cuts the signal fourfold. The paper is transparent about both, so this is a quantitative-reach concern, not an error in the derivations of Eqs. (50) and (88), which check correctly against their stated limits.","agreement_with_reader":"agree"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The paper argues that event-by-event fluctuations of charged-pion multiplicity and mean transverse momentum in heavy ion collisions can be used to locate the second-order critical endpoint E of the QCD first-order transition, because signatures characteristic of freeze-out near E rise and then fall nonmonotonically as a control parameter (collision energy, centrality, ion size) is varied. The authors first build the noncritical baseline: an ideal Bose gas analysis (§2) establishing the microscopic correlator ⟨Δn_p Δn_k⟩ and the relation between event-by-event and inclusive variances, then a resonance-gas simulation (§3) including Bose enhancement, resonance decays and radial flow, compared with preliminary NA49 PbPb data at 160 AGeV. They then treat two critical effects: thermal contact between the observed pions and a heat bath with possibly divergent specific heat (§4), which produces an off-diagonal anti-correlation from energy conservation that is switched off near E; and the direct Gσππ coupling of pions to critical sigma fluctuations (§5), which adds a term ∝ 1/m_σ² to the correlator and yields Eqs. (115)–(118) for the enhancement of p_T, multiplicity and cross-correlation fluctuations. §6 treats the delayed decays of light sigmas after freeze-out, predicting a low-p_T pion excess at ⟨p_T⟩ ~ m_π/2 with enhanced multiplicity fluctuations. An Appendix gives finite-sample estimators, including a multiplicity-weighted estimator for v_ebe (Eq. 145) that removes 1/N bias.","tokens_in":49859,"tokens_out":6494,"duration_ms":118143,"significance":"If the central claim holds, this is a high-impact paper: it converts a qualitative suggestion into a quantitative experimental program with a specific, falsifiable signature — a nonmonotonic bump in √F and in ⟨(ΔN)²⟩ whose maximum must occur at the same control-parameter value across several observables, plus a low-p_T pion population at ⟨p_T⟩ ~ m_π/2 that appears and disappears. The derivations are checkable and the authors check them: Eq. (50) reduces to the canonical result as C_A → ∞ and enforces fixed energy as C_A → 0, and Eq. (88) is cross-checked against the mass-shift interpretation of Eq. (81). The Appendix estimator (Eq. 145) is a concrete, reusable methodological contribution that removes a systematic 1/N bias from published v_ebe values, and the correction (147) is applied transparently to the NA49 numbers. The authors are also candid about their limitations — the finite-time estimate of §5.4 is explicitly called 'suspect for several reasons', the prefactor ambiguity in G_E is stated in §5.3, the lowest-order-in-G treatment is flagged in §5, and the blue-shift approximation is flagged in §3.3 — which makes the paper easy to assess and to build on.","major_comments":[{"comment":"The abstract's conclusion that the effects 'should be easily observable' rests on the adopted ξ_max ≈ 6 fm, but this exceeds the paper's own finite-time scaling estimate by a factor of about 2.4. Eq. (112) gives ξ/l < (τ_H/l)^{1/Δ_t} ~ 2.5, i.e. ξ ≲ 2.5 fm, while Eq. (110) allows ξ ≲ 12 fm; the binding constraint is therefore the finite-time one, yet the paper adopts 6 fm. Because F_σ − 1 scales as ξ² (Eqs. 115–116), taking Eq. (112) at face value reduces F_σ − 1 from 0.35 to about 0.06 at μ_π = 60 MeV (and from 0.14 to about 0.024 at μ_π = 0), and after the factor-of-two dilution from non-direct pions the enhancement of √F falls to roughly 1%, comparable to the Bose (~2%) and thermal-contact (~1%) backgrounds the paper itself computes. The 'compensating good effect' invoked on p. 42 — that the large specific heat makes the system dwell longer near T_c — is qualitatively plausible but is never estimated. Please either quantify the dwell-time enhancement (for instance, estimate how much longer the system spends within |T − T_c| ≲ ξ^{−1/ν} given the divergent C_V) or soften 'easily observable' in the abstract, §5.5 and §7, and state explicitly the minimum ξ at which the signal exceeds the ~2% noncritical background.","section":"§5.4, Eqs. (110)–(112), and abstract:s9","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The inference in §3.4 that 'about 75% of the observed fluctuation is thermodynamic in origin', drawn from comparing v²_ebe(N)/⟨N⟩ = 2.008 ± 0.009 to the resonance-gas value 1.5, is not a decomposition the data support once the initial-state contributions are included. The Note added in proof concedes that Refs. [33,34] reproduce the observed multiplicity fluctuations to within a few percent with no new free parameters by combining geometric, spectator and NN-cross-section contributions with the thermodynamic ones. That being so, the 1.5-versus-2.0 comparison does not isolate a 75% thermodynamic share; the two classes of contribution are not shown to be separately identifiable. §7 nevertheless repeats the 75% figure without qualification. Please qualify both passages, or explain why the resonance-gas prediction and the geometric contributions of [33,34] are not overlapping descriptions of the same variance.","section":"§3.4 and §7, versus the Note added in proof","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The multiplicity enhancement quoted in Eq. (117) is a factor 1 + 1.0 (G/300 MeV)²(ξ/6 fm)² at μ_π = 0 (rising to 2.0 at μ_π = 60 MeV), i.e. a factor of two to three on the direct-pion variance at the central parameter values, yet the correlator (88) from which it follows is derived to lowest order in G. The paper itself states in §5 that 'the multiplicity fluctuations are large enough that in this case a treatment which goes beyond lowest order in G seems called for', but no estimate of the size of the neglected terms is given, so the reader cannot judge whether Eq. (117) is reliable to a factor of two or to an order of magnitude. Since the multiplicity channel is one of the two headline observables, please supply at least an order-of-magnitude estimate of the next-order correction (the σ-bubble resummation mentioned in §5.4 is the natural starting point) or state a quantitative reliability bound on Eq. (117).","section":"§5.1 and §5.5, Eq. (117)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The bookkeeping for √F does not close, and the residual is larger than the statistical error by several times. The data give √F = 1.002 ± 0.002 (Eq. 44); the prediction is √F_B ≈ 1.02 from Bose enhancement, reduced to about 1.01 by the thermal-contact factor √F_T ≈ 0.99 estimated in Eq. (59). The product remains roughly 1% above the measured value, i.e. about five times the quoted statistical error, and the paper attributes the remainder to two-track resolution without estimating it. Because the entire detection strategy is a percent-level nonmonotonic variation of √F on top of this baseline, the baseline should be pinned down more tightly. Please either estimate the two-track contribution or set out an explicit systematic error budget for √F, so that the size of a detectable critical bump can be stated against a known background.","section":"§3.4, Eq. (44), and §4.3, Eq. (59)","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The factor of 1/2 in Eq. (26) and the accompanying sentence that the Bose factor is 'reduced from 1 + n_p to 1 + (n_p/2)' are easy to misread, since n_p here counts both pion charges. A one-line reminder at the equation, and a note in the Fig. 1 caption that F_B is defined per charged-pion pair, would prevent confusion.","section":"§2.2, Eq. (26) and Fig. 1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that |F_res − 1| < 0.01 is quoted without an uncertainty, although it comes from slicing a simulated sample of 1.254 × 10⁶ pions into up to 2500 events. Please state the statistical error on F_res so the reader can see that the bound is meaningful rather than a consequence of limited sampling.","section":"§3.2, Table 1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The multiplicity-weighting correction (147) is applied to v_ebe(p_T) but not to v_inc(p_T), and the two are then compared as a ratio in Eq. (43). Please confirm that no analogous correction is required for the inclusive variance, or state why it cancels in the ratio.","section":"§3.4, Eqs. (42)–(43) and Appendix Eq. (147)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The choice σ²_p = v²_p in Eq. (60) is made 'for convenience', and the resulting fluctuation ⟨(ΔT̃)²⟩ in Eq. (67) is therefore definition-dependent. A sentence making explicit that other choices of σ²_p give the same mean but a different normalization of the variance would help the reader interpret the result.","section":"§4.4, Eqs. (60)–(67)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The vacuum estimate G ≈ 1900 MeV follows from λ = 20 fixed by m_σ = 600 MeV, and the suppression factor (m/Λ_QCD)^{3/5} in Eq. (106) is then applied to that number. Please state how G_E in Eq. (106) shifts if the assumed vacuum m_σ is varied over, say, 500–700 MeV, since the width argument in Eq. (93) is only a plausibility check.","section":"§5.3, Eqs. (92)–(93) and (106)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The linear parametrization m_σ(t) = 2m_π(1 + t/τ) and the assumption that sigmas decay at rest are both approximations, the second justified only qualitatively. A brief remark on how the spectrum (124) and the estimate (126) would change if the sigmas retained thermal momenta of order T_f would strengthen this section, which otherwise carries a distinctive and directly testable prediction.","section":"§6, Eqs. (119)–(126)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The central quantitative comparison rests on preliminary NA49 data cited as 'paper to be submitted' and a private communication. This is clearly labelled and is not a criticism of the analysis, but the reader should be reminded in the caption that the numbers may change, and that the quoted errors are statistical only.","section":"Table 2 and Ref. [7]","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The numerical prefactors are printed with stray decimal points — '1 + 10./(C_A/C_B + 1)' in Eq. (70) and '12.' and '13.' in Eq. (118). These should be cleaned up, and the hyphenation of 'freeze-out' versus 'freezeout' made consistent throughout.","section":"Equations (70) and (118)","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is sound and the qualitative program is likely to be influential; my recommendation turns on the fact that the abstract's 'easily observable' conclusion is stronger than the finite-time scaling estimate in the paper's own §5.4 supports, and on the 75% thermodynamic statement in §3.4 and §7, which the Note added in proof effectively retires. Both are textual fixes rather than new calculations, so I do not think a second round of refereeing is needed if the authors address them. One point for the editor's attention: the manuscript's central data comparison depends on NA49 results that were unpublished at the time of submission, which is acceptable here given the explicit labelling but is worth noting for the record."},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"This is the machinery paper behind the 1998 letter, and it holds up. The two correlators that carry the physics — Eq. (50) for thermal contact with a heat bath of finite heat capacity, and Eq. (88) for sigma exchange — are derived cleanly from Gaussian path integrals and checked against the right limits: C_A to infinity recovers the canonical result, C_A to zero enforces energy conservation, and the 1/m_sigma^2 singularity in (88) matches the Ising susceptibility exponent gamma/nu. The NA49 comparison is presented as broad agreement, not a precision test, which is the right posture given the data were preliminary. The appendix on finite-multiplicity estimators (Eq. 145) is a real, reusable contribution that outlives the critical-point program; the observation that the naive estimator biases v_ebe by O(<(Delta N)^2>/<N>^2) is the kind of thing that should have been in the experimental papers and wasn't. Credit where due: the paper is unusually explicit about what it does not know, flagging the blue-shift approximation, the tricritical extrapolation for G, and the finite-time estimate of xi_max. The classifier's 'limitations: none found' is simply wrong on the text; Section 5.4 calls Eq. (112) 'suspect for several reasons' in so many words. On the soft spot: the stress-test is right that there is an unresolved tension between Eq. (110) (xi <~ 12 fm from finite size) and Eq. (112) (xi <~ 2.5 fm from finite time), and that the paper adopts 6 fm on the strength of an unquantified 'compensating good effect' plus an unknown prefactor. Since F_sigma - 1 scales as xi^2, taking (112) at face value drops the pT signal from ~10% to a few percent after dilution — comparable to the Bose and thermal-contact backgrounds the paper itself estimates, not the 50x statistical error quoted in the abstract. G is a co-equal factor and equally soft. None of this touches the derivations; it is a quantitative-reach concern about the headline, and the nonmonotonic scan argument survives it because the background is roughly constant in the control parameter. But the paper should have made that argument explicitly rather than resting on the 6 fm number. Who should read it: anyone doing event-by-event fluctuation analysis, and anyone thinking about how to estimate finite-time critical scaling in a expanding fireball. It deserves referees, and the derivations deserve to be checked line by line by someone who knows the Ising-to-QCD mapping.","headline":"The analytic core is solid and the NA49 comparison is honest; the 'easily observable' claim leans on a correlation length the authors' own finite-time estimate does not support.","tokens_in":50577,"tokens_out":869,"would_cite":true,"duration_ms":16907,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["25.75.-q","12.38.Mh","05.70.Jk"],"model":"selfhosted/deepseek-v4.1-flash","headline":"The paper argues that event-by-event fluctuations of charged-pion multiplicity and mean transverse momentum should show a nonmonotonic bump when heavy-ion collisions freeze out near the QCD critical endpoint, and that the bump is large en","keywords":["QCD critical point","event-by-event fluctuations","heavy-ion collisions","pion multiplicity","transverse momentum fluctuations","sigma field","finite-size scaling","critical slowing down"],"falsifier":"A concrete test is an energy scan at fixed centrality: measure the scaled event-by-event variance of mean $p_T$, the multiplicity variance, the $Delta NDelta p_T$ cross correlation, and the low-$p_T$ pion yield, and check whether all show a simultaneous nonmonotonic maximum at the collision energy where the freeze-out trajectory is expected to pass the endpoint. A monotonic trend, or a maximum in only one observable with no correlated soft-pion excess, would falsify the central claim as stated.","tokens_in":50767,"feed_emoji":"⚛️","tokens_out":4496,"duration_ms":84210,"temperature":0.7,"pith_summary":"The paper's target is a practical experimental question: can event-by-event fluctuations tell whether a heavy-ion collision has frozen out near the second-order critical endpoint of the QCD phase transition? Its strategy is to first model the noncritical background as a thermal resonance gas, compare that baseline with NA49 data, and then add two critical effects: thermal contact with a heat bath whose heat capacity diverges, and direct coupling of pions to long-wavelength fluctuations of the sigma field. The authors argue that the direct sigma coupling dominates and produces a nonmonotonic increase in the variance of pion multiplicity and mean transverse momentum as a control parameter such as collision energy moves the freeze-out point toward, through, and past the endpoint. They estimate that finite size and finite time cap the sigma correlation length at a few femtometers, but that the resulting signal can still be an order of magnitude above present statistical errors.","feed_headline":"A bump in pion fluctuations could mark QCD's critical point","feed_subtitle":"If the sigma field goes critical, event-by-event pion spreads should rise, then fall, as collision energy is scanned.","key_machinery":"The argument is carried by the microscopic correlator $langleDelta n_pDelta n_krangle$, the covariance between fluctuations of pion occupation numbers in momentum modes $p$ and $k$. In an ideal Bose gas this is diagonal, $v_p^2delta_{pk}$ with $v_p^2=langle n_prangle(1+langle n_prangle)$; resonance decays add small off-diagonal terms; thermal contact with a finite heat capacity $C_A$ adds a negative off-diagonal term from energy conservation; and $\\sigma$ exchange adds the critical term $frac{1}{V m_sigma^2}G^2T,frac{v_p^2}{omega_p}frac{v_k^2}{omega_k}$. Once this correlator is known, any observable built from pion occupation numbers follows by weighting with the single-particle quantity and dividing by $langle Nrangle^2$. The finite-size and finite-time discussion then replaces the divergent $1/m_sigma^2$ by a capped correlation length $xi_{max}$, using scaling arguments $xi_{max}sim R$ and $xi_{max}simtau^{1/Delta_t}$.","core_discovery":"The central claim is that the QCD critical endpoint leaves a measurable imprint on the pions already observed in heavy-ion collisions. Near the endpoint the $\\sigma$ field becomes light and its long-wavelength fluctuations grow; through the $Gsigma pi pi$ coupling these fluctuations add a term to the pion occupation-number correlator proportional to $G^2 T/m_sigma^2$, which diverges as $m_sigma to 0$. The paper turns this into estimates for the scaled event-by-event variance of mean transverse momentum, $F_sigma=1+0.35(G/300,mathrm{MeV})^2(xi/6,mathrm{fm})^2$ at $mu_pi=60,mathrm{MeV}$, and for the multiplicity variance, with both growing near the endpoint. It also argues that the noncritical baseline—Bose enhancement, resonance decays, flow fluctuations, and energy-conservation anticorrelations—is broadly consistent with NA49 data, so that a nonmonotonic departure from that baseline would be evidence for critical freeze-out.","pith_inferences":["A null result in an energy scan would not by itself rule out the endpoint; it could mean the correlation length stayed below the finite-time cap, so the measurement would bound the combination $Gxi$ rather than the existence of the endpoint.","Because the predicted signal depends quadratically on $xi$, an independent interferometric or fluctuation-based estimate of the homogeneous patch size at freeze-out would sharpen the test considerably.","The same correlator logic suggests that a differential scan in pion momentum, rather than a single acceptance-averaged number, would separate the sigma-exchange term from the Bose and resonance backgrounds.","If the soft-pion excess from late sigma decays is seen, its angular and charge correlations could distinguish a critical source from ordinary resonance decays."],"forward_implications":["As collision energy, centrality, or ion size is varied, the scaled $p_T$ variance and the multiplicity variance should first rise, peak at the same setting, and then fall.","The critical increase should be strongest for soft pions, so a low-$p_T$ cut should amplify the nonmonotonic signal relative to the full-acceptance variance.","The cross correlation $langleDelta NDelta p_Trangle$, normally small, should grow in magnitude near the endpoint because it receives a large sigma-exchange contribution.","If freeze-out is near the endpoint, a nonthermal population of pions with mean $p_Tsim m_pi/2$ should appear from late sigma decays and then disappear as the endpoint is passed.","Off-diagonal correlations between $pi^+$ and $pi^-$ fluctuations should grow near the endpoint, since the diagonal Bose term is absent in that channel."],"supporting_citations":[{"why":"Prior letter that proposes the critical-endpoint signatures and the nonmonotonic search strategy this paper develops.","marker":"[1]"},{"why":"NA49 preliminary data used as the noncritical baseline for multiplicity and mean-$p_T$ fluctuations.","marker":"[7]"},{"why":"Supplies the universal critical properties and Ising-like exponents used to characterize the endpoint.","marker":"[2]"},{"why":"Provides the phase-diagram and tricritical-point context for the endpoint.","marker":"[3]"},{"why":"Hydrodynamic freeze-out parameters and pion chemical potential used in the numerical estimates.","marker":"[20]"},{"why":"Chemical freeze-out temperature and baryon chemical potential used in the resonance-gas baseline.","marker":"[23]"},{"why":"Two-pion interferometry patch sizes used to estimate the finite-size cap on the correlation length.","marker":"[25]"},{"why":"Dynamic universality classification and exponent used in the finite-time scaling estimate.","marker":"[32]"},{"why":"Linear sigma model used to estimate the vacuum $sigma pi pi$ coupling $G$.","marker":"[29]"}],"fun_headline_variants":["Pion spreads could bend near the QCD critical endpoint","Sigma-pion coupling adds a critical term to pion correlations","NA49 fits the noncritical baseline, leaving critical pion bumps to find","Finite size and time still leave critical pion fluctuations observable","Multiplicity and mean-pT variances should turn nonmonotonic at the endpoint"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The load-bearing premise is that freeze-out occurs from an equilibrated hadronic system close enough to the critical endpoint that the $\\sigma$ correlation length reaches several femtometers, because the predicted enhancement scales as $xi^2$ and a correlation length of only about $2,mathrm{fm}$ would shrink the $p_T$ signal to roughly the percent level.","fun_headline_variants_meta":{"raw":{"variants":["Pion spreads could bend near the QCD critical endpoint","Sigma-pion coupling adds a critical term to pion correlations","NA49 fits the noncritical baseline, leaving critical pion bumps to find","Finite size and time still leave critical pion fluctuations observable","Multiplicity and mean-pT variances should turn nonmonotonic at the endpoint"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":1790,"prompt_tokens":1062,"completion_tokens":728,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":652}},"tokens_in":1062,"tokens_out":728,"duration_ms":13361,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":652,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-29T18:32:41.105257+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"falsifier":"A concrete test is an energy scan at fixed centrality: measure the scaled event-by-event variance of mean $p_T$, the multiplicity variance, the $Delta NDelta p_T$ cross correlation, and the low-$p_T$ pion yield, and check whether all show a simultaneous nonmonotonic maximum at the collision energy where the freeze-out trajectory is expected to pass the endpoint. A monotonic trend, or a maximum in only one observable with no correlated soft-pion excess, would falsify the central claim as stated.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior letter that proposes the critical-endpoint signatures and the nonmonotonic search strategy this paper develops."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"NA49 preliminary data used as the noncritical baseline for multiplicity and mean-$p_T$ fluctuations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the universal critical properties and Ising-like exponents used to characterize the endpoint."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the phase-diagram and tricritical-point context for the endpoint."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Hydrodynamic freeze-out parameters and pion chemical potential used in the numerical estimates."},{"cited_title":"Braun-Munzinger and J","cited_arxiv_id":null,"evidence_quote":"Chemical freeze-out temperature and baryon chemical potential used in the resonance-gas baseline."},{"cited_title":"Two-Particle Correlations in Relativistic Heavy-Ion Collisions","cited_arxiv_id":"nucl-th/9902020","evidence_quote":"Two-pion interferometry patch sizes used to estimate the finite-size cap on the correlation length."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Dynamic universality classification and exponent used in the finite-time scaling estimate."},{"cited_title":"Gell-Mann and M","cited_arxiv_id":null,"evidence_quote":"Linear sigma model used to estimate the vacuum $sigma pi pi$ coupling $G$."}],"review_version":1}