• 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)$.
  • We present the results of the first lattice QCD calculation of the $K \to \pi$ matrix elements of the chromomagnetic operator $O_{CM} = g\, \bar s\, \sigma_{\mu\nu} G_{\mu\nu} d$, which appears in the effective Hamiltonian describing $\Delta S = 1$ transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been determined non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing. Our result for the B-parameter of the chromomagnetic operator at the physical pion and kaon point is $B_{CMO}^{K \pi} = 0.273 ~ (70)$, while in the SU(3) chiral limit we obtain $B_{CMO} = 0.072 ~ (22)$. Our findings are significantly smaller than the model-dependent estimate $B_{CMO} \sim 1 - 4$, currently used in phenomenological analyses, and improve the uncertainty on this important phenomenological quantity.
  • We present the first lattice $N_f=2+1+1$ determination of the tensor form factor $f_T^{D \pi(K)}(q^2)$ corresponding to the semileptonic $D \to \pi(K) \ell \nu_\ell$ decays as a function of the squared four-momentum transfer $q^2$. Together with our recent determination of the vector and scalar form factors we complete the set of hadronic matrix elements regulating the semileptonic $D \to \pi(K)$ transition within and beyond the Standard Model, where a non-zero tensor coupling is possible. Our analysis is based on the gauge configurations produced by ETMC with $N_f=2+1+1$ flavors of dynamical quarks, which include three different values of the lattice spacing and pion masses as small as 220 MeV. The matrix elements of the tensor current are determined for a plethora of kinematical conditions in which parent and child mesons are either moving or at rest. As in the case of the vector and scalar form factors, Lorentz symmetry breaking due to hypercubic effects is clearly observed also in the data for the tensor form factor and included in the decomposition of the current matrix elements in terms of additional form factors. After the extrapolations to the physical pion mass and to the continuum and infinite volume limits we determine the tensor form factor in the whole kinematical region accessible in the experiments. A set of synthetic data points for $f_T^{D \pi(K)}(q^2)$ at several selected values of $q^2$ is provided and the corresponding covariance matrix is also available. At $q^2=0$ we get $f_T^{D \pi}(0) = 0.506(79)$ and $f_T^{D K}(0) = 0.687(54)$, which correspond to $f_T^{D \pi}(0)/f_+^{D \pi}(0) = 0.827(114)$ and $f_T^{D K}(0)/f_+^{D K}(0)= 0.898(50)$.
  • The leading-order electromagnetic and strong isospin-breaking corrections to the ratio of $K_{\mu 2}$ and $\pi_{\mu 2}$ decay rates are evaluated for the first time on the lattice, following a method recently proposed. The lattice results are obtained using the gauge ensembles produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks. Systematics effects are evaluated and the impact of the quenched QED approximation is estimated. Our result for the correction to the tree-level $K_{\mu 2} / \pi_{\mu 2}$ decay ratio is $-1.22\,(16) \%$ to be compared to the estimate $-1.12\,(21) \%$ based on Chiral Perturbation Theory and adopted by the Particle Data Group.
  • We present an investigation of the electromagnetic pion form factor, $F_\pi(Q^2)$, at small values of the four-momentum transfer $Q^2$ ($\lesssim 0.25$ GeV$^2$), based on the gauge configurations generated by European Twisted Mass Collaboration with $N_f = 2$ twisted-mass quarks at maximal twist including a clover term. Momentum is injected using non-periodic boundary conditions and the calculations are carried out at a fixed lattice spacing ($a \simeq 0.09$ fm) and with pion masses equal to its physical value, 240 MeV and 340 MeV. Our data are successfully analyzed using Chiral Perturbation Theory at next-to-leading order in the light-quark mass. For each pion mass two different lattice volumes are used to take care of finite size effects. Our final result for the squared charge radius is $\langle r^2 \rangle_\pi = 0.443~(29)$ fm$^2$, where the error includes several sources of systematic errors except the uncertainty related to discretization effects. The corresponding value of the SU(2) chiral low-energy constant $\overline{\ell}_6$ is equal to $\overline{\ell}_6 = 16.2 ~ (1.0)$.
  • We present a study of the isospin-breaking (IB) corrections to pseudoscalar (PS) meson masses using the gauge configurations produced by the ETM Collaboration with $N_f=2+1+1$ dynamical quarks at three lattice spacings varying from 0.089 to 0.062 fm. Our method is based on a combined expansion of the path integral in powers of the small parameters $(\widehat{m}_d - \widehat{m}_u)/\Lambda_{QCD}$ and $\alpha_{em}$, where $\widehat{m}_f$ is the renormalized quark mass and $\alpha_{em}$ the renormalized fine structure constant. We obtain results for the pion, kaon and $D$-meson mass splitting; for the Dashen's theorem violation parameters $\epsilon_\gamma(\overline{\mathrm{MS}}, 2~\mbox{GeV})$, $\epsilon_{\pi^0}$, $\epsilon_{K^0}(\overline{\mathrm{MS}}, 2~\mbox{GeV})$; for the light quark masses $(\widehat{m}_d - \widehat{m}_u)(\overline{\mathrm{MS}}, 2~\mbox{GeV})$, $(\widehat{m}_u / \widehat{m}_d)(\overline{\mathrm{MS}}, 2~\mbox{GeV})$; for the flavour symmetry breaking parameters $R(\overline{\mathrm{MS}}, 2~\mbox{GeV})$ and $Q(\overline{\mathrm{MS}}, 2~\mbox{GeV})$ and for the strong IB effects on the kaon decay constants.
  • We present a lattice calculation of the Hadronic Vacuum Polarization (HVP) contribution of the strange and charm quarks to the anomalous magnetic moment of the muon including leading-order electromagnetic (e.m.) corrections. We employ the gauge configurations generated by the European Twisted Mass Collaboration (ETMC) with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) with pion masses in the range $M_\pi \simeq 210 - 450$ MeV. The strange and charm quark masses are tuned at their physical values. Neglecting disconnected diagrams and after the extrapolations to the physical pion mass and to the continuum limit we obtain: $a_\mu^s(\alpha_{em}^2) = (53.1 \pm 2.5) \cdot 10^{-10}$, $a_\mu^s(\alpha_{em}^3) = (-0.018 \pm 0.011) \cdot 10^{-10}$ and $a_\mu^c(\alpha_{em}^2) = (14.75 \pm 0.56) \cdot 10^{-10}$, $a_\mu^c(\alpha_{em}^3) = (-0.030 \pm 0.013) \cdot 10^{-10}$ for the strange and charm contributions, respectively.
  • We present a lattice calculation of the Hadronic Vacuum Polarization (HVP) contribution of the strange and charm quarks to the anomalous magnetic moment of the muon including leading-order electromagnetic corrections. We employ the gauge configurations generated by the European Twisted Mass Collaboration (ETMC) with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) with pion masses in the range $M_\pi \simeq 210 - 450$ MeV. The strange and charm quark masses are tuned at their physical values. Neglecting disconnected diagrams and after the extrapolations to the physical pion mass and to the continuum limit we obtain: $a_\mu^s(\alpha_{em}^2) = (53.1 \pm 2.5) \cdot 10^{-10}$, $a_\mu^s(\alpha_{em}^3) = (-0.018 \pm 0.011) \cdot 10^{-10}$ and $a_\mu^c(\alpha_{em}^2) = (14.75 \pm 0.56) \cdot 10^{-10}$, $a_\mu^c(\alpha_{em}^3) = (-0.030 \pm 0.013) \cdot 10^{-10}$ for the strange and charm contributions, respectively.
  • We present a lattice determination of the vector and scalar form factors of the $D \to \pi(K) \ell \nu$ semileptonic decays, which are relevant for the extraction of the CKM matrix elements $|V_{cd}|$ and $|V_{cs}|$ from experimental data. Our analysis is based on the gauge configurations produced by the European Twisted Mass Collaboration with $N_f = 2+1+1$ flavors of dynamical quarks, at three different values of the lattice spacing and with pion masses as small as 210 MeV. The matrix elements of both vector and scalar currents are determined for a plenty of kinematical conditions in which parent and child mesons are either moving or at rest. Lorentz symmetry breaking due to hypercubic effects is clearly observed in the data and included in the decomposition of the current matrix elements in terms of additional form factors. After the extrapolations to the physical pion mass and to the continuum limit we determine the vector and scalar form factors in the whole kinematical region from $q^2 = 0$ up to $q^2_{max} = (M_D - M_{\pi(K)})^2$ accessible in the experiments, obtaining a good overall agreement with experiments, except in the region at high values of $q^2$ where some deviations are visible. A set of synthetic data points, representing our results for $f_+^{D \pi(K)}(q^2)$ and $f_0^{D \pi(K)}(q^2)$ for several selected values of $q^2$, is provided and also the corresponding covariance matrix is available. At zero 4-momentum transfer we get: $f_+^{D \to \pi}(0) = 0.612 ~ (35)$ and $f_+^{D \to K}(0) = 0.765 ~ (31)$. Using the experimental averages for $|V_{cd}| f_+^{D \to \pi}(0)$ and $|V_{cs}| f_+^{D \to K}(0)$, we extract $|V_{cd}| = 0.2330 ~ (137)$ and $|V_{cs}| = 0.945 ~ (38)$, respectively. The second-row of the CKM matrix is found to be in agreement with unitarity within the current uncertainties: $|V_{cd}|^2 + |V_{cs}|^2 + |V_{cb}|^2 = 0.949 ~ (78)$.
  • We present a lattice calculation of the masses and decay constants of $D_{(s)}^*$ and $B_{(s)}^*$ mesons using the gauge configurations produced by the European Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ dynamical quarks at three values of the lattice spacing $a \sim (0.06 - 0.09)$ fm. Pion masses are simulated in the range $M_\pi \simeq (210 - 450)$ MeV, while the strange and charm sea-quark masses are close to their physical values. We compute the ratios of vector to pseudoscalar masses and decay constants for various values of the heavy-quark mass $m_h$ in the range $0.7 m_c^{\mathrm{phys}} \lesssim m_h \lesssim 3 m_c^{\mathrm{phys}}$. In order to reach the physical b-quark mass, we exploit the Heavy Quark Effective Theory prediction that, in the static limit of infinite heavy-quark mass, the considered ratios are equal to one. At the physical point our results are: $M_{D^*} / M_{D} = 1.0769(79)$, $M_{D^*_{s}} / M_{D_{s}} = 1.0751(56)$, $ f_{D^*} / f_{D} = 1.078(36),$ $f_{D^*_{s}} / f_{D_{s}} = 1.087(20)$, $M_{B^*} / M_{B} = 1.0078(15)$, $M_{B^*_{s}} / M_{B_{s}} = 1.0083(10)$, $f_{B^*} / f_{B} = 0.958(22)$ and $f_{B^*_{s}} / f_{B_{s}} = 0.974(10)$. Combining them with the experimental values of the pseudoscalar meson masses (used as input to fix the quark masses) and the values of the pseudoscalar decay constants calculated by ETMC, we get: $M_{D^*} = 2013(14)~\mathrm{MeV}$, $ M_{D_s^*} = 2116(11)~\mathrm{MeV}$, $f_{D^*} = 223.5(8.4)~\mathrm{MeV}$, $f_{D_s^*} = 268.8(6.6)~\mathrm{MeV}$, $M_{B^*} = 5320.5(7.6)~\mathrm{MeV}$, $M_{B_s^*} = 5411.36(5.3)~\mathrm{MeV}$, $f_{B^*} = 185.9(7.2)~\mathrm{MeV}$ and $f_{B_s^*} = 223.1(5.4)~\mathrm{MeV}$.
  • We present a lattice computation of the isospin-breaking corrections to pseudoscalar meson masses using the gauge configurations produced by the European Twisted Mass collaboration with $N_f = 2 + 1 + 1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082$ and $0.089$ fm) with pion masses in the range $M_\pi \simeq 210 - 450$ MeV. The strange and charm quark masses are tuned at their physical values. We adopt the RM123 method based on the combined expansion of the path integral in powers of the $d$- and $u$-quark mass difference ($\widehat{m}_d - \widehat{m}_u$) and of the electromagnetic coupling $\alpha_{em}$. Within the quenched QED approximation, which neglects the effects of the sea-quark charges, and after the extrapolations to the physical pion mass and to the continuum and infinite volume limits, we provide results for the pion, kaon and (for the first time) charmed-meson mass splittings, for the prescription-dependent parameters $\epsilon_{\pi^0}$, $\epsilon_\gamma(\overline{MS}, 2~\mbox{GeV})$, $\epsilon_{K^0}(\overline{MS}, 2~\mbox{GeV})$, related to the violations of the Dashen's theorem, and for the light quark mass difference $(\widehat{m}_d - \widehat{m}_u)(\overline{MS}, 2~\mbox{GeV})$.
  • We present a precise lattice computation of pseudoscalar and vector heavy-light meson masses for heavy-quark masses ranging from the physical charm mass up to $\simeq 4$ times the physical b-quark mass. We employ the gauge configurations generated by the European Twisted Mass Collaboration (ETMC) with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) with pion masses in the range $M_\pi \simeq 210 - 450$ MeV. The heavy-quark mass is simulated directly on the lattice up to $\simeq 3$ times the physical charm mass. The interpolation to the physical $b$-quark mass is performed using the ETMC ratio method, based on ratios of the meson masses computed at nearby heavy-quark masses, and adopting the kinetic mass scheme. The extrapolation to the physical pion mass and to the continuum limit yields $m_b^{\rm kin}(1~\mbox{GeV}) = 4.61 (20)$ GeV, which corresponds to $\overline{m}_b(\overline{m}_b) = 4.26 (18)$ GeV in the $\overline{MS}$ scheme. The lattice data are analyzed in terms of the Heavy Quark Expansion (HQE) and the matrix elements of dimension-4 and dimension-5 operators are extracted with a good precision, namely: $\overline{\Lambda} = 0.552 (26)$ GeV, $\mu_\pi^2 = 0.321 (32)$ GeV$^2$ and $\mu_G^2(m_b) = 0.253 (25)$ GeV$^2$. The data also allow for a rough estimate of the dimension-6 operator matrix elements. As the HQE parameters play a crucial role in the inclusive determination of the Cabibbo-Kobayashi-Maskawa matrix elements $V_{ub}$ and $V_{cb}$, their precise determination on the lattice may eventually validate and improve the analyses based on fits to the semileptonic moments.
  • Leptonic and semileptonic kaon decays represent till now the golden modes for the extraction of the Cabibbo angle from experiments, provided the relevant hadronic quantities, namely the ratio of the kaon and pion leptonic decay constants, $f_K / f_\pi$, and the semileptonic vector form factor at zero four-momentum transfer, $f_+(0)$, are computed accurately from QCD. In the last years, using large-scale lattice QCD simulations, the determination of both $f_K / f_\pi$ and $f_+(0)$ has reached an impressive level of precision ($\approx 0.3 \%$), which is expected to be further improved in the near future. However, at a permille level of precision both electromagnetic and strong isospin-breaking effects cannot be neglected anymore. In this contribution a new lattice strategy aiming at determining QCD+QED effects in the inclusive leptonic decay rates of charged pseudoscalar mesons is briefly illustrated, and the preliminary results obtained in the case of $\pi_{\ell 2}$ and $K_{\ell 2}$ decays are presented. As for the semileptonic $K_{\ell 3}$ decays, it is pointed out that the determination of the relevant form factors in the full kinematical range covered by the experiments is necessary for a more stringent test of the Standard Model. In this respect the perspectives based on the only existing results from the European Twisted Mass Collaboration (ETMC) are quite encouraging.
  • In Carrasco et al. we have recently proposed a method to calculate $O(e^2)$ electromagnetic corrections to leptonic decay widths of pseudoscalar mesons. The method is based on the observation that the infrared divergent contributions (that appear at intermediate stages of the calculation and that cancel in physical quantities thanks to the Bloch-Nordsieck mechanism) are universal, i.e. depend on the charge and the mass of the meson but not on its internal structure. In this talk we perform a detailed analysis of the finite-volume effects associated with our method. In particular we show that also the leading $1/L$ finite-volume effects are universal and perform an analytical calculation of the finite-volume leptonic decay rate for a point-like meson.
  • We demonstrate that the leading and next-to-leading finite-volume effects in the evaluation of leptonic decay widths of pseudoscalar mesons at $O(\alpha)$ are universal, i.e. they are independent of the structure of the meson. This is analogous to a similar result for the spectrum but with some fundamental differences, most notably the presence of infrared divergences in decay amplitudes. The leading non-universal, structure-dependent terms are of $O(1/L^2)$ (compared to the $O(1/L^3)$ leading non-universal corrections in the spectrum). We calculate the universal finite-volume effects, which requires an extension of previously developed techniques to include a dependence on an external three-momentum (in our case, the momentum of the final state lepton). The result can be included in the strategy proposed in Ref.\,\cite{Carrasco:2015xwa} for using lattice simulations to compute the decay widths at $O(\alpha)$, with the remaining finite-volume effects starting at order $O(1/L^2)$. The methods developed in this paper can be generalised to other decay processes, most notably to semileptonic decays, and hence open the possibility of a new era in precision flavour physics.
  • We present a lattice calculation of the decay constants and masses of $D_{(s)}^*$ and $B_{(s)}^*$ mesons using the gauge configurations produced by the European Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ dynamical quarks and at three values of the lattice spacing $a \sim 0.06 - 0.09$ fm. Pion masses are simulated in the range $m_{\pi} \sim 210 - 450$ MeV, while the strange and charm quark masses are close to their physical values. We computed the ratios of vector to pseudoscalar decay constants or masses for various values of the heavy-quark mass $m_h$ in the range $0.7 m_c^{phys} \lesssim m_h \lesssim 3 m_c^{phys}$. In order to reach the physical b-quark mass, we exploited the HQET prediction that, in the static limit of infinite heavy-quark mass, all the considered ratios are equal to one. We obtain: $ f_{D^*}/f_{D} = 1.078(36),$ $m_{D^*}/m_{D} = 1.0769(79)$, $f_{D^*_{s}}/f_{D_{s}} = 1.087(20)$, $m_{D^*_{s}}m_{D_{s}} = 1.0751(56)$, $f_{B^*}/f_{B} = 0.958(22)$, $m_{B^*}/m_{B} = 1.0078(15)$, $f_{B^*_{s}}/f_{B_{s}} = 0.974(10)$ and $m_{B^*_{s}}/m_{B_{s}} = 1.0083(10)$. Combining them with the corresponding experimental masses from the PDG and the pseudoscalar decay constants calculated by ETMC, we get: $f_{D^*} = 223.5(8.4)~\mathrm{MeV}$, $m_{D^*} = 2013(14)~\mathrm{MeV}$, $f_{D^*_{s}} = 268.8(6.6)~\mathrm{MeV}$, $ m_{D^*_{s}} = 2116(11)~\mathrm{MeV}$, $f_{B^*} = 185.9(7.2)~\mathrm{MeV}$, $m_{B^*} = 5320.5(7.6)~\mathrm{MeV}$, $f_{B^*_{s}} = 223.1(5.4)~\mathrm{MeV}$ and $m_{B^*_{s}}= 5411.36(5.3)~\mathrm{MeV}$.
  • We present a lattice determination of the vector and scalar form factors of the semileptonic $D \to \pi \ell \nu$ decays, which are relevant for the extraction of the CKM matrix element $\lvert V_{cd} \rvert$ from experimental data. Our analysis is based on the gauge configurations produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ flavors of dynamical quarks. We simulated at three different values of the lattice spacing and with pion masses as small as 210 MeV. Quark momenta are injected on the lattice using non-periodic boundary conditions. The matrix elements of both vector and scalar currents are determined for a plenty of kinematical conditions in which parent and child mesons are either moving or at rest. Lorentz symmetry breaking due to hypercubic effects is clearly observed in the data and included in the decomposition of the current matrix elements in terms of additional form factors. Our preliminary estimate for the vector form factor at zero 4-momentum transfer is $f_+^{D \to \pi}(0) = 0.631\,(40)$, which can be compared with the latest FLAG average $f_+^{D \to \pi}(0) = 0.666\,(29)$ available only at $N_f = 2 + 1$.
  • Electromagnetic effects in the leptonic decay rates $\pi^+ \to \mu^+ \nu$ and $K^+ \to \mu^+ \nu$ are evaluated for the first time on the lattice. Following a method recently proposed in Ref. [1] the emission of virtual photons at leading order in the electromagnetic coupling is evaluated on the lattice and the infrared divergence computed for a point-like meson at finite lattice volume is subtracted. The physical decay rate is then obtained by adding the emission of real and virtual photons regularised with a photon mass. Using the gauge ensembles produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks the feasibility of our approach is demonstrated. Preliminary results for the electromagnetic corrections to charged (neutral) pion and kaon masses as well as to the leptonic decay rates of charged pions and kaons are presented.
  • We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.
  • We present precise lattice computations for the b-quark mass, the quark mass ratios mb/mc and mb/ms as well as the leptonic B-decay constants. We employ gauge configurations with four dynamical quark flavors, up/down, strange and charm, at three values of the lattice spacing (a ~ 0.06 - 0.09 fm) and for pion masses as low as 210 MeV. Interpolation in the heavy quark mass to the bottom quark point is performed using ratios of physical quantities computed at nearby quark masses exploiting the fact that these ratios are exactly known in the static quark mass limit. Our results are also extrapolated to the physical pion mass and to the continuum limit and read: mb(MSbar, mb) = 4.26(10) GeV, mb/mc = 4.42(8), mb/ms = 51.4(1.4), fBs = 229(5) MeV, fB = 193(6) MeV, fBs/fB = 1.184(25) and (fBs/fB)/(fK/fpi) = 0.997(17).
  • We present a lattice QCD determination of the vector and scalar form factors of the semileptonic $K \to \pi \ell \nu$ decay which are relevant for the extraction of the CKM matrix element $|V_{us}|$ from experimental data. Our results are based on the gauge configurations produced by the European Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical fermions, which include in the sea, besides two light mass degenerate quarks, also the strange and the charm quarks. We use data simulated at three different values of the lattice spacing and with pion masses as small as $210$ MeV. Our final result for the vector form factor at zero momentum transfer is $f_+(0) = 0.9709 (46)$, where the uncertainty is both statistical and systematic combined in quadrature. Using the latest experimental value of $f_+(0) |V_{us}|$ from $K_{\ell 3}$ decays, we obtain $|V_{us}| = 0.2230 (11)$, which allows to test the unitarity constraint of the Standard Model below the permille level once the determination of $|V_{ud}|$ from superallowed nuclear $\beta$ decays is adopted. A slight tension with unitarity at the level of $\sim 2$ standard deviations is observed. Moreover we present our results for the semileptonic scalar $f_0(q^2)$ and vector $f_+(q^2)$ form factors in the whole range of values of the squared four-momentum transfer $q^2$ measured in $K_{\ell 3}$ decays, obtaining a very good agreement with the momentum dependence of the experimental data. We provide a set of synthetic data points representing our results for the vector and scalar form factors at the physical point for several selected values of $q^2$.
  • We present a lattice QCD determination of the vector and scalar form factors of the semileptonic decays D->\pi l \nu and D -> K l \nu which are relevant for the extraction of the CKM matrix elements |Vcd| and |Vcs| from experimental data. Our analysis is based on the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2+1+1 dynamical fermions. We simulated at three different values of the lattice spacing and with pion masses as small as 210 MeV. Our preliminary estimates for the vector form factor at zero 4-momentum transfer are f+(D -> \pi)(0) = 0.610 (23) and f+(D -> K)(0) = 0.747 (22), where the uncertainties are only statistical. By combining our results with the experimental values of f+(D -> \pi)(0) |Vcd| and f+(D -> K)(0) |Vcs| we obtain |Vcd| = 0.2336 (93) and |Vcs| = 0.975 (30), which together with the PDG determination of |Vcb| are in agreement with the unitarity constraint of the Standard Model.
  • We present a lattice QCD determination of the vector and scalar form factors of the kaon semileptonic decay $K \to \pi \ell \nu$, which is relevant for the determination of the CKM matrix element $|V_{us}|$ from experimental data. Our results are based on the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2+1+1 dynamical fermions. We simulated at three different values of the lattice spacing and with pion masses as small as 210 MeV. Our estimate for the vector form factor at zero 4-momentum transfer is $f_+(0) = 0.9709 (46)$, where the uncertainty is both statistical and systematic. By combining our result with the latest experimental value of $f_+(0)|V_{us}|$ we obtain $|V_{us}| = 0.2230 (11)$, which satisfies the unitarity constraint of the Standard Model at the permille level using the updated determination of $|V_{ud}|$ coming from superallowed nuclear $\beta$ decays. We present also the momentum dependence of the vector and scalar form factors in the whole range of values of the squared 4-momentum transfer measured in $K_{\ell 3}$ decays, obtaining a good agreement with the experimental data.
  • We present unquenched lattice QCD results for the matrix elements of four-fermion operators relevant to the description of the neutral K and D mixing in the Standard Model and its extensions. We have employed simulations with Nf = 2 + 1 + 1 dynamical sea quarks at three values of the lattice spacings in the interval 0.06 - 0.09 fm and pseudoscalar meson masses in the range 210 - 450 MeV. Our results are extrapolated to the continuum limit and to the physical pion mass. Renormalization constants have been determined non-perturbatively in the RI-MOM scheme. In particular, for the Kaon bag-parameter, which is relevant for the \overline{K}^0-K^0 mixing in the Standard Model, we obtain B_K^{RGI} = 0.717(24).
  • We present our study of the renormalization of the chromomagnetic operator,O(CM), which appears in the effective Hamiltonian describing Delta S = 1 transitions in and beyond the Standard Model. We have computed, perturbatively to one-loop, the relevant Green's functions with two (quark-quark) and three (quark-quark-gluon) external fields, at nonzero quark masses, using both the lattice and dimensional regularizations. The perturbative computation on the lattice is carried out using the maximally twisted-mass action for the fermions, while for the gluons we employed the Symanzik improved gauge action for different sets of values of the Symanzik coefficients. We have identified all the operators which can possibly mix with O(CM), including lower dimensional and non gauge invariant operators, and we have calculated those elements of the mixing matrix which are relevant for the renormalization of O(CM). We have also performed numerical lattice calculations to determine non-perturbatively the mixings of the chromomagnetic operator with lower dimensional operators, through proper renormalization conditions. For the first time the 1/a**2-divergent mixing of the chromomagnetic operator with the scalar density has been determined non-perturbatively with high precision. Moreover, the 1/a-divergent mixing with the pseudoscalar density, due to the breaking of parity within the twisted-mass regularization of QCD, has been calculated non-perturbatively and found to be smaller than its one-loop perturbative estimate. The QCD simulations have been carried out using the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 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.