SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Herein we review recent theoretical developments relevant to heavy-ion experiments carried out within the Beam Energy Scan program at the Relativistic Heavy Ion Collider. Our main focus is on the description of the dynamics of systems created in heavy-ion collisions and establishing the necessary connection between the experimental observables and the QCD phase diagram.
In a Quark-Gluon Plasma (QGP), the fundamental building blocks of matter, quarks and gluons, are under extreme conditions of temperature and density. A QGP could exist in the early stages of the Universe, and in various objects and events in the cosmos. The thermodynamic and hydrodynamic properties of the QGP are described by Quantum Chromodynamics (QCD) and can be studied in heavy-ion collisions. Despite being a key thermodynamic parameter, the QGP temperature is still poorly known. Thermal lepton pairs (e + e − and μ + μ − ) are ideal penetrating probes of the true temperature of the emitting source, since their invariant-mass spectra suffer neither from strong final-state interactions nor from blue-shift effects due to rapid expansion. Here we measure the QGP temperature using thermal e+e− production at the Relativistic Heavy Ion Collider (RHIC). The average temperature from the low-mass region (in-medium ρ 0 vector-meson dominant) is (2.01 ± 0.23) × 10 12 K, consistent with the chemical freeze-out temperature from statistical models and the phase transition temperature from Lattice QCD. The average temperature from the intermediate mass region (above the ρ0 mass, QGP dominant) is significantly higher at (3.25 ± 0.60) × 10 12 K. This work provides essential experimental thermodynamic measurements to map out the QCD phase diagram and understand the properties of matter under extreme conditions.
It is demonstrated that the presence of a phase transition in heavy-ion collisions, at beam energies that probe dense quantum chromodynamics (QCD) matter, leads to a significant enhancement of the dilepton yield at low invariant mass and requires the embedding of low beam energies per produced pion due to the extended emission time. In addition, the temperature of low-mass dileptons shows a modest decrease due to the mixed phase. The emission of dileptons in the SIS18–SIS100 beam energy range is studied by augmenting the ultra-relativistic quantum molecular dynamics (UrQMD) transport model with a realistic density-dependent equation of state, as well as two different phase transitions. This is achieved by extending the molecular dynamics interaction part of the UrQMD model to a density-dependent interaction potential with a high-density minimum, leading to a phase transition and metastable coexisting high-density states. Together with a high-precision measurement, these simulations will be able to constrain the existence of a phase transition in QCD up to densities of several times nuclear saturation density.
Relativistic collision of isobaric systems is found to be valuable in differentiating the nucleon distributions for nuclei with the same mass number. In recent contrast experiment of $^{96}_{44}$Ru + $^{96}_{44}$Ru versus $^{96}_{40}$Zr + $^{96}_{40}$Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV, the ratios of multiplicity distribution, elliptic flow, triangular flow, and radial flow are precisely measured and found to be significantly different from unity, indicating the difference in the shapes of the isobar pair. Here, in this work, we investigate the feasibility of nuclear structure reconstruction from heavy-ion collision observables. We perform Bayesian Inference with employing the Monte-Carlo Glauber model as an estimator of the mapping from nuclear structure to the final state observables and to provide the mock data for reconstruction. By varying combination of observables included in the mock data, we find it plausible to infer Woods–Saxon parameters from the observables. We also observe that single-system multiplicity distribution for the isobar system, rather than their ratio, is crucial to simultaneously determine the nuclear structure for the isobar system.
The pseudorapidity distribution of charged hadrons produced in Au + Au collisions at a center-of-mass energy of $\sqrt{s_{NN}}$ = 200 GeV is measured using data collected by the sPHENIX detector. Charged hadron yields are extracted by counting cluster pairs in the inner and outer layers of the Intermediate Silicon Tracker, with corrections applied for detector acceptance, reconstruction efficiency, combinatorial pairs, and contributions from secondary decays. The measured distributions cover |η| < 1.1 across various centralities, and the average pseudorapidity density of charged hadrons at mid-rapidity is compared to predictions from Monte Carlo heavy-ion event generators. This result, featuring full azimuthal coverage at mid-rapidity, is consistent with previous experimental measurements at the Relativistic Heavy Ion Collider, thereby supporting the broader sPHENIX physics program.
Artificial intelligence (AI) generative models, such as generative adversarial networks (GANs), variational autoencoders, and normalizing flows, have been widely used and studied as efficient alternatives for traditional scientific simulations. However, they have several drawbacks, including training instability and inability to cover the entire data distribution, especially for regions where data are rare. This is particularly challenging for whole-event, full-detector simulations in high-energy heavy-ion experiments, such as sPHENIX at the Relativistic Heavy Ion Collider and Large Hadron Collider experiments, where thousands of particles are produced per event and interact with the detector. This work investigates the effectiveness of denoising diffusion probabilistic models (DDPMs) as an AI-based generative surrogate model for the sPHENIX experiment that includes the heavy-ion event generation and response of the entire calorimeter stack. DDPM performance in sPHENIX simulation data is compared with a popular rival, GANs. Results show that both DDPMs and GANs can reproduce the data distribution where the examples are abundant (low-to-medium calorimeter energies). Nonetheless, DDPMs significantly outperform GANs, especially in high-energy regions where data are rare. Additionally, DDPMs exhibit superior stability compared to GANs. The results are consistent between both central and peripheral centrality heavy-ion collision events. Moreover, DDPMs offer a substantial speedup of approximately a factor of 100 compared to the traditional Geant4 simulation method.
High-energy, heavy-ion collisions can create local domains of chirality-imbalanced quarks, reflecting the topological features of quantum chromodynamics. The chiral magnetic effect (CME) predicts an electric charge separation of quarks in such topological domains along the magnetic field ($\vec{B}$) generated by the passing of two high-Z nuclei. Here, we use a correlation observable Δ𝛾$^{112}$ between charged meson pairs to detect the CME-induced charge separation and a novel event shape selection (ESS) method to mitigate the background effects related to elliptic flow (𝑣 2 ). The ESS method classifies events based on the emission pattern of final-state particles and determines Δ𝛾$^{112}_{ESS}$ from the zero-flow limit. We reconstruct the $\vec{B}$ field direction from the spectator nucleons, which minimizes backgrounds unrelated to the collective motion of the system. In this work, we report the measurements of Δ𝛾$^{112}$ and a background indicator Δ𝛾$^{132}$ in Au+Au collisions from the Brookhaven National Laboratory Relativistic Heavy Ion Collider (RHIC) Beam Energy Scan phase II and at the top RHIC energy. After background suppression, Δ𝛾$^{132}_{ESS}$ aligns with zero, and Δ𝛾$^{112}_{ESS}$ is reduced to no more than 20% of Δ𝛾$^{112}$. We observe a finite residual charge separation with 2.5𝜎, 3𝜎, and 3.2𝜎 significance in the 20–50% centrality range of Au + Au collisions at 11.5, 14.6, and 19.6 GeV. The results at 17.3 and 27 GeV also show positive values but with a lower significance of 1.3𝜎 and 1.1𝜎, respectively. The corresponding Δ𝛾$^{112}_{ESS}$ values at 7.7, 9.2, and 200 GeV are consistent with zero within uncertainties.
We combine, for the first time, event-by-event T R ENTo initial conditions with the relativistic viscous hydrodynamic model v-USPhydro and the Monte Carlo event generator JEWEL to make predictions for the nuclear modification factor R AA and jet azimuthal anisotropies v n { 2 } in $\sqrt{s_{NN}}$ = 5.02 TeV PbPb collisions for multiple centralities and values of the jet cone radius R. The R-dependence of R AA and v 2 { 2 } strongly depends on the presence of recoiling scattering centers. We find a small jet v 3 { 2 } in mid-central collisions and consistent results in wide jet p T regions and centralities with ATLAS data.
Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far-reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons is exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges is being developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients are calculated across beam energies. New techniques are being developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project provides essential guidance to the Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.
Exclusive production of dielectron pairs, γγ → e + e – , is studied using $\mathcal{L}_{int}$ = 1.72 nb –1 of data from ultraperipheral collisions of lead nuclei at $\sqrt{s_{NN}}$ = 5.02 TeV recorded by the ATLAS detector at the LHC. The process of interest proceeds via photon–photon interactions in the strong electromagnetic fields of relativistic lead nuclei. Dielectron production is measured in the fiducial region defined by following requirements: electron transverse momentum $p^{e}_{T}$ > 2.5 GeV, absolute electron pseudorapidity |η e | < 2.5, dielectron invariant mass m ee > 5 GeV, and dielectron transverse momentum $p^{ee}_{T}$ < 2 GeV. Differential cross-sections are measured as a function of m ee , average $p^{e}_{T}$, absolute dielectron rapidity |y ee |, and scattering angle in the dielectron rest frame, |cos θ*|, in the inclusive sample, and also with a requirement of no activity in the forward direction. The total integrated fiducial cross-section is measured to be 215 ± 1 (stat) $^{+23}_{-20}$ (syst.) ± 4 (lumi.) μ b. Within experimental uncertainties the measured integrated cross-section is in good agreement with the QED predictions from the Monte Carlo programs STARLIGHT and SUPERCHIC, confirming the broad features of the initial photon fluxes. The differential cross-sections show systematic differences from these predictions which are more pronounced at high |y ee | and |cos θ*| values.
To meet the demands of the Beam Energy Scan phase-II (BES-II) program, the STAR experiment at the Relativistic Heavy Ion Collider (RHIC) developed a dual real-time framework consisting of a High Level Trigger (HLT) and an Express Data Production system (xProduction). The HLT operates online within the Data Acquisition (DAQ) chain on a dedicated multi-core CPU cluster with the option to offload compute-intensive kernels to Xeon Phi coprocessors. It uses parallelized algorithms, such as the Cellular Automaton (CA) Track Finder, to perform rapid tracking, vertexing, and event filtering. This allows it to select events of interest in real time and provide immediate feedback on detector and beam conditions. In contrast, the xProduction workflow runs concurrently and independently of the DAQ loop. It applies near offline-quality calibration and reconstruction within hours of data collection. The xProduction input is the express data stream, whose content can be enriched by HLT trigger/priority selections under DAQ/HLT resource constraints, and it uses the STAR calibration/conditions framework, incorporating online calibration/QA information when available. This enables early preliminary physics analysis, including the reconstruction of rare signals, such as hyperons and hypernuclei. It also provides collaboration-wide access to analysis-ready datasets. Together, the HLT and xProduction systems form a complementary architecture: the HLT performs online event selection while the xProduction chain delivers high-quality results within a short amount of time. This integrated framework has enabled the prompt reconstruction of the $^5_Λ$ He hypernucleus with high statistical significance and the efficient processing of hundreds of millions of heavy-ion collision events. In conclusion, its demonstrated scalability and robustness establish a model for future high-luminosity experiments requiring both online event filtering and rapid access to analysis-quality data.
We report directed flow (v 1 ) of multistrange baryons ($Ξ$ and Ω) and improved v 1 data for K – ,$\overline{p}$, $\overline{Λ}$, and ϕ in Au+Au collisions at $\sqrt{s_{NN}}$ = 27 and 200 GeV from the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). We focus on particles whose constituent quarks are not transported from the incoming nuclei but instead are produced in the collisions. At intermediate impact parameters, we examine quark coalescence behavior for particle combinations with identical quark content, and search for any departure from this behavior (“splitting”) for combinations having non-identical quark content. Under the assumption of quark coalescence for produced quarks, the splitting strength appears to increase with the electric charge difference of the constituent quarks in the combinations, consistent with electromagnetic effect expectations.
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.
Abstract This article summarizes our present knowledge about nuclear matter at the highest energy densities and its formation in relativistic heavy ion collisions. We review what is known about the structure and properties of the quark-gluon plasma and survey the observables that are used to glean information about it from experimental data.
Experimental studies of ultra-relativistic heavy ion collisions at the Large Hadron Collider (LHC) depend crucially on Zero Degree Calorimeters (ZDCs) that measure neutrons produced at near-beam rapidity in nucleus-nucleus collisions. In hadronic nuclear collisions these neutrons are mainly spectator neutrons, those that do not scatter from opposing nucleons during the collision. As a result, the ZDCs provide a vital probe of heavy ion collision geometry. The ZDCs are also essential in the study of ultra-peripheral collisions that are initiated by photons associated with the electric fields of one or both nuclei. Coherent photon emission typically leaves the photon emitter intact, making the observation of no ZDC signal, on one or both sides, a tag of such processes. The ATLAS ZDCs, built prior to Run 1, were substantially upgraded for LHC Run 3. The primary upgrades included replacement of the quartz Cherenkov radiator with H_2-doped fused silica rods; installation of fast air-core signal cables between the ZDC and the ATLAS USA15 cavern; a new LED-based calibration system; and new electronics implemented for readout and fully-digital triggering. The ZDCs were also augmented with new “Reaction Plane Detectors” (RPDs) designed to measure the transverse centroid of multi-neutron showers to allow event-by-event reconstruction of the directed-flow plane in nuclear collisions. The Run 3 ZDC detectors, including the RPDs, are described in detail with emphasis on aspects that are new for Run 3.
We extend a previously constructed T-matrix approach to the quark-gluon plasma (QGP) to include the effects of spin-dependent interactions between partons. Following earlier work within the relativistic quark model, the spin-dependent interactions figure as relativistic corrections to the Cornell potential. When applied to the vacuum spectroscopy of quarkonia, in particular their mass splittings in S- and P-wave states, the issue of the Lorentz structure of the confining potential arises. We confirm that a significant admixture of a vector interaction (to the previously assumed scalar interaction) improves the description of the experimental mass splittings. The temperature corrections to the in-medium potential are constrained by results from thermal lattice-QCD for the equation of state (EoS) and heavy-quark (HQ) free energy in a self-consistent set-up for heavy- and light-parton spectral functions in the QGP. We then deploy the refined in-medium heavy-light T-matrix to compute the charm-quark transport coefficients in the QGP. The vector component of the confining potential, through its relativistic corrections, enhances the friction coefficient for charm quarks in the QGP over previous calculations by tens of percent at low momenta and temperatures, and more at higher momenta. Our results are promising for improving the current phenomenology of open heavy-flavor observables at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC).
We describe a numerical method for simulating stochastic fluid dynamics near a critical point in the Ising universality class. This theory is known as model H, and is expected to govern the nonequilibrium dynamics of quantum chromodynamics (QCD) near a possible critical endpoint of the phase transition between a hadron liquid and the quark-gluon plasma. The numerical algorithm is based on a Metropolis scheme, and automatically ensures that the distribution function of the hydrodynamic variables in equilibrium is independent of the transport coefficients and only governed by the microscopic free energy. We verify dynamic scaling near the critical point of a two and three-dimensional fluid and extract the associated critical exponent z. Here, we find z≃3 in three dimensions, and z≃2 for a two-dimensional fluid. In a finite system, we observe a crossover between the mean field value z=4 and the true critical exponent z≃3 (z≃2 in d=2). This crossover is governed by the values of the correlation length and the renormalized shear viscosity.
State-of-the-art hydrodynamic simulations of the quark-gluon plasma are unable to reproduce the elliptic flow of particles observed at the BNL Relativistic Heavy Ion Collider (RHIC) in relativistic 238 U + 238 U collisions when they rely on information obtained from low-energy experiments for the implementation of deformation in the colliding 238 U ions. Here, we show that this is due to an inappropriate treatment of well-deformed nuclei in the modeling of the initial conditions of the quark-gluon plasma. Past studies have identified the deformation of the nuclear surface with that of the nuclear volume, though these are different concepts. In particular, a volume quadrupole moment can be generated by both a surface hexadecapole and a surface quadrupole moment. This feature was so far neglected in the modeling of heavy-ion collisions, and is particularly relevant for nuclei like 238 U, which is both quadrupole deformed and hexadecapole deformed. With rigorous input from Skyrme density functional calculations, we show that correcting for such effects in the implementation of nuclear deformations in hydrodynamic simulations restores agreement with BNL RHIC data. This brings consistency to the results of nuclear experiments across energy scales, and demonstrates the impact of the hexadecapole deformation of 238 U on high-energy collisions.