
We present results on the light, strange and charm nucleon scalar and tensor
charges from lattice QCD, using simulations with $N_f=2$ flavors of twisted
mass Cloverimproved fermions with a physical value of the pion mass. Both
connected and disconnected contributions are included, enabling us to extract
the isoscalar, strange and charm charges for the first time directly at the
physical point. Furthermore, the renormalization is computed nonperturbatively
for both isovector and isoscalar quantities. We investigate excited state
effects by analyzing several sinksource time separations and by employing a
set of methods to probe ground state dominance. Our final results for the
scalar charges are $g_S^u = 5.20(42)(15)(12)$, $g_S^d = 4.27(26)(15)(12)$,
$g_S^s=0.33(7)(1)(4)$, $g_S^c=0.062(13)(3)(5)$ and for the tensor charges
$g_T^u = 0.782(16)(2)(13)$, $g_T^d = 0.219(10)(2)(13)$,
$g_T^s=0.00319(69)(2)(22)$, $g_T^c=0.00263(269)(2)(37)$ in the $\overline{\rm
MS}$ scheme at 2~GeV. The first error is statistical, the second is the
systematic error due to the renormalization and the third the systematic
arising from possible contamination due to the excited states.

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 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 results for the isospin0 $\pi\pi$ swave scattering length
calculated with OsterwalderSeiler valence quarks on Wilson twisted mass gauge
configurations. We use three $N_f = 2$ ensembles with unitary (valence) pion
mass at its physical value (250$\sim$MeV), at 240$\sim$MeV (320$\sim$MeV) and
at 330$\sim$MeV (400$\sim$MeV), respectively. By using the stochastic Laplacian
Heaviside quark smearing method, all quark propagation diagrams contributing to
the isospin0 $\pi\pi$ correlation function are computed with sufficient
precision. The chiral extrapolation is performed to obtain the scattering
length at the physical pion mass. Our result $M_\pi a^\mathrm{I=0}_0 =
0.198(9)(6)$ agrees reasonably well with various experimental measurements and
theoretical predictions. Since we only use one lattice spacing, certain
systematics uncertainties, especially those arising from unitary breaking, are
not controlled in our result.

We present physics results from simulations of QCD using $N_f = 2$ dynamical
Wilson twisted mass fermions at the physical value of the pion mass. These
simulations were enabled by the addition of the clover term to the twisted mass
quark action. We show evidence that compared to previous simulations without
this term, the pion mass splitting due to isospin breaking is almost completely
eliminated. Using this new action, we compute the masses and decay constants of
pseudoscalar mesons involving the dynamical up and down as well as valence
strange and charm quarks at one value of the lattice spacing, $a \approx 0.09$
fm. Further, we determine renormalized quark masses as well as their
scaleindependent ratios, in excellent agreement with other lattice
determinations in the continuum limit. In the baryon sector, we show that the
nucleon mass is compatible with its physical value and that the masses of the
$\Delta$ baryons do not show any sign of isospin breaking. Finally, we compute
the electron, muon and tau lepton anomalous magnetic moments and show the
results to be consistent with extrapolations of older ETMC data to the
continuum and physical pion mass limits. We mostly find remarkably good
agreement with phenomenology, even though we cannot take the continuum and
thermodynamic limits.

We present results for the isospin0 $\pi\pi$ swave scattering length
calculated in twisted mass lattice QCD. We use three $N_f = 2$ ensembles with
unitary pion mass at its physical value, 240~MeV and 330~MeV respectively. We
also use a large set of $N_f = 2 + 1 +1$ ensembles with unitary pion masses
varying in the range of 230~MeV  510~MeV at three different values of the
lattice spacing. A mixed action approach with the OsterwalderSeiler action in
the valence sector is adopted to circumvent the complications arising from
isospin symmetry breaking of the twisted mass quark action. Due to the
relatively large lattice artefacts in the $N_f = 2 + 1 +1$ ensembles, we do not
present the scattering lengths for these ensembles. Instead, taking the
advantage of the many different pion masses of these ensembles, we
qualitatively discuss the pion mass dependence of the scattering properties of
this channel based on the results from the $N_f = 2 + 1 +1$ ensembles. The
scattering length is computed for the $N_f = 2$ ensembles and the chiral
extrapolation is performed. At the physical pion mass, our result $M_\pi
a^\mathrm{I=0}_0 = 0.198(9)(6)$ agrees reasonably well with various
experimental measurements and theoretical predictions.

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

Taking inspiration from lattice QCD data, we argue that a finite
nonperturbative contribution to the quark mass is generated as a consequence
of the dynamical phenomenon of spontaneous chiral symmetry breaking, in turn
triggered by the explicitly breaking of chiral symmetry induced by the critical
Wilson term in the action. In pure lattice QCD this mass term cannot be
separated from the unavoidably associated linearly divergent contribution.
However, if QCD is enlarged to a theory where also a scalar field is present,
coupled to an SU(2) doublet of fermions via a Yukawa and a Wilsonlike term,
then in the phase where the scalar field takes a nonvanishing expectation
value, a dynamically generated and "naturally" light fermion mass (numerically
unrelated to the expectation value of the scalar field) is conjectured to
emerge at a critical value of the Yukawa coupling where the symmetry of the
model is maximally enhanced. Masses dynamically generated in this way display a
natural hierarchy according to which the stronger is the strongest of the
interactions the fermion is subjected to the larger is its mass.

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.

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.

In this contribution, results from $N_f=2$ lattice QCD simulations at one
lattice spacing using twisted mass fermions with a clover term at the physical
pion mass are presented. The mass splitting between charged and neutral pions
(including the disconnected contribution) is shown to be around
$20(20)~\mathrm{MeV}$. Further, a first measurement using the clover twisted
mass action of the average momentum fraction of the pion is given. Finally, an
analysis of pseudoscalar meson masses and decay constants is presented
involving linear interpolations in strange and charm quark masses. Matching to
meson mass ratios allows the calculation of quark mass ratios:
$\mu_s/\mu_l=27.63(13)$, $\mu_c/\mu_l=339.6(2.2)$ and $\mu_c/\mu_s=12.29(10)$.
From this mass matching the quantities $f_K=153.9(7.5)~\mathrm{MeV}$,
$f_D=219(11)~\mathrm{MeV}$, $f_{D_s}=255(12)~\mathrm{MeV}$ and
$M_{D_s}=1894(93)~\mathrm{MeV}$ are determined without the application of
finite volume or discretization artefact corrections and with errors dominated
by a preliminary estimate of the lattice spacing.

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 study matrix elements of the "chromomagnetic" operator on the lattice.
This operator is contained in the strangenesschanging effective Hamiltonian
which describes electroweak effects in the Standard Model and beyond.
Having dimension 5, the chromomagnetic operator is characterized by a rich
pattern of mixing with other operators of equal and lower dimensionality,
including also non gauge invariant quantities; it is thus quite a challenge to
extract from lattice simulations a clear signal for the hadronic matrix
elements of this operator.
We compute all relevant mixing coefficients to one loop in lattice
perturbation theory; this necessitates calculating both 2point
(quarkantiquark) and 3point (gluonquarkantiquark) Green's functions at
nonzero quark masses. We use the twisted mass lattice formulation, with
Symanzik improved gluon action.
For a comprehensive presentation of our results, along with detailed
explanations and a more complete list of references, we refer to our
forthcoming publication [1].

In this contribution, a first look at simulations using maximally twisted
mass Wilson fermions at the physical point is presented. A lattice action
including clover and twisted mass terms is presented and the Monte Carlo
histories of one run with two massdegenerate flavours at a single lattice
spacing are shown. Measurements from the light and heavylight pseudoscalar
sectors are compared to previous $N_f = 2$ results and their phenomenological
values. Finally, the strategy for extending simulations to $N_f = 2 + 1 + 1$ is
outlined.

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.