⚛️ phenomenology

Gauge-Invariant Longitudinal Modes in the Herwig 7 Electroweak Parton Shower

This paper presents a gauge-invariant scheme for longitudinal electroweak gauge bosons in the Herwig 7 parton shower by incorporating Ward-identity-fixed Goldstone-matching terms, which ensures numerical stability and controlled corrections at lower energy scales while preserving the transverse sector and yielding interpretable shifts in electroweak-sensitive observables.

Original authors: M. R. Masouminia, P. Richardson

Published 2026-03-19

📖 6 min read🧠 Deep dive

CC0 1.0

Original authors: M. R. Masouminia, P. Richardson

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict the chaotic aftermath of a high-speed car crash. In the world of particle physics, this "crash" happens inside the Large Hadron Collider (LHC), where protons smash together at near light speed. To understand what happens next, physicists use a computer program called a Parton Shower. Think of this program as a simulation that tracks how the debris (particles) sprays out, bounces off each other, and splits into smaller pieces, much like a tree branching out or a river delta splitting into smaller streams.

For decades, this simulation has been very good at tracking the "side-to-side" wiggles of particles (called transverse modes). However, there is a tricky, delicate part of the simulation involving particles that move "up and down" or "forward and backward" (called longitudinal modes).

The Problem: The "Ghost" in the Machine

In the Standard Model of physics, some particles (like the W and Z bosons) have mass. When these heavy particles move, they have a special "longitudinal" way of vibrating. Mathematically, describing this vibration is messy. It's like trying to describe the motion of a heavy pendulum using a formula that includes a "ghost" term—a mathematical artifact that shouldn't be there but appears because of how we chose to write the equations.

In the current version of the software (Herwig 7), the programmers had a clever workaround. They simply subtracted this "ghost" term from the equations to make the numbers stable. It worked well in many situations, like when the particles are moving very fast (high energy).

But, there was a catch. In physics, you can't just delete a piece of a puzzle and expect the picture to remain perfect. The "ghost" term is actually connected to a hidden piece of the puzzle called a Goldstone boson (a theoretical particle that gives mass to others). By simply subtracting the ghost, the old method was ignoring this hidden connection. It was like fixing a leaky roof by taping over the hole but forgetting to fix the broken beam underneath. At lower energies or in specific scenarios, this "quick fix" could lead to slightly wrong predictions.

The Solution: The "Gauge-Invariant" Patch

The authors of this paper, Masouminia and Richardson, decided to fix the roof properly. They developed a new, Gauge-Invariant (GI) scheme.

Instead of just deleting the "ghost" term, they:

Kept the remainder of the subtraction (the part that was left over).
Added a new piece that represents the "Goldstone boson" connection, which the laws of physics (specifically something called "Ward identities") demand must be there.

Think of it like this:

Old Method (SL): You have a recipe for a cake. The recipe says "add 1 cup of flour, then subtract the ghost of the flour." You get a cake, but it's missing a secret ingredient that makes it taste right.
New Method (GI): You realize the "ghost of flour" was actually a placeholder for "baking powder." So, you keep the flour, and you explicitly add the baking powder. Now, the cake is mathematically perfect and tastes exactly as nature intended, no matter how you slice it.

What Did They Find?

The team implemented this new "perfect recipe" into the Herwig 7 software and ran thousands of simulations to see what changed.

High Energy is Fine: When the particles are moving incredibly fast (high energy), the old method and the new method give almost the exact same results. The "ghost" was so small it didn't matter.
Low Energy Matters: When the particles are moving slower or are heavy (like the top quark or the W boson), the new method makes a difference. It corrects the predictions for how often these particles split and in what direction they fly.
The "Sudakov" Surprise: One of the most interesting findings is about probability. In these simulations, there's a "no-emission" rule (Sudakov suppression). If the new method says a particle is more likely to split, the computer also says it's less likely to survive without splitting.
- Analogy: Imagine a crowded dance floor. If the new rules say people are more likely to bump into each other (split), the simulation also says fewer people will make it through the night without bumping into anyone. This creates a weird, non-linear effect where a small change in the rules can actually make the final count of "bumpers" go up or down in unexpected ways.

Why Does This Matter?

This isn't just about fixing a math error; it's about precision.

For Discovery: As we look for new physics (like new particles) at the LHC, we need to know the "background noise" (the Standard Model) perfectly. If our simulation of the background is slightly off because of a "ghost" term, we might mistake a simulation error for a new discovery, or miss a real discovery because we thought it was just a glitch.
For Stability: The authors proved that their new method is stable. It doesn't crash the computer, and it doesn't break the other parts of the simulation. It just tweaks the longitudinal parts to be mathematically honest.

The Bottom Line

The authors have upgraded the "engine" of the particle physics simulation. They replaced a "quick fix" for heavy particles with a "proper fix" that respects the deep symmetries of the universe.

For the general public: It's like upgrading a GPS. The old GPS got you to the destination 99% of the time. The new GPS fixes the 1% of cases where the map was slightly wrong, ensuring that when we are looking for the most elusive treasures in the universe, our map is as accurate as possible.
For the scientists: It provides a "Gauge-Invariant" baseline. Now, when they see a difference in their data, they can be more confident that it's real physics, not just an artifact of how they chose to do the math.

The paper concludes that while the changes might seem small, they are crucial for the next generation of precision experiments, ensuring that our understanding of the universe's building blocks is as solid as the laws of physics themselves.

1. Problem Statement

Electroweak (EW) parton showers are essential for precision phenomenology at high-energy colliders (multi-TeV regime), where large logarithms ( $\alpha_W \ln^2(Q^2/M_W^2)$ ) compete with QCD effects. However, implementing these showers in a broken gauge theory (Standard Model) presents a specific technical challenge regarding longitudinal massive vector bosons ( $W^\pm, Z^0$ ).

The Gauge Artifact: The longitudinal polarization vector $\epsilon^\mu_L(p)$ contains a term proportional to $p^\mu/M_V$ . In a broken theory, this term does not vanish diagram-by-diagram; instead, Ward identities relate it to amplitudes involving the associated would-be Goldstone boson ( $G_V$ ) and Yukawa interactions for massive fermions.
The Limitation of Current Methods: The default Herwig 7 implementation uses a Subtraction-Longitudinal (SL) scheme (Dawson's approach). This scheme subtracts the $p^\mu/M_V$ component to ensure finite, numerically stable splitting kernels. While valid in the Equivalent Vector Boson Approximation (EWA) limit, the SL scheme effectively discards the Goldstone contribution.
The Risk: Outside the strict EWA limit, the SL scheme may violate gauge invariance at the level of the splitting operator. It risks mis-modeling longitudinal dynamics in the ultracollinear region (where transverse momenta are $\sim M_V$ ) and failing to capture symmetry-breaking effects (Yukawa couplings) that are crucial for massive fermions (e.g., top/bottom quarks).

2. Methodology

The authors construct and implement a Gauge-Invariant (GI) scheme to replace or supplement the default SL approach in Herwig 7.

Theoretical Framework:
- The GI scheme promotes the subtraction remainder (the SL current) to one component of the physical longitudinal current.
- It adds a Ward-identity-fixed Goldstone-matching term to complete the current. This ensures the longitudinal splitting operator satisfies the broken-theory Ward identities: $p_\mu \mathcal{M}^\mu \propto M_V \mathcal{M}_{Goldstone}$ .
- The construction is applied to two primary branching types:
  1. Fermion emission: $q \to q' V$ (where $V=W, Z$ ).
  2. Vector boson emission: $V \to V' V''$ .
Derivation of Kernels:
- The authors derive helicity-resolved splitting amplitudes ( $\hat{H}$ ) and compact quasi-collinear splitting kernels ( $F$ ) for both SL and GI prescriptions.
- They utilize a spin-density matrix formalism to retain spin correlations within the Markovian branching framework.
- The GI completion modifies only the longitudinal ( $0^*$ ) entries of the splitting functions. The transverse sector remains identical to the SL scheme.
- For $q \to q' V$ , the GI kernel includes interference terms ( $\propto c_G$ ) and pure Goldstone terms ( $\propto c_G^2$ ) involving Yukawa couplings ( $y_L, y_R$ ).
- For $V \to V' V''$ , the completion involves momentum contractions of the triple-gauge vertex with the Goldstone boson.
Implementation:
- The scheme is implemented in Herwig 7 as a switchable option (GaugeInvariant vs. Subtraction).
- It preserves backward compatibility and allows for direct comparison within identical shower setups.
- The implementation includes a parameter $c_G$ (default 1.0) to control the Goldstone-matching coefficient, allowing for stress tests.

3. Key Contributions

Formal Construction: A rigorous derivation of gauge-invariant longitudinal splitting operators that satisfy Ward identities in the broken Standard Model, explicitly linking the vector current to the Goldstone sector.
Analytic Kernels: The provision of closed-form expressions for helicity-dependent splitting functions for both fermion and boson branchings, separating transverse and longitudinal contributions.
Software Integration: A fully functional implementation in Herwig 7, including user-interface switches for the longitudinal scheme and Goldstone couplings.
Validation Framework: A comprehensive validation strategy comparing the SL and GI schemes at three levels:
1. Analytic splitting functions.
2. Single-resummed evolution (isolating one emission).
3. Full-resummed evolution (multi-emission LHC-like environment).

4. Results

The study compares the default SL scheme against the new GI scheme across various observables:

Analytic Level:
- The GI and SL schemes coincide in the transverse sector and in the strict massless/collinear limit.
- Differences arise in the longitudinal sector, driven by symmetry-breaking mass terms ( $m^2/Q^2$ ) and Yukawa couplings.
- Example ( $b \to t W$ ): The GI kernel shows a significant deformation (O(1) shift) compared to SL due to the large top-quark Yukawa coupling enhancing the Goldstone contribution.
- Example ( $b \to b Z$ ): The difference is negligible as Goldstone couplings are suppressed by $m_b/M_W$ .
Single-Resummed Evolution:
- In exclusive single-emission studies, the GI prescription generally yields a larger splitting rate for bosonic channels ( $V \to V' V''$ ) compared to SL.
- For fermionic channels ( $q \to q' V$ ), the observed spectrum can sometimes be smaller in GI than SL despite a larger kernel. This is explained by Sudakov suppression: a locally larger emission probability increases the no-emission probability (Sudakov form factor), reducing the exclusive one-emission cross-section. This highlights the non-monotonic nature of scheme dependence in exclusive observables.
Full-Resummed Evolution (LHC-like):
- The implementation is numerically stable.
- Angular Observables: Distributions like $\Delta R$ and $\cos \theta$ remain unchanged between SL and GI, confirming that the radiation pattern coherence and recoil maps are not destabilized.
- Kinematic Shifts: Differences appear in rate/tilt of $p_\perp$ and the reconstructed lightcone fraction $z$ .
- Hadronization Effects: The study notes that reconstructed variables referencing the event-record "parent" (like $z_W$ ) are sensitive to hadronization (non-perturbative momentum reshuffling), whereas hard, perturbative spectra (like $Z^0$ observables) remain stable.

5. Significance

Theoretical Consistency: The GI scheme resolves the ambiguity of the SL approach by ensuring that longitudinal currents are strictly gauge-invariant, correctly accounting for the interplay between vector bosons and Goldstone modes in the broken phase.
Phenomenological Impact: While inclusive observables dominated by on-shell matrix elements are largely unaffected, splitting-tagged observables and those probing the ultracollinear regime (low $Q^2$ , high symmetry breaking) show controlled, interpretable shifts. This is crucial for precision measurements involving heavy fermions (top quarks) and high-multiplicity EW radiation.
Robustness: The paper demonstrates that the GI completion can be integrated into a general-purpose event generator without disrupting the shower's coherence or stability.
Future Utility: The framework provides a baseline for studying Beyond Standard Model (BSM) scenarios with modified chiral couplings or triple-gauge couplings, where the treatment of longitudinal modes is even more critical.

In conclusion, this work establishes a theoretically sound and practically robust method for handling longitudinal electroweak gauge bosons in parton showers, moving beyond the approximations of the subtraction-only method to a fully gauge-invariant formulation.