Extraction of $|V_{cd}|$ and $|V_{cs}|$ from experimental decay rates using lattice QCD $D \to \pi(K) \ell \nu$ form factors

We present a determination of the Cabibbo-Kobayashi-Maskawa matrix elements $|V_{cd}|$ and $|V_{cs}|$ obtained by combining the momentum dependence of the semileptonic vector form factors $f_+^{D \to \pi}(q^2)$ and $f_+^{D \to K}(q^2)$, recently determined from lattice QCD simulations, with the differential rates measured for the semileptonic $D \to \pi \ell \nu$ and $D \to K \ell \nu$ decays. Our analysis is based on the results for the semileptonic form factors produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ flavors of dynamical quarks in the whole range of values of the squared 4-momentum transfer accessible in the experiments. The statistical and systematic correlations between the lattice data as well as those present in the experimental data are properly taken into account. With respect to the standard procedure based on the use of only the vector form factor at zero 4-momentum transfer, we obtain more precise and consistent results: $|V_{cd} |= 0.2341 ~ (74)$ and $|V_{cs} |= 0.970 ~ (33)$. The second-row CKM unitarity is fulfilled within the current uncertainties: $|V_{cd}|^2 + |V_{cs}|^2 + |V_{cb}|^2 = 0.996 ~ (64)$. Moreover, using for the first time hadronic inputs determined from first principles, we have calculated the ratio of the semileptonic $D \to \pi(K)$ decay rates into muons and electrons, which represent a test of lepton universality within the SM, obtaining in the isospin-symmetric limit of QCD: ${\cal{R}}_{LU}^{D\pi} = 0.985~(2)$ and ${\cal{R}}_{LU}^{DK} = 0.975~(1)$.


G. Salerno
Dipartimento di Matematica e Fisica and INFN, Università Roma Tre, Via della Vasca Navale 84, I-00146 Rome, Italy We present a determination of the Cabibbo-Kobayashi-Maskawa matrix elements |V cd | and |Vcs| obtained by combining the momentum dependence of the semileptonic vector form factors f D→π + (q 2 ) and f D→K + (q 2 ), recently determined from lattice QCD simulations, with the differential rates measured for the semileptonic D → π ν and D → K ν decays. Our analysis is based on the results for the semileptonic form factors produced by the European Twisted Mass Collaboration with N f = 2+1+1 flavors of dynamical quarks in the whole range of values of the squared 4-momentum transfer accessible in the experiments. The statistical and systematic correlations between the lattice data as well as those present in the experimental data are properly taken into account. With respect to the standard procedure based on the use of only the vector form factor at zero 4-momentum transfer, we obtain more precise and consistent results: |V cd | = 0.2341 (74) and |Vcs| = 0.970 (33). The secondrow CKM unitarity is fulfilled within the current uncertainties: |V cd | 2 + |Vcs| 2 + |V cb | 2 = 0.996 (64). Moreover, using for the first time hadronic inputs determined from first principles, we have calculated the ratio of the semileptonic D → π(K) decay rates into muons and electrons, which represent a test of lepton universality within the SM, obtaining in the isospin-symmetric limit of QCD: R Dπ LU = 0.985 (2) and R DK LU = 0.975 (1).

I. INTRODUCTION
Weak decays of hadrons represent an important source of direct information about the Cabibbo-Kobayashi-Maskawa (CKM) [1,2] matrix elements, which are fundamental parameters describing the quark flavor mixing in the electroweak sector of the Standard Model (SM). The CKM matrix is unitary in the SM and this gives rise to unitarity conditions between the elements of its rows and columns, whose investigation has been the focus of much of the experimental and theoretical efforts in Flavor Physics during the recent years and may provide stringent consistency tests of the SM.
The golden modes for the determination of the CKM entries |V cd | and |V cs | are represented by leptonic and semileptonic decays of D and D s mesons, which are sensitive probes of the c → d and c → s quark transitions. According to the V − A structure of the SM weak currents, leptonic and semileptonic D and D s decays should provide consistent results for the CKM elements |V cd | and |V cs |.
As is well known, the extraction of the CKM matrix elements from the experimental decay rates requires theoretical inputs, namely the leptonic decay constants and the semileptonic form factors, which encode the non-perturbative QCD dynamics and can be determined from first principles by simulating QCD on a lattice. The current Lattice QCD (LQCD) determinations of the leptonic decay constants f D and f Ds and of the semileptonic vector form factors at zero four-momentum transfer f D→π + (0) and f D→K + (0) are summarized in the latest FLAG review [3]. Uncertainties equal to 2.4% and 1.7% are quoted respectively for |V cd | and |V cs | extracted from the leptonic decays, while larger uncertainties equal to 4.5% and 2.6%, presently dominated by the theoretical errors, affect the determinations of |V cd | and |V cs | from semileptonic decays.
Within the SM the semileptonic D → P ν differential decay rate is given by is the daughter quark, | p P | is the momentum of the daughter P = π(K) pseudoscalar meson in the D-meson rest frame and q = (p D − p P ) is the 4-momentum of the outgoing lepton pair. Since in Eq. (1.1) the contribution of the scalar form factor f 0 (q 2 ) is proportional to m 2 , in the case of = e(µ) final leptons the differential decay rate simplifies to The analyses of the experimental data from BELLE [4], BABAR [5,6], CLEO [7] and BESIII [8,9] on D → π(K)e(µ)ν decays are based on a variety of parameterizations for the shape of the vector form factor f Dπ(K) + (q 2 ), like effective pole models, inspired by either dispersion relations [10] or heavy-quark expansion arguments [11], the z-expansion method [12] or even relativistic quark model predictions [13]. Thus, each experiment provides results for the products |V cx |f DP + (q 2 ) in various q 2 -bins. The LQCD calculations of the semileptonic D → π(K) form factors, fulfilling the quality criteria of the latest FLAG summary [3], provide the value of the vector form factor at zero 4-momentum transfer, f Dπ(K) + (q 2 = 0). This can be used to extract the relevant CKM entry adopting the experimental (averaged) values for the products |V cd |f D→π + (0) and |V cs |f D→K + (0), like the ones determined by the Heavy Flavor Averaging Group (HFAG) in Ref. [14] by combining all the experiments.
Recently [15] the momentum dependence of both the vector and scalar form factors f Dπ(K) + (q 2 ) and f Dπ(K) 0 (q 2 ) at the physical pion mass and in the continuum and infinite volume limits has been calculated on the lattice by the European Twisted Mass Collaboration (ETMC), obtaining an overall agreement with the momentum dependence of the experimental data from BELLE [4], BABAR [5,6], CLEO [7] and BESIII [8]. Some deviations have been nevertheless observed at high values of q 2 . At variance with the LQCD calculations considered in the FLAG review [3] the ETMC calculations of Ref. [15] employ gauge configurations with N f = 2 + 1 + 1 dynamical quarks (which include in the sea, besides two light mass-degenerate quarks, also the strange and charm quarks with masses close to their physical values), and provide also a set of synthetic (correlated) data points for the vector and scalar form factors in the whole range of values of q 2 accessible in the experiments, i.e. from q 2 = 0 up to q 2 max = (M D − M π(K) ) 2 . It is the aim of this work to combine the new LQCD results of Ref. [15] for the vector form factor f with the experimental data on the differential rates of the semileptonic D → π(K)e(µ)ν decays in order to perform a consistent analysis within the SM for the extraction of the CKM matrix elements |V cd | and |V cs | 1 . We show that a more precise and consistent determination of |V cd | and |V cs |, compared to the one based on the use of only the theoretical form factor values at q 2 = 0, can be obtained.

II. EXTRACTION OF |V cd | AND |Vcs|
The starting point is the partial decay rate provided by each experiment for various bins of values of q 2 (i.e., q 2 i ± ∆q 2 i /2 for i = 1, ..., N bins ). By integrating Eq. (1.2) in each experimental bin one has where the r.h.s. contains the phase-space integral over the vector form factor, viz.
Using the results of Ref. [15] for the vector form factor f DP + (q 2 ), we can combine the theoretical prediction I(q 2 i ) LAT of Eq. (2.2) with the experimental measurement of the partial decay rate (2.1) and get a determination of |V cx | for each experimental q 2 -bin: The values of ∆Γ(q 2 i ) EXP measured for the D → π and D → K semileptonic decays by BABAR [5,6], CLEO [7] and BESIII [8] collaborations are collected in Tables I-II, III-IV and V-VI, respectively. In the case of the BELLE experiment [4] the only available data are the values of |V cx |f DP + (q 2 i ), which result from a variety of model-dependent shapes adopted for describing the vector form factor. The values of |V cx |f DP + (q 2 i ) from the BELLE experiment are listed in Tables VII and VIII for D → π and D → K semileptonic decays, respectively. Experiment q 2 (GeV 2 ) ∆Γ(q 2 ) EXP · 10 16 (GeV) for the D → π transition in the q 2 -bins measured by BABAR [5].
Though the Belle data on |V cx |f DP + (q 2 i ) may contain some model dependence, they are considered for the extraction of |V cx | adopting the lattice determination of f DP + (q 2 ) from Ref. [15] at the center of each q 2 -bin 2 . Combining the BABAR, CLEO and BESIII data from Tables I-VI with the theoretical results for I(q 2 i ) LAT , and the BELLE data from Tables VII-VIII with the theoretical values for f DP + (q 2 i ), we get a determination of the CKM matrix element |V cx | for each experimental bin. The results are shown in Fig. 1 for both |V cd | and |V cs |. They exhibit an approximate constant behavior, except for |V cs | in the high-q 2 region, where some deviations are visible.
The data are strongly correlated, since among the various q 2 -bins the theoretical values of I(q 2 i ) LAT as well as the values of ∆Γ(q 2 i ) EXP coming from the same experiment are correlated. In order to take into account the correlations of the experimental data we have calculated the individual statistical+systematic covariance matrices for each experiments. Then, a global covariance matrix is obtained by combining in a block-diagonal form the separate covariance matrices given by BABAR, CLEO and BESIII collaborations in Refs. [5][6][7][8][9]. No covariance matrix is provided for the BELLE data on |V cx |f + (q 2 i ) in Ref. [4], which are treated as uncorrelated. According to each individual experimental covariance matrix we generate bootstrap events for the values of the experimental partial Experiment q 2 (GeV 2 ) ∆Γ(q 2 ) EXP · 10 16 (GeV) BESIII D 0 (0.00, 0.20) 6.14 ± 0.27 (0. 20  LAT , based on the vector form factor f DP + (q 2 ) evaluated on the lattice in Ref. [15], with the experimental data of the D → π and D → K semileptonic differential decay rates ∆Γ(q 2 i ) EXP , measured by BELLE [4], BABAR [5,6], CLEO [7] and BESIII [8,9] collaborations. The bands are described in the text.
decay rates ∆Γ(q 2 i ) EXP . Then, the bootstrap samplings of Ref. [15] are used to evaluate the lattice values of In this way the bootstrap distribution for V cx (q 2 i ) 2 is obtained and used to evaluate the final covariance matrix (always block-diagonal with respect to different experiments).  [4], collected in Ref. [19]. where the errors include all the various sources of statistical and systematic uncertainties of both the lattice form factors and of the experimental data. We have estimated that ≈ 80% of the errors quoted in Eq. (2.4) comes from the theoretical uncertainties and correlations. In order to check the stability of these results we have also performed a series of constant fits including only the data below each given value of q 2 , i.e., the number of data included in the fitting procedure increases as q 2 increases starting from q 2 0. The orange bands in Fig. 1 illustrate the results of these fits for |V cd | and |V cs | as a function of q 2 . The widths of the bands represent the spread of |V cd | and |V cs | at the level of one standard deviation. It can clearly be seen that: i) the variations of |V cd | and |V cs | are always within the uncertainties, and ii) the uncertainties themselves do not change appreciably when the data at the highest values of q 2 are included.
We have also carried out a series of constant fits including only the data above each given value of q 2 , i.e., the number of data included in the fitting procedure increases as q 2 decreases starting from q 2 q 2 max . The cyan bands in Fig. 1 illustrate the corresponding results. The extracted values of |V cx | differ from those corresponding to the orange bands only when the fits are restricted to the data at high values of q 2 . The tension appears to be at the level of one standard deviation. Once all the data are included in the fitting procedure (which for the cyan bands occurs close to q 2 0 and for the orange bands close to q 2 q 2 max ), the impact of the data at the highest values of q 2 is, however, negligible on the central values and errors given in Eq. (2.4).
Our results (2.4) can be compared with those obtained in Ref. [15] using the values of the vector form factor at q 2 = 0 and the experimental results for |V cd |f D→π + (0) and |V cs |f D→K + (0) provided by HFAG [14], namely |V cd | = 0.2330 (137) , |V cs | = 0.945 (38). (2.5) It turns out that the uncertainty of |V cd | in Eq. (2.4) is smaller by ≈ 50% with respect to the corresponding uncertainty in Eq. (2.5), while for |V cs | the reduction of the uncertainty is marginal. This is partially due to the higher degree of the correlations, found in Ref. [15], among the theoretical values of the vector form factor f DP + (q 2 ) in the various q 2 -bins in the case of the D → K transition with respect to the D → π one. We stress that the same theoretical input from LQCD is used for describing the shape of the vector form factor f Dπ(K) + (q 2 ) in all the experimental data, obtaining in this way a consistent SM analysis. The impact of the above consistency might become more significant as the precision of LQCD calculations of the semileptonic form factors will be improved in the future.
Thus, the theoretical information on f DP + (q 2 ) in the full q 2 -range has allowed not only to guarantee a consistent extraction of |V cd | and |V cs | within the SM, but also to get a more precise determination of |V cd |.
Our semileptonic results (2.4) can be compared with the determinations |V cd | = 0.2221(68) and |V cs | = 1.014 (25), obtained from the experimental D and D s leptonic decay rates [20] adopting the ETMC results [21] for the decay constants f D and f Ds . In Fig. 2 the above results from leptonic and semileptonic decays are reported as ellipses in the (|V cd |, |V cs |) plane corresponding to a 68% probability contour. The ellipses corresponding also to the leptonic and semileptonic FLAG averages [3] for |V cd | and |V cs | are shown as well as the constraint imposed by the second-row unitarity is indicated by a dotted line. Note that in the case of both the leptonic and semileptonic FLAG results the FIG. 2. Results for |V cd | and |Vcs| obtained from leptonic D-and Ds-meson and semileptonic D-meson decays, represented respectively by red and green ellipses corresponding to a 68% probability contour. The solid ellipses are the results of Ref. [21] and of this work, obtained with N f = 2+1+1 dynamical quarks. The correlations between |V cd | and |Vcs| are properly taken into account (see text). The striped ellipses correspond to the latest FLAG results [3], which for the semileptonic decays are based on the LQCD results obtained in Refs. [22,23] with N f = 2 + 1 dynamical quarks. The dashed line indicates the correlation between |V cd | and |Vcs| that follows from the CKM unitarity. correlation between |V cd | and |V cs | is not available, while in the case of the ETMC results the correlation between |V cd | and |V cs | have been properly considered, namely ρ(|V cd |, |V cs |) 0.6 in the case of leptonic decays and −0.1 for semileptonic decays.
Using |V cb | = 0.0360(9) from Ref. [24] and our semileptonic results (2.4) we can test the unitarity of the second row of the CKM matrix, obtaining |V cd | 2 + |V cs | 2 + |V cb | 2 = 0.996 (64), (2.6) which can be compared with the corresponding result |V cd | 2 + |V cs | 2 + |V cb | 2 = 0.949(78) from Ref. [15]. Before closing this Section, we mention that very recently [25] the BESIII collaboration has performed a test of the lepton universality (LU) in the D → π ν decays, obtaining   [15] for both the D → πµ(e)ν and D → Kµ(e)ν decays, we can predict for the first time both R Dπ LU and R DK LU within the SM using hadronic inputs calculated from first principles. Thanks to the strong correlation between the numerator and the denominator we get quite precise values in the isospin-symmetric limit of QCD, namely R DK LU = 0.975 (1) . (2.10) The measured values (2.7)-(2.8) are consistent with our SM prediction (2.9) within 2.3σ and 1.2σ, respectively.

III. CONCLUSIONS
We have presented a determination of the CKM matrix elements |V cd | and |V cs | obtained by combining the momentum dependence of the semileptonic vector form factors f D→π + (q 2 ) and f D→K + (q 2 ), recently determined from lattice QCD simulations, with the differential rates measured for the semileptonic D → π ν and D → K ν decays.
Our analysis is based on the results for the semileptonic form factors produced by ETMC with N f = 2 + 1 + 1 flavors of dynamical quarks in the whole range of values of the squared 4-momentum transfer accessible in the experiments [15]. The statistical and systematic correlations between the lattice data as well as those present in the experimental data are properly taken into account.
With respect to the standard procedure based on the use of only the vector form factor at zero 4-momentum transfer, we obtain more precise and consistent results: |V cd | = 0.2341 (74) and |V cs | = 0.970 (33). The second-row CKM unitarity is fulfilled within the current uncertainties: |V cd | 2 + |V cs | 2 + |V cb | 2 = 0.996 (64).
The results presented in this work depend crucially on the uncertainties and correlations of the input quantities, namely the experimental data available for the differential decay rates and the theoretical calculations of the semileptonic form factors of Ref. [15]. Future improvements of the precision of the above quantities are mandatory in order to assess the significance of any possible discrepancies with the Standard Model predictions, like in the case of the slight tension observed at high values of q 2 (see Fig. 1).
Moreover, using for the first time hadronic inputs determined from first principles, we have calculated the ratio of the semileptonic D → π(K) decay rates into muons and electrons, which represent a test of lepton universality within the SM, obtaining in the isospin-symmetric limit of QCD: R Dπ LU = 0.985 (2) and R DK LU = 0.975 (1).