Resonance- and Width-aware Parton Shower Evolution and NLO Matching

This paper presents a novel technique for next-to-leading order accurate simulation of $e^+e^-\to W^+W^-b\bar{b}$ that incorporates finite width effects beyond the Breit-Wigner approximation in resonance- and width-aware parton shower evolution and NLO matching, offering a publicly available tool for future electron-positron collider studies.

The Gist

Imagine you are trying to take a perfect photograph of a very fast, very fragile event in a particle collider. Specifically, you are watching two heavy particles (top quarks) collide, briefly exist, and then instantly explode into other particles (like a $W$ boson and a $b$-quark).

The problem is that these top quarks are like unstable soap bubbles. They exist for such a tiny fraction of a second that they don't have a single, fixed size or "mass." They are fuzzy, fluctuating, and have a "width" (a range of possible masses) rather than a sharp point.

For decades, physicists have used computer simulations (called "Parton Showers") to predict what happens when these bubbles pop and spray debris everywhere. However, the old simulations had a major flaw: they treated the bubbles like solid, unbreakable billiard balls.

The Problem: The "Recoil" Mistake

Think of the top quark bubble as a dancer spinning on a stage. When the dancer throws a heavy ball (a gluon) into the crowd, they naturally stumble backward (recoil).

In the old simulations, when the dancer threw a ball, the computer would make the entire stage (the other particles) stumble backward to balance the energy.

• The Issue: Because the top quark is a "soap bubble" (a resonance with a finite width), if you make the stage stumble, you accidentally change the size and shape of the bubble itself. You shift it off its natural "mass."
• The Consequence: This creates a ripple effect. The simulation thinks the bubble is heavier or lighter than it actually is, which throws off the entire calculation of how the explosion happens. It's like trying to measure the weight of a balloon while you're constantly squeezing it; your measurement will be wrong.

The Solution: "Resonance-Aware" and "Width-Aware" Dancing

The authors of this paper (Stefan Höche and Daniel Reichelt) invented a new way to simulate this dance. They call it "Resonance- and Width-Aware Parton Shower Evolution."

Here is how their new method works, using simple analogies:

1. Resonance-Aware: The "Family Group" Rule

In the old method, if a particle from the "Top Quark Family" threw a ball, the computer would make the entire universe stumble.
In the new method, the computer realizes: "Wait, this particle belongs to a specific family (the Top Quark decay). If this family member throws a ball, only the other members of that same family should stumble."

• Analogy: Imagine a group of friends holding hands in a circle. If one friend throws a ball, only the friends in that specific circle move to catch the balance. The people standing outside the circle (the other particles) stay perfectly still. This ensures the "bubble" (the top quark) doesn't get squished or stretched by the movement of unrelated particles.

2. Width-Aware: The "Fuzzy Edge" Rule

Even with the "Family Group" rule, there was still a problem near the "threshold" (the exact energy where the top quarks are created).

• The Old Way: The simulation assumed the top quark was a sharp, hard sphere.
• The New Way: The simulation acknowledges the top quark is a fuzzy cloud. It knows that near the edge of the cloud, the rules of physics change slightly. The authors added a special "fuzzy factor" to their math that accounts for the fact that the top quark isn't a point, but a cloud with a specific width.

Why Does This Matter? (The Future of Physics)

The authors tested this new method for a future particle collider called the FCC-ee (Future Circular Collider). This machine is designed to be the "Hubble Telescope" of particle physics, capable of measuring things with incredible precision.

• The Goal: They want to measure the mass of the top quark to within 50 millionths of a gram (50 MeV). That is like measuring the weight of a single grain of sand on a mountain.
• The Result: If you use the old "billiard ball" simulations, your measurement will be off by a huge margin because the "squishing" of the bubble introduces errors. The new "fuzzy, family-aware" simulation removes these errors.

The Takeaway

This paper is essentially a manual for a better camera lens.

For years, physicists were taking photos of the universe with a slightly blurry lens that distorted the shape of the most important objects (the top quarks). This paper introduces a new lens that:

1. Keeps the object steady (Resonance-aware) so it doesn't get distorted by the background.
2. Respects the object's natural fuzziness (Width-aware) so it doesn't try to force a cloud into a hard shape.

With this new tool, scientists can finally take a crystal-clear picture of the top quark, allowing them to test the fundamental laws of the universe with unprecedented accuracy. They have also released the software (based on the ALARIC and SHERPA tools) so other scientists can use this new "lens" immediately.

Technical Summary

1. Problem Statement

The paper addresses a critical limitation in simulating the process $e^+e^- \to W^+W^-b\bar{b}$ (top-quark pair production and decay) near the production threshold ($\sqrt{s} \approx 2m_t$). This process is vital for future high-precision lepton colliders like the FCC-ee, where top-quark mass measurements require precision at the level of 50 MeV.

Key Challenges Identified:

• Finite Width Effects: The top quark decays before hadronization, meaning its finite width ($\Gamma_t$) must be treated dynamically. Traditional parton showers often assume resonances are on-shell or treat them only via Breit-Wigner distributions, neglecting interference and width effects near the threshold.
• Recoil and Virtuality: In standard parton showers, when a parton (e.g., a $b$-quark) radiates a gluon, the recoil is distributed to a color partner (e.g., the $\bar{b}$). This momentum reshuffling alters the virtuality of the intermediate top-quark propagators ($q^2 \neq m_t^2$).
• Inconsistency: Changing the virtuality of the intermediate resonance significantly alters the Born-level cross-section rate (by factors of $O(1)$) due to the steepness of the Breit-Wigner peak. This violates the assumption that parton-shower evolution is a Markovian process with sub-leading power corrections.
• NLO Matching: Existing "resonance-aware" schemes (which preserve resonance virtuality) fail to account for finite width effects in the radiation kernels themselves, specifically the interference between the two top decays near the threshold.

2. Methodology

The authors propose a novel framework combining the ALARIC parton shower with the SHERPA event generator, utilizing S-MC@NLO matching. The solution involves three main technical pillars:

A. Resonance-Aware Evolution (Preserving Virtuality)

To prevent radiation from altering the top-quark virtuality, the authors modify the recoil scheme:

• Recoil Assignment: If a parton originating from a specific top decay splits, the recoil is assigned only to the other decay products of that same top quark, rather than the global color partner.
• Kernel Modification: The standard scalar eikonal kernel (Eq. 5 in the paper), which depends on a global recoil vector, is replaced by the original Catani-Seymour scalar operator (Eq. 7). This operator uses partial fractioning based on scalar invariants, making it independent of the specific recoil vector choice, thus allowing the resonance-aware recoil scheme to function without breaking the dipole structure.

B. Width-Aware Evolution (Threshold Effects)

To account for finite width effects near the $t\bar{t}$ threshold, the authors derive a modified splitting kernel:

• Interference Treatment: They analyze the matrix element squared near threshold, identifying contributions from individual top-bottom dipoles and the interference term between the two top decays.
• Modified Kernel: They introduce a width-dependent factor, $\chi(z)$, into the scalar splitting kernel (Eq. 14). This factor interpolates between the standard dipole radiation pattern and the interference-suppressed pattern based on the gluon energy relative to the top width.
  $$V^{(s,w)} \propto \dots + \chi(z) \left( \dots \right)$$
  where $\chi(z)$ regulates the soft divergence and suppresses radiation that would otherwise violate the resonance structure near the threshold.

C. Fixed-Order Infrared Subtraction

To ensure NLO accuracy, the authors derive the corresponding integrated infrared subtraction counterterms:

• Analytic Derivation: They compute the integrated versions of the new resonance- and width-aware kernels.
• New Expressions: They provide analytic expressions for the finite remainder terms involving the width-dependent factor (Eq. 20), which includes transcendental functions (dilogarithms) and arctangents dependent on the parameter $\gamma \propto \Gamma_t$.
• Validation: These subtraction terms are implemented in SHERPA to cancel infrared singularities in the real-emission corrections.

3. Key Contributions

1. New Algorithm: A fully consistent NLO-matched parton shower algorithm that simultaneously preserves the virtuality of intermediate resonances and incorporates finite width effects near the production threshold.
2. Analytic Results: Derivation of new analytic expressions for integrated infrared subtraction terms specific to width-aware evolution, extending the Catani-Seymour formalism.
3. Public Implementation: The technique is implemented in the ALARIC parton shower within the SHERPA event generator, making it available for public use.
4. Phenomenological Application: First application of this technique to $e^+e^- \to W^+W^-b\bar{b}$ at $\sqrt{s} = 365$ GeV (FCC-ee energy).

4. Results

The paper presents validation and phenomenological results comparing the Default, Resonance-Aware (Res. Aware), and Resonance + Width-Aware (Res. + Width) schemes:

• Cross-Section Validation: The integrated cross-sections for $e^+e^- \to W^+W^-b\bar{b}$ show excellent agreement (sub-permille level) between the different subtraction schemes, confirming the correctness of the new analytic counterterms.
• Radiation Patterns (Lund Plane):
  • The Width-Aware scheme significantly alters the rapidity distribution of the first emitted gluon compared to the standard scheme.
  • It suppresses emissions from the $b$-quark into the $\bar{b}$ hemisphere (and vice versa) in the central rapidity region, reflecting the interference effects near the threshold.
  • The standard scheme overestimates radiation in these regions, while the width-aware scheme correctly recovers the suppression predicted by fixed-order calculations.
• Top Mass Reconstruction:
  • Default/Res. Aware: Show significant distortions in the reconstructed top-mass tail due to clustering mismatches where radiation is incorrectly assigned to the wrong jet.
  • Width-Aware: Closely follows the fixed-order NLO prediction, preserving the shape of the Breit-Wigner distribution.
• Jet Multiplicity: The width-aware scheme predicts fewer additional jets compared to default schemes, consistent with the suppression of soft wide-angle radiation near the threshold.

5. Significance

• Precision Physics: This work is essential for the FCC-ee physics program. The differences between resonance-aware and width-aware evolution are large enough to impact template fits for the top-quark mass, potentially introducing biases if not corrected.
• Theoretical Advancement: It resolves the apparent contradiction between maintaining resonance virtuality (required for NLO matching) and modeling finite width effects (required for threshold physics).
• Generalizability: While focused on $t\bar{t}$ at lepton colliders, the methodology for resonance-aware recoil and width-dependent kernels is general and can be applied to other processes involving internal resonances, such as top production at the LHC or Higgs decays, provided the kinematic constraints are met.

In summary, the paper provides a rigorous, NLO-accurate solution for simulating top-quark pair production near threshold, correcting systematic errors in current Monte Carlo generators that are critical for future high-precision collider experiments.