
We present a determination of the CabibboKobayashiMaskawa 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 4momentum 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 4momentum transfer, we obtain more precise and consistent
results: $V_{cd} = 0.2341 ~ (74)$ and $V_{cs} = 0.970 ~ (33)$. The
secondrow 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 isospinsymmetric 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 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 an investigation of the electromagnetic pion form factor,
$F_\pi(Q^2)$, at small values of the fourmomentum transfer $Q^2$ ($\lesssim
0.25$ GeV$^2$), based on the gauge configurations generated by European Twisted
Mass Collaboration with $N_f = 2$ twistedmass quarks at maximal twist
including a clover term. Momentum is injected using nonperiodic 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
nexttoleading order in the lightquark 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 lowenergy constant $\overline{\ell}_6$ is equal to
$\overline{\ell}_6 = 16.2 ~ (1.0)$.

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})$.

We present a precise lattice computation of pseudoscalar and vector
heavylight meson masses for heavyquark masses ranging from the physical charm
mass up to $\simeq 4$ times the physical bquark 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 heavyquark 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 heavyquark 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 dimension4 and dimension5 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 dimension6 operator matrix
elements. As the HQE parameters play a crucial role in the inclusive
determination of the CabibboKobayashiMaskawa 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
fourmomentum transfer, $f_+(0)$, are computed accurately from QCD. In the last
years, using largescale 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 isospinbreaking 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 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.