Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heavy ion collisions”

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 631 records · Page 35

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↗

Microscopic constraints for the equation of state and structure of neutron stars: A Bayesian model mixing framework

Bayesian model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒⁢EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20⁢𝑛 0 . The uncertainties of the 𝜒⁢EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2⁢𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.

Bayesian methods↗

Back-to-back relative-excess observable to identify the chiral magnetic effect

Background: The chiral magnetic effect (CME) is extensively studied in heavy-ion collisions at RHIC and LHC. In the commonly used reaction-plane (RP)-dependent, charge-dependent azimuthal correlator ($Δ_γ$), both the close and back-to-back pairs are included. Many backgrounds contribute to the close pairs (e.g., resonance decays, jet correlations), whereas the back-to-back pairs are relatively free of those backgrounds. Purpose: In order to reduce those backgrounds, we propose a new observable, which only focuses on the back-to-back pairs, namely, the relative back-to-back opposite-sign over same-sign pair excess ($r_{BB}$) as a function of the pair azimuthal orientation with respect to the RP ($φ_{BB}$). Methods: We use analytical calculations and toy model simulations to demonstrate the sensitivity of $r_{BB}(φ_{BB})$ to the CME and its insensitivity to backgrounds. Results: With finite CME, the $φ_{BB}$ distribution of $r_{BB}$ shows a clear characteristic modulation. Its sensitivity to background is significantly reduced compared to the previous $Δ_γ$ observable. The simulation results are consistent with our analytical calculations. Conclusions: In this work, our studies demonstrate that the $r_{BB}(φ_{BB})$ observable is sensitive to the CME signal and rather insensitive to the resonance backgrounds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Elliptic flow coalescence to identify the $f_0$(980) content

Here we use a simple coalescence model to generate $f_0$(980) particles for three configurations: a $s\bar{s}$ meson, a $u\bar{u}s\bar{s}$ tetraquark, and a $K^+K^–$ molecule. The phase-space information of the coalescing constituents is taken from a multiphase transport (AMPT) simulation of heavy-ion collisions. It is shown that the number of constituent quarks scaling of the elliptic flow anisotropy can be used to discern $s\bar{s}$ from $u\bar{u}s\bar{s}$ and $K^+K^–$ configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GW190814: Impact of a 2.6 solar mass neutron star on the nucleonic equations of state

Is the secondary component of GW190814 the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system? This is the central question animating this paper. Covariant density functional theory provides a unique framework to investigate both the properties of finite nuclei and neutron stars, while enforcing causality at all densities. By tuning existing energy density functionals we were able to: (i) account for a 2.6 M ⊙ neutron star, (ii) satisfy the original constraint on the tidal deformability of a 1.4 M ⊙ neutron star, and (iii) reproduce ground-state properties of finite nuclei. Yet, for the class of models explored in this work, we find that the stiffening of the equation of state required to support supermassive neutron stars is inconsistent with either constraints obtained from energetic heavy-ion collisions or from the low deformability of medium-mass stars. Furthermore, we speculate that the maximum neutron star mass can not be significantly higher than the existing observational limit and that the 2.6 M ⊙ compact object is likely to be the lightest black hole ever discovered.

79 ASTRONOMY AND ASTROPHYSICS↗

Angular correlation between the fission fragment intrinsic spins

The generation of fission fragment (FF) spins and of their relative orbital angular momentum has been debated for more than six decades, and no consensus has been achieved so far. The interpretation of recent experimental results of Wilson et al., has been challenged by several recent theoretical studies, which are not in agreement with one another. According to the interpretation by Wilson et al., the FF spins emerge very long after scission occurrs. Randrup and Vogt, while agreeing that the FF spins are uncorrelated, conclude on the basis of a phenomenological model that these spins are uncorrelated already before scission. Bulgac et al., in a fully microscopic study demonstrate that the primordial FF spin final values are defined before the emission of prompt neutrons and statistical gammas and are strongly correlated with a relative angle between spins close to 2π/3, a result in full agreement with the present independent analysis. The prompt neutrons and statistical gammas carry a significant amount of angular momentum according to the study of Stetcu et al. which can lead to a decorrelation of the FF bandhead spins of the yrast lines measured by Wilson et al., and which also provides arguments why these measured spins are so different from the primordial FF spins evaluated by Bulgac et al.. Here, I show that the unexpected character of the angular correlation between the primordial FF intrinsic spins, recently evaluated by Bulgac et al., which favors FF intrinsic spins pointing predominantly in opposite directions, can be understood by using simple general phase space arguments. The observation by Wilson et al., that the FF spins are uncorrelated follows from both the results of the microscopic calculation of Bulgac et al., and the present analysis of the full correlated probability distribution of the FF spins together with the relative orbital angular momentum of the primordial FFs. These arguments may apply also to heavy-ion collisions and, since there is no use of specifics of the particle interactions, the present results might apply to atomic and molecular systems as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Entanglement entropy, single-particle occupation probabilities, and short-range correlations

For quantum many-body systems with short-range correlations (SRCs), the intimate relationship between their magnitude, the behavior of the single-particle occupation probabilities at momenta larger than the Fermi momentum, and the entanglement entropy is a new qualitative aspect not studied and exploited yet. A large body of recent condensed matter studies indicates that the time evolution of the entanglement entropy describes the nonequilibrium dynamics of isolated and strongly interacting many-body systems, in a manner similar to the Boltzmann entropy, which is strictly defined for dilute and weakly interacting many-body systems. Both theoretical and experimental studies in nuclei and cold atomic gases have shown that the fermion momentum distribution has a generic behavior n(k)=C/k 4 at momenta larger than the Fermi momentum, due to the presence of SRCs, with approximately 20% of the particles having momenta larger than the Fermi momentum. Further, the presence of the long momentum tails in the presence of SRCs changes the textbook relation between the single-particle kinetic energy and occupation probabilities, n mf ⁡(k) = 1/{1+ exp ⁡β[ϵ⁡(k)-μ]} for momenta very different form the Fermi momentum, particularly for dynamics processes. SRCs induced high-momentum tails of the single-particle occupation probabilities increase the entanglement entropy of fermionic systems, which in its turn affects the dynamics of many nuclear reactions, such as heavy-ion collisions and nuclear fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Probing initial baryon stopping and equation of state with rapidity-dependent directed flow of identified particles

In this work, using a (3 + 1)-dimensional hybrid framework with parametric initial conditions, we study the rapidity-dependent directed flow $v_1 (y)$ of identified particles, including pions, kaons, protons, and lambdas in heavy-ion collisions. Cases involving Au + Au collisions are considered, performed at $\sqrt {^SNN}$ ranging from 7.7 to 200 GeV. The dynamics in the beam direction is constrained using the measured pseudo-rapidity distribution of charged particles and the net proton rapidity distribution. Within this framework, the directed flow of mesons is driven by the sideward pressure gradient from the tilted source, and that of baryons mainly due to the initial asymmetric baryon distribution with respect to the beam axis driven by the transverse expansion. Our approach successfully reproduces the rapidity- and beam energy-dependence of $v_1$ for both mesons and baryons. We find that the $v_1 (y)$ of baryons has strong constraining power on the initial baryon stopping, and together with that of mesons, the directed flow probes the equation of state of the dense nuclear matter at finite chemical potentials.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coordinate versus momentum cuts and effects of collective flow on critical fluctuations

We analyze particle number fluctuations in the crossover region near the critical endpoint of a first-order phase transition by utilizing molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our previous study [V. A. Kuznietsov , ] by incorporating longitudinal collective flow. The scaled variance of particle number distribution inside different coordinate and momentum space acceptances is computed through ensemble averaging and found to agree with earlier results obtained using time averaging, validating the ergodic hypothesis for fluctuation observables. Presence of a sizable collective flow is found to be essential for observing large fluctuations from the critical point in momentum space acceptances. We discuss our findings in the context of heavy-ion collisions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effectiveness of denoising diffusion probabilistic models for fast and high-fidelity whole-event simulation in high-energy heavy-ion experiments

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.

42 ENGINEERING↗

In-medium changes of nucleon cross sections tested in neutrino-induced reactions

Historically, studied in the context of heavy-ion collisions, the extent to which free nucleon-nucleon cross sections are modified in-medium remains undetermined by these data sets. Therefore, we investigate the impact of N N in-medium modifications on neutrino-nucleus cross section predictions using the GiBUU transport model. We find that including an in-medium lowering of the N N cross section and density dependence on Δ excitation improves agreement with MicroBooNE neutrino-argon scattering data. This is observed for both proton and neutral pion spectra in charged-current muon neutrino and neutral-current single pion production data sets. The impact of collision broadening of the Δ resonance is also investigated. The absence of Δ broadening is slightly favored, but the larger uncertainties on the pion production data prevent definitive conclusions. Published by the American Physical Society 2024

Bogart, B. (ORCID:0000000305588934)↗

Covariant formulation of spinodal decomposition in rapidly expanding quark gluon plasma

Quantum chromodynamics (QCD) is expected to have a first order phase transition between the confined hadron gas and the deconfined quark gluon plasma at high baryon densities. This will result in phase boundary effects in the metastable and unstable regions. It is important to include these effects in phenomenological models of heavy ion collisions to identify experimental signatures of a phase transition. This requires building intuition on phase separation in rapidly expanding fluids. In this work we present the covariant equations of relativistic hydrodynamics with a phase boundary, provide prescriptions to extend the equation of state to metastable and unstable regions, and show the effects of spinodal separation in a Bjorken flow. Published by the American Physical Society 2024

Kapusta, Joseph I. (ORCID:0000000259429835)↗

Tensor supercurrent in QCD

An external Abelian magnetic field excites in the QCD vacuum a tensor super current that represents the tensor polarization of the chiral condensate. This tensor supercurrent can be deduced from the chiral lagrangian in the presence of anoma lies; a similar tensor supercurrent emerges in rotating systems at finite chemical potential. We discuss the microscopic origin of this supercurrent and argue that it screens the instanton-anti-instanton $I\overline I$ molecules in the QCD vacuum, similarly to the vector supercurrent screening Abrikosov vortices in a superconductor. A number of possible experimental manifestations of the tensor supercurrent are discussed: i) spin alignment of axial-vector and vector mesons in heavy ion collisions; ii) tensor charge of the nucleon; iii) transversity of quark distributions in polarized nucleons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Soft pions and transport near the chiral critical point

During the expansion of a heavy ion collision, the system passes close to the $\textit{O}$(4) critical point of QCD, and thus the fluctuations of the order parameter $(σ, \overrightarrow{π}$) are expected to be enhanced. In this work, our goal is to compute how these enhanced fluctuations modify the transport coefficients of QCD near the pseudocritical point. We also make a phenomenological estimate for how chiral fluctuations could effect the momentum spectrum of soft pions. We first formulate the appropriate stochastic hydrodynamic equations close to the $\textit{O}$(4) critical point. Then, working in mean field, we determine the correlation functions of the stress tensor and the currents which result from this stochastic real-time theory, and use these correlation functions to determine the scaling behavior of the transport coefficients. The hydrodynamic theory also describes the propagation of pion waves, fixing the scaling behavior of the dispersion curve of soft pions. We present scaling functions for the shear viscosity and the charge conductivities near the pseudocritical point and estimate the absolute magnitude of the critical fluctuations to these parameters and the bulk viscosity. Using the calculated pion dispersion curve, we estimate the expected critical enhancement of soft pion yields, and this estimate provides a plausible explanation for the excess seen in experiment relative to ordinary hydrodynamic computations. Our results motivate further phenomenological and numerical work on the implications of chiral symmetry on real-time properties of thermal QCD near the pseudocritical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Constraining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges

We simultaneously incorporate two common extensions of the hadron resonance gas model, namely the addition of extra, unconfirmed resonances to the particle list and the excluded volume repulsive interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow to constrain them separately. In particular, ratios of second-order susceptibilities like $\chi_{11}^{BQ}/\chi_2^B$ and $\chi_{11}^{BS}/\chi_2^B$ are sensitive only to the baryon spectrum, while fourth-to-second order ratios like $\chi_4^B/\chi_2^B$, $\chi_{31}^{BS}/\chi_{11}^{BS}$, or $\chi_{31}^{BQ}/\chi_{11}^{BQ}$ are mainly determined by repulsive interactions. Analysis of the available lattice results suggests the presence of both the extra states in the baryon-strangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. Furthermore, the modified hadron resonance gas model presented here significantly improves the description of lattice QCD susceptibilities at chemical freeze-out and can be used for the analysis of event-by-event fluctuations in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Far-from-equilibrium slow modes and momentum anisotropy in an expanding plasma

The momentum distribution of particle production in heavy-ion collisions encodes information about thermalization processes in the early-stage quark-gluon plasma. We use kinetic theory to study the far-from-equilibrium evolution of an expanding plasma with an anisotropic momentum-space distribution. We identify slow and fast degrees of freedom in the far-from-equilibrium plasma from the evolution of moments of this distribution. At late times, the slow modes correspond to hydrodynamic degrees of freedom and are naturally gapped from the fast modes by the inverse of the relaxation time, τ R − 1 . At early times, however, there are an infinite number of slow modes with a gap inversely proportional to time, τ − 1 . From the evolution of the slow modes we generalize the paradigm of the far-from-equilibrium attractor to vector and tensor components of the energy-momentum tensor, and even to higher moments of the distribution function that are not part of the hydrodynamic evolution. We predict that initial-state momentum anisotropy decays slowly in the far-from-equilibrium phase and may persist until the relaxation time. Published by the American Physical Society 2024

Brewer, Jasmine (ORCID:0000000230840663)↗

Lattice QCD estimates of thermal photon production from the QGP

Thermal photons produced in heavy-ion collision experiments are an important observable for understanding quark-gluon plasma (QGP). The thermal photon rate from the QGP at a given temperature can be calculated from the spectral function of the vector current correlator. Extraction of the spectral function from the lattice correlator is known to be an ill-conditioned problem, as there is no unique solution for a spectral function for a given lattice correlator with statistical errors. The vector current correlator, on the other hand, receives a large ultraviolet contribution from the vacuum, which makes the extraction of the thermal photon rate difficult from this channel. We therefore consider the difference between the transverse and longitudinal part of the spectral function, only capturing the thermal contribution to the current correlator, simplifying the reconstruction significantly. The lattice correlator is calculated for light quarks in quenched QCD at T = 470 MeV ( ∼ 1.5 T c ), as well as in 2 + 1 flavor QCD at T = 220 MeV ( ∼ 1.2 T p c ) with m π = 320 MeV . In order to quantify the nonperturbative effects, the lattice correlator is compared with the corresponding NLO + LPM LO estimate of correlator. The reconstruction of the spectral function is performed in several different frameworks, ranging from physics-informed models of the spectral function to more general models in the Backus-Gilbert method and Gaussian process regression. We find that the resulting photon rates agree within errors. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Strangeness-correlations on the pseudocritical line in ( 2 + 1 )-flavor QCD

We present some lattice QCD results on first ( χ 1 i ) and second ( χ 2 i ) cumulants of and correlations ( χ 11 i j ) among net baryon-number ( B ), strangeness ( S ) and electric charge ( Q ) along the pseudocritical line [ T p c ( μ B ) ] in the temperature ( T )–baryon chemical potential ( μ B ) phase diagram of ( 2 + 1 )-flavor QCD. We point out that violations of sum rules among second order cumulants, which hold in the isospin symmetric limit of vanishing electric charge chemical potential, are small along the T p c ( μ B ) for the entire range of μ B covered in the RHIC beam energy scan. For the strangeness neutral matter produced in heavy-ion collisions this leads to a close relation between χ 11 B S and χ 11 Q S . We compare lattice QCD results for χ 11 B S / χ 2 S along the T p c ( μ B ) line with preliminary experimental measurements of χ 11 B S / χ 2 S for collision energies 7.7 GeV ≤ s N N ≤ 62.4 GeV . While we find good agreements for s N N ≥ 39 GeV , differences are sizeable at smaller values of s N N . Moreover, we compare lattice QCD results for the ratio of the strangeness ( μ S ) to baryon ( μ B ) chemical potentials, which define a strangeness neutral system with fixed electric charge to baryon number density, with experimental results obtained by the STAR collaboration for μ S / μ B using strange baryon yields on the freeze-out line. Finally, we determine the baryon chemical potential at the freeze-out ( μ B f ) by comparing χ 1 B / χ 2 B along the T p c ( μ B ) with the experimentally measured net-proton cumulants χ 1 p / χ 2 p . We find that { μ B f , T p c ( μ B f ) } are consistent with the freeze-out parameters of the statistical-model fits to experimentally measured hadron yields for s N N ≥ 11.5 GeV . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗