C. Urbach

Dynamical Twisted Mass Fermions with Light Quarks: Simulation and Analysis Details

In a recent paper [ETMC, P. Boucaud et al., Phys. Lett. B 650 (2007) 304, 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.

Dynamical twisted mass fermions with light quarks

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.

HMC algorithm with multiple time scale integration and mass preconditioning

Abstract 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 β = 5.6 and at pion masses ranging from 380 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.

Pseudoscalar decay constants of kaon and D-mesons from N f = 2 twisted mass Lattice QCD

We present the results of a lattice QCD calculation of the pseudoscalar meson decay constants fπ, 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 mπ 270 MeV. Our results for the light meson decay constants are fK = 158.1(2.4) MeV and fK/fπ = 1.210(18). From the latter ratio, by using the experimental determination of Γ(K → μμ(γ))/Γ(π → μμ(γ)) 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 σ level.

tmLQCD: A program suite to simulate Wilson twisted mass lattice QCD☆

We discuss a program suite for simulating Quantum Chromodynamics on a 4-dimensional space–time lattice. The basic Hybrid Monte Carlo algorithm is introduced and a number of algorithmic improvements are explained. We then discuss the implementations of these concepts as well as our parallelisation strategy in the actual simulation code. Finally, we provide a user guide to compile and run the program.

Light quark masses and pseudoscalar decay constants from Nf = 2 lattice QCD with twisted mass fermions

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 (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 mud(2 GeV) = 3.85±0.12±0.40 MeV, ms(2 GeV) = 105±3±9 MeV and ms/mud = 27.3±0.3±1.2. We also obtain fK = 161.7±1.2±3.1 MeV and the ratio fK/fπ = 1.227±0.009±0.024. From this ratio, by using the experimental determination of Γ(K → μμ(γ))/Γ(π → μμ(γ)) 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.

Light baryon masses with dynamical twisted mass fermions

We present results on the mass of the nucleon and the

Twisted mass quarks and the phase structure of lattice QCD

\ensuremath{\Delta}

Light hadrons from Nf=2+1+1 dynamical twisted mass fermions

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\char21{}600 MeV on lattices of 2.1\char21{}2.7 fm at three lattice spacings less than 0.1 fm. We check for cutoff effects by evaluating these baryon masses on lattices of spatial size 2.1 fm at

The phase structure of lattice QCD with two flavors of Wilson quarks and renormalization group improved gluons

\ensuremath{\beta}=3.9