Engineering topics

Jia, Jiangyong

Publications and source records attributed to Jia, Jiangyong.

At least 19 records

DOE OSTI2023

Dense nuclear matter equation of state from heavy-ion collisions

The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. Furthermore, this white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Probe nuclear structure using the anisotropic flow at the Large Hadron Collider

Abstract Recent studies have shown that the shape and radial profile of the colliding nuclei have strong influences on the initial condition of the heavy ion collisions and the subsequent development of the anisotropic flow. Using A Multi-Phase Transport model (AMPT) model, we investigated the impact of nuclear quadrupole deformation $$\beta _2$$ β 2 and nuclear diffuseness $$a_0$$ a 0 of $$^{129}$$ 129 Xe on various of flow observables in Xe–Xe collisions at $$\sqrt{s_{\textrm{NN}}}=$$ s NN = 5.44 TeV. We found that $$\beta _2$$ β 2 has a strong influence on central collisions while $$a_0$$ a 0 mostly influences the mid-central collisions. The relative change of flow observables induced by a change in $$\beta _2$$ β 2 and $$a_0$$ a 0 are also found to be insensitive to the values of parameters controlling the strength of the interaction among final state particles. Our study demonstrates the potential for constraining the initial condition of heavy ion collisions using future system scans at the LHC.

Physics↗

DOE OSTI2023

Separating the Impact of Nuclear Skin and Nuclear Deformation in High-Energy Isobar Collisions

Bulk nuclear structure properties, such as radii and deformations, leave distinct signatures in the final state of relativistic heavy-ion collisions. Isobaric collisions offer an easy route to establish explicit connections between the colliding nuclei's structure and the observable outcomes. Here, we investigate the effects of nuclear skin thickness and nuclear deformations on the elliptic flow (v 2 ) and its fluctuations in high-energy 96 Ru + 96 Ru and 96 Zr + 96 Zr collisions. Furthermore, our findings reveal that the difference in skin thickness between these isobars only influences the inherent ellipticity of the collision systems, $v_2^{rp}$. In contrast, differences in nuclear deformations solely impact the fluctuations of v 2 around $v_2^{rp}$. Hence, we have identified a data-driven method to disentangle the effects of nuclear skin and nuclear deformations, marking a significant step toward assessing the consistency of nuclear phenomena across energy scales.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Precision Tests of the Nonlinear Mode Coupling of Anisotropic Flow via High-Energy Collisions of Isobars

Valuable information on dynamics of expanding fluids can be inferred from the response of such systems to perturbations in their initial geometry. We apply this technique in high-energy 96 Ru+ 96 Ru and 96 Zr+ 96 Zr collisions to scrutinize the expansion dynamics of the quark-gluon plasma, where the initial geometry perturbations are sourced by the differences in deformations and radial profiles between 96 Ru and 96 Zr, and the collective response is captured by the change in anisotropic flow V n between the two collision systems. Using a transport model, we analyze how the nonlinear coupling between lower-order flow harmonics V 2 and V 3 to the higher-order flow harmonics V 4 and V 5 , expected to scale as V 4 NL = χ 4 V 2 2 and V 5NL = χ 5 V 2 V 3 , gets modified as one moves from 96 Ru+ 96 Ru to 96 Zr+ 96 Zr systems. We find that these scaling relations are valid to high precision: variations of order 20% in V 4NL and V 5NL due to differences in quadrupole deformation, octupole deformation, and nuclear skin modify χ 4 and χ 5 by about 1–2%. Percent-level deviations are however larger than the expected experimental uncertainties and could be measured. Therefore, collisions of isobars with different nuclear structures are a unique tool to isolate subtle nonlinear effects in the expansion of the quark-gluon plasma that would be otherwise impossible to access in a single collision system.

Physics↗

DOE OSTI2023

Scaling approach to nuclear structure in high-energy heavy-ion collisions

In high-energy heavy-ion collisions, the initial condition of the produced quark-gluon plasma (QGP) and its evolution are sensitive to collective nuclear structure parameters describing the shape and radial profiles of the nuclei. Here, we find a general scaling relation between these parameters and many experimental observables such as elliptic flow, triangular flow, and particle multiplicity distribution. In particular, the ratios of observables between two isobar systems depend only on the differences of these parameters, but not on the details of the final state interactions, hence offering a new way to constrain the QGP initial condition. Using this scaling relation, we show how the structure parameters of $^{96}_{44}$Ru and $^{96}_{40}$Zr conspire to produce the rich centrality dependences of these ratios, as measured by the STAR Collaboration. Our scaling approach demonstrates that isobar collisions are a precision tool to probe the initial condition of heavy-ion collisions, as well as the collective nuclear structures, including the neutron skin, of the atomic nuclei across energy scales.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Free spectator nucleons in ultracentral relativistic heavy-ion collisions as a probe of neutron skin

Besides the yield ratio of free spectator neutrons produced in ultracentral 96 Zr + 96 Zr to 96 Ru + 96 Ru collisions, we propose that the yield ratio N n /N p of free spectator neutrons to protons in a single collision system at the BNL Relativistic Heavy Ion Collider or the CERN Large Hadron Collider can be a more sensitive probe of the neutron-skin thickness Δrnp and the slope parameter L of the symmetry energy. Here, the idea is demonstrated based on the proton and neutron density distributions of colliding nuclei obtained from Skyrme-Hartree-Fock-Bogolyubov calculations, and on a Glauber model that provides information of spectator matter. The final spectator particles are produced from direct emission, clusterization by a minimum spanning tree algorithm or a Wigner function approach, and deexcitation of heavy clusters by gemini. A larger Δr np associated with a larger L value increases the isospin asymmetry of spectator matter and thus leads to a larger N n /N p , especially in ultracentral collisions where the multiplicity of free nucleons are free from the uncertainties of cluster formation and deexcitation. We have further shown that the double ratio of N n /N p in isobaric collision systems or in collisions by isotopes helps to cancel the detecting efficiency for protons. Effects from nuclear deformation and electromagnetic excitation are studied, and they are found to be subdominant compared to the expected sensitivity to Δr np .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Probing neutron-skin thickness with free spectator neutrons in ultracentral high-energy isobaric collisions

We show that the yield ratio of free spectator neutrons produced in high-energy 96 Zr+ 96 Zr to 96 Ru+ 96 Ru collisions is a clean probe of the neutron-skin thickness of colliding nuclei and the slope parameter L of the symmetry energy. The idea is demonstrated based on the proton and neutron density distributions via a state-of-the-art Skyrme-Hartree-Fock-Bogolyubov calculation. Among spectator nucleons given by the Glauber model, free spectator neutrons include those from direct production that survive from clusterization as well as those from deexcitation of heavy clusters described by the popular GEMINI model. More free neutrons are produced in collisions of 96 Zr nucleus due to its larger neutron skin, compared to those produced in collisions of 96 Ru nucleus with a smaller neutron skin. The difference of the free spectator neutron yield is further increased with the increasing difference of the neutron-skin thickness between 96 Zr and 96 Ru with a larger L value, and the increase in ultracentral collisions is particularly insensitive to model details and experimental uncertainties. Since the production of free spectator neutrons is not affected by the complicated dynamics in the mid-rapidity region, the ratio of their multiplicities in ultracentral isobaric collisions is a robust observable for constraining the neutron skin and L value.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Ratios of collective flow observables in high-energy isobar collisions are insensitive to final-state interactions

The ratios of bulk observables, such as harmonic flow v 2 and v 3 , between high-energy 96 Ru + 96 Ru and 96 Zr + 96 Zr collisions were recently argued to be a clean probe of the nuclear structure differences between 96 Ru and 96 Zr . Using a transport model simulation of isobar collisions, we quantify this claim from the dependence of the ratios v 2,Ru /v 2, Zr and v 3,Ru /v 3,Zr on various final-state effects, such as the shear viscosity, hadronization, and hadronic cascade. Although the v 2 and v 3 change by more than 50% when varying the final-state effects, the ratios are unchanged. In addition, these ratios are independent of the transverse momentum p T and hadron species, despite of up to a factor of 2 change in v n . Further, the ratio of mean transverse momentum $\langle$p T $\rangle$ is found to be controlled by the nuclear skin and nuclear radius but is only slightly impacted by the final-state effects. Therefore, these isobar ratios serve as a clean probe of the initial condition of the quark-gluon plasma, which, in turn, is controlled by the collective structure of the colliding nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Phenomenological study on correlation between flow harmonics and mean transverse momentum in nuclear collisions

To assess the properties of the quark-gluon plasma formed in nuclear collisions, the Pearson correlation coefficient between flow harmonics and mean transverse momentum, \rho\left(v_{n}^{2},\left[p_{\mathrm{T}}\right]\right) ρ ( v n 2 , [ p T ] ) , reflecting the overlapped geometry of colliding atomic nuclei, is measured. \rho\left(v_{2}^{2},\left[p_{\mathrm{T}}\right]\right) ρ ( v 2 2 , [ p T ] ) was found to be particularly sensitive to the quadrupole deformation of the nuclei. We study the influence of the nuclear quadrupole deformation on \rho\left(v_{n}^{2},\left[p_{\mathrm{T}}\right]\right) ρ ( v n 2 , [ p T ] ) in Au+Au and U+U collisions at RHIC energy using AMPT transport model, and show that the \rho\left(v_{2}^{2},\left[p_{\mathrm{T}}\right]\right) ρ ( v 2 2 , [ p T ] ) is reduced by the quadrupole deformation \beta_2 β 2 and turns to change sign in ultra-central collisions (UCC).

Zhang, Chunjian↗

DOE OSTI2022

Probing triaxial deformation of atomic nuclei in high-energy heavy ion collisions

Most atomic nuclei are deformed with a quadrupole shape described by its overall strength β 2 and triaxiality γ. The deformation can be accessed in high-energy heavy ion collisions by measuring the collective flow response of the produced quark-gluon plasma to the eccentricity ε 2 and the density gradient d ⊥ in the initial state. Using an analytical estimate and a Glauber model, I show that the variances $\langleε^{2}_{2}\rangle$ or $\langle$(δd ⊥ /d ⊥ ) 2 $\rangle$ and skewnesses $\langleε^{2}_{2}δd_{⊥}/d_{⊥}\rangle$ or $\langle(δd_{⊥}/d_{⊥})^{3}\rangle$ have a simple analytical form of a'+b'β$^{2}_{2}$ and a'+[b'+c' cos(3γ)]β$^{3}_{2}$, respectively. From these, I constructed several normalized skewnesses to isolate the γ dependence from that of β 2 and show that the correlations between a normalized skewness and a variance can constrain simultaneously β 2 and γ. Assuming a linear relation with elliptic flow v 2 and mean-transverse momentum [p T ] of final-state particles, v 2 ∝ ε 2 and δ [p T ] / [p T ] ∝ δd ⊥ /d ⊥ , similar conclusions are also expected for the variances and skewnesses of v 2 and [p T ], i.e., a + bβ$^{2}_{2}$ for $\langle v^{2}_{2}\rangle$ and $\langle(δ[p_{T}] / [p_{T}])^{2}\rangle$ and a + [b + c cos (3γ)] β$^{3}_{2}$ for $\langle v^{2}_{2}δ [p_{T}] / [p_{T}]\rangle$ or $\langle(δ[p_{T}] / [p_{T}])^{3}\rangle$. My findings motivate a dedicated system scan of high-energy heavy ion collisions at RHIC and LHC to measure triaxiality of atomic nuclei: one first determines the coefficients b and c by collisions of isobaric near prolate nuclei, cos (3γ) ≈ 1, and near oblate nuclei, cos (3γ) ≈ –1 , with known β 2 values, followed by collisions of other species of interest with similar mass number. The (β 2 , γ) values for this species can be inferred directly from the measured variance and skewness observables from these collisions. Finally, the results demonstrate the unique opportunities offered by high-energy collisions as a tool to perform interdisciplinary nuclear physics studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Higher-order transverse momentum fluctuations in heavy-ion collisions

In relativistic heavy-ion collisions, the event-by-event mean transverse momentum fluctuations are sensitive to overlap area and energy density fluctuations in the initial state. We present a framework to calculate p T fluctuations up to fourth order using standard and subevent methods, which is validated using the HIJING model. We observe a power-law dependence for cumulants of all orders as a function of charged particle multiplicity N ch , consistent with a simple independent source picture. The fluctuation in pp collisions is observed to be larger than for p+Pb, Pb + Pb, and Xe + Xe collisions at the same N ch due to bias in the number of contributing sources. The short-range correlations are greatly suppressed in the subevent method in comparison to calculations based on the standard method. This paper provides a baseline for transverse momentum fluctuations without the presence of final-state effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Proceedings of RIKEN BNL Research Center Workshop: Physics Opportunities from the RHIC Isobar Run (Volume 138) [Slides]

This workshop “Physics Opportunities from the RHIC Isobar Run” was held virtually due to the COVID-19 pandemic, with 102 registered participants. We gathered experts from both the experimental and theoretical community to advance our understanding of the physics potential of the isobar dataset by discussing analysis strategies and theoretical predictions. The goal was to explore and identify interesting observables and associated physics which the Isobar Runs can provide. We have identified possible backgrounds and determined the required experimental precisions. Finally, we called for theoretical predictions and interesting isobar species for possible future experiments. We appreciate the effort of, and would like to congratulate, the organizers for conducting a very successful workshop. We also welcome and thank the world-wide participants who attended.

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

Evidence of Quadrupole and Octupole Deformations in Zr 96 + Zr 96 and Ru 96 + Ru 96 Collisions at Ultrarelativistic Energies

In the hydrodynamic model description of heavy-ion collisions, the elliptic flow v 2 and triangular flow v 3 are sensitive to the quadrupole deformation β 2 and octupole deformation β 3 of the colliding nuclei. The relations between v n and β n have recently been clarified and were found to follow a simple parametric form. The STAR Collaboration has just published precision v n data from isobaric Ru 96 + Ru 96 and Zr 96 + Zr 96 collisions, where they observe large differences in central collisions v 2 , Ru > v 2 , Zr and v 3 , Ru < v 3 , Zr . Using a transport model simulation, we show that these orderings are a natural consequence of $\beta_{2,Ru} \gg \beta_{2,Zr}$ and $\beta_{3,Ru} \ll \beta_{3,Zr}$ We reproduce the centrality dependence of the v 2 ratio qualitatively and v 3 ratio quantitatively and extract values of β 2 and β 3 that are consistent with those measured at low-energy nuclear structure experiments. STAR data provide the first direct evidence of strong octupole correlations in the ground state of Zr 96 in heavy-ion collisions. Our analysis demonstrates that flow measurements in high-energy, heavy-ion collisions, especially using isobaric systems, are a new precision tool to study nuclear structure physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Probing nuclear quadrupole deformation from correlation of elliptic flow and transverse momentum in heavy ion collisions

In heavy ion collisions, elliptic flow v2 and radial flow, characterized by eventwise average transverse momentum [p T ], are related to the shape and size of the overlap region, which are sensitive to the shape of colliding atomic nuclei. The Pearson correlation coefficient between v 2 and [p T ], ρ 2 , was found to be particularly sensitive to the quadrupole deformation parameter β that is traditionally measured in low energy experiments. Built on earlier insight that the prolate deformation β>0 reduces the ρ 2 in ultracentral collisions (UCC), we show that the prolate deformation β<0 enhances the value of ρ 2 . As β>0 and β<0 are the two extremes of triaxiality, the strength and sign of v 2 2 -[pT] correlation can be used to provide valuable information on the triaxiality of the nucleus. Our study provide further arguments for using the hydrodynamic flow as a precision tool to directly image the deformation of the atomic nuclei at extremely short timescale (<10 -24 s).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Shape of atomic nuclei in heavy ion collisions

In the hydrodynamic model description of heavy ion collisions, the final-state anisotropic flows v n are linearly related to the strengths of the multipole shape of the distribution of nucleons in the transverse plane, ε n : v n ∝ε n . The ε n , for n=1, 2, 3, 4, are sensitive to the shapes of the colliding ions, characterized by the quadrupole β 2 , octupole β 3 , and hexadecapole β 4 deformations. This sensitivity is investigated analytically and also in a Monte Carlo Glauber model. One observes a robust linear relation, $\langleε^{2}_{n}\rangle$ = $a^{'}_{n} + b^{'}_{n}β^{2}_{n}$, for events in a fixed centrality. The $\langleε^{2}_{1}\rangle$ has a contribution from β 3 and β 4 , and $\langleε^{2}_{3}\rangle$ from β 4 . In ultracentral collisions, there are little cross contributions between β 2 and ε 3 and between β 3 and ε 2 , but clear cross contributions are present in noncentral collisions. Additionally, $\langleε^{2}_{n}\rangle$ are insensitive to nonaxial shape parameters such as the triaxiality. This is good news because the measurements of v 2 , v 3 , and v 4 can be used to constrain simultaneously the β 2 , β 3 , and β 4 values. This is best done by comparing two colliding ions with similar mass numbers and therefore nearly identical $a^{'}_{n}$, to obtain a simple equation that relates the β n of the two species. Finally, this opens up the possibility to map the shape of the atomic nuclei at a timescale (<10 –24 s) much shorter than probed by low-energy nuclear structure physics (<10 –21 s), which ultimately may provide information complementary to that obtained in the nuclear structure experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Impact of Nuclear Deformation on Relativistic Heavy-Ion Collisions: Assessing Consistency in Nuclear Physics across Energy Scales

In the hydrodynamic framework of heavy-ion collisions, elliptic flow v 2 is sensitive to the quadrupole deformation β of the colliding ions. This enables one to test whether the established knowledge on the low-energy structure of nuclei is consistent with collider data from high-energy experiments. We derive a formula based on generic scaling laws of hydrodynamics to relate the difference in v 2 measured between collision systems that are close in size to the value of β of the respective species. We validate our formula in simulations of 238 U + 238 U and 197 Au + 197 Au collisions at top Relativistic Heavy Ion Collider (RHIC) energy, and subsequently apply it to experimental data. Using the deformation of 238 U from low-energy experiments, we find that RHIC v 2 data implies 0.16 ≲ | β | ≲ 0.20 for 197 Au nuclei, i.e., significantly more deformed than reported in the literature, posing an interesting issue in nuclear phenomenology.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Investigation of experimental observables in search of the chiral magnetic effect in heavy-ion collisions in the STAR experiment *

The chiral magnetic effect (CME) is a novel transport phenomenon, arising from the interplay between quantum anomalies and strong magnetic fields in chiral systems. In high-energy nuclear collisions, the CME may survive the expansion of the quark-gluon plasma fireball and be detected in experiments. Over the past two decades, experimental searches for the CME have attracted extensive interest at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). The main goal of this study is to investigate three pertinent experimental approaches: the $γ$ correlator, the R correlator, and the signed balance functions. We exploit simple Monte Carlo simulations and a realistic event generator (EBE-AVFD) to verify the equivalence of the core components among these methods and to ascertain their sensitivities to the CME signal and the background contributions for the isobar collisions at the RHIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Accessing the shape of atomic nuclei with relativistic collisions of isobars

Nuclides sharing the same mass number (isobars) are observed ubiquitously along the stability line. While having nearly identical radii, stable isobars can differ in shape, and present different quadrupole deformations. We show that even small differences in these deformations can be probed by relativistic nuclear collisions experiments, where they manifest as deviations from unity in the ratios of elliptic flow coefficients taken between isobaric systems. Collider experiments with isobars represent, thus, a unique means to gain precise knowledge of the geometric shape of atomic nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗