REVIEW 4 major objections 5 minor 27 references
A two-component split of proton collision spectra keeps each piece's average momentum flat as multiplicity grows.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-04 10:47 UTC pith:YO2SQEKC
load-bearing objection A cleanly presented two-component fit whose headline flatness is likely an artifact of the multiplicity-dependent p0 window; worth refereeing but needs major revision. the 4 major comments →
Revisiting the soft-hard separation in the transverse momentum spectra of pp collisions
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that after fitting a single exponential (Boltzmann) function to the low-transverse-momentum region of each multiplicity class and subtracting it, the residual spectra—identified with hard QCD-like fragmentation—show no significant evolution in shape or peak position across multiplicity classes, and both the soft and hard mean transverse momenta are nearly multiplicity-independent. This contradicts the common expectation that the rise of the total mean transverse momentum with multiplicity reflects a change in the soft/thermal sector; instead it reflects an increasing fraction of hard fragments. The decomposition is demonstrated at three collision energies and is reproduc
What carries the argument
The key object is the cut parameter p0, the upper limit of the transverse-momentum window over which a single exponential (Boltzmann) function is fitted. For each multiplicity class, p0 is chosen as the value where the fit's chi-squared per degree of freedom is closest to one, marking the presumed boundary below which hard fragmentation is negligible. The fitted exponential is the soft component, and the residual after subtraction is the hard component. The stability of the extracted components across energies and multiplicity classes, plus agreement with Monte Carlo simulations, carries the argument.
Load-bearing premise
The decomposition assumes there is a low-transverse-momentum window, ending at p0, in which the soft contribution is a single exponential and hard fragmentation contributes negligibly; if jet-fragment particles leak into that window, the fitted 'soft' component absorbs some hard yield and the observed flatness becomes a fitting artifact.
What would settle it
Measure the yield of jet-fragmentation hadrons inside the low-transverse-momentum window [0, p0]—for example by tagging jets or using jet-veto techniques. If that yield is a significant fraction of the spectrum, or if the extracted soft slope changes appreciably when p0 is varied within the region where chi-squared is near one, the separation is not robust.
If this is right
- The rise of the total mean transverse momentum with multiplicity is driven by a growing weight of hard processes, not by a change in the soft component.
- The soft component has an approximately constant effective temperature across multiplicity classes.
- The fragmentation tails keep a stable shape, so no intermediate source such as collective flow or recombination is needed to describe the data.
- The two-component description is a viable alternative to hydrodynamic interpretations of proton-proton spectra.
- A standard Monte Carlo event generator reproduces both the total and separated trends, supporting the decomposition method.
Where Pith is reading between the lines
- If the flat soft mean transverse momentum holds for identified particle species, it would allow a cleaner separation of thermal-like and fragmentation contributions in searches for collective effects in small collision systems.
- The procedure suggests a practical way to define a 'contamination-free' low-transverse-momentum window that could serve as a baseline for calibrating underlying-event models in Monte Carlo generators.
- One could test the decomposition at even higher multiplicities or in proton-nucleus collisions; if the soft mean transverse momentum remains flat, the two-component interpretation would be strengthened.
- Applying the same fitting strategy to identified-particle spectra (pions, kaons, protons) would reveal whether the flatness is universal or species-dependent.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript revisits the two-component (soft + hard) description of charged-particle transverse momentum spectra in pp collisions at ALICE. The soft component is modeled as a Boltzmann exponential A exp(-β pT) fitted over a low-pT window [0, p0]; the residual spectrum is identified as hard QCD fragmentation. The authors apply this decomposition to ALICE data at 2.76, 5.02 and 13 TeV, and compare with Pythia 8 (Monash tune). Their central claim is that after separation both the soft and hard components have nearly multiplicity-independent mean transverse momenta, while the total <pT> rises with multiplicity, supporting a two-component picture over a hydrodynamic interpretation. They also state that a comparison of Pythia simulations with and without color reconnection demonstrates robustness.
Significance. If established, the result would strengthen the two-component interpretation of pp spectra and provide a concrete, falsifiable alternative to blast-wave/hydrodynamic descriptions. The paper is methodologically transparent in stating its operational definitions and in cautioning that 'soft' and 'hard' are labels, not uniquely identified dynamical components. Its strengths include the simplicity of the fit ansatz, the use of published ALICE data, and the attempt to validate trends with Pythia 8. However, the central claim is currently under-determined: the flatness of the soft mean is largely a restatement of the fitted inverse slope, and the hard component is defined as a residual without an independent shape fit. The analysis window p0 is multiplicity-dependent and its values and uncertainties are not reported, so the observed flatness of <pT> in either component may be an artifact of the chosen boundary. The claimed color-reconnection comparison is absent from the body of the paper. These are fixable with additional quantitative analysis and reporting, but they are load-bearing for the paper's main conclusion.
major comments (4)
- [Sec. II, Figs. 1 and 4] The central claim that both separated components have nearly multiplicity-independent <pT> rests on the multiplicity-dependent choice of p0, but the paper does not report p0 values, integration limits, or their uncertainties. The text states that at high multiplicity the Boltzmann fit worsens 'already at lower p0 values', implying p0 decreases with multiplicity. For a power-law tail ~ pT^{-n}, the hard mean over [p0,∞) scales as p0 (n-1)/(n-2); a decrease from ~1.1 to ~0.8 GeV/c would lower the hard mean by roughly 25%, potentially masking a genuine rising trend. Similarly, the soft mean of an exponential is identically 1/β, and a shrinking fit window can bias β. The paper must report p0 for every multiplicity class and energy, propagate its uncertainty, and repeat the separation with a fixed p0 as a robustness test.
- [Eq. (1) and Sec. II] The flatness of the soft <pT> is in large part an algebraic consequence of the exponential ansatz rather than an independent empirical finding. Since <pT>_soft = 1/β for f(pT)=A exp(-β pT), the claim 'the soft component shows a weak dependence on multiplicity' is essentially the statement that the fitted β is nearly constant. This is not circular in a logical sense, but it is definitional and needs to be supported by reporting β(pT, multiplicity) and A values, including their fit uncertainties, and by showing that a different soft ansatz with the same p0 prescription does not alter the conclusion. The hard component, moreover, is not fit independently; it is defined as the residual. The paper should fit the subtracted spectra with an explicit hard form (e.g., a power law or PQCD-like fragmentation function) and report the resulting parameters.
- [Abstract and Sec. II, Fig. 4] The abstract promises that 'the robustness of the decomposition is demonstrated by a comparison between simulations with and without color reconnection, yielding consistent results.' No such comparison appears anywhere in the body of the manuscript. Only the Monash-tune Pythia 8 comparison is shown in Fig. 4. Because the color-reconnection comparison is advertised as the main robustness check, it must either be included with a full description of the simulation and separation procedure, or the abstract and discussion must be revised. Without this, the validation claim is unverifiable.
- [Figs. 2-4 and Sec. II] No uncertainties are shown or reported for the experimental data points or for the fitted soft/hard <pT> values in Fig. 4. The claim that the separated components are 'nearly constant' as a function of multiplicity is not quantitatively assessable without error bars, fit-parameter uncertainties, and systematic uncertainties associated with the p0 choice and the subtraction procedure. In addition, the criterion 'χ2/ndf closest to one' is not a rigorous model-selection rule; the authors should show how p0 and the resulting <pT> values vary under reasonable variations of this criterion, and report the χ2/ndf values around the selected p0.
minor comments (5)
- [Fig. 3 caption] The caption says 'Solid lines correspond to the original Boltzmann fits' for subtracted spectra. This is confusing because the subtracted spectra should not contain the Boltzmann component; please clarify whether the lines indicate the fitted soft component before subtraction or a reference curve.
- [Fig. 4 caption] The lower panel is described as 'the fragmentation and soft contributions from top to bottom', but the ordering is not clear from the caption. Please label the curves directly or use explicit legend entries.
- [Sec. II] The phrase 'the spectra still carry some contribution from the soft QCD component' in the discussion of Fig. 3 seems to contradict the intended interpretation of the subtracted spectra; please rephrase to clarify what is meant.
- [References] Reference [25] (Gardim, Giannini, Ollitrault) is about accessing the speed of sound with mean transverse momentum; the connection to the present exponential-slope interpretation is not explained and may be misleading. Please either elaborate the connection or remove the citation.
- [Throughout] There are several formatting issues, including the running header 'pp collisions' and inconsistent spacing in the abstract. These should be corrected in the final version.
Circularity Check
No significant circularity: the decomposition is a documented fit ansatz, and the reported trends are not forced by an equation-level identity.
full rationale
The paper's derivation chain is an explicit phenomenological fit, not a first-principles claim. It assumes f(pT)=A exp(-β pT) (Eq. 1) below a multiplicity-dependent p0 selected by χ², subtracts the fit, and interprets the residual as hard fragmentation. The reported flatness of the soft mean is algebraically the same information as the fitted β (for an exponential dN/dpT, ⟨pT⟩=1/β), but the paper does not use this relation to predict independent data; it labels the quantity as 'Boltzmann-fitted soft part' and frames the result as consistency, not derivation. The statement that the total ⟨pT⟩ rise is due to a growing hard weight is a weighted-average bookkeeping identity once the decomposition is accepted, but the paper's central claim is only that the two-component picture is 'consistent with' and a 'viable alternative' to hydrodynamical interpretations, with soft/hard explicitly declared operational labels. The Pythia8 Monash comparison is an external simulation anchor for the fitting recipe. Reproducibility concerns remain—p0 values and hard-mean integration limits are not tabulated, and the abstract's promised with/without color-reconnection comparison does not appear in the body—but these are missing-support issues, not circular reductions of the result to its inputs. No equation-level equivalence or fitted-parameter-renamed-as-prediction was found.
Axiom & Free-Parameter Ledger
free parameters (3)
- β (inverse slope of Boltzmann fit)
- A (normalization)
- p0 (soft-hard boundary)
axioms (5)
- domain assumption Spectrum decomposes additively into soft exponential and hard fragmentation components.
- domain assumption There exists a pT range [0,p0] where the hard component is negligible.
- domain assumption A single exponential describes the soft component over the fitted interval.
- domain assumption χ²/ndf closest to one is a valid criterion for the cutoff.
- standard math Standard calculus of exponential integrals for ⟨pT⟩.
read the original abstract
We study the separation of soft and hard components in the transverse momentum spectra of charged particles as measured by ALICE in proton-proton collisions at $\sqrt{s}$ = 2.76, 5.02 and 13 TeV at the LHC. The soft component is described by a Boltzmann fit, while the residual spectra are identified as a hard QCD-like fragmentation contribution. After separation, the subtracted spectra show no significant evolution in shape or peak position with multiplicity, consistent with a two-component interpretation. Mean transverse momenta for both contributions remain nearly constant across multiplicity classes, while Pythia 8 Monte Carlo simulations confirm these trends. The robustness of the decomposition is demonstrated by a comparison between simulations with and without color reconnection, yielding consistent results. The terms `soft' and `hard' are used as operational labels within this framework and should not be interpreted as uniquely identified dynamical components. Our results are consistent with the two-component description as a viable and physically motivated alternative to hydrodynamical interpretations.
Figures
Reference graph
Works this paper leans on
-
[1]
L. Evans and P. Bryant, LHC Machine, JINST3, S08001. [2] S. Acharyaet al.(ALICE Collaboration), The AL- ICE experiment – A journey through QCD, (2022), 4 arXiv:2211.04384 [nucl-ex]
Pith/arXiv arXiv 2022
-
[3]
Strickland, Small system studies: A theory overview, Nucl
M. Strickland, Small system studies: A theory overview, Nucl. Phys. A982, 92 (2019)
2019
-
[4]
Ansariet al.(UA2 Collaboration), Jet Measures and Hadronic Event Shapes at the CERN ¯ppCollider, Z
R. Ansariet al.(UA2 Collaboration), Jet Measures and Hadronic Event Shapes at the CERN ¯ppCollider, Z. Phys. C36, 175 (1987)
1987
-
[5]
Wang and R
X.-n. Wang and R. C. Hwa, The Effect of Jet Production on the Multiplicity Dependence of Average Transverse Momentum, Phys. Rev. D39, 187 (1989)
1989
-
[6]
Ceradini (UA1 Collaboration), Study of minimum bias trigger events at 0.2 TeV to 0.9 TeV with magnetic and calorimetric analysis at the CERN proton-antiproton col- lider, Conf
F. Ceradini (UA1 Collaboration), Study of minimum bias trigger events at 0.2 TeV to 0.9 TeV with magnetic and calorimetric analysis at the CERN proton-antiproton col- lider, Conf. Proc. C850718, 363 (1985)
1985
-
[7]
T. A. Trainor, A two-component model for identified- hadronp t spectra from 5 TeV p-Pb collisions, J. Phys. G47, 045104 (2020), arXiv:1812.01151 [hep-ph]
Pith/arXiv arXiv 2020
-
[8]
T. A. Trainor, Comparing the PYTHIA Monte Carlo to a two-component (soft + hard) model of hadron produc- tion in high-energyp−pcollisions, Int. J. Mod. Phys. E 28, 1950023 (2019), arXiv:1805.09681 [hep-ph]
Pith/arXiv arXiv 2019
-
[9]
T. A. Trainor, How the blast-wave model describes PID hadron spectra from 5 TeV p-Pb collisions, SciPost Phys. Proc.15, 025 (2024), arXiv:2209.09474 [hep-ph]
Pith/arXiv arXiv 2024
-
[10]
Acostaet al.(CDF Collaboration), Soft and Hard Interactions inp¯pCollisions at √s= 1800-GeV and 630- GeV, Phys
D. Acostaet al.(CDF Collaboration), Soft and Hard Interactions inp¯pCollisions at √s= 1800-GeV and 630- GeV, Phys. Rev. D65, 072005 (2002)
2002
-
[11]
Kanetaet al.(NA44 Collaboration), Kaon and proton ratios from central Pb + Pb collisions at the CERN SPS, Nucl
M. Kanetaet al.(NA44 Collaboration), Kaon and proton ratios from central Pb + Pb collisions at the CERN SPS, Nucl. Phys. A638, 419 (1998)
1998
-
[12]
Roland (NA49 Collaboration), Recent results on cen- tral Pb+Pb collisions from experiment NA49, Nucl
G. Roland (NA49 Collaboration), Recent results on cen- tral Pb+Pb collisions from experiment NA49, Nucl. Phys. A654, 635c (1999)
1999
-
[13]
Andersenet al., Enhancement of central Lambda, Xi and Omega yields in Pb - Pb collisions at 158 A-GeV/c, Phys
E. Andersenet al., Enhancement of central Lambda, Xi and Omega yields in Pb - Pb collisions at 158 A-GeV/c, Phys. Lett. B433, 209 (1998)
1998
-
[14]
B. I. Abelevet al.(STAR Collaboration), Systematic Measurements of Identified Particle Spectra inpp, d+ Au and Au+Au Collisions from STAR, Phys. Rev. C79, 034909 (2009), arXiv:0808.2041 [nucl-ex]
Pith/arXiv arXiv 2009
-
[15]
Kisiel, Signatures of collective flow in high multi- plicity pp collisions, Phys
A. Kisiel, Signatures of collective flow in high multi- plicity pp collisions, Phys. Rev. C84, 044913 (2011), arXiv:1012.1517 [nucl-th]
Pith/arXiv arXiv 2011
-
[16]
W. Florkowski, Blast-wave model description of the HBT radii measured inppcollisions at the LHC energies, Acta Phys. Polon. B47, 2241 (2016), arXiv:1603.06418 [nucl- th]
Pith/arXiv arXiv 2016
-
[17]
B ´ ır´ o, L
G. B ´ ır´ o, L. Serkin, G. Pai´ c, and G. G. Barnaf¨ oldi, Inves- tigating the soft and hard limits in transverse momen- tum spectra in pp collisions, Eur. Phys. J. ST234, 2931 (2025)
2025
-
[18]
Tsallis, Possible Generalization of Boltzmann-Gibbs Statistics, J
C. Tsallis, Possible Generalization of Boltzmann-Gibbs Statistics, J. Statist. Phys.52, 479 (1988)
1988
-
[19]
J. Gu, C. Li, Q. Wang, W. Zhang, and H. Zheng, Collec- tive expansion in pp collisions using the Tsallis statistics, J. Phys. G49, 115101 (2022), arXiv:2201.02091 [nucl-th]
Pith/arXiv arXiv 2022
-
[20]
G. B ´ ır´ o, G. G. Barnaf¨ oldi, and T. S. Bir´ o, Tsallis- thermometer: a QGP indicator for large and small collisional systems, J. Phys. G47, 105002 (2020), arXiv:2003.03278 [hep-ph]
Pith/arXiv arXiv 2020
-
[21]
P.-P. Yang, M.-Y. Duan, F.-H. Liu, and R. Sahoo, Anal- ysis of Identified Particle Transverse Momentum Spec- tra Produced in pp, p–Pb and Pb–Pb Collisions at the LHC Using TP-like Function, Symmetry14, 1530 (2022), arXiv:2112.13223 [hep-ph]
Pith/arXiv arXiv 2022
-
[22]
S. Acharyaet al.(ALICE Collaboration), Multiplicity dependence ofπ, K, and p production in pp collisions at √s= 13 TeV, Eur. Phys. J. C80, 693 (2020), arXiv:2003.02394 [nucl-ex]
Pith/arXiv arXiv 2020
-
[23]
S. Acharyaet al.(ALICE Collaboration), Charged- particle production as a function of multiplicity and transverse spherocity in pp collisions at√s= 5.02 and 13 TeV, Eur. Phys. J. C79, 857 (2019), arXiv:1905.07208 [nucl-ex]
Pith/arXiv arXiv 2019
-
[24]
J. Adamet al.(ALICE Collaboration), Pseudorapidity and transverse-momentum distributions of charged par- ticles in proton–proton collisions at √s= 13 TeV, Phys. Lett. B753, 319 (2016), arXiv:1509.08734 [nucl-ex]
Pith/arXiv arXiv 2016
-
[25]
F. G. Gardim, A. V. Giannini, and J.-Y. Ollitrault, Accessing the speed of sound in relativistic ultracen- tral nucleus-nucleus collisions using the mean trans- verse momentum, Phys. Lett. B856, 138937 (2024), arXiv:2403.06052 [nucl-th]
Pith/arXiv arXiv 2024
-
[26]
T. Sj¨ ostrand, S. Ask, J. R. Christiansen, R. Corke, N. De- sai, P. Ilten, S. Mrenna, S. Prestel, C. O. Rasmussen, and P. Z. Skands, An introduction to PYTHIA 8.2, Comput. Phys. Commun.191, 159 (2015), arXiv:1410.3012 [hep- ph]
Pith/arXiv arXiv 2015
-
[27]
P. Skands, S. Carrazza, and J. Rojo, Tuning PYTHIA 8.1: the Monash 2013 Tune, Eur. Phys. J. C74, 3024 (2014), arXiv:1404.5630 [hep-ph]
Pith/arXiv arXiv 2013
-
[28]
A. Mazeliauskas and V. Vislavicius, Temperature and fluid velocity on the freeze-out surface fromπ,K,pspec- tra in pp, p-Pb and Pb-Pb collisions, Phys. Rev. C101, 014910 (2020), arXiv:1907.11059 [hep-ph]. 5
Pith/arXiv arXiv 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.