New Global Fit Sharpens the CKM Matrix Element Vcb and Keeps the R(D) Puzzle Alive

A team of theoretical physicists has delivered one of the most precise exclusive determinations to date of the CKM matrix element |Vcb|, the number that governs how strongly a bottom quark transforms into a charm quark in the weak interaction. Writing in The European Physical Journal C, Xue-Wen Chen, Bo-Yan Cui, Jia-Wei Zhang, Xue-Chen Zhao and Ya-Hui Chen combined the latest Belle II measurements, state-of-the-art lattice-QCD calculations of the B-to-D and B-to-D* form factors, and correlated light-cone sum-rule predictions at large hadronic recoil into a single global fit. Their headline result is |Vcb| = (39.18 ± 0.47) × 10⁻³, a value with a smaller overall uncertainty than many previous exclusive extractions, and one that comes with a full accounting of the correlations among all fitted parameters.

The CKM matrix sits at the heart of the Standard Model’s description of quark flavour mixing. Its element |Vcb| quantifies the charged-current transition strength between the charm and bottom quarks, and it plays a double role in particle physics: it is used to verify the unitarity of the CKM matrix, and it serves as a benchmark against which possible new physics beyond the Standard Model can be tested. A precise determination of |Vcb| has therefore long been one of the central goals of B-physics phenomenology. What makes the problem genuinely difficult is a long-standing discrepancy between exclusive extractions, which rely on specific decays such as B → D(*)ℓν, and inclusive extractions, which sum over many final states. That tension has motivated extensive theoretical and experimental work for years, and it remains unresolved.

The experimental side of the analysis draws on a decade of increasingly sophisticated measurements. The Belle collaboration measured the semileptonic decay B → Dℓν using fully reconstructed B meson tag events, obtained the isospin-averaged branching fraction, and extracted |Vcb| with two different form-factor parameterization frameworks, presenting binned differential decay rates in bins of the recoil variable w. Belle II then performed an analogous measurement adopting the BCL parameterization. For the vector-vector channel B⁰ → D*⁻ℓ⁺ν, independent fits using the CLN and BGL parameterizations have been carried out alongside full binned differential distributions, response matrices, and statistical and systematic correlation matrices, all designed to alleviate the inclusive–exclusive |Vcb| puzzle. Belle II has also provided separate partial decay rates for electron and muon final states as functions of the recoil parameter ω and three decay angular observables.

These experimental fits require form factors evaluated near zero recoil as external theoretical inputs, and lattice QCD has undergone steady methodological and precision improvements over the past decade. The first unquenched lattice gauge simulation of the B → Dℓν form factors came from the Fermilab-MILC collaborations. HPQCD then produced full-q² unquenched calculations of the scalar and vector form factors f₊(q²) and f₀(q²). For the exclusive mode B → Dℓν, unquenched lattice-QCD form factors at nonzero recoil were computed by Fermilab-MILC with 2+1 flavours, and HPQCD later delivered complete vector, axial-vector, and tensor form factors covering the entire q² range. JLQCD introduced Möbius domain-wall fermion discretization for the B → Dℓν transition, offering an independent cross-check against the widely used clover heavy-quark schemes. Collectively, these efforts trace a trajectory from narrow zero-recoil single-channel studies toward full-q², multi-flavour, multi-operator form factor sets with multiple independent discretization frameworks, drastically reducing lattice systematic uncertainties in the |Vcb| extraction.

The theoretical backbone of the new analysis is the BGL dispersive parameterization, a model-independent z-expansion scheme in which the form factors are written as series in a conformal variable z(q²) that maps the physical kinematic region into the unit disk. Blaschke factors remove poles associated with Bc resonances below the pair-production threshold, and outer functions φᵢ, built from perturbative QCD dispersion relations, control the series convergence. The authors impose strong unitarity bounds derived from four heavy-to-heavy transition channels — B → D, B → D*, B* → D and B* → D* — covering twenty form factors in total, with heavy-quark effective theory power corrections retained through O(αs), O(εb), O(εc) and O(εc²) terms. These bounds act as hard constraints in the fit, effectively restricting the space of BGL expansion coefficients and suppressing uncontrolled systematic uncertainties from z-series extrapolation. Sensitivity tests varying the HQET inputs over their quoted ranges shift the extracted |Vcb| by at most about 0.16σ of the nominal uncertainty.

A key innovation is the treatment of the light-cone sum-rule (LCSR) inputs at large hadronic recoil. The team used correlated SCET sum-rule predictions including next-to-leading-order QCD corrections and four sources of next-to-leading-power contributions, evaluated at six representative kinematic points q² ∈ {−3, −2, −1, 0, 1, 2} GeV². Rather than treating these predictions as independent numbers, the authors generated 1000 complete input-parameter sets by sampling the underlying LCSR parameters simultaneously according to their specified distributions, so that each replica is a full correlated prediction vector. This preserves parameter variations that alter the q² dependence, including those associated with the parameter λB. The empirical mean, standard deviations and correlation matrix of the replica ensemble then define the LCSR input to a correlated Gaussian χ². Because the predictions at different q² values are extremely strongly correlated, the raw correlation matrix is nearly singular; the authors regularize it mildly while retaining all correlation directions, and they verified that the regularization choice does not materially control the extracted |Vcb| or its uncertainty.

The numerical analysis proceeds through three nested fit scenarios with progressively enlarged constraint sets. The first relies purely on lattice-QCD form-factor simulations; the second adds the correlated SCET sum-rule predictions; the third supplements both with binned differential decay spectra and angular observables from Belle and Belle II, which is the only scenario in which |Vcb| can be determined. Comparing across the scenarios shows that incorporating the correlated large-recoil LCSR predictions dramatically shrinks the uncertainties of high-order BGL coefficients, while the inclusion of precise Belle II decay observables yields further moderate tightening of the parameter space. The complete global fit yields χ²min = 218.505 for 229 degrees of freedom, corresponding to a healthy p-value of 0.680, with the lattice-QCD, LCSR and experimental contributions to the χ² being 42.91, 14.44 and 161.16 respectively.

From the full global fit the team extracts |Vcb| = (39.18 ± 0.47) × 10⁻³. Individual channel fits give |Vcb| = (40.3 ± 0.8) × 10⁻³ from B → Dℓν and |Vcb| = (38.6 ± 0.6) × 10⁻³ from B → D*ℓν, and the authors document in detail how methodological choices — the exclusion of the Belle 2017 dataset, the retention of only the muonic mode of the Belle 2018 measurement to mitigate well-documented electron–muon tension, and the use of constrained frequentist χ² minimization rather than Bayesian MCMC or nested sampling — shift the central value relative to previous analyses. Robustness checks are thorough: a constrained profile-likelihood scan gives |Vcb| = (39.1801 +0.4611/−0.4577) × 10⁻³, consistent with the Hessian uncertainty, and increasing the BGL truncation from O(z²) to O(z³) shifts the central value by only 0.032 × 10⁻³ while the four strong unitarity bounds become essentially saturated, showing that the extra cubic coefficients are only weakly constrained by the present data.

Using the fitted form-factor coefficients from the combined lattice-QCD and LCSR fit, the team derives updated Standard Model predictions for the lepton-flavour-universality ratios R(D) = 𝓑𝓡(B → Dτντ)/𝓑𝓡(B → Dμνμ) = 0.3069 ± 0.0080 and R(D*) = 0.2548 ± 0.0043, alongside full differential decay distributions in the momentum transfer and angular variables for both the muon and tau channels. The muonic spectrum rises sharply near q² → 0 and decreases monotonically, while the tau channel has a nonvanishing kinematic threshold at finite q² owing to the tauon’s larger mass. Belle II measurements broadly follow the theoretical trend for the muonic channel, though mild local deviations appear in certain angular regions. At present precision, the authors emphasize, these discrepancies do not provide compelling evidence for physics beyond the Standard Model.

The persistent tension, however, remains vivid. Experimental world averages for both R(D) and R(D) lie systematically above the aggregate Standard Model theoretical band, and the team’s two-dimensional correlated 68% confidence contour in the (R(D), R(D)) plane shows little overlap with the experimental average. Their R(D*) central value aligns well with the global theory average, while their R(D) prediction sits slightly above pure lattice averages such as the FLAG 2024 result. The discrepancy may originate from underestimated lattice or LCSR systematic uncertainties, residual shape biases in the BGL z-expansion, missing higher-order or higher-power corrections, or unvalidated experimental reconstruction efficiencies. Improved lattice-QCD calculations, refined sum-rule predictions, and high-statistics Belle II measurements will be essential for determining whether the remaining gaps are theoretical artefacts or genuine signals of new physics in charged-current b → cℓν transitions — and this complete set of correlated Standard Model predictions now provides a robust, up-to-date reference for that decisive comparison.

Subject of Research: Exclusive determination of the CKM matrix element |Vcb| from semileptonic B meson decays

Article Title: Exclusive determination of (|V_{cb}|) from semileptonic decays (B\rightarrow D^{(*)}\ell \nu _{\ell })

Article References: Chen, X.-W., Cui, B.-Y., Zhang, J.-W., Zhao, X.-C., & Chen, Y.-H. (2026). Exclusive determination of $$|V_{cb}|$$ from semileptonic decays $$B\rightarrow D^{(*)}\ell \nu _{\ell }$$. The European Physical Journal C, 86(10), Article 1174. https://doi.org/10.1140/epjc/s10052-026-16486-2

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16486-2

Keywords: CKM matrix, Vcb, semileptonic decays, B mesons, lattice QCD, BGL parameterization, Belle II, R(D), R(D*), lepton flavour universality, Standard Model, flavour physics