Update of the nonlocal sub-leading ${O}_1$-${O}_7$ contribution to $\bar B \to X_s γ$ at LO

This paper presents a revised calculation of the complete non-local sub-leading contribution to the inclusive decay $\bar B \to X_s \gamma$ by removing the previously subtracted local Voloshin term, thereby addressing the high correlation between uncertainties and significantly impacting the estimated range of this contribution.

The Gist

Here is an explanation of the paper using simple language and creative analogies.

The Big Picture: A Noisy Radio Station

Imagine you are trying to listen to a very specific radio station (a particle physics event called $\bar{B} \to X_s\gamma$) where a heavy particle (a Bottom quark) decays into a lighter one (a Strange quark) and emits a flash of light (a photon).

Physicists want to predict exactly how bright this flash should be. They have a "perfect" theoretical model (the Local OPE) that works like a crystal-clear radio signal. However, in the real world, there is static and interference.

In this paper, the authors are focusing on a specific type of interference called "resolved contributions." Think of this as a ghostly echo. Instead of the light coming directly from the main event, it gets "bounced" off other particles (like charm quarks) inside the atom before it reaches you. This echo is messy, hard to calculate, and creates a lot of uncertainty in the final prediction.

The Problem: The "Subtracted" Mistake

For years, physicists tried to calculate this messy echo. However, they made a specific choice: they calculated the "echo" part and then subtracted a specific, well-understood piece of it (called the Voloshin term) because they thought it was already accounted for in their main model.

The Analogy:
Imagine you are trying to measure the weight of a backpack.

1. You weigh the whole backpack.
2. You realize there is a heavy, known rock inside.
3. You decide to subtract the weight of the rock to find the weight of the "rest" of the backpack.
4. The Mistake: You then try to estimate the uncertainty of the "rest" of the backpack without considering that your estimate of the rock's weight was also shaky.

In the old calculations, the authors treated the "rock" (Voloshin term) and the "rest of the backpack" (the non-local shape function) as two completely separate, unrelated things. They calculated the uncertainty for the "rest" and added a fixed number for the "rock."

The Reality:
The authors of this paper realized that the "rock" and the "rest of the backpack" are actually tightly linked. If you change how you estimate the rock, it changes how you estimate the rest. They are correlated. By treating them separately, the old calculations were underestimating the total messiness (uncertainty).

The New Calculation: Putting the Puzzle Back Together

The authors (Benzke, Garzelli, and Hurth) decided to stop subtracting. Instead, they calculated the entire messy echo at once, treating the "rock" and the "rest" as a single, unified system.

How they did it:

1. The Shape Function: They used a mathematical "net" (a set of basis functions like Hermite polynomials) to catch all possible shapes this messy echo could take. They didn't guess the shape; they let the math explore every possibility allowed by the laws of physics.
2. The Variables: They ran thousands of simulations, changing the "knobs" of the universe:
   • How heavy is the charm quark? (The "rock")
   • How heavy is the bottom quark?
   • How "fuzzy" is the internal structure?
3. The Result: When they put the rock and the rest back together, the range of possible outcomes got much wider.

The Findings: The Uncertainty Grows

In the old days, they thought the uncertainty for this specific effect was roughly 2.9% to 8%.

After their new, more honest calculation (which includes the correlation between the parts):

• The range expanded to 2.6% to 13% (when accounting for scale ambiguities).

Why does this matter?
In particle physics, if you want to find "New Physics" (like a new particle that breaks the Standard Model), you need your predictions to be incredibly precise. If your prediction has a huge "fuzzy zone" (uncertainty), you can't tell if a new discovery is real or just a fluctuation in the noise.

This paper says: "Hey, our previous estimate of the noise was too optimistic. The noise is actually much louder and more unpredictable than we thought."

The "Scale" Ambiguity: Tuning the Radio

The paper also mentions a "scale ambiguity."
Analogy: Imagine you are measuring a room with a ruler. If you measure in inches, you get one number. If you measure in centimeters, you get another. In quantum physics, the "ruler" (the energy scale) isn't fixed at the lowest level of calculation.

• The authors showed that if you change the ruler from "hard" (short distance) to "hard-collinear" (slightly longer distance), the uncertainty range gets even bigger.
• They are currently working on a more advanced calculation (Next-to-Leading Order, or NLO) that will fix the ruler, making the measurement more precise in the future.

Summary in One Sentence

This paper fixes a calculation error where physicists were treating two connected parts of a messy quantum effect as separate things, revealing that the true uncertainty in predicting how heavy particles decay is significantly larger than previously thought.

Why Should You Care?

Even though this is about subatomic particles, it represents the scientific method in action. Scientists are constantly refining their tools. By admitting that their previous "best guess" was too confident, they are actually making the field stronger. It ensures that when they do find something new in the future, they won't mistake a calculation error for a discovery.

Technical Summary

Here is a detailed technical summary of the paper "Update of the nonlocal sub-leading O1 - O7 contribution to $\bar{B} \to X_s\gamma$ at LO".

1. Problem Statement

The paper addresses a specific theoretical uncertainty in the inclusive radiative penguin decay $\bar{B} \to X_s\gamma$. This decay is a crucial probe for New Physics, but its theoretical prediction is limited by non-local sub-leading power corrections (resolved contributions).

• The Specific Issue: Previous calculations of the interference between the electroweak operators $O_1^c$ and $O_{7\gamma}$ treated the contribution as two separate parts:
  1. A non-local shape function term ($\Lambda_{Shape}^{17}$), calculated using Soft Collinear Effective Theory (SCET).
  2. A local Voloshin term ($\Lambda_{Voloshin}^{17}$), which was traditionally subtracted from the non-local calculation.
• The Flaw: In prior analyses (specifically Ref. [9]), the Voloshin term was subtracted, and its uncertainty was not included in the final error budget. However, the authors argue that this separation is unnatural from the SCET perspective. The uncertainties of the local Voloshin term and the non-local shape function term are highly correlated. Treating them separately and ignoring the Voloshin uncertainty leads to an underestimation of the total theoretical error.
• Goal: To present a calculation of the complete non-local contribution (combining both terms) to provide a more accurate estimate of the uncertainty range at Leading Order (LO).

2. Methodology

The authors employ a systematic approach based on Soft Collinear Effective Theory (SCET) and Heavy Quark Effective Theory (HQET).

• Factorization Framework: The differential decay rate is described by a factorization formula involving a hard function ($H$), jet functions ($J, \hat{J}$), and a non-perturbative shape function ($g$). The relevant contribution is the convolution of perturbative jet functions with the non-perturbative shape function $h_{17}$.
• Unified Calculation: Instead of separating the Voloshin term, the authors calculate the integral of the complete expression:
  $$\Lambda_{17}^{Complete} = \Lambda_{Shape}^{17} + \Lambda_{Voloshin}^{17}$$
  This corresponds to integrating the full kernel (including the penguin function $F(x)$ and the local term) against the shape function.
• Systematic Shape Function Analysis:
  • They utilize the method proposed in Ref. [14], using a complete set of basis functions (Hermite polynomials multiplied by a Gaussian) to model the unknown functional form of the shape function $h_{17}(\omega_1)$.
  • The model is constrained by known physical properties:
    • PT Invariance: The function is even.
    • Moments: The zeroth moment (related to $\lambda_2$), second moment, and dimensional estimates for the 4th and 6th moments are used as constraints.
  • Grid Scan: A comprehensive numerical scan is performed over all input parameters within their uncertainties:
    • Charm mass ($m_c$): $1.14 - 1.23$ GeV.
    • Bottom mass ($m_b$): $4.55 - 4.61$ GeV.
    • Shape function moments ($m_0, m_2, m_4, m_6$).
    • Model parameter $\sigma$ (Gaussian width).
• Scale Ambiguity: To manifest the LO scale ambiguity, the authors vary the renormalization scale of the Wilson coefficients from the hard scale ($\mu_H \sim m_b$) to the hard-collinear scale ($\mu_J \sim \sqrt{m_b \Lambda_{QCD}}$).

3. Key Contributions

• Removal of Artificial Subtraction: The paper demonstrates that subtracting the Voloshin term and ignoring its uncertainty is theoretically inconsistent. By calculating the complete contribution, they eliminate the artificial separation.
• Correlated Uncertainty Analysis: The study explicitly accounts for the high correlation between the local and non-local terms, which was previously neglected.
• Revised Error Budget: The authors provide a new, expanded range for the resolved contribution, showing that previous estimates were too optimistic regarding precision.

4. Results

The numerical analysis yields the following results for the normalized contribution $F_{17}^{Complete}$:

• Central Range (Fixed Scale):
  • Previous estimate (Shape term only + Voloshin central value): $[2.9\%, 8.0\%]$.
  • New Complete Calculation: $F_{17}^{Complete} \in [2.6\%, 8.9\%]$.
  • This represents a 22.5% increase in the uncertainty range compared to the previous analysis.
• Including Scale Ambiguity:
  • When the scale of the Wilson coefficients is shifted to the hard-collinear scale to reflect LO ambiguity, the range expands further.
  • Final LO Result: $F_{17}^{Complete} \in [2.6\%, 13\%]$.
• Parameter Sensitivity:
  • The extreme values in the range are driven by the smallest charm mass ($m_c \approx 1.14$ GeV) and the largest bottom mass ($m_b \approx 4.61$ GeV).
  • The increase in the range is primarily attributed to the inclusion of the uncertainties of the Voloshin term, which were previously omitted.
  • The range is non-Gaussian, representing a confidence interval within which the true value is expected to lie, rather than a statistical standard deviation.

5. Significance

• Impact on Phenomenology: The resolved contribution due to $O_1^c - O_{7\gamma}$ interference is one of the largest theoretical uncertainties in the inclusive $\bar{B} \to X_s\gamma$ decay rate. The updated range ($2.6% - 13%$) significantly impacts the precision of Standard Model predictions and the sensitivity to New Physics searches.
• Theoretical Consistency: The work corrects a methodological flaw in previous literature by treating the local and non-local parts as a single, correlated entity, adhering to the principles of SCET factorization.
• Future Directions: The authors note that the large scale ambiguity and charm mass dependence are artifacts of the Leading Order calculation. They highlight that ongoing Next-to-Leading Order (NLO) calculations (specifically $\alpha_s$ corrections and RG resummation) are necessary to:
  1. Establish a well-defined scale dependence (reducing the scale ambiguity).
  2. Potentially reduce the large dependence on the charm quark mass.

In summary, this paper provides a mandatory update to the theoretical error budget for $\bar{B} \to X_s\gamma$, demonstrating that the uncertainty is significantly larger than previously thought when the local Voloshin term is treated correctly and correlated with the non-local shape function.