Search

This document presents the initial scientific case for upgrading the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) to 22 GeV. It is the result of a community effort, incorporating insights from a series of workshops conducted between March 2022 and April 2023. With a track record of over 25 years in delivering the world's most intense and precise multi-GeV electron beams, CEBAF's potential for a higher energy upgrade presents a unique opportunity for an innovative nuclear physics program, which seamlessly integrates a rich historical background with a promising future. The proposed physics program encompass a diverse range of investigations centered around the nonperturbative dynamics inherent in hadron structure and the exploration of strongly interacting systems. It builds upon the exceptional capabilities of CEBAF in high-luminosity operations, the availability of existing or planned Hall equipment, and recent advancements in accelerator technology. The proposed program cover various scientific topics, including Hadron Spectroscopy, Partonic Structure and Spin, Hadronization and Transverse Momentum, Spatial Structure, Mechanical Properties, Form Factors and Emergent Hadron Mass, Hadron-Quark Transition, and Nuclear Dynamics at Extreme Conditions, as well as QCD Confinement and Fundamental Symmetries. Each topic highlights the key measurements achievable at a 22 GeV CEBAF accelerator. Furthermore, this document outlines the significant physics outcomes and unique aspects of these programs that distinguish them from other existing or planned facilities. In summary, this document provides an exciting rationale for the energy upgrade of CEBAF to 22 GeV, outlining the transformative scientific potential that lies within reach, and the remarkable opportunities it offers for advancing our understanding of hadron physics and related fundamental phenomena., Comment: Updates to the list of authors; Preprint number changed from theory to experiment; Updates to sections 4 and 6, including additional figures

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 show that a recent interesting idea to circumvent the difficulties with the continuation of parton distribution functions to the Euclidean region, that consists in looking at equal time correlators between proton states of infinite momentum, encounters some problems related to the power divergent mixing pattern of DIS operators, when implemented within the lattice regularization., Comment: 15 pages, no figures, Physical Review D (2017)

We present precise lattice computations for the b-quark mass, the quark mass ratios mb/mc and mb/ms as well as the leptonic B-decay constants. We employ gauge configurations with four dynamical quark flavors, up/down, strange and charm, at three values of the lattice spacing (a ~ 0.06 - 0.09 fm) and for pion masses as low as 210 MeV. Interpolation in the heavy quark mass to the bottom quark point is performed using ratios of physical quantities computed at nearby quark masses exploiting the fact that these ratios are exactly known in the static quark mass limit. Our results are also extrapolated to the physical pion mass and to the continuum limit and read: mb(MSbar, mb) = 4.26(10) GeV, mb/mc = 4.42(8), mb/ms = 51.4(1.4), fBs = 229(5) MeV, fB = 193(6) MeV, fBs/fB = 1.184(25) and (fBs/fB)/(fK/fpi) = 0.997(17)., Comment: Version to appear in PRD. Added comments to simulation setup and error budget discussion. 1+20 pages, 9 figures

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 unquenched lattice QCD results for the matrix elements of four-fermion operators relevant to the description of the neutral K and D mixing in the Standard Model and its extensions. We have employed simulations with Nf = 2 + 1 + 1 dynamical sea quarks at three values of the lattice spacings in the interval 0.06 - 0.09 fm and pseudoscalar meson masses in the range 210 - 450 MeV. Our results are extrapolated to the continuum limit and to the physical pion mass. Renormalization constants have been determined non-perturbatively in the RI-MOM scheme. In particular, for the Kaon bag-parameter, which is relevant for the \overline{K}^0-K^0 mixing in the Standard Model, we obtain B_K^{RGI} = 0.717(24)., Comment: Added comments to error budget discussion; fig.19 corrected. Version to appear in PRD

We present a lattice QCD calculation of the pseudoscalar decay constants fK, fD and fDs 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 values in the real world. The simulations are based on a unitary setup for the two light mass-degenerate quarks and on a mixed action approach for the strange and charm quarks. We use data simulated at three different values of the lattice spacing in the range 0.06 - 0.09 fm and at pion masses in the range 210 - 450 MeV. Our main results are: fK+ / fpi+ = 1.184 (16), fK+ = 154.4 (2.0) MeV, which incorporate the leading strong isospin breaking correction due to the up- and down-quark mass difference, and fK = 155.0 (1.9) MeV, fD = 207.4 (3.8) MeV, fDs = 247.2 (4.1) MeV, fDs / fD = 1.192 (22) and (fDs / fD) / (fK / fpi) = 1.003 (14) obtained in the isospin symmetric limit of QCD. Combined with the experimental measurements of the leptonic decay rates of kaon, pion, D- and Ds-mesons our results lead to the following determination of the CKM matrix elements: |Vus| = 0.2269 (29), |Vcd| = 0.2221 (67) and |Vcs| = 1.014 (24). Using the latest value of |Vud| from superallowed nuclear beta decays the unitarity of the first row of the CKM matrix is fulfilled at the permille level., Comment: 20 pp., 4 figures; revised version to appear in PRD; improved calculation of IB effects for fK+; minor changes in the final values. arXiv admin note: text overlap with arXiv:1403.4504

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 discuss a lattice QCD computation of the $B$-meson decay constants by the ETM collaboration where suitable ratios allow to reach the bottom quark sector by combining simulations around the charm-quark mass with an exactly known static limit. The different steps involved in this ratio method are discussed together with an account of the assessment of various systematic effects. A comparison of results from simulations with two and four flavour dynamical quarks is presented., Comment: 8 pages, 4 figures. Proceedings of the 8th International Workshop on the CKM Unitarity Triangle (CKM 2014), Vienna, Austria, September 8-12, 2014

We present precision lattice calculations of the pseudoscalar decay constants of the charmed sector as well as determinations of the bottom quark mass and its ratio to the charm quark mass. We employ Nf=2+1+1 dynamical quark gauge configurations generated by the European Twisted Mass Collaboration, using data at three values of the lattice spacing and pion masses as low as 210 MeV. Strange and charm sea quark masses are close to their physical values., Comment: 7 pages, Contribution to the Proceedings of the 37th International Conference on High Energy Physics (ICHEP2014), Valencia, Spain, July 2-9, 2014

We present the first unquenched lattice QCD results for the bag parameters controlling the short distance contribution to D meson oscillations in the Standard Model and beyond. We have used the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2 dynamical quarks, at four lattice spacings and light meson masses in the range 280-500 MeV. Renormalization is carried out non-perturbatively with the RI-MOM method. The bag-parameter results have been used to constrain New Physics effects in D-Dbar mixing, to put a lower bound to the generic New Physics scale and to constrain off-diagonal squark mass terms for TeV-scale Supersymmetry., Comment: 18 pages, 4 figures, 10 tables, updated to match the version accepted for publication in Physical Review D. This article draws form arXiv:1308.1851 which is on B-Bbar mixing instead of D-Dbar mixing

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

Taking inspiration from lattice QCD data, we argue that a finite non-perturbative contribution to the quark mass is generated as a consequence of the dynamical phenomenon of spontaneous chiral symmetry breaking, in turn triggered by the explicitly breaking of chiral symmetry induced by the critical Wilson term in the action. In pure lattice QCD this mass term cannot be separated from the unavoidably associated linearly divergent contribution. However, if QCD is enlarged to a theory where also a scalar field is present, coupled to an SU(2) doublet of fermions via a Yukawa and a Wilson-like term, then in the phase where the scalar field takes a non-vanishing expectation value, a dynamically generated and "naturally" light fermion mass (numerically unrelated to the expectation value of the scalar field) is conjectured to emerge at a critical value of the Yukawa coupling where the symmetry of the model is maximally enhanced. Masses dynamically generated in this way display a natural hierarchy according to which the stronger is the strongest of the interactions the fermion is subjected to the larger is its mass., Comment: Added more information in Fig 1. Added Fig. 10. Added an extra Appendix. Restructured a few sentences according to referee suggestions. Corrected a few misprints. All results unchanged. Now 50 pages and 10 Figures

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 a lattice QCD calculation of the pseudoscalar decay constants $f_K$, $f_D$ and $f_{D_s}$ performed by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical fermions. We simulated at three different values of the lattice spacing, the smallest being approximately $0.06fm$, and with pion masses as small as $210$MeV. Our main results are: $f_{K^+}/f_{\pi^+}=1.183(17)$, $f_{K^+}=154.4(2.1)$MeV, $f_{D_s}=242.1(8.3)$MeV, $f_D=201.9(8.0)$MeV, $f_{D_s}/f_D=1.199(25)$ and $(f_{D_s}/f_D) / (f_K/f_\pi) = 1.005(15)$., Comment: 7 pages, 5 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We present a lattice QCD determination of the $b$-quark mass and of the $f_{B_s}$ and $f_B$ decay constants performed with $N_f = 2 + 1 + 1$ twisted mass Wilson fermions. We have used simulations at three values of the lattice spacing generated by ETMC with pion masses ranging from 210 to 440 MeV. To obtain physical quantities we performed a combined chiral and continuum limit and an extrapolation in the heavy quark mass from the charm to the $b$-quark region using suitable ratios calculated at nearby heavy-quark masses having an exactly known static limit. Our results are: $m_b(m_b) = 4.29 (13)$ GeV, $f_B = 196 (9)$ MeV, $f_{B_s} = 235 (9)$ MeV, $f_{B_s} / f_B = 1.201 (25)$, $(f_{B_s}/f_B)/(f_K/f_\pi) = 1.007 (16)$ and $(f_{B_s}/f_B)/(f_{D_s}/f_D) = 1.008 (13)$., Comment: 7 pages, 7 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We present a lattice QCD determination of the average up-down, strange and charm quark masses based on simulations performed by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical fermions. We simulated at three different values of the lattice spacing, the smallest being approximately $0.06fm$, and with pion masses as small as $210 \text{MeV}$. Our results are: $m_{ud}(2\text{GeV})=3.70(17)\text{MeV}$, $m_s(2\text{GeV})=99.2(3.9)\text{MeV}$, $m_c(m_c)=1.350(49)\text{GeV}$, $m_s/m_{ud}=26.64(30)$ and $m_c/m_s=11.65(12)$.

We present an accurate lattice QCD computation of the b-quark mass, the B and Bs decay constants, the B-mixing bag-parameters for the full four-fermion operator basis, as well as estimates for \xi and f_{Bq}\sqrt{B_q} extrapolated to the continuum limit and the physical pion mass. We have used Nf = 2 dynamical quark gauge configurations at four values of the lattice spacing generated by ETMC. Extrapolation in the heavy quark mass from the charm to the bottom quark region has been carried out using ratios of physical quantities computed at nearby quark masses, having an exactly known infinite mass limit., Comment: 7 pages, 4 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We present a fully general derivation of the Laplace--Young formula and discuss the interplay between the intrinsic surface geometry and the extrinsic one ensuing from the immersion of the surface in the ordinary euclidean three-dimensional space. We prove that the (reversible) work done in a general surface deformation can be expressed in terms of the surface stress tensor and the variation of the intrinsic surface metric.

We show how it is possible to define and compute the potential between $q$ and $\bar q$ external sources in the singlet and octet (adjoint) representation of the colour group.

We present a lattice QCD computation of the b-quark mass, the B and B_s decay constants, the B-mixing bag parameters for the full four-fermion operator basis as well as determinations for \xi and f_{Bq}\sqrt{B_i^{(q)}} extrapolated to the continuum limit and to the physical pion mass. We used N_f = 2 twisted mass Wilson fermions at four values of the lattice spacing with pion masses ranging from 280 to 500 MeV. Extrapolation in the heavy quark mass from the charm to the bottom quark region has been carried out on ratios of physical quantities computed at nearby quark masses, exploiting the fact that they have an exactly known infinite mass limit. Our results are m_b(m_b, \overline{\rm{MS}})=4.29(12) GeV, f_{Bs}=228(8) MeV, f_{B}=189(8) MeV and f_{Bs}/f_B=1.206(24). Moreover with our results for the bag-parameters we find \xi=1.225(31), B_1^{(s)}/B_1^{(d)}=1.01(2), f_{Bd}\sqrt{\hat{B}_{1}^{(d)}} = 216(10) MeV and f_{Bs}\sqrt{\hat{B}_{1}^{(s)}} = 262(10) MeV. We also computed the bag parameters for the complete basis of the four-fermion operators which are required in beyond the SM theories. By using these results for the bag parameters we are able to provide a refined Unitarity Triangle analysis in the presence of New Physics, improving the bounds coming from B_{(s)}-\bar B_{(s)} mixing., Comment: version accepted for publication in JHEP, extended discussion on error budget, lattice results unchanged, updated upper bound for C1 in Tables 5 and 6, references added; 55 pages, 25 figures, 14 tables

We present a method to evaluate on the lattice the leading isospin breaking effects due to both the small mass difference between the up and down quarks and the QED interaction. Our proposal is applicable in principle to any QCD+QED gauge invariant hadronic observable which can be computed on the lattice. It is based on the expansion of the path-integral in powers of the small parameters (m_d - m_u)/Lambda_{QCD} and alpha_{em}, where m_f is the renormalized quark mass and alpha_{em} the renormalized fine structure constant. In this paper we discuss in detail the general strategy of the method and the conventional, although arbitrary, separation of QCD from QED isospin breaking corrections. We obtain results for the pion mass splitting, M_{pi+}^2-M_{pi0}^2= 1.44(13)(16) x 10^3 MeV^2, for the Dashen's theorem breaking parameter epsilon_{gamma}= 0.79(18)(18), for the light quark masses, [m_d - m_u](MSbar,2 GeV)= 2.39(8)(17) MeV, [m_u / m_d](MSbar,2 GeV)= 0.50(2)(3) and for the flavour symmetry breaking parameters R and Q. We also update our previous results for the QCD isospin breaking corrections to the Kl2 decay rate and for the QCD contribution to the neutron-proton mass splitting., Comment: 39 pages, 17 figures

We present a precise lattice QCD determination of the b-quark mass, of the B and Bs decay constants and first results for the B-meson bag parameters. For our computation we employ the so-called ratio method and our results benefit from the use of improved interpolating operators for the B-mesons. QCD calculations are performed with Nf = 2 dynamical light-quarks at four values of the lattice spacing and the results are extrapolated to the continuum limit. The preliminary results are mb(mb) = 4.35(12) GeV for the MSbar b-quark mass, fBs = 234(6) MeV and fB = 197(10) MeV for the B-meson decay constants, BBs(mb) = 0.90(5) and BB(mb) = 0.87(5) for the B-meson bag parameters., Comment: 6 pages, 3 figures. Proceedings of the 36th International Conference on High Energy Physics - ICHEP 2012; July 4-11 2012; Melbourne, Australia

We present a precise lattice QCD determination of the b-quark mass, of the B and Bs decay constants and first preliminary results for the B-mesons bag parameter. Simulations are performed with Nf = 2 Wilson twisted mass fermions at four values of the lattice spacing and the results are extrapolated to the continuum limit. Our calculation benefits from the use of improved interpolating operators for the B-mesons and employs the so-called ratio method. The latter allows a controlled interpolation at the b-quark mass between the relativistic data around and above the charm quark mass and the exactly known static limit., Comment: 7 pages, 4 figures, 1 table. Proceedings of the 30th International Symposium on Lattice Field Theory - Lattice 2012; June 24-29, 2012; Cairns, Australia

We present the ETMC results for the bag parameters describing the neutral kaon mixing in the Standard Model and beyond and preliminary results for the bag parameters controlling the short distance contributions in the D^0-\bar{D}^0 oscillations. We also present preliminary results for the B_{Bd}, B_{Bs}, B_{Bs}/B_{Bd} and \xi -parameter controlling B^0_-\bar{B}^0 oscillations in the Standard Model employing the so-called ratio method. Using Nf=2 maximally twisted sea quarks and Osterwalder-Seiler valence quarks we achieve both O(a)-improvement and continuum like renormalization pattern. Simulations are performed at three-values of the lattice spacing and several values of quark masses in the light, strange, charm region and above charm up to ~2.5m_c. Our results are extrapolated to the continuum limit and extrapolated/interpolated to the physical quark masses., Comment: 7 pages, 7 figures. Proceedings of Lattice 2012, the 30th International Symposium on Lattice Field Theory,June 24-29 2012, Cairns, Australia

We present the first unquenched, continuum limit, lattice QCD results for the matrix elements of the operators describing neutral kaon oscillations in extensions of the Standard Model. Owing to the accuracy of our calculation on \Delta S=2 weak Hamiltonian matrix elements, we are able to provide a refined Unitarity Triangle analysis improving the bounds coming from model independent constraints on New Physics. In our non-perturbative computation we use a combination of Nf=2 maximally twisted sea quarks and Osterwalder-Seiler valence quarks in order to achieve both O(a)-improvement and continuum-like renormalization properties for the relevant four-fermion operators. The calculation of the renormalization constants has been performed non-perturbatively in the RI-MOM scheme. Based on simulations at four values of the lattice spacing and a number of quark masses we have extrapolated/interpolated our results to the continuum limit and physical light/strange quark masses., Comment: small corrections to RC matrix Z_Q for i,j=2,3; marginal corrections to bag parameters B_4, B_5, and ratios R_4, R_5; results on NP scales from UTfit analysis unchanged; 66 pages, 12 figures

We present a new method to evaluate with high precision the isospin breaking effects due to the mass difference between the up and down quarks using lattice QCD. Our proposal is applicable in principle to any hadronic observable which can be computed on the lattice. It is based on the expansion of the path-integral in powers of the small parameter $m_d - m_u$. In this talk we discuss how to apply this method to compute the leading isospin breaking effects for several physical quantities of interest: the kaon masses, the kaon decay constants and the neutron-proton mass splitting., Comment: 7 pages, 3 figures. Talk presented at the XXIX International Symposium on Lattice Field Theory - Lattice 2011, July 10-16, 2011, Squaw Valley, Lake Tahoe, California. Typo corrections, accepted version

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 present preliminary results at {\beta} = 1.95 (a = 0.077 fm) on the first unquenched N_f=2+1+1 lattice computation of the B_K parameter which controls the neutral kaon oscillations in the Standard Model. Using N_f=2+1+1 maximally twisted sea quarks and Osterwalder-Seiler valence quarks we achieve O(a) improvement and a continuum-like renormalization pattern for the four-fermion operator. Our results are extrapolated/interpolated to the physical light/strange quark mass but not yet to the continuum limit. The computation of the relevant renormalization constants is performed non perturbatively in the RI'-MOM scheme using dedicated simulations with N_f=4 degenerate sea quark flavours produced by the ETM collaboration. We get B_K^{RGI} (a = 0.077) = 0.747(18), which when compared to our previous unquenched N_f=2 determination and most of the existing results, suggests a rather weak B_K^{RGI} dependence on the number of dynamical flavours. We are at the moment analysing lattice data at two additional {\beta} values which will allow us to perform an extrapolation to the continuum limit., Comment: 7 pages, 8 figures, Proceedings of Lattice 2011, XXIX International Symposium on Lattice Field Theory, Squaw Valley, Lake Tahoe, California

We present a new method to evaluate with high precision isospin breaking effects due to the small mass difference between the up and down quarks using lattice QCD. Our proposal is applicable in principle to any hadronic observable which can be computed on the lattice. It is based on the expansion of the path-integral in powers of the small parameter md-mu. In this paper, we apply this method to compute the leading isospin breaking effects for several physical quantities of interest: the kaon meson masses, the kaon decay constant, the form factors of semileptonic Kl3 decays and the neutron-proton mass splitting., Comment: 30 pages, 17 plots, version accepted for publication in JHEP, added few references and some comment, results unchanged

We present a lattice QCD determination of the b quark mass and of the B and B_s decay constants, performed with N_f=2 twisted mass Wilson fermions, by simulating at four values of the lattice spacing. In order to study the b quark on the lattice, two methods are adopted in the present work, respectively based on suitable ratios with exactly known static limit and on the interpolation between relativistic data, evaluated in the charm mass region, and the static point, obtained by simulating the HQET on the lattice. The two methods provide results in good agreement. For the b quark mass in the MSbar scheme and for the decay constants we obtain m_b(m_b)=4.29(14) GeV, f_B=195(12) MeV, f_Bs=232(10) MeV and f_Bs/f_B=1.19(5). As a byproduct of the analysis we also obtain the results for the f_D and f_Ds decay constants: f_D=212(8) MeV, f_Ds=248(6) MeV and f_Ds/f_D=1.17(5)., Comment: 23 pages, 10 figures, 2 tables. Added appendix showing the agreement of the data for the ratios with the HQE prediction. Matching JHEP published version

We report on an ongoing non-perturbative computation of RI-MOM scheme renormalization constants for the lattice action with four dynamical flavours currently in use by ETMC. For this goal dedicated simulations with four degenerate sea quark flavours are performed at several values of the standard and twisted quark mass parameters. We discuss a method for removing possible O(a) artifacts at all momenta and extrapolating renormalization constant estimators to the chiral limit. We give preliminary results at one lattice spacing., Comment: Presented at "The XXVIII International Symposium on Lattice Field Theory, Lattice2010", June 14-19, 2010, Villasimius, Italy; 7 pages, 8 figures

We present preliminary results on the of neutral kaon oscillations in extensions of the Standard Model. Using Nf=2 maximally twisted sea quarks and Osterwalder-Seiler valence quarks, we achieve both O(a)-improvement and continuum-like renormalization pattern for the relevant four-fermion operators. We perform simulations at three values of the lattice spacing and extrapolate/interpolate our results to the continuum limit and physical light/strange quark mass. The calculation of the renormalization constants of the complete operator basis is performed non- perturbatively in the RI-MOM scheme., Comment: 7 pages, 2 figures; presented at the XXVIII International Symposium on Lattice Field Theory Villasimius, Sardinia, Italy

We present the results of the recent high precision lattice calculation of the average up/down, strange and charm quark masses performed by ETMC with Nf=2 twisted mass Wilson fermions. The analysis includes data at four values of the lattice spacing and pion masses as low as ~270 MeV, allowing for accurate continuum limit and chiral extrapolation. The strange and charm masses are extracted by using several methods, based on different observables: the kaon and the eta_s meson for the strange quark and the D, D_s and eta_c mesons for the charm. The quark mass renormalization is carried out non-perturbatively using the RI-MOM method. The results for the quark masses in the MSbar scheme read: m_ud(2 GeV)= 3.6(2) MeV, m_s(2 GeV)=95(6) MeV and m_c(m_c)=1.28(4) GeV. We have also obtained the ratios m_s/m_ud=27.3(9) and m_c/m_s=12.0(3). Moreover, we provide the updated result for the bottom quark mass, m_b(m_b)=4.3(2) GeV, obtained using the method presented in 0909.3187 [hep-lat]., Comment: 7 pages, 7 figures, talk given at the XXVIII International Symposium on Lattice Field Theory (Lattice 2010), June 14-19 2010, Villasimius, Italy

We present an unquenched $N_f=2$ lattice computation of the $B_{K}$ parameter which controls $K^0-\bar K^0$ oscillations. A partially quenched setup is employed with two maximally twisted dynamical (sea) light Wilson quarks, and valence quarks of both the maximally twisted and the Osterwalder--Seiler variety. Suitable combinations of these two kinds of valence quarks lead to a lattice definition of the $B_{K}$ parameter which is both multiplicatively renormalizable and O($a$) improved. Employing the non-perturbative RI-MOM scheme, in the continuum limit and at the physical value of the pion mass we get $B^{\rm RGI}_K=0.729\pm 0.030$, a number well in line with the existing quenched and unquenched determinations., Comment: Accepted for publication in Phys. Rev. D; systematic error discussion slightly modified; total uncertainty estimate on the BK result unchanged; 42 pages, 25 figures

We present results for the renormalization constants of bilinear quark operators obtained by using the tree-level Symanzik improved gauge action and the Nf=2 twisted mass fermion action at maximal twist, which guarantees automatic O(a)-improvement. Our results are also relevant for the corresponding standard (un-twisted) Wilson fermionic action since the two actions only differ, in the massless limit, by a chiral rotation of the quark fields. The scale-independent renormalization constants ZV, ZA and the ratio ZP/ZS have been computed using the RI-MOM approach, as well as other alternative methods. For ZA and ZP/ZS, the latter are based on both standard twisted mass and Osterwalder-Seiler fermions, while for ZV a Ward Identity has been used. The quark field renormalization constant Zq and the scale dependent renormalization constants ZS, ZP and ZT are determined in the RI-MOM scheme. Leading discretization effects of O(g^2 a^2), evaluated in one-loop perturbation theory, are explicitly subtracted from the RI-MOM estimates., Comment: 31 pages, 11 figures

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 compute non-perturbatively the B-parameters of the complete basis of four-fermion operators needed to study the Kaon oscillations in the SM and in its supersymmetric extension. We perform numerical simulations with two dynamical maximally twisted sea quarks at three values of the lattice spacing on configurations generated by the ETMC. Unwanted operator mixings and O(a) discretization effects are removed by discretizing the valence quarks with a suitable Osterwalder-Seiler variant of the Twisted Mass action. Operators are renormalized non-perturbatively in the RI/MOM scheme. Our preliminary result for BK(RGI) is 0.73(3)(3)., Comment: 7 pages, 3 figures, 1 table, proceedings of the XXVII Int'l Symposyum on Lattice Field Theory (LAT2009), July 26-31 2009, Peking University, Beijing (China)

This work is the result of the ideas developed by Ken Yoshida about the possibility of extending the range of applications of the matrix model approach to the computation of the holomorphic superpotential of the beta-deformed N=4 super Yang-Mills theory both in the presence of a mass term and in the massless limit. Our formulae, while agreeing with all the existing results we can compare with, are valid also in the case of spontaneously broken gauge symmetry. We dedicate this paper to the memory of Ken, an unforgettable friend for all of us and a great scientist., Comment: 25 pages, uses JHEP3.cls; v2: minor corrections, references added; matches published version

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 present work in progress on the computation of the K-meson vector and tensor decay constants, as well as the B-parameter in Kaon oscillations. Our simulations are performed in a partially quenched setup, with two dynamical (sea) Wilson quark flavours, having a maximally twisted mass term. Valence quarks are either of the standard or the Osterwalder-Seiler maximally twisted variety. These two regularizations can be suitably combined in order to obtain a BK parameter which is both multiplicatively renormalizable and O(a) improved., Comment: 7 pages, 4 figures, talk presented at the XXVI International Symposium on Lattice Field Theory, July 14 - 19, 2008, Williamsburg, Virginia, USA

We show that in the massless N=1 supersymmetric Wess-Zumino theory it is possible to devise a computational strategy by which the x-space calculation of the ladder 4-point correlators can be carried out without introducing any regularization. As an application we derive a representation valid at all loop orders in terms of conformal invariant integrals. We obtain an explicit expression of the 3-loop ladder diagram for collinear external points., Comment: LaTeX, 17 pages, 8 figures

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