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

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.

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 recent confirmation of the charged charmonium like resonance Z(4430) by
the LHCb experiment strongly suggests the existence of QCD multi quarks bound
states. Some preliminary results about hypothetical flavored tetraquark mesons
are reported. Such states are particularly amenable to Lattice QCD studies as
their interpolating operators do not overlap with those of ordinary
hiddencharm mesons.

We attempt a description of the recently discovered Z_{c,b} states in terms
of Feshbach resonances arising from the interaction between the `closed'
subspace of hadrocharmonium levels and the `open' one of opencharm/beauty
thresholds. We show how the neutrality of the X(3872) might be understood in
this scheme and provide a preliminary explanation of the pattern of the
measured total widths of X,Z_{c,b}.

The phenomenology of the socalled X, Y and Z hadronic resonances is hard to
reconcile with standard charmonium or bottomonium interpretations. It has been
suggested that some of these new hadrons can possibly be described as tightly
bound tetraquark states and/or as loosely bound twomeson molecules. In the
present paper we focus on the hypothetical existence of flavored, doubly
charmed, tetraquarks. Such states might also carry double electric charge, and
in this case, if discovered, they could univocally be interpreted in terms of
compact tetraquarks. Flavored tetraquarks are also amenable to lattice studies
as their interpolating operators do not overlap with ordinary meson ones. We
show that doubly charmed tetraquarks could significantly be produced at LHC
from B_c or Xi_bc heavy baryons.

We present a method to evaluate on the lattice the leading isospin breaking
effects due to both the small mass difference between the up and down quarks
and the QED interaction. Our proposal is applicable in principle to any QCD+QED
gauge invariant hadronic observable which can be computed on the lattice. It is
based on the expansion of the pathintegral in powers of the small parameters
(m_d  m_u)/Lambda_{QCD} and alpha_{em}, where m_f is the renormalized quark
mass and alpha_{em} the renormalized fine structure constant. In this paper we
discuss in detail the general strategy of the method and the conventional,
although arbitrary, separation of QCD from QED isospin breaking corrections. We
obtain results for the pion mass splitting, M_{pi+}^2M_{pi0}^2= 1.44(13)(16) x
10^3 MeV^2, for the Dashen's theorem breaking parameter epsilon_{gamma}=
0.79(18)(18), for the light quark masses, [m_d  m_u](MSbar,2 GeV)= 2.39(8)(17)
MeV, [m_u / m_d](MSbar,2 GeV)= 0.50(2)(3) and for the flavour symmetry breaking
parameters R and Q. We also update our previous results for the QCD isospin
breaking corrections to the Kl2 decay rate and for the QCD contribution to the
neutronproton mass splitting.

We propose a method to expand correlation functions with respect to the
spatial components of external momenta. From the coefficients of the expansion
it is possible to extract Lorentzinvariant form factors at zero spatial
momentum transfer avoiding model dependent extrapolations. These objects can be
profitably calculated on the lattice. We have explicitly checked the validity
of the proposed procedure by considering twopoint correlators with insertions
of the axial current, the form factors of the semileptonic decay of
pseudoscalar mesons, and the hadronic vacuum polarization tensor entering, for
example, the lattice calculation of the anomalous magnetic moment of the muon.

We present a new method to evaluate with high precision isospin breaking
effects due to the small mass difference between the up and down quarks using
lattice QCD. Our proposal is applicable in principle to any hadronic observable
which can be computed on the lattice. It is based on the expansion of the
pathintegral in powers of the small parameter mdmu. In this paper, we apply
this method to compute the leading isospin breaking effects for several
physical quantities of interest: the kaon meson masses, the kaon decay
constant, the form factors of semileptonic Kl3 decays and the neutronproton
mass splitting.

We present a new method to evaluate with high precision the isospin breaking
effects due to the mass difference between the up and down quarks using lattice
QCD. Our proposal is applicable in principle to any hadronic observable which
can be computed on the lattice. It is based on the expansion of the
pathintegral in powers of the small parameter $m_d  m_u$. In this talk we
discuss how to apply this method to compute the leading isospin breaking
effects for several physical quantities of interest: the kaon masses, the kaon
decay constants and the neutronproton mass splitting.

We propose a preconditioning of the Dirac operator based on the factorisation
of a predefined function related to the decay of the propagator with the
distance. We show that it can improve the accuracy of correlators involving
heavy quarks at large distances and accelerate the computation of light quark
propagators.

One of the major challenges of particle physics has been to gain an indepth
understanding of the role of quark flavor and measurements and theoretical
interpretations of their results have advanced tremendously: apart from masses
and quantum numbers of flavor particles, there now exist detailed measurements
of the characteristics of their interactions allowing stringent tests of
Standard Model predictions. Among the most interesting phenomena of flavor
physics is the violation of the CP symmetry that has been subtle and difficult
to explore. Till early 1990s observations of CP violation were confined to
neutral $K$ mesons, but since then a large number of CPviolating processes
have been studied in detail in neutral $B$ mesons. In parallel, measurements of
the couplings of the heavy quarks and the dynamics for their decays in large
samples of $K, D$, and $B$ mesons have been greatly improved in accuracy and
the results are being used as probes in the search for deviations from the
Standard Model. In the near future, there will be a transition from the current
to a new generation of experiments, thus a review of the status of quark flavor
physics is timely. This report summarizes the results of the current generation
of experiments that is about to be completed and it confronts these results
with the theoretical understanding of the field.

The possibility of revealing new physics by studying the flavor sector of the
Standard Model strongly depends upon the accuracy that will be achieved in
(near) future lattice QCD calculations and, in particular, on heavylight
meson's observables. In turn, handling with heavylight mesons on the lattice
is a challenging problem, because of the presence of two largely separated
energy scales, and at present it is impossible to extract matrix elements
involving B mesons in external states without recurring to some approximation.
In this note I give a fast overview of some of the methods that have been
devised to handle such kind of problems, emphasizing those based on finite
volume techniques, and briefly discuss some recent results obtained by their
application.

We calculate, in the continuum limit of quenched lattice QCD, the form factor
that enters the decay rate of the semileptonic decay B > D* l nu. By using
the step scaling method (SSM), previously introduced to handle two scale
problems in lattice QCD, and by adopting flavor twisted boundary conditions we
extract F(w) at finite momentum transfer and at the physical values of the
heavy quark masses. Our results can be used in order to extract the CKM matrix
element Vcb by the experimental decay rate without model dependent
extrapolations. The value of Vcb agrees with the one obtained from the B > D
l nu channel and makes us confident that the quenched approximation well
applies to these transitions.

We calculate, in the continuum limit of quenched lattice QCD, the form factor
that enters in the decay rate of the semileptonic decay B > D l nu. Making
use of the step scaling method (SSM), previously introduced to handle two scale
problems in lattice QCD, and of flavour twisted boundary conditions we extract
G(w) at finite momentum transfer and at the physical values of the heavy quark
masses. Our results can be used in order to extract the CKM matrix element Vcb
by the experimental decay rate without model dependent extrapolations.

We calculate, in the continuum limit of quenched lattice QCD, the matrix
elements of the heavyheavy vector current between heavylight pseudoscalar
meson states. We present the form factors for different values of the initial
and final meson masses at finite momentum transfer. In particular, we calculate
the nonperturbative correction to the differential decay rate of the process B
> D l nu including the case of a nonvanishing lepton mass.

In this second report on our recent numerical simulations of twoflavour QCD,
we provide further technical details on the simulations and describe the
methods we used to extract the meson masses and decay constants from the
generated ensembles of gauge fields. Among the topics covered are the choice of
the DDHMC parameters, the issue of stability, autocorrelations and the
statistical error analysis. Extensive data tables are included as well as a
short discussion of the quarkmass dependence in partially quenched QCD,
supplementing the physics analysis that was presented in the first paper in
this series.

Recent conceptual, algorithmic and technical advances allow numerical
simulations of lattice QCD with Wilson quarks to be performed at significantly
smaller quark masses than was possible before. Here we report on simulations of
twoflavour QCD at seaquark masses from slightly above to approximately 1/4 of
the strangequark mass, on lattices with up to 64x32^3 points and spacings from
0.05 to 0.08 fm. Physical seaquark effects are clearly seen on these lattices,
while the lattice effects appear to be quite small, even without O(a)
improvement. A striking result is that the dependence of the pion mass on the
seaquark mass is accurately described by leadingorder chiral perturbation
theory up to meson masses of about 500 MeV.

We study the spectral gap of the WilsonDirac operator in twoflavour
lattice QCD as a function of the lattice spacing $a$, the spacetime volume $V$
and the currentquark mass $m$. It turns out that the median of the probability
distribution of the gap scales proportionally to $m$ and that its width is
practically equal to $a/\sqrt{V}$. In particular, numerical simulations are
safe from accidental zero modes in the largevolume regime of QCD.

We present the current status of APENet, our custom 3dimensional
interconnect architecture for PC clusters environment. We report some
microbenchmarks on our recent large installation as well as new developments
on the software and hardware side. The low level device driver has been
reworked by following a custom hardware RDMA architecture, and MPICHVMI, an
implementation of the MPI library, has been ported to APENet.

The adoption of two distinct boundary conditions for two fermions species on
a finite lattice allows to deal with arbitrary relative momentum between the
two particle species, in spite of the momentum quantization rule due to a
limited physical box size. We test the physical significance of this
topological momentum by checking in the continuum limit the validity of the
expected energymomentum dispersion relations.