{"id":"c8535e81-4e1f-4fef-ac72-ead8c252af2d","arxiv_id":"1908.04208","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Charged-particle spectra from pp collisions at 5.02 and 13 TeV are described by Tsallis and blast-wave fits; the extracted temperature and non-extensivity rise with multiplicity while radial flow stays flat.","lead":"This paper fits newly published ALICE data on charged particle production in proton-proton collisions at 5.02 and 13 TeV with two statistical models, extracting temperature and flow parameters. Tsallis statistics describes the full spectrum, a blast-wave model only the bulk, and the extracted parameters shift with multiplicity, feeding the debate on whether small collisions form a quark-gluon plasma.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Tsallis q and T trends are not yet demonstrated to be physical: the paper's own fit-range scan (Figs. 5-7) shows both parameters depend strongly on the arbitrary upper pT cutoff, especially for the high-multiplicity class.","rationale":"Read in good faith, this is a standard phenomenological analysis of public ALICE data. The authors are transparent about their fit-range study, and the algebra of Eqs. (1)-(11) is internally consistent and follows the established thermodynamically consistent Tsallis form. The weak point is evidential rather than algebraic: the paper's own Figs. 5-7 show that q and T vary substantially with the upper pT cutoff, and for V0M1 the fit quality deteriorates steeply above about 2 GeV/c. Since the abstract and summary present q and T trends as physical freeze-out properties, and since no uncertainties accompany any extracted parameter, the reader cannot tell whether those trends are robust or are consequences of always fitting to 20 GeV/c. This is a concrete, testable concern, and it matches the reader's weakest assumption about the arbitrary fit range. It does not rise to a rejection-level flaw because the study is exploratory and the authors partially flag the fit-range sensitivity. A fixed-range re-analysis with uncertainties would settle the issue. Therefore the reader's CONDITIONAL verdict is appropriate and unchanged.","tokens_in":10584,"tokens_out":5950,"duration_ms":63403,"concrete_test":"Re-fit all ten V0M classes at both 5.02 and 13 TeV with Eq. (11), holding the upper pT boundary fixed at 3 GeV/c, where chi2/ndf is still acceptable for all classes, and separately at 20 GeV/c; attach bootstrap or MINUIT uncertainties to each q and T. If the multiplicity and energy orderings of q and T in Figs. 3-4 survive both boundaries within uncertainties, the central trend claims stand. If the orderings reorder or lose significance, the headline claims are fit-range artifacts.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central thermodynamic claims: q rising and saturating with dNch/deta, T increasing with multiplicity, and the energy ordering of q, are read from fits over a single, arbitrarily chosen interval of 0 to 20 GeV/c (Sec. III A). The paper's own fit-range study undermines that choice. Fig. 5 shows that for the high-multiplicity V0M1 class, chi2/ndf is acceptable only up to about pT = 2 GeV/c and then grows strongly, while for low multiplicity it stays acceptable over the full range. Fig. 6 shows q for V0M1 decreases as the upper fit boundary is extended, whereas q for V0M10 increases and saturates above 8 GeV/c. Fig. 7 shows the Tsallis temperature T moving in opposite directions for the two classes. Because no parameter uncertainties are given, the multiplicity ordering in Figs. 3 and 4 could be a consequence of choosing pT,max = 20 GeV/c for all classes rather than a property of the collision system. The BGBW conclusions are similarly vulnerable: no errors are quoted, and the model deviates visibly from data at low pT, so the statements that <beta> is almost independent of multiplicity and that Tkin depends on multiplicity are not statistically supported. This does not disprove the fits, but it makes the thermodynamic interpretation conditional on an unexamined fit-range choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper fits the ALICE inclusive charged-particle transverse momentum spectra in pp collisions at sqrt(s) = 5.02 and 13 TeV, in ten V0M multiplicity classes, using two functional forms: a 'thermodynamically consistent' Tsallis distribution (Sec. II A) and the Boltzmann-Gibbs Blast Wave (BGBW) model (Sec. II B). From the Tsallis fits the authors extract the non-extensive parameter q and the Tsallis temperature T as functions of charged-particle multiplicity and collision energy, and they carry out a fit-range scan for the highest and lowest multiplicity classes. From the BGBW fits to the bulk region (pT up to about 2.5 GeV/c) they extract the average radial flow <beta> and the kinetic freeze-out temperature Tkin. The main claims are that the Tsallis form describes the complete pT range better than BGBW, that q increases with multiplicity and saturates above dNch/deta around 15, that q is larger at 13 TeV than at 5.02 TeV, that T increases with multiplicity, that <beta> is almost independent of multiplicity and collision energy, and that Tkin depends clearly on multiplicity class.","tokens_in":10819,"tokens_out":4840,"duration_ms":54171,"significance":"If the claimed trends were quantitatively established, this would be a useful phenomenological map of freeze-out parameters across multiplicity classes in small systems, directly relevant to current discussions of collectivity and equilibration in high-multiplicity pp collisions. The paper has real strengths: it uses recent ALICE data, shows data/fit ratios in the lower panels, reports chi2/ndf values, and includes a fit-range study that is often missing from such analyses. Its main weaknesses are that the extracted parameters are presented without uncertainties and that the paper's own fit-range scan shows strong variation of q and T with the upper pT cutoff for the high-multiplicity class. The transition from 'we observe a trend' to 'the data demonstrate a property of the collision system' therefore requires error propagation and a defended choice of fit range. The paper does not provide code, parameter tables with uncertainties, or independent validation of the thermodynamic form, which limits reproducibility.","major_comments":[{"comment":"The central multiplicity and energy trends of q and T are read from fits performed with a single upper limit of pT,max = 20 GeV/c for all classes, while the paper's own fit-range scan shows that for the high-multiplicity V0M1 class the extracted q decreases and T increases as the upper boundary is extended (Figs. 6 and 7), and chi2/ndf grows strongly above roughly pT = 2 GeV/c (Fig. 5). Because Figs. 3 and 4 show no parameter uncertainties, the reported ordering of q and T across multiplicity classes could be an artifact of the chosen fit range rather than a property of the collision system. Please report fit uncertainties and either use a common fit range in which the chi2/ndf is acceptable for all classes or demonstrate explicitly that the trends in Figs. 3 and 4 persist when the same physically motivated pT,max is used consistently.","section":"Sec. III A, Figs. 3-7"},{"comment":"The thermodynamic interpretation of q and T rests on the 'addition power of q' in Eq. (1), which is imported from the authors' own earlier papers (refs. [3,4,6]) rather than derived or independently verified here. If that form is not the correct thermostatistical normalization for the system, the fitted q and T are only shape parameters of an ad-hoc function and the abstract's thermodynamic claims do not follow. As a concrete test, please provide the derivation of Eq. (6) or a precise citation to a derivation, and perform a robustness check by fitting the same data with the standard Tsallis form (without the extra power of q) and showing whether the q and T trends and the fit-range behavior survive.","section":"Sec. II A, Eqs. (1)-(6) and (11)"},{"comment":"The BGBW conclusions that <beta> is almost independent of multiplicity and energy while Tkin depends clearly on multiplicity are stated without parameter uncertainties, so 'almost independent' cannot be distinguished from 'consistent within large errors'. In addition, the bottom panels of Figs. 8 and 9 show a systematic low-pT deviation, and the interpretation that this is due to resonance decays is not demonstrated. Please report the fit uncertainties for beta_s, <beta>, and Tkin, provide per-class chi2/ndf values, and state whether the observed Tkin variation is statistically significant relative to those uncertainties.","section":"Sec. III B, Figs. 10-12"}],"minor_comments":[{"comment":"The phrase 'an addition power of q' should read 'an additional power of q'.","section":"Sec. II A, after Eq. (2)"},{"comment":"The ordinate label chi2/NDF is not defined; stating whether NDF is the number of data points minus the number of fit parameters would make the goodness-of-fit values interpretable.","section":"Sec. III A, Fig. 5"},{"comment":"The bullet 'pT ~ 8 GeV/c ... needs a closure look' contains a typo ('closure' should be 'closer') and the claimed threshold at 8 GeV/c is not supported by a quantitative test of, for example, a change in slope or saturation.","section":"Sec. IV, summary bullets"},{"comment":"The pion, kaon, and proton weights (0.8, 0.12, 0.08) are stated to follow experimental yields, but no uncertainties or reference values are given; a short table of weights and the result of the pion-fraction sensitivity test would make this statement more concrete.","section":"Sec. II A, Eq. (8)"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the manuscript is within the journal's scope, but the central phenomenological claims are currently supported only by fits without error bars and with a fit-range dependence that the authors themselves document. The overlap with the authors' prior papers (refs. [3,4,6,12,13]) is substantial, and the 'thermodynamically consistent' form is taken from those papers rather than re-derived; I would ask that the refereeing process ensure those cited derivations are available and that the revision either adds rigorous error propagation and a defended fit-range choice or substantially softens the thermodynamic interpretation. No concerns about research integrity beyond the above are raised."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"For an internal take: this is a solid, conventional fit analysis, not a discovery. It extends established Tsallis and BGBW descriptions to the 2019 ALICE pp spectra at 5.02 and 13 TeV across ten V0M multiplicity classes, and it adds a systematic scan over the upper pT fit boundary. That fit-range study is the most useful part of the paper.\n\nWhat is actually new: the application to these data, plus the explicit demonstration that q and T move with pT,max and move in opposite directions for high and low multiplicity. The fits are honestly presented with ratio panels and chi2/ndf, and the authors are transparent that the high-multiplicity Tsallis fit degrades above about 2 GeV/c. Credit is due for showing that rather than hiding it.\n\nThe soft spots are real but not fatal. First, the extracted parameters have no error bars anywhere, so the multiplicity and energy trends in Figs. 3, 4, 11, and 12 cannot be assessed statistically. Second, the headline values come from a single pT < 20 GeV/c choice for all classes, while the paper's own Fig. 5 shows that this is not a safe range for the high-multiplicity class. The 8 GeV/c threshold is promoted to a summary bullet even though the text says it needs a closer look. Third, the BGBW comparison fits a region where the model visibly underpredicts low-pT data, so claiming <beta> is multiplicity-independent without errors is too strong. On the thermodynamics: the q-dual form is cited to the authors' own papers, not re-derived, but this is a literature assumption rather than a hidden circular step.\n\nDoes the central argument hold? The descriptive claim does: Tsallis describes the full spectra, BGBW only the bulk part. The physical interpretation of q and T as freeze-out thermodynamics is conditional, not demonstrated. The paper is transparent enough that a careful reader can see the distinction.\n\nThis is for practitioners who do Tsallis fits of LHC small-system spectra. It is a useful reference on fitting-range behavior and parameter dependences, less useful for physics conclusions. I would send it to peer review; a decent referee can require uncertainty estimates, a fit-range robustness statement, and softer wording on the threshold and on the BGBW conclusions.","headline":"Competent Tsallis and BGBW fits to new ALICE pp data, but the thermodynamic trends are only as solid as the arbitrary fit range and the missing error bars.","tokens_in":11508,"tokens_out":2471,"would_cite":false,"duration_ms":27469,"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":"This paper reports that a thermodynamically consistent Tsallis distribution describes the full transverse-momentum spectra of charged particles in proton-proton collisions at 5.02 and 13 TeV, while a blast-wave model captures only the…","keywords":["Tsallis non-extensive statistics","transverse momentum spectra","pp collisions","LHC","kinetic freeze-out temperature","radial flow","blast-wave model","charged particle multiplicity"],"falsifier":"Re-fit the $pp$ spectra at $\\sqrt{s}=5.02$ and 13 TeV with a fixed upper boundary, say $p_{T,\\max}=3$ GeV/$c$, for every multiplicity class; if $q$ no longer saturates above $\\mathrm{d}N_{\\rm ch}/\\mathrm{d}\\eta\\sim15$ or $T_{\\rm kin}$ no longer rises with multiplicity, the reported trends are consequences of the changing fitting range rather than properties of the collision system.","tokens_in":10273,"feed_emoji":"📈","tokens_out":14112,"duration_ms":127143,"temperature":0.7,"pith_summary":"Using measured charged-particle transverse-momentum ($p_{T}$) spectra in $pp$ collisions at $\\sqrt{s}=5.02$ and 13 TeV, this paper asks which statistical description governs the produced hadronic system: a non-extensive Tsallis distribution or a Boltzmann-Gibbs blast-wave model. It reports that the thermodynamically consistent Tsallis form fits the complete $p_{T}$ spectra up to 20 GeV/$c$, whereas the blast-wave model covers only the bulk region below about 2.5 GeV/$c$. The fitted non-extensive parameter $q$ rises with event multiplicity and saturates above $\\mathrm{d}N_{\\rm ch}/\\mathrm{d}\\eta \\sim 15$, with larger values at 13 TeV; the fitted temperature also grows with multiplicity. The blast-wave fits show radial flow that is essentially independent of collision energy and multiplicity, alongside a kinetic freeze-out temperature that clearly increases with multiplicity. This matters because it separates genuine freeze-out trends in small collision systems from fitting-function artifacts, and it identifies empirical thresholds ($p_{T}\\sim 8$ GeV/$c$; $\\mathrm{d}N_{\\rm ch}/\\mathrm{d}\\eta\\sim15$) in particle production.","feed_headline":"Tsallis statistics fits full pp momentum spectra to 20 GeV/c","feed_subtitle":"A blast-wave model covers only the bulk; freeze-out temperature grows with multiplicity while radial flow stays flat.","key_machinery":"The load-bearing object is the thermodynamically consistent Tsallis distribution, Eq. (11): a two-parameter $(T,q)$ invariant yield for charged particles in which the power-law exponent carries an extra factor of $q$ so that entropy, energy, pressure, and particle number satisfy standard thermodynamic relations. Fitting this form to the $p_{T}$ spectra yields the Tsallis temperature $T$ and the non-extensive parameter $q$, interpreted as kinetic freeze-out temperature and distance from equilibrium. The comparison mechanism is the BGBW model, Eq. (14), using modified Bessel functions $K_{1}$ and $I_{0}$ with a linear flow profile, which gives the average radial flow $\\langle\\beta\\rangle$ and $T_{\\rm kin}$ from the bulk $p_{T}$ region. The paper's central systematic tool is a scan of the upper $p_{T}$ fitting boundary from 1.2 to 20 GeV/$c$, showing opposite $q$ and $T$ trends for the highest and lowest multiplicity classes and a saturation threshold near 8 GeV/$c$.","core_discovery":"The paper's central claim is that the charged-particle $p_{T}$ spectra in $pp$ collisions at $\\sqrt{s}=5.02$ and 13 TeV separate into two regimes: a bulk component below $\\sim 2.5$ GeV/$c$ captured by the Boltzmann-Gibbs blast-wave model, and a full-range component captured by Tsallis non-extensive statistics. Using the thermodynamically consistent Tsallis form (with an extra power of $q$ in the exponent), the fitted $q$ rises from about 1.12 to 1.18 with multiplicity and saturates above $\\mathrm{d}N_{\\rm ch}/\\mathrm{d}\\eta \\sim 15$, while the fitted temperature $T$ rises into the $\\sim 0.1$ GeV range and is higher at 13 TeV. The BGBW fits yield a radial-flow velocity $\\langle\\beta\\rangle$ that is essentially independent of collision energy and multiplicity, and a kinetic freeze-out temperature $T_{\\rm kin}$ that increases with multiplicity. The paper further claims that $p_{T}\\sim 8$ GeV/$c$ is a threshold: above it, $q$, $T$ and $\\chi^{2}/{\\rm ndf}$ saturate, and high-multiplicity fits worsen when the upper fit boundary exceeds $\\sim 2$ GeV/$c$.","pith_inferences":["Re-fitting the same spectra with an event-shape or jet-veto selection, as the paper itself suggests, would test whether the $q$ rise at high multiplicity comes from jet fragmentation rather than from genuine non-extensive equilibration.","The demonstrated sensitivity of $q$ and $T$ to the upper $p_{T}$ cutoff implies that freeze-out parameters extracted with different fitting windows should not be compared directly; a fixed, pre-registered window would make cross-energy and cross-experiment comparisons cleaner.","If the saturation threshold near $\\mathrm{d}N_{\\rm ch}/\\mathrm{d}\\eta\\sim15$ is real, it gives a concrete target for multi-parton-interaction models and for future LHC Run 3 $pp$ measurements at higher multiplicities.","Applying the BGBW fit separately to identified pions, kaons, and protons would test whether the multiplicity-independent radial flow and multiplicity-dependent $T_{\\rm kin}$ persist per species, rather than only in the charged-particle sum with fixed weight factors."],"forward_implications":["At both collision energies, the Tsallis distribution tracks the measured charged-particle spectra over the full measured $p_{T}$ range in all multiplicity classes, while the BGBW model underpredicts the low-$p_{T}$ region below about 0.3 GeV/$c$, attributed to resonance decays.","The non-extensive parameter $q$ and the Tsallis temperature $T$ both increase with charged-particle multiplicity, with $q$ saturating above $\\mathrm{d}N_{\\rm ch}/\\mathrm{d}\\eta\\sim15$ and taking higher values at 13 TeV than at 5.02 TeV.","The BGBW analysis gives a radial-flow velocity $\\langle\\beta\\rangle$ that is almost constant across collision energy and multiplicity, and a kinetic freeze-out temperature $T_{\\rm kin}$ that rises with multiplicity.","Fitting-range scans show that $q$ and $T$ move in opposite directions for the highest (V0M1) and lowest (V0M10) multiplicity classes as the upper $p_{T}$ cutoff grows, with $p_{T}\\sim8$ GeV/$c$ marking a saturation threshold; for V0M1, $\\chi^{2}/{\\rm ndf}$ degrades once the fit extends beyond about 2 GeV/$c$.","Together these results support a two-component picture: the bulk of the system behaves thermally with a multiplicity-dependent freeze-out temperature, while the high-$p_{T}$ tail requires a non-extensive power-law description."],"supporting_citations":[{"why":"Provides the thermodynamically consistent Tsallis distribution with the extra power of q that is fitted throughout.","marker":"[3]"},{"why":"Companion paper establishing the same thermodynamic consistency conditions for Tsallis statistics.","marker":"[4]"},{"why":"Demonstrates the Tsallis-fit approach to LHC transverse-momentum spectra that this analysis applies.","marker":"[6]"},{"why":"Earlier Tsallis fits to pion spectra at LHC energies, showing the saturation of q with multiplicity that the present charged-particle analysis extends.","marker":"[12]"},{"why":"Supplies the experimental pp spectra at 5.02 and 13 TeV and the V0M multiplicity classes used for every fit.","marker":"[20]"},{"why":"Provides the analytic pT integral that lets the fit normalize by mid-rapidity yield and eliminate the volume parameter.","marker":"[21]"},{"why":"Gives the identified-particle yields used to fix pion, kaon, and proton weight factors of 0.8, 0.12, and 0.08.","marker":"[22]"},{"why":"Defines the Boltzmann-Gibbs blast-wave formula, with modified Bessel functions, used for the bulk-region fits.","marker":"[23]"},{"why":"Supplies the Bjorken rapidity correlation used in the BGBW reduction to the measured pT distribution.","marker":"[24]"}],"fun_headline_variants":["Tsallis outshines blast-wave for full pp pT range","Freeze-out temperature scales with multiplicity, radial flow does not","pp collisions: blast-wave handles bulk, Tsallis covers everything","Tsallis parameters saturate beyond 8 GeV/c in pp spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results stand or fall on the choice of the thermodynamically consistent Tsallis formula as the correct statistics for the system: if that functional form is wrong, the fitted temperatures and $q$-values are properties of the fitting function, not of the collisions.","fun_headline_variants_meta":{"raw":{"variants":["Tsallis outshines blast-wave for full pp pT range","Freeze-out temperature scales with multiplicity, radial flow does not","pp collisions: blast-wave handles bulk, Tsallis covers everything","Tsallis parameters saturate beyond 8 GeV/c in pp spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00026,"raw_usage":{"total_tokens":1645,"prompt_tokens":1057,"completion_tokens":588,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":673,"completion_tokens_details":{"reasoning_tokens":515}},"tokens_in":673,"tokens_out":588,"duration_ms":6786,"temperature":1.0,"reasoning_tokens":515,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:49:16.283302+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the $pp$ spectra at $\\sqrt{s}=5.02$ and 13 TeV with a fixed upper boundary, say $p_{T,\\max}=3$ GeV/$c$, for every multiplicity class; if $q$ no longer saturates above $\\mathrm{d}N_{\\rm ch}/\\mathrm{d}\\eta\\sim15$ or $T_{\\rm kin}$ no longer rises with multiplicity, the reported trends are consequences of the changing fitting range rather than properties of the collision system.","supporting_citations":[{"cited_title":"Cleymans and D","cited_arxiv_id":null,"evidence_quote":"Provides the thermodynamically consistent Tsallis distribution with the extra power of q that is fitted throughout."},{"cited_title":"Cleymans and D","cited_arxiv_id":null,"evidence_quote":"Companion paper establishing the same thermodynamic consistency conditions for Tsallis statistics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates the Tsallis-fit approach to LHC transverse-momentum spectra that this analysis applies."},{"cited_title":"Khuntia, H","cited_arxiv_id":null,"evidence_quote":"Earlier Tsallis fits to pion spectra at LHC energies, showing the saturation of q with multiplicity that the present charged-particle analysis extends."},{"cited_title":"Acharya et al","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental pp spectra at 5.02 and 13 TeV and the V0M multiplicity classes used for every fit."},{"cited_title":"Rybczynski, Z","cited_arxiv_id":null,"evidence_quote":"Provides the analytic pT integral that lets the fit normalize by mid-rapidity yield and eliminate the volume parameter."},{"cited_title":"Acharya et al","cited_arxiv_id":null,"evidence_quote":"Gives the identified-particle yields used to fix pion, kaon, and proton weight factors of 0.8, 0.12, and 0.08."},{"cited_title":"Schnedermann, J","cited_arxiv_id":null,"evidence_quote":"Defines the Boltzmann-Gibbs blast-wave formula, with modified Bessel functions, used for the bulk-region fits."}],"review_version":1}