Search

We present a lattice determination of the inclusive decay rate of the process $\tau\mapsto X_{us} \nu_\tau$ in which the $\tau$ lepton decays into a generic hadronic state $X_{us}$ with $\bar u s$ flavour quantum numbers. Our results have been obtained in $n_f=2+1+1$ iso-symmetric QCD with full non-perturbative accuracy, without any OPE approximation and, except for the presently missing long-distance isospin-breaking corrections, include a solid estimate of all sources of theoretical uncertainties. This has been possible by using the Hansen-Lupo-Tantalo method [1] that we have already successfully applied in [2] to compute the inclusive decay rate of the process $\tau\mapsto X_{ud} \nu_\tau$ in the $\bar u d$ flavour channel. By combining our first-principles theoretical results with the presently-available experimental data we extract the CKM matrix element $\vert V_{us}\vert$, the Cabibbo angle, with a $0.9$\% accuracy, dominated by the experimental error., Comment: 9 pages, 5 figures, 1 table. Version accepted by PRL, expanded technical discussion moved to an appendix, results unchanged

The neutral pion generates the leading pole contribution to the hadronic light-by-light tensor, which is given in terms of the nonperturbative transition form factor $\mathcal{F}_{\pi^0\gamma\gamma}(q_1^2,q_2^2)$. Here we present an ab-initio lattice calculation of this quantity in the continuum and at the physical point using twisted-mass lattice QCD. We report our results for the transition form factor parameterized using a model-independent conformal expansion valid for arbitrary space-like kinematics and compare it with experimental measurements of the single-virtual form factor, the two-photon decay width, and the slope parameter. We then use the transition form factors to compute the pion-pole contribution to the hadronic light-by-light scattering in the muon $g-2$, finding $a_\mu^{\pi^0\text{-pole}} = 56.7(3.2) \times 10^{-11}$., Comment: 21 pages, 17 figures, 4 tables, updated to published version

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{\rm SD}$ and $a_{\mu}^{\rm W}$. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavours of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavours tuned close to their physical values. The simulations are carried out at three values of the lattice spacing equal to $\simeq 0.057, 0.068$ and $0.080$ fm with spatial lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{\rm SD} = 69.27\,(34) \cdot 10^{-10}$, in agreement with the dispersive determination based on experimental $e^+ e^-$ data. For the intermediate window we get instead $a_\mu^{\rm W} = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with recent determinations by other lattice collaborations, but disagrees with the dispersive determination at the level of $3.6\,\sigma$., Comment: 10 pages, 4 figures, 1 table, LATTICE 2022

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{\rm SD}$ and $a_{\mu}^{\rm W}$, defined by the RBC/UKQCD Collaboration [1]. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavors tuned close to their physical values. The simulations are carried out at three values of the lattice spacing equal to $\simeq 0.057, 0.068$ and $0.080$ fm with spatial lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{\rm SD}({\rm ETMC}) = 69.27\,(34) \cdot 10^{-10}$, which is consistent with the recent dispersive value of $a_\mu^{\rm SD}(e^+ e^-) = 68.4\,(5) \cdot 10^{-10}$ [2]. In the case of the intermediate window we get the value $a_\mu^{\rm W}({\rm ETMC}) = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with the result $a_\mu^{\rm W}({\rm BMW}) = 236.7\,(1.4) \cdot 10^{-10}$ [3] by the BMW collaboration as well as with the recent determination by the CLS/Mainz group of $a_\mu^{\rm W}({\rm CLS}) = 237.30\,(1.46) \cdot 10^{-10}$ [4]. However, it is larger than the dispersive result of $a_\mu^{\rm W}(e^+ e^-) = 229.4\,(1.4) \cdot 10^{-10}$ [2] by approximately $3.6$ standard deviations. The tension increases to approximately $4.5$ standard deviations if we average our ETMC result with those by BMW and CLS/Mainz. Our accurate lattice results in the short and intermediate windows point to a possible deviation of the $e^+ e^-$ cross section data with respect to Standard Model predictions in the low and intermediate energy regions, but not in the high energy region., Comment: 78 pages, 22 figures, 14 tables. Analysis improved with more data and fits, presentation reorganized, more material in appendices, fixed typos

We review lattice results related to pion, kaon, $D$-meson, $B$-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $f_+(0)$ arising in the semileptonic $K \to \pi$ transition at zero momentum transfer, as well as the decay constant ratio $f_K/f_\pi$ and its consequences for the CKM matrix elements $V_{us}$ and $V_{ud}$. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $SU(2)_L\times SU(2)_R$ and $SU(3)_L\times SU(3)_R$ Chiral Perturbation Theory. We review the determination of the $B_K$ parameter of neutral kaon mixing as well as the additional four $B$ parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $m_c$ and $m_b$ as well as those for the decay constants, form factors, and mixing parameters of charmed and bottom mesons and baryons. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant $\alpha_s$. We consider nucleon matrix elements, and review the determinations of the axial, scalar and tensor bilinears, both isovector and flavor diagonal. Finally, in this review we have added a new section reviewing determinations of scale-setting quantities., Comment: 418 pages, 53 figures, 200 tables, 1056 references. Minor changes, version as published in EPJC. arXiv admin note: substantial text overlap with arXiv:1902.08191, arXiv:1607.00299, arXiv:1310.8555

We present results for the light, strange and charm quark masses using $N_f=2+1+1$ twisted mass fermion ensembles at three values of the lattice spacing, including two ensembles simulated with the physical value of the pion mass. The analysis is done both in the meson and baryon sectors. The difference in the mean values found in the two sectors is included as part of the systematic error. The presentation is based on the work of Ref. [1], where more details can be found., Comment: Talk presented at the 38th International Symposium on Lattice Field Theory, Lattice 2021, Zoom/Gather@Massachusetts Institute of Technology, 26-30 Jul. 2021; 8 pages, 5 figures

We present a calculation of the up, down, strange and charm quark masses performed within the lattice QCD framework. We use the twisted mass fermion action and carry out simulations that include in the sea two light mass-degenerate quarks, as well as the strange and charm quarks. In the analysis we use gauge ensembles simulated at three values of the lattice spacing and with light quarks that correspond to pion masses in the range from 350 MeV to the physical value, while the strange and charm quark masses are tuned approximately to their physical values. We use several quantities to set the scale in order to check for finite lattice spacing effects and in the continuum limit we get compatible results. The quark mass renormalization is carried out non-perturbatively using the RI'-MOM method converted into the $\overline{\rm MS}$ scheme. For the determination of the quark masses we use physical observables from both the meson and the baryon sectors, obtaining $m_{ud} = 3.636(66)(^{+60}_{-57})$~MeV and $m_s = 98.7(2.4)(^{+4.0}_{-3.2})$~MeV in the $\overline{\rm MS}(2\,{\rm GeV})$ scheme and $m_c = 1036(17)(^{+15}_{-8})$~MeV in the $\overline{\rm MS}(3\,{\rm GeV})$ scheme, where the first errors are statistical and the second ones are combinations of systematic errors. For the quark mass ratios we get $m_s / m_{ud} = 27.17(32)(^{+56}_{-38})$ and $m_c / m_s = 11.48(12)(^{+25}_{-19})$., Comment: 37 pages, 23 figures, 24 tables. One reference added

We present a determination of the ratio of kaon and pion leptonic decay constants in isosymmetric QCD (isoQCD), $f_K / f_\pi$, making use of the gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks, including configurations close to the physical point for all dynamical flavors. The simulations are carried out at three values of the lattice spacing ranging from $\sim 0.068$ to $\sim 0.092$ fm with linear lattice size up to $L \sim 5.5$~fm. The scale is set by the PDG value of the pion decay constant, $f_\pi^{isoQCD} = 130.4~(2)$ MeV, at the isoQCD pion point, $M_\pi^{isoQCD} = 135.0~(2)$ MeV, obtaining for the gradient-flow (GF) scales the values $w_0 = 0.17383~(63)$ fm, $\sqrt{t_0} = 0.14436~(61)$ fm and $t_0 / w_0 = 0.11969~(62)$ fm. The data are analyzed within the framework of SU(2) Chiral Perturbation Theory (ChPT) without resorting to the use of renormalized quark masses. At the isoQCD kaon point $M_K^{isoQCD} = 494.2~(4)$ MeV we get $(f_K / f_\pi)^{isoQCD} = 1.1995~(44)$, where the error includes both statistical and systematic uncertainties. Implications for the Cabibbo-Kobayashi-Maskawa (CKM) matrix element $|V_{us}|$ and for the first-row CKM unitarity are discussed., Comment: 68 pages, 14 figures, 12 tables. Version to appear in PRD

We review lattice results related to pion, kaon, $D$-meson, $B$-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $f_+(0)$ arising in the semileptonic $K \to \pi$ transition at zero momentum transfer, as well as the decay constant ratio $f_K/f_\pi$ and its consequences for the CKM matrix elements $V_{us}$ and $V_{ud}$. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $SU(2)_L\times SU(2)_R$ and $SU(3)_L\times SU(3)_R$ Chiral Perturbation Theory. We review the determination of the $B_K$ parameter of neutral kaon mixing as well as the additional four $B$ parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $m_c$ and $m_b$ 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. We review the status of lattice determinations of the strong coupling constant $\alpha_s$. Finally, in this review we have added a new section reviewing results for nucleon matrix elements of the axial, scalar and tensor bilinears, both isovector and flavor diagonal., Comment: 458 pages, 46 figures, 209 tables, 1146 references. Minor changes, version as published in EPJC. arXiv admin note: substantial text overlap with arXiv:1607.00299, arXiv:1310.8555

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)$., Comment: 23 pages, 9 figures, 7 tables

We present results on the light, strange and charm nucleon scalar and tensor charges from lattice QCD, using simulations with $N_f=2$ flavors of twisted mass Clover-improved fermions with a physical value of the pion mass. Both connected and disconnected contributions are included, enabling us to extract the isoscalar, strange and charm charges for the first time directly at the physical point. Furthermore, the renormalization is computed non-perturbatively for both isovector and isoscalar quantities. We investigate excited state effects by analyzing several sink-source time separations and by employing a set of methods to probe ground state dominance. Our final results for the scalar charges are $g_S^u = 5.20(42)(15)(12)$, $g_S^d = 4.27(26)(15)(12)$, $g_S^s=0.33(7)(1)(4)$, $g_S^c=0.062(13)(3)(5)$ and for the tensor charges $g_T^u = 0.782(16)(2)(13)$, $g_T^d = -0.219(10)(2)(13)$, $g_T^s=-0.00319(69)(2)(22)$, $g_T^c=-0.00263(269)(2)(37)$ in the $\overline{\rm MS}$ scheme at 2~GeV. The first error is statistical, the second is the systematic error due to the renormalization and the third the systematic arising from possible contamination due to the excited states., Comment: 20 pages and 13 figures

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., Comment: 383 pages, 32 figures, 160 tables, 795 references

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., Comment: 45 pages, 15 figures

We present results on the nucleon scalar, axial and tensor charges as well as on the momentum fraction, and the helicity and transversity moments. The pion momentum fraction is also presented. The computation of these key observables is carried out using lattice QCD simulations at a physical value of the pion mass. The evaluation is based on gauge configurations generated with two degenerate sea quarks of twisted mass fermions with a clover term. We investigate excited states contributions with the nucleon quantum numbers by analyzing three sink-source time separations. We find that, for the scalar charge, excited states contribute significantly and to a less degree to the nucleon momentum fraction and helicity moment. Our analysis yields a value for the nucleon axial charge agrees with the experimental value and we predict a value of 1.027(62) in the $\overline{\text{MS}}$ scheme at 2 GeV for the isovector nucleon tensor charge directly at the physical point. The pion momentum fraction is found to be $\langle x\rangle_{u-d}^{\pi^\pm}=0.214(15)(^{+12}_{-9})$ in the $\overline{\rm MS}$ at 2 GeV., Comment: The experimental value of the helicity plotted in Fig.17 is corrected. ETM Collaboration. 25 pages, 19 figures. Version accepted for PRD (volume 92, eid 114513)

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., Comment: 61 pages, 22 figures, 21 tables. Final version to appear in NPB

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., Comment: presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

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., Comment: 7 pages, 5 figures, talk presented at the XXVIII International Symposium on Lattice Field Theory, Lattice2010, June 14-19, 2010 Villasimius, Italy

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., Comment: 8 pages, 4 figures. Talk presented at the XXVII International Symposium on Lattice Field Theory, July 26-31 2009, Peking University, Beijing, China

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., Comment: 40 pages, 6 figures

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., Comment: 20 pages. v2: one author added, label in fig.3 corrected, reference changed

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., Comment: 19 pages, 18 figures, Version as accepted in Phys. Rev. D, Discussion on nucleon sigma-term extended, typos corrected

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., Comment: March 2008. 59pp. Minor changes. Added references

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., Comment: 7 pages, 4 figures, contribution to Lattice 2007, Regensburg, Germany, 30 July - 4 August 2007

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., Comment: 7 pages, 4 figures, Contribution to Lattice 2007, Regensburg, Germany, 30 July - 4 August 2007

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., Comment: 14 pages. Presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4 2007, Regensburg

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., Comment: 20 pages, 5 figures

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., Comment: 7 pages, 5 figures, Poster presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4 2007, Regensburg, Germany

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., Comment: 17 pages, 5 figures, references added, typos corrected

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., Comment: 44 pages, 10 figures

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., Comment: 53 pages, 11 tables, 18 figures; discussion extended

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., Comment: 8 pages. To appear in proceedings of Computational Hadron Physics, Nicosia, 14-17 Sep 2005. Notation made uniform with other publications

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., Comment: 6 pages, Talk presented at Lattice 2005 (machines and algorithms)

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., Comment: Contribution to the proceedings of Lattice 2005 resulting from three talks and one poster (Hadron spectrum and quark masses), 21 pages, 12 figures

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., Comment: 6 pages, poster (hadron spectrum and quark masses) at Lattice 2005,Dublin, July 25-30

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., Comment: 15 pages, 7 figures, typo corrected, web-list of authors corrected

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., Comment: 27 pages, 4 figures, revised version: typos corrected, reference added, minor other changes, accepted for publication in Comp. Phys. Commun

We extend our earlier investigation of the finite temperature deconfinement transition in SU(N) gauge theories, with the emphasis on what happens as N->oo. We calculate the latent heat in the continuum limit, and find the expected quadratic in N behaviour at large N. We confirm that the phase transition, which is second order for SU(2) and weakly first order for SU(3), becomes robustly first order for N>3 and strengthens as N increases. As an aside, we explain why the SU(2) specific heat shows no sign of any peak as T is varied across what is supposedly a second order phase transition. We calculate the effective string tension and electric gluon masses at T=Tc confirming the discontinuous nature of the transition for N>2. We explicitly show that the large-N `spatial' string tension does not vary with T for TTc it increases as T-squared to a good approximation, and the k-string tension ratios closely satisfy Casimir Scaling. Within very small errors, we find a single Tc at which all the k-strings deconfine, i.e. a step-by-step breaking of the relevant centre symmetry does not occur. We calculate the interface tension but are unable to distinguish between linear or quadratic in N variations, each of which can lead to a striking but different N=oo deconfinement scenario. We remark on the location of the bulk phase transition, which bounds the range of our large-N calculations on the strong coupling side, and within whose hysteresis some of our larger-N calculations are performed., Comment: 50 pages, 14 figures

We compare the behavior of overlap fermions, which are chirally invariant, and of Wilson twisted mass fermions at full twist in the approach to the chiral limit. Our quenched simulations reveal that with both formulations of lattice fermions pion masses of O(250 MeV) can be reached in practical applications. Our comparison is done at a fixed value of the lattice spacing a=0.123 fm. A number of quantities are measured such as hadron masses, pseudoscalar decay constants and quark masses obtained from Ward identities. We also determine the axial vector renormalization constants in the case of overlap fermions., Comment: 22 pages, 10 figures

We compare overlap fermions, which are chirally invariant, and Wilson twisted mass fermions in the approach to the chiral limit. Our quenched simulations reveal that with both formulations of lattice fermions pion masses of O(250 MeV) can be reached in practical simulations. Our comparison is done at a fixed lattice spacing a=0.123 fm. Several quantities are measured, such as hadron masses and pseudoscalar decay constants., Comment: Lattice2004(chiral)

We perform a comparative study of the infrared properties of overlap and staggered fermions in QCD. We observe that the infrared spectrum of the APE/HYP improved staggered Dirac operator develops a four-fold near-degeneracy and is in quantitative agreement with the infrared spectrum of the overlap operator. The near-degeneracy allows us to identify the zero modes of the staggered operator and we find that the number of zero modes is in line with the topological index of the overlap operator., Comment: Talk presented at Lattice2004(chiral), Fermilab, June 21-26, 2004, 3 pages, 2 figures

We present a systematic comparison of various iterative methods to obtain the fermion propagator with both overlap and twisted mass fermions at fixed pion mass in the quenched approximation. Taking the best available algorithm in each case we find that calculations with the overlap operator are by a factor of 20-40 more expensive than with the twisted mass operator at the parameter values considered here. For the overlap operator we also compare the efficiency of various methods for calculating the topological index., Comment: Poster presented at Lattice2004(machines), Fermilab, June 21-26, 2004, 3 pages, 4 figures

We review recent lattice results for the large $N$ limit of SU(N) gauge theories. In particular, we focus on glueball masses, topology and its relation to chiral symmetry breaking (relevant for phenomenology), on the tension of strings connecting sources in higher representations of the gauge group (relevant for models of confinement and as a comparative ground for theories beyond the Standard Model) and on the finite temperature deconfinement phase transition (relevant for RHIC-like experiments). In the final part we present open challenges for the future., Comment: 6 pages, 3 figures; summary of the talk given by B. Lucini and the poster presented by U. Wenger at the conference "Confinement 2003"

We investigate spontaneous chiral symmetry breaking in SU(N) gauge theories at large N using overlap fermions. The exact zero modes and the low-lying modes of the Dirac operator provide the tools to gain insight into the interplay between chiral symmetry breaking and topology. We find that topology indeed drives chiral symmetry breaking at N=3 as well as at large N. By comparing the results on various volumes and at different lattice spacings we are able to show that our conclusions are not affected by finite volume effects and also hold in the continuum limit. We then address the question whether the topology can be usefully described in terms of instantons., Comment: Talk at Lattice 2003 (chiral); 3 pages, 2 figures, espcrc2.sty

We study the deconfinement phase transition in SU(N) gauge theories for $N$=2,3,4,6,8. The transition is first order for $N \ge 3$, with the strength increasing as $N$ increases. We extrapolate $T_c/\sqrt{\sigma}$ to the continuum limit for each $N$, and observe a rapid approach to the large $N$ limit. As $N$ increases the phase transition becomes clear-cut on smaller spatial volumes, indicating the absence of (non-singular) finite volume corrections at $N=\infty$ -- reminiscent of large $N$ reduction. The observed rapid increase of the inter-phase surface tension with $N$ may indicate that for $N=\infty$ the deconfinement transition cannot, in practise, occur., Comment: Lattice2003(nonzero)

We address the question of whether the low modes of the Dirac operator are caused by topological objects such as instantons in SU(N) gauge theories. We study the pseudo-scalar density of these modes, finding the size distributions of the instantons, and comparing it with the underlying gauge field. We find that, although the near-zero modes of the Dirac operator depend on topology for all N, their small instanton content decreases as N increases., Comment: lattice2002(chiral)

We report on our ongoing investigation of the deconfining phase transition in SU(4) and SU(6) gauge theories. We calculate the critical couplings while taking care to avoid the influence of a nearby bulk phase transition. We determine the latent heat of the phase transition and investigate the order and the strength of the transition at large N. We also report on our determination of the critical temperature expressed in units of the string tension in the large N limit., Comment: Lattice 2002 (nonzerot), 3 pages, 2 figures

We investigate the properties of the deconfinement transition in SU(4) and SU(6) gauge theories. We find that it is a `normal' first order transition in both cases, from which we conclude that the transition is first order in the N->infinity limit. Comparing our preliminary estimates of the continuum values of Tc/sqrt(K) with existing values for SU(2) and SU(3) demonstrates a weak dependence on N for all values of N., Comment: 18 pages