No Arabic abstract
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.
Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors. Most lattice QCD collaborations advocate for the use of high-statistics calculations at large time separations ($t_{rm sep}gtrsim1$ fm) to combat the signal-to-noise degradation. In this work we demonstrate that, for the nucleon axial charge, $g_A$, the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations ($0.2lesssim t_{rm sep}lesssim1$ fm), combined with a multi-state analysis, provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty, including correlated late-time fluctuations that may bias the ground state. We show that two classes of excited states largely cancel in the ratio of the three-point to two-point functions, leaving the third class, the transition matrix elements, as the dominant source of contamination. On an $m_piapprox310$ MeV ensemble, we observe the expected exponential suppression of excited state contamination in the Feynman-Hellmann correlation function relative to the standard three-point function; the excited states of the regular three-point function reduce to the 1% level for $t_{rm sep} >2$ fm while, for the Feynman-Hellmann correlation function, they are suppressed to 1% at $t_{rm sep}approx1$ fm. Independent analyses of the three-point and Feynman-Hellmann correlators yield consistent results for the ground state. However, a combined analysis allows for a more detailed and robust understanding of the excited state contamination, improving the demonstration that the ground state parameters are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.
We present a lattice QCD calculation of the $Delta(1232)$ matrix elements of the axial-vector and pseudoscalar currents. The decomposition of these matrix elements into the appropriate Lorentz invariant form factors is carried out and the techniques to calculate the form factors are developed and tested using quenched configurations. Results are obtained for 2+1 domain wall fermions and within a hybrid scheme with domain wall valence and staggered sea quarks. Two Goldberger-Treiman type relations connecting the axial to the pseudoscalar effective couplings are derived. These and further relations based on the pion-pole dominance hypothesis are examined using the lattice QCD results, finding support for their validity. Utilizing lattice QCD results on the axial charges of the nucleon and the $Delta$, as well as the nucleon-to-$Delta$ transition coupling constant, we perform a combined chiral fit to all three quantities and study their pion mass dependence as the chiral limit is approached.
Matrix elements of six-quark operators are needed to extract new physics constraints from experimental searches for neutron-antineutron oscillations. This work presents in detail the first lattice quantum chromodynamics calculations of the necessary neutron-antineutron transition matrix elements including calculation methods and discussions of systematic uncertainties. Implications of isospin and chiral symmetry on the matrix elements, power counting in the isospin limit, and renormalization of a chiral basis of six-quark operators are discussed. Calculations are performed with a chiral-symmetric discretization of the quark action and physical light quark masses in order to avoid the need for chiral extrapolation. Non-perturbative renormalization is performed, including a study of lattice cutoff effects. Excited-state effects are studied using two nucleon operators and multiple values of source-sink separation. Results for the dominant matrix elements are found to be significantly larger compared to previous results from the MIT bag model. Future calculations are needed to fully account for systematic uncertainties associated with discretization and finite-volume effects but are not expected to significantly affect this conclusion.
Recent results from lattice QCD simulations provide a realistic picture, based upon first principles, of~$Upsilon$ physics. We combine these results with the experimentally measured mass of the $Upsilon$~meson to obtain an accurate and reliable value for the $b$-quarks pole mass. We use two different methods, each of which yields a mass consistent with $M_b = 5.0(2)$~GeV. This corresponds to a bare mass of $M_b^0 = 4.0(1)$~GeV in our lattice theory and an $msbar$~mass of $M_b^msbar(M_b)=4.0(1)$~GeV. We discuss the implications of this result for the $c$-quark mass. ******************************************************************************* THIS IS THE VERSION WHICH WILL BE PUBLISHED IN PRL. SUBSTANTIAL MATERIAL HAS BEEN ADDED, INCLUDING RESULTS WITH DYNAMICAL FERMIONS AND A CALCULATION OF THE MSBAR MASS. *******************************************************************************
We present lattice results for the non-perturbative Collins-Soper (CS) kernel, which describes the energy-dependence of transverse momentum-dependent parton distributions (TMDs). The CS kernel is extracted from the ratios of first Mellin moments of quasi-TMDs evaluated at different nucleon momenta.The analysis is done with dynamical $N_f=2+1$ clover fermions for the CLS ensemble H101 ($a=0.0854,mathrm{fm}$, $m_{pi}=m_K=422,mathrm{MeV}$). The computed CS kernel is in good agreement with experimental extractions and previous lattice studies.