World Average Value for the Vcs from leptonic Ds+ 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 (cs in the case of Ds+ 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(Ds+→ l+ν) = GF2/(8π)τDsfDs2|Vcs|2MDsMl2(1-Ml2/MDs2)2.

Here, MDs is the Ds meson mass, τDs is its lifetime, Ml is the charged lepton mass, |Vcs| is the magnitude of the CKM matrix element, and GF is the Fermi coupling constant. The parameter fDs is the Ds 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 fDsLQCD=(248.6 ± 2.7) MeV that is used within this section to extract the magnitudes of the |Vcs| CKM matrix element from experimentally measured branching fractions of leptonic Ds+ meson decays.

 

Experimental results and world averages for branching fractions of purely leptonic Ds+ decays and fDs|Vcs|.
The first uncertainty is statistical and the second is experimental systematic. The third uncertainty in the case fDs|Vcs| is due to external inputs (dominated by the uncertainty of Ds meson's lifetime, τDs). We have recalculated Br(Ds+→ τν) quoted by CLEO-c and BaBar using the latest values for branching fractions of τ decays to electron, muon, or pion and neutrinos given in the 2014 PDG edition. CLEO-c and BaBar include statistical uncertainty of number of Ds tags (denominator in the calculation of branching fraction) in the statistical uncertainty of measured branching fraction. We subtract this uncertainty from the statistical one and add it to the systematic uncertainty.
Mode Branching Fraction fDs|Vcs| (MeV) Experiment Source
μ+ ν (5.65 ± 0.44 ± 0.20) × 10-3 250.8 ± 9.8 ± 4.4 ± 1.8 CLEO J. P. Alexander et al. (CLEO Collaboration), Phys. Rev. D 79, 052001 (2009)
(6.02 ± 0.37 ± 0.32) × 10-3 258.9 ± 8.0 ± 6.9 ± 1.8 BaBar P. del Amo Sanchez et al. (BaBar Collaboration), Phys. Rev. D 82, 091103 (2010)
(5.31 ± 0.28 ± 0.20) × 10-3 243.1 ± 6.4 ± 4.6 ± 1.7 Belle A. Zupanc et al. (Belle Collaboration), JHEP 1309, 139 (2013)
(5.57 ± 0.20 ± 0.14) × 10-3 249.0 ± 4.5 ± 3.1 ± 1.7 HFAG Average
τ+(e+) ν (5.31 ± 0.47 ± 0.22) × 10-2 246.1 ± 10.9 ± 5.1 ± 1.7 CLEO P. U. E. Onyisi et al. (CLEO Collaboration), Phys. Rev. D 79, 052002 (2009)
τ++) ν (6.46 ± 0.80 ± 0.23) × 10-2 271.4 ± 16.8 ± 4.8 ± 1.9 J. P. Alexander et al. (CLEO Collaboration), Phys. Rev. D 79, 052001 (2009)
τ++) ν (5.50 ± 0.54 ± 0.24) × 10-2 250.4 ± 12.3 ± 5.5 ± 1.8 P. Naik et al. (CLEO Collaboration), Phys. Rev. D 80, 112004 (2009)
τ+ ν (5.57 ± 0.32 ± 0.15) × 10-2 252.0 ± 7.2 ± 3.4 ± 1.8 CLEO
τ+(e+) ν (5.08 ± 0.52 ± 0.68) × 10-2 240.7 ± 12.3 ± 16.1 ± 1.7 BaBar P. del Amo Sanchez et al. (BaBar Collaboration), Phys. Rev. D 82, 091103 (2010)
τ++) ν (4.90 ± 0.46 ± 0.54) × 10-2 236.4 ± 11.1 ± 13.0 ± 1.7
τ+ ν (4.95 ± 0.36 ± 0.58) × 10-2 237.6 ± 8.6 ± 13.8 ± 1.7 BaBar
τ+(e+) ν (5.37 ± 0.33 +0.35-0.31) × 10-2 247.4 ± 7.6 +8.1-7.1 ± 1.7 Belle A. Zupanc et al. (Belle Collaboration), JHEP 1309, 139 (2013)
τ++) ν (5.86 ± 0.37 +0.34-0.59) × 10-2 258.5 ± 8.2 +7.5-13.0 ± 1.8
τ++) ν (6.04 ± 0.43 +0.46-0.40) × 10-2 262.4 ± 9.3 +10.0-8.7 ± 1.8
τ+ ν (5.70 ± 0.21 ± 0.31) × 10-2 254.9 ± 4.7 ± 6.9 ± 1.8 Belle
τ+ ν (5.55 ± 0.18 ± 0.17) × 10-2 251.5 ± 4.1 ± 3.9 ± 1.8 HFAG average
μ+ ν
τ+ ν
250.6 ± 3.1 ± 2.8 ± 1.8 HFAG average
e+ ν < 8.3 × 10-5 at 90% C.L. Belle A. Zupanc et al. (Belle Collaboration), JHEP 1309, 139 (2013)

 

We use measurements of the branching fraction Br(Ds+→ μ+ν) from CLEO, BaBar and Belle and obtain its world average (WA) value of

Br(Ds+→ μ+ν)WA = (5.57 ± 0.24) × 10-3.

The WA value for Br(Ds+→ τ+ν) is also calculated from CLEO, BaBar and Belle measurements. CLEO made separate measurements for leptonic τ decays and hadronic decay to ρ; BaBar made measurements using the leptonic τ decays only; and Belle made measurements using the leptonic τ decays and the hadronic decay to π. Combining all of them we obtain the WA value of

Br(Ds+→ τ+ν)WA = (5.55 ± 0.24) × 10-2.

The ratio of the branching fractions of leptonic Ds+ decays is thus found to be

R = Br(Ds+→ τ+ν)/Br(Ds+→ μ+ν) = 9.96 ± 0.57,

and is consistent with the value expected in the SM (from lepton universality), RSM = 9.76 ± 0.03.

 

From the average values of branching fractions of muonic and tauonic decays we determine from the equation on top of this page the product of Ds meson decay constant and the |Vcs| CKM matrix element to be (figure)

fDs |Vcs| = (250.6 ± 4.5) MeV,

where the uncertainty includes the uncertainties on Br(Ds+→ μ+ν)WA and Br(Ds+→ τ+ν)WA and external inputs. To obtain the average we have taken into account the correlations within each experiment for the uncertainties related to: normalization, tracking, particle identification, signal and background parameterizations, and peaking background contributions. Using the LQCD value for the decay constant, fDsLQCD=(248.6 ± 2.7) MeV, we finally obtain the CKM matrix element Vcs to be

|Vcs|Ds→ lν = 1.008 ± 0.018(exp.) ± 0.011(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→ Klν decays, |Vcs|D→ Klν = 0.975 ± 0.007(exp.) ± 0.025(LQCD). Averaging both, assuming 100% correlation between the LQCD uncertainties, and none between the experimental uncertainties, gives (figure)

|Vcs|D(s)→ (K)lν = 0.998 ± 0.020.

 

Assuming unitarity of the CKM matrix, the value of the element relevant in the case of leptonic Ds+ decays is known from the global fit of the CKM matrix, |Vcs| = 0.973395+0.000095-0.000176. This value can be used to extract the Ds meson decay constant from the experimentally measured product fDs|Vcs|. This leads to the experimentally measured Ds meson decay constant to be:

fDs = 257.4 ± 4.6 MeV,

which is in agreement with the LQCD determination given above. Using the D meson decay constant determined from the leptonic decays we find the ratio of the Ds and D meson decay constants to be

fDs/fD = 1.264 ± 0.038.

The ratio is in agreement with the average of LQCD calculations, fDs/fD = 1.187 ± 0.012, at the level of 1.9σ.

 


WA value for fDs|Vcs|. 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 |Vcs| determined from the leptonic and semileptonic D meson decays and from W decays 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
            Mτ = (1.77682 ± 0.00016) GeV/c2
            MDs = (1.96850 ± 0.00032) GeV/c2
            τDs = (500 ± 7) × 10-15 s

This page is maintained by A. Zupanc and was last updated