# Lattice QCD at Los Alamos

Lattice QCD (LQCD) is a non-perturbative formulation of Quantum Chromodynamics (QCD), the regnant theory of strong interactions, on a discrete Euclidian space-time grid. Large scale simulations of lattice QCD allow us to calculate the contributions of QCD to the properties, decays and interactions of hadrons composed of quarks and gluons. With the development of highly efficient numerical algorithms and peta-scale computing resources, we are able to calculate, from first principles, many quantities that are being used to both validate QCD and search for new physics beyond the standard model of elementary particles and interactions.

These simulations are extremely computationally intensive and require leadership class computing resources. They are being done at national supercomputer centers at NERSC and Oak Ridge, USQCD clusters at Fermilab and JLab, and institutional computing resources at Los Alamos.

The Los Alamos lattice QCD group has made many pioneering contributions since its inception in 1984. Current staff members include Tanmoy Bhattacharya, Vicenzo Cirigliano, Michael Graesser, Rajan Gupta, Emanuele Mereghetti, and Boram Yoon. They are working with postdoctoral fellows Yong-Chull Jang, Santanu Mondal and Sungwoo Park on five projects:

- Calculations of matrix elements of isovector quark bilinear operators within the neutron and proton states are being carried out to explore signatures of novel scalar and tensor interactions that can arise at the TeV scale. By combining the results for the matrix elements with precision measurements of [ultra]cold neutron decays, one can put constraints on the possible size of these interactions. (see Figure 1) The PNDME (Precision Neutron Decay Matrix Elements) collaboration motivated these calculations and has published a number of state-of-the-art results that have recently been reviewed in the FLAG 2019 report and included in the community averages. Published papers include Probing Novel Scalar and Tensor Interactions from (Ultra)Cold Neutrons to the LHC; Axial, scalar, and tensor charges of the nucleon from 2 + 1 + 1 -flavor Lattice QCD; Isovector charges of the nucleon from 2 + 1 + 1 -flavor lattice QCD
- Calculations of the axial vector, electric and magnetic form factors that are needed to calculate the cross-section of the scattering of neutrinos, electrons and muons off nuclear targets. Our recent work has resolved the violation of the PCAC relation between the form factors. (See Figure 2) Published papers include Axial-vector form factors of the nucleon from lattice QCD; Axial Vector Form Factors from Lattice QCD that Satisfy the PCAC Relation; Nucleon Electromagnetic Form Factors in the Continuum Limit from 2+1+1-flavor Lattice QCD
- Calculations of matrix elements of novel CP violating operators that arise at the TeV scale. Each CP violating interaction makes a contribution to the neutron electric dipole moment (nEDM), and in many BSM theories these are large enough to explain the observed matter-antimatter asymmetry of the universe generated via baryogenesis. Using the tools of effective field theory, these operators are written in terms of quark and gluon fields at the hadronic (~1 GeV) scale. Lattice QCD is then used to calculate their matrix elements within the neutron state. Their contribution to the nEDM is given by
*d*where_{n}= ∑_{i}ε_{i}g_{i}*ε*is the strength (coupling at ~1 GeV scale) of the_{i}*i*operator and^{th}*g*is its matrix element. Our goal is to provide high precision results for_{i}*g*for all operators that make a significant contribution and use future bounds on nEDM to constrain new CP violation, ie, the . (See Figure 3 for current constraints on the split SUSY model.) The operators under investigated include the Q-term, the quark EDM, the quark chromo EDM and Weinberg three gluon operators. Published papers include Neutron Electric Dipole Moment from Beyond the Standard Model; Dimension-5 CP-odd operators: QCD mixing and renormalization; Neutron Electric Dipole Moment and Tensor Charges from Lattice QCD_{i} - Calculations of the matrix elements of flavor diagonal operators with axial, scalar, tensor and vector Lorentz structures. These provide (i) how the spin of the nucleons arises from the spin of the quarks and gluons and their orbital angular momentum; (ii) the nucleon sigma term
*σ*and its strangeness content_{πN}*σ*; (iii) the strength of the interaction of dark matter with nuclear targets; and (iv) the average momentum, transversity and helicity distributions of each quark flavor within nucleons. Published papers include Quark contribution to the proton spin from 2 + 1 + 1 -flavor lattice QCD; Flavor diagonal tensor charges of the nucleon from ( 2 + 1 + 1 )-flavor lattice QCD_{s} - Calculations of matrix elements of light cone operators within the proton state that shed light on the contribution of the angular momentum of quarks to the proton spin. Of particular interest is the extraction of the Sivers and Boer-Mulders distribution functions. (See Figure 5) These calculations are being done in collaboration with Michael Engelhardt (NMSU), and Andreas Schaefer (University of Regensburg, Germany). Published papers include Nucleon transverse momentum-dependent parton distributions in lattice QCD: Renormalization patterns and discretization effects

During 2007—2014, Gupta and Bhattacharya co-led the US wide HotQCD collaboration and investigated the behavior of QCD at finite temperature to elucidate the nature of the chiral transition to the quark-gluon plasma, determine the equation of state, and examine the restoration of U(1) axial symmetry. These results were published in six papers:

(i) The equation of state in (2+1)-flavor QCD;

(ii) QCD Phase Transition with Chiral Quarks and Physical Quark Masses;

(iii) Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model;

(iv) Chiral transition and U(1)_A symmetry restoration from lattice QCD using domain wall fermions;

(v) Chiral and deconfinement aspects of the QCD transition__;__ and

(vi) Equation of state and QCD transition at finite temperature.

Their work showed that, for physical values of the up, down and strange quark masses, the transition from the hadron phase to the quark-gluon plasma is a crossover, with a crossover temperature T=154(9) MeV (See Chiral and deconfinement aspects of the QCD transition). The equation of state is shown in the Figure 6 below and compared with a similar calculation carried out by the Wuppertal-Budapest collaboration, who used a different discretized lattice action (Stout).