
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 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 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 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 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 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 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.

We discuss a lattice QCD computation of the $B$meson decay constants by the
ETM collaboration where suitable ratios allow to reach the bottom quark sector
by combining simulations around the charmquark mass with an exactly known
static limit. The different steps involved in this ratio method are discussed
together with an account of the assessment of various systematic effects. A
comparison of results from simulations with two and four flavour dynamical
quarks is presented.

We present precision lattice calculations of the pseudoscalar decay constants
of the charmed sector as well as determinations of the bottom quark mass and
its ratio to the charm quark mass. We employ Nf=2+1+1 dynamical quark gauge
configurations generated by the European Twisted Mass Collaboration, using data
at three values of the lattice spacing and pion masses as low as 210 MeV.
Strange and charm sea quark masses are close to their physical values.

We present a lattice QCD calculation of the up, down, strange and charm quark
masses 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 physical values. The simulations are
based on a unitary setup for the two light quarks and on a mixed action
approach for the strange and charm quarks. The analysis uses data at three
values of the lattice spacing and pion masses in the range 210  450 MeV,
allowing for accurate continuum limit and controlled chiral extrapolation. The
quark mass renormalization is carried out nonperturbatively using the RIMOM
method. The results for the quark masses converted to the bar{MS} scheme are:
mud(2 GeV) = 3.70(17) MeV, ms(2 GeV) = 99.6(4.3) MeV and mc(mc) = 1.348(46)
GeV. We obtain also the quark mass ratios ms/mud = 26.66(32) and mc/ms =
11.62(16). By studying the mass splitting between the neutral and charged kaons
and using available lattice results for the electromagnetic contributions, we
evaluate mu/md = 0.470(56), leading to mu = 2.36(24) MeV and md = 5.03(26) MeV.

We present the first unquenched lattice QCD results for the bag parameters
controlling the short distance contribution to D meson oscillations in the
Standard Model and beyond. We have used the gauge configurations produced by
the European Twisted Mass Collaboration with Nf = 2 dynamical quarks, at four
lattice spacings and light meson masses in the range 280500 MeV.
Renormalization is carried out nonperturbatively with the RIMOM method. The
bagparameter results have been used to constrain New Physics effects in DDbar
mixing, to put a lower bound to the generic New Physics scale and to constrain
offdiagonal squark mass terms for TeVscale Supersymmetry.

We present a lattice QCD computation of the bquark mass, the B and B_s decay
constants, the Bmixing bag parameters for the full fourfermion operator basis
as well as determinations for \xi and f_{Bq}\sqrt{B_i^{(q)}} extrapolated to
the continuum limit and to the physical pion mass. We used N_f = 2 twisted mass
Wilson fermions at four values of the lattice spacing with pion masses ranging
from 280 to 500 MeV. Extrapolation in the heavy quark mass from the charm to
the bottom quark region has been carried out on ratios of physical quantities
computed at nearby quark masses, exploiting the fact that they have an exactly
known infinite mass limit. Our results are m_b(m_b,
\overline{\rm{MS}})=4.29(12) GeV, f_{Bs}=228(8) MeV, f_{B}=189(8) MeV and
f_{Bs}/f_B=1.206(24). Moreover with our results for the bagparameters we find
\xi=1.225(31), B_1^{(s)}/B_1^{(d)}=1.01(2), f_{Bd}\sqrt{\hat{B}_{1}^{(d)}} =
216(10) MeV and f_{Bs}\sqrt{\hat{B}_{1}^{(s)}} = 262(10) MeV. We also computed
the bag parameters for the complete basis of the fourfermion operators which
are required in beyond the SM theories. By using these results for the bag
parameters we are able to provide a refined Unitarity Triangle analysis in the
presence of New Physics, improving the bounds coming from B_{(s)}\bar B_{(s)}
mixing.

We update the analysis of $D$ meson mixing including the latest experimental
results as of January 2014. We derive constraints on the parameters $M_{12}$,
$\Gamma_{12}$ and $\Phi_{12}$ that describe $D$ meson mixing using all
available data, allowing for CP violation. We also provide posterior
distributions for observable parameters appearing in $D$ physics.

We present a lattice QCD calculation of the pseudoscalar decay constants
$f_K$, $f_D$ and $f_{D_s}$ performed by the European Twisted Mass Collaboration
with $N_f = 2 + 1 + 1$ dynamical fermions. We simulated at three different
values of the lattice spacing, the smallest being approximately $0.06fm$, and
with pion masses as small as $210$MeV. Our main results are:
$f_{K^+}/f_{\pi^+}=1.183(17)$, $f_{K^+}=154.4(2.1)$MeV,
$f_{D_s}=242.1(8.3)$MeV, $f_D=201.9(8.0)$MeV, $f_{D_s}/f_D=1.199(25)$ and
$(f_{D_s}/f_D) / (f_K/f_\pi) = 1.005(15)$.

We present a lattice QCD determination of the average updown, strange and
charm quark masses based on simulations performed by the European Twisted Mass
Collaboration with $N_f = 2 + 1 + 1$ dynamical fermions. We simulated at three
different values of the lattice spacing, the smallest being approximately
$0.06fm$, and with pion masses as small as $210 \text{MeV}$. Our results are:
$m_{ud}(2\text{GeV})=3.70(17)\text{MeV}$,
$m_s(2\text{GeV})=99.2(3.9)\text{MeV}$, $m_c(m_c)=1.350(49)\text{GeV}$,
$m_s/m_{ud}=26.64(30)$ and $m_c/m_s=11.65(12)$.

We present a lattice QCD determination of the $b$quark mass and of the
$f_{B_s}$ and $f_B$ decay constants performed with $N_f = 2 + 1 + 1$ twisted
mass Wilson fermions. We have used simulations at three values of the lattice
spacing generated by ETMC with pion masses ranging from 210 to 440 MeV. To
obtain physical quantities we performed a combined chiral and continuum limit
and an extrapolation in the heavy quark mass from the charm to the $b$quark
region using suitable ratios calculated at nearby heavyquark masses having an
exactly known static limit. Our results are: $m_b(m_b) = 4.29 (13)$ GeV, $f_B =
196 (9)$ MeV, $f_{B_s} = 235 (9)$ MeV, $f_{B_s} / f_B = 1.201 (25)$,
$(f_{B_s}/f_B)/(f_K/f_\pi) = 1.007 (16)$ and $(f_{B_s}/f_B)/(f_{D_s}/f_D) =
1.008 (13)$.

We present an accurate lattice QCD computation of the bquark mass, the B and
Bs decay constants, the Bmixing bagparameters for the full fourfermion
operator basis, as well as estimates for \xi and f_{Bq}\sqrt{B_q} extrapolated
to the continuum limit and the physical pion mass. We have used Nf = 2
dynamical quark gauge configurations at four values of the lattice spacing
generated by ETMC. Extrapolation in the heavy quark mass from the charm to the
bottom quark region has been carried out using ratios of physical quantities
computed at nearby quark masses, having an exactly known infinite mass limit.

We present the first unquenched, continuum limit, lattice QCD results for the
matrix elements of the operators describing neutral kaon oscillations in
extensions of the Standard Model. Owing to the accuracy of our calculation on
\Delta S=2 weak Hamiltonian matrix elements, we are able to provide a refined
Unitarity Triangle analysis improving the bounds coming from model independent
constraints on New Physics. In our nonperturbative computation we use a
combination of Nf=2 maximally twisted sea quarks and OsterwalderSeiler valence
quarks in order to achieve both O(a)improvement and continuumlike
renormalization properties for the relevant fourfermion operators. The
calculation of the renormalization constants has been performed
nonperturbatively in the RIMOM scheme. Based on simulations at four values of
the lattice spacing and a number of quark masses we have
extrapolated/interpolated our results to the continuum limit and physical
light/strange quark masses.

We present a precise lattice QCD determination of the bquark mass, of the B
and Bs decay constants and first results for the Bmeson bag parameters. For
our computation we employ the socalled ratio method and our results benefit
from the use of improved interpolating operators for the Bmesons. QCD
calculations are performed with Nf = 2 dynamical lightquarks at four values of
the lattice spacing and the results are extrapolated to the continuum limit.
The preliminary results are mb(mb) = 4.35(12) GeV for the MSbar bquark mass,
fBs = 234(6) MeV and fB = 197(10) MeV for the Bmeson decay constants, BBs(mb)
= 0.90(5) and BB(mb) = 0.87(5) for the Bmeson bag parameters.