• 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 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})$.
  • 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.
  • We present a lattice QCD calculation of the pseudoscalar decay constants fK, fD and fDs performed 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 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: fK+ / fpi+ = 1.184 (16), fK+ = 154.4 (2.0) MeV, which incorporate the leading strong isospin breaking correction due to the up- and down-quark mass difference, and fK = 155.0 (1.9) MeV, fD = 207.4 (3.8) MeV, fDs = 247.2 (4.1) MeV, fDs / fD = 1.192 (22) and (fDs / fD) / (fK / fpi) = 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 Ds-mesons our results lead to the following determination of the CKM matrix elements: |Vus| = 0.2269 (29), |Vcd| = 0.2221 (67) and |Vcs| = 1.014 (24). Using the latest value of |Vud| from superallowed nuclear beta decays the unitarity of the first row of the CKM matrix is fulfilled at the permille level.
  • In this paper, for the first time a method is proposed to compute electromagnetic effects in hadronic processes using lattice simulations. The method can be applied, for example, to the leptonic and semileptonic decays of light or heavy pseudoscalar mesons. For these quantities the presence of infrared divergences in intermediate stages of the calculation makes the procedure much more complicated than is the case for the hadronic spectrum, for which calculations already exist. In order to compute the physical widths, diagrams with virtual photons must be combined with those corresponding to the emission of real photons. Only in this way do the infrared divergences cancel as first understood by Bloch and Nordsieck in 1937. We present a detailed analysis of the method for the leptonic decays of a pseudoscalar meson. The implementation of our method, although challenging, is within reach of the present lattice technology.
  • The Chromomagnetic operator (CMO) mixes with a large number of operators under renormalization. We identify which operators can mix with the CMO, at the quantum level. Even in dimensional regularization (DR), which has the simplest mixing pattern, the CMO mixes with a total of 9 other operators, forming a basis of dimension-five, Lorentz scalar operators with the same flavor content as the CMO. Among them, there are also gauge noninvariant operators; these are BRST invariant and vanish by the equations of motion, as required by renormalization theory. On the other hand using a lattice regularization further operators with $d \leq 5$ will mix; choosing the lattice action in a manner as to preserve certain discrete symmetries, a minimul set of 3 additional operators (all with $d<5$) will appear. In order to compute all relevant mixing coefficients, we calculate the quark-antiquark (2-pt) and the quark-antiquark-gluon (3-pt) Green's functions of the CMO at nonzero quark masses. These calculations were performed in the continuum (dimensional regularization) and on the lattice using the maximally twisted mass fermion action and the Symanzik improved gluon action. In parallel, non-perturbative measurements of the $K-\pi$ matrix element are being performed in simulations with 4 dynamical ($N_f = 2+1+1$) twisted mass fermions and the Iwasaki improved gluon action.
  • We present preliminary results of the first lattice QCD calculation of the K -> pi matrix elements of the chromomagnetic operator O_{CM}=g sbar sigma_{munu} G_{munu} 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 computed 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 preliminary result for the B-parameter of the chromomagnetic operator is B_{CMO}=0.29(11), which can be compared with the estimate B_{CMO}~1-4 currently used in phenomenological analyses.