Inclusive $B$ decays, lifetimes, and $|V_{cb}|$: a unified heavy-quark expansion analysis

Technical Summary: Inclusive B decays, lifetimes, and |Vcb|: a unified heavy-quark expansion analysis

Problem Statement
The extraction of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element $|V_{cb}|$ from inclusive semileptonic $\bar{B} \to X_c \ell \bar{\nu}_\ell$ decays relies on the Heavy Quark Expansion (HQE). While global fits of spectral moments have historically achieved high precision, recent tensions between different moment measurements (specifically $q^2$, $E_\ell$, and $M_X$ moments) and the " $|V_{cb}|$ puzzle" (the discrepancy between inclusive and exclusive determinations) have raised questions regarding the robustness of current theoretical treatments.

Previous analyses typically addressed missing higher-order corrections (both perturbative $\alpha_s$ and power $1/m_b$) by introducing ad-hoc uncertainties on specific HQE matrix elements (e.g., varying $\rho_3^D$ and $\mu_G^2$ by 20–30%). The authors argue that as perturbative accuracy improves (reaching N$^3$LO), these heuristic uncertainties often dominate the error budget, masking the true sensitivity of data to higher-order power corrections and complicating the assessment of HQE convergence. Furthermore, the treatment of perturbative scale dependence has traditionally involved independent scale variations for each observable, failing to account for the fact that scale variations induce coherent shifts across all inclusive observables derived from the same differential decay rate.

Methodology
The authors propose a unified global analysis of inclusive semileptonic moments, partial branching fractions, and $B^+$/$B^0$ lifetimes using a novel fit strategy implemented in the Python package kolya. The methodology is characterized by three main pillars:

1. Correlated Nuisance Parameters for Scales: Instead of estimating perturbative uncertainties by varying renormalization scales ($\mu_{\alpha_s}, \mu_{m_b}^{kin}, \mu_{m_c}$) independently for each observable, the authors treat these unphysical scales as fully correlated nuisance parameters within the global fit. The scales are profiled over, allowing the data to constrain the preferred scale regions. This approach assumes that scale variations induce coherent shifts across all predictions because they originate from the same underlying triple-differential decay rate.
2. Fixed-Order Convergence Study: Rather than assigning ad-hoc theoretical uncertainties for missing higher-order power corrections, the authors perform global fits at fixed orders of the HQE, specifically from $\mathcal{O}(1/m_b^3)$ up to $\mathcal{O}(1/m_b^5)$. They assess the convergence of the HQE by observing the stability of the extracted parameters (such as $|V_{cb}|$ and HQE matrix elements) as higher-order terms are included. An additional uncertainty is estimated based on the maximum difference between the $\mathcal{O}(1/m_b^3)$ and $\mathcal{O}(1/m_b^{4,5})$ results.
3. Unified Lifetime and Semileptonic Analysis: The analysis extends beyond semileptonic observables to include $B$ meson lifetimes. The authors implement lifetime predictions (including non-leptonic contributions and four-quark operators) into kolya. They perform a simultaneous analysis where semileptonic moments constrain the HQE parameters and scales, which are then used to predict total widths and lifetimes. For non-leptonic contributions, the scale parameters are varied over a broader range to remain agnostic about the specific scale of the non-leptonic process.

Key Contributions

• New Fit Strategy: The introduction of perturbative scales as correlated nuisance parameters determined by the data, rather than varied over a fixed range, represents a shift in how theoretical uncertainties are handled in inclusive B physics.
• Order-by-Order Convergence Analysis: The paper provides a systematic study of HQE convergence up to $\mathcal{O}(1/m_b^5)$ without relying on ad-hoc uncertainty models, offering a direct test of the expansion's validity.
• Implementation in kolya: The authors extend the public kolya package to include:
  • Renormalization Group Evolution (RGE) of $\mu_G^2$.
  • Predictions for $B$ meson lifetimes including non-leptonic decays and four-quark operators.
  • Inclusive semileptonic rates for $\bar{B} \to X_c \tau \bar{\nu}_\tau$ and $\bar{B} \to X_u \ell \bar{\nu}_\ell$.
• Leave-One-Out Validation: A rigorous validation of the experimental dataset is performed to identify sources of tension, specifically highlighting the impact of specific datasets (e.g., Belle $q^2$ moments) on the global fit quality.

Results

• $|V_{cb}|$ Determination: In the default scenario (free scales, $\mathcal{O}(1/m_b^3)$), the authors find:
  $$|V_{cb}| = (42.32 \pm 0.28_{\text{fit}} \pm 0.10_{\text{h.o.}}) \times 10^{-3} = (42.32 \pm 0.30) \times 10^{-3}$$
  The central value is approximately $1\sigma$ higher than recent analyses using similar data but different fit strategies (e.g., Ref. [11]), with a smaller uncertainty.
• Branching Ratio: The inclusive semileptonic branching ratio is determined as:
  $$\mathcal{B}(\bar{B} \to X_c \ell \bar{\nu}_\ell) = (10.91 \pm 0.11_{\text{fit}} \pm 0.21_{\text{h.o.}})\%$$
• Fit Quality and Tensions: The fit to the full dataset yields a poor $p$-value ($8.4 \times 10^{-5}$). A leave-one-out study reveals that the tension is largely driven by the Belle $q^2$ moment measurements; excluding these moments improves the $p$-value to 0.2257. This suggests internal inconsistencies in the current experimental dataset or missing higher-order effects.
• HQE Parameters:
  • $|V_{cb}|$ and $\mu_\pi^2$ are found to be stable across orders $\mathcal{O}(1/m_b^{3,4,5})$.
  • $\mu_G^2(2.3 \text{ GeV})$ shows large variations (up to 43%) between orders, and the extracted value is significantly lower than in previous works, creating tension with the leading-order relation to the $B^*-B$ mass splitting.
  • $\rho_3^D(1 \text{ GeV})$ increases monotonically by about 13% per order.
• Lifetime Predictions: Using the extracted parameters, the authors predict:
  $$\tau_{B^+} = 1.89^{+0.04}_{-0.22} \text{ ps}, \quad \tau_{B^0} = 1.74^{+0.03}_{-0.23} \text{ ps}$$
  These predictions lie $1.1\text{--}1.2\sigma$ below experimental measurements. The large, asymmetric uncertainties are dominated by the perturbative scale variation of the non-leptonic width. However, the lifetime ratio $\tau_{B^+}/\tau_{B^0} = 1.082^{+0.031}_{-0.013}$ agrees well with experiment.

Significance and Claims
The paper claims that its fixed-order approach, free of ad-hoc power-correction uncertainties, provides a more rigorous test of the HQE's convergence and the consistency of the experimental data. The authors emphasize that the poor fit quality of the full dataset, particularly the tension with Belle $q^2$ moments, indicates that the current dataset cannot be described consistently within a fixed-order HQE analysis without additional theoretical inputs or new correlated measurements.

The authors assert that their treatment of scales as fitted nuisance parameters absorbs the perturbative truncation uncertainty into the central value, resulting in uncertainties that are not directly comparable to previous works where scales were varied. They conclude that new, correlated measurements of moments (especially $q^2$) and updated branching ratio measurements are essential to resolve the persisting $|V_{cb}|$ puzzle and the observed tensions in the data. The work also highlights the importance of using a consistent set of HQE parameters across different applications (semileptonic rates vs. lifetimes) due to the strong correlations between fitted parameters and scales.