Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quark-gluon plasma”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics↗
DOE OSTI2026

Practical considerations for measuring global spin density matrix elements of vector mesons in heavy-ion collisions

The STAR Collaboration has reported a significant 𝜙-meson global spin alignment (𝜌 00 ) signal in Au+Au collisions at $\sqrt{s_{NN}}$ ≤ 62 GeV by measuring the polar angle distribution of 𝜙-meson daughters with respect to the orbital angular momentum direction of the collision system. Here, in this paper, a new method is explored for studying vector-meson global spin alignment in heavy-ion collisions by examining the two-dimensional (2D) polar and azimuthal angle distribution. This method allows simultaneous extraction of 𝜌 00 and off-diagonal spin density matrix elements (SDMEs), providing unique access to local quark-antiquark spin correlations and spin hydrodynamics in quark-gluon plasma. The new 2D method also removes potential biases from nonzero off-diagonal SDMEs on 𝜌 00 with the 1D method. A detailed procedure to correct for detector acceptance and resolution effects is also presented and validated by simulation studies.

Wilks, Gavin [University of Illinois, Chicago, IL ↗
DOE OSTI2024

Λ polarization from vortex rings as the medium response for jet thermalization

We performed a systematic study on the formation of vorticity rings as the process for jet thermalization in the medium created in high-energy nuclear collisions. In this work, we expanded our previous analysis to a more realistic framework by considering noncentral events and fluctuations in the initial condition. We simulate the formation and evolution of the flow vortex structure in a relativistic viscous hydrodynamic model and study the sensitivity of the proposed “ring observable” (ℛ$^{𝑡}_{Λ}$) that can be measured experimentally through the polarization of Λ hyperons. We show that this observable is robust with respect to fluctuating initial conditions to capture the jet-induced vortex flow signal and further study its dependence on different model parameters, such as the jet's velocity, position, the fluid's shear viscosity, and the collision centrality. The proposed observable is associated with the formation of vorticity in a quark-gluon plasma, showing that the measurement of particle polarization can be a powerful tool to probe different properties of jet-medium interactions and to understand better the polarization induced by the transverse and longitudinal expansions of the medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Impact of QCD energy evolution on observables in heavy-ion collisions

Here, we study how the inclusion of energy dependence as dictated by quantum chromodynamic (QCD) small-x evolution equations affects key observables in ultra-relativistic heavy-ion collisions. Specifically, we incorporate JIMWLK evolution into the IP-Glasma framework, which serves as the initial condition for a simulation pipeline that includes viscous relativistic hydrodynamics and a hadronic afterburner. This approach enables a consistent modeling of highly energetic nuclei across varying Bjorken-x values, which are relevant for different collision energies and rapidity regions. In comparison to the standard IP-Glasma setup without small-x evolution, we observe pronounced changes in particle multiplicities and spectral distributions, especially in smaller systems and at the highest available energies. We further explore effects on anisotropic flow observables and correlations between mean transverse momentum and elliptic flow. Our findings underscore the critical role of nonlinear QCD evolution in accurately modeling the early stages of heavy-ion collisions, as well as its implications for extracting transport properties of the quark-gluon plasma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

New approaches to Bayesian uncertainty quantification for Nuclear Science (Final Technical Report)

Inverse problems play a central role in experimentation and theory/data comparisons for many areas of modern Nuclear Physics (NP) and High-Energy Physics (HEP). Bayes’s Theorem is a powerful tool for solving Inverse Problems, providing conceptually transparent and unbiased constraints on theoretical parameters and their uncertainties (“Bayesian Inference”) and enabling the quantification of agreement or tension between models and data. However, analyses based on Bayesian Inference are often challenging for NP and HEP applications, either because of the large number of parameters in the problem, the high computational cost, or both. We propose a multi-institutional collaboration to develop and deploy novel Bayesian analysis tools that advance the scientific scope of a broad range of current and future NP experiments. This project brings together NP domain scientists working on several high-profile NP projects for which new, high-performance Bayesian Uncertainty Quantification (“Bayesian UQ”) methods are essential to carry out the science, and data scientists who are developing state-of-the-art methods applicable to these problems. The NP projects in this proposal comprise measurements of the mass and fundamental nature of the neutrino; study of the Quark-Gluon Plasma that filled the early universe; and mapping of natural and anthropogenic radiation environments. While these NP projects have very different scientific goals, with datasets and analysis approaches that differ significantly, they share common requirements for improving computationally intensive Bayesian analyses using advanced Machine Learning algorithms and will benefit strongly from a coherent effort to develop general solutions. This proposal brings together these projects and forefront ML-based data science algorithms to develop such general solutions. The methods developed in this project will also be more widely applicable, thereby advancing science in the larger Nuclear Physics portfolio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Probing early parton emissions in heavy ion collisions using the Lund jet plane

In scattering experiments, high-virtuality partons, i.e., quarks and gluons, initiate a series of additional parton emissions to create collimated sprays of particles known as jets. This paper presents a measurement of the Lund jet plane (LJP) of high-energy jets produced in lead-lead (PbPb) collisions and compares the results to data for proton-proton (pp) collisions. The LJP is formed by iteratively declustering the constituents of a jet into consecutive emissions and recording the relative transverse momentum ($k_\mathrm{T}$) and angle of the resulting emission with respect to its emitter. The angular distributions of two different $k_\mathrm{T}$ slices of the LJP are investigated for jets with radius parameter of 0.4 and transverse momentum in the range 200$-$1000 GeV. The PbPb (pp) data were recorded by the CMS experiment in 2018 (2017) and correspond to an integrated luminosity of 1.7 nb$^{-1}$ (301 pb$^{-1}$) at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The measurement was designed to test whether the earliest jet emissions are produced before the formation of the quark-gluon plasma (QGP) in PbPb collisions. Within the experimental uncertainties, no significant difference is observed between the angular distribution of high-$k_\mathrm{T}$ emissions in \pp and PbPb collisions, which is consistent with these emissions occurring early in the jet evolution, before substantial interaction with the QGP.

Belyaev, Andrey [Yerevan Phys. Inst.] (ORCID:00000↗
DOE OSTI2025

Polarization of the 𝜙 meson in the hadronic phase with nucleon scatterings and a viscous hydrodynamic background

We extend our previous work on the spin alignment of the ϕ vector meson in the hadronic phase of the quark gluon plasma, by including effects of nucleon scatterings. The emission rates are calculated in a realistic hydrodynamic background simulated with the code FLUIDuM , for different beam energies. We find that all the effects taken into account cannot explain the out-of-plane spin alignment of the ϕ meson observed experimentally.

Polarization phenomena↗
DOE OSTI2026

Open quantum system approach to inclusive jet production in heavy-ion collisions

We derive a factorization formula for inclusive jet production in heavy-ion collisions using the tools of Effective Field Theory (EFT). We show how physics at widely separated scales in this process can be systematically separated by matching to EFTs at successively lower virtualities. Owing to a strong scale separation, we recover a vacuum-like DGLAP evolution above the jet scale, while the additional low-energy scales induced by the medium effectively probe the internal structure of the jet. As a result, the cross section can be written as a series with an increasing number of subjets characterized by perturbative matching coefficients each of which is convolved with a distinct function. These functions encode broadening, medium-induced radiations as well as quantum interference such as the Landau-Pomeranchuk-Migdal effect and color coherence dynamics to all orders in perturbation theory. As a first application of this EFT framework, we investigate the case of an unresolved jet and show how the cross section can be factorized and fully separate the jet dynamics from the universal physics of the medium. To compare to the existing literature, we explicitly compute the medium jet function at next-to-leading order in the coupling and leading order in medium opacity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Exploiting ν-dependence of projected energy correlators in HICs

We extend the recently derived factorization formula for energy-energy correlators to study the analytic structure of general ν-point projected energy correlators in heavy ion collisions. The ν-point projected energy correlators (or, ν-correlators) are an analytically continued family of the integer N-point projected energy correlators, which probe correlations between N final-state particles. By tracking the largest separation (χ) between the N particles, in vacuum, their structure is closely related to the DGLAP splitting functions and exhibits a classical scaling behavior ∼ 1/χ which is modified by resummation through the anomalous dimensions. We show that, in a thermal medium, the ν-correlators display non-trivial angular scaling already at the leading order in perturbation theory. We find that for non-integer values, particularly ν < 1, medium-induced jet function is enhanced compared to ν > 1. This is particularly manifested in the ratios of ν-correlators with respect to the two-point energy correlator, which encodes intrinsic angular information for ν < 1 when compared to large ν values. Moreover, for small-ν values, the ν-correlators appear to saturate at ν = 0.01 . We further confirm our leading-order numerical computations against simulated events from JEWEL for the parton-level production cross-section. Finally, we qualitatively discuss the effect of BFKL resummation for various values of ν.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Measurement of inclusive $J/ψ$ production in Au+Au collisions at $\sqrt{s_{\textrm{NN}}}$ = 54.4 GeV at STAR

This article presents measurements of inclusive J/ψ production at midrapidity (|y| < 1.0) in Au+Au collisions at $\sqrt{s_{\textrm{NN}}}$ = 54.4 GeV with the STAR detector at the Relativistic Heavy Ion Collider. A suppression of the J/ψ yield, quantified using the nuclear modification factors (R AA , R CP ), is observed with respect to the scaled production p + p in collisions. The dependence of R AA on collision centrality and J/ψ transverse momentum is measured with improved precision compared to previous measurements at 39 and 62.4 GeV, while the centrality dependence of R CP is measured and compared to the same results at 39, 62.4, and 200 GeV. In central collisions, no significant collision energy dependence of R AA is found within uncertainties for collision energies between 17.3 and 200 GeV. Two transport model calculations that include dissociation and regeneration contributions are consistent with the experimental results within uncertainties. Although no significant collision energy dependence of the J/ψ suppression in high energy heavy-ion collisions up to $\sqrt{s_{\textrm{NN}}}$ = 200 GeV is observed within uncertainties, the newly measured results at 54.4 GeV Au+Au collisions provide additional constraints on theoretical calculations of the hot medium evolution and cold nuclear matter effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

System-size dependence of charged-particle suppression in ultrarelativistic nucleus-nucleus collisions

High-energy partons lose energy while propagating through the hot, strongly interacting medium produced in ultrarelativistic nucleus-nucleus collisions, leading to a suppression of particle production at high transverse momentum (p T ). The dependence of this energy loss on the size of the colliding nuclear system has yet to be firmly established experimentally. This Letter presents a systematic study of charged-particle suppression across four different nucleus-nucleus collision systems using nuclear modification factors (R AA ) measured by the CMS Collaboration at the CERN LHC. Previous CMS measurements of R AA in oxygen-oxygen, xenon-xenon, and lead-lead collisions are recast with identical p T intervals and are complemented by the first measurement of the charged-particle R AA in neon-neon collisions at $\sqrt{S_{NN}}$ = 5.36 Te V . The neon-neon data correspond to an integrated luminosity of 0.76 nb −1 . The R AA in all collision systems examined show similar qualitative trends as a function of p T , but have a magnitude which is ordered with the nucleon number A. The R AA feature a downward slope at low p T , a local minimum at around 5–7 Ge V , and an upward slope with increasing p T . The R AA are also compared in terms of A 1/3 , which is proportional to the nuclear radius. Models including only initial-state nuclear effects fail to reproduce the observed trends, whereas energy loss models reproduce the trends in the region p T > 9.6 Ge V .

CMS↗
DOE OSTI2025

Universality of scaled particle spectra in ultrarelativistic heavy-ion collisions

We study the transverse momentum spectra of identified particles in ultrarelativistic collisions of large and small collision systems. In order to isolate information contained in the momentum dependence, we propose to scale the spectra by the total particle number and mean transverse momentum, global quantities which are already well studied. We observe an interesting, nearly universal, centrality-independent shape in the scaled spectra, similar to scalings that have been studied previously. This scaling behavior breaks down at large transverse momentum and for very small systems, such as those produced in p − p collisions. We perform hybrid hydrodynamic simulations and show that, in these simulations, a centrality-independent shape is a consequence of an event-by-event independence. Here, our results motivate further theoretical and experimental investigations of the regime of validity of this scaling phenomenon and their physical interpretation at different collision energies and systems.

Collective flow↗
DOE OSTI2026

Charm quark evolution in the early stages of heavy-ion collisions

Heavy quarks are predominantly generated at the initial stage of relativistic heavy-ion collisions such that heavy flavor observables have the potential to provide information on the pre-equilibrium medium dynamics. In this study, we investigate the sensitivity of D-meson $R_{AA}$ and $v_2$ to early-time charm quark dynamics in Pb+Pb collisions at $\sqrt{s_{NN}}$=5.02 TeV. We employ the IP-Glasma+MUSIC+UrQMD framework to model the evolution of the bulk medium. Charm quarks are generated using PYTHIA with nuclear parton distribution functions and evolved using Langevin dynamics within MARTINI. We observe that even though there is significant momentum broadening in the earliest stage, D-meson $R_{AA}$ and $v_2$ are only weakly sensitive to pre-equilibrium interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Locating the QCD critical point through contours of constant entropy density

We propose a new method to investigate the existence and location of the conjectured high-temperature critical point of strongly interacting matter via contours of constant entropy density. By approximating these lines as a power series in the baryon chemical potential 𝜇 𝐵 , one can extrapolate them from first-principle results at zero net-baryon density, and use them to locate the quantum chromodynamics (QCD) critical point, including the associated first-order and spinodal lines. As a proof of principle, we employ currently available continuum-extrapolated lattice data from the Wuppertal-Budapest collaboration to find a critical point at a temperature and a baryon chemical potential of 𝑇 𝑐 = 114.3 ± 6.9 MeV and 𝜇 𝐵,𝑐 = 602.1 ± 62.1 MeV, respectively, at expansion order 𝒪⁡(𝜇$^{2}_{𝐵}$). We advocate for a more precise determination of the required expansion coefficients via lattice QCD simulations as a means of pinpointing the location of the critical endpoint in the phase diagram of strongly interacting matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Nonlinear causality of Israel-Stewart theory with diffusion

We present the first fully nonlinear causality constraints in D = 3 + 1 dimensions for Israel-Stewart theory in the presence of energy and number diffusion in the Eckart and Landau hydrodynamic frames, respectively. These constraints are algebraic inequalities that make no assumption on the underlying geometry of the spacetime or the equation of state. In order to highlight the distinct physical and structural behavior of the two hydrodynamic frames, we discuss the special ultrarelativistic ideal gas equation of state considered in earlier literature in D = 1 + 1 dimensions, and show that our general D = 3 + 1 constraints reduce to their results upon an appropriate choice of angles. For this equation of state in both D = 1 + 1 and D = 3 + 1 dimensions one can show that: (i) there exists a region allowed by nonlinear causality in which the baryon current transitions into a spacelike vector in the Landau frame, and (ii) an analogous argument shows that the solutions of the Eckart frame equations of motion never violate the dominant energy condition, assuming nonlinear causality holds. Furthermore, we then compare our results with those from linearized Israel-Stewart theory and show that the linear causality bounds fail to capture the new physical constraints on energy and number diffusion that are successfully obtained through our nonlinear causality approach.

Quark-gluon plasma↗
DOE OSTI2026

Observation of Charmonium Sequential Suppression in Heavy-Ion Collisions at the Relativistic Heavy Ion Collider

We report measurements of charmonium sequential suppression in Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ =200 GeV with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The inclusive yield ratio of 𝜓⁡(2⁢S) to J/𝜓 as a function of transverse momentum is reported, along with the centrality dependence of the double ratio, defined as the 𝜓⁡(2⁢S) to J/𝜓 ratio in heavy-ion collisions relative to that in 𝑝 +𝑝 collisions. In the 0–80% centrality class, the double ratio is found to be 0.41±0.10 (stat)±0.03 (syst)±0.02 (ref), lower than unity with a significance of 5.6 standard deviations. This provides experimental evidence that 𝜓⁡(2⁢S) is significantly more suppressed than J/𝜓 in heavy-ion collisions at RHIC. This sequential suppression pattern seems to increase from peripheral to central collisions, but with no significant dependence on the transverse momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Measurement of fifth- and sixth-order fluctuations of (net-)proton number in Au + Au collisions from phase II of the beam energy scan program at RHIC

We report high-statistics measurements of fifth- and sixth-order factorial cumulants and cumulant ratios of (net-)proton multiplicity distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7–27GeV, using data from the STAR experiment collected during the Beam Energy Scan Phase II at RHIC. Protons and antiprotons are identified at midrapidity (|𝑦| < 0.5) with transverse momentum 0.4 < 𝑝 𝑇 < 2.0GeV/𝑐. The proton factorial cumulants 𝜅 4 , 𝜅 5 , and 𝜅 6 increase with order but exhibit no sign alternation within current uncertainties, offering no evidence for a two-component structure in the proton multiplicity distribution, as might be expected near a first-order phase transition. The cumulant ratios 𝐶 5 /𝐶 1 and 𝐶 6 /𝐶 2 fluctuate around zero in collisions at 0–40% centrality. The results are consistent with both the negative predictions from lattice QCD and the positive trends obtained from the ultrarelativistic quantum molecular dynamics (UrQMD) model. At $\sqrt{s_{NN}}$ ≳ 27GeV, the 𝐶 4 /𝐶 2 and 𝐶 5 /𝐶 1 results are compatible with predictions from lattice QCD, functional renormalization group (FRG), and hadron resonance gas (HRG) models, while UrQMD describes the data better at lower energies. Here, these measurements place constraints on baryon number fluctuations and offer valuable insights into the QCD phase structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗