Leptonic decay constants f(K), f(D), and f(Ds) with N-f=2+1+1 twisted-mass lattice QCD

We present a lattice QCD calculation of the pseudoscalar decay constants f(K), f(D) and f(Ds) performed using the gauge configurations produced by the European Twisted Mass Collaboration 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 values in the real world. The simulations are based on a unitary setup for the two light mass-degenerate quarks and on a mixed action approach for the strange and charm quarks. We use data simulated at three different values of the lattice spacing in the range 0.06-0.09 fm and at pion masses in the range 210-450 MeV. Our main results are f(K+)/f(pi+) = 1.184(16), f(K+) = 154.4(2.0) MeV, which incorporate the leading strong isospin breaking correction due to the up and down quark mass difference, and f(K) = 155.0(1.9) MeV, f(D) = 207.4(3.8) MeV, f(Ds) = 247.2(4.1) MeV, f(Ds)/f(D) = 1.192(22) and (f(Ds)/f(D))/(f(K)/f(pi)) = 1.003(14) obtained in the isospin symmetric limit of QCD. Combined with the experimental measurements of the leptonic decay rates of kaon, pion, D and D-s mesons our results lead to the following determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements: vertical bar V-us vertical bar = 0.2269(29), vertical bar V-cd vertical bar = 0.2221(67) and vertical bar V-cs vertical bar = 1.014(24). Using the latest value of vertical bar V-ud vertical bar from superallowed nuclear beta decays the unitarity of the first row of the CKM matrix is fulfilled at the per mill level.