# World Average Value for the Vcd from leptonic D+ decays (updated 10 December 2014)

People working on this:   Anze Zupanc

The branching fraction of leptonic decays of pseudo-scalar mesons that proceed via the annihilation of the initial quark-antiquark pair (cd in the case of D+ meson) into a virtual W+ that finally materializes as an antilepton-neutrino pair (lν, where l = e, μ, or τ) is given in the Standard Model by

Br(D+→ l+ν) = GF2/(8π)τDfD2|Vcd|2MDMl2(1-Ml2/MD2)2.

Here, MD is the D meson mass, τD is its lifetime, Ml is the charged lepton mass, |Vcd| is the magnitude of the CKM matrix element, and GF is the Fermi coupling constant. The parameter fD is the D meson decay constant and is related to the wave-function overlap of the meson's constituent quark and anti-quark. Within the SM, the decay constants have been predicted using several methods, the most precise being the lattice gauge theory (LQCD) calculations. The
Flavor Lattice Averaging Group combines all LQCD calculations and provides averaged value for fDLQCD=(209.2 ± 3.3) MeV that is used within this section to extract the magnitudes of the |Vcd| CKM matrix element from experimentally measured branching fractions of leptonic D+ meson decays.

Experimental results and world averages for branching fractions of purely leptonic D+ decays and fD|Vcd|.
The first uncertainty is statistical and the second is experimental systematic. The third uncertainty in the case fD|Vcd| is due to external inputs (dominated by the uncertainty of D meson's lifetime, τD).
Mode Branching Fraction fD|Vcd| (MeV) Experiment Source
μ+ ν (3.82 ± 0.32 ± 0.09) × 10-4 46.4 ± 1.9 ± 0.5 ± 0.2 CLEO B. I. Eisenstein et al. (CLEO Collaboration), Phys. Rev. D 78, 052003 (2008)
(3.71 ± 0.19 ± 0.06) × 10-4 45.7 ± 1.2 ± 0.4 ± 0.2 BESIII M. Ablikim et al. (BESIII Collaboration), Phys. Rev. D 89, 051104 (2014)
(3.74 ± 0.16 ± 0.05) × 10-4 45.9 ± 1.0 ± 0.3 ± 0.2 HFAG Average
e+ ν < 8.8 × 10-6 at 90% C.L. CLEO B. I. Eisenstein et al. (CLEO Collaboration), Phys. Rev. D 78, 052003 (2008)
τ+ ν < 1.2 × 10-3 at 90% C.L.

We use measurements of the branching fraction Br(D+→ μ+ν) from CLEO and BESIII to calculate its world average (WA) value. Uncertainties of the two measurements are considered to be uncorrelated. We obtain

Br(D+→ μ+ν)WA = (3.74 ± 0.17) × 10-4,

from which we determine the product of the decay constant and the CKM matrix element using the equation from the top of the page to be (figure)

fD|Vcd| = (45.9 ± 1.1) MeV,

where the uncertainty includes the uncertainty on Br(D+→ μ+ν)WA and external inputs. Using the LQCD value for the decay constant, fDLQCD=(209.2 ± 3.3) MeV, we finally obtain the CKM matrix element Vcd to be

|Vcd|D→ lν = 0.219 ± 0.005(exp.) ± 0.003(LQCD),

where the uncertainties are from the experiments and lattice calculations, respectively. The value is found to be consistent with the one determined from the semileptonc D→ πlν decays, |Vcd|D→ πlν = 0.214 ± 0.003(exp.) ± 0.009(LQCD). Averaging both, assuming 100% correlation between the LQCD uncertainties, and none between the experimental uncertainties, gives (figure)

|Vcd|D→ (π)lν = 0.219 ± 0.006.

Assuming unitarity of the CKM matrix, the value of the element relevant in the case of leptonic D+ decays is known from the global fit of the CKM matrix, |Vcd| = 0.22537+0.00068-0.00035. This value can be used to extract the D meson decay constant from the experimentally measured product fD|Vcd|. This leads to the experimentally measured D meson decay constant to be:

fD = 203.7 ± 4.9 MeV,

which is in agreement with the LQCD determination given above.

WA value for fD|Vcd|. For each point, the first error listed is the statistical
and the second error is the systematic error. (Click on figure for higher resolution.)

Comparison of magnitudes of the CKM matrix element |Vcd| determined from the leptonic and semileptonic D meson decays and from neutrino scattering data and indirect determination from the global fit assuming CKM unitarity. (Click on figure for higher resolution.)

External parameter values as taken from the 2014 Particle Data Book:
GF /(ℏc)3 = (1.1663787 ± 0.0000006) × 10-5 GeV-2
Mμ = (0.1056583715 ± 0.0000000035) GeV/c2
MD = (1.86962 ± 0.00015) GeV/c2
τD = (1040 ± 7) × 10-15 s