Exploration Of Nucleon Structure In Lattice Qcd With Chiral Quarks

In this work, we calculate various nucleon structure observables using the fundamental theory of quarks and gluons, QCD, simulated on a lattice. In our simulations, we use the full QCD action including Nf = 2+ 1 dynamical quarks in the SU(2) isospin limit. We compute the nucleon vector and axial vector form factors as well as the generalized form factors, and analyze the nucleon charge, magnetization, and axial radii, anomalous magnetic moment, and axial charge. In addition, we compute quark contributions to the nucleon momentum and spin. Our calculation is novel for three reasons. It is a first full QCD calculation using both sea and valence chiral quarks with pion masses as low as m[pi] = 300 MeV. We develop a method to keep systematic effects in the lattice nucleon matrix elements under control, which helps us to obtain a better signal-to-noise ratio, to achieve higher precision and to test the applicability of low-energy effective theories. Finally, we compare the results from lattice QCD calculations with two different discretization methods and lattice spacings, with the rest of the calculation technique kept equal. The level of agreement between these results indicates that our calculations are not significantly affected by discretization effects.

Lattice gauge theory is a valuable tool for understanding how properties of the nucleon arise from the fundamental interactions of QCD. Numerical computations on the lattice can be used not only for first principles calculations of experimentally accessible quantities, but also for calculations of quantities that are not (yet) known from experiment. This thesis presents two lattice studies of the quark substructure of nucleons. The first study used overlaps calculated on the lattice to evaluate the goodness of trial nucleon sources. A variational study was performed to find the trial source that best approximated the true nucleon ground state. In this exploratory work with relatively simple trial sources on quenched lattices, we obtained overlaps as high as 80%. The second study was performed using domain wall valence fermions on Asqtad improved staggered lattices provided by the MILC collaboration, with pion masses as low as 290 MeV. We compute nucleon matrix elements of local quark operators: (F', S'l@P(0) F{Il Dt12 ... i D 0 (0)P, S), where F" E {y", -y"-y, -io*}. These operators are parameterized by generalized form factors, which in the infinite momentum frame can be unambiguously interpreted in terms of Fourier transforms of the transverse spatial distributions of quarks in a nucleon. By calculating the local operators at many different values of nucleon momentum, we extract a complete set of generalized form factors for the lowest two moments of the vector, axial and tensor operators. From the form factors, we compute a variety of quantities characterizing the internal structure of the nucleon. Finally, we explore chiral extrapolations of the lattice results to the physical pion mass.

It is emphasized in the 2015 NSAC Long Range Plan [1] that "understanding the structure of hadrons in terms of QCD's quarks and gluons is one of the central goals of modern nuclear physics." Over the last three decades, lattice QCD has developed into a powerful tool for ab initio calculations of strong-interaction physics. Up until now, it is the only theoretical approach to solving QCD with controlled statistical and systematic errors. Since 1985, we have proposed and carried out first-principles calculations of nucleon structure and hadron spectroscopy using lattice QCD which entails both algorithmic development and large scale computer simulation. We started out by calculating the nucleon form factors − electromagnetic [2], axial-vector [3], ? NN [4], and scalar [5] form factors, the quark spin contribution [6] to the proton spin, the strangeness magnetic moment [7], the quark orbital angular momentum [8], the quark momentum fraction [9], and the quark and glue decomposition of the proton momentum and angular momentum [10]. These first round of calculations were done with Wilson fermions in the q̀uenched' approximation where the dynamical effects of the quarks in the sea are not taken into account in the Monte Carlo simulation to generate the background gauge configurations. Beginning in 2000, we have started implementing the overlap fermion formulation into the spectroscopy and structure calculations [11, 12]. This is mainly because the overlap fermion honors chiral symmetry as in the continuum. It is going to be more and more important to take the symmetry into account as the simulations move closer to the physical point where the u and d quark masses are as light as a few MeV only. We began with lattices which have quark masses in the sea corresponding to a pion mass at ̃300 MeV and obtained the strange form factors [13], charm and strange quark masses, the charmonium spectrum and the Ds meson decay constant fDs [14], the strangeness and charmness [15], the meson mass decomposition [16] and the strange quark spin from the anomalous Ward identity [17]. Recently, we have started to include multiple lattices with different lattice spacings and different volumes including large lattices at the physical pion mass point. We are getting quite close to being able to calculate the hadron structure at the physical point and to do the continuum and large volume extrapolations which is our ultimate aim. We have now finished several projects which have included these systematic corrections. They include the leptonic decay width of the [18], the N sigma and strange sigma terms [19], and the strange quark magnetic moment [20]. Over the years, we have also studied hadron spectroscopy with lattice calculations and in phenomenology. These include Roper resonance [21, 22], pentaquark state [23], charmonium spectrum [24, 14], glueballs [25, 26, 27, 28], scalar mesons a0(1450) and (600) [29] and other scalar mesons [30], and the 1−+ meson [31]. In addition, we have employed the canonical approach to explore the first order phase transition and the critical point at finite density and finite temperature [32, 33]. We have also discovered a new parton degree of freedom − the connected sea partons, from the path-integral formulation of the hadronic tensor [34, 35] which explains the experimentally observed Gottfried sum rule violation [34]. Combining experimental result on the strange parton distribution, the CT10 global fitting results of the total u and d anti-partons and the lattice result of the ratio of the momentum fraction of the strange vs that of u or d in the disconnected insertion, we have shown that the connected sea partons can be isolated [36]. In this final technical report, we shall present a few representative highlights that have been achieved in the project.

We report RBC and RBC/UKQCD lattice QCD numerical calculations of nucleon electroweak matrix elements with dynamical domain-wall fermions (DWF) quarks. The first, RBC, set of dynamical DWF ensembles employs two degenerate flavors of DWF quarks and the DBW2 gauge action. Three sea quark mass values of 0.04, 0.03 and 0.02 in lattice units are used with about 200 gauge configurations each. The lattice cutoff is about 1.7 GeV and the spatial volume is about (1.9 fm)3. Despite the small volume, the ratio of the isovector vector and axial charges g{sub A}/g{sub V} and that of structure function moments x{sub u-d}/x{sub [Delta] u - [Delta] d} are in agreement with experiment, and show only very mild quark mass dependence. The second, RBC/UK, set of ensembles employs one strange and two degenerate (up and down) dynamical DWF quarks and Iwasaki gauge action. The strange quark mass is set at 0.04, and three up/down mass values of 0.03, 0.02 and 0.01 in lattice units are used. The lattice cutoff is about 1.6 GeV and the spatial volume is about (3.0 fm)3. Even with preliminary statistics of 25-30 gauge configurations, the ratios g{sub A}/g{sub V} and x{sub u-d}/x{sub {Delta} u - {Delta} d} are consistent with experiment and show only very mild quark mass dependence. Another structure function moment, d1, though yet to be renormalized, appears small in both sets.

The formalism of Quantum Chromodynamics on the lattice (or Lattice QCD) allows to perform ab-initio non-perturbative studies of strong-interaction driven processes, as it provides both a covariant regularisation of the theory of QCD and a natural framework for numerical computations. In this work, after a review of the main features of QCD and a step-by-step presentation of our discretization of QCD on a lattice, we undertake detailed studies of two problems of hadronic physics: the phenomenon of resonant scattering and the structure of the nucleon. The lattice calculations are performed with the Budapest-Marseille-Wuppertal Collaboration's 2+1-flavour gauge configurations, which give access to a wide range of lattice spacings, volumes and quarks masses, thereby allowing to study the sensibility of our results on these parameters, and to perform a complete continuum extrapolation. These configurations include dynamical quarks, and use a clover-improved Wilson QCD action. To investigate the scattering of particles on the lattice, we set up a Lüscher analysis for the emblematic case of pion-pion scattering in the channel of the rho meson resonance. We analyse our data with a variational generalized eigenvalue method, and give an in-depth calculation of the scattering phase-shifts and the corresponding resonance parameters, with a full control of the systematic errors. Our results provide an important step for lattice studies of scattering states, as they are the first to be performed at the physical pion mass, where one can see the actual decay of the rho into two pions. The obtained rho meson parameters are in good agreement with the experimental values, and consistent with a weak pion mass dependence of the coupling between the rho and two pions. As for our probe of the structure of the nucleon, we present a complete extraction of the electroweak isovector form factors, with a comprehensive study of the electric charge squared radius and of the axial charge. Our analysis also feature data at the physical pion mass, which turns out to be crucial in order to perform safe extrapolations to the physical point, as the chiral perturbation theory predicts violent variations of these quantities near the massless-quarks point. Our calculation includes source and sink projections onto the nucleon state, as well as a combined fit method between the two-point and three-point correlation functions to control the contamination of our data by excited states. Although one would need more data to perform a high-accuracy determination of the nucleon radius and axial charge at the physical point with a relevant estimation of the systematic errors, the results we obtain are in good agreement with the experiment and suggest that the excited-state effects are under control. Our analysis also highlights that gauge configurations ensembles near the physical pion mass and with large volumes must be used in order to extract accurate information about the nucleon structure from lattice calculations.

We report RBC and RBC/UKQCD lattice QCD numerical calculations of nucleon electroweak matrix elements with dynamical domain-wall fermions (DWF) quarks. The first, RBC, set of dynamical DWF ensembles employs two degenerate flavors of DWF quarks and the DBW2 gauge action. Three sea quark mass values of 0.04, 0.03 and 0.02 in lattice units are used with 220 gauge configurations each. The lattice cutoff is a{sup -1} {approx} 1.7GeV and the spatial volume is about (1.9fm){sup 3}. Despite the small volume, the ratio of the isovector vector and axial charges g{sub A}/g{sub V} and that of structure function moments x{sub u-d}/x{sub {Delta}u-{Delta}d} are in agreement with experiment, and show only very mild quark mass dependence. The second, RBC/UK, set of ensembles employs one strange and two degenerate (up and down) dynamical DWF quarks and Iwasaki gauge action. The strange quark mass is set at 0.04, and three up/down mass values of 0.03, 0.02 and 0.01 in lattice units are used. The lattice cutoff is a{sup -1} {approx} 1.6GeV and the spatial volume is about (3.0fm){sup 3}. Even with preliminary statistics of 25-30 gauge configurations, the ratios g{sub A}/g{sub V} and x{sub u-d}/x{sub {Delta}u-{Delta}d} are consistent with experiment and show only very mild quark mass dependence. Another structure function moment, d{sub 1}, though yet to be renormalized, appears small in both sets.

With ever increasing computational resources and improvements in algorithms, new opportunities are emerging for lattice gauge theory to address key questions in strongly interacting systems, such as nuclear matter. Calculations today use dynamical gauge-field ensembles with degenerate light up/down quarks and the strange quark and it is possible now to consider including charm-quark degrees of freedom in the QCD vacuum. Pion masses and other sources of systematic error, such as finite-volume and discretization effects, are beginning to be quantified systematically. Altogether, an era of precision calculation has begun and many new observables will be calculated at the new computational facilities. The aim of this set of lectures is to provide graduate students with a grounding in the application of lattice gauge theory methods to strongly interacting systems and in particular to nuclear physics. A wide variety of topics are covered, including continuum field theory, lattice discretizations, hadron spectroscopy and structure, many-body systems, together with more topical lectures in nuclear physics aimed a providing a broad phenomenological background. Exercises to encourage hands-on experience with parallel computing and data analysis are included.

This volume is centered on recent developments in the exploration of hadronic structure through lepton scattering, in the description of hadron physics directly from lattice QCD and non-perturbative QCD models, and in efforts to strengthen the links between these activities. Specific topics that are covered include: parton distribution functions, polarized structure functions, generalized structure functions, nuclear effects, quark-hadron duality, electromagnetic form factors, structure functions and hadron properties from lattice QCD, and QCD models based on the Dyson-Schwinger equations.