Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Equilibrium expectations for non-Gaussian fluctuations near a QCD critical point

With the highly anticipated results from the Beam Energy Scan II program at RHIC being recently revealed, an understanding of particle-number fluctuations and their significance as a potential signature of a possible QCD critical point is crucial. Early works that embarked on this endeavor sought to estimate the fluctuations due to the presence of a critical point assuming they stay in equilibrium. From these results came the proposal to focus efforts on higher, non-Gaussian, moments of the event-by-event distributions, in particular of the number of protons. These non-Gaussian moments are especially sensitive to critical fluctuations, as their magnitudes are proportional to high powers of the critical correlation length. As the equation of state provides key input for hydrodynamical simulations of heavy-ion collisions, we estimate equilibrium fluctuations from the BEST equation of state (EoS) that includes critical features from the 3D Ising Model. In particular, the proton factorial cumulants and their dependence on non-universal mapping parameters is investigated within the BEST EoS. Furthermore, the correlation length, as a central quantity for the assessment of fluctuations in the vicinity of a critical point, is also calculated in a consistent manner with the scaling equation of state. An understanding of the equilibrium estimates of proton factorial cumulants will be useful for further comparison to estimates of out-of-equilibrium fluctuations in order to determine the magnitude of the observable fluctuations to be expected in heavyion collision experiments, in which the time spent near a critical point is short.

Karthein, Jamie M. [Massachusetts Institute of Tec↗
DOE OSTI2025

Theoretical tools for neutrino scattering: interplay between lattice QCD, EFTs, nuclear physics, phenomenology, and neutrino event generators

Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neutrino scattering. Higher-energy interactions involve a variety of reaction mechanisms including quasi-elastic scattering, resonance production, and deep inelastic scattering that must all be included to reliably predict cross sections for energies relevant to DUNE and other accelerator neutrino experiments. Refined nuclear interaction models in these energy regimes will also be valuable for other applications, such as measurements of reactor, solar, and atmospheric neutrinos. This manuscript discusses the theoretical status, challenges, required resources, and path forward for achieving precise predictions of neutrino-nucleus scattering and emphasizes the need for a coordinated theoretical effort involved lattice QCD, nuclear effective theories, phenomenological models of the transition region, and event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Variance reduction in lattice QCD observables via normalizing flows

Normalizing flows can be used to construct unbiased, reduced-variance estimators for lattice field theory observables that are defined by a derivative with respect to action parameters. This work implements the approach for observables involving gluonic operator insertions in the SU(3) Yang-Mills theory and two-flavor QCD in four space-time dimensions. Variance reduction by factors of 10–60 is achieved in glueball correlation functions and in gluonic matrix elements related to hadron structure, with demonstrated computational advantages. The observed variance reduction is found to be approximately independent of the lattice volume, so volume transfer can be utilized to minimize training costs.

Abbott, Ryan [Columbia U.; MIT, Cambridge, CTP; IA↗
DOE OSTI2025

String-based axial and helicity-flip GPDs: A comparison to lattice QCD

We construct an analytic, string-based representation of the nucleon’s axial and helicity-flip conformal moments of generalized parton distributions (GPDs) that holds for skewness and for both the quark and gluon channels. The starting point is the Mellin-Barnes resummation of the conformal partial-wave expansion, where the moments are parametrized by open- (Reggeon) and closed-string (Pomeron) trajectories with slopes determined by experimental form factors and meson/glueball spectroscopy. The forward limits are fixed by the empirical unpolarized and polarized parton distributions. Polynomiality, crossing symmetry, and support are satisfied by construction. After next-to-leading order Dokshitzer-Gribov-Lipatov-Altarelli-Parisi/Efremov-Radyushkin-Brodsky-Lepage evolution to μ = 2 GeV our analytic framework (i) reproduces some of the currently available lattice moments of E and H˜ in the nonsinglet sector, (ii) predicts sea-quark and gluon polarized moments that will be testable by forthcoming simulations and experiments at Jefferson Lab and the future Electron-Ion Collider, and (iii) yields axial and helicity-flip GPDs in x-space in reasonable agreement with lattice QCD.

Gauge-gravity dualities↗
DOE OSTI2026

Enhancement of the photon emission rate near the QCD critical point

We compute photon emission rate enhancement near the QCD critical point using an effective theory of dynamic critical phenomena and derive a universal photon spectrum. The emission rate scales similar to conductivity, increasing with the correlation length (𝜉), diverging at the critical point. The spectrum exhibits 𝜔⁢𝑑⁢𝑁 𝛾 /𝑑 3 ⁢𝑘 ∝ 𝜔 −1/2 in the scaling regime, with the transition occurring at a frequency comparable to shear damping rate 𝜔 ∼ 𝛾 𝜂 /𝜉 2 , reflecting the nonequilibrium properties of the near-critical liquid.

dynamic critical phenomena↗
DOE OSTI2025

Collision-Energy Dependence in Heavy-Ion Collisions from Nonlinear QCD Evolution

We explore the effects of including the energy dependence determined from evolution equations within the color glass condensate framework on observables in ultrarelativistic heavy-ion collisions. This amounts to integrating the JIMWLK evolution equations into the impact parameter-dependent glasma model, which is then coupled to viscous relativistic hydrodynamics. This methodology allows for a systematic representation of nuclei at specific Bjorken-x values, which are probed at different center-of-mass energies of the collision and rapidities of final state particles. Comparing the methodology to the conventional impact parameter-dependent glasma model, we find significant effects on multiplicity distributions and particle spectra, especially in smaller collision systems at the highest center-of-mass energies. Our results highlight the importance of incorporating nonlinear QCD evolution in the description of heavy-ion collisions at varying center-of-mass energies, as the precise extraction of transport coefficients will be affected. This Letter establishes a robust framework for understanding the quark gluon plasma and nuclear structure at high energy, integrating small-x physics into the initial conditions of heavy-ion collisions.

Mäntysaari, Heikki↗
DOE OSTI2025

Ω 3⁢𝑐 ⁢𝑁⁢𝑁 and Ω 3⁢𝑐 ⁢Ω 3⁢𝑐 ⁢𝑁 systems with HAL QCD potentials

Here, this study employs the Faddeev formalism in configuration space to investigate the Ω 3⁢𝑐 ⁢𝑁⁢𝑁 cluster containing a triply charmed Omega baryon (Ω 3⁢𝑐 ). Using the recently reported HAL QCD 𝑆-wave Ω 3⁢𝑐 ⁢𝑁 potentials in the 3 𝑆 1 and 5 𝑆 2 channels, together with the MT-I–III nucleon-nucleon potential and neglecting the Coulomb force, we find no bound state for the Ω 3⁢𝑐 ⁢𝑛⁢𝑝 system. We predict near-threshold resonances in the 𝐽 𝜋 =5/2 + (maximal total spin) and 𝐽 𝜋 =1/2 + (minimal total spin) states, with resonance energies of 1.1 MeV below and 0.0 MeV at the three-body breakup threshold, respectively, at Euclidean time 𝑡/𝑎 =16. A similar analysis of the Ω 3⁢𝑐⁢ Ω 3⁢𝑐⁢ 𝑁 system likewise reveals no bound states, though a possible resonance is indicated.

few-body systems↗
DOE OSTI2025

String theory and grand unification suggest a submicroelectronvolt QCD axion

Axions, grand unification, and string theory are each compelling extensions of the Standard Model. We show that combining these frameworks imposes strong constraints on the QCD axion mass. Using perturbative unitarity arguments and explicit string compactifications—such as those from the Kreuzer-Skarke (KS) type IIB ensemble—we find that the axion mass is favored to lie within the range 10 −11 eV ≲ 𝑚 𝑎 ≲ 10 −8 eV. This range is directly relevant for near-future axion dark matter searches, including ABRACADABRA/DMRadio and CASPEr. We argue that grand unification and the absence of proton decay suggest a compactification volume that keeps the string scale above the unification scale ( ∼ 10 16 GeV), which in turn limits how heavy the axion can be. The same requirements limit the KS axiverse to have at most ∼ 47 axions. As an additional application of our methodology, we search for axions in the KS axiverse that could explain the recent Dark Energy Spectroscopic Instrument hints of evolving dark energy but find none with high enough decay constant (𝑓 𝑎 ≳ 2.5 ×10 17 GeV); we comment on why such high decay constants and low axion masses are difficult to obtain in string compactifications more broadly.

Axions↗
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↗
DOE OSTI2025

QCD jets in a hot and dense medium: A study of shower formation time and collision kernels

With the use of , a model which considers evolving QCD jets against a fluid dynamical background, it is shown that the introduction of formation time to the parton shower after the initial hard scattering is essential for a simultaneous description of charged hadron and jet R A A . This inclusion also improves jet shape ratios at small angle and jet fragmentation function ratios of leading charged hadrons. The framework is then aimed at a study of the leading-order, next-to-leading-order, and nonperturbative collision kernels. Sizable differences in the modification of jet substructure observables, i.e., jet shape and fragmentation functions, are observed. Such differences are caused by the difference in the radiation rates of relatively soft gluons that survive in the evolution in medium. Published by the American Physical Society 2025

Modarresi-Yazdi, Rouzbeh (ORCID:0000000290377447)↗
DOE OSTI2024

Applications of flow models to the generation of correlated lattice QCD ensembles

Machine-learned normalizing flows can be used in the context of lattice quantum field theory to generate statistically correlated ensembles of lattice gauge fields at different action parameters. This work demonstrates how these correlations can be exploited for variance reduction in the computation of observables. Three different proof-of-concept applications are demonstrated using a novel residual flow architecture: continuum limits of gauge theories, the mass dependence of QCD observables, and hadronic matrix elements based on the Feynman–Hellmann approach. In all three cases, it is shown that statistical uncertainties are significantly reduced when machine-learned flows are incorporated as compared with the same calculations performed with uncorrelated ensembles or direct reweighting. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Charmonium χ c 0 and χ c 2 resonances in coupled-channel scattering from lattice QCD

In order to explore the spectrum of hidden-charm scalar and tensor resonances, we study meson-meson scattering with J P C = 0 + + , 2 + + in the charmonium energy region using lattice QCD. Employing a light-quark mass corresponding to m π ≈ 391 MeV , we determine coupled-channel scattering amplitudes up to around 4100 MeV considering all kinematically relevant channels consisting of a pair of open-charm mesons or a charmonium meson with a light meson. A single isolated scalar resonance near 4000 MeV is found with large couplings to D D ¯ , D s D ¯ s and the kinematically closed D * D ¯ * channel. A single tensor resonance at a similar mass couples strongly to D D ¯ , D D ¯ * and D * D ¯ * . We compare the extracted resonances to contemporary experimental candidate states, previous lattice results and theoretical modeling. In contrast to several other studies, we do not find any significant feature in the scalar amplitudes between the ground state χ c 0 ( 1 P ) and the resonance found around 4000 MeV. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Unpolarized transverse momentum dependent parton distributions of the nucleon from lattice QCD

We present a first lattice QCD calculation of the unpolarized nucleon’s isovector transverse-momentum-dependent parton distribution functions (TMDPDFs), which are essential to predict observables of multiscale, semi-inclusive processes in the standard model. We use a N f = 2 + 1 + 1 MILC ensemble with valence clover fermions on a highly improved staggered quark (HISQ) sea to compute the quark momentum distributions in a large-momentum nucleon on the lattice. The state-of-the-art techniques in renormalization and extrapolation in the correlation distance on the lattice are adopted. The perturbative kernel up to next-to-next-to-leading order is taken into account, and the dependence on the pion mass and the hadron momentum is explored. Our results are qualitatively comparable with phenomenological TMDPDFs, which provide an opportunity to predict high energy scatterings from first principles. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Glue in hadrons at medium resolution and the QCD instanton vacuum

We discuss a general framework for the evaluation of the gluonic form factors in light hadrons at low momentum transfer, in the QCD instanton vacuum. At medium resolution of the order of the inverse mean instanton size, the glue is mostly localized in single or pairs of pseudoparticles and globally constrained by the fluctuations of their topological charges. These pseudoparticles trap light quarks, giving rise to emerging multiflavor ’t Hooft interactions. We explicitly evaluate the gluonic scalar, pseudoscalar, energy-momentum tensor (EMT), and the leading C -odd and C -even three-gluon hadronic form factors, at next to leading order in the instanton density, including molecular clusters of like and unlike instantons. We use the results for the EMT to address the contribution of the gluons in Ji’s mass and spin sum rules, at low resolution. When evolved, our results for the mass and spin composition of the nucleon are shown to be in good agreement with the recently reported lattice results at higher resolution. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Constraints on the U ( 1 ) B − L model from global QCD analysis

We perform the first global QCD analysis of electron-nucleon deep-inelastic scattering and related high-energy data including the beyond the Standard Model U ( 1 ) B − L gauge boson, Z ′ . Contrary to the dark photon case, we find no improvement in the χ 2 relative to the baseline result. The finding allows us to place exclusion limits on the coupling constant of the Z ′ with mass in the range M Z ′ = 2 to 160 GeV. 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 OSTI2025

QCD factorization with multihadron fragmentation functions

Important aspects of quantum chromodynamics (QCD) factorization theorems are the properties of the objects involved that can be identified as universal. One example is that the definitions of parton densities and fragmentation functions for different types of hadrons differ only in the identity of the nonperturbative states that form the matrix elements, but are otherwise the same. This leads to independence of perturbative calculations on nonperturbative details of external states. It also lends support to interpretations of correlation functions as encapsulations of intrinsic nonperturbative properties. These characteristics have usually been presumed to still hold true in fragmentation functions even when the observed nonperturbative state is a small-mass cluster of n hadrons rather than simply a single isolated hadron. However, the multidifferential aspect of cross sections that rely on these latter types of fragmentation functions complicates the treatment of kinematical approximations in factorization derivations. That has led to recent claims that the operator definitions for fragmentation functions need to be modified from the single hadron case with nonuniversal prefactors. With such concerns as our motivation, we retrace the steps for factorizing the unpolarized semi-inclusive e + e − annihilation cross section and confirm that they do apply without modification to the case of a small-mass multihadron observed in the final state. In particular, we verify that the standard operator definition from single hadron fragmentation, with its usual prefactor, remains equally valid for the small-mass n -hadron case with the same hard parts and evolution kernels, whereas the more recently proposed definitions with nonuniversal prefactors do not. Our results reaffirm the reliability of most past phenomenological applications of dihadron fragmentation functions. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗