{"id":"5a39e58e-cd32-4dac-bea0-3bd966b62778","arxiv_id":"2412.05221","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The sqrt(s)-dependent intrinsic-kT width in PYTHIA 8 is traced to soft gluon emissions controlled by the ISR cutoff pT0Ref, with the fitted width increasing nearly linearly with pT0Ref.","lead":"By tuning PYTHIA 8 to Drell-Yan data at four collision energies, this paper shows that the fitted width of the intrinsic transverse momentum distribution grows with both collision energy and the parton shower's soft-gluon cutoff. The result suggests that the energy dependence long seen in Monte Carlo generators comes from non-perturbative soft gluon emissions, not from the partons' internal motion alone.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Attribution of the sigma(pT0Ref) trend to non-perturbative soft gluons is not established: pT0Ref smoothly suppresses all low-pT shower emissions, so fitted sigma may absorb perturbative shifts.","rationale":"The reader's weakest-assumption analysis and my stress-test pass converge on the same load-bearing point: the paper interprets the fitted Gaussian width sigma as cleanly separating intrinsic transverse motion from soft non-perturbative gluon emissions, with pT0Ref controlling only the latter. That separation is not demonstrated in the manuscript. In Pythia, pT0Ref enters as a smooth regulator of the ISR splitting probability, so varying it also changes the rate of emissions that are normally regarded as perturbative. The fitted sigma is therefore a global compensation parameter, and without a direct decomposition of which emitted partons drive the sigma shift, the causal attribution to non-perturbative soft gluons is under-supported. The absence of error bars in the key figure (Fig. 6) and the use of only three pT0Ref values per energy further weaken the specific linearity claim, though the qualitative upward trend is plausible. The paper is otherwise internally consistent and makes a useful observation about a longstanding Monte Carlo tuning issue; the independent support from the chi^2 fits and the comparison with data is real. My concern does not overturn the verdict: the conditional status remains appropriate, because the requested uncertainty propagation and causal decomposition would settle the strongest caveat. I recommend no change to the reader's CONDITIONAL verdict.","tokens_in":13176,"tokens_out":6057,"duration_ms":64803,"concrete_test":"Modify the Pythia ISR implementation to replace the smooth suppression factor pT^2/(pT^2+pT0^2) by a hard cutoff at pT0Ref, i.e., no suppression for emissions with pT > pT0Ref and zero splitting below pT0Ref, keeping all other settings of the intrinsic-kT-optimized tune fixed. Re-run the chi^2 extraction of sigma for pT0Ref = 0.5, 1, and 2 GeV at sqrt(s) = 13 TeV. If the fitted sigma values shift by more than the reported delta-chi^2 uncertainty, the default pT0Ref dependence is not a pure soft-gluon effect; the smooth suppression of the perturbative region contributes to the observed sigma(pT0Ref) trend.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central causal claim (Section 4.1, Conclusions) is that increasing pT0Ref removes non-perturbative soft gluon emissions, so the fitted intrinsic-kT width sigma grows to compensate. But in the Pythia shower the cutoff is implemented not as a hard separation between non-perturbative and perturbative emissions, but as the smooth factor pT^2/(pT^2+pT0^2) multiplying the splitting probability in the ISR evolution (Section 4, Ref [9]; see also Eq. (2) for the Sudakov form factor). Lowering pT0Ref from 2 to 0.5 GeV therefore changes the emission probability over the whole pT spectrum: even at pT = 2 GeV the suppression factor changes from 0.5 to 0.94. The paper provides no decomposition showing that the sigma shift is driven by emissions below a genuinely non-perturbative scale rather than by this smooth modification of the perturbative region. The comparison in Fig. 5(right) only shows that the fitted sigma is wrong at the other pT0Ref; it does not identify which emissions are responsible. Because an alpha_s adjustment was also needed to obtain the 'intrinsic kT optimized' tune (Section 2.2), other model parameters may be absorbing part of the pT0Ref effect. Figure 6 also shows no uncertainty bars (the sigma uncertainty from delta-chi^2 = 1 is quoted in Section 3 but not plotted), and the linear trend rests on only three pT0Ref points per energy.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper investigates the origin of the collision-energy dependence of the Gaussian intrinsic transverse-momentum width σ (PYTHIA 8 parameter BeamRemnants:primodialkThard) that is required when fitting low-pT Drell-Yan data. Using an 'intrinsic kT-optimized tune' (a Monash-2013-based setting with an adjusted αs for ISR/FSR), the authors fit σ to measurements at √s = 38.8, 200, 1960, and 13000 GeV for three values of the ISR cutoff SPACESHOWER:pT0Ref (0.5, 1.0, 2.0 GeV) via χ² minimization (Section 3, Eq. 1). Figure 6 shows that σ increases approximately linearly with pT0Ref, with a stronger slope at TeV-scale energies; the paper attributes the √s-dependence of σ to soft, non-perturbative gluon emissions controlled by the Sudakov form factor (Section 4, Eq. 2) and supports this with the observed mDY-dependence of σ (Section 4.2, Figs. 7–9). The central conclusion is that σ reflects the interplay of two non-perturbative processes: intrinsic parton motion and soft gluon emission.","tokens_in":13512,"tokens_out":9574,"duration_ms":93305,"significance":"If the central attribution holds, this paper provides a physical explanation for a long-standing generator-tuning puzzle: parton-shower generators such as PYTHIA and HERWIG require an intrinsic-kT width that grows with √s, whereas PB-TMD fits yield an essentially energy-independent σ. The study contains several useful falsifiable statements, including a small dσ/dpT0Ref slope at low energies with a growing slope at high energy, and a weak mDY-dependence of σ at 13 TeV with a visible dependence at low mDY. Credit is due for the transparent methodology (public data, MCatNLO+PYTHIA, Rivet, explicit χ² definition) and for the strong diagnostic in Fig. 5(right), where the σ tuned at pT0Ref = 0.5 GeV fails catastrophically at pT0Ref = 2 GeV (χ²/n = 70.14), demonstrating that the pT0Ref effect is large and not a small correction. The main weakness is that the Sudakov form factor is discussed qualitatively but never evaluated to make a quantitative prediction, and the key evidence (Fig. 6) is presented without uncertainties. The paper is therefore best read as a well-motivated interpretive study of a Monte-Carlo artifact rather than a quantitative derivation of the energy dependence.","major_comments":[{"comment":"The σ values in Fig. 6 are plotted without any uncertainty bars, even though Section 3 defines a Δχ² = 1 uncertainty on each tuned value (Fig. 4, left). Consequently the 'approximately linear' increase of σ with pT0Ref, and in particular the comparison of the slopes at 13 TeV and 1.96 TeV, cannot be quantitatively assessed. The problem is compounded by the fact that the fit at pT0Ref = 2 GeV for the CMS 13 TeV data has χ²/n = 2.34 (Fig. 5, left), which is not a good description even at the best-fit σ. Please provide the numerical values σ ± δσ for all points (as a table or as error bars in the figure) and explicitly comment on the reliability of the pT0Ref = 2 GeV extraction.","section":"Section 4.1 / Fig. 6"},{"comment":"The attribution of the σ(pT0Ref) trend specifically to non-perturbative soft gluon emission is not established by the presented analysis. In the PYTHIA shower, the pT0Ref cutoff enters not as a hard separation but as the smooth suppression factor pT²/(pT²+pT0²) multiplying the splitting probability at all transverse momenta (Section 4, Ref. [9]); lowering pT0Ref from 2 GeV to 0.5 GeV changes that factor from 0.5 to 0.94 even at pT = 2 GeV, so perturbative-region emissions are also modified. The manuscript provides no decomposition showing that the σ shift is driven by emissions below a genuinely non-perturbative scale. A concrete test would be to evaluate the Sudakov form factor of Eq. (2) with and without the pT0-suppression and compare the resulting low-pT(ℓℓ) predictions with the fitted σ values, thereby predicting σ(pT0Ref) and σ(√s) rather than reading them from the fits.","section":"Section 4.1 / Eq. (2)"},{"comment":"The definition of the 'intrinsic kT-optimized tune' is incomplete: the text states that 'by changing the value of αs(mZ) = 0.130 in the Monash tune, an excellent description is achieved. In the following, we take this value,' but the adopted αs value (and whether it applies to ISR, FSR, or both) is never specified. Since αs strongly affects the 2–5 GeV region shown in Fig. 2, and σ is subsequently fitted with αs held fixed, a possible degeneracy between σ and αs should be checked, for example by a two-dimensional χ² scan or by refitting αs at each pT0Ref. Without such a check, part of the σ shift in Fig. 6 could be absorbing the αs adjustment rather than the soft-gluon cutoff.","section":"Section 2.2"}],"minor_comments":[{"comment":"The caption reads '√s = 13 GeV' but should read '√s = 13 TeV' for the proton-proton collision energy.","section":"Figure 1 caption"},{"comment":"The parameter name 'primodialkThard' is misspelled throughout the manuscript (it should be 'primordialkThard'); please standardize the spelling.","section":"Section 3"},{"comment":"The text defines NDF as equal to the number of bins, which is non-standard (normally NDF = number of bins − number of fitted parameters); please clarify this convention or rename the quantity.","section":"Section 3"},{"comment":"The footnote contains a typo: 'higher order αs contribtuins' should be 'higher order αs contributions'.","section":"Section 2.2, footnote"},{"comment":"The sentence 'a significantly larger fraction of soft gluon contribution is removed when when the ISR cutoff parameter is increased' contains a duplicated 'when'.","section":"Section 4.1"},{"comment":"The phrase 'soft gluon effects become increasingly prominent with rising collision energy—contrary to initial expectations' is confusing, because an increasing role of soft emissions with increasing phase space is precisely the standard Sudakov expectation; please rephrase or justify what expectation is being contrasted.","section":"Abstract and Section 1"},{"comment":"The fit range for the χ² calculation is stated as '8 GeV when applicable', but the actual fitted pT(ℓℓ) ranges for the CDF, PHENIX, and E605 datasets are not given; please list them for each dataset.","section":"Section 3"},{"comment":"The description of the χ² scan does not state the number of σ-grid points or the step size used for the cubic spline interpolation; please provide these details for reproducibility.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a question of real interest to the MC-generator community, and the qualitative conclusion is likely correct, but the quantitative evidence is currently too weak: the central figure contains no uncertainties, the tune definition omits the αs value, and the attribution to non-perturbative soft gluons needs a sharper test. I would ask the authors to supply the numerical σ ± δσ table, the precise tune settings, and a Sudakov-based or decomposition-based check of the attribution; these are within the scope of a major revision. The citation pattern is reasonable, with the authors' previous PB-method work [1,2,26] used as a contrast rather than over-cited."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about arXiv:2412.05221. First, it pins down something concrete: in PYTHIA 8 the fitted Gaussian width sigma of the intrinsic-kT distribution rises almost linearly with the ISR cutoff pT0Ref, across four collision energies. That is a new observation, and it directly speaks to a known puzzle—why sigma appears to grow with sqrt(s) in parton-shower MCs. Second, the paper's interpretation of that trend, that sigma is compensating for suppressed non-perturbative soft gluon emissions, is plausible but not proven. The evidence supports the trend; the causal story is stronger than the evidence.\n\nThe paper does several things well. It builds an 'intrinsic kT-optimized tune' with pT0Ref = 0.5 GeV and shows it describes DY data from 38.8 GeV to 13 TeV with reasonable chi2. The fitting procedure is transparent: chi2 minimization in pT(ll), with the sigma uncertainty from delta-chi2 = 1 quoted. It also checks the mDY dependence and shows that the pT0Ref effect on sigma is larger at low mDY, which is qualitatively consistent with soft emissions at low scales. That is a nice cross-check.\n\nThe soft spots are real but not load-bearing. The stress-test concern is valid: pT0Ref enters as a smooth suppression factor pT^2/(pT^2+pT0^2), not a hard cut, so lowering it changes emission rates across the whole pT spectrum. The paper itself notes that 'soft emissions can appear at all scales Q2.' That means the sigma(pT0Ref) trend does not cleanly isolate non-perturbative soft gluons; it could partly absorb shifted perturbative emissions or other model deficiencies. The optimized tune also adjusts alpha_s, so the attribution is not unique. Figure 6 has no error bars on sigma, and the linear trend uses only three points per energy. These are presentation and evidence issues, not fatal flaws: the qualitative pattern—steeper sigma(pT0Ref) at higher sqrt(s), flatter at low energy—would likely survive better uncertainty handling.\n\nWho benefits: people who tune PYTHIA or interpret intrinsic-kT in TMD and DY contexts. It is a practical, reproducible study, but not a first-principles derivation. It deserves a serious referee. If I were handling it, I would send it out and ask for explicit error bars in Fig. 6, a test of whether the sigma shift is driven by emissions below a genuine non-perturbative scale, and a discussion of how much the alpha_s adjustment alters the conclusion. My bottom line: engage with this paper, but treat the soft-gluon claim as a working hypothesis rather than a settled result.","headline":"The paper pins down a clear empirical trend—sigma rises with pT0Ref—but the soft-gluon interpretation is plausible rather than proven; send to review.","tokens_in":727,"tokens_out":801,"would_cite":true,"duration_ms":33313,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The energy growth of the intrinsic-$k_T$ width in PYTHIA is caused by soft gluon emissions, not by a change in partons' internal motion.","keywords":["intrinsic transverse momentum","Drell-Yan production","parton shower","soft gluon emissions","Sudakov form factor","PYTHIA","transverse momentum spectrum","non-perturbative QCD"],"falsifier":"Extract the intrinsic-kT width from low-pT Drell-Yan data using a TMD calculation whose Sudakov form factor is computed analytically (rather than through a shower cutoff), and compare the fitted sigma across sqrt(s) = 38.8 GeV, 200 GeV, 1.96 TeV and 13 TeV. The paper's claim implies the TMD-extracted width is roughly energy-independent; if it grows with energy, the claim is falsified.","tokens_in":12959,"feed_emoji":"⚛️","tokens_out":8776,"duration_ms":77678,"temperature":0.7,"pith_summary":"This paper asks why the intrinsic transverse momentum width fitted in PYTHIA grows with collision energy, reaching values that seem unphysical. The answer it argues for is that the width is not a clean measure of partons' internal motion: it also absorbs non-perturbative soft gluon emissions. Using Drell-Yan pair production at low transverse momentum, where these effects dominate, the authors fit the Gaussian width sigma at several center-of-mass energies and values of the initial-state-radiation cutoff pT0Ref. They find that sigma increases approximately linearly with pT0Ref in the range 0.5-2.0 GeV, so the energy dependence enters through the Sudakov no-emission probability rather than through intrinsic motion. If correct, this resolves a longstanding puzzle in parton-shower modeling and tells generators how to separate the two non-perturbative contributions.","feed_headline":"Soft gluons, not internal motion, explain PYTHIA's energy puzzle","feed_subtitle":"Fits show the intrinsic-kT width grows with the shower cutoff, tracing its energy dependence to soft gluons.","key_machinery":"The load-bearing object is the PYTHIA initial-state parton shower with its non-perturbative regulator pT0Ref, which suppresses small-transverse-momentum emissions by a factor $pT^{2}$/($pT^{2}$ + $pT0Ref^{2}$). The Sudakov form factor, the probability of no radiation between two scales, depends on this cutoff, so changing pT0Ref changes how many soft gluons are emitted and thus how much transverse momentum the intrinsic Gaussian distribution must provide to describe the measured low-pT Drell-Yan spectrum. The argument is carried by the fitted $\\sigma$ as a function of pT0Ref and $\\sqrt$(s), together with the dependence of $\\sigma$ on the DY invariant mass, which probes the scale dependence of the Sudakov form factor.","core_discovery":"The central claim is that the sqrt(s)-dependence of the intrinsic-kT width sigma observed in PYTHIA 8 arises from the interplay between two non-perturbative processes: the internal transverse motion of partons and soft gluon emissions. By fitting sigma to measured Drell-Yan transverse momentum spectra at sqrt(s) = 38.8 GeV, 200 GeV, 1.96 TeV and 13 TeV, and repeating the fit for three values of the ISR cutoff pT0Ref (0.5, 1.0 and 2.0 GeV), the paper shows that at fixed energy sigma rises nearly linearly with pT0Ref, and at fixed pT0Ref sigma rises with energy. Since lowering pT0Ref unsuppresses soft gluon emissions, the fitted width must be understood as a combined observable: the part of the low-pT spectrum that the shower does not generate through soft radiation is compensated by a larger intrinsic Gaussian width. The paper therefore concludes that the energy scaling in PYTHIA is a compensation mechanism for soft gluon emissions, not a property of intrinsic transverse momentum.","pith_inferences":["If the paper is right, other event generators with different soft-gluon regulators should show a different sigma-versus-energy pattern, so comparing the same data across generators would test whether the effect is universal or regulator-specific.","A practical extension would be to replace the Gaussian intrinsic-kT model with a TMD distribution that already contains the Sudakov suppression, in which case the fitted sigma should become independent of pT0Ref; this is a testable prediction of the paper's interpretation.","Because sigma absorbs soft-gluon effects, the unphysical high-energy widths are not a failure of the generator but a sign that the Gaussian ansatz is being pushed to cover two distinct physics effects; a two-component fit separating both contributions would be more physical.","The near-linear sigma-pT0Ref relation could serve as a tuning rule: once pT0Ref is chosen, sigma is determined by the requirement that the total low-pT spectrum matches data, turning a puzzling energy dependence into a known compensation."],"forward_implications":["Fitted intrinsic-kT widths from parton-shower generators should not be interpreted as direct measurements of intrinsic transverse momentum; they are effective parameters that also encode soft-gluon physics.","At fixed collision energy, sigma can be changed by nearly 1 GeV by moving pT0Ref from 0.5 to 2.0 GeV, so the low-pT DY spectrum alone does not fix both parameters; the two must be tuned together.","The energy dependence reported by previous studies disappears once the soft-gluon contribution is accounted for, meaning intrinsic transverse motion could be roughly energy-independent as the parton-branching method suggests.","Soft gluon effects grow with collision energy and are largest at low DY mass, so low-mass Drell-Yan measurements are the most sensitive place to constrain the non-perturbative Sudakov form factor."],"supporting_citations":[{"why":"The work showing that with the parton-branching method the intrinsic-kT width is nearly sqrt(s)-independent, providing the baseline the paper must explain.","marker":"[1]"},{"why":"The study showing that excluding the non-perturbative region in the parton-branching method reproduces a sqrt(s)-dependence similar to parton-shower generators.","marker":"[2]"},{"why":"The CMS measurement confirming the energy-dependent intrinsic-kT width in PYTHIA and HERWIG tunes, the observation this paper explains.","marker":"[3]"},{"why":"The description of the PYTHIA 8.3 initial-state shower whose pT0Ref suppression mechanism is varied throughout the analysis.","marker":"[9]"},{"why":"The default tune used as the baseline for the simulations before the intrinsic-kT-optimized tune is defined.","marker":"[17]"},{"why":"The CMS 13 TeV Drell-Yan transverse momentum data used to fit sigma at the highest energy.","marker":"[18]"},{"why":"The PHENIX 200 GeV Drell-Yan cross section used for the intermediate-energy fit.","marker":"[21]"},{"why":"The E605 38.8 GeV dimuon data used for the lowest-energy fit.","marker":"[22]"},{"why":"The CDF 1.96 TeV Z-region data used for the Tevatron-energy fit.","marker":"[25]"}],"fun_headline_variants":["Soft gluons, not motion, drive PYTHIA's kT energy scaling","PYTHIA's kT energy puzzle is soft gluon compensation","Soft gluons, not intrinsic motion, set PYTHIA's kT width","Drell-Yan data trace PYTHIA's kT width scaling to soft gluons","Energy scaling of PYTHIA's kT width comes from soft gluons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that changing PYTHIA's pT0Ref parameter alters only the non-perturbative soft-gluon part of the shower, so that the fitted Gaussian width sigma cleanly separates intrinsic motion from soft-gluon effects and does not absorb other model deficiencies.","fun_headline_variants_meta":{"raw":{"variants":["Soft gluons, not motion, drive PYTHIA's kT energy scaling","PYTHIA's kT energy puzzle is soft gluon compensation","Soft gluons, not intrinsic motion, set PYTHIA's kT width","Drell-Yan data trace PYTHIA's kT width scaling to soft gluons","Energy scaling of PYTHIA's kT width comes from soft gluons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00152,"raw_usage":{"total_tokens":6104,"prompt_tokens":977,"completion_tokens":5127,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":5023}},"tokens_in":593,"tokens_out":5127,"duration_ms":31615,"temperature":1.0,"reasoning_tokens":5023,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T20:48:38.645565+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Extract the intrinsic-kT width from low-pT Drell-Yan data using a TMD calculation whose Sudakov form factor is computed analytically (rather than through a shower cutoff), and compare the fitted sigma across sqrt(s) = 38.8 GeV, 200 GeV, 1.96 TeV and 13 TeV. The paper's claim implies the TMD-extracted width is roughly energy-independent; if it grows with energy, the claim is falsified.","supporting_citations":[{"cited_title":"Dimuon production in proton - copper collisions at √s = 38.8 GeV","cited_arxiv_id":null,"evidence_quote":"The E605 38.8 GeV dimuon data used for the lowest-energy fit."}],"review_version":1}