$B$-meson decay width up to $1/m_b^3$ corrections within and beyond the Standard Model

This paper presents a complete calculation of $B$-meson decay widths up to $1/m_b^3$ corrections within and beyond the Standard Model by deriving analytic expressions for all matching coefficients of two-quark operators and previously missing weak-annihilation contributions, thereby finalizing the theoretical framework for non-leptonic $b$-quark decays relevant to $B$-meson lifetimes and addressing tensions in QCD factorization.

The Gist

Imagine you have a very heavy, unstable ball (a B-meson) sitting in a box. Eventually, this ball breaks apart into smaller pieces. Physicists want to know exactly how long it takes for this ball to break apart (its "lifetime").

For a long time, scientists have had a very good rulebook (the Standard Model) for predicting this. However, when they look at real experiments, the predictions are sometimes slightly off, like a clock that gains or loses a few seconds a day. This paper is about sharpening that rulebook to see if the clock is actually broken or if we just needed a better way to read it.

Here is a breakdown of what the authors did, using simple analogies:

1. The "Heavy Quark Expansion" (The Recipe Book)

To predict how long the ball lasts, the authors use a method called the Heavy Quark Expansion (HQE).

• The Analogy: Imagine trying to predict the exact path of a bowling ball rolling down a lane.
  • The Big Picture (Leading Order): First, you just look at the ball rolling straight. This is the easiest part and gives you a rough idea of the time.
  • The Details (Power Corrections): But the ball isn't perfect. It wobbles, it spins, and the lane isn't perfectly smooth. To get a precise prediction, you have to add corrections for these wobbles and spins.
  • The Paper's Job: The authors calculated the math for these "wobbles" and "spins" up to a very high level of detail (specifically, up to the third layer of corrections). Before this paper, some of these detailed corrections were missing or incomplete.

2. The "New Ingredients" (Beyond the Standard Model)

The Standard Model is like a standard recipe for a cake. But sometimes, the cake tastes a little different than the recipe says it should. Scientists suspect there might be "secret ingredients" (New Physics or BSM) mixed in that we haven't discovered yet.

• The Analogy: Imagine you are baking a cake, but you suspect someone might have secretly added a pinch of salt or a drop of vanilla that isn't in the official recipe.
• The Paper's Job: Instead of guessing what that secret ingredient is, the authors wrote down a Master Recipe. This Master Recipe includes every possible ingredient (Standard and non-Standard) that could theoretically be added. They then calculated exactly how each of these ingredients would change the baking time. This allows future scientists to look at the real cake and say, "Aha! The time is off by exactly this much, which means the secret ingredient must be this specific one."

3. Fixing the "Glitches" (Infrared Divergences)

When doing these complex calculations, the math sometimes hits a "glitch" where numbers blow up to infinity. In physics, this is called an infrared divergence.

• The Analogy: Imagine you are counting the number of people in a room, but the door is open and people are walking in and out so fast that your counter breaks.
• The Paper's Job: The authors found a specific type of glitch caused by "soft gluons" (tiny particles of force) being emitted by the lighter pieces of the broken ball. They realized that to fix the counter, they had to also account for a specific interaction called Weak Annihilation (where two particles inside the ball destroy each other).
  • The Result: They calculated this missing piece (the "Weak Annihilation" contribution) for the first time in this specific context. By adding this missing piece, the "glitch" disappears, and the math works perfectly. They even double-checked their work using two completely different mathematical tools (like measuring a room with a tape measure and then with a laser) to ensure the numbers matched.

4. The "Penguin" Surprise

In the world of particle physics, there are special particles called "Penguins" (named because of a joke, not because they look like birds). These are rare interactions that usually happen very quietly.

• The Analogy: Most of the time, the ball breaks apart because of the main ingredients. But sometimes, a tiny, rare "Penguin" interaction happens in the background.
• The Paper's Job: The authors also calculated how these "Penguin" interactions affect the lifetime, including how they mix with the main ingredients. While these effects are usually very small, the authors provided the precise math for them, ensuring that even the tiniest whispers of these interactions are accounted for in the final prediction.

Summary of the Achievement

Think of the prediction for the B-meson's lifetime as a high-precision clock.

• Before this paper: The clock was accurate to the minute, but the "seconds" and "milliseconds" were a bit fuzzy because some of the internal gears (the math for the wobbles and the "Weak Annihilation" piece) were missing or uncalculated.
• After this paper: The authors have built the missing gears and polished the existing ones. They have provided a complete, mathematically rigorous set of instructions (analytic expressions) for how the clock ticks, whether it follows the standard rules or if there are secret "New Physics" ingredients mixed in.

What they did NOT do:
They did not build a new machine, they did not find the secret ingredient yet, and they did not change the physical laws. They simply provided the perfectly detailed mathematical map that allows others to compare real-world experiments against theory with much higher precision. If the real clock still doesn't match this new, sharper map, then we will know for sure that there is a "secret ingredient" (New Physics) at play.

Technical Summary

Technical Summary: B-meson decay width up to $1/m_b^3$ corrections within and beyond the Standard Model

Problem Statement
The lifetimes of heavy hadrons containing a $b$ quark are fundamental observables for testing the Standard Model (SM) and constraining Beyond the Standard Model (BSM) physics. While experimental precision for $B$-meson lifetimes and lifetime ratios (e.g., $\tau(B^0_s)/\tau(B^0_d)$) has reached the per-mille level, theoretical predictions based on the Heavy Quark Expansion (HQE) currently possess uncertainties roughly 20 times larger. A specific motivation for this work arises from observed tensions between experimental branching ratios and theoretical predictions (based on QCD factorisation) for several colour-allowed non-leptonic $B$-meson decays. To address these tensions and constrain potential BSM effects in tree-level $b$-quark decays ($b \to q_1 \bar{q}_2 q_3$), a complete calculation of the HQE coefficients up to dimension-six is required. Prior to this work, the matching coefficients for two-quark operators at dimension-five and six, particularly those induced by QCD-penguin operators and generalised BSM operators, were incomplete. Furthermore, the treatment of infrared (IR) divergences in the Darwin operator coefficient, which mixes with four-quark operators, lacked a complete subtraction of weak-annihilation (WA) contributions.

Methodology
The authors compute the total decay width of a $B$ meson using the optical theorem, expressed as the imaginary part of the forward matrix element of the time-ordered product of two effective Hamiltonians. The calculation proceeds within the framework of the HQE, expanding the non-local operator in inverse powers of the heavy $b$-quark mass ($m_b$).

1. Effective Hamiltonian: The study starts from the most general model-independent effective Hamiltonian for non-leptonic $b \to q_1 \bar{q}_2 q_3$ decays, including all possible Dirac structures (current-current, scalar, pseudoscalar, tensor, and their chiral counterparts) and both SM and BSM Wilson coefficients.
2. Operator Expansion: The internal quark propagators are expanded in the external background gluon field using the Fock-Schwinger gauge. This expansion generates a tower of local operators. The calculation retains terms up to three covariant derivatives, corresponding to contributions up to dimension-six in the HQE.
3. Regularisation Schemes: To handle IR divergences arising from the emission of soft gluons from light-quark propagators (specifically affecting the Darwin operator at dimension-six), the authors employ two independent schemes:
   • Mass Regularisation: Working in $D=4$ dimensions with a finite light-quark mass ($m_q$).
   • Dimensional Regularisation (DR): Working in $D=4-2\epsilon$ with massless light quarks.
     The agreement between these two approaches serves as a crucial cross-check.
4. IR Divergence Cancellation: The calculation explicitly includes the matching of four-quark operators at dimension-six. Crucially, the authors compute the previously missing weak-annihilation (WA) contributions. These contributions are necessary to cancel the IR divergences appearing in the Darwin operator coefficients, ensuring finite physical results as $m_q \to 0$.
5. Computational Tools: The calculation utilizes tools such as qgraf for diagram generation, tapir for integral family identification, LiteRed for IBP reductions, and FeynCalc for Dirac algebra. The results are expressed in terms of master integrals with known imaginary parts.

Key Contributions and Results
The paper provides analytic expressions at Leading Order (LO) in QCD for all matching coefficients of two-quark operators in the HQE up to mass-dimension-six for the most general effective Hamiltonian.

• BSM Completeness: The results complete the calculation of BSM effects in non-leptonic, tree-level $b$-quark decays relevant for $B$-meson lifetimes. This includes the leading-power results at dimension-three and the coefficients of the chromomagnetic operator (dimension-five) and the Darwin operator (dimension-six).
• Weak Annihilation: The authors derive new analytic expressions for the WA contributions at dimension-six, which were previously missing and are essential for the consistent subtraction of IR divergences.
• QCD-Penguin Operators: As a by-product, the paper determines the SM matching coefficients up to dimension-six originating from QCD-penguin operators. This includes:
  • Interference terms between current-current and penguin operators.
  • Contributions quadratic in penguin operators.
  • New results for the chromomagnetic and Darwin operator coefficients induced by penguins, which are typically suppressed by $\alpha_s$ and $\alpha_s^2$ in the SM.
• Specific Channels: Explicit analytic results are provided for the CKM-dominant transitions $b \to c\bar{u}d$ and $b \to c\bar{c}s$, as well as the CKM-suppressed modes $b \to u\bar{c}s$ and $b \to u\bar{u}d$ (listed in the appendices).
• Consistency Checks: The calculation reproduces known dimension-three expressions and provides a non-trivial cross-check via the agreement between mass regularisation and dimensional regularisation.

Significance and Claims
The paper claims that these results "complete the calculation of BSM effects in non-leptonic, tree-level, $b$-quark decays relevant for $B$ meson lifetimes and lifetime ratios." By providing the full set of matching coefficients up to dimension-six, the work enables:

• Global Phenomenological Analyses: The results allow for the combination of lifetime observables with other precision flavour data (mixing parameters, CP asymmetries) to improve constraints on generic BSM Wilson coefficients.
• Resolution of Tensions: The framework offers a tool to investigate the origin of tensions observed in colour-allowed non-leptonic $B$-meson decays, distinguishing between missing SM contributions (e.g., QED corrections, rescattering) and potential BSM effects.
• Model Mapping: The generic operator basis can be mapped onto specific BSM scenarios (e.g., $W'$ models, diquarks, 2HDM) to test consistency with collider bounds.

The authors note that while their work provides the necessary perturbative coefficients, further theoretical improvements—specifically NLO-QCD corrections to the Darwin and chromomagnetic contributions and refined lattice QCD determinations of non-perturbative matrix elements—are required to fully enhance the sensitivity of lifetime observables to New Physics.