Structure functions from the Compton amplitude

We have initiated a program to compute the Compton amplitude from lattice QCD with the Feynman-Hellman method. This amplitude is related to the structure function via a Fredholm integral equation of the first kind. It is known that these types of equations are inherently ill--posed - they are, e.g., extremely sensitive to perturbations of the system. We discuss two methods which are candidates to handle these problems: the model free inversion based on singular value decomposition and one Bayesian type approach. We apply the Bayesian method to currently available lattice data for the Compton amplitude

A Feynman-Hellmann approach to the spin structure of hadrons

We perform a Nf = 2 + 1 lattice QCD simulation to determine the quark spin fractions of hadrons using the Feynman-Hellmann theorem. By introducing an external spin operator to the fermion action, the matrix elements relevant for quark spin fractions are extracted from the linear response of the hadron energies. Simulations indicate that the Feynman-Hellmann method offers statistical precision that is comparable to the standard three-point function approach, with the added benefit that it is less susceptible to excited state contamination. This suggests that the Feynman-Hellmann technique offers a promising alternative for calculations of quark line disconnected contributions to hadronic matrix elements. At the SU(3)-flavour symmetry point, we find that the connected quark spin fractions are universally in the range 55-70% for vector mesons and octet and decuplet baryons. There is an indication that the amount of spin suppression is quite sensitive to the strength of SU(3) breaking.Comment: 13 pages, 7 figure

Determining the glue component of the nucleon

Computing the gluon component of momentum in the nucleon is a difficult and computationally expensive problem, as the matrix element involves a quark-line-disconnected gluon operator which suffers from ultra-violet fluctuations. But also necessary for a successful determination is the non-perturbative renormalisation of this operator. As a first step we investigate here this renormalisation in the RI-MOM scheme. Using quenched QCD as an example, a statistical signal is obtained in a direct calculation using an adaption of the Feynman-Hellmann technique.Comment: 7 pages, Proceedings of the 37th Annual International Symposium on Lattice Field Theory (Lattice 2019), 16-22 June 2019, Wuhan, Chin

Feynman--Hellmann approach to transition matrix elements and quasi-degenerate energy states

The Feynman--Hellmann approach to computing matrix elements in lattice QCD by first adding a perturbing operator to the action is described using the transition matrix and the Dyson expansion formalism. This perturbs the energies in the two-point baryon correlation function, from which the matrix element can be obtained. In particular at leading order in the perturbation we need to diagonalise a matrix of near-degenerate energies. While the method is general for all hadrons, we apply it here to a study of a Sigma to Nucleon baryon transition vector matrix element.Comment: 50 pages. Minor typos fixed. Published versio

Determining the glue component of the nucleon

Computing the gluon component of momentum in the nucleon is a difficult and computationally expensive problem, as the matrix element involves a quark-line-disconnected gluon operator which suffers from ultra-violet fluctuations. But also necessary for a successful determination is the non-perturbative renormalisation of this operator. As a first step we investigate here this renormalisation in the RI − MOM scheme. Using quenched QCD as an example, a statistical signal is obtained in a direct calculation using an adaption of the Feynman-Hellmann technique

Electromagnetic form factors at large momenta from lattice QCD

Accessing hadronic form factors at large momentum transfers has traditionally presented a challenge for lattice QCD simulations. Here, we demonstrate how a novel implementation of the Feynman-Hellmann method can be employed to calculate hadronic form factors in lattice QCD at momenta much higher than previously accessible. Our simulations are performed on a single set of gauge configurations with three flavors of degenerate mass quarks corresponding to mπ≈470  MeV. We are able to determine the electromagnetic form factors of the pion and nucleon up to approximately 6  GeV2, with results for the ratio of the electric and magnetic form factors of the proton at our simulated quark mass agreeing well with experimental results.A.J. Chambers, J. Dragos, R. Horsley, Y. Nakamura, H. Perlt, D. Pleiter, P.E.L. Rakow, G. Schierholz, A. Schiller, K. Somfleth, H. Stüben, R.D. Young and J.M. Zanott

Moments and power corrections of longitudinal and transverse proton structure functions from lattice QCD

We present a simultaneous extraction of the moments of $F_2$ and $F_L$ structure functions of the proton at a range of photon virtuality, $Q^2$. This is achieved by computing the forward Compton amplitude via an application of the second-order Feynman-Hellmann method. We find the moments of $F_{2,L}$ in good agreement with experimental values. By studying the $Q^2$ dependence of $F_2$ moments, we estimate the power corrections

Anomalous magnetic moment of the muon with dynamical QCD+QED

The current $3.5\sigma$ discrepancy between experimental and Standard Model determinations of the anomalous magnetic moment of the muon $a_\mu=(g-2)/2$ can only be extended to the discovery $5\sigma$ regime through a reduction of both experimental and theoretical uncertainties. On the theory side, this means a determination of the hadronic vacuum polarisation (HVP) contribution to better than 0.5%, a level of precision that demands the inclusion of QCD + QED effects to properly understand how the behaviour of quarks are modified when their electric charges are turned on. The QCDSF collaboration has generated an ensemble of configurations with dynamical QCD and QED fields with the specific aim of studying flavour breaking effects arising from differences in the quark masses and charges in physical quantities. Here we study these effects in a calculation of HVP around the SU(3) symmetric point. Furthermore, by performing partially-quenched simulations we are able to cover a larger range of quark masses and charges on these configurations and then fit the results to an SU(3) flavour breaking expansion. Subsequently, this allows for an extrapolation to the physical point

Patterns of flavor symmetry breaking in hadron matrix elements involving u, d, and s quarks

By considering a flavor expansion about the SU(3) flavor symmetric point, we investigate how flavor blindness constrains octet baryon matrix elements after SU(3) is broken by the mass difference between quarks. Similarly to hadron masses we find the expansions to be constrained along a mass trajectory where the singlet quark mass is held constant, which provides invaluable insight into the mechanism of flavor symmetry breaking and proves beneficial for extrapolations to the physical point. Expansions are given up to third order in the expansion parameters. Considering higher orders would give no further constraints on the expansion parameters. The relation of the expansion coefficients to the quark-line-connected and quark-line-disconnected terms in the three-point correlation functions is also given. As we consider Wilson cloverlike fermions, the addition of improvement coefficients is also discussed and shown to be included in the formalism developed here. As an example of the method we investigate this numerically via a lattice calculation of the flavor-conserving matrix elements of the vector first-class form factors

Feynman-Hellmann approach to transition matrix elements and quasi-degenerate energy states

The Feynman-Hellmann approach to computing matrix elements in lattice QCD by first adding a perturbing operator to the action is described using the transition matrix and the Dyson expansion formalism. This perturbs the energies in the two-point baryon correlation function, from which the matrix element can be obtained. In particular at leading order in the perturbation we need to diagonalise a matrix of near-degenerate energies. While the method is general for all hadrons, we apply it here to a study of a Sigma to Nucleon baryon transition vector matrix element