Model-Independent Extraction of $|V_{cb}|$ from $\bar{B}\rightarrow D^* \ell \overline{\nu}$

We fit the unfolded data of $\bar{B}^0\rightarrow D^{*+} \ell \overline{\nu}$ from the Belle experiment, where $\ell \equiv e, \mu$, using a method independent of heavy quark symmetry to extrapolate to zero-recoil and extract the value of $|V_{cb}|$. This results in $|V_{cb}| = (41.9^{~+2.0}_{~-1.9})\times 10^{-3}$, which is robust to changes in the theoretical inputs and very consistent with the value extracted from inclusive semileptonic $B$ decays.


I. INTRODUCTION
The discrepancy between the measured values of the CKM matrix element |V cb | from inclusive versus exclusive semileptonic B decays has been an ongoing dilemma for a few decades; for a recent review, see Ref. [1].Currently, the world averages for |V cb | are [2,3]:  [4,5], which utilizes not only dispersion relations to bound the hadronic form factors, but also relations at 1/m Q in the heavy quark expansion, and claims to describe the full differential decay to within a few percent given only 4 parameters.Interestingly, it has now become clear that the CLN parameterization is no longer a good fit to the B → D ν data from current experiments [6].One would naively expect the same to be true, in general, for B → D * ν as well.
A recent analysis by Belle [7] used the CLN parameterization to extract a value for |V cb | at zero recoil given a dataset of B0 → D * + ν, resulting in |V cb | = (37.4± 1.3) × 10 −3 .This is one of the first times an experiment with a large dataset has published unfolded data for B → D * ν, so one can easily explore if a different parameterization of the hadronic form factors gives rise to a different value of |V cb |.We use the method developed by Boyd, Grinstein, and Lebed (BGL) in Refs.[8][9][10][11][12], from which the original CLN parameterization was based, but which does not utilize any assumptions from heavy quark symmetry. 1We find that while both the CLN and BGL parameterizations give a good fit to the Belle data, the BGL parameterization gives a larger value for |V cb | when extrapolating to zero recoil, consistent with the value of |V cb | measured in inclusive analyses [2,3].As we prepared this manuscript, a similar analysis was made public, using the same BGL parameterization [14].We are pleased to see the results match nearly identically.Another analysis was recently made public, which highlighted the importance of including theoretical uncertainties from Λ QCD /m Q corrections in the relations between the form factors in the heavy quark expansion [15], which suggests that the uncertainties associated with the extracted value of |V cb | using the CLN method may be underestimated.Moreover, in their analysis of B → D ν lattice data, the Fermilab Lattice and MILC collaborations decide not to quote results of their CLN fits because they are more confident in the errors obtained from the BGL parameterization which, they state, can be used to obtain |V cb | even as the uncertainties become arbitrarily more precise [16].
In all, the current evidence points to the strong possibility that the tension between inclusive and exclusive measurements of |V cb | in semileptonic B decays may be due to the use of the CLN parameterization of B → D ( * ) ν, and a different parameterization should be employed.

II. DIFFERENTIAL DECAY RATE AND BGL PARAMETERIZATION
Using the notation in Ref. [12], the B → D * matrix elements are defined as where ε µ is the polarization tensor of the vector D * meson.In the limit when the final-state leptons are massless, the full differential decay rate for B → D * ν is where q µ is the 4-momentum of the lepton system, r ≡ m D * /m B , and |p D * | is the magnitude of the D * 3-momentum in the rest frame of the B: Here, H + , H − , and H 0 are form factors associated with each of the three helicity states of the D * , all of which are functions of q 2 .Also, θ is the angle between the anti-neutrino and the direction antiparallel to the D * in the rest frame of the leptonic system, θ v is the angle between the D * momentum and its daughter D meson, and χ is the angle between the planes defined by the the leptonic system and the D * system.The factor η ew incorporates the leading electroweak corrections [17], In terms of the form factors in Eqs. ( 5) and ( 6), A detailed discussion about the BGL method for parameterizing the form factors f , g, and F 1 can be found in Ref. [12].The final result gives a parametrization of each form factor in terms of N + 1 coefficients: where the conformal variable z is defined as Here, a = 1 can be chosen such that z = 0 corresponds to zero recoil, and the coefficients a n , b n , and c n are bounded by unitarity [10], From Eq. ( 11), F 1 (0) = (m B − m D * )f (0); hence b 0 and c 0 are not independent, i.e., The Blaschke factors P (z) remove poles for q 2 < (m B + m D * ) 2 associated with on-shell production of B * c bound states: For the form factor g, the index i runs over the 4 vector B * c states, and for the f and F 1 , it runs over the 4 axial vector states, as listed in Table I.The outer functions φ are defined as: where n I is the effective number of light quarks, u ≡ m c /m b , and χ T (±u) is related to a perturbative calculation at q 2 = 0. We use the result for χ T (±u) in the pole mass scheme, and take numerical values to be those listed in Table II, as done in Ref. [12], ignoring the small contributions from condensates:   We do not incorporate any uncertainties associated with the numerical inputs to the BGL theory parameterization.Specifically, the uncertainties associated with m B and m D * are negligible in this analysis.Also, varying the inputs that only affect the overall normalization of the outer function, namely n I and those used in the perturbative calculation of χ T , amounts to changing the right-hand side of the bounds in Eq. ( 14) by a corresponding fractional amount.We impose the bounds Eq. ( 14) as constraints, however, the resultant best-fit values are orders of magnitude from saturating the unitarity bounds, and these constraints effectively play no role.Because of this, changing the values of, say, α s , m b , etc., will change the best-fit values for the free parameters in the theory, but will not change the extracted value |V cb |, to a very good approximation.The remaining inputs are B * c pole masses.They affect both the effectiveness of the unitarity bound and the shape of form factors.For example, for F ≡ g, f, F 1 , the unitarity bound implies that the truncation error for finite N in Eq. ( 12) is bounded over the physical region Since P F (z) varies by less than 15% when the input pole masses are individually varied by 1%, the ability to control the expansion is unaffected by uncertainties in pole masses.Also, the precise location of the pole masses affects the shape of the form factors in the physical region, but the parameters of the fit accomodate small changes in pole masses easily.For example, consider the N = 0 parametrization of the form factor f .Varying the value of the lowest pole mass by 1%, the change in f is compensated by a multiplicative factor of 1.07 − 0.17z to better than 1% over the whole physical region.For these reasons, we expect that while the values of the fit parameters may depend somewhat sensitively on the inputs, the extracted value of |V cb | does not.We illustrate this by using the same input values as Ref. [12], as listed in Tables I and II, which are about 20 years old, but find that our extracted value of |V cb | is nearly identical to that of a recent analysis by the authors of Ref. [14], which uses the current state-of-the-art theoretical inputs.Extracting the value of |V cb | using the BGL parameterization is robust to changes in the numerical inputs.We discuss this more in Section III.

III. FITTING RESULTS AND EXTRACTING |V cb |
We fit to the unfolded data provided by Belle in Ref. [7], using the numerical inputs found in Tables I and II.We choose to truncate the series in Eq. ( 12) at O(z 2 ) for f and g, and at order O(z 3 ) for F 1 , as well as define the value c 0 in terms of b 0 , as in Eq. ( 15).This results in 6 free parameters, defined as ãi ≡ |V cb |η ew a i , bi ≡ |V cb |η ew b i , where i = 0, 1, and ci ≡ |V cb |η ew c i , where i = 1, 2. We choose to truncate to these orders because it results in a good fit to the data.A Markov Chain Monte Carlo (MCMC) is used to estimate the shape of the likelihood and the shape around its extremum.The result of the fit is shown in Fig. 1.We find χ 2 min /dof 27.7/34.Comparing to the Belle data, our best fit value for the BGL parameterization can be found in Fig. 2, compared to the fit performed by Belle using the CLN parameterization.
One requires input from the lattice to extract the value of |V cb |.The B → D * ν rate is given in terms of a function F(w) defined as [20]: One can estimate on the lattice the value of F(1).Given the definition of the form factors in Eqs. ( 7), (9), and (12), the relationship between |V cb |η ew F(1) and our fitting parameters is: Given the statistical ensemble of b0 provided by the MCMC, we can estimate the likelihood, and thus the variance in the χ 2 about its minimum, i.e., ∆χ 2 , for the value of |V cb |η ew F(1) extracted from the Belle data, as shown in Fig. 3.This results in: Using the η ew = 1.0066 and the FNAL/MILC value of F(1) = 0.906 ± 0.013 from Ref. [20], the final result for |V cb | is This result is nearly identical to that of the recent analysis by the authors of Ref. [14], and systematically larger than the value extracted by the Belle experiment using the CLN parameterization [7].Our result is comparable to a recent analysis by the authors of Ref. [15], which reevaluate the CLN parameterization, using relations from heavy quark symmetry, but not ignoring important sources of theoretical uncertainty when fitting to the same data.As explained above, extracting the value of |V cb | using the BGL parameterization is robust to changes in the The result of the fit to the form factors defined in Eq. ( 12), where ãi, bi, and ci are defined as ηew|V cb |ai, ηew|V cb |bi, and ηew|V cb |ci, respectively.The blue, yellow, and red lines correspond to the 68.27%, 95%, and 99% CL, respectively.Any jaggedness in the contours is from the finite size of the MCMC ensemble.
numerical inputs.We explore further the effect of Blaschke factors on the extraction of |V cb | by squaring the Blaschke factors for f and F 1 .One might, a priori, expect to produce radical variations in the z dependence of the form factors. Surprisingly, while we fit very different values for the free parameters, the extracted value of |V cb |η ew F(1) is nearly unchanged (after updating the expression in Eq. ( 23) to account for the different definition of the form factors).This gives further evidence that the BGL parameterization is robust when extrapolating to zero recoil.On the other hand, the precise values of the theoretical inputs may be essential when extracting the functional dependence of the form factors over the entire kinematic range.FIG.2: The Belle data (black points) compared to the results of our fit using the BGL parameterization (blue), and the results of the Belle analysis using the CLN parameterization (red).

IV. SUMMARY AND CONCLUSIONS
The long-standing tension between the exclusive and inclusive determinations of |V cb | in semileptonic B decays could have been exacerbated by underestimated uncertainties in the exclusive analyses when using the CLN parameterization [6,14,15].As it stands, the data has become too precise to use the original CLN proposal in the case of B → D ν, as discussed in Ref. [6], and one might expect the same could be true for B → D * ν.
The Belle experiment recently released an unfolded dataset of B0 → D * + ν, which resulted in the following value of |V cb | using the CLN parameterization [7]: Using this data, we fit to the BGL parameterization, which does not rely on any assumptions regarding heavy quark symmetry, and obtain: |V cb | = (41.9+2.0 −1.9 ) × 10 This result is consistent with the value of |V cb | measured using inclusive semileptonic B decays.We are able to corroborate the results of the analysis in Ref. [14], which appeared as our manuscript was in preparation.One can see by eye in Fig. 2 that the BGL parametrization gives a different fit to the data, and a systematically larger value of the B → D * ν rate near zero-recoil, than when using the CLN parameterization.Furthermore, we find that the BGL method is robust when extrapolating to zero-recoil, i.e., large changes in the inputs associated with the form factors in Eq. ( 12) result in almost no change in the extracted value of |V cb |.

ã0 [ 10 − 2 ] 11 −
FIG.1:The result of the fit to the form factors defined in Eq.(12), where ãi, bi, and ci are defined as ηew|V cb |ai, ηew|V cb |bi, and ηew|V cb |ci, respectively.The blue, yellow, and red lines correspond to the 68.27%, 95%, and 99% CL, respectively.Any jaggedness in the contours is from the finite size of the MCMC ensemble.

TABLE I :
[18]c masses used in this analysis from[18].Only states that have a mass less than mB + mD * are included.

TABLE II :
[19]rical inputs for this analysis.The meson masses are taken from Ref.[19].See the text at the end of Section II for a discussion regarding the choices of the other parameters.