Nucleon axial-vector form factor from lattice QCD, MINERvA hydrogen data, and future neutrino experiments
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
We present some lattice QCD results on first ( χ 1 i ) and second ( χ 2 i ) cumulants of and correlations ( χ 11 i j ) among net baryon-number ( B ), strangeness ( S ) and electric charge ( Q ) along the pseudocritical line [ T p c ( μ B ) ] in the temperature ( T )–baryon chemical potential ( μ B ) phase diagram of ( 2 + 1 )-flavor QCD. We point out that violations of sum rules among second order cumulants, which hold in the isospin symmetric limit of vanishing electric charge chemical potential, are small along the T p c ( μ B ) for the entire range of μ B covered in the RHIC beam energy scan. For the strangeness neutral matter produced in heavy-ion collisions this leads to a close relation between χ 11 B S and χ 11 Q S . We compare lattice QCD results for χ 11 B S / χ 2 S along the T p c ( μ B ) line with preliminary experimental measurements of χ 11 B S / χ 2 S for collision energies 7.7 GeV ≤ s N N ≤ 62.4 GeV . While we find good agreements for s N N ≥ 39 GeV , differences are sizeable at smaller values of s N N . Moreover, we compare lattice QCD results for the ratio of the strangeness ( μ S ) to baryon ( μ B ) chemical potentials, which define a strangeness neutral system with fixed electric charge to baryon number density, with experimental results obtained by the STAR collaboration for μ S / μ B using strange baryon yields on the freeze-out line. Finally, we determine the baryon chemical potential at the freeze-out ( μ B f ) by comparing χ 1 B / χ 2 B along the T p c ( μ B ) with the experimentally measured net-proton cumulants χ 1 p / χ 2 p . We find that { μ B f , T p c ( μ B f ) } are consistent with the freeze-out parameters of the statistical-model fits to experimentally measured hadron yields for s N N ≥ 11.5 GeV . Published by the American Physical Society 2024
We present a first machine learning-based decoding of the gluonic structure of the proton from lattice QCD using a variational autoencoder inverse mapper (VAIM). Harnessing the power of generative AI, we predict the parton distribution function (PDF) of the gluon given information on the reduced pseudo-Ioffe-time distributions (RpITDs) as calculated from an ensemble with lattice spacing a ≈ 0.09 fm and a pion mass of M π ≈ 310 MeV. The resulting gluon PDF is consistent with phenomenological global fits within uncertainties, particularly in the intermediate-to-high-x region where lattice data are most constraining. A subsequent correlation analysis confirms that the VAIM learns a meaningful latent representation, highlighting the potential of generative AI to bridge lattice QCD and phenomenological extractions within a unified analysis framework.
Recent lattice QCD calculations show strong indications that the crossover of QCD at zero baryon chemical potential (µB) is a remnant of the second order chiral phase transition. The non-universal parameters needed to map temperature T and µ B to the universal properties of the second order chiral phase transition were determined by lattice QCD calculations. Motivated by these advances, first, we discuss the analytic structure of the partition function - the so-called Yang-Lee edge singularity - in the QCD crossover regime, solely based on universal properties. Then, utilizing the lattice calculated non-universal parameters, we map this singularity to the real T and complex µ B plane, in order to find the radius of convergence for a Taylor series expansion of QCD partition function around µ B = 0 in the QCD crossover regime. Our most important findings are: (i) An universality-based estimate of the radius of convergence around µ B = 0; (ii) Universality and lattice QCD based constraints on the location of the QCD critical point in the T - µ B plane.
An accurate calculation of the leading-order hadronic vacuum polarization (LOHVP) contribution to the anomalous magnetic moment of the muon ( a μ ) is key to determining whether a discrepancy, suggesting new physics, exists between the Standard Model and experimental results. This calculation can be expressed as an integral over Euclidean time of a current-current correlator G ( t ) , where G ( t ) can be calculated using lattice QCD or, with dispersion relations, from experimental data for e + e − → hadrons . The BMW/DMZ collaboration recently presented a hybrid approach in which G ( t ) is calculated using lattice QCD for most of the contributing t range, but using experimental data for the largest t (lowest energy) region. Here we study the advantages of varying the position t = t 1 separating lattice QCD from data-driven contributions. The total LOHVP contribution should be independent of t 1 , providing both a test of the experimental input and the robustness of the hybrid approach. We use this criterion and a correlated fit to show that Fermilab/HPQCD/MILC lattice QCD results from 2019 strongly favor the CMD-3 cross section data for e + e − → π + π − over a combination of earlier experimental results for this channel. Further, the resulting total LOHVP contribution obtained is consistent with the result obtained by BMW/DMZ, and supports the scenario in which there is no significant discrepancy between the experimental value for a μ and that expected in the Standard Model. We then discuss how improved lattice results in this hybrid approach could provide a more accurate total LOHVP across a wider range of t 1 values with an uncertainty that is smaller than that from either lattice QCD or data-driven approaches on their own. Published by the American Physical Society 2025
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.
This report was prepared on behalf of three Energy Frontier Topical Groups of the Snowmass 2021 Community Planning Exercise. It summarizes the status and implications of studies of strong interactions in high-energy experiments and QCD theory. We emphasize the rich landscape and broad impact of these studies in the decade ahead. Hadronic interactions play a central role in the high-luminosity Large Hadron Collider (LHC) physics program, and strong synergies exist between the (HL-)LHC and planned or proposed experiments at the U.S. Electron-Ion Collider, CERN forward physics experiments, high-intensity facilities, and future TeV-range lepton and hadron colliders. Prospects for precision determinations of the strong coupling and a variety of nonperturbative distribution and fragmentation functions are examined. We also review the potential of envisioned tests of new dynamical regimes of QCD in high-energy and high-density scattering processes with nucleon, ion, and photon initial states. The important role of the high-energy heavy-ion program in studies of nuclear structure and the nuclear medium, and its connections with QCD involving nucleons are summarized. We address ongoing and future theoretical advancements in multi-loop QCD computations, lattice QCD, jet substructure, and event generators. Cross-cutting connections between experimental measurements, theoretical predictions, large-scale data analysis, and high-performance computing are emphasized.
We present a lattice QCD calculation of the nucleon electric polarizabilities at the physical pion mass. Our findings reveal the substantial contributions of the N π states to these polarizabilities. Without considering these contributions, the lattice results fall significantly below the experimental values, consistent with previous lattice studies. This observation has motivated us to compute both the parity-negative N π scattering up to a nucleon momentum of ∼ 0.5 GeV in the center-of-mass frame and corresponding N γ * → N π matrix elements using lattice QCD. Our results confirm that incorporating dynamic N π contributions is crucial for a reliable determination of the polarizabilities from lattice QCD. This methodology lays the groundwork for future lattice QCD investigations into various other polarizabilities. Published by the American Physical Society 2024
Lattice QCD simulations with chirally symmetric quarks have recently established approximate SU(2) CS and SU(2N F ) symmetries of the quantum effective action in a temperature range above the chiral crossover T ch , in which color-electric interactions between quarks dominate the dynamics. We show that such an intermediate temperature range between the chirally broken and plasma regimes is fully consistent with published screening mass spectra, which demonstrate the breakdown of thermal perturbation theory at the crossover between the partonic and the chiral spin symmetric regime at T s ~(2–3)T ch . From the known behavior of screening masses with baryon chemical potential, we deduce qualitatively how this chiral spin symmetric band extends into the QCD phase diagram. In the cold and dense region, we propose parity doubled baryons as possible candidates for chiral spin symmetric matter. This represents a special case of quarkyonic matter with confinement and restored chiral symmetry, and can smoothly transform to quark matter at sufficiently high densities. Finally, we discuss the potential of dilepton spectra to identify such matter forms.
Parton distributions encode the momentum-space structure and, in their generalizations, the spatial tomography of quarks and gluons inside hadrons, the building blocks of visible matter. We present a unified neural-network approach that learns these distributions directly from matrix elements calculated via numerical simulations of quantum chromodynamics (QCD) on the lattice by fitting two complementary inputs simultaneously: data matched to physical quantities via known momentum-space and coordinate-space formalisms. Utilizing data from both methods stabilizes the extraction and mitigates biases that can arise when either is used alone. We validate the method on controlled mock data and apply it to lattice-QCD matrix elements to extract parton distribution functions (PDFs). We show benefits of such an approach for determining the physical quantities. We further extend the framework to zero-skewness generalized parton distributions and demonstrate nucleon tomography within the same neural-network parameterization. Our results provide an adaptable and systematically improvable approach for extracting partonic distributions from Euclidean correlators. It can incorporate polarization, additional channels, and future experimental constraints from current and future facilities, such as the Electron-Ion Collider.
Decay constants of pseudoscalar mesons D, D s , η c , and vector mesons D*, D$^{*}_{s}$, J/ψ are determined from the N f = 2 + 1 lattice QCD at a lattice spacing a ~ 0.08 fm. For vector mesons, the decay constants defined by tensor currents are given in the $\overline{MS}$ scheme at 2 GeV. The calculation is performed on domain wall fermion configurations generated by the RBC-UKQCD collaborations and the overlap fermion action is used for the valence quarks. Comparing the current results with our previous results at a coarser lattice spacing a ~ 0.11 fm provides a better understanding of the discretization error. We obtain $f$$^{T}_{D^*_s}$($\overline{MS}$, 2 GeV)/f D$^{*}_{s}$ = 0.909(18) with a better precision than our previous result. Combining our f D$^{*}_{s}$ = 277(11) MeV with the total width of D$^{*}_{s}$ determined in a recent study gives a branching fraction 4.26(52) x 10 –5 for D$^{*}_{s}$ leptonic decay.
The three-particle K-matrix, K df,3 , is a scheme-dependent quantity that parametrizes short-range three-particle interactions in the relativistic-field-theory three-particle finite-volume formalism. In this work, we compute its value for systems of three pions at maximal isospin through next-to-leading order (NLO) in Chiral Perturbation Theory (ChPT). We compare the values to existing lattice QCD results and find that the agreement between lattice QCD data and ChPT in the first two coefficients of the threshold expansion of K df,3 is significantly improved with respect to leading order once NLO effects are incorporated.
The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.
We develop a comprehensive framework for extracting the pole position and properties of the doubly-charmed tetraquark $T$$^{+}_{cc}$(3875) from lattice QCD data using the relativistic three-particle formalism. This approach incorporates the effect of the one-pion exchange diagram in DDπ and DD ∗ scattering, making it applicable at energies coinciding with the left-hand cut in the partial-wave projected DD ∗ amplitude. We present an example application of this framework to existing lattice QCD data at mπ = 280 MeV. We solve the integral equations describing the DDπ reaction, use LSZ reduction to determine the corresponding DD ∗ amplitude, and find the values of the infinite-volume two- and three-body K matrices that lead to agreement with lattice DD ∗ phase shifts within their uncertainties. Using these K matrices in the three-particle quantization condition, we describe the finite- volume DD ∗ spectrum and find good agreement with the lattice QCD energies. Our results suggest that, at this pion mass, the tetraquark appears as a pair of subthreshold complex poles whose precise location strongly depends on the value of the DDπ three-particle K matrix.
We present the current Standard Model (SM) prediction for the muon anomalous magnetic moment, a μ , updating the first White Paper (WP20) [1]. The pure QED and electroweak contributions have been further consolidated, while hadronic contributions continue to be responsible for the bulk of the uncertainty of the SM prediction. Significant progress has been achieved in the hadronic light-by-light scattering contribution using both the data-driven dispersive approach as well as lattice-QCD calculations, leading to a reduction of the uncertainty by almost a factor of two. The most important development since WP20 is the change in the estimate of the leading-order hadronic-vacuum-polarization (LO HVP) contribution. A new measurement of the e + e - → π + π - cross section by CMD-3 has increased the tensions among data-driven dispersive evaluations of the LO HVP contribution to a level that makes it impossible to combine the results in a meaningful way. At the same time, the attainable precision of lattice-QCD calculations has increased substantially and allows for a consolidated lattice-QCD average of the LO HVP contribution with a precision of about 0.9%. Adopting the latter in this update has resulted in a major upward shift of the total SM prediction, which now reads $a^{SM}_{μ}$ = 116 592 033 (62) x $10^{-11}$ (530 ppb). When compared against the current experimental average based on the E821 experiment and runs 1–6 of E989 at Fermilab, one finds $a^{exp}_{μ} -a^{SM}_{μ}= 38 (63)$ x $10^{-11}$, which implies that there is no tension between the SM and experiment at the current level of precision. The final precision of E989 (127 ppb) is the target of future efforts by the Theory Initiative. The resolution of the tensions among data-driven dispersive evaluations of the LO HVP contribution will be a key element in this endeavor.
Calculation of many-body correlation functions is one of the critical kernels utilized in many scientific computing areas, especially in Lattice Quantum Chromodynamics (Lattice QCD). It is formalized as a sum of a large number of contraction terms each of which can be represented by a graph consisting of vertices describing quarks inside a hadron node and edges designating quark propagations at specific time intervals. Due to its computation- and memory-intensive nature, real-world physics systems (e.g., multi-meson or multi-baryon systems) explored by Lattice QCD prefer to leverage multi-GPUs. Different from general graph processing, many-body correlation function calculations show two specific features: a large number of computation-/data-intensive kernels and frequently repeated appearances of original and intermediate data. The former results in expensive memory operations such as tensor movements and evictions. The latter offers data reuse opportunities to mitigate the data-intensive nature of many-body correlation function calculations. However, existing graph-based multi-GPU schedulers cannot capture these data-centric features, thus resulting in a sub-optimal performance for many-body correlation function calculations. To address this issue, this paper presents a multi-GPU scheduling framework, MICCO, to accelerate contractions for correlation functions particularly by taking the data dimension (e.g., data reuse and data eviction) into account. This work first performs a comprehensive study on the interplay of data reuse and load balance, and designs two new concepts: local reuse pattern and reuse bound to study the opportunity of achieving the optimal trade-off between them. Based on this study, MICCO proposes a heuristic scheduling algorithm and a machine-learning-based regression model to generate the optimal setting of reuse bounds. Specifically, MICCO is integrated into a real-world Lattice QCD system, Redstar, for the first time running on multiple GPUs. The evaluation demonstrates MICCO outperforms other state-of-art works, achieving up to 2.25× speedup in synthesized datasets, and 1.49× speedup in real-world correlation functions.
Transverse-momentum-dependent parton distributions (TMDs) can be calculated from first principles by computing a related set of Euclidean lattice observables and connecting them via a factorization formula. This work focuses on the leading-power factorization formula connecting the lattice quasi-TMD and continuum Collins TMD for gluons. We calculate the one-loop gluon matching coefficient, which is known to be independent of spin and exhibits no mixing with quarks. We demonstrate that this coefficient satisfies Casimir scaling with respect to the quark matching coefficient at one-loop order. Our result facilitates reliable lattice QCD calculations of gluon TMDs.
We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.