
We present results of calculations of the hadronic vacuum polarisation
contribution to the muon anomalous magnetic moment. Specifically, we focus on
controlling the infrared regime of the vacuum polarisation function. Our
results are corrected for finitesize effects by combining the GounarisSakurai
parameterisation of the timelike pion form factor with the L\"uscher formalism.
The impact of quarkdisconnected diagrams and the precision of the scale
determination is discussed and included in our final result in twoflavour QCD,
which carries an overall uncertainty of 6%. We present preliminary results
computed on ensembles with $N_f=2+1$ dynamical flavours and discuss how the
longdistance contribution can be accurately constrained by a dedicated
spectrum calculation in the isovector channel.

We present a calculation of the hadronic vacuum polarization contribution to
the muon anomalous magnetic moment, $a_\mu^{\mathrm hvp}$, in lattice QCD
employing dynamical up and down quarks. We focus on controlling the infrared
regime of the vacuum polarization function. To this end we employ several
complementary approaches, including Pad\'e fits, time moments and the
timemomentum representation. We correct our results for finitevolume effects
by combining the GounarisSakurai parameterization of the timelike pion form
factor with the L\"uscher formalism. On a subset of our ensembles we have
derived an upper bound on the magnitude of quarkdisconnected diagrams and
found that they decrease the estimate for $a_\mu^{\mathrm hvp}$ by at most 2%.
Our final result is $a_\mu^{\mathrm hvp}=(654\pm32\,{}^{+21}_{23})\cdot
10^{10}$, where the first error is statistical, and the second denotes the
combined systematic uncertainty. Based on our findings we discuss the prospects
for determining $a_\mu^{\mathrm hvp}$ with subpercent precision.

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 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 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 determine a basis of dimension7 operators which arise at O($a$) in the
Symanzik expansion of the $\Delta B=2$ operators with static heavy quarks. We
consider both Wilsonlike and GinspargWilson light quarks. Exact chiral
symmetry reduces the number of these O($a$) counterterms by a factor of two.
Only a subset of these operators has previously appeared in the literature. We
then extend the analysis to the O($1/m$) operators contributing beyond the
static approximation.

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

We present the ETMC results for the bag parameters describing the neutral
kaon mixing in the Standard Model and beyond and preliminary results for the
bag parameters controlling the short distance contributions in the
D^0\bar{D}^0 oscillations. We also present preliminary results for the B_{Bd},
B_{Bs}, B_{Bs}/B_{Bd} and \xi parameter controlling B^0_\bar{B}^0
oscillations in the Standard Model employing the socalled ratio method. Using
Nf=2 maximally twisted sea quarks and OsterwalderSeiler valence quarks we
achieve both O(a)improvement and continuum like renormalization pattern.
Simulations are performed at threevalues of the lattice spacing and several
values of quark masses in the light, strange, charm region and above charm up
to ~2.5m_c. Our results are extrapolated to the continuum limit and
extrapolated/interpolated to the physical quark masses.

We present a lattice QCD determination of the b quark mass and of the B and
B_s decay constants, performed with N_f=2 twisted mass Wilson fermions, by
simulating at four values of the lattice spacing. In order to study the b quark
on the lattice, two methods are adopted in the present work, respectively based
on suitable ratios with exactly known static limit and on the interpolation
between relativistic data, evaluated in the charm mass region, and the static
point, obtained by simulating the HQET on the lattice. The two methods provide
results in good agreement. For the b quark mass in the MSbar scheme and for the
decay constants we obtain m_b(m_b)=4.29(14) GeV, f_B=195(12) MeV, f_Bs=232(10)
MeV and f_Bs/f_B=1.19(5). As a byproduct of the analysis we also obtain the
results for the f_D and f_Ds decay constants: f_D=212(8) MeV, f_Ds=248(6) MeV
and f_Ds/f_D=1.17(5).

We present preliminary results of the nonperturbative computation of the
RIMOM renormalisation constants in a massindependent scheme for the action
with Iwasaki glue and four dynamical Wilson quarks employed by ETMC. Our
project requires dedicated gauge ensembles with four degenerate sea quark
flavours at three lattice spacings and at several values of the standard and
twisted quark mass parameters. The RIMOM renormalisation constants are
obtained from appropriate O(a) improved estimators extrapolated to the chiral
limit.

We report on an ongoing nonperturbative computation of RIMOM scheme
renormalization constants for the lattice action with four dynamical flavours
currently in use by ETMC. For this goal dedicated simulations with four
degenerate sea quark flavours are performed at several values of the standard
and twisted quark mass parameters. We discuss a method for removing possible
O(a) artifacts at all momenta and extrapolating renormalization constant
estimators to the chiral limit. We give preliminary results at one lattice
spacing.

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.

A method to extract Bphysics parameters (bquark mass and fB, fBs decay
constants) from currently available lattice data is presented and tested. The
approach is based on the idea of constructing appropriate ratios of heavylight
meson masses and decay constants, respectively, possessing a precisely known
static limit, and evaluating them at various pairs of heavy quark masses around
the charm. Via a smooth interpolation in the heavy quark mass from the easily
accessible charm region to the asymptotic point, Bphysics parameters are
computed with a few percent (statistical + systematic) error using recently
produced Nf=2 maximally twisted Wilson fermions data.

We present results for the renormalization constants of bilinear quark
operators obtained by using the treelevel Symanzik improved gauge action and
the Nf=2 twisted mass fermion action at maximal twist, which guarantees
automatic O(a)improvement. Our results are also relevant for the corresponding
standard (untwisted) Wilson fermionic action since the two actions only
differ, in the massless limit, by a chiral rotation of the quark fields. The
scaleindependent renormalization constants ZV, ZA and the ratio ZP/ZS have
been computed using the RIMOM approach, as well as other alternative methods.
For ZA and ZP/ZS, the latter are based on both standard twisted mass and
OsterwalderSeiler fermions, while for ZV a Ward Identity has been used. The
quark field renormalization constant Zq and the scale dependent renormalization
constants ZS, ZP and ZT are determined in the RIMOM scheme. Leading
discretization effects of O(g^2 a^2), evaluated in oneloop perturbation
theory, are explicitly subtracted from the RIMOM estimates.

We study lattice QCD determinations of the pion mass and decay constant by
means of chiral perturbation theory (ChPT). The lattice data are obtained from
large scale simulations with N_f=2 flavours of twisted mass fermions at maximal
twist. We perform a scaling test to the continuum limit of these data and use
ChPT to perform both the chiral and the infinite volume extrapolations.

We study lightquark observables by means of dynamical lattice QCD
simulations using two flavours of twisted mass fermions at maximal twist. We
employ chiral perturbation theory to describe our data for the pion mass and
decay constant. In this way, we extract precise determinations for the
lowenergy constants of the effective theory as well as for the lightquark
mass and the chiral condensate.

We present first results from runs performed with Nf=2+1+1 flavours of
dynamical twisted mass fermions at maximal twist: a degenerate light doublet
and a mass split heavy doublet. An overview of the input parameters and tuning
status of our ensembles is given, together with a comparison with results
obtained with Nf=2 flavours. The problem of extracting the mass of the K and
Dmesons is discussed, and the tuning of the strange and charm quark masses
examined. Finally we compare two methods of extracting the lattice spacings to
check the consistency of our data and we present some first results of ChiPT
fits in the light meson sector.

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 a lattice QCD calculation of the heavylight decay constants f_B
and f_Bs performed with Nf=2 maximally twisted Wilson fermions, at four values
of the lattice spacing. The decay constants have been also computed in the
static limit and the results are used to interpolate the observables between
the charm and the infinitemass sectors, thus obtaining the value of the decay
constants at the physical b quark mass. Our preliminary results are f_B=191(14)
MeV, f_Bs=243(14) MeV, f_Bs/f_B=1.27(5). They are in good agreement with those
obtained with a novel approach, recently proposed by our Collaboration (ETMC),
based on the use of suitable ratios having an exactly known static limit.

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 an update of the results for pion mass and pion decay constant as
obtained by the ETM collaboration in large scale simulations with maximally
twisted mass fermions and two mass degenerate flavours of light quarks. We
discuss the continuum, chiral and infinite volume extrapolation of these
quantities as well as the extraction of low energy constants, and investigate
possible systematic uncertainties.

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.