Search

Studies of the $\Delta$ baryon resonance and the $K_0^\ast(700)$ and $a_0(980)$ meson resonances using $N_f=2+1$ lattice QCD for pion masses near 200 MeV are presented. The $s$-wave scattering lengths for both the $I=1/2$ $N \pi$ and $I=3/2$ $N \pi$ channels and properties of the $\Delta$ resonance are identified from the finite-volume energy levels of the lattice simulation. The importance of a three-quark $\Delta$-operator in the $N\pi$ system and tetraquark operators in the mesonic systems is investigated., Comment: 10 pages, 5 figures, 3 tables

This work presents technical details of determining the finite-volume energy spectra for the scattering amplitude of the coupled-channel $\pi\Sigma - \bar{K}N$ from lattice QCD data. The importance of reliably extracting such spectra lies in the crucial dependence of the hadronic scattering amplitudes analysis on the energy spectrum when using L\"{u}scher's formalism. Results of the methods used are presented and the final finite-volume spectra are shown. The analysis of the scattering amplitude based on these results, exhibits a two-pole structure for the $\Lambda(1405)$, a virtual bound state below the $\pi\Sigma$ threshold and a resonance pole right below the $\bar{K}N$ threshold.

Recent results studying the masses and widths of low-lying baryon resonances in lattice QCD are presented. The $S$-wave $N\pi$ scattering lengths for both total isospins $I = 1/2$ and $I = 3/2$ are inferred from the finite-volume spectrum below the inelastic threshold together with the $I = 3/2$ $P$-wave containing the $\Delta(1232)$ resonance. A lattice QCD computation employing a combined basis of three-quark and meson-baryon interpolating operators with definite momentum to determine the coupled channel $\Sigma\pi$-$N\overline{K}$ scattering amplitude in the $\Lambda(1405)$ region is also presented. Our results support the picture of a two-pole structure suggested by theoretical approaches based on $SU(3)$ chiral symmetry and unitarity., Comment: 5 pages, 3 figures, talk presented at the 20th International Conference on Hadron Spectroscopy and Structure (HADRON 2023) in Genova, Italy, June 5th to 9th 2023

A lattice QCD computation of the coupled channel $\pi\Sigma-\bar{K}N$ scattering amplitudes in the $\Lambda(1405)$ region is detailed. Results are obtained using a single ensemble of gauge field configurations with $N_{\rm f} = 2+1$ dynamical quark flavors and $m_{\pi} \approx 200$ MeV and $m_K\approx487$ MeV. Hermitian correlation matrices using both single baryon and meson-baryon interpolating operators for a variety of different total momenta and irreducible representations are used. Several parametrizations of the two-channel scattering $K$-matrix are utilized to obtain the scattering amplitudes from the finite-volume spectrum. The amplitudes, continued to the complex energy plane, exhibit a virtual bound state below the $\pi\Sigma$ threshold and a resonance pole just below the $\bar{K}N$ threshold., Comment: 18 pages, 11 figures, 14 tables. Final version accepted for publication in PRD with post-publication corrections

This letter presents the first lattice QCD computation of the coupled channel $\pi\Sigma-\bar{K}N$ scattering amplitudes at energies near $1405\,{\rm MeV}$. These amplitudes contain the resonance $\Lambda(1405)$ with strangeness $S=-1$ and isospin, spin, and parity quantum numbers $I(J^P)=0(1/2^-)$. However, whether there is a single resonance or two nearby resonance poles in this region is controversial theoretically and experimentally. Using single-baryon and meson-baryon operators to extract the finite-volume stationary-state energies to obtain the scattering amplitudes at slightly unphysical quark masses corresponding to $m_\pi\approx200$ MeV and $m_K\approx487$ MeV, this study finds the amplitudes exhibit a virtual bound state below the $\pi\Sigma$ threshold in addition to the established resonance pole just below the $\bar{K}N$ threshold. Several parametrizations of the two-channel $K$-matrix are employed to fit the lattice QCD results, all of which support the two-pole picture suggested by $SU(3)$ chiral symmetry and unitarity., Comment: 7 pages, 4 figures, and 1 table. Includes various corrections made during the last stages before publication

Lattice gauge theory continues to be a powerful theoretical and computational approach to simulating strongly interacting quantum field theories, whose applications permeate almost all disciplines of modern-day research in High-Energy Physics. Whether it is to enable precision quark- and lepton-flavor physics, to uncover signals of new physics in nucleons and nuclei, to elucidate hadron structure and spectrum, to serve as a numerical laboratory to reach beyond the Standard Model, or to invent and improve state-of-the-art computational paradigms, the lattice-gauge-theory program is in a prime position to impact the course of developments and enhance discovery potential of a vibrant experimental program in High-Energy Physics over the coming decade. This projection is based on abundant successful results that have emerged using lattice gauge theory over the years: on continued improvement in theoretical frameworks and algorithmic suits; on the forthcoming transition into the exascale era of high-performance computing; and on a skillful, dedicated, and organized community of lattice gauge theorists in the U.S. and worldwide. The prospects of this effort in pushing the frontiers of research in High-Energy Physics have recently been studied within the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021), and the conclusions are summarized in this Topical Report., Comment: 57 pages, 1 figure. Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Topical Group Report for TF05 - Lattice Gauge Theory

Elastic nucleon-pion scattering amplitudes are computed using lattice QCD on a single ensemble of gauge field configurations with $N_{\rm f} = 2+1$ dynamical quark flavors and $m_{\pi} = 200~{\rm MeV}$. The $s$-wave scattering lengths with both total isospins $I=1/2$ and $I=3/2$ are inferred from the finite-volume spectrum below the inelastic threshold together with the $I=3/2$ $p$-wave containing the $\Delta(1232)$ resonance. The amplitudes are well-described by the effective range expansion with parameters constrained by fits to the finite-volume energy levels enabling a determination of the $I=3/2$ scattering length with statistical errors below $5\%$, while the $I = 1/2$ is somewhat less precise. Systematic errors due to excited states and the influence of higher partial waves are controlled, providing a pathway for future computations down to the physical light quark masses with multiple lattice spacings and physical volumes., Comment: 33 pages, 16 figures; version accepted for publication

The observation of neutrino oscillations and hence non-zero neutrino masses provided a milestone in the search for physics beyond the Standard Model. But even though we now know that neutrinos are massive, the nature of neutrino masses, i.e., whether they are Dirac or Majorana, remains an open question. A smoking-gun signature of Majorana neutrinos is the observation of neutrinoless double-beta decay, a process that violates the lepton-number conservation of the Standard Model. This white paper focuses on the theoretical aspects of the neutrinoless double-beta decay program and lays out a roadmap for future developments. The roadmap is a multi-scale path starting from high-energy models of neutrinoless double-beta decay all the way to the low-energy nuclear many-body problem that needs to be solved to supplement measurements of the decay rate. The path goes through a systematic effective-field-theory description of the underlying processes at various scales and needs to be supplemented by lattice quantum chromodynamics input. The white paper also discusses the interplay between neutrinoless double-beta decay, experiments at the Large Hadron Collider and results from astrophysics and cosmology in probing simplified models of lepton-number violation at the TeV scale, and the generation of the matter-antimatter asymmetry via leptogenesis. This white paper is prepared for the topical groups TF11 (Theory of Neutrino Physics), TF05 (Lattice Gauge Theory), RF04 (Baryon and Lepton Number Violating Processes), NF03 (Beyond the Standard Model) and NF05 (Neutrino Properties) within the Theory Frontier, Rare Processes and Precision Frontier, and Neutrino Physics Frontier of the U.S. Community Study on the Future of Particle Physics (Snowmass 2021)., Comment: Contribution to Snowmass 2021

This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges ahead, and new ideas and strategies to make progress in nuclear structure and reaction physics, effective field theory, lattice QCD, quantum information, and quantum computing. The contributed pieces solely reflect the perspectives of the respective authors and do not represent the viewpoints of the Institute for Nuclear theory or the organizers of the program., Comment: Perspective pieces of the virtual INT program 21-1b "Nuclear Forces for Precision Nuclear Physics", 107 pages

Hyperon decays present a promising alternative for extracting $\vert V_{us} \vert$ from lattice QCD combined with experimental measurements. Currently $\vert V_{us} \vert$ is determined from the kaon decay widths and a lattice calculation of the associated form factor. In this proceeding, I will present preliminary work on a lattice determination of the hyperon mass spectrum. I will additionally summarize future goals in which we will calculate the hyperon transition matrix elements, which will provide an alternative means for accessing $\vert V_{us} \vert$. This work is based on a particular formulation of SU(2) chiral perturbation theory for hyperons; determining the extent to which this effective field theory converges is instrumental in understanding the limits of its predictive power, especially since some hyperonic observables are difficult to calculate near the physical pion mass (e.g., hyperon-to-nucleon form factors), and thus the use of heavier than physical pion masses is likely to yield more precise results when combined with extrapolations to the physical point.}, Comment: 8 pages, 3 figures, presented at the 38th International Symposium on Lattice Field Theory (LATTICE2021)

Lattice QCD calculations of two-nucleon interactions have been underway for about a decade, but still haven't reached the pion mass regime necessary for matching onto effective field theories and extrapolating to the physical point. Furthermore, results from different methods, including the use of the Luscher formalism with different types of operators, as well as the HALQCD potential method, do not agree even qualitatively at very heavy pion mass. We investigate the role that different operators employed in the literature may play on the extraction of spectra for use within the Luscher method. We first explore expectations from Effective Field Theory solved within a finite volume, for which the exact spectrum may be computed given different physical scenarios. We then present preliminary lattice QCD results for two-nucleon spectra calculated using different operators on a common lattice ensemble., Comment: 11 pages, 7 figures, proceeding for The 38th International Symposium on Lattice Field Theory, LATTICE2021, 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Progress in computing various meson-baryon scattering amplitudes is presented on a single ensemble from the Coordinated Lattice Simulations (CLS) consortium with $m_\pi=200$ MeV and $N_f=2+1$ dynamical fermions. The finite-volume L\"uscher approach is employed to determine the lowest few partial waves from ground- and excited-state energies computed from correlation matrices rotated in a single pivot using a generalized eigenvector solution. This analysis requires evaluating matrices of correlation functions between single- and two-hadron interpolating operators which are projected onto definite spatial momenta and finite-volume irreducible representations. The stochastic LapH method is used to estimate all needed quark propagators. Preliminary results are presented for $I=\frac{1}{2}, \frac{3}{2}$ $N\pi$ amplitudes including the $\Delta(1232)$ resonance and the $I=0$ $S$-wave amplitude with unit strangeness relevant for the $\Lambda(1405)$., Comment: The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology 9 pages, 4 figures

The Deep Underground Neutrino Experiment (DUNE) is an upcoming neutrino oscillation experiment that is poised to answer key questions about the nature of neutrinos. Lattice QCD has the ability to make significant impact upon DUNE, beginning with computations of nucleon-neutrino interactions with weak currents. Nucleon amplitudes involving the axial form factor are part of the primary signal measurement process for DUNE, and precise calculations from LQCD can significantly reduce the uncertainty for inputs into Monte Carlo generators. Recent calculations of the nucleon axial charge have demonstrated that sub-percent precision is possible on this vital quantity. In these proceedings, we discuss preliminary results for the CalLat collaboration's calculation of the axial form factor of the nucleon. These computations are performed with M\"obius domain wall valence quarks on HISQ sea quark ensembles generated by the MILC and CalLat collaborations. The results use a variety of ensembles including several at physical pion mass., Comment: The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology 9 pages, 2 figures

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: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations $\gtrsim1$ fm are commonly used to combat these issues. In this work, focusing on $g_A$, we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2--1 fm), combined with a multi-state analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis 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., Comment: v2: updates based on referee comments and some community response, consistent with published version; v1: 13 pages plus appendices. The correlation function data and analysis code accompanying this publication can be accessed at this github repository: https://github.com/callat-qcd/project_fh_vs_3pt

We report on a sub-percent scale determination using the omega baryon mass and gradient-flow methods. The calculations are performed on 22 ensembles of $N_f=2+1+1$ highly improved, rooted staggered sea-quark configurations generated by the MILC and CalLat Collaborations. The valence quark action used is M\"obius Domain-Wall fermions solved on these configurations after a gradient-flow smearing is applied with a flowtime of $t_{\rm gf}=1$ in lattice units. The ensembles span four lattice spacings in the range $0.06 \lesssim a \lesssim 0.15$ fm, six pion masses in the range $130 \lesssim m_\pi \lesssim 400$ MeV and multiple lattice volumes. On each ensemble, the gradient-flow scales $t_0/a^2$ and $w_0/a$ and the omega baryon mass $a m_\Omega$ are computed. The dimensionless product of these quantities is then extrapolated to the continuum and infinite volume limits and interpolated to the physical light, strange and charm quark mass point in the isospin limit, resulting in the determination of $\sqrt{t_0}=0.1422(14)$ fm and $w_0 = 0.1709(11)$ fm with all sources of statistical and systematic uncertainty accounted for. The dominant uncertainty in this result is the stochastic uncertainty, providing a clear path for a few-per-mille uncertainty, as recently obtained by the Budapest-Marseille-Wuppertal Collaboration., Comment: v3: Published version; v2: Added determination of t_0 as well as w_0; v1: 13 pages plus appendices. The correlation function data, mass results and analysis code accompanying this publication can be found at this github repository: https://github.com/callat-qcd/project_scale_setting_mdwf_hisq

We report on the first application of the stochastic Laplacian Heaviside method for computing multi-particle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a positive definite set of two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance in order to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the $SU(3)$-flavor symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of $a\sim0.086$~fm, lattice volume of $L=48a\simeq4.1$~fm and pion mass $m_\pi\simeq714$ MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel., Comment: v2: version to be published in Phys. Rev. C.; v1: 13 pages plus figures and appendices

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

Second order beta-decay processes with and without neutrinos in the final state are key probes of nuclear physics and of the nature of neutrinos. Neutrinoful double-beta decay is the rarest Standard Model process that has been observed and provides a unique test of the understanding of weak nuclear interactions. Observation of neutrinoless double-beta decay would reveal that neutrinos are Majorana fermions and that lepton number conservation is violated in nature. While significant progress has been made in phenomenological approaches to understanding these processes, establishing a connection between these processes and the physics of the Standard Model and beyond is a critical task as it will provide input into the design and interpretation of future experiments. The strong-interaction contributions to double-beta decay processes are non-perturbative and can only be addressed systematically through a combination of lattice Quantum Chromoodynamics (LQCD) and nuclear many-body calculations. In this review, current efforts to establish the LQCD connection are discussed for both neutrinoful and neutrinoless double-beta decay. LQCD calculations of the hadronic contributions to the neutrinoful process $nn\to pp e^- e^- \bar\nu_e\bar\nu_e$ and to various neutrinoless pionic transitions are reviewed, and the connections of these calculations to the phenomenology of double-beta decay through the use of effective field theory (EFTs) is highlighted. At present, LQCD calculations are limited to small nuclear systems, and to pionic subsystems, and require matching to appropriate EFTs to have direct phenomenological impact. However, these calculations have already revealed qualitatively that there are terms in the EFTs that can only be constrained from double-beta decay processes themselves or using inputs from LQCD. Future prospects for direct calculations in larger nuclei are also discussed., Comment: 53 pages, 17 figures

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

Exascale computing could soon enable a predictive theory of nuclear structure and reactions rooted in the Standard Model, with quantifiable and systematically improvable uncertainties. Such a predictive theory will help exploit experiments that use nucleons and nuclei as laboratories for testing the Standard Model and its limitations. Examples include direct dark matter detection, neutrinoless double beta decay, and searches for permanent electric dipole moments of the neutron and atoms. It will also help connect QCD to the properties of cold neutron stars and hot supernova cores. We discuss how a quantitative bridge between QCD and the properties of nuclei and nuclear matter will require a synthesis of lattice QCD (especially as applied to two- and three-nucleon interactions), effective field theory, and ab initio methods for solving the nuclear many-body problem. While there are significant challenges that must be addressed in developing this triad of theoretical tools, the rapid advance of computing is accelerating progress. In particular, we focus this review on the anticipated advances from lattice QCD and how these advances will impact few-body effective theories of nuclear physics by providing critical input, such as constraints on unknown low-energy constants of the effective (field) theories. We also review particular challenges that must be overcome for the successful application of lattice QCD for low-energy nuclear physics. We describe progress in developing few-body effective (field) theories of nuclear physics, with an emphasis on HOBET, a non-relativistic effective theory of nuclear physics, which is less common in the literature. We use the examples of neutrinoless double beta decay and the nuclear-matter equation of state to illustrate how the coupling of lattice QCD to effective theory might impact our understanding of symmetries and exotic astrophysical environments., Comment: v3: version published in Progress in Particle and Nuclear Physics; v2: updated manuscript based upon community feedback and referee comments. Also, substantially re-written section on two-nucleon lattice QCD controversy. 53.5 pages plus a "few more" references; v1: Contribution to: The tower of effective (field) theories and the emergence of nuclear phenomena; 47 pages plus a "few" references

The Feynman-Hellmann method, as implemented by Bouchard et al. [1612.06963], was recently employed successfully to determine the nucleon axial charge. A limitation of the method was the restriction to a single operator and a single momentum during the computation of each "Feynman- Hellmann" propagator. By using stochastic techniques to estimate the all-to-all propagator, we relax this constraint and demonstrate the successful implementation of this new method. We show reproduction of the axial charge on a test ensemble and non-zero momentum transfer points of the axial and vector form factors., Comment: 7 pages, 5 figures, Proceedings for The 36th International Symposium on Lattice Field Theory

The search for neutrinoless double beta decay of nuclei is believed to be one of the most promising means to search for new physics. Observation of this very rare nuclear process, which violates Lepton Number conservation, would imply the neutrino sector has a Majorana mass component and may also provide an explanation for the universe matter-antimatter asymmetry of the universe. In the case where a heavy intermediate particle is exchanged in this process, QCD contributions from short range interactions become relevant and the calculation of matrix elements with four-quark operators becomes necessary. In these proceedings we will discuss our current progress in the calculation of these four-quark operators from LQCD., Comment: Proceedings of the 36th Annual International Symposium on Lattice Field Theory, 22-28 July, 2018, Michigan State University, East Lansing, Michigan

Anchoring the nuclear interaction in QCD is a long-outstanding problem in nuclear physics. While the lattice community has made enormous progress in mesonic physics and single nucleon physics, continuum-limit physical-point multi-nucleon physics has remained out of reach. I will review CalLat's strategy for multi-nucleon spectroscopy and our latest results., Comment: Proceedings of the 36th Annual International Symposium on Lattice Field Theory, 22-28 July, 2018, Michigan State University, East Lansing, Michigan

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

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