Next-to-next-to-leading QCD corrections to the $\mathbf{B^+}$-$\mathbf{B_d^0}$, $\mathbf{D^+}$-$\mathbf{D^0}$, and $\mathbf{D_s^+}$-$\mathbf{D^0}$ lifetime ratios

This paper presents next-to-next-to-leading QCD corrections to the lifetime ratios of $B^+$, $D^+$, $D_s^+$, and $D^0$ mesons, combining three-loop perturbative calculations with hadronic matrix elements to produce theoretical predictions that show good agreement with experimental data.

The Gist

Imagine the universe as a giant, complex machine made of tiny, invisible building blocks called quarks. Some of these blocks are heavy and slow, like a bowling ball (the "b" and "c" quarks), while others are light and fast, like ping-pong balls. When these heavy blocks form particles called "mesons" (like the $B$ and $D$ mesons), they don't last forever; they eventually decay, or fall apart, into lighter particles.

The main question this paper answers is: How long do these heavy particles live, and why do some live slightly longer than their "twins"?

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

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

To predict how long a particle lives, physicists use a method called the Heavy Quark Expansion (HQE). Think of this like a recipe for a cake.

• The Main Ingredient: The most important part of the recipe is the heavy quark itself. If you just look at this, all heavy particles should have the exact same "lifetime" (how long the cake lasts before it crumbles).
• The Secret Spices: However, in reality, some particles live a tiny bit longer or shorter. This is because of the "spices" mixed in—interactions with the other, lighter quarks inside the particle.
• The Hierarchy: The recipe says the main ingredient is the biggest factor. The spices are smaller factors. The paper focuses on the third layer of spices (mathematically called terms suppressed by $1/m^3$). These are the specific interactions that cause the differences in lifetimes between particles that look almost identical.

2. The Problem: The "Three-Loop" Puzzle

Calculating these "spice" interactions is incredibly hard. It involves solving complex mathematical puzzles involving quantum mechanics.

• Previous Attempts: Before this paper, scientists had calculated the first and second layers of complexity (called Leading Order and Next-to-Leading Order). It was like trying to bake a cake with a blurry recipe; the results were close, but not precise enough to match the ultra-accurate measurements taken in modern labs.
• The New Achievement: This team calculated the third layer of complexity (Next-to-Next-to-Leading Order, or NNLO). In the language of Feynman diagrams (the maps physicists use to draw particle interactions), this required solving three-loop calculations.
  • Analogy: If the previous calculations were like drawing a map with a pencil, this paper drew the map with a laser, accounting for every tiny twist and turn in the quantum world that was previously ignored.

3. The Twins: $B$ and $D$ Mesons

The authors looked at two specific pairs of "twins":

• The $B$ Mesons: A charged one ($B^+$) and a neutral one ($B^0_d$).
• The $D$ Mesons: A charged one ($D^+$), a neutral one ($D^0$), and a strange one ($D^+_s$).

In the world of particle physics, these twins are almost identical, but they have different "flavors" of light quarks attached to them. The paper calculates exactly how much longer the charged version lives compared to the neutral version.

4. The Results: A Perfect Match

The team combined their new, ultra-precise mathematical "recipe" with data from other methods (like "Lattice QCD," which is like running a supercomputer simulation of the particle's interior).

• For the $B$ Mesons: They predicted the ratio of lifetimes to be 1.072. The actual experiment measured 1.076.
  • The Verdict: This is a perfect match. The difference is so small it's within the margin of error. This proves that their "recipe" (the Heavy Quark Expansion) is working correctly and that the "spices" they calculated are the right ones.
• For the $D$ Mesons: They predicted ratios of 2.344 and 1.289. The experimental values are 2.510 and 1.222.
  • The Verdict: These are also in good agreement, though the $D$ mesons are a bit trickier because they are lighter and the "spices" are a bit messier. The small differences between their prediction and the experiment help scientists estimate how much "noise" comes from even smaller, higher-order effects they haven't calculated yet.

5. Why This Matters

Think of this paper as a calibration check for the entire field of heavy particle physics.

• Validation: By showing that their complex math matches the real-world measurements so well, they confirmed that the Heavy Quark Expansion is a reliable tool.
• The "Unknowns": Because their prediction matches the experiment so well, they can now confidently say that any remaining tiny differences must come from effects they haven't calculated yet (like the "fourth layer of spices"). This helps them estimate the size of those unknown effects without needing to calculate them immediately.
• Future Safety: Since this method works so well for these "boring" particles (where we know the answer), scientists can now use this same method to study "exotic" particles where we don't know the answer yet, looking for signs of new physics beyond our current understanding.

In short: The authors built a super-precise mathematical model to explain why heavy particles live slightly different amounts of time. They tested it against real data, and it passed with flying colors, proving their model is solid and ready to be used for even more complex mysteries in the universe.

Technical Summary

1. Problem Statement

The lifetimes of heavy hadrons ($H_Q$) are calculated using the Heavy Quark Expansion (HQE), an operator product expansion where the total decay width $\Gamma(H_Q)$ is expanded in powers of $1/m_Q$ and the strong coupling $\alpha_s(m_Q)$.

• The Challenge: While the leading term (dimension-3) is universal for all hadrons containing a specific heavy quark, lifetime differences (splittings) arise from higher-dimensional operators (dimension-6 and above) where the spectator quark participates in the weak decay.
• Current Status: For $B$-mesons, the expansion parameter $\Lambda_{\text{QCD}}/m_b \approx 0.1$ suggests good convergence. For $D$-mesons, $\Lambda_{\text{QCD}}/m_c \approx 0.3$ implies slower convergence and larger theoretical uncertainties.
• The Gap: Previous theoretical predictions for lifetime ratios (e.g., $\tau(B^+)/\tau(B^0_d)$) were limited to Next-to-Leading Order (NLO) in QCD. Given the high precision of modern experimental data and the availability of new non-perturbative inputs (Lattice QCD and QCD Sum Rules), NLO precision is insufficient to rigorously test the HQE or constrain Beyond-Standard-Model (BSM) physics. There was a lack of Next-to-Next-to-Leading Order (NNLO) calculations for the $\Delta B = 0$ Wilson coefficients governing these lifetime splittings.

2. Methodology

The authors perform a systematic calculation of the NNLO QCD corrections to the Wilson coefficients of the dimension-6 operators in the HQE.

• Theoretical Framework:

  • Optical Theorem: The total decay width is related to the forward scattering matrix element $\langle H_Q | \text{Im} \, i \int d^4x \, T \{ \mathcal{H}_{\text{eff}}(x) \mathcal{H}_{\text{eff}}(0) \} | H_Q \rangle$.
  • Effective Theories: The calculation involves matching the full theory (with $|\Delta B|=1$ effective Hamiltonian) to an effective theory with $\Delta B = 0$ local operators.
  • Operators: The dominant contributions to lifetime splittings come from Weak Annihilation (WA) and Pauli Interference (PI) diagrams, represented by four-quark operators of dimension-6.
  • Symmetry Limits: Calculations are performed in the limit of exact Isospin symmetry (for $B$ and $D$ systems) and V-spin symmetry (for $D_s$ vs $D^0$), neglecting isospin-breaking effects in matrix elements.

• Computational Techniques:

  • Loop Order: The calculation involves a three-loop computation on the $|\Delta B|=1$ side and a two-loop computation on the $\Delta B = 0$ side.
  • Tools: The authors used independent computer codes in FORM and Mathematica. Diagram generation was done with qgraf, amplitude processing with tapir, and Integration-by-Parts (IBP) reduction to master integrals using Kira.
  • Renormalization:
    • Dimensional regularization with anticommuting $\gamma_5$.
    • $\overline{\text{MS}}$ scheme for $\alpha_s$, charm mass, and matching coefficients; Pole scheme for the bottom mass (later converted).
    • Careful handling of evanescent operators to ensure Fierz symmetry and correct renormalization group evolution.
  • CKM Structure: The calculation includes NNLO corrections to the CKM-favored contributions and NLO corrections to CKM-suppressed terms (involving $|V_{ub}|^2$, $|V_{cd}|^2$, etc.).

• Non-Perturbative Input:

  • The theoretical predictions are combined with hadronic matrix elements (Bag parameters) calculated via:
    1. QCD Sum Rules (HQET sum rules).
    2. Lattice QCD (recent complete calculations).

3. Key Contributions

1. First NNLO Calculation for $\Delta B = 0$: This is the first computation of NNLO QCD corrections to the Wilson coefficients governing heavy meson lifetime ratios.
2. CKM-Suppressed Terms: The authors provide new NLO corrections to the Cabibbo-suppressed contributions in the $B^+-B^0_d$ ratio, which were previously neglected or approximated.
3. Application to Charm System: The results are adapted to the $D$-meson system, providing NNLO predictions for $\tau(D^+)/\tau(D^0)$ and $\tau(D^+_s)/\tau(D^0)$, including V-spin breaking effects.
4. Validation of HQE: By comparing high-precision theory with experiment, the paper provides a stringent test of the convergence of the HQE.

4. Results

A. $B$-Meson System ($\tau(B^+)/\tau(B^0_d)$)

• Prediction: Combining NNLO coefficients with Bag parameters from QCD Sum Rules:
  $$\tau(B^+)/\tau(B^0_d) = 1.072 \pm 0.024$$
• Comparison: This agrees excellently with the experimental value of $1.076 \pm 0.004$.
• Implications:
  • The agreement validates the HQE framework.
  • It suggests that contributions from dimension-7 operators ($1/m_b^4$ terms) are small.
  • The reduction in renormalization scale dependence from NLO to NNLO is significant, though the scale uncertainty remains the dominant theoretical error source.

B. $D$-Meson System

Due to the larger $\Lambda_{\text{QCD}}/m_c$, the theoretical uncertainties are larger, and the choice of non-perturbative input (Sum Rules vs. Lattice) has a dramatic effect.

1. Ratio $\tau(D^+)/\tau(D^0)$:

• Using Sum Rules: Prediction $\approx 2.344 \pm 0.170$.
• Using Lattice QCD: Prediction $\approx 2.344 \pm 0.170$ (Note: The central value is similar, but the error budget shifts).
• Comparison: Experimental value is $2.510 \pm 0.015$.
• Analysis: The discrepancy between theory and experiment implies that dimension-7 ($1/m_c^4$) terms contribute roughly 10% (using Lattice inputs) to 50% (using Sum Rules) of the total splitting. This highlights the sensitivity of the charm system to higher-order HQE terms.

2. Ratio $\tau(D^+_s)/\tau(D^0)$:

• Prediction (Lattice inputs): $\tau(D^+_s)/\tau(D^0) = 1.289 \pm 0.042$.
• Comparison: Experimental value is $1.222 \pm 0.006$.
• Analysis: The difference is attributed to:
  1. Unaccounted $1/m_c^4$ terms.
  2. V-spin breaking: The $D^+_s$ and $D^0$ are V-spin partners, but the strange quark mass breaks this symmetry ($m_s - m_u \sim 0.3 \Lambda_{\text{QCD}}$). The calculation accounts for this via decay constant ratios ($f_{D_s}/f_D$), but residual breaking effects remain.

5. Significance

• Precision Physics: This work represents a major step forward in precision flavor physics, bringing theoretical predictions for heavy hadron lifetimes to the same order of accuracy as experimental measurements.
• Validation of HQE: The successful agreement for $B$-mesons confirms that the HQE is a robust tool for calculating lifetime differences, even for the charm system where convergence is slower.
• BSM Constraints: Since lifetime ratios are largely insensitive to BSM physics (dominated by SM weak decays), the agreement between theory and experiment sets a baseline. Any future deviation could signal new physics.
• Future Directions: The paper provides the necessary perturbative coefficients to allow future improvements in the determination of hadronic matrix elements (via Lattice QCD) to directly reduce the total theoretical uncertainty. It also quantifies the size of neglected higher-order terms ($1/m_Q^4$), guiding future theoretical efforts.

In summary, the paper establishes a new standard for theoretical precision in heavy meson lifetime calculations, confirming the validity of the HQE while quantifying the specific limitations and higher-order effects in the charm sector.