Pith. sign in

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 →

arxiv 2510.09692 v2 pith:YO2SQEKC submitted 2025-10-09 hep-ph nucl-th

Revisiting the soft-hard separation in the transverse momentum spectra of pp collisions

classification hep-ph nucl-th
keywords transverse momentum spectrasoft-hard separationtwo-component modelBoltzmann fitproton-proton collisionsmultiplicity dependencemean transverse momentumfragmentation
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper argues that the transverse momentum spectra of charged particles from proton-proton collisions can be cleanly divided into a soft piece (an exponential or Boltzmann-like form) and a hard piece (a QCD-like fragmentation tail). Once this division is made, the mean transverse momentum of each component stays nearly constant as the charged-particle multiplicity increases, while the mean transverse momentum of the total spectrum rises. The rise is therefore not caused by a change in the soft sector but by the growing weight of hard fragmentation. The paper presents this as support for a two-component interpretation of the spectra, an alternative to hydrodynamic or blast-wave descriptions, and reproduces the trends with a standard Monte Carlo event generator. The authors stress that 'soft' and 'hard' are operational labels, not uniquely identified dynamical components.

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.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 5 minor

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)
  1. [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.
  2. [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.
  3. [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.
  4. [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)
  1. [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.
  2. [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.
  3. [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.
  4. [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.
  5. [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

0 steps flagged

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

3 free parameters · 5 axioms · 0 invented entities

The paper introduces no new entities, but it relies on the existence and identifiability of a purely soft low-pT region, a single-exponential shape, and a χ²-based cut. These are domain assumptions that carry the central interpretation.

free parameters (3)
  • β (inverse slope of Boltzmann fit)
    Fitted separately for each multiplicity class; the paper's claim that the soft component has stable ⟨pT⟩ is essentially a statement that β is stable.
  • A (normalization)
    Fitted per spectrum; scale only.
  • p0 (soft-hard boundary)
    Chosen per spectrum as min χ²/ndf; acts as a free parameter defining the two components; no systematic variation is shown.
axioms (5)
  • domain assumption Spectrum decomposes additively into soft exponential and hard fragmentation components.
    Introduced in Sec. I-II as the premise to be tested; never independently verified.
  • domain assumption There exists a pT range [0,p0] where the hard component is negligible.
    Sec. II: p0 'marks the region where the contribution from hard QCD processes is still negligible.' Load-bearing for the inference that the residual spectrum is fragmentation.
  • domain assumption A single exponential describes the soft component over the fitted interval.
    Sec. II: 'Our goal here was not to assign a physical temperature but to test how well the Boltzmann form can represent the data.'
  • domain assumption χ²/ndf closest to one is a valid criterion for the cutoff.
    Used in Sec. II Fig. 1 to define p0; this choice is data-driven and affects all downstream results.
  • standard math Standard calculus of exponential integrals for ⟨pT⟩.
    Used implicitly when converting fitted β to ⟨pT⟩.

pith-pipeline@v1.3.0-alltime-deepseek · 5830 in / 13899 out tokens · 124857 ms · 2026-08-04T10:47:36.608579+00:00 · methodology

0 comments
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

Figures reproduced from arXiv: 2510.09692 by G\'abor B\'ir\'o, Gergely G\'abor Barnaf\"oldi, Guy Pai\'c, Leonid Serkin.

Figure 1
Figure 1. Figure 1: FIG. 1. The [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Subtracted hadron spectra obtained by removing the [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Boltzmann-distribution fits shown on linear scale at [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The mean transverse momenta, [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

27 extracted references · 17 linked inside Pith

  1. [1]

    Evans and P

    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]

  2. [3]

    Strickland, Small system studies: A theory overview, Nucl

    M. Strickland, Small system studies: A theory overview, Nucl. Phys. A982, 92 (2019)

  3. [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)

  4. [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)

  5. [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)

  6. [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]

  7. [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]

  8. [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]

  9. [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)

  10. [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)

  11. [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)

  12. [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)

  13. [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]

  14. [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]

  15. [16]

    Florkowski, Blast-wave model description of the HBT radii measured inppcollisions at the LHC energies, Acta Phys

    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]

  16. [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)

  17. [18]

    Tsallis, Possible Generalization of Boltzmann-Gibbs Statistics, J

    C. Tsallis, Possible Generalization of Boltzmann-Gibbs Statistics, J. Statist. Phys.52, 479 (1988)

  18. [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]

  19. [20]

    B ´ ır´ o, G

    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]

  20. [21]

    Yang, M.-Y

    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]

  21. [22]

    Acharyaet al.(ALICE Collaboration), Multiplicity dependence ofπ, K, and p production in pp collisions at √s= 13 TeV, Eur

    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]

  22. [23]

    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

    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]

  23. [24]

    Adamet al.(ALICE Collaboration), Pseudorapidity and transverse-momentum distributions of charged par- ticles in proton–proton collisions at √s= 13 TeV, Phys

    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]

  24. [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]

  25. [26]

    Sj¨ ostrand, S

    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]

  26. [27]

    Skands, S

    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]

  27. [28]

    Mazeliauskas and V

    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