SEARCH · Engineering Papers

Results for “nonperturbative 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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Spatial String Tension and Its Effects on Screening Correlators in a Thermal QCD Plasma

We calculate the spatial Wilson line correlator for 2 +1 flavor QCD using highly improved staggered quark discretization for fermions and in quenched QCD for a wide range of temperatures, from the chiral crossover temperature 𝑇 pc ≃ 156 MeV or the deconfinement temperature ≃300 MeV, respectively, up to 2 GeV. Extracting the spatial string tension for different lattice cutoffs and by performing a continuum extrapolation of this observable, we show that the soft (magnetic) gluons interact nonperturbatively even at temperatures ≳1 GeV. We provide incriminating evidences to demonstrate that dimensionally reduced effective theories can describe these soft quark and gluon quasi-particles for both quenched and 2 +1 flavor QCD, at temperatures 𝑇 ≳ 5⁢𝑇 pc . We also show for the first time the imprints of the nonperturbative pseudopotential in the properties of mesonic screening masses for temperatures ranging from 0.8 to 164 GeV in the quark-gluon plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Tomography of the rho meson in the QCD instanton vacuum: Transverse momentum dependent parton distribution functions

We analyze the rho meson unpolarized and polarized transverse momentum dependent parton distribution functions (TMDPDFs) in the instanton liquid model (ILM). The corresponding TMDs in ILM are approximated by a constituent quark beam function in the leading Fock state multiplied by a rapidity-dependent factor resulting from the staple-shaped Wilson lines, for fixed longitudinal momentum, transverse separation, and rapidity. At the resolution of the ILM, all of the rho meson TMDs are symmetric in parton 𝑥 for fixed transverse momentum, and Gaussian-like in the transverse momentum for fixed parton 𝑥. The latter is a direct consequence of the profiling of the quark zero modes in the ILM. The evolved TMDs at higher rapidity using the Collins-Soper kernel, and higher resolution using the renormalization group, show substantial skewness towards low parton 𝑥.

Gluons↗

DOE OSTI2025

What is the Quark-Gluon Plasma made of?

This article surveys our present understanding of the internal structure of the fully developed quark-gluon plasma at temperatures outside the crossover region. The theoretical part of the review covers perturbative and nonperturbative approaches to quark-gluon plasma structure, in particular, hard-thermal loop effective theory, lattice QCD and the functional renormalization group. The phenomenological part of the review scrutinizes the information that has been derived from bulk observables and hard probes in relativistic heavy ion collisions in terms of how it informs our knowledge about the structure of the quark-gluon plasma. The final section lists possible avenues for future progress.

FOS: Physical sciences↗

DOE OSTI2024

Gauge-invariant renormalization of four-quark operators

We study the renormalization of four-quark operators in one-loop perturbation theory. We employ a coordinate-space gauge-invariant renormalization scheme (GIRS), which can be advantageous compared to other schemes, especially in nonperturbative lattice investigations. From our perturbative calculations, we extract the conversion factors between GIRS and the modified minimal subtraction scheme ( MS ¯ ) at the next-to-leading order. As a by-product, we also obtain the relevant anomalous dimensions in the GIRS scheme. A formidable issue in the study of the four-quark operators is that operators with different Dirac matrices mix among themselves upon renormalization. We focus on both parity-conserving and parity-violating four-quark operators, which change flavor numbers by two units ( Δ F = 2 ). The extraction of the conversion factors entails the calculation of two-point Green’s functions involving products of two four-quark operators, as well as three-point Green’s functions with one four-quark and two bilinear operators. The significance of our results lies in their potential to refine our understanding of QCD phenomena, offering insights into the precision of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and shedding light on the nonperturbative treatment of complex mixing patterns associated with four-quark operators. Published by the American Physical Society 2024

Constantinou, M. (ORCID:0000000269881745)↗

DOE OSTI2026

First constraints on the nonperturbative gluon Collins-Soper kernel

The gluon Collins-Soper kernel, which encodes the rapidity evolution of transverse-momentum-dependent gluon distributions, is constrained for the first time in the nonperturbative regime, for transverse momentum scales $q_{T} \in [ 300\text{ MeV}, 1.3\text{ GeV}]$. The constraints are determined in lattice QCD at a close-to-physical pion mass $M_π= 172(3)\text{ MeV}$, a single lattice spacing $a=0.15\text{ fm}$, and next-to-next-to-leading logarithmic matching in Large-Momentum Effective Theory. These results represent the first step toward a controlled determination of the gluon Collins-Soper kernel in QCD, with eventual phenomenological import and relevance to present and future experiments sensitive to the gluon structure of hadronic matter.

Avkhadiev, Artur [Argonne; MIT, Cambridge, CTP]↗

DOE OSTI2026

Quantum Scaling in Energy Correlators beyond the Confinement Transition

We study the QCD scaling behavior of the small-angle energy-energy correlator (EEC), focusing on the transition between its perturbative preconfinement and nonperturbative postconfinement regimes. Applying the light-ray operator product expansion (OPE), we develop a formalism that describes the scaling of the EEC with the input energy 𝑄 in the transition and the postconfinement region, where the latter quantum scaling is determined by the 𝐽 = 5 DGLAP anomalous dimension. A key result of our Letter is a novel connection between the light-ray OPE and the dihadron fragmentation function (DFF), where we show that the nonperturbative OPE coefficients correspond to moments of the DFF. This finding establishes a new paradigm for studying hadronization. Our theoretical predictions are validated against Monte Carlo simulations for both 𝑒 + ⁢𝑒 − and 𝑝⁢𝑝 collisions, showing excellent agreement. The potential role of the quantum scaling in the precision determination of 𝛼 𝑠 is also discussed.

Chang, Cyuan-Han [University of Chicago, IL (Unite↗

DOE OSTI2025

Understanding baryon stopping at the BNL Relativistic Heavy Ion Collider top energies

The nucleon exhibits a rich internal structure governed by quantum chromodynamics (QCD), where its electric charge arises from valence quarks, while its spin and mass emerge from complex interactions among valence quarks, sea (anti)quarks, and gluons. At the advent of QCD, an alternative hypothesis emerged suggesting, at high energies, the transport of a nucleon's baryon number could be traced by a nonperturbative configuration of gluon fields connecting its three valence quarks, forming a 𝑌-shaped topology known as the gluon junction. Recent measurements by the STAR experiment are compatible with this scenario. In light of these measurements, this study aims to explore the mechanisms of baryon transport in high-energy nuclear collisions using the pythia-8 framework, which incorporates a state-of-the-art hadronization model with advanced color flow (CF) and color reconnection (CR) mechanisms that mimic signatures of a baryon junction. Within this model setup, we investigate (i) the rapidity slope of the net-baryon distributions in photon-included processes (𝛾 + 𝑝) and (ii) baryon over charge transport in the isobaric (Ru + Ru and Zr + Zr) collisions. Our study highlights the importance of the CF and CR mechanisms in pythia-8, which play a crucial role in baryon transport. The results show that the CF and CR schemes significantly affect the isobaric baryon-to-charge ratio, leading to different predictions for baryon stopping and underscoring the need to account for CF and CR effects in comparisons with experimental measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Toward a microscopic picture of hadronization and multi-parton processes

This project advanced the understanding of how quarks and gluons produced in high-energy collisions transform into the hadrons observed in particle detectors, a fundamental process known as quantum chromodynamics (QCD) hadronization. By combining theoretical calculations, quantum simulation methods, and modern AI techniques, the research developed new tools to study multi-parton dynamics and nonperturbative effects that are essential for interpreting data from current and future nuclear physics experiments. Key outcomes include new theoretical frameworks for jet and hadron measurements, pioneering quantum simulation algorithms for real-time dynamics in field theories, and the development of advanced machine-learning models, such as diffusion models and explainable classifiers, to simulate and analyze collider events. These results are directly relevant to experiments at Jefferson Lab, Brookhaven National Laboratory, and the future Electron-Ion Collider, and they also have a broader impact in areas such as quantum information science and data-driven modeling of complex systems. The project supported the training of graduate students and postdoctoral fellows and contributed to the broader scientific community through publications, workshops, and collaborative activities. Overall, this work provides new insights into the microscopic mechanisms of hadron formation and establishes a foundation for future studies at the intersection of nuclear physics, artificial intelligence, and quantum computing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Baryon-baryon, meson-meson, and meson-baryon interactions in nonrelativistic QCD

Van der Waals potentials describing interactions between color-singlet mesons and/or baryons vanish at leading order in potential nonrelativistic quantum chromodynamics (pNRQCD). This result and constraints from Gauss’s law are used to prove that weakly coupled pNRQCD van der Waals potentials in generic non-Abelian gauge theories with only heavy quarks are too weak to form bound states whose color state is a product of color singlets. Quantum Monte Carlo calculations of four, five, and six quarks with equal masses provide numerical evidence that exotic color configurations are higher energy than products of color-singlet hadrons, suggesting that equal-mass fully heavy tetraquark, pentaquark, and hexaquark bound states do not exist at next-to-leading order in pNRQCD and at all orders in QCD-like theories in which all quark masses are asymptotically large. Mechanisms for generating hadron-hadron bound states are identified, which necessarily involve large quark-mass hierarchies, relativistic effects arising from the presence of sufficiently light quarks, or nonperturbative effects outside the scope of weakly coupled pNRQCD.

Effective field theory↗

DOE OSTI2025

Quantum entanglement correlations in double quark PDFs

Methods from Quantum Information Theory are used to scrutinize quantum correlations encoded in the two-quark density matrix over light-cone momentum fractions x 1 and x 2 . A nonperturbative three quark model light-cone wavefunction predicts significant nonclassical correlations associated with the “entanglement negativity” measure for asymmetric and small quark momentum fractions. We perform one step of QCD scale evolution of the entire density matrix, not just its diagonal (dPDF), by computing collinearly divergent corrections due to the emission of a gluon. Finally, we present first qualitative numerical results for single-step scale evolution of quantum entanglement correlations in double quark PDFs. At a higher Q 2 scale, the nonclassical correlations manifest in the dPDF for nearly symmetric momentum fractions.

Entanglement in field theory↗

DOE OSTI2026

Nuclear cold QCD: Review and future strategy

This review examines data from hadron-nucleus collisions, primarily focusing on hard processes like Drell-Yan, heavy flavor, and quarkonium production. It highlights observed modifications of particle yields as functions of momentum and rapidity, aiming to clarify the underlying QCD effects on cold nuclear matter. It outlines strategies for future experiments, including the Electron-Ion Collider, to distinguish between these effects. Key questions address the universality of suppression mechanisms and the role of nonperturbative physics, providing a road map for upcoming measurements of hadrons on nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Collins-Soper kernel and reduced soft function in lattice QCD

We evaluate the Collins-Soper kernel and the reduced soft function in lattice QCD, incorporating 𝒪⁡(𝛼 𝑠 ) matching corrections. The calculation relies on the evaluation of the quasitransverse momentum–dependent wave function with asymmetric staple-shaped quark bilinear operators and four-point meson form factors. These quantities are computed nonperturbatively using two 𝑁 𝑓 =2 + 1 + 1 twisted-mass fermion ensembles with the same lattice spacing of 𝑎 = 0.093 fm: the first ensemble has a lattice size of 24 3 × 48 and a pion mass of 346 MeV, and the second one has a lattice size of 32 3 × 64 and a pion mass of 261 MeV. The Collins-Soper kernel and the soft function are needed for the determination of the transverse momentum–dependent parton distribution functions.

Lattice QCD↗

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

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 OSTI2025

Bottomonium properties in the quark-gluon plasma from a lattice-QCD informed 𝑇-matrix approach

Recent computations of bottomonium correlation functions with extended sources in lattice-discretized quantum chromodynamics (lQCD) provide new insights into heavy-quark dynamics at distance scales which are of the order of the inverse temperature. We analyze these results employing the thermodynamic 𝑇-matrix framework, in a continued effort to interpret lQCD data for quarkonium correlation functions in a nonperturbative and self-consistently solved quantum-many-body approach to a strongly coupled quark-gluon plasma (QGP). Its key inputs are the in-medium driving kernel (potential) of the scattering equation and an interference function which implements 3-body effects in the quarkonium coupling to the thermal medium. A simultaneous description of lQCD results for the bottomonium correlators with extended operators and the previously analyzed Wilson line correlators only requires minor refinements of the potential but calls for stronger interference effects at larger separation of the bottom quark and antiquark. We then analyze the poles of the self-consistent 𝑇 matrices on the real axis to assess the survival of the various bound states. Here, we estimate the pertinent temperatures where the poles disappear for the various bottomonium states and discuss the relation to the corresponding peaks in the bottomonium spectral functions. We also recalculate the spatial diffusion coefficient of the QGP and find it to be similar to that in our previous study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Lattice calculation of short-range contributions to neutrinoless double-beta decay 𝜋 − → 𝜋 + ⁢𝑒⁢𝑒 at physical pion mass

Neutrinoless double-beta (0⁢𝜈⁢𝛽⁢𝛽) decays provide an excellent probe for determining whether neutrinos are Dirac or Majorana fermions. The short-range matrix elements associated with the 𝜋 − → 𝜋 + ⁢𝑒⁢𝑒 process contribute at leading order in the 0⁢𝜈⁢𝛽⁢𝛽 decay channel 𝑛⁢𝑛 → 𝑝⁢𝑝⁢𝑒⁢𝑒 through pion exchange between nucleons. However, current lattice calculations show notable discrepancies in predicting these short-range contributions. To address this issue, we perform a lattice QCD calculation of the 𝜋 − → 𝜋 + ⁢𝑒⁢𝑒 matrix elements using domain wall fermion ensembles at the physical pion mass generated by the RBC/UKQCD Collaboration. To mitigate contamination from around-the-world effects, we develop a new method to reconstruct and subtract them directly from lattice data. We then perform nonperturbative renormalization in the regularization-independent symmetric momentum-subtraction scheme (RI/SMOM), using the (𝛾 𝜇 , 𝛾 𝜇 ) and ($\not{𝑞}$, $\not{𝑞}$) projectors. Compared with previous studies, this work reduces the uncertainties in the matrix elements and provides an independent cross-check that helps to reconcile the discrepancies among previous lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Kaon gluon parton distribution and momentum fraction from 2+1+1 lattice QCD with high statistics

We present a high-statistics lattice-QCD determination of the kaon gluon parton distribution function and gluon momentum fraction. We use clover valence fermion action to take 1,296,640 kaon-correlator measurements on a highly improved staggered quark ensemble with 𝑎 ≈ 0.12 fm and 310-MeV pion mass generated by the MILC Collaboration. A detailed investigation into the impact of gauge-link smearing on the gluonic matrix elements indicates that five steps of hypercubic smearing offer an effective balance between signal quality and preservation of long-distance physics. We report a nonperturbatively renormalized kaon gluon momentum fraction of Math output error at 𝜇 = 2 GeV in the Math output error scheme. Using reduced pseudo-Ioffe-time-distribution matrix elements and pseudo-parton-distribution-function (PDF) matching, we extract the kaon gluon PDF and compare with the prediction from the Dyson-Schwinger equation and with the pion PDF obtained from the same ensemble.

First-principles calculations↗

DOE OSTI2025

Resummation for lattice QCD calculation of generalized parton distributions at nonzero skewness

Large-momentum effective theory (LaMET) provides an approach to directly calculate the x-dependence of generalized parton distributions (GPDs) on a Euclidean lattice through power expansion and a perturbative matching. When a parton’s momentum becomes soft, the corresponding logarithms in the matching kernel become non-negligible at higher orders of perturbation theory, which requires a resummation. But the resummation for the off-forward matrix elements at nonzero skewness ξ is difficult due to their multi-scale nature. In this work, we demonstrate that these logarithms are important only in the threshold limit, and derive the threshold factorization formula for the quasi-GPDs in LaMET. We then propose an approach to resum all the large logarithms based on the threshold factorization, which is implemented on a GPD model. We demonstrate that the LaMET prediction is reliable for [−1 + x 0 , −ξ − x 0 ] ∪ [−ξ + x 0 , ξ − x 0 ] ∪ [ξ + x 0 , 1 − x 0 ], where x 0 is a cutoff depending on hard parton momenta. Through our numerical tests with the GPD model, we demonstrate that our method is self-consistent and that the inverse matching does not spread the nonperturbative effects or power corrections to the perturbatively calculable regions.

hadronic spectroscopy↗