Search

Some calculations of parton distributions from first principles only give access to a limited range of Fourier modes of the function to reconstruct. We present a physically motivated procedure to regularize the inverse integral problem using a Gaussian process as a Bayesian prior. We propose to fix the hyperparameters of the prior in a meaningful physical fashion, offering a simple implementation, great numerical efficiency, and allowing us to understand and keep control easily of the uncertainty of the reconstruction., Comment: 14 pages, 13 figures

We present an exploration of the unpolarized isovector proton generalized parton distributions (GPDs) $H^{u-d}(x, \xi, t)$ and $E^{u-d}(x, \xi, t)$ in the pseudo-distribution formalism using distillation. Taking advantage of the large kinematic coverage made possible by this approach, we present results on the moments of GPDs up to the order $x^3$ -- including their skewness dependence -- at a pion mass $m_\pi = 358$ MeV and a lattice spacing $a = 0.094$ fm., Comment: added discussion on isovector D-term and stability of fits with respect to the range in non-local separation of the operator

This paper presents a comprehensive review of both the theory and experimental successes of Quantum Chromodynamics, starting with its emergence as a well defined theory in 1972-73 and following developments and results up to the present day. Topics include a review of the earliest theoretical and experimental foundations; the fundamental constants of QCD; an introductory discussion of lattice QCD, the only known method for obtaining exact predictions from QCD; methods for approximating QCD, with special focus on effective field theories; QCD under extreme conditions; measurements and predictions of meson and baryon states; a special discussion of the structure of the nucleon; techniques for study of QCD at high energy, including treatment of jets and showers; measurements at colliders; weak decays and quark mixing; and a section on the future, which discusses new experimental facilities or upgrades currently funded. The paper is intended to provide a broad background for Ph.D. students and postdocs starting their career. Some contributions include personal accounts of how the ideas or experiments were developed., Comment: Invited volume for the EJPC; 567 pages if text/figures and 4783 references occupying about 160 additional pages. arXiv abstract abridged, for the complete abstract please see the full text

The calculation of disconnected diagram contributions to physical signals is a computationally expensive task in Lattice QCD. To extract the physical signal, the trace of the inverse Lattice Dirac operator, a large sparse matrix, must be stochastically estimated. Because the variance of the stochastic estimator is typically large, variance reduction techniques must be employed. Multilevel Monte Carlo (MLMC) methods reduce the variance of the trace estimator by utilizing a telescoping sequence of estimators. Frequency Splitting is one such method that uses a sequence of inverses of shifted operators to estimate the trace of the inverse lattice Dirac operator, however there is no a priori way to select the shifts that minimize the cost of the multilevel trace estimation. In this article, we present a sampling and interpolation scheme that is able to predict the variances associated with Frequency Splitting under displacements of the underlying space time lattice. The interpolation scheme is able to predict the variances to high accuracy and therefore choose shifts that correspond to an approximate minimum of the cost for the trace estimation. We show that Frequency Splitting with the chosen shifts displays significant speedups over multigrid deflation, and that these shifts can be used for multiple configurations within the same ensemble with no penalty to performance.

The non-singlet helicity quark parton distribution functions (PDFs) of the nucleon are determined from lattice QCD, by jointly leveraging pseudo-distributions and the distillation spatial smearing paradigm. A Lorentz decomposition of appropriately isolated space-like matrix elements reveals pseudo-distributions that contain information on the leading-twist helicity PDFs, as well as an invariant amplitude that induces an additional $z^2$ contamination of the leading-twist signal. An analysis of the short-distance behavior of the space-like matrix elements using matching coefficients computed to next-to-leading order (NLO) exposes the desired PDF up to this additional $z^2$ contamination. Due to the non-conservation of the axial current, we elect to isolate the helicity PDFs normalized by the nucleon axial charge at the same scale $\mu^2$. The leading-twist helicity PDFs as well as several sources of systematic error, including higher-twist effects, discretization errors, and the aforementioned $z^2$ contaminating amplitude are jointly determined by characterizing the computed pseudo-distribution in a basis of Jacobi polynomials. The Akaike Information Criterion is exploited to effectively average over distinct model parameterizations and cuts on the pseudo-distribution. Encouraging agreement is observed with recent global analyses of each non-singlet quark helicity PDF, notably a rather small non-singlet anti-quark helicity PDF for all quark momentum fractions., Comment: 32 pages, 38 figures

We present the first exploratory lattice quantum chromodynamics (QCD) calculation of the polarized gluon Ioffe-time pseudo-distribution in the nucleon. The Ioffe-time pseudo-distribution provides a frame-independent and gauge-invariant framework to determine the gluon helicity in the nucleon from first principles. We employ a high-statistics computation using a $32^3\times 64$ lattice ensemble characterized by a $358$ MeV pion mass and a $0.094$ fm lattice spacing. We establish the pseudo-distribution approach as a feasible method to address the proton spin puzzle with successive improvements in statistical and systematic uncertainties anticipated in the future. Within the statistical precision of our data, we find a good comparison between the lattice determined polarized gluon Ioffe-time distribution and the corresponding expectations from the state-of-the-art global analyses. We find a hint for a nonzero gluon spin contribution to the proton spin from the model-independent extraction of the gluon helicity pseudo-distribution over a range of Ioffe-time, $\nu\lesssim 9$., Comment: Version to be published in Phys. Rev. D, 24 pages, 13 figures

In this document, we summarize the status and challenges of calculating parton physics in lattice QCD for the US Particle Physics Community Planning Exercise (a.k.a. "Snowmass"). While PDF-moments calculations have been very successful and been continuously improved, new methods have been developed to calculate distributions directly in $x$-space. Many recent lattice studies have been focused on calculating isovector PDFs of the pion and nucleon, learning to control systematics associated with excited-state contamination, renormalization and continuum extrapolations, pion-mass and finite-volume effects, etc. Although in some cases, the lattice results are already competitive with experimental data, to reach the level of precision in a wide range of $x$ for unpolarized nucleon PDFs impactful for future collider physics remains a challenge, and may require exascale supercomputing power. The new theoretical methods open the door for calculating other partonic observables which will be the focus of the experimental program in nuclear physics, including generalized parton distributions and transverse-momentum dependent PDFs. A fruitful interplay between experimental data and lattice-QCD calculations will usher in a new era for parton physics and hadron structure., Comment: 14 pages

Advances in machine learning methods provide tools that have broad applicability in scientific research. These techniques are being applied across the diversity of nuclear physics research topics, leading to advances that will facilitate scientific discoveries and societal applications. This Review gives a snapshot of nuclear physics research which has been transformed by machine learning techniques., Comment: Comments are welcome

We present a determination of the non-singlet transversity parton distribution function (PDF) of the nucleon, normalized with respect to the tensor charge at $\mu^2=2$ GeV$^2$ from lattice quantum chromodynamics. We apply the pseudo-distribution approach, using a gauge ensemble with a lattice spacing of 0.094 fm and the light quark mass tuned to a pion mass of 358 MeV. We extract the transversity PDF from the analysis of the short-distance behavior of the Ioffe-time pseudo-distribution using the leading-twist next-to-leading order (NLO) matching coefficients calculated for transversity. We reconstruct the $x$-dependence of the transversity PDF through an expansion in a basis of Jacobi polynomials in order to reduce the PDF ansatz dependence. Within the limitations imposed by a heavier-than-physical pion mass and a fixed lattice spacing, we present a comparison of our estimate for the valence transversity PDF with the recent global fit results based on single transverse spin asymmetry. We find the intrinsic nucleon sea to be isospin symmetric with respect to transversity., Comment: 25 pages, 13 figures

In this study, we present a determination of the unpolarized gluon Ioffe-time distribution in the nucleon from a first principles lattice quantum chromodynamics calculation. We carry out the lattice calculation on a $32^3\times 64$ ensemble with a pion mass of $358$ MeV and lattice spacing of $0.094$ fm. We construct the nucleon interpolating fields using the distillation technique, flow the gauge fields using the gradient flow, and solve the summed generalized eigenvalue problem to determine the glounic matrix elements. Combining these techniques allows us to provide a statistically well-controlled Ioffe-time distribution and unpolarized gluon PDF. We obtain the flow time independent reduced Ioffe-time pseudo-distribution, and calculate the light-cone Ioffe-time distribution and unpolarized gluon distribution function in the $\overline{\rm MS}$ scheme at $\mu = 2$ GeV, neglecting the mixing of the gluon operator with the quark singlet sector. Finally, we compare our results to phenomenological determinations., Comment: 23 pages, 12 figures, Phys Rev D accepted version

We apply the Distillation spatial smearing program to the extraction of the unpolarized isovector valence PDF of the nucleon. The improved volume sampling and control of excited-states afforded by distillation leads to a dramatically improved determination of the requisite Ioffe-time Pseudo-distribution (pITD). The impact of higher-twist effects is subsequently explored by extending the Wilson line length present in our non-local operators to one half the spatial extent of the lattice ensemble considered. The valence PDF is extracted by analyzing both the matched Ioffe-time Distribution (ITD), as well as a direct matching of the pITD to the PDF. Through development of a novel prescription to obtain the PDF from the pITD, we establish a concerning deviation of the pITD from the expected DGLAP evolution of the pseudo-PDF. The presence of DGLAP evolution is observed once more following introduction of a discretization term into the PDF extractions. Observance and correction of this discrepancy further highlights the utility of distillation in such structure studies., Comment: 34 pages, 26 figures

In this work, we study probing for the more general problem of computing the trace of a permutation of $A^{-1}$, say $PA^{-1}$. The motivation comes from Lattice QCD where we need to construct "disconnected diagrams" to extract flavor-separated Generalized Parton functions. In Lattice QCD, where the matrix has a 4D toroidal lattice structure, these non-local operators correspond to a $PA^{-1}$ where $P$ is the permutation relating to some displacement $\vec{p}$ in one or more dimensions. We focus on a single dimension displacement ($p$) but our methods are general. We show that probing on $A^k$ or $(PA)^k$ do not annihilate the largest magnitude elements. To resolve this issue, our displacement-based probing works on $PA^k$ using a new coloring scheme that works directly on appropriately displaced neighborhoods on the lattice. We prove lower bounds on the number of colors needed, and study the effect of this scheme on variance reduction, both theoretically and experimentally on a real-world Lattice QCD calculation. We achieve orders of magnitude speedup over the unprobed or the naively probed methods.

In this study, we present continuum limit results for the unpolarized parton distribution function of the nucleon computed in lattice QCD. This study is the first continuum limit using the pseudo-PDF approach with Short Distance Factorization for factorizing lattice QCD calculable matrix elements. Our findings are also compared with the pertinent phenomenological determinations. Inter alia, we are employing the summation Generalized Eigenvalue Problem (sGEVP) technique in order to optimize our control over the excited state contamination which can be one of the most serious systematic errors in this type of calculations. A crucial novel ingredient of our analysis is the parameterization of systematic errors using Jacobi polynomials to characterize and remove both lattice spacing and higher twist contaminations, as well as the leading twist distribution. This method can be expanded in further studies to remove all other systematic errors., Comment: 56 pages, 29 figures

In this work we present an efficient construction of baryon interpolating fields for lattice QCD computations of two and three point functions. These are essential building blocks of computations of nucleon parton distribution functions (PDFs), generalized parton distribution functions (GPDs) and transverse momentum dependent distributions functions (TMDs). Lattice QCD computations of these quantities can provide additional input to assist with the global fits on experimental data for determining TMDs, GPDs and PDFs.

We extract two nonsinglet nucleon Parton Distribution Functions from lattice QCD data for reduced Ioffe-time pseudodistributions. We perform such analysis within the NNPDF framework, considering data coming from different lattice ensembles and discussing in detail the treatment of the different source of systematics involved in the fit. We introduce a recipe for taking care of systematics and use it to perform our extraction of light-cone PDFs., Comment: 23 pages, 30 figures

The interactions between two octet baryons are studied at low energies using lattice QCD (LQCD) with larger-than-physical quark masses corresponding to a pion mass of $m_{\pi}\sim 450$ MeV and a kaon mass of $m_{K}\sim 596$ MeV. The two-baryon systems that are analyzed range from strangeness $S=0$ to $S=-4$ and include the spin-singlet and triplet $NN$, $\Sigma N$ ($I=3/2$), and $\Xi\Xi$ states, the spin-singlet $\Sigma\Sigma$ ($I=2$) and $\Xi\Sigma$ ($I=3/2$) states, and the spin-triplet $\Xi N$ ($I=0$) state. The $s$-wave scattering phase shifts, low-energy scattering parameters, and binding energies when applicable, are extracted using L\"uscher's formalism. While the results are consistent with most of the systems being bound at this pion mass, the interactions in the spin-triplet $\Sigma N$ and $\Xi\Xi$ channels are found to be repulsive and do not support bound states. Using results from previous studies at a larger pion mass, an extrapolation of the binding energies to the physical point is performed and is compared with experimental values and phenomenological predictions. The low-energy coefficients in pionless EFT relevant for two-baryon interactions, including those responsible for $SU(3)$ flavor-symmetry breaking, are constrained. The $SU(3)$ symmetry is observed to hold approximately at the chosen values of the quark masses, as well as the $SU(6)$ spin-flavor symmetry, predicted at large $N_c$. A remnant of an accidental $SU(16)$ symmetry found previously at a larger pion mass is further observed. The $SU(6)$-symmetric EFT constrained by these LQCD calculations is used to make predictions for two-baryon systems for which the low-energy scattering parameters could not be determined with LQCD directly in this study, and to constrain the coefficients of all leading $SU(3)$ flavor-symmetric interactions, demonstrating the predictive power of two-baryon EFTs matched to LQCD., Comment: 69 pages, 31 figures and 25 tables; published version

Extraction of hadronic observables at finite-momenta from Lattice QCD (LQCD) is constrained by the well-known signal-to-noise problems afflicting all such LQCD calculations. Traditional quark smearing algorithms are commonly used tools to improve the statistical quality of hadronic $n$-point functions, provided operator momenta are small. The momentum smearing algorithm of Bali et al. extends the range of momenta that are cleanly accessible, and has facilitated countless novel lattice calculations. Momentum smearing has, however, not been explicitly demonstrated within the framework of distillation. In this work we extend the momentum-smearing idea, by exploring a few modifications to the distillation framework. Together with enhanced time slice sampling and expanded operator bases engendered by distillation, we find ground-state nucleon energies can be extracted reliably for $\left|\vec{p}\right|\lesssim3\text{ GeV}$ and matrix elements featuring a large momentum dependence can be resolved., Comment: 19 pages, 30 figures

The fraction of the longitudinal momentum of ${}^3\text{He}$ that is carried by the isovector combination of $u$ and $d$ quarks is determined using lattice QCD for the first time. The ratio of this combination to that in the constituent nucleons is found to be consistent with unity at the few-percent level from calculations with quark masses corresponding to $m_\pi\sim 800$ MeV, extrapolated to the physical quark masses. This constraint is consistent with, and significantly more precise than, determinations from global nuclear parton distribution function fits. Including the lattice QCD determination of the momentum fraction in the nNNPDF global fitting framework results in the uncertainty on the isovector momentum fraction ratio being reduced by a factor of 2.5, and thereby enables a more precise extraction of the $u$ and $d$ parton distributions in ${}^3\text{He}$.

Over the last decade, numerical solutions of Quantum Chromodynamics (QCD) using the technique of lattice QCD have developed to a point where they are beginning to connect fundamental aspects of nuclear physics to the underlying degrees of freedom of the Standard Model. In this review, the progress of lattice QCD studies of nuclear matrix elements of electroweak currents and beyond-Standard-Model operators is summarized, and connections with effective field theories and nuclear models are outlined. Lattice QCD calculations of nuclear matrix elements can provide guidance for low-energy nuclear reactions in astrophysics, dark matter direct detection experiments, and experimental searches for violations of the symmetries of the Standard Model, including searches for additional CP violation in the hadronic and leptonic sectors, baryon-number violation, and lepton-number or flavor violation. Similarly, important inputs to neutrino experiments seeking to determine the neutrino-mass hierarchy and oscillation parameters, as well as other electroweak and beyond-Standard-Model processes can be determined. The phenomenological implications of existing studies of electroweak and beyond-Standard-Model matrix elements in light nuclear systems are discussed, and future prospects for the field toward precision studies of these matrix elements are outlined.

This report is an outcome of the workshop "AI for Nuclear Physics" held at Thomas Jefferson National Accelerator Facility on March 4-6, 2020. The workshop brought together 184 scientists to explore opportunities for Nuclear Physics in the area of Artificial Intelligence. The workshop consisted of plenary talks, as well as six working groups. The report includes the workshop deliberations and additional contributions to describe prospects for using AI across Nuclear Physics research., Comment: This version includes reference updates, improved figures and minor clarifications in the text

We report the results of a lattice quantum chromodynamics calculation of $F_K/F_\pi$ using M\"{o}bius domain-wall fermions computed on gradient-flowed $N_f=2+1+1$ highly-improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from $130 \lesssim m_\pi \lesssim 400$ MeV, four lattice spacings of $a\sim 0.15, 0.12, 0.09$ and $0.06$ fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the $a\sim0.06$ fm ensemble is helpful, but not necessary to achieve a subpercent determination of $F_K/F_\pi$. We also include an estimate of the strong isospin breaking corrections and arrive at a final result of $F_{K^\pm}/F_{\pi^\pm} = 1.1942(45)$ with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of $|V_{us}|/|V_{ud}| = 0.2311(10)$., Comment: v3: published version; v2: version submitted to journal; v1: 26 pages including 13 figures, appendices, and references. See https://github.com/callat-qcd/project_fkfpi for the analysis and data

We present results for the unpolarized parton distribution function of the nucleon computed in lattice QCD at the physical pion mass. This is the first study of its kind employing the method of Ioffe time pseudo-distributions. Beyond the reconstruction of the Bjorken-$x$ dependence we also extract the lowest moments of the distribution function using the small Ioffe time expansion of the Ioffe time pseudo-distribution. We compare our findings with the pertinent phenomenological determinations., Comment: v2, 8 pages, 6 figures, matches published version, to appear in Phys.Rev.Lett

We extract the pion valence quark distribution $q^\pi_{\rm v}(x)$ from lattice QCD (LQCD) calculated matrix elements of spacelike correlations of one vector and one axial vector current analyzed in terms of QCD collinear factorization, using a new short-distance matching coefficient calculated to one-loop accuracy. We derive the Ioffe time distribution of the two-current correlations in the physical limit by investigating the finite lattice spacing, volume, quark mass, and higher-twist dependencies in a simultaneous fit of matrix elements computed on four gauge ensembles. We find remarkable consistency between our extracted $q^\pi_{\rm v}(x)$ and that obtained from experimental data across the entire $x$-range. Further, we demonstrate that the one-loop matching coefficient relating the LQCD matrix computed in position space to the $q_{\rm v}^{\pi}(x)$ in momentum space has well-controlled behavior with Ioffe time. This justifies that LQCD calculated current-current correlations are good observables for extracting partonic structures by using QCD factorization, which complements to the global effort to extract partonic structure from experimental data., Comment: Published version, Physical Review D

The nucleon axial coupling, $g_A$, is a fundamental property of protons and neutrons, dictating the strength with which the weak axial current of the Standard Model couples to nucleons, and hence, the lifetime of a free neutron. The prominence of $g_A$ in nuclear physics has made it a benchmark quantity with which to calibrate lattice QCD calculations of nucleon structure and more complex calculations of electroweak matrix elements in one and few nucleon systems. There were a number of significant challenges in determining $g_A$, notably the notorious exponentially-bad signal-to-noise problem and the requirement for hundreds of thousands of stochastic samples, that rendered this goal more difficult to obtain than originally thought. I will describe the use of an unconventional computation method, coupled with "ludicrously'" fast GPU code, access to publicly available lattice QCD configurations from MILC and access to leadership computing that have allowed these challenges to be overcome resulting in a determination of $g_A$ with 1% precision and all sources of systematic uncertainty controlled. I will discuss the implications of these results for the convergence of $SU(2)$ Chiral Perturbation theory for nucleons, as well as prospects for further improvements to $g_A$ (sub-percent precision, for which we have preliminary results) which is part of a more comprehensive application of lattice QCD to nuclear physics. This is particularly exciting in light of the new CORAL supercomputers coming online, Sierra and Summit, for which our lattice QCD codes achieve a machine-to-machine speed up over Titan of an order of magnitude., Comment: Plenary presentation at The 9th International workshop on Chiral Dynamics

Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work summarizes our efforts to reduce the computational cost of computing the deflation space while achieving the desired variance reduction for Lattice QCD applications. Previous efforts computed the lower part of the singular spectrum of the Dirac operator by using an eigensolver preconditioned with a multigrid linear system solver. Despite the improvement in performance in those applications, as the problem size grows the runtime and storage demands of this approach will eventually dominate the stochastic estimation part of the computation. In this work, we propose to compute the deflation space in one of the following two ways. First, by using an inexact eigensolver on the Hermitian, but maximally indefinite, operator $A \gamma_5$. Second, by exploiting the fact that the multigrid prolongator for this operator is rich in components toward the lower part of the singular spectrum. We show experimentally that the inexact eigensolver can approximate the lower part of the spectrum even for ill-conditioned operators. Also, the deflation based on the multigrid prolongator is more efficient to compute and apply, and, despite its limited ability to approximate the fine level spectrum, it obtains similar variance reduction on the trace estimator as deflating with approximate eigenvectors from the fine level operator.

We present a calculation of the pion valence quark distribution extracted using the formalism of reduced Ioffe time pseudo-distributions or more commonly known as pseudo-PDFs. Our calculation is carried out on two different 2+1 flavor QCD ensembles using the isotropic-clover fermion action, with lattice dimensions $24^3\times 64$ and $32^3\times 96$ at the lattice spacing of $a=0.127$ fm, and with the quark mass equivalent to a pion mass of $m_\pi \simeq 415$ MeV. We incorporate several combinations of smeared-point and smeared-smeared pion source-sink interpolation fields in obtaining the lattice QCD matrix elements using the summation method. After one-loop perturbative matching and combining the pseudo-distributions from these two ensembles, we extract the pion valence quark distribution using a phenomenological functional form motivated by the global fits of parton distribution functions. We also calculate the lowest four moments of the pion quark distribution through the "OPE without OPE". We present a qualitative comparison between our lattice QCD extraction of the pion valence quark distribution with that obtained from global fits and previous lattice QCD calculations., Comment: Version accepted by Phys. Rev. D , 18 pages, 14 figures

In this paper, we present a detailed study of the unpolarized nucleon parton distribution function (PDF) employing the approach of parton pseudo-distribution functions. We perform a systematic analysis using three lattice ensembles at two volumes, with lattice spacings $a=$ 0.127 fm and $a=$ 0.094 fm, for a pion mass of roughly 400 MeV. With two lattice spacings and two volumes, both continuum limit and infinite volume extrapolation systematic errors of the PDF are estimated. In addition to the $x$ dependence of the PDF, we compute their first two moments and compare them with the pertinent phenomenological determinations., Comment: 1+46 pages, 21 figures

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD calculations related to the structure and spectroscopy of hadrons and nuclei. An overview of recent lattice calculations of the structure of the proton and other hadrons is presented along with prospects for future extensions. Progress and prospects of hadronic spectroscopy and the study of resonances in the light, strange and heavy quark sectors is summarized. Finally, recent advances in the study of light nuclei from lattice QCD are addressed, and the scope of future investigations that are currently envisioned is outlined., Comment: 45 pages

The computation of the parton distribution functions (PDF) or distribution amplitudes (DA) of hadrons from first principles lattice QCD constitutes a central open problem. In this study, we present and evaluate the efficiency of a selection of methods for inverse problems to reconstruct the full $x$-dependence of PDFs. Our starting point are the so called Ioffe time PDFs, which are accessible from Euclidean time calculations in conjunction with a matching procedure. Using realistic mock data tests, we find that the ill-posed incomplete Fourier transform underlying the reconstruction requires careful regularization, for which both the Bayesian approach as well as neural networks are efficient and flexible choices., Comment: 1+41 pages, 20 figures

We present the first exploratory lattice QCD calculation of the pion valence quark distribution extracted from spatially separated current-current correlations in coordinate space. We show that an antisymmetric combination of vector and axial-vector currents provides direct information on the pion valence quark distribution. Using the collinear factorization approach, we calculate the perturbative tree-level kernel for this current combination and extract the pion valence distribution. The main goal of this article is to demonstrate the efficacy of this general lattice QCD approach in the reliable extraction of parton distributions. With controllable power corrections and a good understanding of the lattice systematics, this method has the potential to serve as a complementary to the many efforts to extract parton distributions in global analyses from experimentally measured cross sections. We perform our calculation on an ensemble of 2+1 flavor QCD using the isotropic-clover fermion action, with lattice dimensions $32^3\times 96$ at a lattice spacing \mbox{$a=0.127$ fm} and the quark mass equivalent to a pion mass $m_\pi \simeq 416$ MeV., Comment: Version to appear in Physical Review D, 13 pages, 7 figures

The fundamental particle theory called Quantum Chromodynamics (QCD) dictates everything about protons and neutrons, from their intrinsic properties to interactions that bind them into atomic nuclei. Quantities that cannot be fully resolved through experiment, such as the neutron lifetime (whose precise value is important for the existence of light-atomic elements that make the sun shine and life possible), may be understood through numerical solutions to QCD. We directly solve QCD using Lattice Gauge Theory and calculate nuclear observables such as neutron lifetime. We have developed an improved algorithm that exponentially decreases the time-to solution and applied it on the new CORAL supercomputers, Sierra and Summit. We use run-time autotuning to distribute GPU resources, achieving 20% performance at low node count. We also developed optimal application mapping through a job manager, which allows CPU and GPU jobs to be interleaved, yielding 15% of peak performance when deployed across large fractions of CORAL., Comment: 2018 Gordon Bell Finalist: 9 pages, 9 figures; v2: fixed 2 typos and appended acknowledgements

We examine the relation of moments of parton distribution functions to matrix elements of non-local operators computed in lattice quantum chromodynamics. We argue that after the continuum limit is taken, these non-local matrix elements give access to moments that are finite and can be matched to those defined in the $\overline{MS}$ scheme. We demonstrate this fact with a numerical computation of moments through non-local matrix elements in the quenched approximation and we find that these moments are in excellent agreement with the moments obtained from direct computations of local twist-2 matrix elements in the quenched approximation., Comment: 1+11 pages, 1 figure, version to appear in JHEP

The $\textit{axial coupling of the nucleon}$, $g_A$, is the strength of its coupling to the $\textit{weak}$ axial current of the Standard Model of particle physics, in much the same way as the electric charge is the strength of the coupling to the electromagnetic current. This axial coupling dictates the rate at which neutrons decay to protons, the strength of the attractive long-range force between nucleons and other features of nuclear physics. Precision tests of the Standard Model in nuclear environments require a quantitative understanding of nuclear physics rooted in Quantum Chromodynamics, a pillar of the Standard Model. The prominence of $g_A$ makes it a benchmark quantity to determine theoretically - a difficult task because quantum chromodynamics is non-perturbative, precluding known analytical methods. Lattice Quantum Chromodynamics provides a rigorous, non-perturbative definition of quantum chromodynamics that can be implemented numerically. It has been estimated that a precision of two percent would be possible by 2020 if two challenges are overcome: contamination of $g_A$ from excited states must be controlled in the calculations and statistical precision must be improved markedly. Here we report a calculation of $g_A^{QCD} = 1.271\pm0.013$, using an unconventional method inspired by the Feynman-Hellmann theorem that overcomes these challenges., Comment: Published in Nature. 46 pages total: Main text 4 pages, Extended Data 8 pages, Supplemental 34 pages. Supporting data and code at https://github.com/callat-qcd/project_gA or https://zenodo.org/record/1241374

Complete flavour decompositions of the scalar, axial and tensor charges of the proton, deuteron, diproton and $^3$He at SU(3)-symmetric values of the quark masses corresponding to a pion mass $m_\pi\sim806$ MeV are determined using lattice QCD. At the physical quark masses, the scalar charges constrain mean-field models of nuclei and the low-energy interactions of nuclei with potential dark matter candidates. The axial and tensor charges of nuclei constrain their spin content, integrated transversity and the quark contributions to their electric dipole moments. External fields are used to directly access the quark-line connected matrix elements of quark bilinear operators, and a combination of stochastic estimation techniques is used to determine the disconnected sea-quark contributions. Significant nuclear modifications are found, with particularly large, O(10%), effects in the scalar charges. Typically, these nuclear effects reduce the effective charge of the nucleon (quenching), although in some cases an enhancement is not excluded. Given the size of the nuclear modifications of the scalar charges resolved here, contributions from correlated multi-nucleon effects should be quantified in the analysis of dark matter direct-detection experiments using nuclear targets., Comment: published version

In the framework of quantum chromodynamics (QCD), parton distribution functions (PDFs) quantify how the momentum and spin of a hadron are divided among its quark and gluon constituents. Two main approaches exist to determine PDFs. The first approach, based on QCD factorization theorems, realizes a QCD analysis of a suitable set of hard-scattering measurements, often using a variety of hadronic observables. The second approach, based on first-principle operator definitions of PDFs, uses lattice QCD to compute directly some PDF-related quantities, such as their moments. Motivated by recent progress in both approaches, in this document we present an overview of lattice-QCD and global-analysis techniques used to determine unpolarized and polarized proton PDFs and their moments. We provide benchmark numbers to validate present and future lattice-QCD calculations and we illustrate how they could be used to reduce the PDF uncertainties in current unpolarized and polarized global analyses. This document represents a first step towards establishing a common language between the two communities, to foster dialogue and to further improve our knowledge of PDFs., Comment: 80 pages, 15 figures, 28 tables, minor typos corrected, published in Progress in Particle and Nuclear Physics

Ioffe-time distributions, which are functions of the Ioffe-time $\nu$, are the Fourier transforms of parton distribution functions with respect to the momentum fraction variable $x$. These distributions can be obtained from suitable equal time, quark bilinear hadronic matrix elements which can be calculated from first principles in lattice QCD, as it has been recently argued. In this talk I present the first numerical calculation of the Ioffe-time distributions of the nucleon in the quenched approximation., Comment: 8 pages, 10 figures. arXiv admin note: text overlap with arXiv:1706.05373

We present state-of-the-art results from a lattice QCD calculation of the nucleon axial coupling, $g_A$, using M\"obius Domain-Wall fermions solved on the dynamical $N_f = 2 + 1 + 1$ HISQ ensembles after they are smeared using the gradient-flow algorithm. Relevant three-point correlation functions are calculated using a method inspired by the Feynman-Hellmann theorem, and demonstrate significant improvement in signal for fixed stochastic samples. The calculation is performed at five pion masses of $m_\pi\sim \{400, 350, 310, 220, 130\}$~MeV, three lattice spacings of $a\sim\{0.15, 0.12, 0.09\}$~fm, and we do a dedicated volume study with $m_\pi L\sim\{3.22, 4.29, 5.36\}$. Control over all relevant sources of systematic uncertainty are demonstrated and quantified. We achieve a preliminary value of $g_A = 1.285(17)$, with a relative uncertainty of 1.33\%., Comment: 18 pages, 8 figures, Lattice 2017 Proceedings

We present a new approach to extracting continuum quasi distributions from lattice QCD. Quasi distributions are defined by matrix elements of a Wilson-line operator extended in a spatial direction, evaluated between nucleon states at finite momentum. We propose smearing this extended operator with the gradient flow to render the corresponding matrix elements finite in the continuum limit. This procedure provides a nonperturbative method to remove the power-divergence associated with the Wilson line and the resulting matrix elements can be directly matched to light-front distributions via perturbation theory., Comment: Eight pages, two figures. Proceedings of the 35th International Symposium on Lattice Field Theory

The workshop on Parton Distributions and Lattice Calculations in the LHC era (PDFLattice2017) was hosted at Balliol College, Oxford (UK), from 22$^{\rm nd}$ to 24$^{\rm th}$ March 2017. The workshop brought together the lattice-QCD and the global-fit physicists who devote their efforts to determine the parton distribution functions (PDFs) of the proton. The goals were to make the two communities more familiar between each other, review developments from both sides, and set precision targets for lattice calculations so that they can contribute, together with the forthcoming experimental input, to the next generation of PDF determinations. This contribution summarises the relevant outcome of the workshop, in anticipation of a thorough white paper., Comment: 5 pages, 1 figure, contribution to the proceedings of the XXV Workshop on Deep-Inelastic Scattering and Related Subjects (DIS2017)

Magnetic moments of the octet baryons are computed using lattice QCD in background magnetic fields, including the first treatment of the magnetically coupled Sigma-Lambda system. Although the computations are performed for relatively large values of the up and down quark masses, we gain new insight into the symmetries and relations between magnetic moments by working at a three-flavor mass-symmetric point. While the spin-flavor symmetry in the large Nc limit of QCD is shared by the naive constituent quark model, we find instances where quark model predictions are considerably favored over those emerging in the large Nc limit. We suggest further calculations that would shed light on the curious patterns of baryon magnetic moments., Comment: 8 pages, 5 figures, talk given at Lattice 2017, the 35th International Symposium on Lattice Field Theory, Granada, Spain, 18-24 June 2017

The role of gluons in the structure of the nucleon and light nuclei is investigated using lattice quantum chromodynamics (QCD) calculations. The first moment of the unpolarised gluon distribution is studied in nuclei up to atomic number $A=3$ at quark masses corresponding to pion masses of $m_\pi\sim 450$ and $806$ MeV. Nuclear modification of this quantity defines a gluonic analogue of the EMC effect and is constrained to be less than $\sim 10$% in these nuclei. This is consistent with expectations from phenomenological quark distributions and the momentum sum rule. In the deuteron, the combination of gluon distributions corresponding to the $b_1$ structure function is found to have a small first moment compared with the corresponding momentum fraction. The first moment of the gluon transversity structure function is also investigated in the spin-1 deuteron, where a non-zero signal is observed at $m_\pi \sim 806$ MeV. This is the first indication of gluon contributions to nuclear structure that can not be associated with an individual nucleon.

Lattice quantum chromodynamics is used to constrain the interactions of two octet baryons at the SU(3) flavor-symmetric point, with quark masses that are heavier than those in nature (equal to that of the physical strange quark mass and corresponding to a pion mass of $\approx 806~\tt{MeV}$). Specifically, the S-wave scattering phase shifts of two-baryon systems at low energies are obtained with the application of L\"uscher's formalism, mapping the energy eigenvalues of two interacting baryons in a finite volume to the two-particle scattering amplitudes below the relevant inelastic thresholds. The values of the leading-order low-energy scattering parameters in the irreducible representations of SU(3) are consistent with an approximate SU(6) spin-flavor symmetry in the nuclear and hypernuclear forces that is predicted in the large-$N_c$ limit of QCD. The two distinct SU(6)-invariant interactions between two baryons are constrained at this value of the quark masses, and their values indicate an approximate accidental SU(16) symmetry. The SU(3) irreducible representations containing the $NN~({^1}S_0)$, $NN~({^3}S_1)$ and $\frac{1}{\sqrt{2}}(\Xi^0n+\Xi^-p)~({^3}S_1)$ channels unambiguously exhibit a single bound state, while the irreducible representation containing the $\Sigma^+ p~({^3}S_1)$ channel exhibits a state that is consistent with either a bound state or a scattering state close to threshold. These results are in agreement with the previous conclusions of the NPLQCD collaboration regarding the existence of two-nucleon bound states at this value of the quark masses., Comment: 49 pages, 21 figures, 14 tables. v2: Minor wording improvements and updated references

We demonstrate a new method of extracting parton distributions from lattice calculations. The starting idea is to treat the generic equal-time matrix element ${\cal M} (Pz_3, z_3^2)$ as a function of the Ioffe time $\nu = Pz_3$ and the distance $z_3$. The next step is to divide ${\cal M} (Pz_3, z_3^2)$ by the rest-frame density ${\cal M} (0, z_3^2)$. Our lattice calculation shows a linear exponential $z_3$-dependence in the rest-frame function, expected from the $Z(z_3^2)$ factor generated by the gauge link. Still, we observe that the ratio ${\cal M} (Pz_3 , z_3^2)/{\cal M} (0, z_3^2)$ has a Gaussian-type behavior with respect to $z_3$ for 6 values of $P$ used in the calculation. This means that $Z(z_3^2)$ factor was canceled in the ratio. When plotted as a function of $\nu$ and $z_3$, the data are very close to $z_3$-independent functions. This phenomenon corresponds to factorization of the $x$- and $k_\perp$-dependence for the TMD ${\cal F} (x, k_\perp^2)$. For small $z_3 \leq 4a$, the residual $z_3$-dependence is explained by perturbative evolution, with $\alpha_s/\pi =0.1$., Comment: 12 pages, 18 figures, to appear in Phys. Rev. D

In this comment, we address a number of erroneous discussions and conclusions presented in a recent preprint by the HALQCD collaboration, arXiv:1703.07210. In particular, we demonstrate that lattice QCD determinations of bound states at quark masses corresponding to a pion mass of $m_\pi = 806$ MeV are robust, and that the phases shifts extracted by the NPLQCD collaboration for these systems pass all of the 'sanity checks' introduced in arXiv:1703.07210., Comment: Clarifications added

We report on a lattice QCD calculation of the nucleon axial charge, $g_A$, using M\"{o}bius Domain-Wall fermions solved on the dynamical $N_f=2+1+1$ HISQ ensembles after they are smeared using the gradient-flow algorithm. The calculation is performed with three pion masses, $m_\pi\sim\{310,220,130\}$ MeV. Three lattice spacings ($a\sim\{0.15,0.12,0.09\}$ fm) are used with the heaviest pion mass, while the coarsest two spacings are used on the middle pion mass and only the coarsest spacing is used with the near physical pion mass. On the $m_\pi\sim220$ MeV, $a\sim0.12$ fm point, a dedicated volume study is performed with $m_\pi L \sim \{3.22,4.29,5.36\}$. Using a new strategy motivated by the Feynman-Hellmann Theorem, we achieve a precise determination of $g_A$ with relatively low statistics, and demonstrable control over the excited state, continuum, infinite volume and chiral extrapolation systematic uncertainties, the latter of which remains the dominant uncertainty. Our final determination at 2.6\% total uncertainty is $g_A = 1.278(21)(26)$, with the first uncertainty including statistical and systematic uncertainties from fitting and the second including model selection systematics related to the chiral and continuum extrapolation. The largest reduction of the second uncertainty will come from a greater number of pion mass points as well as more precise lattice QCD results near the physical pion mass., Comment: 17 pages + 11 pages of references and appendices. 15 figures. Interested readers can download the Python analysis scripts and an hdf5 data file at https://github.com/callat-qcd/project_gA_v0