Showing 1–50 of 54 results for author: Du, L

Probing the equilibration of the QCD matter created in heavy-ion collisions with dileptons

Authors: Xiang-Yu Wu, Lipei Du, Charles Gale, Sangyong Jeon

Abstract: A systematic study of intermediate invariant mass dilepton production in Pb+Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV is performed, using next-to-leading-order (NLO) thermal QCD dilepton emission rates with a multistage dynamical approach which includes event-by-event IP-Glasma initial conditions, relativistic viscous fluid dynamics, and a hadronic afterburner. Considering dilepton yield and ani… ▽ More A systematic study of intermediate invariant mass dilepton production in Pb+Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV is performed, using next-to-leading-order (NLO) thermal QCD dilepton emission rates with a multistage dynamical approach which includes event-by-event IP-Glasma initial conditions, relativistic viscous fluid dynamics, and a hadronic afterburner. Considering dilepton yield and anisotropic flow, special attention is paid to the out-of-equilibrium aspects, both thermal and chemical, and to the contribution of the Drell-Yan process. The relative contribution of each of those different channels to dilepton observables is calculated and discussed. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 13 pages, 10 figures

The QCD phase diagram and Beam Energy Scan physics: a theory overview

Authors: Lipei Du, Agnieszka Sorensen, Mikhail Stephanov

Abstract: 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. 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. △ Less

Submitted 6 May, 2024; v1 submitted 15 February, 2024; originally announced February 2024.

Comments: v1: 105 pages + 33 pages references, 24 figures; v2: 142 pages. A chapter contribution to book Quark-Gluon Plasma 6

Report number: INT-PUB-24-017

Journal ref: Int. J. Mod. Phys. E (2024) 2430008

Photon-triggered jets as probes of multi-stage jet modification

Authors: C. Sirimanna, Y. Tachibana, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kordell II, A. Kumar, R. Kunnawalkam-Elayavalli, J. Latessa, S. Lee , et al. (28 additional authors not shown)

Abstract: Prompt photons are created in the early stages of heavy ion collisions and traverse the QGP medium without any interaction. Therefore, photon-triggered jets can be used to study the jet quenching in the QGP medium. In this work, photon-triggered jets are studied through different jet and jet substructure observables for different collision systems and energies using the JETSCAPE framework. Since t… ▽ More Prompt photons are created in the early stages of heavy ion collisions and traverse the QGP medium without any interaction. Therefore, photon-triggered jets can be used to study the jet quenching in the QGP medium. In this work, photon-triggered jets are studied through different jet and jet substructure observables for different collision systems and energies using the JETSCAPE framework. Since the multistage evolution used in the JETSCAPE framework is adequate to describe a wide range of experimental observables simultaneously using the same parameter tune, we use the same parameters tuned for jet and leading hadron studies. The same isolation criteria used in the experimental analysis are used to identify prompt photons for better comparison. For the first time, high-accuracy JETSCAPE results are compared with multi-energy LHC and RHIC measurements to better understand the deviations observed in prior studies. This study highlights the importance of multistage evolution for the simultaneous description of experimental observables through different collision systems and energies using a single parameter tune. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 4 pages, 5 figures. Proceedings of Quark Matter 2023 - XXXth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions, Houston, TX, 3-9 September 2023

Measuring jet quenching with a Bayesian inference analysis of hadron and jet data by JETSCAPE

Authors: R. Ehlers, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, L. Du, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kordell II, A. Kumar, R. Kunnawalkam-Elayavalli, J. Latessa, S. Lee, Y. -J. Lee, D. Liyanage , et al. (28 additional authors not shown)

Abstract: The JETSCAPE Collaboration reports the first multi-messenger study of the QGP jet transport parameter $\hat{q}$ using Bayesian inference, incorporating all available hadron and jet inclusive yield and jet substructure data from RHIC and the LHC. The theoretical model utilizes virtuality-dependent in-medium partonic energy loss coupled to a detailed dynamical model of QGP evolution. Tension is obse… ▽ More The JETSCAPE Collaboration reports the first multi-messenger study of the QGP jet transport parameter $\hat{q}$ using Bayesian inference, incorporating all available hadron and jet inclusive yield and jet substructure data from RHIC and the LHC. The theoretical model utilizes virtuality-dependent in-medium partonic energy loss coupled to a detailed dynamical model of QGP evolution. Tension is observed when constraining $\hat{q}$ for different kinematic cuts of the inclusive hadron data. The addition of substructure data is shown to improve the constraint on $\hat{q}$, without inducing tension with the constraint due to inclusive observables. These studies provide new insight into the mechanisms of jet interactions in matter, and point to next steps in the field for comprehensive understanding of jet quenching as a probe of the QGP. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 4 pages, 2 figures. Proceedings of Quark Matter 2023 - XXXth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions, Houston, TX, 3-9 September 2023

Bulk medium properties of heavy-ion collisions from the beam energy scan with a multistage hydrodynamic model

Authors: Lipei Du

Abstract: We introduce a method to reconstruct full rapidity distributions of charged particle multiplicity and net proton yields, crucial for constraining the longitudinal dynamics of nuclear matter created in the beam energy scan program. Employing rapidity distributions within a multistage hydrodynamic model calibrated for Au+Au collisions at $\sqrt{s_\mathrm{NN}}=7.7-200\,$GeV, we estimate the total ene… ▽ More We introduce a method to reconstruct full rapidity distributions of charged particle multiplicity and net proton yields, crucial for constraining the longitudinal dynamics of nuclear matter created in the beam energy scan program. Employing rapidity distributions within a multistage hydrodynamic model calibrated for Au+Au collisions at $\sqrt{s_\mathrm{NN}}=7.7-200\,$GeV, we estimate the total energy and baryon number deposited into the collision fireball, offering insights into initial dynamics and the identification of nuclear remnants. We explore the potential of rapidity-dependent measurements in probing equations of state at finite chemical potentials. Furthermore, we compare the freeze-out parameters derived from both hydrodynamics and thermal models, highlighting that the parameters extracted via thermal models represent averaged properties across rapidities. △ Less

Submitted 8 July, 2024; v1 submitted 31 December, 2023; originally announced January 2024.

Comments: v1: 15 pages, 12 figures; v2: figs. 2 and 4 slightly adjusted; v3: Title slightly changed. A few short discussions were added, but there are no changes to the results or conclusions. Published in PRC

Journal ref: Phys. Rev. C 110, 014904 (2024)

3D Multi-system Bayesian Calibration with Energy Conservation to Study Rapidity-dependent Dynamics of Nuclear Collisions

Authors: Andi Mankolli, Aaron Angerami, Ritu Arora, Steffen Bass, Shanshan Cao, Yi Chen, Lipei Du, Raymond Ehlers, Hannah Elfner, Wenkai Fan, Rainer J. Fries, Charles Gale, Yayun He, Ulrich Heinz, Barbara Jacak, Peter Jacobs, Sangyong Jeon, Yi Ji, Lauren Kasper, Michael Kordell II, Amit Kumar, R. Kunnawalkam-Elayavalli, Joseph Latessa, Sook H. Lee, Yen-Jie Lee , et al. (26 additional authors not shown)

Abstract: Considerable information about the early-stage dynamics of heavy-ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems -- large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynam… ▽ More Considerable information about the early-stage dynamics of heavy-ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems -- large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynamic simulations initialized by a parameterized model for rapidity-dependent energy deposition, which we calibrate on the hadron multiplicity and anisotropic flow coefficients. We utilize Bayesian inference to constrain properties of the early- and late- time dynamics of the system, and highlight the impact of enforcing global energy conservation in our 3D model. △ Less

Submitted 31 December, 2023; originally announced January 2024.

Thermal dilepton production in heavy-ion collisions at beam-energy-scan (BES) energies

Authors: Jessica Churchill, Lipei Du, Bailey Forster, Han Gao, Greg Jackson, Sangyong Jeon, Charles Gale

Abstract: The lepton pair production rate at finite temperature and at next-to-leading-order (NLO) is calculated for quark-gluon plasma at non-zero baryon density. Yields are obtained using a (3+1)D multicomponent simulation capable of reproducing hadronic observables measured in the RHIC Beam Energy Scan. Spectra of intermediate invariant mass dileptons are compared with measurements from the STAR collabor… ▽ More The lepton pair production rate at finite temperature and at next-to-leading-order (NLO) is calculated for quark-gluon plasma at non-zero baryon density. Yields are obtained using a (3+1)D multicomponent simulation capable of reproducing hadronic observables measured in the RHIC Beam Energy Scan. Spectra of intermediate invariant mass dileptons are compared with measurements from the STAR collaboration. The result of a study where the temperature information obtained from the dilepton spectrum from heavy-ion collisions performed at different energies and centralities is presented. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 4 pages, 3 figures. Contribution to the proceedings of Quark Matter 2023 (Houston, TX, 3-9 Sep. 2023)

Exploring the freeze-out hypersurface of relativistic nuclear collisions with a rapidity-dependent thermal model

Authors: Han Gao, Lipei Du, Sangyong Jeon, Charles Gale

Abstract: Considering applications to relativistic heavy-ion collisions, we develop a rapidity-dependent thermal model that includes thermal smearing effect and longitudinal boost. We calibrate the model with thermal yields obtained from a multistage hydrodynamic simulation. Through Bayesian analysis, we find that our model extracts freeze-out thermodynamics with better precision than rapidity-independent m… ▽ More Considering applications to relativistic heavy-ion collisions, we develop a rapidity-dependent thermal model that includes thermal smearing effect and longitudinal boost. We calibrate the model with thermal yields obtained from a multistage hydrodynamic simulation. Through Bayesian analysis, we find that our model extracts freeze-out thermodynamics with better precision than rapidity-independent models. The potential of such a model to constrain longitudinal flow from data is also demonstrated. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: 4 pages, 2 figures. Contribution to the proceedings of Quark Matter 2023 (Houston, TX, 3-9 Sep. 2023)

Constraints on initial baryon stopping and equation of state from directed flow

Authors: Lipei Du, Chun Shen, Sangyong Jeon, Charles Gale

Abstract: Our investigation focuses on the rapidity-dependent directed flow, $v_1(y)$, of identified hadrons in Au+Au collisions across a broad range of $\sqrt{s_{\rm NN}}$ from 7.7 to 200 GeV. Employing a (3+1)-dimensional hybrid framework, our study successfully reproduces the characteristic features of the measured $v_1(y)$ for both mesons and baryons across the considered beam energies. Notably, our ana… ▽ More Our investigation focuses on the rapidity-dependent directed flow, $v_1(y)$, of identified hadrons in Au+Au collisions across a broad range of $\sqrt{s_{\rm NN}}$ from 7.7 to 200 GeV. Employing a (3+1)-dimensional hybrid framework, our study successfully reproduces the characteristic features of the measured $v_1(y)$ for both mesons and baryons across the considered beam energies. Notably, our analysis reveals the constraining power of baryonic $v_1(y)$ on the initial baryon stopping mechanism. Together with mesonic $v_1(y)$, the directed flow serves as a crucial tool for probing the equation of state governing dense nuclear matter at finite chemical potentials. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 4 pages, 1 figure. Contribution to the proceedings of Quark Matter 2023 (Houston, TX, 3-9 Sep. 2023)

Rapidity scan approach for net-baryon cumulants with a statistical thermal model

Authors: Jianing Li, Lipei Du, Shuzhe Shi

Abstract: Utilizing rapidity-dependent measurements to map the QCD phase diagram provides a complementary approach to traditional beam-energy-dependent measurements around midrapidity. The changing nature of thermodynamic properties of QCD matter along the beam axis in heavy-ion collisions at low collision energies both motivates and poses challenges for this method. In this study, we derive the analytical… ▽ More Utilizing rapidity-dependent measurements to map the QCD phase diagram provides a complementary approach to traditional beam-energy-dependent measurements around midrapidity. The changing nature of thermodynamic properties of QCD matter along the beam axis in heavy-ion collisions at low collision energies both motivates and poses challenges for this method. In this study, we derive the analytical cumulant-generating function for subsystems within distinct rapidity windows, while accounting for global net-baryon charge conservation of the full system. Rapidity-dependent net-baryon cumulants are then calculated for a system exhibiting inhomogeneity along the beam axis, and their sensitivity to finite acceptances through changing rapidity bin widths is explored. We highlight the nontrivial behaviors exhibited by these cumulants, underscoring their importance in establishing a noncritical baseline for interpreting net-proton cumulants in the search for the QCD critical point. Finally, we discuss the implications of the rapidity scan for mapping the QCD phase diagram within the current context. △ Less

Submitted 11 March, 2024; v1 submitted 19 November, 2023; originally announced November 2023.

Comments: 16 pages, 9 figures. Content, figures, and references are updated to the publication version

Journal ref: Phys. Rev. C 109, 034906 (2024)

Virtual Photons Shed Light on the Early Temperature of Dense QCD Matter

Authors: Jessica Churchill, Lipei Du, Charles Gale, Greg Jackson, Sangyong Jeon

Abstract: Dileptons produced during heavy-ion collisions represent a unique probe of the QCD phase diagram, and convey information about the state of the strongly interacting system at the moment their preceding off-shell photon is created. In this study, we compute thermal dilepton yields from Au+Au collisions performed at different beam energies, employing a (3+1)-dimensional dynamic framework combined wi… ▽ More Dileptons produced during heavy-ion collisions represent a unique probe of the QCD phase diagram, and convey information about the state of the strongly interacting system at the moment their preceding off-shell photon is created. In this study, we compute thermal dilepton yields from Au+Au collisions performed at different beam energies, employing a (3+1)-dimensional dynamic framework combined with emission rates accurate at next-to-leading order in perturbation theory and which include baryon chemical potential dependencies. By comparing the effective temperature extracted from the thermal dilepton invariant mass spectrum with the average temperature of the fluid, we offer a robust quantitative validation of dileptons as effective probe of the early quark-gluon plasma stage. △ Less

Submitted 18 February, 2024; v1 submitted 12 November, 2023; originally announced November 2023.

Comments: v1: 6+2 pages, 4+2 figures; v2: refs added, a table added in the Supplemental Material. Companion paper: arXiv:2311.06675

Journal ref: Phys. Rev. Lett. 132, 172301 (2024)

Dilepton production at next-to-leading order and intermediate invariant-mass observables

Authors: Jessica Churchill, Lipei Du, Charles Gale, Greg Jackson, Sangyong Jeon

Abstract: The thermal QCD dilepton production rate is calculated at next-to-leading order in the strong coupling and at finite baryon chemical potential. The two-loop virtual photon self-energy is evaluated using finite temperature field theory and combined consistently with the self-energy in the Landau-Pomeranchuk-Migdal regime. We present new results for a dense baryonic plasma. The rates are then integr… ▽ More The thermal QCD dilepton production rate is calculated at next-to-leading order in the strong coupling and at finite baryon chemical potential. The two-loop virtual photon self-energy is evaluated using finite temperature field theory and combined consistently with the self-energy in the Landau-Pomeranchuk-Migdal regime. We present new results for a dense baryonic plasma. The rates are then integrated using (3+1)-dimensional fluid-dynamical simulations calibrated to reproduce hadronic experimental results obtained at RHIC at energies ranging from those of the Beam Energy Scan to $\sqrt{s_{_{\rm NN}}} = 200$ GeV. We elaborate on the ability for dileptons to relay information about the plasma baryonic content and temperature. △ Less

Submitted 6 May, 2024; v1 submitted 11 November, 2023; originally announced November 2023.

Comments: v1: 19 pages, 15 figures; v2: minor revision; v3: published version. Companion paper: Phys. Rev. Lett. 132, 172301 (2024) [arXiv:2311.06951]

Journal ref: Phys. Rev. C 109, 044915 (2024)

A new metric improving Bayesian calibration of a multistage approach studying hadron and inclusive jet suppression

Authors: W. Fan, G. Vujanovic, S. A. Bass, A. Angerami, R. Arora, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kordell II, A. Kumar, J. Latessa, Y. -J. Lee , et al. (30 additional authors not shown)

Abstract: We study parton energy-momentum exchange with the quark gluon plasma (QGP) within a multistage approach composed of in-medium DGLAP evolution at high virtuality, and (linearized) Boltzmann Transport formalism at lower virtuality. This multistage simulation is then calibrated in comparison with high $p_T$ charged hadrons, D-mesons, and the inclusive jet nuclear modification factors, using Bayesian… ▽ More We study parton energy-momentum exchange with the quark gluon plasma (QGP) within a multistage approach composed of in-medium DGLAP evolution at high virtuality, and (linearized) Boltzmann Transport formalism at lower virtuality. This multistage simulation is then calibrated in comparison with high $p_T$ charged hadrons, D-mesons, and the inclusive jet nuclear modification factors, using Bayesian model-to-data comparison, to extract the virtuality-dependent transverse momentum broadening transport coefficient $\hat{q}$. To facilitate this undertaking, we develop a quantitative metric for validating the Bayesian workflow, which is used to analyze the sensitivity of various model parameters to individual observables. The usefulness of this new metric in improving Bayesian model emulation is shown to be highly beneficial for future such analyses. △ Less

Submitted 27 October, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 24 pages, 16 figures

Multiscale evolution of heavy flavor in the QGP

Authors: G. Vujanovic, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, M. Kordell II, A. Kumar, J. Latessa, Y. -J. Lee , et al. (30 additional authors not shown)

Abstract: Shower development dynamics for a jet traveling through the quark-gluon plasma (QGP) is a multiscale process, where the heavy flavor mass is an important scale. During the high virtuality portion of the jet evolution in the QGP, emission of gluons from a heavy flavor is modified owing to heavy quark mass. Medium-induced radiation of heavy flavor is sensitive to microscopic processes (e.g. diffusio… ▽ More Shower development dynamics for a jet traveling through the quark-gluon plasma (QGP) is a multiscale process, where the heavy flavor mass is an important scale. During the high virtuality portion of the jet evolution in the QGP, emission of gluons from a heavy flavor is modified owing to heavy quark mass. Medium-induced radiation of heavy flavor is sensitive to microscopic processes (e.g. diffusion), whose virtuality dependence is phenomenologically explored in this study. In the lower virtuality part of shower evolution, i.e. when the mass is comparable to the virtuality of the parton, scattering and radiation processes of heavy quarks differ from light quarks. The effects of these mechanisms on shower development in heavy flavor tagged showers in the QGP is explored here. Furthermore, this multiscale study examines dynamical pair production of heavy flavor (via virtual gluon splittings) and their subsequent evolution in the QGP, which is not possible otherwise. A realistic event-by-event simulation is performed using the JETSCAPE framework. Energy-momentum exchange with the medium proceeds using a weak coupling recoil approach. Using leading hadron and open heavy flavor observables, differences in heavy versus light quark energy-loss mechanisms are explored, while the importance of heavy flavor pair production is highlighted along with future directions to study. △ Less

Submitted 27 October, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 6 pages, 3 figures, contribution to the Hard Probes 2023 proceedings

Effects of multi-scale jet-medium interactions on jet substructures

Authors: JETSCAPE Collaboration, Y. Tachibana, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, K. Kauder, L. Kasper, W. Ke, M. Kelsey , et al. (35 additional authors not shown)

Abstract: We utilize event-by-event Monte Carlo simulations within the JETSCAPE framework to examine scale-dependent jet-medium interactions in heavy-ion collisions. The reduction in jet-medium interaction during the early high-virtuality stage, where the medium is resolved at a short distance scale, is emphasized as a key element in explaining multiple jet observables, particularly substructures, simultane… ▽ More We utilize event-by-event Monte Carlo simulations within the JETSCAPE framework to examine scale-dependent jet-medium interactions in heavy-ion collisions. The reduction in jet-medium interaction during the early high-virtuality stage, where the medium is resolved at a short distance scale, is emphasized as a key element in explaining multiple jet observables, particularly substructures, simultaneously. By employing the MATTER+LBT setup, which incorporates this explicit reduction of medium effects at high virtuality, we investigate jet substructure observables, such as Soft Drop groomed observables. When contrasted with existing data, our findings spotlight the significant influence of the reduction at the early high-virtuality stages. Furthermore, we study the substructure of gamma-tagged jets, providing predictive insights for future experimental analyses. This broadens our understanding of the various contributing factors involved in modifying jet substructures. △ Less

Submitted 16 July, 2023; originally announced July 2023.

Comments: 6 pages, 4 figures, contribution to the Hard Probes 2023 proceedings

Hot QCD White Paper

Authors: M. Arslandok, S. A. Bass, A. A. Baty, I. Bautista, C. Beattie, F. Becattini, R. Bellwied, Y. Berdnikov, A. Berdnikov, J. Bielcik, J. T. Blair, F. Bock, B. Boimska, H. Bossi, H. Caines, Y. Chen, Y. -T. Chien, M. Chiu, M. E. Connors, M. Csanád, C. L. da Silva, A. P. Dash, G. David, K. Dehmelt, V. Dexheimer , et al. (149 additional authors not shown)

Abstract: Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the… ▽ More Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the temperature dependence of the transport properties of quark-gluon plasma, the phase diagram of nuclear matter, the interaction of quarks and gluons at different scales and much more. This document, as part of the 2023 nuclear science long range planning process, was written to review the progress in hot QCD since the 2015 Long Range Plan for Nuclear Science, as well as highlight the realization of previous recommendations, and present opportunities for the next decade, building on the accomplishments and investments made in theoretical developments and the construction of new detectors. Furthermore, this document provides additional context to support the recommendations voted on at the Joint Hot and Cold QCD Town Hall Meeting, which are reported in a separate document. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 190 pages, 69 figures

Rapidity scan with multistage hydrodynamic and statistical thermal models

Authors: Lipei Du, Han Gao, Sangyong Jeon, Charles Gale

Abstract: We calibrate a (3+1)-dimensional multistage hybrid framework using the measured pseudo-rapidity distribution of charged particles and rapidity distribution of net protons for central Au+Au collisions at $\sqrt{s_{\rm NN}}=7.7,\,19.6,\,62.4,\,200$ GeV. We then study the thermodynamic properties of the nuclear matter along the beam direction, and the phase diagram regions probed by the hadronization… ▽ More We calibrate a (3+1)-dimensional multistage hybrid framework using the measured pseudo-rapidity distribution of charged particles and rapidity distribution of net protons for central Au+Au collisions at $\sqrt{s_{\rm NN}}=7.7,\,19.6,\,62.4,\,200$ GeV. We then study the thermodynamic properties of the nuclear matter along the beam direction, and the phase diagram regions probed by the hadronization process near the chemical freeze-out. Using the rapidity-dependent thermal yields of identified particles with full rapidity coverage from the hybrid framework, we apply different scenarios of the statistical thermal model to extract the thermodynamic parameters at the freeze-out, with the known system properties from the hybrid model as a closure test. We find significant theoretical uncertainties in the thermal models when applied to regions away from midrapidity. We also propose a thermal model inspired by the hybrid approach that includes thermal smearing and longitudinal flow for the nuclear matter created at low beam energies. △ Less

Submitted 22 January, 2024; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: v1: 17 pages, 12 figures; v2: structure mainly adjusted in Sec. IIB, results remain unchanged

Journal ref: Phys. Rev. C 109, 014907 (2024)

Dense Nuclear Matter Equation of State from Heavy-Ion Collisions

Authors: Agnieszka Sorensen, Kshitij Agarwal, Kyle W. Brown, Zbigniew Chajęcki, Paweł Danielewicz, Christian Drischler, Stefano Gandolfi, Jeremy W. Holt, Matthias Kaminski, Che-Ming Ko, Rohit Kumar, Bao-An Li, William G. Lynch, Alan B. McIntosh, William G. Newton, Scott Pratt, Oleh Savchuk, Maria Stefaniak, Ingo Tews, ManYee Betty Tsang, Ramona Vogt, Hermann Wolter, Hanna Zbroszczyk, Navid Abbasi, Jörg Aichelin , et al. (111 additional authors not shown)

Abstract: 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 mu… ▽ More 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. 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. △ Less

Submitted 25 January, 2024; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

Report number: INT-PUB-23-001, LA-UR-23-20514, LLNL-TR-844629

Journal ref: Prog.Part.Nucl.Phys. 134 (2024) 104080

Hard Jet Substructure in a Multi-stage Approach

Authors: Y. Tachibana, A. Kumar, A. Majumder, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, K. Kauder, L. Kasper, W. Ke , et al. (34 additional authors not shown)

Abstract: We present predictions and postdictions for a wide variety of hard jet-substructure observables using a multi-stage model within the JETSCAPE framework. The details of the multi-stage model and the various parameter choices are described in [A. Kumar et al., arXiv:2204.01163]. A novel feature of this model is the presence of two stages of jet modification: a high virtuality phase (modeled using MA… ▽ More We present predictions and postdictions for a wide variety of hard jet-substructure observables using a multi-stage model within the JETSCAPE framework. The details of the multi-stage model and the various parameter choices are described in [A. Kumar et al., arXiv:2204.01163]. A novel feature of this model is the presence of two stages of jet modification: a high virtuality phase (modeled using MATTER), where coherence effects diminish medium-induced radiation, and a lower virtuality phase (modeled using LBT), where parton splits are fully resolved by the medium as they endure multiple scattering induced energy loss. Energy loss calculations are carried out on event-by-event viscous fluid dynamic backgrounds constrained by experimental data. The uniformed and consistent descriptions of multiple experimental observables demonstrate the essential role of coherence effects and the multi-stage modeling of the jet evolution. Using the best choice of parameters from [A. Kumar et al., arXiv:2204.01163], and with no further tuning, we present calculations for the medium modified jet fragmentation function, the groomed jet momentum fraction $z_g$ and angular separation $r_g$ distributions, as well as the nuclear modification factor of groomed jets. These calculations provide accurate descriptions of published and preliminary data from experiments at RHIC and LHC. Furthermore, we provide predictions from the multi-stage model for future measurements at RHIC. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 20 pages, 12 figures

Comprehensive Study of Multi-scale Jet-medium Interaction

Authors: Y. Tachibana, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, W. Ke, M. Kelsey, M. Kordell II, A. Kumar , et al. (33 additional authors not shown)

Abstract: We explore jet-medium interactions at various scales in high-energy heavy-ion collisions using the JETSCAPE framework. The physics of the multi-stage modeling and the coherence effect at high virtuality is discussed through the results of multiple jet and high-$p_{\mathrm{T}}$ particle observables, compared with experimental data. Furthermore, we investigate the jet-medium interaction involved in… ▽ More We explore jet-medium interactions at various scales in high-energy heavy-ion collisions using the JETSCAPE framework. The physics of the multi-stage modeling and the coherence effect at high virtuality is discussed through the results of multiple jet and high-$p_{\mathrm{T}}$ particle observables, compared with experimental data. Furthermore, we investigate the jet-medium interaction involved in the hadronization process. △ Less

Submitted 23 December, 2022; originally announced December 2022.

Comments: 6 pages, 5 figures, contribution to the Quark Matter 2022 proceedings

Journal ref: Acta Phys. Pol. B Proc. Suppl. 16, 1-A50 (2023)

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

Authors: Lipei Du, Chun Shen, Sangyong Jeon, Charles Gale

Abstract: 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{s_{\rm NN}}$ ranging from 7.7 to 200 GeV. The dynamics in the beam direction is constrained using the… ▽ More 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{s_{\rm NN}}$ 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. △ Less

Submitted 11 November, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: v1: 6+3 pages, 4 figures; v2: add discussion on the Lund string model

Journal ref: Phys. Rev. C 108, L041901 (2023)

Bayesian analysis of QGP jet transport using multi-scale modeling applied to inclusive hadron and reconstructed jet data

Authors: R. Ehlers, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, L. Du, T. Dai, H. Elfner, W. Fan, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, L. Kasper, W. Ke, M. Kelsey, M. Kordell II, A. Kumar, J. Latessa , et al. (33 additional authors not shown)

Abstract: The JETSCAPE Collaboration reports a new determination of jet transport coefficients in the Quark-Gluon Plasma, using both reconstructed jet and hadron data measured at RHIC and the LHC. The JETSCAPE framework incorporates detailed modeling of the dynamical evolution of the QGP; a multi-stage theoretical approach to in-medium jet evolution and medium response; and Bayesian inference for quantitati… ▽ More The JETSCAPE Collaboration reports a new determination of jet transport coefficients in the Quark-Gluon Plasma, using both reconstructed jet and hadron data measured at RHIC and the LHC. The JETSCAPE framework incorporates detailed modeling of the dynamical evolution of the QGP; a multi-stage theoretical approach to in-medium jet evolution and medium response; and Bayesian inference for quantitative comparison of model calculations and data. The multi-stage framework incorporates multiple models to cover a broad range in scale of the in-medium parton shower evolution, with dynamical choice of model that depends on the current virtuality or energy of the parton. We will discuss the physics of the multi-stage modeling, and then present a new Bayesian analysis incorporating it. This analysis extends the recently published JETSCAPE determination of the jet transport parameter $\hat{q}$ that was based solely on inclusive hadron suppression data, by incorporating reconstructed jet measurements of quenching. We explore the functional dependence of jet transport coefficients on QGP temperature and jet energy and virtuality, and report the consistency and tensions found for current jet quenching modeling with hadron and reconstructed jet data over a wide range in kinematics and $\sqrt{s_{\text{NN}}}$. This analysis represents the next step in the program of comprehensive analysis of jet quenching phenomenology and its constraint of properties of the QGP. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 6 pages, 2 figures, contribution to the Quark Matter 2022 proceedings

Multi-scale evolution of charmed particles in a nuclear medium

Authors: JETSCAPE collaboration, W. Fan, G. Vujanovic, S. A. Bass, A. Majumder, A. Angerami, R. Arora, S. Cao, Y. Chen, T. Dai, L. Du, R. Ehlers, H. Elfner, R. J. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, Y. Ji, K. Kauder, L. Kasper, W. Ke , et al. (35 additional authors not shown)

Abstract: Parton energy-momentum exchange with the quark gluon plasma (QGP) is a multi-scale problem. In this work, we calculate the interaction of charm quarks with the QGP within the higher twist formalism at high virtuality and high energy using the MATTER model, while the low virtuality and high energy portion is treated via a (linearized) Boltzmann Transport (LBT) formalism. Coherence effect that reduc… ▽ More Parton energy-momentum exchange with the quark gluon plasma (QGP) is a multi-scale problem. In this work, we calculate the interaction of charm quarks with the QGP within the higher twist formalism at high virtuality and high energy using the MATTER model, while the low virtuality and high energy portion is treated via a (linearized) Boltzmann Transport (LBT) formalism. Coherence effect that reduces the medium-induced emission rate in the MATTER model is also taken into account. The interplay between these two formalisms is studied in detail and used to produce a good description of the D-meson and charged hadron nuclear modification factor RAA across multiple centralities. All calculations were carried out utilizing the JETSCAPE framework. △ Less

Submitted 13 May, 2023; v1 submitted 1 August, 2022; originally announced August 2022.

Journal ref: 107(5), 2023, 054901

Far-from-equilibrium attractor in non-conformal plasmas

Authors: Sunil Jaiswal, Subrata Pal, Chandrodoy Chattopadhyay, Lipei Du, Ulrich Heinz

Abstract: We study the far-off-equilibrium dynamics of a Bjorken expanding non-conformal system within kinetic theory and hydrodynamics. We show that, in contrast to the conformal case, neither shear nor bulk viscous pressure relax quickly to a non-equilibrium attractor. In kinetic theory an early-time, far-from-equilibrium attractor exists for the scaled longitudinal pressure, driven by the rapid longitudi… ▽ More We study the far-off-equilibrium dynamics of a Bjorken expanding non-conformal system within kinetic theory and hydrodynamics. We show that, in contrast to the conformal case, neither shear nor bulk viscous pressure relax quickly to a non-equilibrium attractor. In kinetic theory an early-time, far-from-equilibrium attractor exists for the scaled longitudinal pressure, driven by the rapid longitudinal expansion of the medium. Second-order dissipative hydrodynamics fails to accurately describe this attractor, but a modified anisotropic hydrodynamic formulation reproduces it and provides excellent agreement with kinetic theory. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: Contribution to the proceedings of Quark Matter 2022. 6 pages, 3 figures

Response of a baryon-charged medium to energetic partons

Authors: Lipei Du, Ulrich Heinz

Abstract: We explore the response to energetic partons of a baryon-charged medium produced in low energy heavy-ion collisions in which the partonic energy loss rate is expected to depend on both temperature and baryon chemical potential. The energy and momentum deposited by the partons are described by dynamical sources added to hydrodynamic equations of motion. We study the distortions of various hydrodyna… ▽ More We explore the response to energetic partons of a baryon-charged medium produced in low energy heavy-ion collisions in which the partonic energy loss rate is expected to depend on both temperature and baryon chemical potential. The energy and momentum deposited by the partons are described by dynamical sources added to hydrodynamic equations of motion. We study the distortions of various hydrodynamic quantities, especially the energy and baryon densities, induced by energetic partons plowing through the medium. By studying the contribution from this medium response to the emission spectra of several identified hadron species we identify qualitative differences between the jet-induced modifications for mesons and net protons. △ Less

Submitted 11 September, 2022; v1 submitted 20 May, 2022; originally announced May 2022.

Comments: v3: title changed, published on PRC

Journal ref: Phys. Rev. C 106, 034903 (2022)

Inclusive jet and hadron suppression in a multistage approach

Authors: A. Kumar, Y. Tachibana, C. Sirimanna, G. Vujanovic, S. Cao, A. Majumder, Y. Chen, L. Du, R. Ehlers, D. Everett, W. Fan, Y. He, J. Mulligan, C. Park, A. Angerami, R. Arora, S. A. Bass, T. Dai, H. Elfner, R. J. Fries, C. Gale, F. Garza, M. Heffernan, U. Heinz, B. V. Jacak , et al. (35 additional authors not shown)

Abstract: We present a new study of jet interactions in the quark-gluon plasma created in high-energy heavy-ion collisions, using a multistage event generator within the JETSCAPE framework. We focus on medium-induced modifications in the rate of inclusive jets and high transverse momentum (high-$p_{\mathrm{T}}$) hadrons. Scattering-induced jet energy loss is calculated in two stages: A high virtuality stage… ▽ More We present a new study of jet interactions in the quark-gluon plasma created in high-energy heavy-ion collisions, using a multistage event generator within the JETSCAPE framework. We focus on medium-induced modifications in the rate of inclusive jets and high transverse momentum (high-$p_{\mathrm{T}}$) hadrons. Scattering-induced jet energy loss is calculated in two stages: A high virtuality stage based on the MATTER model, in which scattering of highly virtual partons modifies the vacuum radiation pattern, and a second stage at lower jet virtuality based on the LBT model, in which leading partons gain and lose virtuality by scattering and radiation. Coherence effects that reduce the medium-induced emission rate in the MATTER phase are also included. The TRENTo model is used for initial conditions, and the (2+1)dimensional VISHNU model is used for viscous hydrodynamic evolution. Jet interactions with the medium are modeled via 2-to-2 scattering with Debye screened potentials, in which the recoiling partons are tracked, hadronized, and included in the jet clustering. Holes left in the medium are also tracked and subtracted to conserve transverse momentum. Calculations of the nuclear modification factor ($R_{\mathrm{AA}}$) for inclusive jets and high-$p_{\mathrm{T}}$ hadrons are compared to experimental measurements at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC). Within this framework, we find that with one extra parameter which codifies the transition between stages of jet modification -- along with the typical parameters such as the coupling in the medium, the start and stop criteria etc. -- we can describe these data at all energies for central and semicentral collisions without a rescaling of the jet transport coefficient $\hat{q}$. △ Less

Submitted 16 April, 2023; v1 submitted 3 April, 2022; originally announced April 2022.

Comments: 34 pages, 23 figures

Journal ref: Phys. Rev. C 107, no.3, 034911 (2023)

Role of bulk viscosity in deuteron production in ultrarelativistic nuclear collisions

Authors: D. Everett, D. Oliinychenko, M. Luzum, J. -F. Paquet, G. Vujanovic, S. A. Bass, L. Du, C. Gale, M. Heffernan, U. Heinz, L. Kasper, W. Ke, D. Liyanage, A. Majumder, A. Mankolli, C. Shen, D. Soeder, J. Velkovska, A. Angerami, R. Arora, S. Cao, Y. Chen, T. Dai, R. Ehlers, H. Elfner , et al. (31 additional authors not shown)

Abstract: We use a Bayesian-calibrated multistage viscous hydrodynamic model to explore deuteron yield, mean transverse momentum and flow observables in LHC Pb-Pb collisions. We explore theoretical uncertainty in the production of deuterons, including (i) the contribution of thermal deuterons, (ii) models for the subsequent formation of deuterons (hadronic transport vs coalescence) and (iii) the overall sen… ▽ More We use a Bayesian-calibrated multistage viscous hydrodynamic model to explore deuteron yield, mean transverse momentum and flow observables in LHC Pb-Pb collisions. We explore theoretical uncertainty in the production of deuterons, including (i) the contribution of thermal deuterons, (ii) models for the subsequent formation of deuterons (hadronic transport vs coalescence) and (iii) the overall sensitivity of the results to the hydrodynamic model -- in particular to bulk viscosity, which is often neglected in studies of deuteron production. Using physical parameters set by a comparison to only light hadron observables, we find good agreement with measurements of the mean transverse momentum $\langle p_T \rangle$ and elliptic flow $v_2$ of deuterons; however, tension is observed with experimental data for the deuteron multiplicity in central collisions. The results are found to be sensitive to each of the mentioned theoretical uncertainties, with a particular sensitivity to bulk viscosity, indicating that the latter is an important ingredient for an accurate treatment of deuteron production. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: 11 pages, 8 figures

Baryon diffusion near the QCD critical point

Authors: Lipei Du, Xin An, Ulrich Heinz

Abstract: Fireballs created in relativistic heavy-ion collisions at different beam energies have been argued to follow different trajectories in the QCD phase diagram in which the QCD critical point serves as a landmark. Using a (1+1)-dimensional model setting with transverse homogeneity, we study the complexities introduced by the fact that the evolution history of each fireball cannot be characterized by… ▽ More Fireballs created in relativistic heavy-ion collisions at different beam energies have been argued to follow different trajectories in the QCD phase diagram in which the QCD critical point serves as a landmark. Using a (1+1)-dimensional model setting with transverse homogeneity, we study the complexities introduced by the fact that the evolution history of each fireball cannot be characterized by a single trajectory but rather covers an entire swath of the phase diagram, with the finally emitted hadron spectra integrating over contributions from many different trajectories. Studying the phase diagram trajectories of fluid cells at different space-time rapidities, we explore how baryon diffusion shuffles them around, and how they are affected by critical dynamics near the QCD critical point. We find a striking insensitivity of baryon diffusion to critical effects. Its origins are analyzed and possible implications discussed. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 6 pages, 3 figures. New discussion on the viscous effects (shear and bulk viscosity) on phase diagram trajectories. Proceedings of CPOD 2021. Contributed talk on arXiv:2107.02302. arXiv admin note: text overlap with arXiv:2107.02302

The BEST framework for the search for the QCD critical point and the chiral magnetic effect

Authors: Xin An, Marcus Bluhm, Lipei Du, Gerald V. Dunne, Hannah Elfner, Charles Gale, Joaquin Grefa, Ulrich Heinz, Anping Huang, Jamie M. Karthein, Dmitri E. Kharzeev, Volker Koch, Jinfeng Liao, Shiyong Li, Mauricio Martinez, Michael McNelis, Debora Mroczek, Swagato Mukherjee, Marlene Nahrgang, Angel R. Nava Acuna, Jacquelyn Noronha-Hostler, Dmytro Oliinychenko, Paolo Parotto, Israel Portillo, Maneesha Sushama Pradeep , et al. (18 additional authors not shown)

Abstract: The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. W… ▽ More The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Finally, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world. △ Less

Submitted 22 November, 2021; v1 submitted 31 August, 2021; originally announced August 2021.

Comments: 103 pages

Journal ref: Nuclear Physics A, Volume 1017, January 2022, 122343

On non-conformal kinetic theory and hydrodynamics for Bjorken flow

Authors: Sunil Jaiswal, Chandrodoy Chattopadhyay, Lipei Du, Ulrich Heinz, Subrata Pal

Abstract: Using and comparing kinetic theory and second-order Chapman-Enskog hydrodynamics, we study the non-conformal dynamics of a system undergoing Bjorken expansion. We use the concept of `free-streaming fixed lines' for scaled shear and bulk stresses in non-conformal kinetic theory and hydrodynamics, and show that these `fixed lines' behave as early-time attractors and repellors of the evolution. In th… ▽ More Using and comparing kinetic theory and second-order Chapman-Enskog hydrodynamics, we study the non-conformal dynamics of a system undergoing Bjorken expansion. We use the concept of `free-streaming fixed lines' for scaled shear and bulk stresses in non-conformal kinetic theory and hydrodynamics, and show that these `fixed lines' behave as early-time attractors and repellors of the evolution. In the conformal limit, the free-streaming fixed lines reduce to the well-known fixed points of conformal Bjorken dynamics. A new fixed point in the free streaming regime is identified which lies at the intersection of these fixed lines. Contrary to the conformal scenario, both kinetic theory and hydrodynamics predict the absence of attractor behavior in the normalised shear stress channel. In kinetic theory a far-off-equilibrium attractor is found for the normalised effective longitudinal pressure, driven by rapid longitudinal expansion. Second-order viscous hydrodynamics fails to accurately describe this attractor. From a thorough analysis of the free-streaming dynamics in Chapman-Enskog hydrodynamics we conclude that this failure results from an inaccurate approximation of the fixed lines and a related incorrect description of the nature of the fixed point. A modified anisotropic hydrodynamic description is presented that provides excellent agreement with kinetic theory results and reproduces the far-from-equilibrium attractor for the scaled longitudinal pressure. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Comments: 28 pages, 22 figures

Hydrodynamic description of the baryon-charged quark-gluon plasma

Authors: Lipei Du

Abstract: One of the primary goals of nuclear physics is studying the phase diagram of Quantum Chromodynamics, where a hypothetical critical point serves as a landmark. A systematic model-data comparison of heavy-ion collisions at center-of-mass energies between 1 and 100 GeV per nucleon is essential for locating the critical point and the phase boundary between the deconfined quark-gluon plasma and the con… ▽ More One of the primary goals of nuclear physics is studying the phase diagram of Quantum Chromodynamics, where a hypothetical critical point serves as a landmark. A systematic model-data comparison of heavy-ion collisions at center-of-mass energies between 1 and 100 GeV per nucleon is essential for locating the critical point and the phase boundary between the deconfined quark-gluon plasma and the confined hadron resonance gas. At these energies the net baryon density of the system can be high and critical fluctuations can become essential in the presence of the critical point. Simulating their dynamical evolution thus becomes an indispensable part of theoretical modeling. In this thesis we first present the (3+1)-dimensional relativistic hydrodynamic code BEShydro, which solves the equations of motion of second-order Denicol-Niemi-Molnar-Rischke theory, including bulk and shear viscous components as well as baryon diffusion current. We then study the effects caused by the baryon diffusion on the longitudinal dynamics and on the phase diagram trajectories of fluid cells at different space-time rapidities of the system, and how they are affected by critical dynamics near the critical point. We finally explore the evolution of non-hydrodynamic slow processes describing long wavelength critical fluctuations near the critical point, by extending the conventional hydrodynamic description by coupling it to additional explicitly evolving slow modes, and their back-reaction to the bulk matter properties. △ Less

Submitted 18 July, 2021; originally announced July 2021.

Comments: PhD thesis, 203 pages

Non-conformal attractor in boost-invariant plasmas

Authors: Chandrodoy Chattopadhyay, Sunil Jaiswal, Lipei Du, Ulrich Heinz, Subrata Pal

Abstract: We study the dissipative evolution of (0+1)-dimensionally expanding media with Bjorken symmetry using the Boltzmann equation for massive particles in relaxation-time approximation. Breaking conformal symmetry by a mass induces a non-zero bulk viscous pressure in the medium. It is shown that even a small mass (in units of the local temperature) drastically modifies the well-known attractor for the… ▽ More We study the dissipative evolution of (0+1)-dimensionally expanding media with Bjorken symmetry using the Boltzmann equation for massive particles in relaxation-time approximation. Breaking conformal symmetry by a mass induces a non-zero bulk viscous pressure in the medium. It is shown that even a small mass (in units of the local temperature) drastically modifies the well-known attractor for the shear Reynolds number previously observed in massless systems. For generic nonzero particle mass, neither the shear nor the bulk viscous pressure relax quickly to a non-equilibrium attractor; they approach the hydrodynamic limit only late, at small values of the inverse Reynolds numbers. Only the longitudinal pressure, which is a combination of thermal, shear and bulk viscous pressures, continues to show early approach to a far-off-equilibrium attractor, driven by the rapid longitudinal expansion at early times. Second-order dissipative hydrodynamics based on a gradient expansion around locally isotropic thermal equilibrium fails to reproduce this attractor. △ Less

Submitted 8 July, 2021; originally announced July 2021.

Comments: 7 pages, 4 figures

Baryon transport and the QCD critical point

Authors: Lipei Du, Xin An, Ulrich Heinz

Abstract: Fireballs created in relativistic heavy-ion collisions at different beam energies have been argued to follow different trajectories in the QCD phase diagram in which the QCD critical point serves as a landmark. Using a (1+1)-dimensional model setting with transverse homogeneity, we study the complexities introduced by the fact that the evolution history of each fireball cannot be characterized by… ▽ More Fireballs created in relativistic heavy-ion collisions at different beam energies have been argued to follow different trajectories in the QCD phase diagram in which the QCD critical point serves as a landmark. Using a (1+1)-dimensional model setting with transverse homogeneity, we study the complexities introduced by the fact that the evolution history of each fireball cannot be characterized by a single trajectory but rather covers an entire swath of the phase diagram, with the finally emitted hadron spectra integrating over contributions from many different trajectories. Studying the phase diagram trajectories of fluid cells at different space-time rapidities, we explore how baryon diffusion shuffles them around, and how they are affected by critical dynamics near the QCD critical point. We find a striking insensitivity of baryon diffusion to critical effects. Its origins are analyzed and possible implications discussed. △ Less

Submitted 16 January, 2022; v1 submitted 5 July, 2021; originally announced July 2021.

Comments: 26 pages, 17 figures; v2: Fig.14 moved to Fig.1, Sec.IIIC included more clarifications; published on PRC

Determining the jet transport coefficient $\hat{q}$ of the quark-gluon plasma using Bayesian parameter estimation

Authors: J. Mulligan, A. Angerami, R. Arora, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, F. Garza, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, W. Ke, B. Kim , et al. (24 additional authors not shown)

Abstract: We present a new determination of $\hat{q}$, the jet transport coefficient of the quark-gluon plasma. Using the JETSCAPE framework, we use Bayesian parameter estimation to constrain the dependence of $\hat{q}$ on the jet energy, virtuality, and medium temperature from experimental measurements of inclusive hadron suppression in Au-Au collisions at RHIC and Pb-Pb collisions at the LHC. These result… ▽ More We present a new determination of $\hat{q}$, the jet transport coefficient of the quark-gluon plasma. Using the JETSCAPE framework, we use Bayesian parameter estimation to constrain the dependence of $\hat{q}$ on the jet energy, virtuality, and medium temperature from experimental measurements of inclusive hadron suppression in Au-Au collisions at RHIC and Pb-Pb collisions at the LHC. These results are based on a multi-stage theoretical approach to in-medium jet evolution with the MATTER and LBT jet quenching models. The functional dependence of $\hat{q}$ on jet energy, 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. The correlation of experimental systematic uncertainties is accounted for in the parameter extraction. These results provide state-of-the-art constraints on $\hat{q}$ and lay the groundwork to extract additional properties of the quark-gluon plasma from jet measurements in heavy-ion collisions. △ Less

Submitted 21 June, 2021; originally announced June 2021.

Comments: contribution to the 2021 QCD session of the 55th Recontres de Moriond

Determining the jet transport coefficient $\hat{q}$ from inclusive hadron suppression measurements using Bayesian parameter estimation

Authors: S. Cao, Y. Chen, J. Coleman, J. Mulligan, P. M. Jacobs, R. A. Soltz, A. Angerami, R. Arora, S. A. Bass, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. J. Fries, C. Gale, F. Garza, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, S. Jeon, W. Ke , et al. (22 additional authors not shown)

Abstract: 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 multi-stage 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 su… ▽ More 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 multi-stage 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 RHIC and Pb+Pb collisions at the 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 multi-stage 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_\mathrm{T}^\mathrm{jet}>40$ GeV/c. The value of $Q_0$, determined here for the first time, is in the range 2.0-2.7 GeV. △ Less

Submitted 28 July, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: Published in Phys Rev C

Journal ref: Phys. Rev. C 104, 024905 (2021)

Multi-system Bayesian constraints on the transport coefficients of QCD matter

Authors: D. Everett, W. Ke, J. -F. Paquet, G. Vujanovic, S. A. Bass, L. Du, C. Gale, M. Heffernan, U. Heinz, D. Liyanage, M. Luzum, A. Majumder, M. McNelis, C. Shen, Y. Xu, A. Angerami, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, R. Ehlers, H. Elfner, W. Fan, R. J. Fries , et al. (23 additional authors not shown)

Abstract: We study the properties of the strongly-coupled quark-gluon plasma with a multistage model of heavy ion collisions that combines the T$_\mathrm{R}$ENTo initial condition ansatz, free-streaming, viscous relativistic hydrodynamics, and a relativistic hadronic transport. A model-to-data comparison with Bayesian inference is performed, revisiting assumptions made in previous studies. The role of param… ▽ More We study the properties of the strongly-coupled quark-gluon plasma with a multistage model of heavy ion collisions that combines the T$_\mathrm{R}$ENTo initial condition ansatz, free-streaming, viscous relativistic hydrodynamics, and a relativistic hadronic transport. A model-to-data comparison with Bayesian inference is performed, revisiting assumptions made in previous studies. The role of parameter priors is studied in light of their importance towards the interpretation of results. We emphasize the use of closure tests to perform extensive validation of the analysis workflow before comparison with observations. Our study combines measurements from the Large Hadron Collider and the Relativistic Heavy Ion Collider, achieving a good simultaneous description of a wide range of hadronic observables from both colliders. The selected experimental data provide reasonable constraints on the shear and the bulk viscosities of the quark-gluon plasma at $T\sim$ 150-250 MeV, but their constraining power degrades at higher temperatures $T \gtrsim 250$ MeV. Furthermore, these viscosity constraints are found to depend significantly on how viscous corrections are handled in the transition from hydrodynamics to the hadronic transport. Several other model parameters, including the free-streaming time, show similar model sensitivity while the initial condition parameters associated with the T$_\mathrm{R}$ENTo ansatz are quite robust against variations of the particlization prescription. We also report on the sensitivity of individual observables to the various model parameters. Finally, Bayesian model selection is used to quantitatively compare the agreement with measurements for different sets of model assumptions, including different particlization models and different choices for which parameters are allowed to vary between RHIC and LHC energies. △ Less

Submitted 6 November, 2020; v1 submitted 2 November, 2020; originally announced November 2020.

Comments: 51 pages, including 35 figures and 8 appendices, long companion paper to arXiv:2010.03928. A useful visualization tool to see the effect of varying individual model parameters on physical observables can be found at jetscape.org/sims-widget. Some references and discussion added. This version submitted for publication

Journal ref: Phys. Rev. C 103, 054904 (2021)

Phenomenological constraints on the transport properties of QCD matter with data-driven model averaging

Authors: D. Everett, W. Ke, J. -F. Paquet, G. Vujanovic, S. A. Bass, L. Du, C. Gale, M. Heffernan, U. Heinz, D. Liyanage, M. Luzum, A. Majumder, M. McNelis, C. Shen, Y. Xu, A. Angerami, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, R. Ehlers, H. Elfner, W. Fan, R. J. Fries , et al. (23 additional authors not shown)

Abstract: Using combined data from the Relativistic Heavy Ion and Large Hadron Colliders, we constrain the shear and bulk viscosities of quark-gluon plasma (QGP) at temperatures of ${\sim\,}150{-}350$ MeV. We use Bayesian inference to translate experimental and theoretical uncertainties into probabilistic constraints for the viscosities. With Bayesian Model Averaging we account for the irreducible model amb… ▽ More Using combined data from the Relativistic Heavy Ion and Large Hadron Colliders, we constrain the shear and bulk viscosities of quark-gluon plasma (QGP) at temperatures of ${\sim\,}150{-}350$ MeV. We use Bayesian inference to translate experimental and theoretical uncertainties into probabilistic constraints for the viscosities. With Bayesian Model Averaging we account for the irreducible model ambiguities in the transition from a fluid description of the QGP to hadronic transport in the final evolution stage, providing the most reliable phenomenological constraints to date on the QGP viscosities. △ Less

Submitted 8 October, 2020; originally announced October 2020.

Comments: 8 pages, 3 figures

Journal ref: Phys. Rev. Lett. 126, 242301 (2021)

Probing the multi-scale dynamical interaction between heavy quarks and the QGP using JETSCAPE

Authors: W. Fan, G. Vujanovic, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, R. Fries, C. Gale, F. Garza, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, W. Ke, E. Khalaj, B. Kim , et al. (25 additional authors not shown)

Abstract: The dynamics of shower development for a jet traveling through the QGP involves a variety of scales, one of them being the heavy quark mass. Even though the mass of the heavy quarks plays a subdominant role during the high virtuality portion of the jet evolution, it does affect longitudinal drag and diffusion, stimulating additional radiation from heavy quarks. These emissions partially compensate… ▽ More The dynamics of shower development for a jet traveling through the QGP involves a variety of scales, one of them being the heavy quark mass. Even though the mass of the heavy quarks plays a subdominant role during the high virtuality portion of the jet evolution, it does affect longitudinal drag and diffusion, stimulating additional radiation from heavy quarks. These emissions partially compensate the reduction in radiation from the dead cone effect. In the lower virtuality part of the shower, when the mass is comparable to the transverse momenta of the partons, scattering and radiation processes off heavy quarks differ from those off light quarks. All these factors result in a different nuclear modification factor for heavy versus light flavors and thus for heavy-flavor tagged jets. In this study, the heavy quark shower evolution and the fluid dynamical medium are modeled on an event by event basis using the JETSCAPE Framework. We present a multi-stage calculation that explores the differences between various heavy quark energy-loss mechanisms within a realistically expanding quark-gluon plasma (QGP). Inside the QGP, the highly virtual and energetic portion of the shower is modeled using the MATTER generator, while the LBT generator models the showers induced by energetic and close-to-on-shell heavy quarks. Energy-momentum exchange with the medium, essential for the study of jet modification, proceeds using a weak coupling recoil approach. The JETSCAPE framework allows for transitions, on the level of individual partons, from one energy-loss prescription to the other depending on the parton's energy and virtuality and the local density. This allows us to explore the effect and interplay between the different regimes of energy loss on the propagation and radiation from hard heavy quarks in a dense medium. △ Less

Submitted 9 August, 2022; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: arXiv admin note: substantial text overlap with arXiv:2002.06643

Photon-jet correlations in p-p and Pb-Pb collisions using JETSCAPE framework

Authors: C. Sirimanna, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, F. Garza, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, W. Ke, E. Khalaj, B. Kim , et al. (25 additional authors not shown)

Abstract: It is now well established that jet modification is a multistage effect; hence a single model alone cannot describe all facets of jet modification. The JETSCAPE framework is a multistage framework that uses several modules to simulate different stages of jet propagation through the QGP medium. These simulations require a set of parameters to ensure a smooth transition between stages. We fine tune… ▽ More It is now well established that jet modification is a multistage effect; hence a single model alone cannot describe all facets of jet modification. The JETSCAPE framework is a multistage framework that uses several modules to simulate different stages of jet propagation through the QGP medium. These simulations require a set of parameters to ensure a smooth transition between stages. We fine tune these parameters to successfully describe a variety of observables, such as the nuclear modification factors of leading hadrons and jets, jet shape, and jet fragmentation function. Photons can be produced in the hard scattering or as radiation from quarks inside jets. In this work, we study photon-jet transverse momentum imbalance and azimuthal correlation for both $p-p$ and $Pb-Pb$ collision systems. All the photons produced in each event, including the photons from hard scattering, radiation from the parton shower, and radiation from hadronization are considered with an isolation cut to directly compare with experimental data. The simulations are conducted using the same set of tuned parameters as used for the jet analysis. No new parameters are introduced or tuned. We demonstrate a significantly improved agreement with photons from $Pb-Pb$ collisions compared to prior efforts. This work provides an independent, parameter free verification of the multistage evolution framework. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 4 pages, 7 figures

First results from Hybrid Hadronization in small and large systems

Authors: M. Kordell II, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, F. Garza, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, W. Ke, E. Khalaj, B. Kim , et al. (25 additional authors not shown)

Abstract: "Hybrid Hadronization" is a new Monte Carlo package to hadronize systems of partons. It smoothly combines quark recombination applicable when distances between partons in phase space are small, and string fragmentation appropriate for dilute parton systems, following the picture outlined by Han et al. [PRC 93, 045207 (2016)]. Hybrid Hadronization integrates with PYTHIA 8 and can be applied to a va… ▽ More "Hybrid Hadronization" is a new Monte Carlo package to hadronize systems of partons. It smoothly combines quark recombination applicable when distances between partons in phase space are small, and string fragmentation appropriate for dilute parton systems, following the picture outlined by Han et al. [PRC 93, 045207 (2016)]. Hybrid Hadronization integrates with PYTHIA 8 and can be applied to a variety of systems from $e^++e^-$ to $A+A$ collisions. It takes systems of partons and their color flow information, for example from a Monte Carlo parton shower generator, as input. In addition, if for $A+A$ collisions a thermal background medium is provided, the package allows sampling thermal partons that contribute to hadronization. Hybrid Hadronization is available for use as a standalone code and is also part of JETSCAPE since the 2.0 release. In these proceedings we review the physics concepts underlying Hybrid Hadronization and demonstrate how users can use the code with various parton shower Monte Carlos. We present calculations of hadron chemistry and fragmentation functions in small and large systems when Hybrid Hadronization is combined with parton shower Monte Carlos MATTER and LBT. In particular, we discuss observable effects of the recombination of shower partons with thermal partons. △ Less

Submitted 11 September, 2020; v1 submitted 8 September, 2020; originally announced September 2020.

Comments: 4 pages, 3 figures, Proceedings of Hard Probes 2020, 1-6 June 2020, Austin, Texas; Updated Author list

Constraints on jet quenching from a multi-stage energy-loss approach

Authors: C. Park, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, F. Garza, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, W. Ke, E. Khalaj, B. Kim , et al. (25 additional authors not shown)

Abstract: We present a multi-stage model for jet evolution through a quark-gluon plasma within the JETSCAPE framework. The multi-stage approach in JETSCAPE provides a unified description of distinct phases in jet shower contingent on the virtuality. We demonstrate a simultaneous description of leading hadron and integrated jet observables as well as jet $v_n$ using tuned parameters. Medium response to the j… ▽ More We present a multi-stage model for jet evolution through a quark-gluon plasma within the JETSCAPE framework. The multi-stage approach in JETSCAPE provides a unified description of distinct phases in jet shower contingent on the virtuality. We demonstrate a simultaneous description of leading hadron and integrated jet observables as well as jet $v_n$ using tuned parameters. Medium response to the jet quenching is implemented based on a weakly-coupled recoil prescription. We also explore the cone-size dependence of jet energy loss inside the plasma. △ Less

Submitted 11 September, 2020; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: 4 pages, 4 figures, Proceedings of Hard Probes 2020, 1-6 June 2020, Austin, Texas

Fluctuation dynamics near the QCD critical point

Authors: Lipei Du, Ulrich Heinz, Krishna Rajagopal, Yi Yin

Abstract: The evolution of non-hydrodynamic slow processes near the QCD critical point is explored with the novel Hydro+ framework, which extends the conventional hydrodynamic description by coupling it to additional explicitly evolving slow modes describing long wavelength fluctuations. Their slow relaxation is controlled by critical behavior of the correlation length and is independent from gradients of m… ▽ More The evolution of non-hydrodynamic slow processes near the QCD critical point is explored with the novel Hydro+ framework, which extends the conventional hydrodynamic description by coupling it to additional explicitly evolving slow modes describing long wavelength fluctuations. Their slow relaxation is controlled by critical behavior of the correlation length and is independent from gradients of matter density and pressure that control the evolution of the hydrodynamic quantities. In this exploratory study we follow the evolution of the slow modes on top of a simplified QCD matter background, allowing us to clearly distinguish, and study both separately and in combination, the main effects controlling the dynamics of critical slow modes. In particular, we show how the evolution of the slow modes depend on their wave number, the expansion of and advection by the fluid background, and the behavior of the correlation length. Non-equilibrium contributions from the slow modes to bulk matter properties that affect the bulk dynamics (entropy, pressure, temperature and chemical potential) are discussed and found to be small. △ Less

Submitted 29 November, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: 29 pages, 18 figures. Some discussions are added; Fig. 1 is separated into two figures; in some figures, axis labels and line color styles are updated. Published as Phys. Rev. C 102, 054911 (2020)

Journal ref: Phys. Rev. C 102, 054911 (2020)

Hydrodynamic response to jets with a source based on causal diffusion

Authors: Y. Tachibana, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, K. Kauder, W. Ke, E. Khalaj, M. Kordell II , et al. (25 additional authors not shown)

Abstract: We study the medium response to jet evolution in the quark-gluon plasma within the JETSCAPE framework. Recoil partons' medium response in the weakly coupled description is implemented in the multi-stage jet energy-loss model in the framework. As a further extension, the hydrodynamic description is rearranged to include in-medium jet transport based on a strong-coupling picture. To interface hydrod… ▽ More We study the medium response to jet evolution in the quark-gluon plasma within the JETSCAPE framework. Recoil partons' medium response in the weakly coupled description is implemented in the multi-stage jet energy-loss model in the framework. As a further extension, the hydrodynamic description is rearranged to include in-medium jet transport based on a strong-coupling picture. To interface hydrodynamics with jet energy-loss models, the hydrodynamic source term is modeled by a causal formulation employing the relativistic diffusion equation. The jet shape and fragmentation function are studied via realistic simulations with weakly coupled recoils. We also demonstrate modifications in the medium caused by the hydrodynamic response. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: 4 pages, 2 figures, contribution to the Quark Matter 2019 proceedings

Jet quenching in a multi-stage Monte Carlo approach

Authors: A. Kumar, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, K. Kauder, W. Ke, E. Khalaj, M. Kordell II , et al. (25 additional authors not shown)

Abstract: We present a jet quenching model within a unified multi-stage framework and demonstrate for the first time a simultaneous description of leading hadrons, inclusive jets, and elliptic flow observables which spans multiple centralities and collision energies. This highlights one of the major successes of the JETSCAPE framework in providing a tool for setting up an effective parton evolution that inc… ▽ More We present a jet quenching model within a unified multi-stage framework and demonstrate for the first time a simultaneous description of leading hadrons, inclusive jets, and elliptic flow observables which spans multiple centralities and collision energies. This highlights one of the major successes of the JETSCAPE framework in providing a tool for setting up an effective parton evolution that includes a high-virtuality radiation dominated energy loss phase (MATTER), followed by a low-virtuality scattering dominated (LBT) energy loss phase. Measurements of jet and charged-hadron $R_{AA}$ set strong constraints on the jet quenching model. Jet-medium response is also included through a weakly-coupled transport description. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 4 pages, 4 figures, contribution to the Quark Matter 2019 proceedings

Multi-stage evolution of heavy quarks in the quark-gluon plasma

Authors: G. Vujanovic, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, K. Kauder, W. Ke, E. Khalaj, M. Kordell II , et al. (25 additional authors not shown)

Abstract: The interaction of heavy flavor with the quark-gluon plasma (QGP) in relativistic heavy-ion collisions is studied using JETSCAPE, a publicly available software package containing a framework for Monte Carlo event generators. Multi-stage (and multi-model) evolution of heavy quarks within JETSCAPE provides a cohesive description of heavy flavor quenching inside the QGP. As the parton shower develops… ▽ More The interaction of heavy flavor with the quark-gluon plasma (QGP) in relativistic heavy-ion collisions is studied using JETSCAPE, a publicly available software package containing a framework for Monte Carlo event generators. Multi-stage (and multi-model) evolution of heavy quarks within JETSCAPE provides a cohesive description of heavy flavor quenching inside the QGP. As the parton shower develops, a model becomes active as soon as its kinematic region of validity is reached. Two combinations of heavy-flavor energy-loss models are explored within a realistic QGP medium, using parameters which were tuned to describe {\it light-flavor} partonic energy-loss. △ Less

Submitted 16 February, 2020; originally announced February 2020.

Comments: 4 pages, 2 figures, contribution to the Quark Matter 2019 proceedings

Revisiting Bayesian constraints on the transport coefficients of QCD

Authors: J. -F. Paquet, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, K. Kauder, W. Ke, E. Khalaj, M. Kordell II , et al. (25 additional authors not shown)

Abstract: Multistage models based on relativistic viscous hydrodynamics have proven successful in describing hadron measurements from relativistic nuclear collisions. These measurements are sensitive to the shear and the bulk viscosities of QCD and provide a unique opportunity to constrain these transport coefficients. Bayesian analyses can be used to obtain systematic constraints on the viscosities of QCD,… ▽ More Multistage models based on relativistic viscous hydrodynamics have proven successful in describing hadron measurements from relativistic nuclear collisions. These measurements are sensitive to the shear and the bulk viscosities of QCD and provide a unique opportunity to constrain these transport coefficients. Bayesian analyses can be used to obtain systematic constraints on the viscosities of QCD, through methodical model-to-data comparisons. In this manuscript, we discuss recent developments in Bayesian analyses of heavy ion collision data. We highlight the essential role of closure tests in validating a Bayesian analysis before comparison with measurements. We discuss the role of the emulator that is used as proxy for the multistage theoretical model. We use an ongoing Bayesian analysis of soft hadron measurements by the JETSCAPE Collaboration as context for the discussion. △ Less

Submitted 12 February, 2020; originally announced February 2020.

Comments: 4 pages, 1 figure, contribution to the Quark Matter 2019 proceedings

The JETSCAPE framework: p+p results

Authors: A. Kumar, Y. Tachibana, D. Pablos, C. Sirimanna, R. J. Fries, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, H. Elfner, D. Everett, W. Fan, C. Gale, Y. He, U. Heinz, B. V. Jacak, P. M. Jacobs, 15 S. Jeon, K. Kauder, W. Ke, E. Khalaj , et al. (21 additional authors not shown)

Abstract: The JETSCAPE framework is a modular and versatile Monte Carlo software package for the simulation of high energy nuclear collisions. In this work we present a new tune of JETSCAPE, called PP19, and validate it by comparison to jet-based measurements in $p+p$ collisions, including inclusive single jet cross sections, jet shape observables, fragmentation functions, charged hadron cross sections, and… ▽ More The JETSCAPE framework is a modular and versatile Monte Carlo software package for the simulation of high energy nuclear collisions. In this work we present a new tune of JETSCAPE, called PP19, and validate it by comparison to jet-based measurements in $p+p$ collisions, including inclusive single jet cross sections, jet shape observables, fragmentation functions, charged hadron cross sections, and dijet mass cross sections. These observables in $p+p$ collisions provide the baseline for their counterparts in nuclear collisions. Quantifying the level of agreement of JETSCAPE results with $p+p$ data is thus necessary for meaningful applications of JETSCAPE to A+A collisions. The calculations use the JETSCAPE PP19 tune, defined in this paper, based on version 1.0 of the JETSCAPE framework. For the observables discussed in this work calculations using JETSCAPE PP19 agree with data over a wide range of collision energies at a level comparable to standard Monte Carlo codes. These results demonstrate the physics capabilities of the JETSCAPE framework and provide benchmarks for JETSCAPE users. △ Less

Submitted 6 November, 2019; v1 submitted 12 October, 2019; originally announced October 2019.

Comments: 23 pages, 23 figures; v1.1: minor bug fixes in author information

Journal ref: Phys. Rev. C 102, 054906 (2020)

(3+1)-dimensional dissipative relativistic fluid dynamics at non-zero net baryon density

Authors: Lipei Du, Ulrich Heinz

Abstract: Heavy-ion collisions at center-of-mass energies between 1 and 100 GeV/nucleon are essential to understand the phase diagram of QCD and search for its critical point. At these energies the net baryon density of the system can be high, and simulating its evolution becomes an indispensable part of theoretical modeling. We here present the (3+1)-dimensional diffusive relativistic hydrodynamic code BES… ▽ More Heavy-ion collisions at center-of-mass energies between 1 and 100 GeV/nucleon are essential to understand the phase diagram of QCD and search for its critical point. At these energies the net baryon density of the system can be high, and simulating its evolution becomes an indispensable part of theoretical modeling. We here present the (3+1)-dimensional diffusive relativistic hydrodynamic code BEShydro which solves the equations of motion of second-order Denicol-Niemi-Molnar-Rischke (DNMR) theory, including bulk and shear viscous currents and baryon diffusion currents. BEShydro features a modular structure that allows to easily turn on and off baryon evolution and different dissipative effects and thus to study their physical effects on the dynamical evolution individually. An extensive set of test protocols for the code, including several novel tests of the precision of baryon transport that can also be used to test other such codes, is documented here and supplied as a permanent part of the code package. △ Less

Submitted 8 December, 2019; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: 62 pages, 8 figures; small changes, accepted version

Journal ref: Comput.Phys.Commun. 251 (2020) 107090

Universal scaling of meson and baryon spectra in p-Pb collisions at 5.02 TeV

Authors: N. Liu, X. L. Du, L. Y. Qiao, G. R. Che, W. C. Zhang

Abstract: We systematically investigate the scaling property of mesons (pions and kaons) and baryons (protons, ${\rm Λらむだ}$, ${\rm Ξくしー}$ and ${\rm Ωおめが}$) transverse momentum ($p_{\rm T}$) spectra at different centrality classes (0-5$\%$, 5-10$\%$, 10-20$\%$, 20-40$\%$, 40-60$\%$, 60-80$\%$ and 80-100$\%$) in proton-lead collisions with center of mass energy per nucleon pair 5.02 TeV. In the low $p_{\rm T}$ region… ▽ More We systematically investigate the scaling property of mesons (pions and kaons) and baryons (protons, ${\rm Λらむだ}$, ${\rm Ξくしー}$ and ${\rm Ωおめが}$) transverse momentum ($p_{\rm T}$) spectra at different centrality classes (0-5$\%$, 5-10$\%$, 10-20$\%$, 20-40$\%$, 40-60$\%$, 60-80$\%$ and 80-100$\%$) in proton-lead collisions with center of mass energy per nucleon pair 5.02 TeV. In the low $p_{\rm T}$ region with $p_{\rm{T}}\leq $ 3.9 (3.1, 2.5, 2.7, 2.4 and 2.8) GeV/c, a universal scaling independent of the centrality is observed in the pion (kaon, proton, ${\rm Λらむだ}$, ${\rm Ξくしー}$ and ${\rm Ωおめが}$) spectra when a dilatation, $p_{\rm T}\rightarrow p_{\rm T}/K$, is applied. Here $K$ is a scaling parameter depending on the centrality class. We find that the rates at which ln$K$ changes with the logarithmic value of the average value of the number of participating nucleons, ln$\langle N_{\rm{part}}\rangle$, are stronger for baryons than those for mesons. In the high $p_{\rm T}$ region, there is a deviation from the scaling. The more peripheral the collisions are, the more obvious the violation of the scaling is. In the framework of the colour string percolation (CSP) model, we show that mesons and baryons are generated from the decay of clusters formed by strings overlapping in the transverse plane with the same size dispersion but with different mean size. The mean size of clusters for baryons is smaller than that of mesons. For the same hadrons at different centrality classes, the mean size of clusters decreases with the increase of centrality. The fragmentation functions for cluster decay are different for different hadrons, while they are universal for the same hadrons at different centrality classes. The universal scaling of the meson and baryon spectra in the low $p_{\rm T}$ region can be quantitatively understood with the CSP model at the same time. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comments: 17 pages, 8 figs.,6 tables. Accepted by International Journal of Modern Physics E. arXiv admin note: text overlap with arXiv:1805.02300

Journal ref: International Journal of Modern Physics E, Vol. 28, No. 4 (2019) 1950030

The JETSCAPE framework

Authors: J. H. Putschke, K. Kauder, E. Khalaj, A. Angerami, S. A. Bass, S. Cao, J. Coleman, L. Cunqueiro, T. Dai, L. Du, H. Elfner, D. Everett, W. Fan, R. J. Fries, C. Gale, Y. He, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, W. Ke, M. Kordell II, A. Kumar, T. Luo, A. Majumder , et al. (20 additional authors not shown)

Abstract: The JETSCAPE simulation framework is an overarching computational envelope for developing complete event generators for heavy-ion collisions. It allows for modular incorporation of a wide variety of existing and future software that simulates different aspects of a heavy-ion collision. The default JETSCAPE package contains both the framework, and an entire set of indigenous and third party routine… ▽ More The JETSCAPE simulation framework is an overarching computational envelope for developing complete event generators for heavy-ion collisions. It allows for modular incorporation of a wide variety of existing and future software that simulates different aspects of a heavy-ion collision. The default JETSCAPE package contains both the framework, and an entire set of indigenous and third party routines that can be used to directly compare with experimental data. In this article, we outline the algorithmic design of the JETSCAPE framework, define the interfaces and describe the default modules required to carry out full simulations of heavy-ion collisions within this package. We begin with a description of the various physics elements required to simulate an entire event in a heavy-ion collision, and distribute these within a flowchart representing the event generator and statistical routines for comparison with data. This is followed by a description of the abstract class structure, with associated members and functions required for this flowchart to work. We then define the interface that will be required for external users of JETSCAPE to incorporate their code within this framework and to modify existing elements within the default distribution. We conclude with a discussion of some of the physics output for both $p$-$p$ and $A$-$A$ collisions from the default distribution, and an outlook towards future releases. In the appendix, we discuss various architectures on which this code can be run and outline our benchmarks on similar hardware. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: 93 pages, 13 figures