Improved Standard-Model predictions for $\eta^{(\prime)}\to \ell^+ \ell^-$

This paper presents improved Standard Model predictions for the rare dilepton decays $\eta^{(\prime)}\to \ell^+\ell^-$ by leveraging recent advances in transition form factor calculations and a robust dispersive evaluation of subleading contributions, yielding precise branching fractions that reveal a mild tension with experimental data for $\eta\to\mu^+\mu^-$ and provide new constraints on physics beyond the Standard Model.

The Gist

Imagine the universe as a giant, intricate machine where tiny particles constantly bump into each other, transform, and vanish. In this paper, a team of physicists acts like master mechanics trying to predict exactly how often a specific, rare event happens: a heavy, unstable particle (called an eta meson) spontaneously splitting into a pair of lighter particles (an electron or a muon and their anti-particles).

Here is the story of their work, broken down into simple concepts and analogies.

1. The Rare Event: The "Ghostly" Split

In the Standard Model (our current best rulebook for physics), this split is incredibly rare. It's like trying to win a lottery where the odds are one in a billion.

Why is it so rare?

• The Helicity Trap: Imagine the eta meson is a spinning top. When it tries to split into two lighter particles, the laws of physics say those new particles must spin in a specific way. But because they are so light (especially the electron), they are "lazy" and don't want to spin that way. This makes the event very difficult to happen.
• The Detour: The particle can't just split directly. It has to take a detour. It briefly turns into two "ghostly" flashes of light (virtual photons) before turning into the final particles. This extra step makes the event even rarer.

2. The Map: The "Transition Form Factor"

To predict how often this happens, the physicists needed a map. In their world, this map is called the Transition Form Factor (TFF).

Think of the TFF as a detailed blueprint of the eta meson's internal structure. It tells us exactly how the meson interacts with those "ghostly" flashes of light.

• Previous Maps: Before this paper, scientists had rough sketches of this blueprint. They were okay for the lightest particle (the pion), but for the heavier eta mesons, the sketches were blurry.
• The New Map: The authors used a technique called dispersion theory. Imagine you are trying to guess the shape of a hidden object by listening to the echoes of sound bouncing off it. They used data from many different experiments (the echoes) to reconstruct the blueprint with incredible precision. They didn't just guess; they built the map piece by piece using real-world data.

3. The Twist: The "Imaginary" Parts

In physics, calculations often have a "real" part and an "imaginary" part. Don't let the name scare you; "imaginary" here just means a specific mathematical component related to probability and timing.

• The Old View: Scientists used to think the only thing that mattered was the "two-photon" path (the main detour).
• The New Discovery: The authors found that for the heavier eta mesons, there are other hidden paths. Imagine the meson doesn't just turn into two flashes of light; sometimes it briefly turns into a pion and a photon, or other combinations, before settling down.
• The Analogy: If you are walking from Point A to Point B, you usually take the highway (two photons). But for the heavier mesons, there are also scenic backroads (other particle combinations) that you can take. The authors calculated exactly how much these backroads affect the final arrival time. They found these backroads are much more important for the heavy eta meson than for the light pion.

4. The Prediction: The "Crystal Ball"

Using their new, ultra-precise map and accounting for those hidden backroads, the team made a prediction for how often these splits happen.

They calculated the "Branching Fraction," which is just a fancy way of saying: "Out of every 10 billion eta mesons that decay, how many will turn into electrons or muons?"

• The Result: They gave very specific numbers with tiny error bars (like saying "5.37 out of 10 billion, give or take 0.04").
• The Tension: When they compared their crystal ball prediction to what experiments have actually measured so far, they found a slight mismatch.
  • For the muon version, the experiment sees slightly more events than the theory predicts. It's like the lottery is happening slightly more often than the math says it should.
  • The mismatch is about 1.6 standard deviations. In the world of physics, this isn't a "smoking gun" (which usually requires 5 sigma), but it's a "raised eyebrow." It's a hint that maybe there is something new we don't understand yet.

5. The Search for New Physics

Why does this matter?
If the experiment keeps seeing more events than the Standard Model predicts, it could mean New Physics.

• The Analogy: Imagine you are a detective. You know exactly how a safe should open based on the blueprints. But every time you try to open it, it clicks open a little faster than expected.
  • Scenario A: Maybe your blueprint was slightly wrong (which is what this paper is trying to fix).
  • Scenario B: Maybe there is a hidden key (a new particle or force) helping the safe open.

The authors used their new, precise prediction to set limits on these "hidden keys." They calculated: "If there is a new particle causing this, it must be heavier than X, or it must interact in a very specific way."

Summary

This paper is a masterclass in precision.

1. The authors built the most accurate map (TFF) ever for these specific particles.
2. They accounted for hidden detours (subleading channels) that previous models ignored.
3. They provided a new, ultra-precise prediction for a rare event.
4. They found a small, intriguing hint that the universe might be doing something slightly different than our current rulebook says, specifically for the muon version of the decay.

It's a reminder that even in the well-trodden territory of the Standard Model, there are still tiny cracks in the pavement where new discoveries might be hiding, waiting for a more precise map to reveal them.

Technical Summary

1. Problem Statement

The rare decays of neutral pseudoscalar mesons ($P = \pi^0, \eta, \eta'$) into lepton pairs ($P \to \ell^+\ell^-$) are highly suppressed in the Standard Model (SM) due to:

• Chirality suppression: The decay rate scales with $m_\ell^2/M_P^2$.
• Loop suppression: The process proceeds via a two-photon exchange loop ($P \to \gamma^*\gamma^* \to \ell^+\ell^-$).

These decays are sensitive probes for Physics Beyond the Standard Model (BSM), as BSM scenarios (e.g., new light particles or effective operators) can lift the helicity suppression. While the SM prediction for $\pi^0 \to e^+e^-$ has reached high precision, predictions for the heavier $\eta$ and $\eta'$ mesons have historically suffered from larger uncertainties. The primary challenge lies in the precise determination of the Transition Form Factor (TFF), $F_{P\gamma^*\gamma^*}(q_1^2, q_2^2)$, which describes the $P \to \gamma^*\gamma^*$ vertex. Previous approaches often relied on model-dependent approximations (like Vector Meson Dominance) or dispersion relations restricted to the two-photon cut, which become insufficient for the heavier $\eta$ and $\eta'$ where other intermediate states and mass effects are significant.

2. Methodology

The authors employ a dispersive approach to reconstruct the TFFs, leveraging recent high-precision studies of $\eta^{(\prime)}$ poles in the context of Hadronic Light-by-Light (HLbL) scattering for the muon $g-2$.

• TFF Decomposition: The normalized TFF is decomposed into isovector ($I=1$), isoscalar ($I=0$), effective pole, and asymptotic contributions:
  $$\tilde{F}_{P\gamma^*\gamma^*} = \tilde{F}^{(I=1)} + \tilde{F}^{(I=0)} + \tilde{F}^{\text{eff}} + \tilde{F}^{\text{asym}}$$
• Isovector Contribution: Reconstructed via a double-spectral representation using $\pi\pi$ intermediate states. This involves solving an inhomogeneous Muskhelishvili–Omnès problem, incorporating $\pi\pi$ P-wave phase shifts and $a_2(1320)$ tensor meson exchange (left-hand cuts).
• Isoscalar Contribution: Described using a narrow-resonance Vector Meson Dominance (VMD) ansatz for $\omega$ and $\phi$ mesons, with finite-width effects checked and found negligible.
• Effective Poles: Introduced to satisfy normalization sum rules and match high-energy space-like data ($e^+e^- \to e^+e^- \eta^{(\prime)}$).
• Asymptotic Contributions & Mass Corrections: A key methodological advancement is the inclusion of pseudoscalar mass corrections in the asymptotic region. The authors generalize the standard light-cone expansion to include terms proportional to $M_P^2$, which are crucial for the $\eta'$ where $M_{\eta'}$ is comparable to the matching scale.
• Loop Integration: The reduced amplitude $A_\ell$ is calculated by integrating the TFF over the loop momentum. The authors derive new analytic expressions for the loop integrals (using Passarino–Veltman reduction) that consistently include the massive pseudoscalar corrections in the asymptotic kernel.
• Imaginary Parts: Unlike the $\pi^0$, where the two-photon cut dominates the imaginary part, the $\eta$ and $\eta'$ receive significant contributions from subleading channels (e.g., $2\pi\gamma$, $3\pi\gamma$). The dispersive framework allows for a robust evaluation of these imaginary parts beyond the unitarity bound.

3. Key Contributions

1. First Robust Evaluation of Subleading Imaginary Parts: For the first time, the dispersive representation allows for a rigorous calculation of the imaginary parts of the amplitude arising from channels other than the two-photon cut. This is particularly important for $\eta'$, where these corrections are large.
2. Inclusion of Pseudoscalar Mass Corrections: The paper derives and implements a generalized asymptotic form of the TFF that accounts for $M_P^2$ effects. This corrects previous approximations that assumed massless pseudoscalars, which were inadequate for the $\eta'$.
3. Unified Framework: The work unifies the TFF analysis for $\eta^{(\prime)}$ dilepton decays with the high-precision HLbL scattering calculations used for the muon $g-2$, ensuring consistency across different low-energy QCD observables.
4. BSM Constraints: The paper translates the improved SM predictions into constraints on effective operators (axial-vector, pseudoscalar, gluonic) and light new particles ($Z'$, axion-like particles).

4. Results

The authors provide Standard Model predictions for the branching ratios with uncertainties at the few-percent level. The results are:

Decay Mode	Branching Ratio Prediction	Comparison to Experiment
$\eta \to e^+e^-$	$5.37(4)(2)[4] \times 10^{-9}$	Limit: $< 7.0 \times 10^{-7}$
$\eta \to \mu^+\mu^-$	$4.54(4)(2)[4] \times 10^{-6}$	Exp: $5.8(8) \times 10^{-6}$
$\eta' \to e^+e^-$	$1.80(2)(3)[3] \times 10^{-10}$	Limit: $< 5.6 \times 10^{-9}$
$\eta' \to \mu^+\mu^-$	$1.22(2)(2)[3] \times 10^{-7}$	No limit available

• Uncertainty Breakdown: The errors are categorized into dispersive systematics, uncertainties from space-like data (Brodsky-Lepage), asymptotic contributions, and isoscalar weights.
• Tension with Experiment: The prediction for $\eta \to \mu^+\mu^-$ exhibits a $1.6\sigma$ tension with the current experimental average. The authors note that a new measurement (e.g., by REDTOP) is needed to clarify this.
• Imaginary Part Corrections: The inclusion of subleading cuts reduces the magnitude of the imaginary part for $\eta'$ significantly (e.g., from $-23.68$ to $-18.87$ for the electron channel), leading to lower total branching fractions compared to naive unitarity bounds.

5. Significance

• Precision Benchmark: The paper establishes the most precise SM predictions to date for $\eta^{(\prime)} \to \ell^+\ell^-$, reducing theoretical uncertainties to a level where experimental improvements of up to three orders of magnitude (for electron modes) could be fully utilized before theory becomes the limiting factor.
• BSM Sensitivity: The improved precision tightens constraints on BSM physics. Specifically, pseudoscalar and gluonic operators benefit from chiral enhancement ($1/m_\ell$), allowing the $\eta^{(\prime)} \to e^+e^-$ channels to probe energy scales up to $\sim 1.6$ TeV, despite the current experimental limits being far above the SM prediction.
• Theoretical Consistency: By resolving the tension between different intermediate state contributions and mass effects, the work provides a consistent picture of pseudoscalar TFFs, bridging the gap between low-energy phenomenology and high-energy constraints.
• Future Directions: The authors highlight the need for new experimental data on $\eta \to \mu^+\mu^-$ and independent validation of the $\eta \to \gamma\gamma$ normalization to resolve the current $1.6\sigma$ discrepancy. They also suggest that future lattice QCD calculations of the doubly-virtual TFF would provide a valuable cross-check.

In summary, this work represents a major step forward in the precision phenomenology of rare meson decays, utilizing advanced dispersive techniques to control theoretical uncertainties and providing a solid foundation for testing the Standard Model and searching for new physics.