Constraints on SMEFT operators from $Z \to \mu \mu bb$ decay

This paper presents the first process-specific constraints on flavor-resolved four-fermion operators involving muons and bottom quarks by analyzing $Z \to \mu\mu bb$ decays within the Standard Model Effective Field Theory framework, utilizing state-of-the-art Monte Carlo simulations and profile likelihood methods to derive complementary limits on relevant Wilson coefficients.

The Gist

Imagine the universe is a giant, high-stakes game of billiards. The Standard Model is the rulebook that physicists have written down over the last 50 years. It tells us exactly how the balls (particles) should bounce, spin, and collide. For the most part, the game plays out exactly as the rulebook predicts.

But physicists suspect there's a "ghost player" or a hidden rule we haven't discovered yet—something heavy and invisible that influences the game from the shadows. We can't see this new player directly because they are too heavy to be caught in our current net.

This paper is about a clever new way to catch a glimpse of these hidden rules by watching a very specific, rare trick shot.

The Setup: The "Z-Boson" Trick Shot

In our billiard game, there's a special ball called the Z-boson. Usually, when this ball breaks apart, it splits into simple pairs, like two muons (a type of electron cousin) or two quarks. Physicists have watched these simple splits millions of times and checked the rulebook. Everything looks perfect.

But this paper focuses on a much rarer, more complex trick shot: The Z-boson splitting into four pieces at once—two muons and two bottom quarks (heavy cousins of the up/down quarks).

Think of it like this:

• The Standard Shot: A cue ball hits a target and splits into two balls. Easy to predict.
• The Rare Shot: The cue ball hits the target and explodes into four balls flying in a specific pattern. This is the $Z \to \mu\mu bb$ decay.

The Detective Work: The "SMEFT" Magnifying Glass

Since we can't see the "ghost player" (New Physics) directly, the authors use a tool called SMEFT (Standard Model Effective Field Theory).

Think of SMEFT as a magnifying glass with a grid overlay. It doesn't show you the ghost directly, but it allows you to measure tiny, subtle distortions in how the four balls fly apart. If the balls fly slightly differently than the rulebook predicts, it suggests the ghost player is nudging them.

The authors are looking for specific "nudges" caused by invisible forces that connect the light muons to the heavy bottom quarks.

The Simulation: A Digital Time Machine

Since we can't run this experiment a billion times in a real lab right now, the authors built a super-advanced video game simulation.

1. The Engine: They used software (MadGraph) to simulate the collision of protons at the Large Hadron Collider (LHC).
2. The Physics: They programmed the "ghost rules" (the new physics) into the game.
3. The Camera: They added a "camera" (Delphes) that mimics the real detectors at CERN, including the fact that the camera isn't perfect (it sometimes misses a ball or confuses a heavy ball for a light one).

They ran this simulation to see what the "perfect" game looks like, and then what the game looks like if the ghost player is interfering.

The Challenge: Finding a Needle in a Haystack

The problem is that this rare trick shot is incredibly hard to spot.

• The Haystack: There are billions of other collisions happening that look almost the same. For example, top quarks decaying can mimic this signal.
• The Needle: The real $Z \to \mu\mu bb$ signal.

To find the needle, the authors set up a series of filters (like a bouncer at a club):

• "Only let in events with exactly two muons."
• "Only let in events with two heavy bottom-quark jets."
• "Throw out anything with too much missing energy (which would mean a neutrino is hiding)."

After filtering out the noise, they looked at the invariant mass—essentially, the "weight" of the four particles combined. If the ghost player isn't there, the weight should cluster perfectly around the known mass of the Z-boson (like a bell curve). If the ghost player is there, the curve might get wider, shift, or change shape.

The Results: The First "Flavor-Specific" Clue

The authors found something exciting:

1. Validation: They confirmed that this rare four-body decay is a viable way to test the rules of physics.
2. New Limits: They placed the first-ever strict limits on how much "ghost interference" can happen specifically between muons and bottom quarks.
   • Previous studies looked at "average" effects across all particles.
   • This study says, "We are specifically watching the relationship between these two types of particles, and here is exactly how much they can deviate from the rulebook."

It's like saying, "We don't just know the average speed of all cars on the highway; we know exactly how fast the red sports cars are allowed to go before they break the law."

The Bottom Line

This paper is a proof of concept. It shows that by looking at these messy, complex, four-particle explosions, we can find new clues about the universe that simpler, two-particle experiments miss.

Even though this was a simulation (a digital experiment), it paves the way for real physicists at the LHC to start hunting for these specific signals. If they find a deviation in the future, it could be the first crack in the Standard Model, revealing the heavy, invisible physics that lies just beyond our current reach.

In short: They built a digital microscope to watch a rare, four-piece particle breakup, proving that this specific angle is a powerful new way to hunt for the universe's hidden secrets.

Technical Summary

1. Problem Statement and Motivation

The Standard Model Effective Field Theory (SMEFT) is the primary framework for probing indirect effects of heavy new physics through precision measurements. While SMEFT constraints have been extensively studied using purely leptonic $Z$ decays (e.g., $Z \to 4\ell$) and inclusive $Z$ production, mixed leptonic–hadronic decay modes remain largely unexplored.

• The Gap: Existing global fits often assume flavor universality or focus on specific sectors (purely leptonic or purely hadronic), potentially missing flavor-resolved effects.
• The Target: The decay channel $Z \to \mu\mu b\bar{b}$ is of special interest because:
  • The presence of bottom quarks enhances sensitivity to operators coupling to the third quark generation.
  • The dimuon final state allows for precise reconstruction of leptonic kinematics.
  • This channel probes four-fermion contact interactions between leptons and quarks ($\ell\ell bb$) and flavor-dependent modifications of $Z$–fermion couplings, which are weakly constrained by inclusive observables.

2. Methodology

A. Theoretical Framework

The authors analyze the process within the SMEFT framework, truncating the Lagrangian expansion at dimension-six operators:
$$\mathcal{L}_{\text{SMEFT}} = \mathcal{L}_{\text{SM}} + \sum_i \frac{C_i}{\Lambda^2} \mathcal{O}_i$$
They focus on six specific dimension-six operators that contribute to $Z \to \mu\mu b\bar{b}$:

1. Four-fermion operators: Specifically those involving second-generation leptons ($\ell_2$) and third-generation quarks ($q_3$), denoted as $C^{(1,3)}_{\ell q, 2233}$. These generate direct tree-level $\mu\mu b\bar{b}$ interactions.
2. Coupling modification operators: Operators modifying $Z$ couplings to leptons ($C^{(1,3)}_{H\ell, 22}$) and bottom quarks ($C^{(1,3)}_{Hq, 33}$).

The analysis retains both interference (linear in Wilson Coefficients, $C_i$) and quadratic ($C_i^2$) terms.

B. Simulation and Event Generation

• Process Generation: Signal ($Z \to \mu\mu b\bar{b}$) and backgrounds are generated at Leading Order (LO) using MadGraph5_aMC@NLO.
• Parton Showering & Hadronization: Performed using Pythia8.
• Detector Simulation: Delphes is used with a CMS-like detector configuration.
  • b-tagging: A medium working point (discriminator > 0.5) is applied, yielding ~70% efficiency for true $b$-jets and low mis-tag rates for $c$ and light jets.
• Reweighting Technique: Instead of generating new samples for every Wilson Coefficient (WC) value, the authors use the SMEFTsim 3.0 UFO model. Precomputed weights for different WC values are applied to SM events, allowing efficient exploration of the parameter space $N(\theta) = c + b\theta + a\theta^2$.

C. Event Selection and Backgrounds

The analysis simulates proton-proton collisions at $\sqrt{s} = 13$ TeV with an integrated luminosity of 138 fb$^{-1}$.

• Signal Definition: Two isolated muons and two $b$-tagged jets.

• Trigger: Single-muon ($p_T > 25$ GeV) or double-muon triggers ($p_{T1} > 18$ GeV, $p_{T2} > 7$ GeV).

• Selection Cuts:

  • $\ge 2$ muons (subleading $p_T > 10$ GeV) and $\ge 2$ jets ($p_T > 25$ GeV).
  • At least 2 $b$-tagged jets.
  • Missing Transverse Energy (MET) < 30 GeV (crucial for suppressing $t\bar{t}$ backgrounds which contain neutrinos).
  • Angular separation requirements ($\Delta R > 0.4$) between muons, jets, and muon-jet pairs.
  • Signal Region: Invariant mass of the $\mu\mu b\bar{b}$ system ($m_{\mu\mu b\bar{b}}$) restricted to 80–115 GeV (centered on the $Z$ mass).

• Dominant Backgrounds:

  1. $t\bar{t}$ (semileptonic): Suppressed significantly by the MET cut and $b$-tagging requirements.
  2. $Z + \text{jets}$: Includes light-flavor misidentification and heavy-flavor irreducible backgrounds. Modeled with MLM matching for up to 2 extra jets.
  3. Diboson ($ZZ$) and $ZH$: Subleading but topologically similar.

3. Key Results

A. Statistical Analysis

The authors construct an Asimov dataset assuming the Standard Model (SM) is the null hypothesis ($N_{\text{obs}} = N_{\text{SM}}$). They use a profile likelihood ratio test statistic to derive 95% Confidence Level (C.L.) intervals on the Wilson Coefficients ($C_i/\Lambda^2$).

B. Constraints on Wilson Coefficients

The study derives the first process-specific limits on flavor-resolved four-fermion operators involving muons and bottom quarks. The 95% C.L. intervals (in TeV$^{-2}$) for the current dataset (138 fb$^{-1}$) and the projected High-Luminosity LHC (HL-LHC, 3000 fb$^{-1}$) are:

Operator	138 fb$^{-1}$ Limit	HL-LHC (3000 fb$^{-1}$) Limit
$C^{(1)}_{\ell q, 2233}$	$[-0.023, 0.014]$	$[-0.005, 0.003]$
$C^{(3)}_{\ell q, 2233}$	$[-0.023, 0.014]$	$[-0.005, 0.003]$
$C^{(1)}_{Hq, 33}$	$[-0.029, 0.026]$	$[-0.006, 0.006]$
$C^{(3)}_{Hq, 33}$	$[-0.030, 0.026]$	$[-0.006, 0.006]$
$C^{(1)}_{H\ell, 22}$	$[-0.025, 0.009]$	$[-0.005, 0.002]$
$C^{(3)}_{H\ell, 22}$	$[-0.028, 0.025]$	$[-0.006, 0.005]$

• Quadratic Behavior: The analysis confirms that for four-fermion operators ($C_{\ell q}$), the constraints are dominated by the quadratic term (symmetric limits), whereas coupling modifiers ($C_{H\ell}, C_{Hq}$) show significant interference effects (asymmetric limits).
• Comparison: The constraints on coupling modifiers are comparable in magnitude to existing global SMEFT fits (e.g., CMS combined analyses), despite being derived from a single decay channel. Crucially, this work provides the first direct constraints on the flavor-specific four-fermion operators $C_{\ell q, 2233}$, which are typically marginalized in global fits.

4. Significance and Contributions

1. First Flavor-Resolved Constraints: This work provides the first process-specific limits on four-fermion operators coupling second-generation leptons to third-generation quarks ($\mu\mu b\bar{b}$), filling a gap in global SMEFT interpretations.
2. Complementarity: It demonstrates that mixed leptonic–hadronic $Z$ decays offer a unique probe for new physics directions that are weakly constrained by inclusive $Z$ observables or purely leptonic modes.
3. Methodological Validation: The study validates the use of the $Z \to \mu\mu b\bar{b}$ channel for SMEFT analyses, showing that despite QCD backgrounds and detector effects (like $b$-tagging), the signal can be reliably reconstructed and used to set competitive limits.
4. Future Prospects: The projected HL-LHC sensitivity shows a significant improvement (roughly a factor of 4-5 tighter constraints), establishing this channel as a vital component for future global SMEFT fits.

Conclusion

The paper establishes that the $Z \to \mu\mu b\bar{b}$ decay is a powerful and viable channel for probing SMEFT operators. By combining precise kinematic reconstruction with advanced reweighting techniques, the authors derive competitive constraints on both $Z$-coupling modifiers and, for the first time, flavor-specific four-fermion contact interactions, highlighting the necessity of including mixed leptonic–hadronic modes in future global new physics searches.