Decay matrix of B-\bar{B} mixing: Mixing of dimension-seven operators into dimension-six operators under renormalization

This paper establishes a methodology for handling the renormalization of dimension-seven operators in $B_s$-$\bar{B}_s$ mixing by demonstrating that their infrared-finite one-loop corrections in the $\overline{MS}$ scheme can be absorbed into finite counterterms proportional to dimension-six operators, thereby preserving power counting and constraining the definition of evanescent operators.

The Gist

The Big Picture: A Tug-of-War in the Subatomic World

Imagine the universe is filled with tiny, unstable particles called B-mesons. These particles are like a pair of dancers who are constantly switching places with their mirror-image partners (antiparticles). Sometimes they are the "light" dancer, sometimes the "heavy" dancer.

Physicists want to know exactly how long these dancers live before they disappear. There are two main ways they can vanish:

1. The Main Act: They vanish quickly and predictably (this is the "mass difference" part, which we already understand very well).
2. The Subtle Twist: They vanish at a slightly different rate depending on whether they are the "light" or "heavy" version. This difference is called the width difference ($\Delta\Gamma_s$).

The Problem: We can measure this "Subtle Twist" very precisely in the lab. But when physicists try to calculate it using math, the numbers don't quite match. The error bars on the math are too big (about 40% uncertainty).

Why? Because the math relies on a "recipe" called the Heavy Quark Expansion. Think of this recipe as a cake.

• Layer 1 (The Cake): The main ingredients (Dimension-6 operators). We know how to bake this perfectly.
• Layer 2 (The Frosting): Smaller corrections (Dimension-7 operators). These are the "power-suppressed" terms. They are supposed to be tiny, like a pinch of salt.

The Glitch: When we try to calculate the "Frosting" (the Dimension-7 terms) using Quantum Chromodynamics (QCD), the math goes haywire. The calculations produce "ghost" ingredients that shouldn't be there. Specifically, the math accidentally mixes the tiny "Frosting" ingredients back into the main "Cake" ingredients.

In physics terms, this is a renormalization problem. The "ghost" ingredients are mathematical artifacts that violate the rules of the recipe (power counting). They make the tiny corrections look huge, ruining the prediction.

The Solution: The "Ghost Catcher"

The authors of this paper, Artyom and Ulrich, act like culinary detectives or accountants. Their job is to find these "ghost" ingredients and cancel them out so the recipe works again.

Here is how they did it, using a simple analogy:

1. The "Evanescent" Ghosts

In their calculations, the physicists use a mathematical trick called Dimensional Regularization. Imagine trying to measure a 3D object, but you are forced to draw it on a 2D piece of paper. To make the math work, you have to pretend the paper has a tiny, invisible extra dimension (let's call it "epsilon").

In this invisible dimension, some operators (mathematical rules) appear that don't exist in our real 4D world. These are called Evanescent Operators. They are like "ghosts" that only exist in the math.

• The Issue: When the math is finished and we return to our real 4D world, these ghosts leave behind a residue. This residue looks like a "Dimension-6" ingredient (the main cake) mixed into our "Dimension-7" frosting.
• The Fix: The authors realized that if you define these ghosts very carefully (using specific rules called Fierz symmetry), you can calculate exactly how much "residue" they leave behind.

2. The "Counter-Term" (The Eraser)

Once they calculated the size of the ghost residue, they invented a Counter-Term.

• Think of the ghost residue as a smudge of ink on a clean white shirt.
• The Counter-Term is a special eraser designed to remove exactly that smudge.
• By adding this eraser to the recipe, the "Dimension-6" ghost disappears, and the "Dimension-7" term is left pure and tiny, just as it should be.

3. The "Large Color" Test

How do they know their eraser works? They used a super-powerful test called the Large $N_c$ Limit.

• Imagine the universe where the number of "colors" (a property of quarks) is infinite. In this imaginary world, the math becomes simple enough to solve exactly, like a puzzle with a clear solution.
• The authors calculated the "ghost residue" in this simple world. They found that their eraser (the counter-term) perfectly cleaned up the mess, restoring the correct order of magnitude.
• If the eraser didn't work, the "frosting" would still look like "cake," and the math would fail. Since it worked, they know their definitions of the "ghosts" (the evanescent operators) are correct.

Why Does This Matter?

This paper is a technical manual for fixing a broken part of the Standard Model of physics.

1. Precision: By fixing the "ghost" mixing, they provide the exact numbers (counter-terms) needed to calculate the B-meson decay width with much higher precision.
2. New Physics: If we can calculate the Standard Model prediction perfectly, and the experimental measurement still doesn't match, it means we have found New Physics. It could be a sign of undiscovered particles (like Supersymmetry) lurking in the shadows.
3. The Foundation: Before we can find new particles, we must be sure our old map is perfect. This paper draws the missing lines on that map.

Summary in One Sentence

The authors found a way to mathematically "erase" invisible errors that were making tiny corrections to particle decay look too big, ensuring that our predictions for how these particles behave are accurate enough to potentially reveal new secrets of the universe.

Technical Summary

1. Problem Statement

The paper addresses a critical bottleneck in the theoretical prediction of the width difference $\Delta\Gamma_s$ in the $B_s-\bar{B}_s$ system.

• Context: While the mass difference $\Delta M_s$ is well-understood, the width difference $\Delta\Gamma_s$ requires a precise calculation of the decay matrix element $\Gamma_{12}^s$.
• The Challenge: The Heavy Quark Expansion (HQE) for $\Gamma_{12}^s$ includes power-suppressed corrections of order $\mathcal{O}(\Lambda_{\text{QCD}}/m_b)$. These corrections involve matrix elements of dimension-7 four-quark operators ($R_j, \tilde{R}_j$).
• The Specific Issue: In the $\overline{\text{MS}}$ renormalization scheme, one-loop QCD corrections to dimension-7 operators generate terms that violate the expected power counting. Specifically, loop integrals involving derivatives on the light quark field ($s$) can yield contributions scaling as $m_b^0$ (dimension-6 behavior) rather than the expected $1/m_b$.
• Gap: There was no general methodology to handle the renormalization of these dimension-7 operators, specifically how to absorb these "wrong-power" terms into finite counterterms proportional to dimension-6 operators ($Q, \tilde{Q}_S$) to restore the correct HQE power counting.

2. Methodology

The authors employ perturbative QCD within the $\overline{\text{MS}}$ scheme using dimensional regularization ($D = 4 - 2\epsilon$).

• Operator Basis: The study focuses on the dimension-7 operators $R_2, R_3$ and their color-rearranged counterparts $\tilde{R}_2, \tilde{R}_3$, which mix into the dimension-6 basis operators $Q$ and $\tilde{Q}_S$.
• Renormalization Strategy:
  • They calculate the one-loop matrix elements of the dimension-7 operators.
  • They identify terms proportional to dimension-6 operators that arise from the UV region (hard loop momenta). These terms are infrared (IR) finite.
  • They introduce finite counterterms ($\delta R_j, \delta \tilde{R}_j$) proportional to dimension-6 operators to cancel these unwanted $m_b^0$ contributions, ensuring the renormalized matrix elements scale correctly as $\mathcal{O}(1/m_b)$.
• Evanescent Operators: A crucial part of the methodology involves the definition of evanescent operators (operators that vanish in $D=4$ but contribute in $D \neq 4$). The authors define evanescent operators for the dimension-7 basis ($E_1[R_2], E_2[R_2]$, etc.) with specific coefficients for the $\epsilon$-terms.
• Consistency Check ($N_c \to \infty$): To verify the correctness of their counterterms and the choice of evanescent operator definitions, the authors calculate the hadronic matrix elements in the limit of a large number of colors ($N_c \to \infty$). In this limit, matrix elements factorize into products of currents, allowing for an exact analytical check.

3. Key Contributions

1. Derivation of Finite Counterterms: The paper provides the first complete calculation of the finite renormalization constants ($Z^{(10)}$) required to mix dimension-7 operators into dimension-6 operators at Next-to-Leading Order (NLO) in $\alpha_s$.
2. Resolution of Power Counting Violation: It demonstrates that the power-counting violating terms in the $\overline{\text{MS}}$ scheme are indeed IR-finite and can be systematically removed via finite counterterms, thereby validating the consistency of the HQE at NLO.
3. Constraints on Evanescent Operators: The authors show that the choice of coefficients for the $\epsilon$-terms in the definitions of evanescent operators is not arbitrary.
   • For color-rearranged operators ($\tilde{R}_2, \tilde{R}_3$), the requirement that the renormalized matrix elements respect Fierz symmetry (invariance under Fierz transformations) uniquely fixes the parameters of the evanescent operators ($e a_2 = 12, e b_2 = 4, e b_3 = 10$).
   • Without this specific choice, the large-$N_c$ factorization check fails.
4. Calculation of Current Matrix Elements: As a by-product, the authors calculated the one-loop matrix elements of several $\bar{b}-s$ current operators involving derivatives and gluon field strength tensors, which are necessary for the large-$N_c$ checks.

4. Key Results

The paper presents explicit analytical expressions for the counterterms $\delta R_j$ and $\delta \tilde{R}_j$. These are summarized in Table 1 of the paper.

• Structure of Counterterms: The counterterms take the form:
  $$\delta R_j = \left[ \frac{Z^{(11)}}{\epsilon} + Z^{(10)} \right] Q + \left[ \frac{\tilde{Z}^{(11)}}{\epsilon} + \tilde{Z}^{(10)} \right] \tilde{Q}_S$$
  where $Z^{(10)}$ and $\tilde{Z}^{(10)}$ are the finite parts derived in the paper.
• Dependence on Scheme: The finite parts depend on the renormalization scale $\mu$, the number of colors $N_c$, and the Casimir factor $C_F$. Crucially, they depend on the parameters ($a_2, b_2, b_3$) defining the evanescent operators.
• Large-$N_c$ Verification:
  • For $R_2$ and $R_3$, the large-$N_c$ calculation confirms that the renormalized matrix elements scale as $\mathcal{O}(1/m_b)$ after applying the derived counterterms.
  • For $\tilde{R}_2$ and $\tilde{R}_3$, the verification only holds if the Fierz-symmetric definitions of the evanescent operators are used. This proves that Fierz symmetry is a necessary condition for the consistency of the renormalization scheme for these specific operators.

5. Significance

• Enabling NLO Precision: This work is a prerequisite for calculating the NLO QCD corrections to the $1/m_b$ terms in $\Gamma_{12}^s$. Without these counterterms, the theoretical prediction for $\Delta\Gamma_s$ cannot be improved beyond Leading Order (LO), limiting the ability to test the Standard Model against the increasingly precise experimental data (e.g., from LHCb).
• Lattice QCD Matching: The derived renormalization constants are essential for matching lattice QCD calculations (which compute matrix elements at a specific lattice spacing) to the continuum $\overline{\text{MS}}$ scheme. The mixing of dimension-7 into dimension-6 operators is a major source of uncertainty in lattice calculations; this paper provides the necessary perturbative input to control it.
• Theoretical Framework: It establishes a rigorous methodology for handling the renormalization of higher-dimensional operators in effective field theories where power counting is threatened by UV divergences, a problem relevant beyond just $B$-physics.

In summary, Hovhannisyan and Nierste have solved a long-standing technical hurdle in the Heavy Quark Expansion, providing the necessary tools to push the theoretical precision of $B_s$ mixing observables to the NLO level.