Search NASA⌕ Search

SEARCH · Search NASA

Results for “lattice quantum chromodynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hadronic Vacuum Polarization for the Muon 𝑔 − 2 from Lattice QCD: Long-Distance and Full Light-Quark Connected Contribution

We present results for the dominant light-quark connected contribution to the long-distance window of the hadronic vacuum polarization (HVP) contribution to the muon 𝑔 − 2 from lattice quantum chromodynamics. Specifically, with a new determination of the lattice scale on MILC’s physical-mass HISQ ensembles, using the Ω− baryon mass, we obtain a result of 𝑎$^{𝑙⁢𝑙,LD}_{𝜇}$⁡(conn) = 400.2⁢(2.3) stat ⁢(3.7) syst ⁢[4.3] total × 10 −10 . Summing this result with our recent determinations of the light-quark connected contributions to the short- and intermediate-distance windows, we obtain a subpercent precision determination of the light-quark-connected contribution to HVP of 𝑎$^{𝑙⁢𝑙}_{𝜇}$⁡(conn) = 655.2⁢(2.3) stat ⁢(3.9) syst ⁢[4.5]total×10 −10 . Finally, as a consistency check, we verify that an independent analysis of the full contribution is in agreement with the sum of individual windows. We discuss our future plans for improvements of our HVP calculations to meet the target precision of the Fermilab 𝑔 − 2 experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Aspects of scaling and scalability for flow-based sampling of lattice QCD

Recent applications of machine-learned normalizing flows to sampling in lattice field theory suggest that such methods may be able to mitigate critical slowing down and topological freezing. However, these demonstrations have been at the scale of toy models, and it remains to be determined whether they can be applied to state-of-the-art lattice quantum chromodynamics calculations. Assessing the viability of sampling algorithms for lattice field theory at scale has traditionally been accomplished using simple cost scaling laws, but as we discuss in this work, their utility is limited for flow-based approaches. We conclude that flow-based approaches to sampling are better thought of as a broad family of algorithms with different scaling properties, and that scalability must be assessed experimentally.

Abbott, Ryan↗
DOE OSTI2026

Comment on “Nonlocal nucleon matrix elements in the rest frame”

In a recent paper, “Nonlocal nucleon matrix elements in the rest frame” [1], it was observed that the next-to-leading order calculations of the renormalization factor can describe, to a few percent accuracy, the logarithm of the lattice quantum chromodynamics (QCD) rest frame matrix elements with separations up to distances of 0.6 fm on multiple lattice spacings. We argue that perturbative QCD breaks down at such a distance scale after resumming the associated large logarithms, while the Ansatz used in the analysis there is not justified in perturbation theory. Besides, we explain the observation in Ref. [1] and demonstrate that the Ansatz fails to describe the data for 𝑧 > 0.3 fm, showing an opposite trend. Finally, although Ref. [1] proposes multiplying the Ansatz by a Gaussian correction model, which is shown to reduce the discrepancy with the data, this does not legitimize the use of perturbative QCD at such distance scales.

Lattice QCD↗
DOE OSTI2025

Polarized and unpolarized gluon PDFs: Generative machine learning applications for lattice QCD matrix elements at short distance and large momentum

Lattice quantum chromodynamics (QCD) calculations share a defining challenge by requiring a small finite range of spatial separation z between quark/gluon bilinears for controllable power corrections in the perturbative QCD factorization, and a large hadron boost p z for a successful determination of collinear parton distribution functions (PDFs). However, these two requirements make the determination of PDFs from lattice data very challenging. We present the application of generative machine learning algorithms to estimate the polarized and unpolarized gluon correlation functions utilizing short-distance data and extending the correlation up to z p z ≲ 14 , surpassing the current capabilities of lattice QCD calculations. We train physics-informed machine learning algorithms to learn from the short-distance correlation at z ≲ 0.36 fm and take the limit, p z → ∞ , thereby minimizing possible contamination from the higher-twist effects for a successful reconstruction of the polarized gluon PDF. We also expose the bias and problems with underestimating uncertainties associated with the use of model-dependent and overly constrained functional forms, such as x α ( 1 − x ) β and its variants to extract PDFs from the lattice data. We propose the use of generative machine learning algorithms to mitigate these issues and present our determination of the polarized and unpolarized gluon PDFs in the nucleon. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Near-threshold states in coupled D D * − D * D * scattering from lattice QCD

The first determination of doubly charmed isospin-0 coupled-channel D D * − D * D * scattering amplitudes from lattice quantum chromodynamics (QCD) is presented. The finite-volume spectrum is computed for three lattice volumes with a light-quark mass corresponding to m π ≈ 391 MeV and is used to extract the scattering amplitudes in J P = 1 + via the Lüscher quantization condition. By analytically continuing the scattering amplitudes to complex energies, a T c c pole corresponding to a virtual bound state is found below D D * threshold. We also find a second pole, T c c ′ , corresponding to a resonance pole below the kinematically closed D * D * channel, to which it has a strong coupling. A nonzero coupling is robustly found between the S -wave D D * and D * D * channels producing a clear cusp in the D D * amplitude at the D * D * threshold energy. This suggests that the experimental T c c ′ should be observable in D D * and D * D * final states at ongoing experiments. Published by the American Physical Society 2025

Whyte, Travis (ORCID:000000027661301X)↗
DOE OSTI2024

U(1) fields from qubits: An approach via D-theory algebra

A new quantum link microstructure was proposed for the lattice quantum chromodynamics (QCD) Hamiltonian, replacing the Wilson gauge links with a bilinear of fermionic qubits, later generalized to D-theory. This formalism provides a general framework for building lattice field theory algorithms for quantum computing. We focus mostly on the simplest case of a quantum rotor for a single compact U(1) field. We also make some progress for non-Abelian setups, making it clear that the ideas developed in the U(1) case extend to other groups. These in turn are building blocks for 1 + 0 -dimensional ( 1 + 0 -D) matrix models, 1 + 1 -D sigma models and non-Abelian gauge theories in 2 + 1 and 3 + 1 dimensions. By introducing multiple flavors for the U(1) field, where the flavor symmetry is gauged, we can efficiently approach the infinite-dimensional Hilbert space of the quantum O(2) rotor with increasing flavors. The emphasis of the method is on preserving the symplectic algebra exchanging fermionic qubits by sigma matrices (or hard bosons) and developing a formal strategy capable of generalization to a SU ( 3 ) field for lattice QCD and other non-Abelian 1 + 1 -D sigma models or 3 + 1 -D gauge theories. For U(1), we discuss briefly the qubit algorithms for the study of the discrete 1 + 1 -D sine-Gordon equation. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Systematic Uncertainties from Gribov Copies in Lattice Calculation of Parton Distributions in the Coulomb Gauge

Recently, a new method has been proposed to compute parton distributions using boosted correlators fixed in the Coulomb gauge (CG) within the framework of large-momentum effective theory. This approach, which does not involve Wilson lines, could greatly improve the efficiency and precision of lattice quantum chromodynamics calculations. However, concerns remain regarding whether systematic uncertainties from Gribov copies, which correspond to ambiguities in lattice gauge-fixing, are adequately controlled. This work assesses the effects of Gribov copies on Coulomb-gauge-fixed quark correlators. We utilize different strategies for Coulomb-gauge fixing, selecting two different groups of Gribov copies based on lattice gauge configurations. We examine the differences in the resulting spatial quark correlators in both vacuum and pion states. Our findings indicate that the statistical errors of the matrix elements from both Gribov copies, regardless of the correlation range, decrease proportionally to the square root of the number of gauge configurations. The difference between the strategies does not show statistical significance compared to the gauge noise, demonstrating that the effect of the Gribov copies can be neglected in practical lattice calculations of quark parton distributions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Probing for the Trace Estimation of a Permuted Matrix Inverse Corresponding to a Lattice Displacement

We report thatpProbing is a general technique that is used to reduce the variance of the Hutchinson stochastic estimator for the trace of the inverse of a large, sparse matrix A. The variance of the estimator is the sum of the squares of the off-diagonal elements of A -1 . Therefore, this technique computes probing vectors that when used in the estimator annihilate the largest off-diagonal elements. For matrices that display decay of the magnitude of |A$^{-1}_{ij}$| with the graph distance between nodes i and j, this is achieved through graph coloring of increasing powers A k . Equivalently, when a matrix stems from a lattice discretization, it is computationally beneficial to find a distance-k coloring of the lattice. Previously, a hierarchical coloring was proposed so that k can be increased at runtime as needed without discarding previous work. 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 quantum chromodynamics (QCD), where we need to construct “disconnected diagrams” to extract flavor-separated generalized parton functions. In lattice QCD, where the matrix has a four-dimensional toroidal lattice structure, these nonlocal 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 does 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.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Domain decomposition efforts in Chroma

Multilevel integration is a technique used in lattice quantum chromodynamics (LQCD) simulations, where the correlation functions are decomposed into factors that depend only on fields localized within specific lattice subdomains. The approach significantly improves the efficiency of the calculations because the factors can be independently integrated. The main challenge in implementing multilevel integration is that the standard formulation of the LQCD action is not local in the presence of fermions. However, past studies have shown that approximations of the quark propagator and the fermionic determinant, which depend on the gauge fields within specific subdomains, can still be effective. We will present the current progress in supporting domain decomposition within the Chroma software, which is necessary for implementing multilevel integration approaches. Most of the efforts have focused on extending the Chroma code base to efficiently manipulate subdomains for lattice fields and operators, as well as performing inversions and eigendecompositions within these domains.

Alcalde, Eloy Romero↗
DOE OSTI2025

Constraints on the finite volume two-nucleon spectrum at 𝑚𝜋 ≈806 MeV

The low-energy, finite-volume spectrum of the two-nucleon system at a quark mass corresponding to a pion mass of 𝑚𝜋≈806 MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Variational study of two-nucleon systems with lattice QCD, Phys. Rev. D 107, 094508 (2023).] are extended by including a complete basis of local hexaquark operators, as well as plane-wave dibaryon operators built from products of both positive- and negative-parity nucleon operators. Results are presented for the isosinglet and isotriplet two-nucleon channels. In both channels, noticeably weaker variational bounds on the lowest few energy eigenvalues are obtained from operator sets which contain only hexaquark operators or operators constructed from the product of two negative-parity nucleons, while other operator sets produce low-energy variational bounds which are consistent within statistical uncertainties. The consequences of these studies for the LQCD understanding of the two-nucleon spectrum are investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Nonlocal nucleon matrix elements in the rest frame

Extracting parton structure from lattice quantum chromodynamics (QCD) calculations requires studying the coordinate scale 𝑧 3 dependence of the matrix elements of bilocal operators. The most significant contribution comes from the 𝑧 3 dependence induced by ultraviolet (UV) renormalization of the Wilson line. We demonstrate that the next-to-leading order perturbative calculations of the renormalization factor can describe, to a few percent accuracy, the logarithm of the lattice QCD rest frame matrix elements with separations up to distances of 0.6 fm on multiple lattice spacings. The residual discrepancies can be modeled by a leading effect from the structure of the nucleon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Lattice QCD calculation of 𝑥-dependent meson distribution amplitudes at physical pion mass with threshold logarithm resummation

We present a lattice quantum chromodynamics (QCD) calculation of the 𝑥-dependent pion and kaon distribution amplitudes (DA) in the framework of large momentum effective theory. This calculation is performed on a fine lattice of 𝑎 = 0.076 fm at physical pion mass, with the pion boosted to 1.8 GeV and kaon boosted to 2.3 GeV. We renormalize the matrix elements in the hybrid scheme and match to Math output error with a subtraction of the leading renormalon in the Wilson-line mass. The perturbative matching is improved by resumming the large logarithms related to the small quark and gluon momenta in the soft-gluon limit. After resummation, we demonstrate that we are able to calculate a range of 𝑥 ∈[𝑥 0 ,1 − 𝑥 0 ] with 𝑥 0 = 0.25 for pion and 𝑥 0 = 0.2 for kaon with theoretical systematic errors under control. The kaon DA is shown to be slighted skewed, and narrower than pion DA. Although the 𝑥-dependence cannot be direct calculated beyond these ranges, we estimate higher moments of the pion and kaon DAs by complementing our calculation with short-distance factorization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Moments of parton distribution functions of the pion from lattice QCD using gradient flow

We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.

lattice QCD↗
DOE OSTI2024

Determination of the Collins-Soper Kernel from Lattice QCD

This Letter presents a determination of the quark Collins-Soper kernel, which relates transverse-momentum-dependent parton distributions (TMDs) at different rapidity scales, using lattice quantum chromodynamics (QCD). This is the first such determination with systematic control of quark mass, operator mixing, and discretization effects. Next-to-next-to-leading logarithmic matching is used to match lattice-calculable distributions to the corresponding TMDs. The continuum-extrapolated lattice QCD results are consistent with several recent phenomenological parametrizations of the Collins-Soper kernel and are precise enough to disfavor other parametrizations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

𝐹 𝐾 ⁡/𝐹 𝜋 from Möbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of 𝐹 𝐾 ⁡/𝐹 𝜋 using Möbius domain-wall fermions computed on gradient-flowed 𝑁 𝑓 =2 +1 +1 highly improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from 130 ≲ 𝑚 𝜋 ≲ 400 MeV , four lattice spacings of 𝑎 ∼ 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 𝑎 ∼ 0.06 fm ensemble is helpful, but not necessary to achieve a subpercent determination of 𝐹 𝐾 ⁡/𝐹 𝜋 . We also include an estimate of the strong isospin breaking corrections and arrive at a final result of 𝐹 $\hat{K}$ + ⁡ /𝐹 $\hat{𝜋}$ + = 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 |𝑉 𝑢⁢𝑠 |/|𝑉 𝑢⁢𝑑 | = 0.2311⁢(10).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A lattice QCD perspective on weak decays of $b$ and $c$ quarks Snowmass 2022 White Paper

Lattice quantum chromodynamics has proven to be an indispensable method to determine nonperturbative strong contributions to weak decay processes. In this white paper for the Snowmass community planning process we highlight achievements and future avenues of research for lattice calculations of weak $b$ and $c$ quark decays, and point out how these calculations will help to address the anomalies currently in the spotlight of the particle physics community. With future increases in computational resources and algorithmic improvements, percent level (and below) lattice determinations will play a central role in constraining the standard model or identifying new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Binding energy of the 𝑇 𝑏⁢𝑏 tetraquark from lattice QCD with relativistic and nonrelativistic heavy-quark actions

We present a new determination of the $b\bar{b}$𝑢⁢𝑑 (𝐽 𝑃 = 1 + , 𝐼 = 0) tetraquark binding energy using lattice quantum chromodynamics (QCD) with domain-wall light quarks and a nonperturbatively tuned three-parameter anisotropic-clover “relativistic” action for the 𝑏 quarks. We also perform a direct comparison with a reanalysis of data generated in prior work using a lattice-nonrelativistic QCD (NRQCD) action for the 𝑏 quarks and otherwise identical parameters. Using the new data with relativistic 𝑏 quarks from seven different ensembles with multiple lattice spacings and pion masses, we perform combined chiral and continuum extrapolations and obtain (𝑚 𝑇 𝑏⁢𝑏 −𝑚 𝐵 −𝑚 𝐵* ) RHQ =(−76 ±23) MeV. For the NRQCD data from five ensembles, we perform chiral-only extrapolations and obtain (𝑚 𝑇 𝑏⁢𝑏 −𝑚 𝐵 −𝑚 𝐵* ) NRQCD = (−74 ±17 ±10) MeV. The lower magnitude of the results obtained here, compared to the original analysis in [Phys. Rev. D 100, 014503 (2019)], is due to the use of the symmetric parts of the correlation matrices with local four-quark operators only.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Extraction of low-energy constants of single- and double-β decays from lattice QCD: A sensitivity analysis

Lattice quantum chromodynamics (LQCD) has the promise of constraining low-energy constants (LECs) of nuclear effective field theories (EFTs) from first-principles calculations that incorporate the dynamics of quarks and gluons. Given the Euclidean and finite-volume nature of LQCD outputs, complex mappings are developed in recent years to obtain the Minkowski and infinite-volume counterparts of LQCD observables. In particular, as LQCD is moving toward computing a set of important few-nucleon matrix elements at the physical values of the quark masses, it is important to investigate whether the anticipated precision of LQCD spectra and matrix elements will be sufficient to guarantee tighter constraints on the relevant LECs than those already obtained from phenomenology, considering the nontrivial mappings involved. With a focus on the leading-order LECs of the pionless EFT, L 1,A and $g$$^{NN}_{v}$, which parametrize, respectively, the strength of the isovector axial two-body current in a single-β decay (and other related processes such pp fusion), and of the isotensor contact two-body operator in the neutrinoless double-β decay within the light neutrino exchange scenario, the expected uncertainty on future extractions of L 1,A and $g$$^{NN}_{v}$ are examined using synthetic data at the physical values of the quark masses. It is observed that achieving small uncertainties in L 1,A will be challenging, and (sub)percent-level precision in the two-nucleon spectra and matrix elements is essential in reducing the uncertainty on this LEC compared to the existing constraints. On the other hand, the short-distance coupling of the neutrinoless double-β decay, $g$$^{NN}_{v}$, is shown to be less sensitive to uncertainties on both LQCD energies and the matrix element, and can likely be constrained with percent-level precision in the upcoming LQCD calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗