
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 nonperturbatively, 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 Bparameter 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 modeldependent 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 fourmomentum 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 nonzero 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 leadingorder electromagnetic and strong isospinbreaking 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 treelevel $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 isospinbreaking (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 leadingorder 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 leadingorder 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 4momentum 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
secondrow 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 seaquark masses are close to their physical values. We
compute the ratios of vector to pseudoscalar masses and decay constants for
various values of the heavyquark 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 bquark mass, we exploit the Heavy Quark Effective Theory
prediction that, in the static limit of infinite heavyquark 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 isospinbreaking 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 seaquark 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)
charmedmeson mass splittings, for the prescriptiondependent 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 BlochNordsieck 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
finitevolume effects associated with our method. In particular we show that
also the leading $1/L$ finitevolume effects are universal and perform an
analytical calculation of the finitevolume leptonic decay rate for a
pointlike meson.

We demonstrate that the leading and nexttoleading finitevolume 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 nonuniversal, structuredependent terms are of
$O(1/L^2)$ (compared to the $O(1/L^3)$ leading nonuniversal corrections in the
spectrum). We calculate the universal finitevolume effects, which requires an
extension of previously developed techniques to include a dependence on an
external threemomentum (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 finitevolume 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 heavyquark 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 bquark mass, we
exploited the HQET prediction that, in the static limit of infinite heavyquark
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 nonperiodic 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 4momentum
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 pointlike 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 Bmeson physics with
the aim of making them easily accessible to the particle physics community.
More specifically, we report on the determination of the lightquark 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 lowenergy 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 heavyquark sector, we provide results for mc and mb (also new
compared to the previous review), as well as those for D and Bmeson decay
constants, form factors, and mixing parameters. These are the heavyquark
quantities most relevant for the determination of CKM matrix elements and the
global CKM unitaritytriangle fit. Finally, we review the status of lattice
determinations of the strong coupling constant alpha_s.

We present precise lattice computations for the bquark mass, the quark mass
ratios mb/mc and mb/ms as well as the leptonic Bdecay 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 fourmomentum 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 4momentum 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 4momentum 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 4momentum 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
fourfermion 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
nonperturbatively in the RIMOM scheme. In particular, for the Kaon
bagparameter, which is relevant for the \overline{K}^0K^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 oneloop, the relevant Green's functions with two
(quarkquark) and three (quarkquarkgluon) 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
twistedmass 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 nonperturbatively the mixings of the chromomagnetic operator with
lower dimensional operators, through proper renormalization conditions. For the
first time the 1/a**2divergent mixing of the chromomagnetic operator with the
scalar density has been determined nonperturbatively with high precision.
Moreover, the 1/adivergent mixing with the pseudoscalar density, due to the
breaking of parity within the twistedmass regularization of QCD, has been
calculated nonperturbatively and found to be smaller than its oneloop
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 massdegenerate
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 downquark 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 Dsmesons 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 dimensionfive, 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 quarkantiquark (2pt) and the
quarkantiquarkgluon (3pt) 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, nonperturbative
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 oneloop in perturbation theory. The power divergent
coefficients controlling the mixing with operators of lower dimension have been
computed nonperturbatively, 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 Bparameter of the chromomagnetic operator is
B_{CMO}=0.29(11), which can be compared with the estimate B_{CMO}~14 currently
used in phenomenological analyses.