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

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

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

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

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

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

Towards high partial waves in lattice QCD with a dumbbell-like operator

An extended two-hadron operator is developed to extract the spectra of irreducible representations (irreps) in the finite volume. The irreps of the group for the finite volume system are projected using a coordinate-space operator. The correlation function of this operator is computationally efficient to extract lattice spectra of the specific irrep. In particular, this new formulation only requires propagators to be computed from two distinct source locations, at fixed spatial separation. We perform a proof-of-principle study on a $24^3 \times 48$ lattice volume with $m_\pi\approx 900$ MeV by isolating various spectra of the $\pi\pi$ system with isospin-2 including a range of total momenta and irreps. By applying the Lüscher formalism, the phase shifts of $S$-, $D$- and $G$-wave $\pi\pi$ scattering with isospin-2 are extracted from the spectra

Towards high partial waves in lattice QCD with a dumbbell-like operator

An extended two-hadron operator is developed to extract the spectra of irreducible representations (irreps) in the finite volume. The irreps of the group for the finite volume system are projected using a coordinate-space operator. The correlation function of this operator is computationally efficient to extract lattice spectra of the specific irrep. In particular, this new formulation only requires propagators to be computed from two distinct source locations, at fixed spatial separation. We perform a proof-of-principle study on a $24^3 \times 48$ lattice volume with $m_\pi\approx 900$ MeV by isolating various spectra of the $\pi\pi$ system with isospin-2 including a range of total momenta and irreps. By applying the L\"uscher formalism, the phase shifts of $S$-, $D$- and $G$-wave $\pi\pi$ scattering with isospin-2 are extracted from the spectra.Comment: 29 Pages, 5 figure

Towards high partial waves in lattice QCD with an extended two-hadron operator

An extended two-hadron operator is developed to extract the spectra of irreducible representations (irreps) in the finite volume. The irreps of the group for the finite volume system are projected using a coordinate-space operator. The correlation function of this operator is computationally efficient to extract lattice spectra of the specific irrep. In particular, this new formulation only requires propagators to be computed from two distinct source locations, at fixed spatial separation. We perform a proof-of-principle study on a 243×48 lattice volume with mπ≈900 MeV by isolating various spectra of the ππ system with isospin-2 including a range of total momenta and irreps. By applying the Lüscher formalism, the phase shifts of S- and D-wave ππ scattering with isospin-2 are extracted from the spectra, with a tentative look at the role and influence of the G-wave