• Nucleon electromagnetic form factors in two-flavour QCD(1504.04628)

April 26, 2015 hep-ph, nucl-th, hep-lat
We present results for the nucleon electromagnetic form factors, including the momentum transfer dependence and derived quantities (charge radii and magnetic moment). The analysis is performed using O(a) improved Wilson fermions in Nf=2 QCD measured on the CLS ensembles. Particular focus is placed on a systematic evaluation of the influence of excited states in three-point correlation functions, which lead to a biased evaluation, if not accounted for correctly. We argue that the use of summed operator insertions and fit ans\"atze including excited states allow us to suppress and control this effect. We employ a novel method to perform joint chiral and continuum extrapolations, by fitting the form factors directly to the expressions of covariant baryonic chiral effective field theory. The final results for the charge radii and magnetic moment from our lattice calculations include, for the first time, a full error budget. We find that our estimates are compatible with experimental results within their overall uncertainties.
• Nucleon axial form factors from two-flavour Lattice QCD(1411.5828)

Nov. 21, 2014 hep-lat
We present preliminary results on the axial form factor $G_A(Q^2)$ and the induced pseudoscalar form factor $G_P(Q^2)$ of the nucleon. A systematic analysis of the excited-state contributions to form factors is performed on the CLS ensemble `N6' with $m_\pi = 340 \ \text{MeV}$ and lattice spacing $a \sim 0.05 \ \text{fm}$. The relevant three-point functions were computed with source-sink separations ranging from $t_s \sim 0.6 \ \text{fm}$ to $t_s \sim \ 1.4 \ \text{fm}$. We observe that the form factors suffer from non-trivial excited-state contributions at the source-sink separations available to us. It is noted that naive plateau fits underestimate the excited-state contributions and that the method of summed operator insertions correctly accounts for these effects.
• Suppression of excited-state effects in lattice determination of nucleon electromagnetic form factors(1411.4804)

Nov. 18, 2014 hep-lat
We study the ability of a variety of fitting techniques to extract the ground state matrix elements of the vector current from ratios of nucleon three- and two-point functions that contain contaminations from excited states. Extending our high-statistics study of nucleon form factors, we are able to demonstrate that the treatment of excited-state contributions in conjunction with approaching the physical pion mass has a significant impact on the $Q^2$-dependence of the form factors.
• A high-statistics study of the nucleon EM form factors, axial charge and quark momentum fraction(1311.5804)

April 8, 2014 hep-lat
We present updated results for the nucleon axial charge and electromagnetic (EM) form factors, which include a significant increase in statistics for all ensembles (up to 4000 measurements), as well as the addition of ensembles with pion masses down to $m_\pi\sim195$ MeV. We also present results for the average quark momentum fraction. The new data allows us to perform a thorough study of the systematic effects encountered in the lattice extraction. We concentrate on systematic effects due to excited-state contaminations for each of the quantities, which we check using several different time separations between the operators at the source and sink through a comparison of plateau fits and the summed operator insertion method (which provides a mechanism to suppress the excited-state contamination). We confirm our earlier finding that a reliable extraction of the axial charge must be based on a method which eliminates excited-state contaminations. Similar conclusions apply to our EM form factor calculations . The measurements are calculated using the CLS ensembles with non-perturbatively O(a) improved Wilson fermions in $N_f=2$ QCD.
• The nucleon axial charge from lattice QCD with controlled errors(1205.0180)

Oct. 2, 2012 hep-ph, nucl-th, hep-lat
We report on our calculation of the nucleon axial charge gA in QCD with two flavours of dynamical quarks. A detailed investigation of systematic errors is performed, with a particular focus on contributions from excited states to three-point correlation functions. The use of summed operator insertions allows for a much better control over such contamination. After performing a chiral extrapolation to the physical pion mass, we find gA=1.223 +/- 0.063 (stat) +0.035 -0.060 (syst), in good agreement with the experimental value.
• Form factors in lattice QCD(1106.1554)

June 8, 2011 hep-ph, nucl-th, hep-lat
Lattice simulations of QCD have produced precise estimates for the masses of the lowest-lying hadrons which show excellent agreement with experiment. By contrast, lattice results for the vector and axial vector form factors of the nucleon show significant deviations from their experimental determination. We present results from our ongoing project to compute a variety of form factors with control over all systematic uncertainties. In the case of the pion electromagnetic form factor we employ partially twisted boundary conditions to extract the pion charge radius directly from the linear slope of the form factor near vanishing momentum transfer. In the nucleon sector we focus specifically on the possible contamination from contributions of higher excited states. We argue that summed correlation functions offer the possibility of eliminating this source of systematic error. As an illustration of the method we discuss our results for the axial charge, gA, of the nucleon.
• Wilson fermions at fine lattice spacings: scale setting, pion form factors and (g-2)_mu(1010.2390)

Oct. 12, 2010 hep-ph, nucl-th, hep-lat
We present an update on our on-going project to compute hadronic observables for Nf=2 flavours of O(a) improved Wilson fermions at small lattice spacings. The procedure to determine the lattice scale via the mass of the Omega baryon is described. Furthermore we present preliminary results for the pion form factor computed using twisted boundary conditions, and report on the implementation of a novel approach to determine the contribution of the hadronic vacuum polarisation to the anomalous magnetic moment of the muon.
• Mesonic and baryonic correlation functions at fine lattice spacings(0910.5578)

Oct. 29, 2009 hep-lat
We report on our on-going project to compute mesonic and baryonic two- and three-point correlation functions in simulations using Nf=2 flavours of O(a) improved Wilson quarks and the Wilson plaquette action. We present performance figures for the DD-HMC algorithm on commodity cluster hardware and discuss the issue of critical slowing down, which is particularly pronounced for the topological charge. The effectiveness of stochastic noise sources and Jacobi smearing are investigated. Our preliminary results obtained at three quark masses on 96x48^3 at beta=5.5 imply that the lattice spacing is about 0.06 fm, while the smallest pion mass in the current runs is around 360 MeV, which corresponds to m_pi*L=5.3.
• Parton Distribution Functions with Twisted Mass Fermions(hep-lat/0511013)

Nov. 4, 2005 hep-lat
We present a first Wilson twisted mass fermion calculation of the matrix element between pion states of the twist-2 operator, which is related to the the lowest moment < x > of the valence quark parton distribution function in a pion. Using Wilson twisted mass fermions in the quenched approximation we demonstrate that < x > can be computed at small pseudoscalar meson masses down to values of order 250 MeV. We investigate the scaling behaviour of this physically important quantity by applying two definitions of the critical mass and observe a scaling compatible with the expected O(a^2) behaviour in both cases. A combined continuum extrapolation allows to obtain reliable results for < x > at very small pseudoscalar meson masses, which previously could not be explored by lattice QCD simulations.
• 3-point functions from twisted mass lattice QCD at small quark masses(hep-lat/0510023)

Oct. 4, 2005 hep-lat
We show at the example of the matrix element between pion states of a twist-2, non-singlet operator that Wilson twisted mass fermions allow to compute this phenomenologically relevant quantitiy at small pseudo scalar masses of O(270 MeV). In the quenched approximation, we investigate the scaling behaviour of this observable that is derived from a 3-point function by applying two definitions of the critical mass and find a scaling compatible with the expected O(a^2) behaviour in both cases. A combined continuum extrapolations allows to obtain reliable results at small pion masses, which previously could not be explored by lattice QCD simulations.
• Comparison between overlap and twisted mass fermions towards the chiral limit(hep-lat/0409109)

Sept. 17, 2004 hep-lat
We compare overlap fermions, which are chirally invariant, and Wilson twisted mass fermions in the approach to the chiral limit. Our quenched simulations reveal that with both formulations of lattice fermions pion masses of O(250 MeV) can be reached in practical simulations. Our comparison is done at a fixed lattice spacing a=0.123 fm. Several quantities are measured, such as hadron masses and pseudoscalar decay constants.
• Extracting physics from an unphysical situation: light mesons in a small box(physics/0309072)

Sept. 16, 2003 hep-ph, hep-lat, physics.comp-ph
Quantum Chromo Dynamics is considered in a setup where the light mesons are squeezed into unphysically small boxes. We show how such a situation can be used to determine the couplings of the low energy chiral Lagrangian from lattice simulations, applying chirally invariant formulations of lattice fermions.
• Continuum Extrapolation of Moments of Nucleon Quark Distributions in Full QCD(hep-lat/0211021)

Nov. 11, 2002 hep-lat
Moments of light cone quark density, helicity, and transversity distributions are calculated in unquenched lattice QCD at $\beta = 5.5$ and $\beta = 5.3$ using Wilson fermions on $16^3 \times 32$ lattices. These results are combined with earlier calculations at $\beta = 5.6$ using SESAM configurations to study the continuum limit.
• Moments of Nucleon Light Cone Quark Distributions Calculated in Full Lattice QCD(hep-lat/0201021)

May 29, 2002 hep-ph, nucl-th, hep-lat
Moments of the quark density, helicity, and transversity distributions are calculated in unquenched lattice QCD. Calculations of proton matrix elements of operators corresponding to these moments through the operator product expansion have been performed on $16^3 \times 32$ lattices for Wilson fermions at $\beta = 5.6$ using configurations from the SESAM collaboration and at $\beta = 5.5$ using configurations from SCRI. One-loop perturbative renormalization corrections are included. At quark masses accessible in present calculations, there is no statistically significant difference between quenched and full QCD results, indicating that the contributions of quark-antiquark excitations from the Dirac Sea are small. Close agreement between calculations with cooled configurations containing essentially only instantons and the full gluon configurations indicates that quark zero modes associated with instantons play a dominant role. Naive linear extrapolation of the full QCD calculation to the physical pion mass yields results inconsistent with experiment. Extrapolation to the chiral limit including the physics of the pion cloud can resolve this discrepancy and the requirements for a definitive chiral extrapolation are described.
• Progress towards a lattice determination of (moments of) nucleon structure functions(hep-lat/0111012)

Nov. 8, 2001 hep-ph, hep-lat
Using unimproved and non-perturbatively O(a) improved Wilson fermions, results are given for the three lowest moments of unpolarised nucleon structure functions. Renormalisation, chiral extrapolation and the continuum limit of the matrix elements are briefly discussed. The simulations are performed for both quenched and two flavours of unquenched fermions. No obvious sign of deviation from linearity in the chiral extrapolations are found. (This is most clearly seen in our quenched unimproved data, which extends to lighter quark mass.) Possible quenching effects also seem to be small. The lowest moment thus remains too large, so it seems to be necessary to reach smaller quark masses in numerical simulations.
• Analysis of the DeltaI=1/2 Rule and e'/e with Overlap Fermions(hep-lat/0011070)

April 4, 2001 hep-ph, hep-lat
We study the renormalization of the DeltaS=1 effective weak Hamiltonian with overlap fermions. The mixing coefficients among dimension-six operators are computed at one loop in perturbation theory. As a consequence of the chiral symmetry at finite lattice spacing and of the GIM mechanism, which turns out to be quadratic in the masses, the K-->pi pi and K -->pi matrix elements relevant for the Delta I =1/2 rule can be computed without any power subtractions. The analogous amplitudes for e'/e require one divergent subtraction only, which can be performed non-perturbatively using K-->0 matrix elements.
• Moments of Structure Functions in Full QCD(hep-lat/0011010)

Nov. 1, 2000 hep-lat
Moments of the quark density distribution, moments of the quark helicity distribution, and the tensor charge are calculated in full QCD. Calculations of matrix elements of operators from the operator product expansion have been performed on $16^3 \times 32$ lattices for Wilson fermions at $\beta = 5.6$ using configurations from the SESAM collaboration and at $\beta = 5.5$ using configurations from SCRI. One-loop perturbative renormalization corrections are included. Selected results are compared with corresponding quenched calculations and with calculations using cooled configurations.

Oct. 26, 2000 hep-lat
We present calculations of matrix elements of 4-quark operators in the pion and in the nucleon extracted from quenched Monte Carlo simulations at beta = 6.0 using Wilson fermions. These operators are relevant for higher-twist effects. We are particularly careful to avoid mixing with lower-dimensional operators by choosing appropriate flavour structures.
• Perturbative Renormalization of Weak-Hamiltonian Four-Fermion Operators with Overlap Fermions(hep-lat/0007011)

July 10, 2000 hep-ph, hep-lat
The renormalization of the most general dimension-six four-fermion operators without power subtractions is studied at one loop in lattice perturbation theory using overlap fermions. As expected, operators with different chirality do not mix among themselves and parity-conserving and parity-violating multiplets renormalize in the same way. The renormalization constants of unimproved and improved operators are also the same. These mixing factors are necessary to determine physical matrix elements relevant to many phenomenological applications of weak interactions. The most important are the K0-K0bar and B0-B0bar mixings in the Standard Model and beyond, the Delta I =1/2 rule and epsilon'/epsilon.
• Renormalisation and off-shell improvement in lattice perturbation theory(hep-lat/0007004)

July 5, 2000 hep-lat
We discuss the improvement of flavour non-singlet point and one-link lattice quark operators, which describe the quark currents and the first moment of the DIS structure functions respectively. Suitable bases of improved operators are given, and the corresponding renormalisation factors and improvement coefficients are calculated in one-loop lattice perturbation theory, using the Sheikholeslami-Wohlert (clover) action. To this order we achieve off-shell improvement by eliminating the effect of contact terms. We use massive fermions, and our calculations are done keeping all terms up to first order in the lattice spacing, for arbitrary m^2/p^2, in a general covariant gauge. We also compare clover fermions with fermions satisfying the Ginsparg-Wilson relation, and show how to remove O(a) effects off-shell in this case too, and how this is in many aspects simpler than for clover fermions. Finally, tadpole improvement is also considered.
• Higher-Twist Contribution to Pion Structure Function: 4-Fermi Operators(hep-lat/9908011)

Dec. 30, 1999 hep-ph, hep-lat
We present quenched lattice QCD results for the contribution of higher-twist operators to the lowest non-trivial moment of the pion structure function. To be specific, we consider the combination $F_2^{\pi^+} + F_2^{\pi^-} - 2 F_2^{\pi^0}$ which has $I = 2$ and receives contributions from 4-Fermi operators only. We introduce the basis of lattice operators. The renormalization of the operators is done perturbatively in the $\bar{\rm{MS}}$ scheme using the 't Hooft-Veltman prescription for $\gamma_5$, taking particular care of mixing effects. The contribution is found to be of $O(f_\pi^2/Q^2)$, relative to the leading contribution to the moment of $F_2^{\pi^+}$.
• On-shell and off-shell improvement for Ginsparg-Wilson fermions(hep-lat/9909167)

Sept. 29, 1999 hep-lat
We discuss the improvement of bilinear fermionic operators for Ginsparg-Wilson fermions. We present explicit formulae for improved Green's functions, which apply both on-shell and off-shell.
• Higher-twist contributions to the Structure Functions coming from 4-fermion operators(hep-lat/9908038)

Aug. 24, 1999 hep-lat
We evaluate the contribution of a class of higher-twist operators to the lowest moment of the Structure Functions, by computing appropriate matrix elements of six four-fermion operators in the quenched approximation. Their perturbative renormalization constants and mixing coefficients are calculated in the 't Hooft-Veltman scheme of dimensional regularization, using codes written in the algebraic manipulation computer language FORM.
• Operator improvement for Ginsparg-Wilson fermions(hep-lat/9908029)

Aug. 20, 1999 hep-lat
The improvement of fermionic operators for Ginsparg-Wilson fermions is investigated. We present explicit formulae for improved Green's functions, which apply both on-shell and off-shell.
• Higher-twist corrections to nucleon structure functions from lattice QCD(hep-ph/9906320)

June 10, 1999 hep-ph, hep-lat
A genuinely non-perturbative evaluation of higher-twist contributions to the structure functions of the nucleon, with all mixing effects and renormalon ambiguities taken care of, is presented. Higher-twist corrections turn out to be significant at moderate values of $q^2$.