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 finite-size effects by combining the Gounaris-Sakurai parameterisation of the timelike pion form factor with the L\"uscher formalism. The impact of quark-disconnected diagrams and the precision of the scale determination is discussed and included in our final result in two-flavour 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 long-distance contribution can be accurately constrained by a dedicated spectrum calculation in the iso-vector 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 time-momentum representation. We correct our results for finite-volume effects by combining the Gounaris-Sakurai 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 quark-disconnected 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 sub-percent 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 charm-quark 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 non-perturbatively using the RI-MOM 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 b-quark mass, the B and B_s decay constants, the B-mixing bag parameters for the full four-fermion 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 bag-parameters 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 four-fermion 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 dimension-7 operators which arise at O($a$) in the Symanzik expansion of the $\Delta B=2$ operators with static heavy quarks. We consider both Wilson-like and Ginsparg-Wilson 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 mass-degenerate flavours at a single lattice spacing are shown. Measurements from the light and heavy-light 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 b-quark mass, the B and Bs decay constants, the B-mixing bag-parameters for the full four-fermion 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 b-quark mass, of the B and Bs decay constants and first results for the B-meson bag parameters. For our computation we employ the so-called ratio method and our results benefit from the use of improved interpolating operators for the B-mesons. QCD calculations are performed with Nf = 2 dynamical light-quarks 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 b-quark mass, fBs = 234(6) MeV and fB = 197(10) MeV for the B-meson decay constants, BBs(mb) = 0.90(5) and BB(mb) = 0.87(5) for the B-meson 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 so-called ratio method. Using Nf=2 maximally twisted sea quarks and Osterwalder-Seiler valence quarks we achieve both O(a)-improvement and continuum like renormalization pattern. Simulations are performed at three-values 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 non-perturbative computation of the RI-MOM renormalisation constants in a mass-independent 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 RI-MOM renormalisation constants are obtained from appropriate O(a) improved estimators extrapolated to the chiral limit.

We report on an ongoing non-perturbative computation of RI-MOM 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 mass-degenerate up and down and mass-split 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 B-physics parameters (b-quark 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 heavy-light 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, B-physics 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 tree-level 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 (un-twisted) Wilson fermionic action since the two actions only differ, in the massless limit, by a chiral rotation of the quark fields. The scale-independent renormalization constants ZV, ZA and the ratio ZP/ZS have been computed using the RI-MOM approach, as well as other alternative methods. For ZA and ZP/ZS, the latter are based on both standard twisted mass and Osterwalder-Seiler 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 RI-MOM scheme. Leading discretization effects of O(g^2 a^2), evaluated in one-loop perturbation theory, are explicitly subtracted from the RI-MOM 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 light-quark 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 low-energy constants of the effective theory as well as for the light-quark 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 D-mesons 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 mass-degenerate 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 heavy-light 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 infinite-mass 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 tree-level 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 [hep-lat/0701012] we presented precise lattice QCD results of our European Twisted Mass Collaboration (ETMC). They were obtained by employing two mass-degenerate 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 tree-level Symanzik improved gauge action. The evaluation is performed at four quark masses corresponding to a pion mass in the range of about 300-600 MeV on lattices of 2.1-2.7 fm. We check for cut-off 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 Delta-system 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.