Multiparticle production in electron-positron annihilation

This paper revisits the analysis of multiplicity in electron-positron annihilation to refine the Gluon Dominance Model, which combines perturbative QCD quark-gluon cascades with phenomenological hadronization schemes to better describe multiparticle production in high-energy collisions.

The Gist

Imagine you are watching a high-speed collision between two tiny particles, an electron and a positron (the electron's antimatter twin). When they smash into each other, they vanish and turn into pure energy, which then instantly explodes into a shower of new particles. This is called multiparticle production.

Physicists have been trying to predict exactly how many particles come out of this explosion and what they look like for decades. They use complex computer programs (called Monte Carlo generators) to simulate these crashes before building real machines. But often, these simulations get the numbers wrong, especially when a huge number of particles are created.

This paper by E.S. Kokoulina proposes a new, simpler way to understand this explosion using a model called the Gluon Dominance Model (GDM). Here is the breakdown using everyday analogies:

1. The Two-Stage Explosion

The author suggests that this particle explosion happens in two distinct phases, like a two-step cooking recipe.

• Stage 1: The Gluon Cascade (The "Branching Tree")
  Imagine a single drop of water hitting a surface and splashing. But instead of just splashing, the droplets keep splitting into smaller droplets, which split again, and again.
  In physics, the initial collision creates a quark and an antiquark. These immediately start shooting out "gluons" (the glue that holds particles together). These gluons split into more gluons, creating a massive, branching tree of energy. The paper uses math (Markov branching processes) to describe how this tree grows.

  • The Analogy: Think of a rumor starting in a room. One person tells two others, who tell two more, and soon the whole room is buzzing. This is the "cascade."

• Stage 2: Hadronization (The "Freezing")
  Eventually, the energy gets too low to keep splitting. The invisible "glue" particles (gluons and quarks) have to turn into real, visible particles (like pions or protons) that we can detect. This is called hadronization.

  • The Analogy: Imagine the hot, chaotic steam from a kettle suddenly hitting a cold window. It instantly condenses into distinct water droplets. The invisible energy "freezes" into solid, countable particles.

2. The Secret Ingredient: The "Second Correlation Moment"

How do the authors know their recipe is right? They look at a specific statistical number called the second correlation moment ($f_2$).

• The Analogy: Imagine you are counting the number of people in different families in a town.
  • At low energy (small town), families are small and predictable. The math shows a negative number here.
  • At high energy (huge city), families are chaotic, with huge variations in size. The math flips to a positive number.
  • The paper shows that as the collision energy increases, this number flips from negative to positive. This "flip" tells the authors exactly when the "splitting tree" (Stage 1) becomes more important than the "freezing" (Stage 2).

3. The "Gluon Dominance" Discovery

The most interesting finding is about who is doing the work.

• At Low Energies: The explosion is mostly driven by the original quarks. It's like a small fire started by a single match.
• At High Energies: The "gluons" take over. The paper finds that at high speeds, the number of gluons becomes massive (dozens of them), while the original quarks are just a tiny fraction of the total.
  • The Analogy: Imagine a small campfire (low energy) vs. a massive forest fire (high energy). In the forest fire, the original spark (the quark) is barely noticeable; the fire is now dominated by the thousands of burning trees (gluons) spreading everywhere. The model is called "Gluon Dominance" because these invisible glue particles are the main actors in the high-energy drama.

4. Changing the Rules of the Game

The paper also noticed something strange happening as energy gets really high (above 100 GeV).

• The "Fragmentation" vs. "Recombination" Switch:
  • At lower energies, one gluon turns into one particle (like a single seed growing into one plant).
  • At very high energies, the math suggests that one gluon might turn into more than one particle.
  • The Analogy: It's like switching from a factory where one machine makes one car, to a factory where the machines start merging and recombining to build cars faster and differently. This suggests the "glue" behaves differently when the crowd gets too dense.

5. Predicting the Future

Using this model, the author predicts what will happen in future, super-powerful particle accelerators (like those planned for 500 GeV or 1 TeV).

• The Prediction: They estimate that at these massive energies, a single collision could produce between 32 and 60 particles on average.
• Why it matters: Current computer simulations often fail to predict the "tail end" of the distribution (the rare, massive explosions). This model handles those rare, high-multiplicity events much better, giving physicists a clearer map for future experiments.

Summary

In short, this paper says: "We used to think the original particles (quarks) were the main drivers of these explosions. But we found that at high speeds, the invisible 'glue' (gluons) takes over, splitting wildly and creating a chaotic storm of particles. By understanding this two-step process (splitting then freezing), we can finally predict the results of these cosmic collisions with much higher accuracy."

This helps physicists build better experiments and understand the fundamental rules of how the universe builds matter from pure energy.

Technical Summary

1. Problem Statement

The paper addresses the challenge of accurately modeling multiparticle production (MP) in high-energy physics, specifically in electron-positron ($e^+e^-$) annihilation.

• The Limitation of Current Tools: While Monte Carlo event generators based on Quantum Chromodynamics (QCD) are standard for planning experiments, they often fail to accurately predict experimental data, particularly in the high-multiplicity tail of distributions. This is because perturbative QCD (pQCD) only describes the hard scattering (quark-gluon cascade) well, while the non-perturbative hadronization stage (conversion of partons to hadrons) relies on phenomenological models that often require ad-hoc adjustments.
• The Specific Gap: Although the Gluon Dominance Model (GDM), also known as the Two-Stage Model (TSM), has been successful in describing proton collisions and heavy ion interactions, a comprehensive, end-to-end re-analysis of $e^+e^-$ annihilation data across the entire available energy range (14 GeV to 189 GeV) using this unified framework was lacking. The authors aim to verify if the GDM can consistently describe $e^+e^-$ data without the discrepancies seen in other models, and to investigate the energy dependence of the underlying parameters to understand the transition in hadronization mechanisms.

2. Methodology

The authors employ the Gluon Dominance Model (GDM), which treats multiparticle production as a convolution of two distinct stages:

Stage 1: Quark-Gluon (qg) Cascade (Perturbative QCD)

• Mechanism: Modeled as a Markov branching process.
• Elementary Events:
  1. Quark bremsstrahlung: $q \to q + g$
  2. Gluon fission: $g \to g + g$
• Mathematical Formulation: The cascade is described using differential-difference equations for generating functions (GF).
  • Quark Jets: Follow a Negative Binomial Distribution (NBD) (Polya-Egenberger distribution).
  • Gluon Jets: Follow the Furry distribution.
• Key Parameter ($k_p$): The ratio of probabilities between quark bremsstrahlung and gluon fission.

Stage 2: Hadronization (Phenomenological)

• Mechanism: The transformation of partons into observable hadrons.
• Distribution: A Binomial (Bernoulli) distribution is used, chosen based on the behavior of the second correlation moment ($f_2$).
  • $f_2 = D^2 - \bar{n} = \overline{n(n-1)} - \bar{n}^2$.
  • The model accounts for the sign change of $f_2$ from negative (low energy) to positive (high energy).
• Parameters:
  • $\bar{n}_p$: Average number of hadrons formed from one parton $p$ (quark or gluon).
  • $N_p$: Maximum possible number of hadrons from one parton.
  • $\alpha$: A ratio parameter linking gluon and quark hadronization parameters ($\bar{n}_g = \alpha \bar{n}_q$).

The Convolution (The Model Equation)

The final multiplicity distribution $P_n(s)$ is obtained by convolving the cascade distribution with the hadronization binomial law:
$$P_n(s) = \Omega \sum_{m=0}^{M_g} P^q_m \binom{(2+\alpha m)N}{n} \left(\frac{\bar{n}_h}{N}\right)^n \left(1-\frac{\bar{n}_h}{N}\right)^{(2+\alpha m)N - n}$$
Where $P^q_m$ is the probability of having $m$ gluons in the quark jet, and $M_g$ is the limit of active gluons.

3. Key Contributions

1. Unified Description: The paper provides a successful, unified description of $e^+e^-$ annihilation multiplicity distributions from 14 GeV up to 189 GeV (including the $Z^0$ resonance region) using a single set of physical principles (GDM).
2. Resolution of High-Multiplicity Discrepancies: The model accurately predicts the "high-multiplicity tail" of the distribution, a region where standard Monte Carlo generators often fail or require significant tuning.
3. Physical Interpretation of Parameters: The authors extract and analyze the energy dependence of six key model parameters, linking them to physical mechanisms:
   • $k_p$ (Cascade ratio): Shows a sharp drop at low energies, indicating suppressed gluon fission below 22 GeV, followed by a stable value ($\approx 5.4$) at higher energies.
   • $\bar{n}_g$ (Gluon hadronization): Reveals a transition in the hadronization mechanism.
4. Low-Energy Validation: The model is successfully applied to 9.4 GeV data (PLUTO collaboration), where gluon fission is negligible, reducing the cascade to a Poisson distribution, further validating the model's robustness across energy scales.

4. Results

• Fit Quality: The GDM fits experimental data from collaborations (TASSO, AMY, OPAL) with $\chi^2$ values of only a few units. The agreement is particularly strong in the high-multiplicity region.
• Second Correlation Moment ($f_2$):
  • The model correctly reproduces the experimental observation that $f_2$ is negative at low energies (dominated by hadronization/binomial statistics) and becomes positive at high energies (dominated by the cascade/NBD statistics).
  • The zero-crossing of $f_2$ coincides with the energy threshold where gluon fission becomes significant.
• Parameter Behavior & Mechanism Transition:
  • Low Energy ($\sqrt{s} < 100$ GeV): The average number of hadrons per gluon ($\bar{n}_g$) is $\le 1$ (approx. 0.9). This supports the fragmentation mechanism (one gluon $\to$ one hadron) and the Local Parton-Hadron Duality (LPHD) hypothesis.
  • High Energy ($\sqrt{s} > 130$ GeV): $\bar{n}_g$ rises slightly above 1 (approx. 1.2). This suggests a shift toward a recombination mechanism, similar to that observed in dense quark-gluon media (proton collisions), implying a change in the hadronization dynamics at very high energies.
• Predictions for Future Colliders:
  • Using the fitted energy dependencies of $\bar{m}$ (gluon multiplicity) and $\bar{n}_q$ (quark hadronization), the authors predict the average charged multiplicity for future $e^+e^-$ accelerators:
    • 500 GeV: $\bar{n} \approx 32 - 39$
    • 1 TeV: $\bar{n} \approx 37 - 60$

5. Significance

• Theoretical Insight: The study demonstrates that a two-stage model combining perturbative QCD cascades with a phenomenological hadronization scheme based on correlation moments is sufficient to describe complex multiparticle production without needing complex, energy-dependent tuning of Monte Carlo generators.
• Mechanism Discovery: The analysis of the parameter $\bar{n}_g$ provides evidence for a change in the hadronization mechanism from pure fragmentation to a recombination-like process as energy increases, a finding that bridges the gap between $e^+e^-$ and hadronic collision physics.
• Future Physics: The results are directly relevant for the physics programs of upcoming facilities like the NICA collider (SPD facility) and the proposed eRHIC, where understanding gluon dominance and multiparticle production is crucial for studying the gluon structure of the nucleon and the formation of baryonic matter.
• Methodological Value: The paper highlights the sensitivity of multiparticle production studies to the choice of phenomenological models, advocating for the use of correlation moments ($f_2$) as a robust tool for constraining hadronization schemes.