• We present an investigation of the electromagnetic pion form factor, $F_\pi(Q^2)$, at small values of the four-momentum transfer $Q^2$ ($\lesssim 0.25$ GeV$^2$), based on the gauge configurations generated by European Twisted Mass Collaboration with $N_f = 2$ twisted-mass quarks at maximal twist including a clover term. Momentum is injected using non-periodic 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 next-to-leading order in the light-quark 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 low-energy 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 scale-independent 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 B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark 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 low-energy 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 heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle 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 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.
  • 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 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.
  • 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 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 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.
  • 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 eight-flavour 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 Gasser-Leutwyler 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 up-down 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 non-perturbative RI-MOM 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 30-120 more expensive than with the twisted mass operator.
  • Discretisation errors in two-flavour lattice QCD with Wilson-quarks and DBW2 gauge action are investigated by comparing numerical simulation data at two values of the bare gauge coupling. Both non-zero and zero twisted mass values are considered. The results, including also data from simulations using the Wilson plaquette gauge action, are compared to next-to-leading 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 Gasser-Leutwyler Low Energy Constants (LECs). A wide set of data with different lattice spacings (a ~ 0.2-0.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 meta-stabilities.
  • 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 tree-level 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.