Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “lattice QCD”

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

Moments of axial-vector GPD from lattice QCD: quark helicity, orbital angular momentum, and spin-orbit correlation

In this work, we present a lattice QCD calculation of the Mellin moments of the twist-2 axial-vector generalized parton distribution (GPD), $\overset{\sim }{H}\left(x,\xi, t\right)$ , at zero skewness, ξ, with multiple values of the momentum transfer, t. Our analysis employs the short-distance factorization framework on ratio-scheme renormalized quasi-GPD matrix elements. The calculations are based on an N f = 2 + 1 + 1 twisted mass fermions ensemble with clover improvement, a lattice spacing of a = 0.093 fm, and a pion mass of m π = 260 MeV. We consider both the iso-vector and iso-scalar cases, utilizing next-to-leading-order perturbative matching while omitting the disconnected contributions and gluon mixing in the iso-scalar case. For the first time, we determine the Mellin moments of $\overset{\sim }{H}$ up to the fifth order. From these moments, we discuss the quark helicity and orbital angular momentum contributions to the nucleon spin, as well as the spin-orbit correlations of the quarks. Additionally, we perform a Fourier transform over the momentum transfer, which allows us to explore the spin structure in the impact-parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Hadronic vacuum polarization for the muon g - 2 from lattice QCD: Complete short and intermediate windows

We present complete results for the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment a μ in the short- and intermediate-distance window regions, which account for roughly 10% and 35% of the total HVP contribution to a μ , respectively. In particular, we perform lattice-QCD calculations for the isospin-symmetric connected and disconnected contributions, as well as corrections due to strong-isospin breaking. For the short-distance window observables, we investigate the so-called log-enhancement effects as well as the significant oscillations associated with staggered quarks in this region. For the dominant, isospin-symmetric light-quark-connected contribution, we obtain a μ l l , SD ( conn ) = 48.139 ( 11 ) stat ( 91 ) syst [ 92 ] total × 10 - 10 and a μ l l , W ( conn ) = 206.90 ( 14 ) stat ( 61 ) syst [ 63 ] total × 10 - 10 . We use Bayesian model averaging to fully estimate the covariance matrix between the individual contributions. Our determinations of the complete window contributions are a μ SD = 69.05 ( 1 ) stat ( 21 ) syst [ 21 ] total × 10 - 10 and a μ W = 236.45 ( 17 ) stat ( 83 ) syst [ 85 ] total × 10 - 10 . This work is part of our ongoing effort to compute all contributions to HVP with an overall uncertainty at the few-permille level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

$Ξ$ 𝑏 → $Ξ$ form factors from lattice QCD and standard-model predictions for $Ξ$ 𝑏 → $Ξ$⁢𝜇 + ⁢𝜇 − and $Ξ$ 𝑏 → $Ξ$⁢𝛾 decays

We present the first lattice QCD determination of the $Ξ$ 𝑏 → $Ξ$ vector, axial-vector, and tensor form factors, which are relevant for the theory of rare decays including $Ξ$ 𝑏 → $Ξ$⁢ℓ + ⁢ℓ − and $Ξ$ 𝑏 → $Ξ$⁢𝛾. The calculation is performed with 2+1 flavors of domain-wall fermions at three different lattice spacings and pion masses in the range from approximately 430 to 230 MeV. The bottom quark is implemented using an anisotropic clover action. Three-point functions with a wide range of source-sink separations and model averaging are used to extract the ground-state contributions. We fit the dependence of the form factors on the momentum transfer, the pion mass, and the lattice spacing using modified 𝑧 expansions that account for subthreshold branch cuts, and apply dispersive bounds and asymptotic behavior constraints to achieve controlled uncertainties in the full semileptonic kinematic region. Using our form factor results, we present standard model predictions for the $Ξ$$^{−}_{𝑏}$ → $Ξ$ − ⁢𝛾 and $Ξ$$^{−}_{𝑏}$ → $Ξ$ − ⁢𝜇 + ⁢𝜇 − branching fractions and two angular observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Exploratory calculation of 𝐾 L → 𝜇 + ⁢𝜇 − decay from lattice QCD at physical pion mass

We compute the complex, long-distance two-photon-exchange amplitude which contributes to the rare 𝐾 L → 𝜇 + ⁢𝜇 − decay from lattice QCD. We use a 24 3 × 64 physical-pion-mass gauge field ensemble at an inverse lattice spacing of 1.023 GeV and a QED ∞ -based formalism. Our implementation strategies for all five non-SU(3)-flavor-suppressed diagram topologies are given in detail. We achieve a 25% statistical precision on the dispersive part of this long-distance amplitude. This calculation is carried out with 2+1 quark flavors and therefore requires the addition of counterterms to compensate for the absence of the Glashow-Iliopoulos-Maiani mechanism. These counterterms are not included in the current calculation and will be the subject of a second paper. Although a direct comparison to experiment cannot yet be made because of those omitted counterterms, the present exploratory calculation allows one to identify principal sources of statistical uncertainty in this calculation. The precision of our results is limited by the reconstruction of the physical contribution of the 𝜂 intermediate state, for which various strategies are tested and compared.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Lattice QCD Benchmark of Proton Helicity and Flavor-Dependent Unpolarized Transverse Momentum-Dependent Parton Distribution Functions at Physical Quark Masses

We present the first lattice QCD calculations of the isovector helicity transverse momentum-dependent parton distribution function (TMDPDF) and the flavor-dependent unpolarized TMDPDFs for up and down quarks in the proton. Our computations utilize domain-wall fermion discretization with physical quark masses. Employing Coulomb-gauge-fixed quark correlation functions within the large-momentum effective theory framework, we access nonperturbative transverse quark separations 𝑏 𝑇 up to approximately 1 fm, corresponding to transverse momenta as low as 200 MeV. Based on the quasi-TMD factorization theorem, we construct renormalization-group–invariant ratios that are equal to the corresponding light-cone TMDPDF ratios. At moderate 𝑥, our results reveal that the isovector helicity and unpolarized TMDPDFs exhibit nearly identical transverse structure up to a normalization factor, and the unpolarized distributions display only mild flavor dependence. These findings not only support key trends observed in recent global analyses but also provide robust, nonperturbative constraints that can distinguish between different parameterizations. This Letter establishes a first-principles benchmark for TMDPDFs, offering valuable input for ongoing and future experimental efforts to map the proton’s three-dimensional structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Elastic and resonance structures of the nucleon from the hadronic tensor in lattice QCD: Implications for neutrino-nucleon scattering and hadron physics

We compute the Euclidean hadronic tensor from charge density operators and extract elastic and resonance structures by employing exponential fits to the four-point function correlator, as well as a Bayesian reconstruction inverse algorithm to obtain the corresponding spectral density for qualitative comparison. We present the determination of the nucleon’s Sachs electric form factor using the hadronic tensor formalism and verify that it is consistent with that from the conventional three-point function calculation. Beyond the elastic peak, we observe a structure located approximately 0.5–0.7 GeV above the nucleon mass in the Bayesian reconstruction. This structure is interpreted as a mixture of the Roper resonance [𝑁⁡(1440)], and states with both positive and negative parities in this mass region, as well as multihadron states. Assuming the observed structure is dominated by 𝐽 𝑃 = 1/2 ± states, we extract the transition electric form factor 𝐺$^*_𝐸$⁡(𝑄 2 ) and the corresponding longitudinal helicity amplitude 𝑆 1/2 ⁡(𝑄 2 ), and compare them with those determined from the CLAS experimental data of nucleon-to-Roper transition. Although fitting to the four-point correlation function or using the inverse algorithm does not resolve individual resonances, it nevertheless enables the determination of total inclusive lepton–nucleon scattering cross sections in appropriate energy bins. This lattice QCD calculation presents the first major step toward studying the inclusive 𝑁 → 𝑋 contributions within the hadronic tensor formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Hadronic light-by-light contribution to the muon anomaly from lattice QCD with infinite volume QED at physical pion mass

The hadronic light-by-light scattering contribution to the muon anomalous magnetic moment, (g–2)⁢/2, is computed in the infinite volume QED framework with lattice QCD. We report $a^{HLbL}_μ$ = 12.47⁢(1.15)⁢(0.95) ×10 –10 where the first error is statistical and the second systematic. The result is mainly based on the 2+1 flavor Möbius domain wall fermion ensemble with inverse lattice spacing a –1 = 1.73 GeV, lattice size L = 5.5 fm, and m π = 139 MeV, generated by the RBC-UKQCD collaborations. The leading systematic error of this result comes from the lattice discretization. This result is consistent with previous determinations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Third moments of nucleon unpolarized, polarized, and transversity parton distribution functions from physical-point lattice QCD

Using forward matrix elements of local leading-twist operators, we present a determination of the isovector third Mellin moments ⟨𝑥 2 ⟩ of nucleon unpolarized, polarized, and transversity parton distribution functions. Two lattice QCD ensembles at the physical pion mass are used, which were generated using a tree-level Symanzik-improved gauge action and 2+1 flavor tree-level improved Wilson Clover fermions coupling via 2-level HEX-smearing. Leveraging a wide set of operators, two extraction methods for the matrix elements, and the automatic inclusion of model uncertainties via bootstrapped model averages, we extract values of the third Mellin moments. Furthermore, this is the first direct calculation of these observables performed at the physical pion mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Accelerating multigrid with streaming chiral SVD for Wilson fermions in lattice QCD

A modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD is presented. A larger basis of test vectors than that used in regular multigrid is calculated by the smoother and truncated by singular value decomposition on the chiral components of the test vectors. The truncated basis is used to form the prolongation and restriction matrices of the multigrid hierarchy. This modification of the setup method is demonstrated to increase the convergence of linear solvers on an anisotropic lattice with m π ≈ 239 MeV from the Hadron Spectrum Collaboration and an isotropic lattice with m π ≈ 220 MeV from the MILC Collaboration. The lattice volume dependence of the method is also examined. Increasing the number of test vectors improves speedup up to a point, but storing these vectors becomes impossible in limited memory resources such as GPUs. To address storage cost, we implement a streaming singular value decomposition of the basis of test vectors on the chiral components and demonstrate a decrease in the number of fine level iterations by a factor of 1.7 for m q ≈ m crit

Iterative methods↗
DOE OSTI2025

Generalized parton distributions from lattice QCD with asymmetric momentum transfer: Tensor case

The calculation of generalized parton distributions (GPDs) in lattice QCD was traditionally done by calculating matrix elements in the symmetric frame. Recent advancements have significantly reduced computational costs by calculating these matrix elements in the asymmetric frame, allowing us to choose the momentum transfer to be in either the initial or final states only. The theoretical methodology requires a new parametrization of the matrix element to obtain Lorentz-invariant amplitudes, which are then related to the GPDs. The formulation and implementation of this approachaveh already been established for the unpolarized and helicity GPDs. Building upon this idea, we extend this formulation to the four leading-twist quark transversity GPDs ($𝐻_𝑇$, $𝐸_𝑇$, $\tilde{𝐻}_𝑇$, $\tilde{𝐸}_𝑇$). We also present numerical results for zero skewness using an 𝑁 𝑓 = 2 +1 +1 ensemble of twisted mass fermions with a clover improvement. The light quark masses employed in these calculations correspond to a pion mass of about 260 MeV. Furthermore, we include a comparison between the symmetric and asymmetric frame calculations to demonstrate frame independence of the Lorentz-invariant amplitudes. Analysis of the matrix elements in the asymmetric frame is performed at several values of the momentum transfer squared, −𝑡, ranging from 0.17 to 2.29 GeV 2 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Lattice-QCD Computable Quark Correlation Functions at Three-Loop Order and Extraction of Splitting Functions

We present the first complete next-to-next-to-next-to-leading-order calculation of the matching coefficients that link unpolarized flavor nonsinglet parton distribution functions with lattice QCD computable correlation functions. By using this high-order result, we notice a reduction in theoretical uncertainties compared to relying solely on previously known lower-order matching coefficients. Furthermore, based on this result we have extracted the three-loop unpolarized flavor nonsinglet splitting function, which is in agreement with the state-of-the-art result. Because of the simplicity of our method, it has the potential to advance the calculation of splitting functions to the desired four-loop order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

EIC Physics from Lattice QCD: investigations beyond leading twist

This project proposes theoretical studies of Quantum Chromodynamics (QCD), the theory describing the strong nuclear force among the building blocks (quarks and gluons) of the visible matter. These appear only confined within hadrons, that make up more than 99% of the mass of the matter. Understanding QCD will significantly advance many aspects of science, from the sub-nuclear interactions to astrophysics, and a quantitative theoretical description is imperative. However, this is a challenging task because QCD is a highly nonlinear theory. We propose hadron structure calculations within lattice QCD (LQCD), an ideal ab initio approach based on space-time discretization, which allows the study of the properties of fundamental particles numerically. This is done by defining the continuous equations on a discrete four-dimensional lattice, which results in equations with hundreds of billions of degrees of freedom, and must be simulated in powerful computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Chiral rank-$k$ truncations for the multigrid preconditioner of Wilson fermions in lattice QCD

We present a modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD. A larger number of test vectors than that used in conventional multigrid is generated by the smoother. This set of test vectors is then truncated by a singular value decomposition on the chiral components of the test vectors, which are subsequently used to form the prolongation and restriction matrices of the multigrid hierarchy. This modification is demonstrated to improve the convergence of linear equations on an anisotropic lattice with 𝑚𝜋 ≈ 280 MeV from the Hadron Spectrum Collaboration and an isotropic lattice with 𝑚𝜋 ≈ 220 MeV from the MILC Collaboration. The lattice volume dependence of the method is also examined.

Whyte, Travis [Jülich Supercomputing Center, Jülic↗
DOE OSTI2025

Lattice QCD Calculation of the Subtraction Function in Forward Compton Amplitude

The subtraction function plays a pivotal role in calculations involving the forward Compton amplitude, which is crucial for predicting the Lamb shift in muonic atoms, as well as the proton-neutron mass difference. In this Letter, we present a lattice QCD calculation of the subtraction function using two domain wall fermion gauge ensembles near the physical pion mass. We utilize a recently proposed subtraction point, demonstrating its advantage in mitigating statistical and systematic uncertainties by eliminating the need for ground-state subtraction. Our results reveal significant contributions from N⁢π intermediate states to the subtraction function. Incorporating these contributions, we compute the proton, neutron, and nucleon isovector subtraction functions at photon momentum transfer Q 2 ϵ[0,2] GeV 2 . For the proton subtraction function, we compare our lattice results with chiral perturbation theory prediction at low Q 2 and with the results from the perturbative operator-product expansion at high Q 2 . Finally, using these subtraction functions as input, we determine their contribution to two-photon exchange effects in the Lamb shift and isovector nucleon electromagnetic self-energy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

𝜂 and 𝜂′ meson production in 𝐽/𝜓 radiative decays from lattice QCD

We report on the computation of amplitudes describing the radiative decay processes 𝐽/𝜓 →𝛾⁢𝜂 and 𝐽/𝜓 →𝛾⁢𝜂′ from first principles via lattice quantum chromodynamics (QCD). Using lattices with two heavier than physical degenerate flavors of light quark where 𝑚 𝜋 ∼391 MeV, and a strange quark tuned to approximately the physical value, we compute three-point correlation functions using optimized meson operators with a range of momenta. The use of optimized operators allows access to the 𝜂′, even though it appears as an excited state, lying above the ground-state 𝜂 in the isoscalar pseudoscalar channel. Statistically precise signals are obtained in this disconnected process by averaging over large numbers of kinematically equivalent correlators. We determine transition form factors as a function of photon virtuality across the timelike region, which describe also the Dalitz processes 𝐽/𝜓 →𝑒 + ⁢𝑒 − ⁢𝜂 (′) . Significantly lower magnitudes of transition form factor for both the 𝜂 and 𝜂′ are found than those extracted from experimental data, and possible explanations for this observation are presented.

form factors↗
DOE OSTI2025

Bootstrap-determined p values in lattice QCD

We present a general method to determine the probability that stochastic Monte Carlo data, in particular those generated in a lattice QCD calculation, would have been obtained were that data drawn from the distribution predicted by a given theoretical hypothesis. Such a probability, or p -value, is often used as an important heuristic measure of the validity of that hypothesis. The proposed method offers the benefit that it remains usable in cases where the standard Hotelling T 2 methods based on the conventional χ 2 statistic do not apply, such as for uncorrelated fits. Specifically, we analyze q 2 , defined as the correlated χ 2 statistic obtained using an arbitrary covariance matrix estimator, and show how to use the bootstrap as a data-driven method to determine the expected distribution of q 2 for a given hypothesis with minimal assumptions. This distribution can then be used to determine the p -value for a fit to the data. We also describe a bootstrap approach for quantifying the impact upon this p -value of estimating population parameters from a single ensemble of N samples. The overall method is accurate up to a 1 / N bias which we do not attempt to quantify. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Optimizing stochastic algorithms for hadron correlation function computations in lattice QCD using a localized distillation basis

Distillation is a quark-smearing method for the construction of a broad class of hadron operators useful in lattice QCD computations and defined via a projection operator into a vector space of smooth gauge-covariant fields. A new orthonormal basis for this space is constructed which builds in locality. This basis is useful for the construction of stochastic methods to estimate the correlation functions computed in Monte Carlo calculations relevant for hadronic physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗