Probing Lepton-Flavor-Violating Four-Lepton Operators at a Muon Collider

Technical Summary: Probing Lepton-flavour-Violating Four-Lepton Operators at a Muon Collider

Problem Statement
Charged lepton flavour violation (LFV) in the charged lepton sector remains unobserved, constituting a sensitive probe for physics beyond the Standard Model (SM). While neutrino oscillations confirm LFV in the neutral sector, charged LFV processes in the minimally extended SM (with neutrino masses) are suppressed by loop factors and $(\Delta m_\nu^2/M_W^2)^2$, rendering them unobservable. Consequently, any observation of charged LFV would signal new physics. While various Beyond-Standard-Model (BSM) scenarios (e.g., heavy neutrinos, multi-Higgs models, supersymmetry) predict observable LFV, a model-independent description is provided by the Standard Model Effective Field Theory (SMEFT).

This work addresses the sensitivity of a future high-energy muon collider to dimension-six four-lepton operators within the SMEFT framework. Specifically, it investigates the processes $\mu^+\mu^- \to e^\pm\tau^\mp$, $\mu^+\mu^- \to e^\pm\mu^\mp$, and $\mu^+\mu^- \to \mu^\pm\tau^\mp$. These processes are forbidden in the SM and are mediated at tree level by four-fermion contact interactions. The study aims to determine the projected sensitivity to the Wilson coefficients of these operators, accounting for beam polarization, hadronic $\tau$ reconstruction, and the evolution of constraints from low-energy scales to multi-TeV collider energies.

Methodology
The analysis employs a "bottom-up" strategy, starting with experimentally allowed Wilson coefficients defined near the electroweak scale and evolving them upward to the characteristic energy of a multi-TeV muon collider using Renormalisation Group Equations (RGE).

1. Effective Lagrangian and Operator Basis:
   The study focuses on the Warsaw basis of dimension-six four-lepton operators ($O_{\ell\ell}, O_{ee}, O_{\ell e}$). Due to the fixed initial state ($\mu^+\mu^-$), only operators with flavor indices $(p,r)=(\mu,\mu)$ contribute at tree level. The operators are mapped to a low-energy vector and scalar contact-interaction basis using Fierz rearrangements. The analysis assumes real Wilson coefficients to avoid CP-violating phases.

2. Constraints and RGE Evolution:

   • Low-Energy Bounds: Constraints are derived from current experimental upper limits on LFV decays (e.g., $\tau \to \mu\mu\mu$, $\tau \to \mu\mu e$) and $\mu \to e\gamma$. For processes involving $\tau$ in the final state, bounds are extracted directly from $\tau$ decay limits. For the $\mu^\pm e^\mp$ channel, where direct tree-level constraints on purely leptonic operators are weak or non-existent, the analysis relies on the lack of strong low-energy bounds to justify the collider probe.
   • RGE: The Wilson coefficients are evolved from the low-energy scale ($\mu_{RS} = m_\tau$) to the collider scale ($\mu_{RS} = \sqrt{s}$) using one-loop RGEs. The evolution is dominated by the gauge sector, with Yukawa contributions being negligible at one-loop order. This ensures the parameter space remains within the perturbative reliability of the EFT.

3. Event Generation and Simulation:

   • Processes: Signal events for $\mu^+\mu^- \to e^\pm\tau^\mp$, $\mu^+\mu^- \to \mu^\pm\tau^\mp$, and $\mu^+\mu^- \to e^\pm\mu^\mp$ are simulated at center-of-mass energies $\sqrt{s} = 3, 10,$ and $14$ TeV.
   • Backgrounds: Dominant SM backgrounds include $\mu^+\mu^- \to W^+W^-$, $\mu^+\mu^- \to \tau^+\tau^-$, $\mu^+\mu^- \to Z\nu\bar{\nu}$, and $\mu^+\mu^- \to \mu^+\mu^-$ (for the $e\mu$ channel). Beam-induced backgrounds (BIB) are considered but found to be suppressed by hard kinematic cuts.
   • Detector Simulation: Events are generated using MadGraph5_aMC@NLO, showered with Pythia8, and processed through Delphes3 with a muon collider-specific detector card. Hadronic $\tau$ decays are reconstructed using the Valencia algorithm in FastJet. The analysis focuses exclusively on hadronically decaying $\tau$ leptons to maximize signal yield and utilize the larger branching fraction (~65%).

4. Analysis Strategy:

   • Kinematic Selection: A hard transverse momentum ($p_T$) cut is applied to the prompt lepton ($e$ or $\mu$) to exploit the distinct kinematic features of the contact interaction signal (peaked near $p_T \sim \sqrt{s}/2$) versus the falling SM background spectrum.
   • Optimal Observable Technique (OOT): To extract Wilson coefficients and resolve correlations, the study employs an optimal-observable analysis of the angular distribution ($\cos\theta$) of the prompt charged lepton. This method utilizes the differential cross-section to construct a covariance matrix, allowing for a global fit of the Wilson coefficients ($C_{\ell\ell}, C_{\ell e}, C_{ee}$).
   • Significance: Statistical significance is estimated using the profile-likelihood formalism, incorporating a 1% systematic uncertainty on the background.

Key Contributions and Results

• Projected Sensitivity: The study projects sensitivities for the Wilson coefficients scaled by the new physics scale, $C/\Lambda^2$, reaching the range of $(0.6\text{--}1.6) \times 10^{-11} \text{ GeV}^{-2}$. This represents an improvement over current experimental limits by up to an order of magnitude, depending on the specific flavor and chiral structure of the operator.
• Energy and Polarization Dependence: The analysis demonstrates that the signal-to-background ratio improves significantly at higher energies (10 and 14 TeV) due to the energy-enhanced behavior of contact interactions, whereas SM backgrounds decrease. The inclusion of beam polarization (specifically $P_{\mu^-} = \pm 30\%$ and $\pm 80\%$) is shown to significantly enhance sensitivity and improve the resolution of correlations among Wilson coefficients.
• Correlation Resolution: A combined analysis across multiple center-of-mass energies and polarization states is shown to effectively disentangle the correlations between the different Wilson coefficients, which is difficult to achieve with a single energy setting.
• Kinematic Discrimination: The study confirms that hard $p_T$ cuts are the primary discriminator, suppressing SM backgrounds by more than an order of magnitude while retaining the majority of the LFV signal.

Significance
The paper establishes that a future multi-TeV muon collider possesses strong sensitivity to charged lepton flavour violating four-lepton interactions. By utilizing the energy growth of contact interactions and the clean environment of a muon collider, the proposed facility can probe the SMEFT parameter space for four-lepton operators with precision exceeding current low-energy bounds. The work highlights the utility of the optimal-observable technique in extracting Wilson coefficients and resolving their correlations, positioning the muon collider as a powerful probe for model-independent searches of new physics in the lepton sector.