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

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

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 results of lattice QCD simulations with massdegenerate up and
down and masssplit strange and charm (Nf=2+1+1) dynamical quarks using Wilson
twisted mass fermions at maximal twist. The tuning of the strange and charm
quark masses is performed at three values of the lattice spacing a~0.06 fm,
a~0.08 fm and a~0.09 fm with lattice sizes ranging from L~1.9 fm to L~3.9 fm.
We perform a preliminary study of SU(2) chiral perturbation theory by combining
our lattice data from these three values of the lattice spacing.

We present a comprehensive investigation of light meson physics using
maximally twisted mass fermions for two massdegenerate quark flavours. By
employing four values of the lattice spacing, spatial lattice extents ranging
from 2.0 fm to 2.5 fm and pseudo scalar masses in the range 280 MeV to 650 MeV
we control the major systematic effects of our calculation. This enables us to
confront our data with chiral perturbation theory and extract low energy
constants of the effective chiral Lagrangian and derived quantities, such as
the light quark mass, with high precision.

We present the results of a lattice QCD calculation of the pseudoscalar meson
decay constants fpi, fK, fD and fDs, performed with Nf=2 dynamical fermions.
The simulation is carried out with the treelevel improved Symanzik gauge
action and with the twisted mass fermionic action at maximal twist. We have
considered for the final analysis three values of the lattice spacing, a~0.10
fm, 0.09 fm and 0.07 fm, with pion masses down to mpi~270 MeV. Our results for
the light meson decay constants are fK=158.1(2.4) MeV and fK/fpi=1.210(18).
From the latter ratio, by using the experimental determination of
Gamma(K>mu nu_mu (gamma))/ Gamma(pi> mu nu_mu (gamma)) and the average
value of Vud from nuclear beta decays, we obtain Vus=0.2222(34), in good
agreement with the determination from semileptonic Kl3 decays and the unitarity
constraint. For the D and Ds meson decay constants we obtain fD=197(9) MeV,
fDs=244(8) MeV and fDs/fD=1.24(3). Our result for fD is in good agreement with
the CLEO experimental measurement. For fDs our determination is smaller than
the PDG 2008 experimental average but in agreement with a recent improved
measurement by CLEO at the 1.4 sigma level.

In a recent paper [heplat/0701012] we presented precise lattice QCD results
of our European Twisted Mass Collaboration (ETMC). They were obtained by
employing two massdegenerate flavours of twisted mass fermions at maximal
twist. In the present paper we give details on our simulations and the
computation of physical observables. In particular, we discuss the problem of
tuning to maximal twist, the techniques we have used to compute correlators and
error estimates. In addition, we provide more information on the algorithm
used, the autocorrelation times and scale determination, the evaluation of
disconnected contributions and the description of our data by means of chiral
perturbation theory formulae.

We present results on the mass of the nucleon and the Delta using two
dynamical degenerate twisted mass quarks and the treelevel Symanzik improved
gauge action. The evaluation is performed at four quark masses corresponding to
a pion mass in the range of about 300600 MeV on lattices of 2.12.7 fm. We
check for cutoff effects by evaluating these baryon masses on lattices of
spatial size 2.1 fm at beta=3.9 and beta=4.05 and on a lattice of 2.4 fm at
beta=3.8. The values we find are compatible within our statistical errors.
Lattice results are extrapolated to the physical limit using continuum chiral
perturbation theory. Performing a combined fit to our lattice data at beta=3.9
and beta=4.05 we find a nucleon mass of 964\pm 28 (stat.) \pm 8 (syst.) MeV.
The nucleon mass at the physical point provides an independent determination of
the lattice spacing. Using heavy baryon chiral perturbation theory at O(p^3) we
find a_{\beta=3.9}=0.0890\pm 0.0039(stat.) \pm 0.0014(syst.) fm, and
a_{\beta=4.05}= 0.0691\pm 0.0034(stat.) \pm 0.0010(syst.) fm, in good agreement
with the values determined from the pion decay constant. Isospin violating
lattice artifacts in the Deltasystem are found to be compatible with zero for
the values of the lattice spacings used in this work. Performing a combined fit
to our lattice data at beta=3.9 and beta=4.05 we find for the masses of the
Delta^{++,} and Delta^{+,0} 1316 \pm 60 (stat.) MeV and 1330 \pm 74 (stat.)
MeV respectively. We confirm that in the continuum limit they are also
degenerate.

In this proceeding contribution we report on a first study in order to
explore the so called epsilon regime with Wilson twisted mass (Wtm) fermions.
To show the potential of this approach we give a preliminary determination of
the chiral condensate.

Using a canonical formalism, we determine the equation of state and the phase
diagram of eightflavour QCD, as a function of temperature and isospin density.
Two mechanisms are at work: Bose condensation of pions at high density, and
deconfinement at high temperature. We study their interplay and find that on
our small and coarse lattice the first order deconfinement transition appears
to end at a critical point at finite density. We investigate the strength of
the overlap and of the sign problems and discuss implications for the baryonic
density case.

We report on the scaling of basic hadronic observables in lattice QCD with
N_f=2 maximally twisted Wilson dynamical quarks. We give preliminary results
for some of the GasserLeutwyler low energy constants, the chiral condensate
and the average mass of u and d quarks.

We present the results of a lattice QCD calculation of the average updown
and strange quark masses and of the light meson pseudoscalar decay constants
with Nf=2 dynamical fermions. The simulation is carried out at a single value
of the lattice spacing with the twisted mass fermionic action at maximal twist,
which guarantees automatic O(a)improvement of the physical quantities. Quark
masses are renormalized by implementing the nonperturbative RIMOM
renormalization procedure. Our results for the light quark masses are
m_ud^{msbar}(2 GeV)= 3.85 + 0.12 + 0.40 MeV, m_s^{msbar}(2 GeV) = 105 + 3 +
9 MeV and m_s/m_ud = 27.3 + 0.3 + 1.2. We also obtain fK = 161.7 + 1.2 +
3.1 MeV and the ratio fK/fpi=1.227 + 0.009 + 0.024. From this ratio, by using
the experimental determination of Gamma(K> mu nu (gamma))/Gamma(pi > mu nu
(gamma)) and the average value of Vud from nuclear beta decays, we obtain
Vus=0.2192(5)(45), in agreement with the determination from Kl3 decays and
the unitarity constraint.

The effect of Stout smearing is investigated in numerical simulations with
twisted mass Wilson quarks. The phase transition near zero quark mass is
studied on 12^3x24, 16^3x32 and 24^3x48 lattices at lattice spacings a = 0.1 
0.125 fm.

We present results of dynamical simulations with 2 flavours of degenerate
Wilson twisted mass quarks at maximal twist in the range of pseudo scalar
masses from 300 to 550 MeV. The simulations are performed at one value of the
lattice spacing a \lesssim 0.1 fm. In order to have O(a) improvement and aiming
at small residual cutoff effects, the theory is tuned to maximal twist by
requiring the vanishing of the untwisted quark mass. Precise results for the
pseudo scalar decay constant and the pseudo scalar mass are confronted with
chiral perturbation theory predictions and the low energy constants F,
\bar{l}_3 and \bar{l}_4 are evaluated with small statistical errors.

We present a comparison of a number of iterative solvers of linear systems of
equations for obtaining the fermion propagator in lattice QCD. In particular,
we consider chirally invariant overlap and chirally improved Wilson (maximally)
twisted mass fermions. The comparison of both formulations of lattice QCD is
performed at four fixed values of the pion mass between 230MeV and 720MeV. For
overlap fermions we address adaptive precision and low mode preconditioning
while for twisted mass fermions we discuss even/odd preconditioning. Taking the
best available algorithms in each case we find that calculations with the
overlap operator are by a factor of 30120 more expensive than with the twisted
mass operator.

Discretisation errors in twoflavour lattice QCD with Wilsonquarks and DBW2
gauge action are investigated by comparing numerical simulation data at two
values of the bare gauge coupling. Both nonzero and zero twisted mass values
are considered. The results, including also data from simulations using the
Wilson plaquette gauge action, are compared to nexttoleading order chiral
perturbation theory formulas.

By employing a twisted mass term, we compare recent results from lattice
calculations of Nf=2 dynamical Wilson fermions with Wilson Chiral Perturbation
Theory (WChPT). The final goal is to determine some combinations of
GasserLeutwyler Low Energy Constants (LECs). A wide set of data with different
lattice spacings (a ~ 0.20.12 fm), different gauge actions (Wilson plaquette,
DBW2) and different quark masses (down to the lowest pion mass allowed by
lattice artifacts and including negative quark masses) provide a strong check
of the applicability of WChPT in this regime and the scaling behaviours in the
continuum limit.

We describe a new HMC algorithm variant we have recently introduced and
extend the published results by preliminary results of a simulation with a
pseudo scalar mass value of about 300 MeV. This new run confirms our
expectation that simulations with such pseudo scalar mass values become
feasible and affordable with our HMC variant. In addition we discuss
simulations from hot and cold starts at a pseudo scalar mass value of about 300
MeV, which we performed in order to test for possible metastabilities.

We summarize four contributions about dynamical twisted mass fermions. The
resulting report covers results for N_f=2 obtained from three different gauge
actions, namely the standard Wilson plaquette gauge action, the DBW2 and the
treelevel Symanzik improved gauge action. In addition, first results for
N_f=2+1+1 flavours of twisted mass fermions are discussed.

Twisted mass fermions allow light quarks to be explored but with the
consequence that there are mass splittings, such as between the neutral and
charged pion. Using a direct calculation of the connected neutral pion
correlator and stochastic methods to evaluate the disconnected correlations, we
determine the neutral pion mass. We explore the dependence on lattice spacing
and quark mass in quenched QCD. For dynamical QCD, we determine the sign of the
splitting which is linked, via chiral PT, to the nature of the phase transition
at small quark mass.

We present a variant of the HMC algorithm with mass preconditioning
(Hasenbusch acceleration) and multiple time scale integration. We have tested
this variant for standard Wilson fermions at beta=5.6 and at pion masses
ranging from 380 MeV to 680 MeV. We show that in this situation its performance
is comparable to the recently proposed HMC variant with domain decomposition as
preconditioner. We give an update of the ``Berlin Wall'' figure, comparing the
performance of our variant of the HMC algorithm to other published performance
data. Advantages of the HMC algorithm with mass preconditioning and multiple
time scale integration are that it is straightforward to implement and can be
used in combination with a wide variety of lattice Dirac operators.

Lattice QCD with Wilson fermions generically shows the phenomenon of a first
order phase transition. We study the phase structure of lattice QCD using
Wilson twisted mass fermions and the Wilson plaquette gauge action are used in
a range of beta values where such a first order phase transition is observed.
In particular, we investigate the dependence of the first order phase
transition on the value of the lattice spacing. Using only data in one phase
and neglecting possible problems arising from the phase transition we are able
to perform a first scaling test for physical quantities using this action.