Search NASA⌕ Search

SEARCH · Search NASA

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

Direct nonparametric multimessenger constraints on the equation of state of cold dense nuclear matter

We utilize the now substantial amount of astrophysical observations of neutron stars (NSs), along with perturbative quantum chromodynamics (pQCD) calculations at high density, to directly constrain the NS equation of state (EOS). To this end, we construct nonparametric EOS priors by using Gaussian processes trained on 75 EOSs, which include models with either hadrons, hyperons, or quarks at high densities. We create a prior using the full EOS sample (model agnostic), and one prior for each EOS family to test model discrimination. We introduce a novel inference approach, which allows the simultaneous sampling of intrinsic and extrinsic parameters of binary NS mergers, as well as a nonparametric equation of state. We showcase this method in a Bayesian updating scheme by first performing a complete analysis of the binary NS merger event GW170817 with minimal assumptions, and sequentially adding information from x-ray and radio NS observations, along with pQCD calculations. Besides providing standard constraints, such as the pressure at twice nuclear saturation density 𝑝⁡(2⁢𝜌 sat ) = 4.3$^{+0.6}_{−0.6}$ × 10 34 dyne/cm 2 , at 95% confidence level, for the model agnostic prior, our methodology shows how the choice of EOS families used in conditioning changes the inferred astrophysical properties of the EOS, namely tidal deformability and maximum supported NS mass. We find hyperonic priors predicting higher tidal deformabilities for a 1.4⁢𝑀 ⊙ NS, and hadronic priors being preferred by the considered astrophysical data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Kinematically enhanced interpolating operators for boosted hadrons

We propose to use interpolating operators for lattice quantum chromodynamics calculations of highly boosted pions and nucleons with kinematically enhanced ground-state overlap factors at large momentum. Because this kinematic enhancement applies to the signal but not the variance of the correlation function, these interpolating operators can achieve better signal-to-noise ratios at large momentum. We perform proof-of-principle calculations with boosted pions and nucleons using close-to-physical and larger quark masses to explore the utility of our proposal. Results for effective energies and matrix elements, as well as Lanczos ground-state energy estimators, are consistent with theoretical expectations for signal-to-noise improvement at large momenta.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Axion Mass Prediction from Adaptive Mesh Refinement Cosmological Lattice Simulations

The quantum chromodynamics (QCD) axion arises as the pseudo-Goldstone mode of a spontaneously broken Abelian Peccei-Quinn (PQ) symmetry. If the scale of PQ symmetry breaking occurs below the inflationary reheat temperature and the domain wall number is unity, then there is a unique axion mass that gives the observed dark matter (DM) abundance. Computing this mass has been the subject of intensive numerical simulations for decades since the mass prediction informs laboratory experiments. Axion strings develop below the PQ symmetry-breaking temperature, and as the string network evolves, it emits axions that go on to become the DM. A key ingredient in the axion mass prediction is the spectral index of axion radiation emitted by the axion strings. We compute this index in this Letter using the most precise and accurate large-scale simulations to date of the axion-string network leveraging adaptive mesh refinement to achieve the precision that would, otherwise, require a static lattice with 262,144 3 lattice sites. We find a scale-invariant axion radiation spectrum to within 1% precision and find no evidence that the spectral index of radiation evolves with time. Accounting for axion production from strings prior to the QCD phase transition leads us to predict that the axion mass should be approximately 𝑚 𝑎 ∈ (45, 65) μ ⁢eV. However, we provide preliminary evidence that axions are produced in greater quantities from the string-domain-wall network collapse during the QCD phase transition, potentially increasing the mass prediction to as much as 300 μ ⁢eV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Empirical investigation of nuclear correlation function distributions in lattice QCD

Two-point correlation functions of systems with baryon number 𝐵 ∈ {1,2,3,4} are investigated using lattice quantum chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the 𝑂⁡(𝑁) model. The theoretical 𝑂⁡(𝑁) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behavior of the baryon number 𝐵 QCD correlation function at large temporal separation is well reproduced by the 𝑂⁡(𝑁 ∼ 2/𝐵) model distribution.

Lattice field theory↗
DOE OSTI2025

Shearing approach to gauge-invariant Trotterization

Universal quantum simulations of gauge field theories are exposed to the risk of gauge symmetry violations when it is not known how to compile the desired operations exactly using the available gate set. In this article, we show how time evolution can be compiled in an Abelian gauge theory—if only approximately—without compromising gauge invariance, by graphically motivating a block-diagonalization procedure. When gauge-invariant interactions are associated with a “spatial network” in the space of discrete quantum numbers, it is seen that cyclically shearing the spatial network converts simultaneous updates to many quantum numbers into conditional updates of a single quantum number; ultimately, this eliminates any need to pass through (and acquire overlap onto) unphysical intermediate configurations. Shearing is explicitly applied to gauge-matter and magnetic interactions of lattice quantum electrodynamics. The features that make shearing successful at preserving Abelian gauge symmetry may also be found in non-Abelian theories, bringing one closer to gauge-invariant simulations of quantum chromodynamics.

Gauge theories↗
DOE OSTI2025

Measurements of ϒ states production in 𝑝 + 𝑝 collisions at $\sqrt{s}$ = 500 GeV with STAR: Cross sections, ratios, and multiplicity dependence

We report measurements of ϒ⁡(1⁢𝑆), ϒ⁡(2⁢𝑆) and ϒ⁡(3⁢𝑆) production in 𝑝 + 𝑝 collisions at $\sqrt{s}$ =500 GeV by the STAR experiment in year 2011, corresponding to an integrated luminosity ℒ int = 13 pb −1 . The results provide precise cross sections, transverse momentum (𝑝 T ) and rapidity (𝑦) spectra, as well as cross section ratios for 𝑝 T < 10 GeV/c and |𝑦| < 1. The dependence of the ϒ yield on charged particle multiplicity has also been measured, offering new insights into the mechanisms of quarkonium production. The data are compared to various theoretical models: the color evaporation model (CEM) accurately describes the ϒ⁡(1⁢𝑆) production, while the color glass condensate+nonrelativistic quantum chromodynamics (CGC+NRQCD) model overestimates the data, particularly at low 𝑝 T . Conversely, the color singlet model (CSM) underestimates the rapidity dependence. These discrepancies highlight the need for further development in understanding the production dynamics of heavy quarkonia in high-energy hadronic collisions. The trend in the multiplicity dependence is consistent with CGC/saturation and string percolation models or ϒ production happening in multiple parton interactions modeled by PYTHIA 8.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

String-based model with Hagedorn temperature of 𝑇 𝐻 ∼ 300 MeV describes the spectrum of mesons and glueballs

We consider the thermodynamics of a color-confined phase of quantum chromodynamics (QCD) and pure gauge theory within a string-inspired corresponding to a physical spatial dimension, 𝑑 = 3. We show that the physical mass spectrum of massive mesons—in both the strange and nonstrange sectors separately—is reasonably well described and extended by the exponential mass spectrum of open strings, 𝜌⁡(𝑚), characterized by a unique Hagedorn temperature 𝑇 𝐻 = $\sqrt3⁢𝜎/2⁢𝜋$, expressed by the string tension 𝜎. This 𝑇 𝐻 is the value appropriate for 𝑑 = 3 spatial dimensions and is of the order of 𝑇 𝐻 ∼ 300 MeV for typical values of the string tension. It is much larger than the values of 𝑇 𝐻 which have been phenomenologically extracted so far to describe the meson spectrum. Glueball states in pure gauge theory, modeled by closed strings, exhibit a similarly large Hagedorn temperature, highlighting a universal feature of the exponential spectrum. We further analyze the thermodynamic properties of the equation of state at finite temperature and demonstrate that, in the confined phase, the string models agree with lattice QCD results. This lends further support to the recent interpretation of the QCD phase diagram that incorporates strings as relevant degrees of freedom.

Phase transitions↗
DOE OSTI2025

Observation of 𝑊⁢𝑍⁢𝛾 production and constraints on new physics scenarios in proton-proton collisions at $\sqrt{s}$ =13 TeV

A measurement of the 𝑊⁢𝑍⁢𝛾 triboson production cross section is presented. The analysis is based on a data sample of proton-proton collisions at a center-of-mass energy of $\sqrt{s}$ =13 TeV recorded with the CMS detector at the LHC, corresponding to an integrated luminosity of 138 fb −1 . The analysis focuses on the final state with three charged leptons, ℓ ± ⁢𝜈⁢ℓ + ⁢ℓ − , where ℓ = 𝑒 or 𝜇, accompanied by an additional photon. The observed (expected) significance of the 𝑊⁢𝑍⁢𝛾 signal is 5.4 (3.8) standard deviations. The cross section is measured in a fiducial region, where events with an ℓ originating from a tau lepton decay are excluded, to be 5.48 ± 1.11 fb, which is compatible with the prediction of 3.69 ± 0.24 fb at next-to-leading order in quantum chromodynamics. Exclusion limits are set on anomalous quartic gauge couplings and on the production cross sections of massive axionlike particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2026

Schwinger-Keldysh effective action for hydrodynamics with approximate symmetries

We study hydrodynamic theories with approximate symmetries in the recently developed effective action approach on the Schwinger-Keldysh contour. We employ the method of spurious symmetry transformation for small explicit symmetry-breaking parameters to systematically constrain symmetry-breaking effects in the nonequilibrium effective action for hydrodynamics. We apply our method to the hydrodynamic theory of chiral symmetry in quantum chromodynamics at finite temperature and density and its explicit breaking by quark masses. We show that the spurious symmetry and the Kubo-Martin-Schwinger relation dictate that the Ward-Takahashi identity for the axial symmetry, i.e., the partial conservation of axial vector current (PCAC) relation, contains a relaxational term proportional to the axial chemical potential, whose kinetic coefficient is at least of the second order in the quark mass. In the phase where the chiral symmetry is spontaneously broken, and the pseudo-Nambu-Goldstone pions appear as hydrodynamic variables, this relaxation effect is subleading compared to the conventional pion mass term in the PCAC relation, which is of the first order in the quark mass. On the other hand, in the chiral symmetry restored phase, we show that our relaxation term, which is of the second order in the quark mass, becomes the leading contribution to the axial charge relaxation. Therefore, the leading axial charge relaxation mechanism is parametrically different in the quark mass across a chiral phase transition.

Goldstone bosons↗
DOE OSTI2026

Real-time estimators for scattering observables: A full account of finite-volume errors for quantum simulation

The real-time correlators of quantum field theories can be directly probed through new approaches to simulation, such as quantum computing and tensor networks. This provides a new framework for computing scattering observables in lattice formulations of strongly interacting theories, such as lattice quantum chromodynamics. In this paper, we prove that the proposal of real-time estimators of scattering observables is universally applicable to all scattering observables of gapped quantum field theories. All finite-volume errors are exponentially suppressed, and the rate of this suppression is controlled by the regulator considered, namely, a displacement of the spectrum of the theory into the complex plane. A partial restoration of Lorentz symmetry by averaging over different boosts gives an additional suppression of finite volume errors. Our results also apply to the simulation of wave packet scattering, where a similar averaging is performed to construct the wave packets that regulate the finite volume effects. This result represents a necessary key step toward determining a broad class of scattering observables via quantum computing that are currently inaccessible via classical computing. Such observables are relevant for various applications, including hadron spectroscopy, hadron structure, and precision tests of the Standard Model. We also comment on potential applications of our results to traditional computational schemes.

Burbano, Ivan M. [University of California, Berkel↗
DOE OSTI2026

Nucleon axial-vector form factor and radius from radiatively corrected antineutrino scattering data

The nucleon axial-vector form factor, 𝐺 𝐴 , is critical to determine the electroweak interactions of leptons with nucleons. Important examples of processes influenced by 𝐺 𝐴 are elastic (anti)neutrino-nucleon scattering and muon capture by the proton. Sparse experimental data results in a large uncertainty on the momentum dependence of 𝐺 𝐴 and has motivated the consideration of new experimental probes and first-principles lattice quantum chromodynamics (QCD) evaluations. The comparison of new and precise theoretical predictions for 𝐺 𝐴 with future experimental data necessitates the application of radiative corrections to experimentally observable processes. We apply these corrections in the extraction of 𝐺 𝐴 and the associated axial-vector radius from the recent MINERvA antineutrino-hydrogen data, compare the effects from radiative corrections to other uncertainties in neutrino scattering experiments, and discuss the comparison of lattice QCD evaluations to experimental measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Block Lanczos algorithm for lattice QCD spectroscopy and matrix elements

Recent work introduced a new framework for analyzing correlation functions with improved convergence and signal-to-noise properties, as well as rigorous quantification of excited-state effects, based on the Lanczos algorithm and spurious eigenvalue filtering with the Cullum-Willoughby test. Here, we extend this framework to the analysis of correlation-function matrices built from multiple interpolating operators in lattice quantum chromodynamics (QCD) by constructing an oblique generalization of the block Lanczos algorithm, as well as a new physically motivated reformulation of the Cullum-Willoughby test that generalizes to block Lanczos straightforwardly. The resulting block Lanczos method directly extends generalized eigenvalue problem (GEVP) methods, which can be viewed as applying a single iteration of block Lanczos. Block Lanczos provides qualitative and quantitative advantages over GEVP methods analogous to the benefits of Lanczos over the standard effective mass, including faster convergence to ground- and excited-state energies, explicitly computable two-sided error bounds, straightforward extraction of matrix elements of external currents, and asymptotically constant signal-to-noise. No fits or statistical inference are required. Proof-of-principle calculations are performed for noiseless mock-data examples as well as two-by-two proton correlation-function matrices in lattice QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

QCD Theory Meets Information Theory

We present a novel technique to incorporate precision calculations from quantum chromodynamics into fully differential particle-level Monte Carlo simulations. By minimizing an information-theoretic quantity subject to constraints, our reweighted Monte Carlo incorporates systematic uncertainties absent in individual Monte Carlo predictions, achieving consistency with the theory input in precision and its estimated systematic uncertainties. Our method can be applied to arbitrary observables known from precision calculations, including multiple observables simultaneously. It generates strictly positive weights, thus offering a clear path to statistically powerful and theoretically precise computations for current and future collider experiments. As a proof of concept, we apply our technique to event-shape observables at electron-positron colliders, leveraging existing precision calculations of thrust. Our analysis highlights the importance of logarithmic moments of event shapes, which have not been previously studied in the collider physics literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2025

New state of matter between the hadronic phase and the quark-gluon plasma?

Lattice-quantum chromodynamics (QCD) simulations and theoretical arguments hint at the existence of an intermediate phase of strongly interacting matter between a confined hadron gas and a deconfined quark-gluon plasma (QGP). We qualitatively and semiquantitatively explore and differentiate the phase structures in the temperature window from the QCD pseudocritical temperature 𝑇 c ≃ 160 MeV to the pure gluonic deconfinement temperature 𝑇 d ≃ 285 MeV. We propose a three-regime picture using a hadron resonance gas description augmented with the exponential spectrum of strings, corresponding to highly excited mesons and glueballs, based on the analysis of a large number 𝑁 c of colors. We estimate the entropy density from our model to confirm that the lattice-QCD data are bracketed with three regimes, i.e., a hadron gas, a QGP, and a new phase for 𝑇 c ≲ 𝑇 ≲ 𝑇 d . In this new phase, which we name a spaghetti of quarks with glueballs (SQGBs), thermal degrees of freedom of quarks are liberated, yet gluons remain confined in glueballs. Since the Hagedorn temperature 𝑇 H ∼ 285 MeV is universal in the meson and the glueball sectors, in the infinite-𝑁 c limit, the phase diagram in the plane of the baryon chemical potential and the temperature is reduced to one with the confined and deconfined phases and quarkyonic matter at high density. At large but finite 𝑁 c , an SQGB window may open between these phases. We point out that the SQGB has interesting similarities with quarkyonic matter and that this matter in the large-𝑁 c limit is confined as measured by the interaction between heavy quarks, but behaves in other respects like a quasifree gas of quarks. As a result of the extrapolation to 𝑁 c = 3, we present a revised phase diagram with the SQGB phase bounded by thermal crossovers. Finally, we give a quantitative analysis of chiral-symmetry restoration in the SQGB phase.

Color confinement↗
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 OSTI2025

Transverse-momentum-dependent pion structures from lattice QCD: Collins-Soper kernel, soft factor, TMDWF, and TMDPDF

We present the first lattice quantum chromodynamics (QCD) calculation of the pion valence-quark transverse-momentum-dependent parton distribution function (TMDPDF) within the framework of large-momentum effective theory (LaMET). Using correlators fixed in the Coulomb gauge (CG), we computed the quasi-TMD beam function for a pion with a mass of 300 MeV, a fine lattice spacing of 𝑎 =0.06 fm, and multiple large momenta up to 3 GeV. The intrinsic soft functions in the CG approach are extracted from form factors with large momentum transfer, and as a byproduct, we also obtain the corresponding Collins-Soper (CS) kernel. Our determinations of both the soft function and the CS kernel agree with perturbation theory at small transverse separations (𝑏 ⊥ ) between the quarks. At larger 𝑏 ⊥ , the CS kernel remains consistent with recent results obtained using both CG and gauge-invariant TMD correlators in the literature. By combining next-to-leading logarithmic factorization of the quasi-TMD beam function and the soft function, we obtain an 𝑥-dependent pion valence-quark TMDPDF for transverse separations 𝑏 ⊥ ≳1 fm. Interestingly, we find that the 𝑏 ⊥ dependence of the phenomenological parametrizations of TMDPDF for moderate values of 𝑥 are in reasonable agreement with our QCD determinations. In addition, we present results for the transverse-momentum-dependent wave function for a heavier pion with 670 MeV mass.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗