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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Here, we report a systematic measurement of cumulants, C n , for net-proton, proton, and antiproton multiplicity distributions, and correlation functions, κ n , for proton and antiproton multiplicity distributions up to the fourth order in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 54.4, 62.4, and 200 GeV. The C n and κ n are presented as a function of collision energy, centrality and kinematic acceptance in rapidity, y, and transverse momentum, p T . The data were taken during the first phase of the Beam Energy Scan (BES) program (2010–2017) at the BNL Relativistic Heavy Ion Collider (RHIC) facility. The measurements are carried out at midrapidity (|y| < 0.5) and transverse momentum 0.4 < p T < 2.0 GeV/c, using the STAR detector at RHIC. We observe a nonmonotonic energy dependence ($\sqrt{s_{NN}}$ = 7.7–62.4 GeV) of the net-proton C 4 /C 2 with the significance of 3.1σ for the 0–5% central Au+Au collisions. This is consistent with the expectations of critical fluctuations in a QCD-inspired model. Thermal and transport model calculations show a monotonic variation with $\sqrt{s_{NN}}$. For the multiparticle correlation functions, we observe significant negative values for a two-particle correlation function, κ 2 , of protons and antiprotons, which are mainly due to the effects of baryon number conservation. Furthermore, it is found that the four-particle correlation function, κ 4 , of protons plays a role in determining the energy dependence of proton C 4 /C 1 below 19.6 GeV, which cannot be understood by the effect of baryon number conservation.
Not Available
Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.
Standard relativistic hydrodynamics has been successful in describing the properties of the strongly interacting matter produced in the heavy-ion collision experiments. Recently, there has been a significant theoretical advancement in this field to explain spin polarization of hadrons emitted in these processes. Although current models have successfully explained some of the experimental data based on the coupling between spin polarization and vorticity of the medium, they still lack a clear understanding of the differential measurements. This is commonly interpreted as an indication that the spin needs to be treated as an independent degree of freedom whose dynamics is not entirely bound to flow circulation. In particular, if the spin is a macroscopic property of the system, in equilibrium its dynamics should follow hydrodynamic laws. Here, we develop a framework of relativistic hydrodynamics which includes spin degrees of freedom from the quantum kinetic theory for Dirac fermions and use it for modeling the dynamics of matter. Following experimental observations, we assume that the polarization effects are small and derive conservation laws for the net baryon current, the energy–momentum tensor and the spin tensor based on the de Groot–van Leeuwen–van Weert definitions of these currents. We present various properties of the spin polarization tensor and its components, analyze the propagation properties of the spin polarization components, and derive the spin-wave velocity for arbitrary statistics. We find that only the transverse spin components propagate, analogously to the electromagnetic waves. Finally, using our framework, we study the space–time evolution of the spin polarization for the systems respecting certain space–time symmetries and calculate the mean spin polarization per particle, which can be compared to the experimental data. We find that, for some observables, our spin polarization results agree qualitatively with the experimental findings and other model calculations.
Here, the multiphase transport model (AMPT) is used to study the final-state effects on the symmetric correlations (SC), asymmetric correlations (ASC), normalized symmetric correlations (NSC), and normalized asymmetric correlations (NASC) in Au+Au collisions at 200 GeV. The correlators' sensitivity to nonflow effects associated with long- and short-range nonflow correlations are also discussed using the HIJING model. The results indicate that SC, ASC, NSC, and NASC can give accompanying constraints for initial- and final-state effects. In addition, conducting further detailed experimental measurements spanning a broad range of collision systems and beam energies will serve as an additional constraint for the theoretical models' calculations.
We present calculations of bulk properties and multiparticle correlations in a large variety of collision systems within a hybrid formalism consisting of IP-Glasma initial conditions, Music viscous relativistic hydrodynamics, and UrQMD microscopic hadronic transport. In particular, we study heavy ion collisions at the Large Hadron Collider (LHC), including Pb + Pb, Xe + Xe, and O + O collisions, and Au + Au, U + U, Ru + Ru, Zr + Zr, and O + O collisions at the Relativistic Heavy Ion Collider (RHIC). We further study asymmetric systems, including p + Au, d + Au, 3 He + Au, and p + Pb collisions at various energies as well as p+p collisions at 0.5 and 13 TeV. Here, we describe experimental observables in all heavy ion systems well with one fixed set of parameters, validating the energy and system dependence of the framework. As a result, many observables in the smaller systems are also well described, although they test the limits of the model.
Recently studies of the differential nature of the flow angle fluctuations, known as event plane angular decorrelation, indicated that measurements that assume a common symmetry plane may need to consider the flow angle fluctuations effect. Using the HIJING and AMPT models, it is shown that the flow angle fluctuations measurements, obtained with the two- and three-subevents correlation method, could have significant nonflow effects associated with long- and short-range nonflow correlations. The current study demonstrates the four-subevents cumulant methods' ability to reduce the nonflow effects. Finally, it is further argued that the measurements using the four-subevents correlation method can be used to provide accurate quantification of the differential flow angle fluctuations.
Proxies for cumulants of baryon number 𝐵, electric charge 𝑄, and strangeness 𝑆 are usually measured in heavy-ion collisions via moments of net-number distribution of given hadronic species. Since these cumulants of conserved charges are expected to be sensitive to the existence of a critical point in the phase diagram of nuclear matter, it is crucial to ensure that the proxies used as substitutes are as close to them as possible. Hence, we use the EPOS 4 framework to generate Au + Au collisions at several collision energies of the BNL Relativistic Heavy Ion Collider beam energy scan. We compute second-order net cumulants of 𝜋, 𝐾, and 𝑝, for which experimental data have been published as well as the corresponding conserved charge cumulants. We then compare them with proxies, defined in previous lattice QCD and hadron resonance gas model studies, which are shown to reproduce more accurately their associated conserved charge cumulants. We investigate the impact of hadronic rescatterings occurring in the late evolution of the system on these quantities, as well as the amount of signal actually originating from the bulk medium which endures a phase transition.
Here, we introduce additional coalescence factors for the production of strange baryons in a multiphase transport (AMPT) model in order to describe the enhanced production of multistrange hadrons observed in Pb + Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV at the Large Hadron Collider (LHC) and Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV at Relativistic Heavy-Ion Collider (RHIC). This extended AMPT model is found to also give a reasonable description of the multiplicity dependence of the strangeness enhancement observed in high multiplicity events in $\textit{pp}$ collisions at $\sqrt{s}$ = 7 TeV and $\textit{p}$-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. We find that the coalescence factors depend on the system size but not much on whether the system is produced from A + A or $\textit{p}$ + A collisions. Here, the extended AMPT model thus provides a convenient way to model the mechanism underlying the observed strangeness enhancement in collisions of both small and large systems at RHIC and LHC energies.
This work presents a model for generating nonthermal power-law tails of particles’ energy probability density functions in turbulent collisionless plasmas, applicable to both nonrelativistic and relativistic scenarios. We propose that strong Alfvénic turbulence energizes plasma particles through curvature acceleration, particularly for particles with Larmor radii comparable to the scales of turbulence. When the energy density of the energized particles increases, the efficiency of the energy exchange process diminishes. As a result, the acceleration process saturates, leading to power-law distributions of particle momentum and energy. In the nonrelativistic case, the momentum probability density function scales as f(p)dp ∝ p −3 dp, while in the ultrarelativistic case, the energy probability density function scales as f(γ)dγ ∝ γ −3 dγ, where γ is the Lorentz factor. This model provides a unified framework for understanding particle acceleration in both energy regimes, complementing existing analytical approaches. The predicted scalings are consistent with available observations of energetic ion distributions in the heliosphere and with the findings from numerical simulations of ultrarelativistic particle acceleration in magnetically dominated plasma turbulence.
Measurements from the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) can be used to study the properties of quark-gluon plasma. Systematic constraints on these properties must combine measurements from different collision systems and methodically account for experimental and theoretical uncertainties. Such studies require a vast number of costly numerical simulations. While computationally inexpensive surrogate models (“emulators”) can be used to efficiently approximate the predictions of heavy ion simulations across a broad range of model parameters, training a reliable emulator remains a computationally expensive task. We use transfer learning to map the parameter dependencies of one model emulator onto another, leveraging similarities between different simulations of heavy ion collisions. By limiting the need for large numbers of simulations to only one of the emulators, this technique reduces the numerical cost of comprehensive uncertainty quantification when studying multiple collision systems and exploring different models.
Particle correlations are a powerful tool to study the properties of the bulk nuclear matter produced in relativistic heavy ion collisions. The momentum correlations between identical particles originating from the same particle-emitting source, referred to as the Bose-Einstein correlations, measure scales that are related to the geometrical size of the source. The two-particle azimuthal angular correlations measure the spatial anisotropy of produced particles, providing information on collective phenomena arising in the dense nuclear medium. This contribution will discuss new LHCb measurements of Bose- Einstein correlations and, for the first time, the collective flow coefficients in the far forward rapidity region.
The chiral magnetic effect (CME) refers to charge separation along a strong magnetic field due to imbalanced chirality of quarks in local parity and charge-parity violating domains in quantum chromodynamics. The experimental measurement of the charge separation is made difficult by the presence of a major background from elliptic azimuthal anisotropy. This background and the CME signal have different sensitivities to the spectator and participant planes, and could thus be determined by measurements with respect to these planes. We report such measurements in Au+Au collisions at a nucleon-nucleon center-of-mass energy of 200 GeV at the Relativistic Heavy-Ion Collider. It is found that the charge separation, with the flow background removed, is consistent with zero in peripheral (large impact parameter) collisions. Some indication of finite CME signals is seen with a significance of 1–3 standard deviations in mid-central (intermediate impact parameter) collisions. Furthermore, significant residual background effects may, however, still be present.
The yields of hadrons and light nuclei in relativistic collisions of heavy-nuclei at a center of mass energy of $\sqrt{s_{NN}}$ = 2.6 TeV can be described remarkably well by a thermal distribution of an ideal gas of hadrons and light nuclei interacting only via the decay of resonances. Given the particularly small binding energy of hypertritons relative to the temperature describing the yields (about 156 MeV), one might naturally expect hypertrions to dissociate in medium, making the agreement of hypertriton yields with thermal predictions highly puzzling. The puzzle is compounded by the fact that small binding energy is associated with the large size of the hypertriton. This size is on a similar scale to the overall size of the fireball and much larger than the length scale over which temperatures in the fireball vary over phenomenologically relevant amounts. Here, this paper quantifies the tension this effect causes and shows that it is sufficiently large to render the thermal model inconsistent: its natural assumptions are in conflict with its outputs. The possibility that hypertritons are formed at freeze out as compact objects, quark droplets, that subsequently evolve into hypertritons is considered as a way to resolve the puzzle. It is noted that beyond making the assumption that compact quark droplets form, additional detailed dynamical assumptions which have not been justified are needed to make the thermal model work. The issue of why, despite these issues, the hypertriton is well described by a simple statistical description at freeze out is unresolved. Resolving the hypertriton puzzle is important as it may clarify whether the phenomenological success of the simple thermal model for yields accurately reflects the simple picture of the underlying physics on which it is based.
We report a new determination of $\hat{q}$, the jet transport coefficient of the quark-gluon plasma. We use the JETSCAPE framework, which incorporates a novel multistage theoretical approach to in-medium jet evolution and Bayesian inference for parameter extraction. The calculations, based on the Matter and Lbt jet quenching models, are compared to experimental measurements of inclusive hadron suppression in Au + Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and Pb + Pb collisions at the CERN Large Hadron Collider (LHC). The correlation of experimental systematic uncertainties is accounted for in the parameter extraction. The functional dependence of $\hat{q}$ on jet energy or virtuality and medium temperature is based on a perturbative picture of in-medium scattering, with components reflecting the different regimes of applicability of Matter and Lbt. In the multistage approach, the switch between Matter and Lbt is governed by a virtuality scale Q 0 . Comparison of the posterior model predictions to the RHIC and LHC hadron suppression data shows reasonable agreement, with moderate tension in limited regions of phase space. The distribution of $\hat{q}/T^3$ extracted from the posterior distributions exhibits weak dependence on jet momentum and medium temperature T, with 90% credible region (CR) depending on the specific choice of model configuration. The choice of Matter + Lbt, with switching at virtuality Q 0 , has 90% CR of $2<\hat{q}/T^3<4$ for p T jet > 40 GeV=c. The value of Q 0 , determined here for the first time, is in the range 2.0-2.7 GeV.
In this study, the beam energy scan (BES) program at the Relativistic Heavy Ion Collider (RHIC) was extended to energies below $\sqrt{s_{NN}}$ = 7.7 GeV in 2015 by successful implementation of the fixed-target mode of operation in the STAR (Solenoidal Tracker At RHIC) experiment. In this mode, ions circulate in one ring of the collider and interact with a stationary target at the entrance of the STAR Time Projection Chamber. The first results for Au + Au collisions at $\sqrt{s_{NN}}$ = 4.5 GeV are presented, demonstrating good performance of all the relevant detector subsystems in fixed-target mode. Results presented here include directed and elliptic flow of identified hadrons, and radii from pion femtoscopy. The latter, together with recent HADES results, reveal a long-sought peak structure that may be caused by the system evolving through a first-order phase transition from quark-gluon plasma to the hadronic phase. Directed and elliptic flow for pions are presented for the first time at this beam energy. Pion and proton elliptic flow show behavior which hints at constituent quark scaling, and demonstrate that a definitive conclusion will be achievable using the full statistics of the on-going second phase of BES (BES-II). In particular, BES-II to date has recorded fixed-target data sets with two orders of magnitude more events at each of nine energies between $\sqrt{s_{NN}}$ = 3.0 and 7.7 GeV.