Search

Quark orbital angular momentum in the proton is evaluated via a Lattice QCD calculation of the second Mellin moment of the twist-3 generalized parton distribution $\widetilde{E}_{2T} $ in the forward limit. The connection between this approach to quark orbital angular momentum and approaches previously utilized in Lattice QCD calculations, via generalized transverse momentum-dependent parton distributions and via Ji's sum rule, is reviewed. This connection can be given in terms of Lorentz invariance and equation of motion relations. The calculation of the second Mellin moment of $\widetilde{E}_{2T} $ proceeds via a finite-momentum proton matrix element of a quark bilocal operator with a straight-line gauge connection and separation in both the longitudinal and transverse directions. The dependence on the former component serves to extract the second Mellin moment, whereas the dependence on the latter component provides a transverse momentum cutoff for the matrix element. Furthermore, a derivative of the matrix element with respect to momentum transfer in the forward limit is required, which is obtained using a direct derivative method. The calculation utilizes a clover fermion ensemble at pion mass 317 MeV. The resulting quark orbital angular momentum is consistent with previous evaluations through alternative approaches, albeit with greater statistical uncertainty using a comparable number of samples., Comment: 10 pages, 3 figures, to appear in the proceedings of the 25th International Spin Physics Symposium (SPIN 2023), PoS SPIN2023 (2024) 075

Previous Lattice QCD calculations of nucleon transverse momentum-dependent parton distributions (TMDs) focused on the case of transversely polarized nucleons, and thus did not encompass two leading-twist TMDs associated with longitudinal polarization, namely, the helicity TMD and the worm-gear TMD corresponding to transversely polarized quarks in a longitudinally polarized nucleon. Based on a definition of TMDs via hadronic matrix elements of quark bilocal operators containing staple-shaped gauge connections, TMD observables characterizing the aforementioned two TMDs are evaluated, utilizing an RBC/UKQCD domain wall fermion ensemble at the physical pion mass., Comment: 9 pages, 4 figures, to appear in the proceedings of the 39th International Symposium on Lattice Field Theory, LATTICE2022, PoS LATTICE2022, 103

Generalized transverse momentum-dependent parton distributions (GTMDs) provide a comprehensive framework for imaging the internal structure of the proton. In particular, by encoding the simultaneous distribution of quark transverse positions and momenta, they allow one to directly access longitudinal quark orbital angular momentum, and, moreover, to correlate it with the quark helicity. The relevant GTMD is evaluated through a lattice calculation of a proton matrix element of a quark bilocal operator (the separation in which is Fourier conjugate to the quark momentum) featuring a momentum transfer (which is Fourier conjugate to the quark position), as well as the Dirac structure appropriate for capturing the quark helicity. The weighting by quark transverse position requires a derivative with respect to momentum transfer, which is obtained in unbiased fashion using a direct derivative method. The lattice calculation is performed directly at the physical pion mass, using domain wall fermions to mitigate operator mixing effects. Both the Jaffe-Manohar as well as the Ji quark spin-orbit correlations are extracted, yielding evidence for a strong quark spin-orbit coupling in the proton., Comment: 8 pages, 2 figures, to appear in the proceedings of the 38th International Symposium on Lattice Field Theory, LATTICE2021, PoS LATTICE2021, 413

A Lattice QCD approach to quark orbital angular momentum in the proton based on generalized transverse momentum-dependent parton distributions (GTMDs) is enhanced methodologically by incorporating a direct derivative technique. This improvement removes a significant numerical bias that had been seen to afflict results of a previous study. In particular, the value obtained for Ji quark orbital angular momentum is reconciled with the one obtained independently via Ji's sum rule, validating the GMTD approach. Since GTMDs simultaneously contain information about the quark impact parameter and transverse momentum, they permit a direct evaluation of the cross product of the latter. They are defined through proton matrix elements of a quark bilocal operator containing a Wilson line; the choice in Wilson line path allows one to continuously interpolate from Ji to Jaffe-Manohar quark orbital angular momentum. The latter is seen to be significantly enhanced in magnitude compared to Ji quark orbital angular momentum, confirming previous results., Comment: 12 pages, 6 figures

Quark orbital angular momentum (OAM) in the proton can be calculated directly given a Wigner function encoding the simultaneous distribution of quark transverse positions and momenta. This distribution can be accessed via proton matrix elements of a quark bilocal operator (the separation in which is Fourier conjugate to the quark momentum) featuring a momentum transfer (which is Fourier conjugate to the quark position). To generate the weighting by quark transverse position needed to calculate OAM, a derivative with respect to momentum transfer is consequently required. This derivative is evaluated using a direct derivative method, i.e., a method in which the momentum derivative of a correlator is directly sampled in the lattice calculation, as opposed to extracting it a posteriori from the numerical correlator data. The method removes the bias stemming from estimating the derivative a posteriori that was seen to afflict a previous exploratory calculation. Data for Ji OAM generated on a clover ensemble at pion mass $m_{\pi } = 317\, \mbox{MeV} $ are seen to agree with the result obtained via the traditional Ji sum rule method. By varying the gauge connection in the quark bilocal operator, also Jaffe-Manohar OAM is extracted, and seen to be enhanced significantly compared to Ji OAM., Comment: 7 pages, 3 figures, to appear in the proceedings of the 23rd International Spin Physics Symposium (SPIN2018), 10-14 September 2018, Ferrara, Italy, and in the proceedings of the 36th Annual International Symposium on Lattice Field Theory (LATTICE2018), 22-28 July 2018, Michigan State University, East Lansing, Michigan, USA

We compute the leading, strong-interaction contribution to the anomalous magnetic moment of the electron, muon and tau using lattice quantum chromodynamics (QCD) simulations. Calculations include the effects of $u$, $d$, $s$ and $c$ quarks and are performed directly at the physical values of the quark masses and in volumes of linear extent larger than $6\,\mathrm{fm}$. All connected and disconnected Wick contractions are calculated. Continuum limits are carried out using six lattice spacings. We obtain $a_e^\mathrm{LO-HVP}=189.3(2.6)(5.6)\times 10^{-14}$, $a_\mu^\mathrm{LO-HVP}=711.1(7.5)(17.4)\times 10^{-10}$ and $a_\tau^\mathrm{LO-HVP}=341.0(0.8)(3.2)\times 10^{-8}$, where the first error is statistical and the second is systematic., Comment: 17 pages, 8 figures (in 13 PDF files), RevTeX 4.1. Minor changes to results and to text. References updated. Matches version published in Physical Review Letters

We present a determination of the corrections to Dashen's theorem and of the individual up and down quark masses from a lattice calculation based on quenched QED and $N_f=2+1$ QCD simulations with 5 lattice spacings down to 0.054 fm. The simulations feature lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. For the parameter which quantifies violations to Dashens's theorem we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4) \, MeV$ and $m_d=4.67(6)(5)(4) \, MeV$ in the $\overline{MS}$ scheme at $2 \, GeV$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than 24 standard deviations., Comment: 7 pages, 4 figures, Proceedings of the 34th annual International Symposium on Lattice Field Theory, 24-30 July 2016, University of Southampton, UK

We compute the slope and curvature, at vanishing four-momentum transfer squared, of the leading order hadron vacuum polarization function, using lattice QCD. Calculations are performed with 2+1+1 flavors of staggered fermions directly at the physical values of the quark masses and in volumes of linear extent larger than 6fm. The continuum limit is carried out using six different lattice spacings. All connected and disconnected contributions are calculated, up to and including those of the charm., Comment: 6 pages, 3 figures

In a previous letter (arXiv:1306.2287) we determined the isospin mass splittings of the baryon octet from a lattice calculation based on quenched QED and $N_f{=}2{+}1$ QCD simulations with 5 lattice spacings down to $0.054~\mathrm{fm}$, lattice sizes up to $6~\mathrm{fm}$ and average up-down quark masses all the way down to their physical value. Using the same data we determine here the corrections to Dashen's theorem and the individual up and down quark masses. For the parameter which quantifies violations to Dashens's theorem, we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4)~\mathrm{MeV}$ and $m_d=4.67(6)(5)(4)~\mathrm{MeV}$ in the $\bar{\mathrm{MS}}$ scheme at $2~\mathrm{GeV}$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than $24$ standard deviations.

We present a QCD calculation of the $u$, $d$ and $s$ scalar quark contents of nucleons based on $47$ lattice ensembles with $N_f = 2+1$ dynamical sea quarks, $5$ lattice spacings down to $0.054\,\text{fm}$, lattice sizes up to $6\,\text{fm}$ and pion masses down to $120\,\text{MeV}$. Using the Feynman-Hellmann theorem, we obtain $f^N_{ud} = 0.0405(40)(35)$ and $f^N_s = 0.113(45)(40)$, which translates into $\sigma_{\pi N}=38(3)(3)\,\text{MeV}$, $\sigma_{sN}=105(41)(37)\,\text{MeV}$ and $y_N=0.20(8)(8)$ for the sigma terms and the related ratio, where the first errors are statistical and the second are systematic. Using isospin relations, we also compute the individual up and down quark contents of the proton and neutron (results in the main text)., Comment: 6 pages, 2 figures

We report initial nucleon structure results computed on lattices with 2+1 dynamical M\"obius domain wall fermions at the physical point generated by the RBC and UKQCD collaborations. At this stage, we evaluate only connected quark contributions. In particular, we discuss the nucleon vector and axial-vector form factors, nucleon axial charge and the isovector quark momentum fraction. From currently available statistics, we estimate the stochastic accuracy of the determination of $g_A$ and $_{u-d}$ to be around 10%, and we expect to reduce that to 5% within the next year. To reduce the computational cost of our calculations, we extensively use acceleration techniques such as low-eigenmode deflation and all-mode-averaging (AMA). We present a method for choosing optimal AMA parameters., Comment: 7 pages, 11 figures; talk presented at the 32nd International Symposium on Lattice Field Theory, 23-28 June, 2014, Columbia University, New York, USA

The existence and stability of atoms rely on the fact that neutrons are more massive than protons. The measured mass difference is only 0.14\% of the average of the two masses. A slightly smaller or larger value would have led to a dramatically different universe. Here, we show that this difference results from the competition between electromagnetic and mass isospin breaking effects. We performed lattice quantum-chromodynamics and quantum-electrodynamics computations with four nondegenerate Wilson fermion flavors and computed the neutron-proton mass-splitting with an accuracy of $300$ kilo-electron volts, which is greater than $0$ by $5$ standard deviations. We also determine the splittings in the $\Sigma$, $\Xi$, $D$ and $\Xi_{cc}$ isospin multiplets, exceeding in some cases the precision of experimental measurements., Comment: 57 pages, 15 figures, 6 tables, revised version

We present lattice QCD calculations of nucleon electromagnetic form factors using pion masses $m_\pi$ = 149, 202, and 254 MeV and an action with clover-improved Wilson quarks coupled to smeared gauge fields, as used by the Budapest-Marseille-Wuppertal collaboration. Particular attention is given to removal of the effects of excited state contamination by calculation at three source-sink separations and use of the summation and generalized pencil-of-function methods. The combination of calculation at the nearly physical mass $m_\pi$ = 149 MeV in a large spatial volume ($m_\pi L_s$ = 4.2) and removal of excited state effects yields agreement with experiment for the electric and magnetic form factors $G_E(Q^2)$ and $G_M(Q^2)$ up to $Q^2$ = 0.5 GeV$^2$., Comment: v2: published version; 30 pages, 25 figures, 6 tables

Recent results for moments of multiplicity distributions of net-protons and net-electric charge from the STAR collaboration are compared to lattice QCD results for higher order fluctuations of baryon number and electric charge by the Wuppertal-Budapest collaboration, with the purpose of extracting the freeze-out temperature and chemical potential. All lattice simulations are performed for a system of 2+1 dynamical quark flavors, at the physical mass for light and strange quarks; all results are continuum extrapolated. We show that it is possible to extract an upper value for the freeze-out temperature, as well as precise baryo-chemical potential values corresponding to the four highest collision energies of the experimental beam energy scan. Consistency between the freeze-out parameters obtained from baryon number and electric charge fluctuations is found. The freeze-out chemical potentials are now in agreement with the statistical hadronization model., Comment: 5 pages, 4 figures, references added, discussion added to the introduction, results unchanged

We present continuum extrapolated lattice results for the higher order fluctuations of conserved charges in high temperature Quantum Chromodynamics. Through the matching of the grand canonical ensemble on the lattice to the net charge and net baryon distribution realized in heavy ion experiments the temperature and the chemical potential may be estimated at the time of chemical freeze-out, Comment: 7 pages, 5 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

In lattice QCD computations a substantial amount of work is spent in solving linear systems arising in Wilson's discretization of the Dirac equations. We show first numerical results of the extension of the two-level DD-\alpha AMG method to a true multilevel method based on our parallel MPI-C implementation. Using additional levels pays off, allowing to cut down the core minutes spent on one system solve by a factor of approximately 700 compared to standard Krylov subspace methods and yielding another speed-up of a factor of 1.7 over the two-level approach., Comment: 12 pages, 5 figures

While electromagnetic and up-down quark mass difference effects on octet baryon masses are very small, they have important consequences. The stability of the hydrogen atom against beta decay is a prominent example. Here we include these effects by adding them to valence quarks in a lattice QCD calculation based on $N_f=2+1$ simulations with 5 lattice spacings down to 0.054 fm, lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. This allows us to gain control over all systematic errors, except for the one associated with neglecting electromagnetism in the sea. We compute the octet baryon isomultiplet mass splittings, as well as the individual contributions from electromagnetism and the up-down quark mass difference. Our results for the total splittings are in good agreement with experiment., Comment: 5 pages, 3 figures. Version accepted for publication by Phys. Rev. Lett

We present our results for ratios of higher order fluctuations of electric charge as functions of the temperature. These results are obtained in a system of 2+1 quark flavors at physical quark masses and continuum extrapolated. We compare them to preliminary data on higher order moments of the net electric charge distribution from the STAR collaboration. This allows us to determine the freeze-out temperature and chemical potential from first principles. We also show continuum-extrapolated results for ratios of higher order fluctuations of baryon number. These will allow to test the consistency of the approach, by comparing them to the corresponding experimental data (once they become available) and thus extracting the freeze-out parameters in an independent way., Comment: 5 pages, 7 figures, revtex

We consider the effect of strong external electromagnetic fields on thermodynamic observables in QCD, through lattice simulations with 1+1+1 flavors of staggered quarks at physical quark masses. Continuum extrapolated results are presented for the light quark condensates and for their tensor polarizations, as functions of the temperature and the magnetic field. We find the light condensates to undergo inverse magnetic catalysis in the transition region, in a manner that the transition temperature decreases with growing magnetic field. We also compare the results to other approaches and lattice simulations. Furthermore, we relate the tensor polarization to the spin part of the magnetic susceptibility of the QCD vacuum, and show that this contribution is diamagnetic., Comment: 13 pages, 9 figures, talks presented by FB and GE at Xth Quark Confinement and the Hadron Spectrum, 8-12 October 2012, TUM Campus Garching, Munich, Germany

We report the first Lattice QCD calculation using the almost physical pion mass mpi=149 MeV that agrees with experiment for four fundamental isovector observables characterizing the gross structure of the nucleon: the Dirac and Pauli radii, the magnetic moment, and the quark momentum fraction. The key to this success is the combination of using a nearly physical pion mass and excluding the contributions of excited states. An analogous calculation of the nucleon axial charge governing beta decay has inconsistencies indicating a source of bias at low pion masses not present for the other observables and yields a result that disagrees with experiment., Comment: journal version; 15 pages, 6 figures

We present 2+1 flavor Lattice QCD calculations of the nucleon scalar and tensor charges. Using the BMW clover-improved Wilson action with pion masses between 150 and 350 MeV and three source-sink separations between 0.9 and 1.4 fm, we achieve good control over excited-state contamination and extrapolation to the physical pion mass. As a consistency check, we also present results from calculations using unitary domain wall fermions with pion masses between 300 and 400 MeV, and using domain wall valence quarks and staggered sea quarks with pion masses between 300 and 600 MeV., Comment: 10 pages, 6 figures

We use the Wilson flow to define the gauge anisotropy at a given physical scale. We demonstrate the use of the anisotropic flow by performing the tuning of the bare gauge anisotropy in the tree-level Symanzik action for several lattice spacings and target anisotropies. We use this method to tune the anisotropy parameters in full QCD, where we also exploit the diminishing effect of a well chosen smearing on the renormalization of the fermion anisotropy., Comment: 21 pages, 7 figures, 2 tables

QCD thermodynamics is considered using Wilson fermions in the fixed scale approach. The temperature dependence of the renormalized chiral condensate, quark number susceptibility and Polyakov loop is measured at four lattice spacings allowing for a controlled continuum limit. The light quark masses are fixed to heavier than physical values in this first study. Finite volume effects are ensured to be negligible by using approriately large box sizes. The final continuum results are compared with staggered fermion simulations performed in the fixed N_t approach. The same continuum renormalization conditions are used in both approaches and the final results agree perfectly., Comment: 15 pages, 9 figures

We determine the equation of state of QCD for nonzero chemical potentials via a Taylor expansion of the pressure. The results are obtained for N_f=2+1 flavors of quarks with physical masses, on various lattice spacings. We present results for the pressure, interaction measure, energy density, entropy density, and the speed of sound for small chemical potentials. At low temperatures we compare our results with the Hadron Resonance Gas model. We also express our observables along trajectories of constant entropy over particle number. A simple parameterization is given (the Matlab/Octave script parameterization.m, submitted to the arXiv along with the paper), which can be used to reconstruct the observables as functions of T and mu, or as functions of T and S/N., Comment: 14 pages, 15 figures, version accepted for publication in JHEP

We study QCD thermodynamics using two flavors of dynamical overlap fermions with quark masses corresponding to a pion mass of 350 MeV. We determine several observables on N_t=6 and 8 lattices. All our runs are performed with fixed global topology. Our results are compared with staggered ones and a nice agreement is found., Comment: 14 pages, 6 figures, 1 table

Scale setting is of central importance in lattice QCD. It is required to predict dimensional quantities in physical units. Moreover, it determines the relative lattice spacings of computations performed at different values of the bare coupling, and this is needed for extrapolating results into the continuum. Thus, we calculate a new quantity, $w_0$, for setting the scale in lattice QCD, which is based on the Wilson flow like the scale $t_0$ (M. Luscher, JHEP 1008 (2010) 071). It is cheap and straightforward to implement and compute. In particular, it does not involve the delicate fitting of correlation functions at asymptotic times. It typically can be determined on the few per-mil level. We compute its continuum extrapolated value in 2+1-flavor QCD for physical and non-physical pion and kaon masses, to allow for mass-independent scale setting even away from the physical mass point. We demonstrate its robustness by computing it with two very different actions (one of them with staggered, the other with Wilson fermions) and by showing that the results agree for physical quark masses in the continuum limit., Comment: 15 pages, 7 figures, 2 tables; Version published in JHEP

The effect of an external magnetic field on the finite temperature transition of QCD is studied. We measure thermodynamic observables including the quark condensates and susceptibilities and the strange quark number susceptibility. We generate configurations at various values of the quantized magnetic flux with N_f=2+1 flavors of stout smeared staggered quarks at physical quark masses. We perform the renormalization of our observables and approach the continuum limit with N_t=6,8 and 10 lattices. We also check finite volume effects using various lattice volumes. Our main result is that the transition temperature significantly decreases with growing magnetic field, and that the transition remains an analytic crossover up to our largest external field eB ~ 1 GeV^2., 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, USA

The effect of an external (electro)magnetic field on the finite temperature transition of QCD is studied. We generate configurations at various values of the quantized magnetic flux with $N_f=2+1$ flavors of stout smeared staggered quarks, with physical masses. Thermodynamic observables including the chiral condensate and susceptibility, and the strange quark number susceptibility are measured as functions of the field strength. We perform the renormalization of the studied observables and extrapolate the results to the continuum limit using $N_t=6,8$ and 10 lattices. We also check for finite volume effects using various lattice volumes. We find from all of our observables that the transition temperature $T_c$ significantly decreases with increasing magnetic field. This is in conflict with various model calculations that predict an increasing $T_c(B)$. From a finite volume scaling analysis we find that the analytic crossover that is present at B=0 persists up to our largest magnetic fields $eB \approx 1 \textmd{GeV}^2$, and that the transition strength increases mildly up to this $eB\approx1 \textmd{GeV}^2$., Comment: 22 pages, 13 figures

By using lattice QCD computations we determine the sigma terms and strangeness content of all octet baryons by means of an application of the Hellmann-Feynman theorem. In addition to polynomial and rational expressions for the quark mass dependence of octet members, we use SU(3) covariant baryon chiral perturbation theory to perform the extrapolation to the physical up and down quark masses. Our N_f=2+1 lattice ensembles include pion masses down to about 190 MeV in large volumes (M_\pi L > 4), and three values of the lattice spacing. Our main results are the nucleon sigma term \sigma_{\pi N} = 39(4)(^{+18}_{-7}) and the strangeness content y_{N} = 0.20(7)(^{+13}_{-17}). Under the assumption of validity of covariant baryon \chi PT in our range of masses one finds y_{N} = 0.276(77)(^{+90}_{-62})., Comment: LaTeX. 15 pages, 5 figures. Text and results corrected according to Erratum [Phys. Rev. D 93, 039905(E) (2016)]

We study the spectra of heavy-light and heavy-heavy mesons containing charm quarks, including higher spin states. We use two sets of $N_f = 2 + 1$ gauge configurations, one set from QCDSF using the SLiNC action, and the other configurations from the Budapest-Marseille-Wuppertal collaboration, using the HEX smeared clover action. To extract information about the excited states, we choose a suitable basis of operators to implement the variational method., Comment: 7 pages, 5 figures, Talk presented at the XXIX International Symposium on Lattice Field Theory, Lattice2011, July 11-16, 2011, The Village at Squaw Valley, California, USA

Indirect CP violation in K \rightarrow {\pi}{\pi} decays plays a central role in constraining the flavor structure of the Standard Model (SM) and in the search for new physics. For many years the leading uncertainty in the SM prediction of this phenomenon was the one associated with the nonperturbative strong interaction dynamics in this process. Here we present a fully controlled lattice QCD calculation of these effects, which are described by the neutral kaon mixing parameter B_K . We use a two step HEX smeared clover-improved Wilson action, with four lattice spacings from a\approx0.054 fm to a\approx0.093 fm and pion masses at and even below the physical value. Nonperturbative renormalization is performed in the RI-MOM scheme, where we find that operator mixing induced by chiral symmetry breaking is very small. Using fully nonperturbative continuum running, we obtain our main result B_K^{RI}(3.5GeV)=0.531(6)_{stat}(2)_{sys}. A perturbative 2-loop conversion yields B_K^{MSbar-NDR}(2GeV)=0.564(6)_{stat}(3)_{sys}(6)_{PT}, which is in good agreement with current results from fits to experimental data., Comment: 10 pages, 7 figures. v2: Added one reference and one figure, replaced 2 figures for better readability and updated ensembles, conclusions unchanged. Final, published version

Based on a series of lattice calculations we determine the ratio FK/Fpi in QCD. With experimental data from kaon decay and nuclear double beta decay, we obtain a precise determination of |Vus|. Our simulation includes 2+1 flavours of sea quarks, with three lattice spacings, large volumes and a simulated pion mass reaching down to about 190 MeV for a full control over the systematic uncertainties., Comment: 6 pages, 2 figures. Proceedings of CKM2010, the 6th International Workshop on the CKM Unitarity Triangle, University of Warwick, UK, 6-10 September 2010

A status report is given for a joint project of the Budapest-Marseille-Wuppertal collaboration and the Regensburg group to study the quark mass-dependence of octet baryons in SU(3) Baryon XPT. This formulation is expected to extend to larger masses than Heavy-Baryon XPT. Its applicability is tested with 2+1 flavor data which cover three lattice spacings and pion masses down to about 190 MeV, in large volumes. Also polynomial and rational interpolations in M_\pi^2 and M_K^2 are used to assess the uncertainty due to the ansatz. Both frameworks are combined to explore the precision to be expected in a controlled determination of the nucleon sigma term and strangeness content., Comment: Lattice 2010

At the precision reached in current lattice QCD calculations, electromagnetic effects are becoming numerically relevant. We will present preliminary results for electromagnetic corrections to light hadron masses, based on simulations in which a $\mathrm{U}(1)$ degree of freedom is superimposed on $N_f=2+1$ QCD configurations from the BMW collaboration., Comment: 7 pages, 2 figures, The XXVIII International Symposium on Lattice Field Theory, June 14-19,2010, Villasimius, Sardinia Italy

While the masses of light hadrons have been extensively studied in lattice QCD simulations, there exist only a few exploratory calculations of the strong decay widths of hadronic resonances. We will present preliminary results of a computation of the rho meson width obtained using $N_f=2+1$ flavor simulations. The work is based on L\"uscher's formalism and its extension to moving frames., Comment: The XXVIII International Symposium on Lattice Field Theory, June 14-19,2010, Villasimius, Sardinia Italy

We give details of our precise determination of the light quark masses m_{ud}=(m_u+m_d)/2 and m_s in 2+1 flavor QCD, with simulated pion masses down to 120 MeV, at five lattice spacings, and in large volumes. The details concern the action and algorithm employed, the HMC force with HEX smeared clover fermions, the choice of the scale setting procedure and of the input masses. After an overview of the simulation parameters, extensive checks of algorithmic stability, autocorrelation and (practical) ergodicity are reported. To corroborate the good scaling properties of our action, explicit tests of the scaling of hadron masses in N_f=3 QCD are carried out. Details of how we control finite volume effects through dedicated finite volume scaling runs are reported. To check consistency with SU(2) Chiral Perturbation Theory the behavior of M_\pi^2/m_{ud} and F_\pi as a function of m_{ud} is investigated. Details of how we use the RI/MOM procedure with a separate continuum limit of the running of the scalar density R_S(\mu,\mu') are given. This procedure is shown to reproduce the known value of r_0m_s in quenched QCD. Input from dispersion theory is used to split our value of m_{ud} into separate values of m_u and m_d. Finally, our procedure to quantify both systematic and statistical uncertainties is discussed., Comment: 45 pages

Ordinary matter is described by six fundamental parameters: three couplings (gravitational, electromagnetic and strong) and three masses: the electron's (m_e) and those of the up (m_u) and down (m_d) quarks. An additional mass enters through quantum fluctuations: the strange quark mass (m_s). The three couplings and m_e are known with an accuracy of better than a few per mil. Despite their importance, $m_u$, $m_d$ (their average m_{ud}) and m_s are relatively poorly known: e.g. the Particle Data Group quotes them with conservative errors close to 25%. Here we determine these quantities with a precision below 2% by performing ab initio lattice quantum chromodynamics (QCD) calculations, in which all systematics are controlled. We use pion and quark masses down to (and even below) their physical values, lattice sizes of up to 6 fm, and five lattice spacings to extrapolate to continuum spacetime. All necessary renormalizations are performed nonperturbatively., Comment: 9 pages, 2 figures, 1 table. v2: published version, one reference added

We determine the ratio FK/Fpi in QCD with Nf=2+1 flavors of sea quarks, based on a series of lattice calculations with three different lattice spacings, large volumes and a simulated pion mass reaching down to about 190 MeV. We obtain FK/Fpi=1.192 +/- 0.007(stat) +/- 0.006(syst). This result is then used to give an updated value of the CKM matrix element |Vus|. The unitarity relation for the first row of this matrix is found to be well observed., Comment: 15 pages, 4 figures, 2 tables

More than 99% of the mass of the visible universe is made up of protons and neutrons. Both particles are much heavier than their quark and gluon constituents, and the Standard Model of particle physics should explain this difference. We present a full ab-initio calculation of the masses of protons, neutrons and other light hadrons, using lattice quantum chromodynamics. Pion masses down to 190 mega electronvolts are used to extrapolate to the physical point with lattice sizes of approximately four times the inverse pion mass. Three lattice spacings are used for a continuum extrapolation. Our results completely agree with experimental observations and represent a quantitative confirmation of this aspect of the Standard Model with fully controlled uncertainties., Comment: 22 pages, 3 Tables, 8 Figures. Published in Science (21 November 2008) with Supporting Online Material. Submission to arXiv has been delayed by 6 months to respect the journal's embargo policy

We extend our previous study [Phys. Lett. B643 (2006) 46] of the cross-over temperatures (T_c) of QCD. We improve our zero temperature analysis by using physical quark masses and finer lattices. In addition to the kaon decay constant used for scale setting we determine four quantities (masses of the \Omega baryon, K^*(892) and \phi(1020) mesons and the pion decay constant) which are found to agree with experiment. This implies that --independently of which of these quantities is used to set the overall scale-- the same results are obtained within a few percent. At finite temperature we use finer lattices down to a <= 0.1 fm (N_t=12 and N_t=16 at one point). Our new results confirm completely our previous findings. We compare the results with those of the 'hotQCD' collaboration., Comment: 19 pages, 8 figures, 3 tables

We present a framework for phenomenological lattice QCD calculations which makes use of a tree level Symanzink improved action for gluons and stout-link Wilson fermions. We give details of our efficient HMC/RHMC algorithm and present a scaling study of the low-lying N_f=3 baryon spectrum. We find a scaling region that extends to a~<0.16fm and conclude that our action and algorithm are suitable for large scale phenomenological investigations of N_f=2+1 QCD. We expect this conclusion to hold for other comparable actions., Comment: 11 pages, 7 figures, aps style. Final published version

Some algorithmic details of our $N_f=2+1$ QCD mixed action simulations with overlap valence and improved Wilson sea quarks are presented., Comment: 7 pages, 6 figures. Presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4 2007, Regensburg, Germany

We present preliminary results for the chiral behavior of charged pseudo-Goldstone-boson masses and decay constants. These are obtained in simulations with N_f=2+1 flavors of tree-level, O(a)-improved Wilson sea quarks. In these simulations, mesons are composed of either valence quarks discretized in the same way as the sea quarks (unitary simulations) or of overlap valence quarks (mixed-action simulations). We find that the chiral behavior of the pseudoscalar meson masses in the mixed-action calculations cannot be explained with continuum, partially-quenched chiral perturbation theory. We show that the inclusion of O(a^2) unitarity violations in the chiral expansion resolves this discrepancy and that the size of the unitarity violations required are consistent with those which we observe in the zero-momentum, scalar-isotriplet-meson propagator., Comment: 7 pages, 3 figures, talk by L. Lellouch at the XXV International Symposium on Lattice Field Theory (LATTICE 2007), 30 July - 4 August 2007, Regensburg, Germany

We present two improvements to our previous dynamical overlap HMC algorithm. We introduce a new method of differentiating the eigenvectors of the Kernel operator, which removes an instability in the fermionic force. Secondly, by simulating part of the fermion determinant exactly, without pseudo-fermions, we are able to increase the rate of topological tunnelling by a factor of more than ten, reducing the auto-correlation., Comment: 7 pages, contribution to lattice 2007 (algorithms)

We present algorithmic improvements to the overlap Hybrid Monte Carlo algorithm, including preconditioning techniques and improvements to the correction step, used when one of the eigenvalues of the Kernel operator changes sign, which is now O(\Delta t^2) exact., Comment: 6 pages, 3 figures; poster contribution at Lattice 2005(Algorithms and machines)