Search

We present a lattice calculation of the Hadronic Vacuum Polarization (HVP) contribution of the strange and charm quarks to the anomalous magnetic moment of the muon in isospin symmetric QCD. We employ the gauge configurations generated by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks at five lattice spacings and at values of the quark mass parameters that are close and/or include the isospin symmetric QCD point of interest. After computing the small corrections necessary to precisely match this point, and carrying out an extrapolation to the continuum limit based on the data at lattice spacings $a \simeq 0.049, 0.057, 0.068, 0.080$~fm and spatial lattice sizes up to $L \simeq 7.6$~fm, we obtain $a_\mu^{\rm HVP}(s) = (53.57 \pm 0.63) \times 10^{-10}$ and $a_\mu^{\rm HVP}(c) = (14.56 \pm 0.13) \times 10^{-10}$, for the quark-connected strange and charm contributions, respectively. Our findings agree well with the corresponding results by other lattice groups., Comment: 33 pages, 18 figures, 3 tables

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

The Hamiltonian limit of lattice gauge theories can be found by extrapolating the results of anisotropic lattice computations, i.e., computations using lattice actions with different temporal and spatial lattice spacings ($a_t\neq a_s$), to the limit of $a_t\to 0$. In this work, we present a study of this Hamiltonian limit for a Euclidean $U(1)$ gauge theory in 2+1 dimensions (QED3), regularized on a toroidal lattice. The limit is found using the renormalized anisotropy $\xi_R=a_t/a_s$, by sending $\xi_R \to 0$ while keeping the spatial lattice spacing constant. We compute $\xi_R$ in $3$ different ways: using both the ``normal'' and the ``sideways'' static quark potential, as well as the gradient flow evolution of gauge fields. The latter approach will be particularly relevant for future investigations of combining quantum computations with classical Monte Carlo computations, which requires the matching of lattice results obtained in the Hamiltonian and Lagrangian formalisms.

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 present a calculation of the Mellin moments of the transverse quark spin densities in the nucleon using lattice QCD. The densities are extracted from the unpolarized and transversity generalized form factors extrapolated to the continuum limit using three $N_f=2+1+1$ twisted mass fermion gauge ensembles simulated with physical quark masses and spanning three lattice spacings. The first moment of transversely polarized quarks in an unpolarized nucleon shows an interesting distortion, which can be traced back to the sharp falloff of the transversity generalized form factor $\bar{B}_{Tn0}(t)$. The isovector tensor anomalous magnetic moment is determined to be $\kappa_T=1.051(94)$, which confirms a negative and large Boer-Mulders function, $h_1^{\perp}$, in the nucleon.

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 compute the nucleon axial and induced pseudoscalar form factors using three ensembles of gauge configurations, generated with dynamical light quarks with mass tuned to approximately their physical value. One of the ensembles also includes the strange and charm quarks with their mass close to physical. The latter ensemble has large statistics and finer lattice spacing and it is used to obtain final results, while the other two are used for assessing volume effects. The pseudoscalar form factor is also computed using these ensembles. We examine the momentum dependence of these form factors as well as relations based on pion pole dominance and the partially conserved axial-vector current hypothesis., Comment: 25 pages and 27 figures

We present results for the moments of nucleon isovector vector and axial generalised parton distribution functions computed within lattice QCD. Three ensembles of maximally twisted mass clover-improved fermions simulated with a physical value of the pion mass are analyzed. Two of these ensembles are generated using two degenerate light quarks. A third ensemble is used having, in addition to the light quarks, strange and charm quarks in the sea. A careful analysis of the convergence to the ground state is carried out that is shown to be essential for extracting the correct nucleon matrix elements. This allows a controlled determination of the unpolarised, helicity and tensor second Mellin moments. The vector and axial-vector generalised form factors are also computed as a function of the momentum transfer square up to about 1 GeV$^2$. The three ensembles allow us to check for unquenching effects and to assess lattice finite volume effects., Comment: Published version. Typos corrected, references and figures added. 22 pages, 17 figures

Using Nf=2+1+1 lattice QCD, we determine the fermionic connected contributions to the second and third moment of the pion PDF. Based on gauge configurations from the European Twisted Mass Collaboration, chiral and continuum extrapolations are performed using pion masses in the range of 230 to 500 MeV and three values of the lattice spacing. Finite volume effects are investigated using different volumes. In order to avoid mixing under renormalisation for the third moment, we use an operator with two non-zero spatial components of momentum. Momenta are injected using twisted boundary conditions. Our final values read $\langle x\rangle=0.2075(106)$ and $\langle x^2\rangle=0.163(33)$, determined at 2 GeV in the $\overline{MS}$-scheme and with systematic and statistical uncertainties summend in quadrature., Comment: 30 pages, matching the published version

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 for the isospin-0 $\pi\pi$ s-wave scattering length calculated in twisted mass lattice QCD. We use three $N_f = 2$ ensembles with unitary pion mass at its physical value, 240~MeV and 330~MeV respectively. We also use a large set of $N_f = 2 + 1 +1$ ensembles with unitary pion masses varying in the range of 230~MeV - 510~MeV at three different values of the lattice spacing. A mixed action approach with the Osterwalder-Seiler action in the valence sector is adopted to circumvent the complications arising from isospin symmetry breaking of the twisted mass quark action. Due to the relatively large lattice artefacts in the $N_f = 2 + 1 +1$ ensembles, we do not present the scattering lengths for these ensembles. Instead, taking the advantage of the many different pion masses of these ensembles, we qualitatively discuss the pion mass dependence of the scattering properties of this channel based on the results from the $N_f = 2 + 1 +1$ ensembles. The scattering length is computed for the $N_f = 2$ ensembles and the chiral extrapolation is performed. At the physical pion mass, our result $M_\pi a^\mathrm{I=0}_0 = 0.198(9)(6)$ agrees reasonably well with various experimental measurements and theoretical predictions., Comment: 7 pages, 2 figures, talk presented at the 34th annual International Symposium on Lattice Field Theory (Lattice2016), 24-30 July 2016, University of Southampton, UK

We present results for the isospin-0 $\pi\pi$ s-wave scattering length calculated with Osterwalder-Seiler valence quarks on Wilson twisted mass gauge configurations. We use three $N_f = 2$ ensembles with unitary (valence) pion mass at its physical value (250$\sim$MeV), at 240$\sim$MeV (320$\sim$MeV) and at 330$\sim$MeV (400$\sim$MeV), respectively. By using the stochastic Laplacian Heaviside quark smearing method, all quark propagation diagrams contributing to the isospin-0 $\pi\pi$ correlation function are computed with sufficient precision. The chiral extrapolation is performed to obtain the scattering length at the physical pion mass. Our result $M_\pi a^\mathrm{I=0}_0 = 0.198(9)(6)$ agrees reasonably well with various experimental measurements and theoretical predictions. Since we only use one lattice spacing, certain systematics uncertainties, especially those arising from unitary breaking, are not controlled in our result., Comment: 21 pages, 5 figures, 6 tables

We present first results on the $\rho$ resonance parameters obtained with $N_f=2+1+1$ Wilson twisted mass fermions at maximal twist. Using ensembles of the ETM collaboration, we provide results for two values of the lattice spacing and a range of pion mass values., Comment: 7 pages, 3 figures. Poster presented at the 33rd International Symposium on Lattice Field Theory, 14 -18 July 2015 (Kobe, Japan)

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)

We present results for the $I=2$ $\pi\pi$ scattering length using $N_f=2+1+1$ twisted mass lattice QCD for three values of the lattice spacing and a range of pion mass values. Due to the use of Laplacian Heaviside smearing our statistical errors are reduced compared to previous lattice studies. A detailed investigation of systematic effects such as discretisation effects, volume effects, and pollution of excited and thermal states is performed. After extrapolation to the physical point using chiral perturbation theory at NLO we obtain $M_\pi a_0=-0.0442(2)_\mathrm{stat}(^{+4}_{-0})_\mathrm{sys}$., Comment: Edited for typos, overhauled figures, more detailed comparison to existing lattice results

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

An overview is given of the lessons learned from the introduction of multi-threading using OpenMP in tmLQCD. In particular, programming style, performance measurements, cache misses, scaling, thread distribution for hybrid codes, race conditions, the overlapping of communication and computation and the measurement and reduction of certain overheads are discussed. Performance measurements and sampling profiles are given for different implementations of the hopping matrix computational kernel., Comment: presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

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 determine mass and mixing angles of eta and eta' states using Nf=2+1+1 Wilson twisted mass lattice QCD. We describe how those flavour singlet states need to be treated in this lattice formulation. Results are presented for three values of the lattice spacing, a=0.061 fm, a=0.078 fm and a=0.086 fm, with light quark masses corresponding to values of the charged pion mass in a range of 230 to 500 MeV and fixed bare strange and charm quark mass values. We obtain 557(15)(45) MeV for the eta mass (first error statistical, second systematic) and 44(5) degrees for the mixing angle in the quark flavour basis, corresponding to -10(5) degrees in the octet-singlet basis., Comment: 28 pages, 9 figures, version to appear in JHEP, extended discussion of autocorrelation times and comparison to results available in the literature, added a comment for FS-effects and clarified the description of our blocking procedure

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., Comment: 7 pages, 8 figures, Talk presented at the XXIX International Symposium on Lattice Field Theory (Lattice 2011), July 10-16, 2011, Squaw Valley, Lake Tahoe, California, USA

We investigate the thermal QCD transition with two flavors of maximally twisted mass fermions for a set of pion masses, 300 MeV \textless $m_\pi$ \textless 500 MeV, and lattice spacings $a$ \textless 0.09 fm. We determine the pseudo-critical temperatures and discuss their extrapolation to the chiral limit using scaling forms for different universality classes, as well as the scaling form for the magnetic equation of state. For all pion masses considered we find resonable consistency with O(4) scaling plus leading corrections. However, a true distinction between the O(4) scenario and a first order scenario in the chiral limit requires lighter pions than are currently in use in simulations of Wilson fermions., Comment: 11 pages, 11 figures

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 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., Comment: 8 pages, 4 figures, presented at International Workshop on Effective Field Theories: From the Pion to the Upsilon (EFT 09), Valencia, Spain, 2-6 Feb 2009

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

In this talk we give an overview of the 3D phase diagram of two-flavour non-zero temperature lattice QCD with twisted-mass Wilson fermions and a tree-level Symanzik-improved gauge action. We present a first feasibility study at maximal twist and, for the quenched case, we demonstrate automatic O(a)-improvement to work., Comment: 7 pages, 12 figures. To appear in the Proceedings of 27th International Symposium on Lattice Field Theory (Lattice 2009), Beijing, China, July 26-31, 2009

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 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., Comment: Proceedings of the 27th International Symposium on Lattice Field Theory (Lattice 2009), Beijing, China, 2009 July 26-31. 8 pages, 3 figures

We review various B meson decays that require knowledge of the transverse decay constant of the b1(1235) meson. We report on an exploratory lattice QCD calculation of the transverse decay constant of the b1 meson. The lattice QCD calculations used unquenched gauge configurations, at two lattice spacings, generated with two flavours of sea quarks. The twisted mass formalism is used., Comment: 11 pages, 3 figures. Fixed some typos, updated value of one of the sum rule results

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., Comment: 23 pages, 8 figures; version to appear in JHEP

In this paper we propose to interpret the large discretization artifacts affecting the neutral pion mass in maximally twisted lattice QCD simulations as O(a^2) effects whose magnitude is roughly proportional to the modulus square of the (continuum) matrix element of the pseudoscalar density operator between vacuum and one-pion state. The numerical size of this quantity is determined by the dynamical mechanism of spontaneous chiral symmetry breaking and turns out to be substantially larger than its natural magnitude set by the value of Lambda_QCD., Comment: 38 pages, 1 figure, 2 tables

We report results for the masses of the flavour non-singlet light 0++, 1--, and 1+- mesons from unquenched lattice QCD at two lattice spacings. The twisted mass formalism was used with two flavours of sea quarks. For the 0++ and 1+- mesons we look for the effect of decays on the mass dependence. For the light vector mesons we study the chiral extrapolations of the mass. We report results for the leptonic and transverse decay constants of the rho meson. We test the mass dependence of the KRSF relations., Comment: 35 pages, 14 figures. Added 2 references. Fixed small typos

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., Comment: accepted for publication in the 40th anniversary issue of CPC, 44 pages

We present results on the mass of the baryon octet and decuplet using two flavors of light dynamical twisted mass fermions. The strange quark mass is fixed to its physical value from the kaon sector in a partially quenched set up. Calculations are performed for light quark masses corresponding to a pion mass in the range 270-500 MeV and lattice sizes of 2.1 fm and 2.7 fm. We check for cut-off effects and isospin breaking by evaluating the baryon masses at two different lattice spacings. We carry out a chiral extrapolation for the octet baryons and discuss results for the Omega., Comment: 7 pages, 9 figures, Talk presented at The XXVI International Symposium on Lattice Field Theory, LATTICE 2008, July 14 - 19, 2008, Williamsburg, Virginia, USA

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 compare flavour singlet and non-singlet vector mesons from first principles using lattice QCD. With N_f=2 flavours of light quark, this addresses the omega-rho mass difference. Using maximally twisted-mass lattice QCD, we are able for the first time to determine this mass difference precisely and we compare with experiment. We also discuss omega-rho mixing effects arising within QCD through the u-d quark mass difference., Comment: 12 pages 4 figures

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., Comment: 7 pages, 2 figures, talk presented at the XXVI International Symposium on Lattice Field Theory, July 14-19, 2008, Williamsburg, Virginia, USA

We present results for chiral extrapolations of the mass and decay constants of the rho meson. The data sets used are the nf=2 unquenched gauge configurations generated with twisted mass fermions by the European Twisted Mass Collaboration. We describe a calculation of three decay constants in charmonium and explain why they are required., Comment: 7 pages, 2 figures. Talk given at the XXVI International Symposium on Lattice Field Theory, July 14 - 19 2008, Williamsburg, Virginia, USA

The behaviour of the one--dimensional random--forced Burgers equation is investigated in the path integral formalism, using a discrete space--time lattice. We show that by means of Monte Carlo methods one may evaluate observables, such as structure functions, as ensemble averages over different field realizations. The regularization of shock solutions to the zero--viscosity limit (Hopf-eq.) eventually leads to constraints on lattice parameters, required for the stability of the simulations. Insight into the formation of localized structures (shocks) and their dynamics is obtained., Comment: 4 pages, 3 figures submitted to EPL

We investigate the behaviour of stochastic differential equations, especially Burgers' eq., by means of Monte-Carlo-techniques., Comment: To be published in proceedings to "Path Integral - New Trends and Perspectives", Dresden 2007 6 pages, 2 figures

We study the flavour singlet pseudoscalar mesons from first principles using lattice QCD. With N_f=2 flavours of light quark, this is the so-called eta_2 meson and we discuss the phenomenological status of this. Using maximally twisted-mass lattice QCD, we extract the mass of the eta_2 meson at two values of the lattice spacing for lighter quarks than previously discussed in the literature. We are able to estimate the mass value in the limit of light quarks with their physical masses., Comment: 16 pages: version accepted for publication