pith. sign in

arxiv: 2604.18743 · v1 · submitted 2026-04-20 · ✦ hep-ph

Multiparticle production in electron-positron annihilation

E. S. Kokoulina This is my paper

Pith reviewed 2026-05-10 03:41 UTC · model grok-4.3

classification ✦ hep-ph
keywords multiparticle productionelectron-positron annihilationgluon dominance modelquark-gluon cascadehadronizationmultiplicity distributionQCD soft processes
0
0 comments X

The pith

The gluon dominance model describes multiplicity in electron-positron annihilation by convoluting quark-gluon cascades with hadronization tuned from hadron collisions.

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

The paper reanalyzes multiplicity data from earlier electron-positron annihilation experiments in light of new hadronization results from proton and heavy-ion collisions. It presents the gluon dominance model as a framework that combines a perturbative Markovian quark-gluon cascade with a subsequent hadronization stage whose parameters are fixed by experimental fits. This convolution is claimed to account for observed multiplicities across both lepton and hadron high-energy processes. A reader would care because the approach offers a consistent alternative to Monte Carlo generators that often mismatch data in soft regimes. The review concludes that fresh hadron-collision insights now justify revisiting the older e+e- measurements.

Core claim

Convolution of the qg-cascade and hadronization in the gluon dominance model describes the multiplicity in practically all processes of multiple production in both lepton and hadron high-energy collisions, making a reanalysis of e+e- annihilation data necessary after new hadronization findings from other reactions.

What carries the argument

The gluon dominance model, which treats multiparticle production as a Markovian branching quark-gluon cascade followed by a phenomenological hadronization stage selected from data.

If this is right

  • Multiplicity distributions in e+e- annihilation can be calculated using the same hadronization parameters extracted from proton and heavy-ion data.
  • The model supplies a unified description of multiplicities for lepton and hadron collisions without separate tuning for each initial state.
  • Soft processes in high-energy interactions can be handled by this convolution rather than relying exclusively on Monte Carlo event generators.
  • Updated hadronization knowledge from current collider programs can be imported to refine predictions for older and future annihilation measurements.

Where Pith is reading between the lines

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

  • A universal hadronization mechanism may operate across different collision initial states once the cascade stage is properly separated.
  • The framework could be tested by predicting multiplicities in new e+e- runs at facilities like Belle II or future colliders using only parameters from LHC heavy-ion data.
  • It suggests that discrepancies between Monte Carlo generators and data often stem from incomplete modeling of the hadronization transition rather than the perturbative cascade itself.

Load-bearing premise

The hadronization stage with parameters chosen from hadron collision data remains unchanged in its essential features when applied back to electron-positron annihilation.

What would settle it

A calculation using the gluon dominance model with hadron-collision-tuned hadronization parameters that deviates substantially from measured multiplicity distributions in e+e- annihilation at several center-of-mass energies would falsify the model’s applicability.

Figures

Figures reproduced from arXiv: 2604.18743 by E. S. Kokoulina.

Figure 1
Figure 1. Figure 1: Experimental values of f2 as a function of the mean multiplicity of negative particles n− for annihilation and non-annihilation processes [17] [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The description of multiplicity distributions, [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The description of multiplicity distributions, [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The description of multiplicity distributions, [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Parameters of GDM, kp (on the left) and m (on the right). With following increasing of energy it takes on a value that already exceeds one (∼ 1.2). Such behavior we observed in hadron interactions [9]. Mechanism of hadronization in pp collisions is recombination. In this case it occurs in dense quark-gluon medium. From this [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Parameters of hadronization for a quark jet: [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Parameters of hadronization for a gluon jet: [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The GDM parameter, α (on the left), the description of MD in e +e − annihilation by GDM and data at 9.4 GeV [30] (on the right) [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The description of linear growth with energy of a parameter [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The description of average gluon multiplicity, [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
read the original abstract

Multiparticle production in hadron and lepton interactions still attracts our attention. Simulation by using Monte Carlo event generators is performed before planning any experiment. But it often overestimates (or underestimates) experimental data. These generators are based on the theory of strong interactions, quantum chromodynamics (QCD), which is capable of performing calculations only in the perturbation theory. Soft processes that make up a significant contribution in high-energy interactions are forced to involve phenomenological models. Of all multiparticle production processes, electron-positron annihilation is the theoretically cleanest, proceeding via an intermediate virtual photon or $Z^0$-boson followed by quark-antiquark pair creation. QCD describes well the development of quark-gluon ($qg$) cascade as marcovian branching process, that is called first stage. The transformation of quarks and gluons produced in the $qg$-cascade into observable hadrons occurs in the second stage, hadronization, to which perturbation theory is no longer applicable. The choice of a scheme for it is based on experimental data. Convolution of $qg$-cascade and hadronization allowed us to describe the multiplicity in practice all processes of multiple production in both lepton and hadron high-energy collisions. This model is called the gluon dominance model. Several decades have passed since a series of $e^+e^-$ annihilation experiments were carried out. Now, the main interests of high energy physicists are focused on the study of multiparticle production in proton and heavy ion collisions. Their research revealed many new results in the theory of strong interactions, including the hadronization. That is why it appeared necessary to analyze multiplicity n $e^+e^-$-annihilation again.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 3 minor

Summary. The manuscript revisits the gluon dominance model for multiparticle production in electron-positron annihilation. The process is described as a two-stage convolution: a perturbative Markovian quark-gluon cascade calculable in QCD, followed by a non-perturbative hadronization stage whose scheme is selected from experimental data. The authors claim this convolution successfully describes multiplicities across lepton and hadron collisions and argue that recent insights from proton and heavy-ion collisions require a re-analysis of older e+e- data.

Significance. If the model can be shown to yield stable, non-circular predictions with explicit parameter comparisons and quantitative fits to data, it would offer a useful phenomenological framework unifying cascade and hadronization across collision types, potentially aiding Monte Carlo generators. The identification of e+e- annihilation as a theoretically clean testbed is appropriate, but the absence of new calculations or invariance tests in the manuscript limits its immediate contribution.

major comments (3)
  1. [Abstract] Abstract: The central claim that the convolution 'allowed us to describe the multiplicity in practice all processes of multiple production in both lepton and hadron high-energy collisions' is unsupported by any equations for the cascade or hadronization, data tables, fits, chi-squared values, or error analysis. This renders the success assertion unverifiable from the text.
  2. [Motivation for re-analysis] Motivation paragraph: The paper asserts that new hadron-collision results on hadronization make re-analysis of e+e- multiplicities necessary, yet provides no re-fits, before/after parameter comparisons, updated multiplicity predictions, or invariance tests for the hadronization parameters when new insights are incorporated. Without such demonstration, the re-analysis motivation is unsubstantiated and the prior success claim cannot be preserved.
  3. [Hadronization stage] Hadronization stage description: The scheme is explicitly 'chosen based on experimental data' while the overall model is said to describe 'practically all processes.' If the hadronization parameters were tuned on the same or overlapping e+e- and hadron datasets, the description risks circularity; the manuscript supplies no test separating the cascade prediction from the fitted hadronization component.
minor comments (3)
  1. [Abstract] Typo: 'marcovian' should read 'Markovian'.
  2. [Abstract] Typo: 'practice all' should read 'practically all'.
  3. [Abstract] The phrase 'multiplicity n $e^+e^-$' is incomplete or mistyped; it appears to intend 'multiplicity in e+e- annihilation'.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, indicating where revisions will be made to improve clarity, verifiability, and substantiation while preserving the note's concise scope as a motivation for re-analysis.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim that the convolution 'allowed us to describe the multiplicity in practice all processes of multiple production in both lepton and hadron high-energy collisions' is unsupported by any equations for the cascade or hadronization, data tables, fits, chi-squared values, or error analysis. This renders the success assertion unverifiable from the text.

    Authors: We acknowledge that the abstract makes a broad claim without supporting details in this short note. The gluon dominance model, including the Markovian qg-cascade equations from perturbative QCD and the hadronization convolution, along with quantitative fits, data comparisons, and error analyses, are fully detailed in our prior publications. In the revised version, we will modify the abstract to reference those specific earlier works where the equations, tables, and statistical measures are presented, thereby making the claim traceable within the current text. revision: yes

  2. Referee: [Motivation for re-analysis] Motivation paragraph: The paper asserts that new hadron-collision results on hadronization make re-analysis of e+e- multiplicities necessary, yet provides no re-fits, before/after parameter comparisons, updated multiplicity predictions, or invariance tests for the hadronization parameters when new insights are incorporated. Without such demonstration, the re-analysis motivation is unsubstantiated and the prior success claim cannot be preserved.

    Authors: The motivation arises from qualitative new insights on hadronization obtained in recent proton and heavy-ion studies. We agree that the manuscript does not yet demonstrate explicit re-fits or parameter shifts. In revision, we will expand the motivation paragraph to cite the specific recent hadron-collision results and outline their expected impact on e+e- multiplicity analysis. A full quantitative re-fit with before/after comparisons and invariance tests lies beyond the scope of this concise note and is reserved for a follow-up study. revision: partial

  3. Referee: [Hadronization stage] Hadronization stage description: The scheme is explicitly 'chosen based on experimental data' while the overall model is said to describe 'practically all processes.' If the hadronization parameters were tuned on the same or overlapping e+e- and hadron datasets, the description risks circularity; the manuscript supplies no test separating the cascade prediction from the fitted hadronization component.

    Authors: The first stage (qg-cascade) is computed independently as a perturbative Markovian branching process using QCD. Hadronization parameters are fixed from data but applied uniformly across collision types to test universality, as shown in our earlier cross-validation studies. To mitigate the circularity concern, we will insert a clarifying paragraph that explicitly distinguishes the perturbative cascade calculation from the non-perturbative hadronization stage and references the separation tests performed in prior work. revision: yes

Circularity Check

1 steps flagged

Hadronization scheme selected from data enables claimed description of multiplicities, but no explicit reduction of new e+e- results to prior fits shown.

specific steps
  1. fitted input called prediction [Abstract]
    "The choice of a scheme for it is based on experimental data. Convolution of qg-cascade and hadronization allowed us to describe the multiplicity in practice all processes of multiple production in both lepton and hadron high-energy collisions. This model is called the gluon dominance model."

    The hadronization scheme is selected from experimental data precisely so that the convolution reproduces observed multiplicities. The statement that the model 'allowed us to describe' the data across processes therefore reduces to the fitting step by construction; the description is not an independent output of the cascade alone.

full rationale

The paper's core assertion is that the gluon dominance model (qg-cascade convolved with hadronization) describes multiplicities across lepton and hadron collisions. The abstract explicitly states the hadronization scheme is chosen based on experimental data and that this convolution 'allowed us to describe' the multiplicities. This is a standard phenomenological construction rather than a first-principles derivation, but the text provides no equations, no parameter tables, and no demonstration that the re-analysis of e+e- data uses unchanged parameters or performs an invariance test against the new hadron-collision insights. The central claim therefore rests on a data-tuned stage without shown independence from the inputs it is said to describe. No self-citation chain or uniqueness theorem is invoked in the supplied text, so the circularity is partial and limited to the descriptive success claim.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on two stages whose interface is not derived from first principles: perturbative QCD cascade (standard) and a data-driven hadronization scheme (phenomenological). No new entities are postulated.

free parameters (1)
  • hadronization parameters
    Scheme for transforming quarks and gluons into hadrons is chosen based on experimental data; specific values not given in abstract but required for the convolution.
axioms (2)
  • domain assumption QCD describes the qg-cascade as a Markovian branching process
    Invoked in the description of the first stage; standard in perturbative QCD but assumes the branching approximation holds down to the hadronization scale.
  • ad hoc to paper Hadronization can be modeled by a convolution with the cascade output
    The model definition itself; no derivation from QCD is provided.

pith-pipeline@v0.9.0 · 5597 in / 1469 out tokens · 22399 ms · 2026-05-10T03:41:16.784151+00:00 · methodology

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

31 extracted references · 31 canonical work pages · 1 internal anchor

  1. [1]

    Two elementary events are taken into account: q-bremsstrahlung (q→q+g) and gluon fission (g→g+g)

    is taken as an evolutionary parameter. Two elementary events are taken into account: q-bremsstrahlung (q→q+g) and gluon fission (g→g+g). The probabilities of these events are determined in QCD (PT). The formation of aqq-pair from a gluon at the qg-cascade stage is suppressed according to QCD estimates, so it is neglected. Differential-difference equations...

    work page

  2. [2]

    Yndurain

    F.J. Yndurain. Quantum Chromosynamics. An introduction to the theory of quarks and gluons. Springer-Verlag. 1983

    work page 1983

  3. [3]

    Yu. L. Dokshitzer, D. I. D’yakonov, and S. I. Troyan. Hard semi-inclusive processes in QCD. Physics Letters. B78(1978) 290

    work page 1978

  4. [4]

    E. Fermi. Multiparticle production of pions in nucleon-nucleon collisions at cosmotron- energies. Physical Review.92(1953) 452

    work page 1953

  5. [5]

    Hagedorn

    R. Hagedorn. Statistical thermodynamics of strong interactions at high-energies. Nuovo Cimento Supplimento.3(1965)

    work page 1965

  6. [6]

    Henley and Alejandro Garcia

    Ernst M. Henley and Alejandro Garcia. Subatomic physics. World Scientific. (2007) 620 p

    work page 2007

  7. [7]

    Kokoulina

    E.C. Kokoulina. Analysis of multiparticle dynamics ine+e− - annihilation into hadrons by two-stage model. The proceeding of XXXII ISMD, Alushta, Ukraine. 2002. World Scientific. 340

    work page 2002

  8. [8]

    E. S. Kokoulina. Multiplicity distribution studies ofe+e−-annihilation at 50-61.4 GeV and 172-189 GeV by two stage model hadronization. 11th Annual Seminar on Nonlinear Phe- nomena in Complex Systems: Chaos, Fractals, Phase Transitions, Selforganization. Minsk, Belarus, 13-16 May 2002. Institute of Physics, Minsk.11(2003) 123

    work page 2002

  9. [9]

    Kokoulina

    E. Kokoulina. Description of pp-interactions with very high multiplicity at 70 GeV/c. Acta Phys. Polon.35(2004) 295

    work page 2004

  10. [10]

    Kokoulina

    E.S. Kokoulina. High multiplicity study and gluon dominance model. Phys.Part.Nucl.Lett. 13(2016) 74

    work page 2016

  11. [11]

    Kokoulina

    E. Kokoulina. and A. Kutov. Gluon dominance model. European Phys. J. Web of Conf.138 (2017) 07007

    work page 2017

  12. [12]

    Riadovikov

    V.N. Riadovikov. Topological cross sections in proton-proton interactions at 50 GeV. Physics of Atomic Nuclei.75(2012) 315

    work page 2012

  13. [13]

    Ammosov, et al

    V.V. Ammosov, et al. Average charged particle multiplicity and topological cross sections in 5-GeV/c and 69-GeV/c pp interactions Physics LettersB42(1972) 519

    work page 1972

  14. [14]

    V. N. Ryadovikov.Fluctuations of the Number of Neutral Pions at High Multiplicity in pp Interactions at 50 GeV. Physics of Atomic Nuclei.75(2012) 1050

    work page 2012

  15. [15]

    Barlykov, V

    N. Barlykov, V. Dudin, E. Kokoulina, A.Kutov, V. Nikitin, V. Popov, R. Shulyakovsky. Collective Phenomena in pp Interactions with High Multiplicity. Nonlinear Dynamics and Applications.28(2022) 169

    work page 2022

  16. [16]

    Konishi, A

    K. Konishi, A. Ukawa, G. Veneziano. A simple algorithm for resolving QCD jets. Nuclear Physics. B157(1979) 45

    work page 1979

  17. [17]

    Giovannini

    A. Giovannini. QCD jets as Markov branching processes. Nuclear Physics. B¯161(1979) 429. 14

    work page 1979

  18. [18]

    Rushbrooke and B.R

    J.G. Rushbrooke and B.R. Webber. High energy antiparticle-particle reaction differences and annihilations. Physics Reports.44(1978) 1

    work page 1978

  19. [19]

    New Discoveries at RHIC

    B. Muller. Hadronic signals of deconfinement at RHIC. Contribution to RBRC Scientific Articles Proceedings Series "New Discoveries at RHIC". Nuclear Physics A750(2005)2005 84

    work page 2005

  20. [20]

    Kokoulina and V.I

    E.S. Kokoulina and V.I. Kuvshinov. Hadron multiplicity distribution in the three gluon decay of quarkonia. Sov.Phys.J. 28 (1985) 757-760

    work page 1985

  21. [21]

    Kokoulina, Y

    E. Kokoulina, Y. Caballero Duran, T. Gorelkina, A. Kutov, V. Nikitin, V. Popov, K. Sarkar, E. Shakhvorostova, and Y. Shousha. Annihilation in the Gluon Dominance Model. Nonlinear Phenomena in Complex Systems.27(2024) 246

    work page 2024

  22. [22]

    Abramov et al

    V.V. Abramov et al. Possible Studies at the First Stage of the NICA Collider Operation with Polarized and Unpolarized Proton and Deuteron Beams. SPD Collab. Phys.Part.Nucl. 52(2021) 1044

    work page 2021

  23. [23]

    Page, Xiaouxuan Chu, E.C

    B.S. Page, Xiaouxuan Chu, E.C. Aschenauer. Experimental aspects of jet physics at a future EIC. Phys. Rev. D101(2020) 072003

    work page 2020

  24. [24]

    Althoff et al

    M. Althoff et al. TASSO Coll. Jet production and fragmentation ine+e− annihilation at 12-43 GeV. Z.Physics22(1984) 307

    work page 1984

  25. [25]

    H. W. Zgeng et al. AMY Coll. Charged particle multiplicities ine+e− annihilation at 50 - 61.4 GeV. Phys. Rev.42(1990) 737

    work page 1990

  26. [26]

    P. D. Acton et al. OPAL Coll. A study of charged particle multiplicities in hadronic decays of theZ 0. Z. Physics.C52(1991) 271

    work page 1991

  27. [27]

    Alexander et al

    G. Alexander et al. OPAL Coll. QCD studies withe+e− annihilation data at 130 and 136 GeV. Z. Physics.72(1996) 191

    work page 1996

  28. [28]

    Ackerstaff et al

    K. Ackerstaff et al. OPAL Coll. QCD studies withe+e− annihilation at 161 GeV. Z. Physics. 75(1997) 193

    work page 1997

  29. [29]

    Abiendi et al

    G. Abiendi et al. OPAL Coll. QCD studies withe +e− annihilation data at 172-189 GeV.Eur.Phys. J.C16(2000) 185

    work page 2000

  30. [30]

    Kokoulina,V

    E. Kokoulina,V. Nikitin. The description of pp-interactions with very high multiplicity at 70 GeV/c by Two Stage Gluon Model. Proceeding of the International School-Seminar Actual Problems of Microworld Physics, Gomel, Belarus, 2003. arXiv:hep-ph/0308139 [hep-ph]

    work page internal anchor Pith review arXiv 2003

  31. [31]

    Berger etal

    Ch. Berger etal. PLUTO Coll. Multiplicity distributions ine+e− annihilation at PETRA energies. Phys. Lett.B95(1980) 313

    work page 1980