New Physics in inclusive $\bar{B} \to X_c \ell \bar{\nu}$ decays

This paper presents the first global phenomenological fit of inclusive $\bar{B} \to X_c \ell \bar{\nu}$ decays to all available experimental data within the Weak Effective Theory, incorporating generic dimension-six New Physics and higher-order corrections, and finds no significant evidence for New Physics while establishing competitive bounds on the relevant Wilson coefficients.

The Gist

Imagine you are a detective trying to solve a mystery inside a tiny, chaotic factory called a B-meson. Inside this factory, heavy particles are constantly breaking apart and transforming into lighter ones. Specifically, we are watching a heavy "beauty" quark turn into a "charm" quark, while shooting out a lepton (like an electron) and a neutrino.

In the world of particle physics, this process is called $\bar{B} \to X_c \ell \bar{\nu}$.

The Mystery: Is the Standard Model Broken?

For decades, scientists have had a "rulebook" called the Standard Model that predicts exactly how these factories should behave. However, there are some clues in the data that don't quite fit the rulebook. It's like if you expected a car to get 30 miles per gallon, but it keeps getting 25. Is the car broken? Or is there a hidden mechanic (New Physics) tweaking the engine?

One of the biggest puzzles is a number called $|V_{cb}|$, which is essentially the "speed limit" for how fast this transformation happens. Different ways of measuring it give slightly different answers, creating a headache for physicists.

The Old Way vs. The New Way

Previously, scientists tried to solve this by looking at exclusive decays. Imagine this like watching a specific, single car drive out of the factory. You can see exactly where it goes, but you only see a tiny fraction of the traffic.

This paper introduces a new approach: inclusive decays. Instead of watching one specific car, the scientists look at the entire traffic jam leaving the factory. They count every possible way the beauty quark can turn into a charm quark, regardless of what other debris is left behind.

• The Analogy: If the exclusive method is counting specific cars, the inclusive method is weighing the entire pile of exhaust and scrap metal coming out of the factory. It's a much bigger sample size, which should make the measurements more precise.

The Detective Work: The "Global Fit"

The authors of this paper performed a massive statistical experiment called a global fit. Think of this as a giant puzzle where they tried to fit two different sets of pieces at the same time:

1. The "Engine" Pieces (QCD Parameters): These are the messy, internal details of the factory (how the quarks interact). These are hard to calculate because they involve the strong nuclear force, which is like trying to predict the weather inside a hurricane.
2. The "New Mechanic" Pieces (New Physics): These are the potential hidden rules (New Physics) that might be breaking the Standard Model.

They used a super-computer to adjust all these pieces simultaneously to see which combination best matched the real-world data collected by experiments like Belle, BaBar, and CDF.

The Twist: Power Corrections

To make their prediction accurate, the scientists had to account for "noise." In physics, this is called power corrections.

• The Analogy: Imagine you are trying to measure the speed of a car on a highway. The "main" speed is easy to calculate. But there are small bumps in the road, wind resistance, and tire friction that slightly alter the speed. The scientists calculated these tiny effects up to a very high level of precision (up to the third power of the noise), ensuring their "theoretical car" was as realistic as possible.

The Verdict: No New Mechanic Found (Yet)

After crunching the numbers, the results were surprisingly calm:

• The Standard Model Wins: The data fits the "Standard Model" rulebook almost perfectly. There is no strong evidence that a "New Mechanic" is secretly tweaking the engine.
• The "Flat" Spot: Interestingly, when they tried to find the "New Physics" numbers, the math got a bit wobbly. It's like trying to find the exact center of a very wide, flat valley; you know you're in the valley, but you can't pinpoint the exact spot. This means the data allows for some New Physics, but it has to be very small.
• The Limits: Even though they didn't find New Physics, they set very strict "speed limits" on how big any hidden New Physics could be. These limits are now just as tight as the ones found by looking at single cars (exclusive decays).

Why This Matters

This paper is a milestone because:

1. It's the First: It's the first time scientists have done this massive "inclusive" puzzle with New Physics included from the start.
2. It's Precise: They provided new, exact formulas for the "noise" (mathematical corrections), which makes future calculations easier and faster for everyone.
3. It's a Safety Net: By proving that the "traffic jam" method agrees with the "single car" method, they have strengthened our confidence in the Standard Model. If New Physics does exist, it's hiding very well, and we now know exactly how well it's hiding.

In short: The scientists looked at the entire mess of particle debris, accounted for every tiny bump and wind gust, and found that the universe is still behaving exactly as the old rulebook predicted. No new mechanics found, but we now have a much better map of the factory floor.

Technical Summary

1. Problem Statement

Semileptonic $b \to c \ell \bar{\nu}$ transitions are central to precision flavor physics. Despite their theoretical tractability via the Heavy Quark Expansion (HQE), long-standing discrepancies exist within the Standard Model (SM), specifically in the extraction of the CKM matrix element $|V_{cb}|$ and the ratio $R(D^{(*)})$. These anomalies, alongside others in the $b \to s \mu \mu$ sector, have fueled speculation regarding New Physics (NP).

While inclusive $\bar{B} \to X_c \ell \bar{\nu}$ observables have reached $\mathcal{O}(1\%)$ accuracy in the SM, previous phenomenological fits have not fully accounted for the impact of generic New Physics on the extraction of non-perturbative HQE parameters. Furthermore, inclusive and exclusive modes probe different combinations of Wilson Coefficients (WCs), yet a simultaneous global fit allowing for generic dimension-six NP interactions in inclusive modes had not been performed.

2. Methodology

The authors present the first global phenomenological fit of inclusive $\bar{B} \to X_c \ell \bar{\nu}$ observables to all available experimental data (from BaBar, CLEO, DELPHI, CDF, Belle, and Belle II) within a generic New Physics scenario.

Theoretical Framework:

• Effective Hamiltonian: The analysis utilizes the Weak Effective Theory with all relevant dimension-6 operators for $b \to c \ell \bar{\nu}$: Vector Left ($O_{VL}$), Vector Right ($O_{VR}$), Scalar ($O_S$), Pseudoscalar ($O_P$), and Tensor ($O_T$).
• Differential Rate: The decay rate is calculated using the optical theorem, relating the hadronic tensor to the forward-scattering amplitude.
• Operator Product Expansion (OPE): The forward-scattering matrix element is expanded systematically in powers of $\Lambda_{\text{QCD}}/m_b$. The calculation includes:
  • Power corrections: Up to $\mathcal{O}(\Lambda_{\text{QCD}}^3/m_b^3)$.
  • Perturbative QCD: Up to $\mathcal{O}(\alpha_s)$ (and $\mathcal{O}(\alpha_s^2)$ for specific SM moments).
  • HQE Parameters: Non-perturbative parameters ($\mu_\pi^2, \mu_G^2, \rho_{LS}^3, \rho_D^3$) are treated as free parameters to be fitted alongside NP Wilson Coefficients.

Fit Procedure:

• Observables: The fit utilizes normalized moments of the decay rate ($M_n = \langle X^n \rangle$) for lepton energy ($E_\ell$), hadronic mass squared ($m_{X_c}^2$), and momentum transfer squared ($q^2$), as well as the branching ratio fraction $R^*$. Central moments are used to minimize correlations.
• Parameterization:
  • In the SM, the fit extracts HQE parameters and $|V_{cb}|$ in a two-step process.
  • In the NP scenario, a parametric degeneracy arises between $|V_{cb}|$ and the Wilson Coefficient $C_{VL}$. To resolve this, the authors define a redefined parameter $\tilde{V}_{cb} \equiv V_{cb} C_{VL}$.
  • NP effects are parameterized using 7 variables: magnitudes ($a_{VR}, a_T, a_P, a_S$) and relative phases of the normalized WCs ($C_X/C_{VL} = a_X e^{i\delta_X}$).

3. Key Contributions

1. First Global NP Fit: This is the first simultaneous fit of NP Wilson Coefficients and HQE parameters using inclusive $\bar{B} \to X_c \ell \bar{\nu}$ data in a generic NP scenario.
2. Analytical $\mathcal{O}(\alpha_s)$ Corrections: The authors provide the first fully analytical expressions for the $\mathcal{O}(\alpha_s)$ QCD corrections to the first three kinematic moments ($E_\ell, m_{X_c}^2, q^2$) in the presence of New Physics. This significantly simplifies future phenomenological applications and improves computational efficiency.
3. Comprehensive Uncertainty Handling: The calculation includes power corrections up to $\mathcal{O}(\Lambda_{\text{QCD}}^3/m_b^3)$ and perturbative corrections up to $\mathcal{O}(\alpha_s)$ (and $\mathcal{O}(\alpha_s^2)$ for $q^2$ moments in the SM), ensuring state-of-the-art theoretical precision.

4. Results

• Goodness of Fit: The data is well-described by both the SM and the NP hypothesis.
  • SM fit: $\chi^2/\text{d.o.f.} = 42.6/74$.
  • NP fit: $\chi^2/\text{d.o.f.} = 36.1/67$.
• New Physics Preference: The NP hypothesis is compatible with the SM at the $0.7\sigma$ level. The fit finds no significant preference for New Physics; all fitted NP parameters are compatible with zero.
• Constraints on Wilson Coefficients: Despite the lack of NP signal, the fit establishes competitive upper bounds on the magnitude of the relevant Wilson Coefficients. These bounds are comparable to those derived from exclusive decay modes, validating inclusive decays as a robust probe for NP.
• Extraction of $|V_{cb}|$:
  • SM: The authors extract $|V_{cb}| = 41.64(47) \times 10^{-3}$, incorporating $\mathcal{O}(\alpha_s^2)$ corrections in $q^2$-moments and full $\mathcal{O}(\alpha_{em})$ corrections in the total rate.
  • NP: The redefined parameter is found to be $|\tilde{V}_{cb}| = (35.3^{+4.4}_{-3.6}) \times 10^{-3}$. The large uncertainty and low central value are attributed to approximate flat directions in the kinematic moments, which are absent in the total decay rate.

5. Significance and Outlook

This work establishes inclusive $\bar{B} \to X_c \ell \bar{\nu}$ decays as a powerful, model-independent probe for New Physics, capable of providing constraints competitive with exclusive modes. The provision of analytical $\mathcal{O}(\alpha_s)$ formulas for kinematic moments in the NP sector removes a significant barrier for future theoretical and phenomenological studies.

The authors identify a "flat direction" in the inclusive fit regarding $|\tilde{V}_{cb}|$ when NP is allowed. They suggest that this degeneracy can be resolved by:

1. Combining inclusive and exclusive decay modes in a global fit.
2. Measuring specific observables optimized for this purpose at Belle II.

This study represents a critical step toward a unified understanding of $b \to c \ell \bar{\nu}$ transitions and the potential discovery of physics beyond the Standard Model.