Search

To study the bulk properties of quark-gluon plasma produced at beam-energy scan (BES) energies at the Relativistic Heavy Ion Collider (RHIC), we extend the (3+1)-dimensional viscous hydrodynamics clvisc to include net baryon number conservation and Israel-Stewart-like equations for baryon diffusion with the neos-bqs equation of state, fluctuating initial conditions from the Monte Carlo Glauber model and the afterburner smash. This integrated framework is shown to provide a good description of identified particle spectra, mean transverse momenta, and anisotropic flows for different centralities and over a wide range of collision energies (7.7-62.4 GeV). It is found that the mean momenta of identified particles and anisotropic flows increase mildly with the collision energy due to larger radial flow. We further compute the multiple-particle cumulant ratio v2{4}/v2{2} of elliptic flow across BES energies, and find that the relative fluctuations of elliptic flow are insensitive to the collision energy, consistent with the preliminary STAR data. Our model provides a benchmark for understanding the RHIC-BES data and studying the critical properties and phase structure of hot and dense QCD matter.

We develop a new heavy quark transport model, QLBT, to simulate the dynamical propagation of heavy quarks inside the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions. Our QLBT model is based on the linear Boltzmann transport (LBT) model with the ideal QGP replaced by a collection of quasi-particles to account for the non-perturbative interactions among quarks and gluons of the hot QGP. The thermal masses of quasi-particles are fitted to the equation of state from lattice QCD simulations using the Bayesian statistical analysis method. Combining QLBT with our advanced hybrid fragmentation-coalescence hadronization approach, we calculate the nuclear modification factor RAA and the elliptic flow v2 of D mesons at the Relativistic Heavy-Ion Collider and the Large Hadron Collider. By comparing our QLBT calculation to the experimental data on the D meson RAA and v2, we extract the heavy quark transport parameter q^ and diffusion coefficient Ds in the temperature range of 1-4Tc, and compare them with the lattice QCD results and other phenomenological studies.

Systems with different interactions could develop the same critical behavior due to the underlying symmetry and universality. Using this principle of universality, we can embed critical correlations modeled on the 3D Ising model into the simulated data of heavy-ion collisions, hiding weak signals of a few inter-particle correlations within a large particle cloud. Employing a point cloud network with dynamical edge convolution, we are able to identify events with critical fluctuations through supervised learning, and pick out a large fraction of signal particles used for decision-making in each single event.

An extraordinary strong magnetic field, eB0≈1018 Gauss, is expected to be generated in non-central ultrarelativistic heavy ion collisions and it is envisaged to induce several effects on hot QCD matter including the possibility of local parity and local charge conjugation and parity symmetry violations. A direct signature of such e.m. fields and a first quantitative measurement of its strength and lifetime are still missing. We point out that both the mean value of leptonic invariant mass of Z0 boson, reconstructed by its decaying lepton pairs, and the relative width are modified in relativistic heavy ion collisions due to the presence of strong initial e.m. fields. We propose a measurement of the leptonic invariant mass of Z0 as a novel probe of the strength of the By. Both shifts could be up to about few hundred MeV and are found to depend on the integral of By over the time duration quadratically (approximate). Hence it provides a novel and clear probe of electromagnetic fields, which can be tested experimentally.

Advances in machine learning methods provide tools that have broad applicability in scientific research. These techniques are being applied across the diversity of nuclear physics research topics, leading to advances that will facilitate scientific discoveries and societal applications. This Colloquium provides a snapshot of nuclear physics research, which has been transformed by machine learning techniques.

Large angle gluon radiations induced by multiple parton scatterings contribute to dijet production in deeply inelastic scattering off a large nucleus at the electron-ion collider. Within the generalized high-twist approach to multiple parton scattering, contributions at the leading order and twist four in perturbative QCD and large Bjorken momentum fraction xB can be expressed as a convolution of the multiple parton scattering amplitudes and the transverse momentum dependent (TMD) two-parton correlation matrix elements. We study this medium-induced dijet spectrum and its azimuthal angle correlation under the approximation of small longitudinal momentum transfer in the secondary scattering and the factorization of two-parton correlation matrix elements as a product of the quark and gluon TMD parton distribution function (PDF). Contributions to dijet cross section from double scattering are power suppressed and only become sizable for minijets at small transverse momentum. We find that the total dijet correlation for these minijets, which also includes the contribution from single scattering, is sensitive to the transverse momentum broadening in the quark TMD PDF at large x and saturation in the gluon TMD PDF at small x inside the nucleus. The correlation from double scattering is also found to increase with the dijet rapidity gap and have a quadratic nuclear-size dependence because of the Landau-Pomeranchuk-Migdal interference in gluon emission induced by multiple scattering. Experimental measurements of such unique features in the dijet correlation can shed light on the Landau-Pomeranchuk-Migdal interference in strong interaction and gluon saturation in large nuclei.

An extraordinary strong magnetic field, eB0≈1018 Gauss, is expected to be generated in non-central ultrarelativistic heavy ion collisions and it is envisaged to induce several effects on hot QCD matter including the possibility of local parity and local charge conjugation and parity symmetry violations. A direct signature of such e.m. fields and a first quantitative measurement of its strength and lifetime are still missing. We point out that both the mean value of leptonic invariant mass of Z0 boson, reconstructed by its decaying lepton pairs, and the relative width are modified in relativistic heavy ion collisions due to the presence of strong initial e.m. fields. We propose a measurement of the leptonic invariant mass of Z0 as a novel probe of the strength of the By. Both shifts could be up to about few hundred MeV and are found to depend on the integral of By over the time duration quadratically (approximate). Hence it provides a novel and clear probe of electromagnetic fields, which can be tested experimentally.

Systems with different interactions could develop the same critical behavior due to the underlying symmetry and universality. Using this principle of universality, we can embed critical correlations modeled on the 3D Ising model into the simulated data of heavy-ion collisions, hiding weak signals of a few inter-particle correlations within a large particle cloud. Employing a point cloud network with dynamical edge convolution, we are able to identify events with critical fluctuations through supervised learning, and pick out a large fraction of signal particles used for decision-making in each single event.

To study the bulk properties of quark-gluon plasma produced at beam-energy scan (BES) energies at the Relativistic Heavy Ion Collider (RHIC), we extend the (3+1)-dimensional viscous hydrodynamics clvisc to include net baryon number conservation and Israel-Stewart-like equations for baryon diffusion with the neos-bqs equation of state, fluctuating initial conditions from the Monte Carlo Glauber model and the afterburner smash. This integrated framework is shown to provide a good description of identified particle spectra, mean transverse momenta, and anisotropic flows for different centralities and over a wide range of collision energies (7.7-62.4 GeV). It is found that the mean momenta of identified particles and anisotropic flows increase mildly with the collision energy due to larger radial flow. We further compute the multiple-particle cumulant ratio v2{4}/v2{2} of elliptic flow across BES energies, and find that the relative fluctuations of elliptic flow are insensitive to the collision energy, consistent with the preliminary STAR data. Our model provides a benchmark for understanding the RHIC-BES data and studying the critical properties and phase structure of hot and dense QCD matter.

We develop a new heavy quark transport model, QLBT, to simulate the dynamical propagation of heavy quarks inside the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions. Our QLBT model is based on the linear Boltzmann transport (LBT) model with the ideal QGP replaced by a collection of quasi-particles to account for the non-perturbative interactions among quarks and gluons of the hot QGP. The thermal masses of quasi-particles are fitted to the equation of state from lattice QCD simulations using the Bayesian statistical analysis method. Combining QLBT with our advanced hybrid fragmentation-coalescence hadronization approach, we calculate the nuclear modification factor RAA and the elliptic flow v2 of D mesons at the Relativistic Heavy-Ion Collider and the Large Hadron Collider. By comparing our QLBT calculation to the experimental data on the D meson RAA and v2, we extract the heavy quark transport parameter q^ and diffusion coefficient Ds in the temperature range of 1-4Tc, and compare them with the lattice QCD results and other phenomenological studies.

In this proceeding, we review our recent work using deep convolutional neural network (CNN) to identify the nature of the QCD transition in a hybrid modeling of heavy-ion collisions. Within this hybrid model, a viscous hydrodynamic model is coupled with a hadronic cascade “after-burner”. As a binary classification setup, we employ two different types of equations of state (EoS) of the hot medium in the hydrodynamic evolution. The resulting final-state pion spectra in the transverse momentum and azimuthal angle plane are fed to the neural network as the input data in order to distinguish different EoS. To probe the effects of the fluctuations in the event-by-event spectra, we explore different scenarios for the input data and make a comparison in a systematic way. We observe a clear hierarchy in the predictive power when the network is fed with the event-by-event, cascade-coarse-grained and event-fine-averaged spectra. The carefully-trained neural network can extract high-level features from pion spectra to identify the nature of the QCD transition in a realistic simulation scenario.

We report a new determination of q, the jet transport coefficient of the quark-gluon plasma. We use the JETSCAPE framework, which incorporates a novel multistage theoretical approach to in-medium jet evolution and Bayesian inference for parameter extraction. The calculations, based on the Matter and Lbt jet quenching models, are compared to experimental measurements of inclusive hadron suppression in Au+Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and Pb+Pb collisions at the CERN Large Hadron Collider (LHC). The correlation of experimental systematic uncertainties is accounted for in the parameter extraction. The functional dependence of q on jet energy or virtuality and medium temperature is based on a perturbative picture of in-medium scattering, with components reflecting the different regimes of applicability of Matter and Lbt. In the multistage approach, the switch between Matter and Lbt is governed by a virtuality scale Q0. 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 q/T3 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 Q0, has 90% CR of 240 GeV/c. The value of Q0, determined here for the first time, is in the range 2.0-2.7 GeV.

We study the properties of the strongly coupled quark-gluon plasma with a multistage model of heavy-ion collisions that combines the TRENTo 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 for 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 (LHC) and the Relativistic Heavy Ion Collider (RHIC), 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≈ 150-250 MeV, but their constraining power degrades at higher temperatures, T 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 TRENTo 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.

We report a new determination of q, the jet transport coefficient of the quark-gluon plasma. We use the JETSCAPE framework, which incorporates a novel multistage theoretical approach to in-medium jet evolution and Bayesian inference for parameter extraction. The calculations, based on the Matter and Lbt jet quenching models, are compared to experimental measurements of inclusive hadron suppression in Au+Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and Pb+Pb collisions at the CERN Large Hadron Collider (LHC). The correlation of experimental systematic uncertainties is accounted for in the parameter extraction. The functional dependence of q on jet energy or virtuality and medium temperature is based on a perturbative picture of in-medium scattering, with components reflecting the different regimes of applicability of Matter and Lbt. In the multistage approach, the switch between Matter and Lbt is governed by a virtuality scale Q0. 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 q/T3 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 Q0, has 90% CR of 240 GeV/c. The value of Q0, determined here for the first time, is in the range 2.0-2.7 GeV.

The large values and constituent-quark-number scaling of the elliptic flow of low- D mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through strong interactions with quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low- charm quarks, we compare the meson spectra and elliptic flow from a hydrodynamic model to experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in the local comoving frame of the medium before freeze-out. mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both the hydrodynamic and transport models can describe the elliptic flow of D0 mesons at GeV/c as measured in Au+Au collisions at sqrts_{NN}} = 200 GeV. Though the experimental data on spectra are consistent with the hydrodynamic result at small p_T\sim 1 GeV/c, they deviate from the hydrodynamic model at high transverse momentum, GeV/c. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Our comparative study indicates that charm quarks only approach local thermal equilibrium at small, even though they acquire sizable elliptic flow comparable to light-quark hadrons at both small and intermediate.

It is widely acknowledged that heavy flavor probes are sensitive to the properties of the quark-gluon plasma and are often considered an important tool for the plasma tomography studies. Forward rapidity observables can provide further insight on the dynamics of the medium due to the interplay between the medium size and the differences in the production spectra of heavy quark probes. In this proceedings we present the nuclear modification factor RAA's for B and D mesons, as well as heavy flavor leptons, in the rapidity range −4.0 < y < 4.0 obtained from relativistic Langevin equation with gluon radiation coupled with a (3+1)-dimensional viscous hydrodynamics medium background. We present comparison with experimental data at mid-rapidity as well as predictions for different rapidity ranges.

The CLVisc (3+1)D viscous hydrodynamic model is extended to include the equation of net baryon conservation and the Israel-Stewart-like equations for dissipative baryon current. Using the NEOSB equation of state, we simulate the dynamical evolution and collectivity of the quark-gluon plasma with finite chemical potential, assuming smooth energy density and net baryon density distributions at the initial proper time. Numerical results are shown for the impact of net-baryon dissipation on particle yields and pT spectra in heavy-ion collisions at beam energy scan energies.

We calculate the radiative parton energy loss in the deeply inelastic scattering (DIS) off a large nucleus within a generalized high-twist approach. The final gluon radiation spectrum is a convolution of the hard partonic part and the transverse momentum dependent (TMD) quark-gluon correlation function, without the twist expansion in the transverse momentum of the initial gluons used in original high twist approach. The TMD quark-gluon correlation function can be factorized approximately as the product of initial quark distribution and TMD gluon distribution which can be used to define the generalized or TMD jet transport coefficient. The radiation spectrum will recover the Gylassy-Levai-Vitev (GLV) result in the first order of the opacity expansion, under the static scattering center and soft gluon radiation approximation. We also investigate numerically the difference as a result of the soft gluon radiation approximation, under the static scattering center approximation.

We study the properties of the strongly coupled quark-gluon plasma with a multistage model of heavy-ion collisions that combines the TRENTo 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 for 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 (LHC) and the Relativistic Heavy Ion Collider (RHIC), 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≈ 150-250 MeV, but their constraining power degrades at higher temperatures, T 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 TRENTo 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.

We report a new determination of q, the jet transport coefficient of the quark-gluon plasma. We use the JETSCAPE framework, which incorporates a novel multistage theoretical approach to in-medium jet evolution and Bayesian inference for parameter extraction. The calculations, based on the Matter and Lbt jet quenching models, are compared to experimental measurements of inclusive hadron suppression in Au+Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and Pb+Pb collisions at the CERN Large Hadron Collider (LHC). The correlation of experimental systematic uncertainties is accounted for in the parameter extraction. The functional dependence of q on jet energy or virtuality and medium temperature is based on a perturbative picture of in-medium scattering, with components reflecting the different regimes of applicability of Matter and Lbt. In the multistage approach, the switch between Matter and Lbt is governed by a virtuality scale Q0. 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 q/T3 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 Q0, has 90% CR of 240 GeV/c. The value of Q0, determined here for the first time, is in the range 2.0-2.7 GeV.

The CLVisc (3+1)D viscous hydrodynamic model is extended to include the equation of net baryon conservation and the Israel-Stewart-like equations for dissipative baryon current. Using the NEOSB equation of state, we simulate the dynamical evolution and collectivity of the quark-gluon plasma with finite chemical potential, assuming smooth energy density and net baryon density distributions at the initial proper time. Numerical results are shown for the impact of net-baryon dissipation on particle yields and pT spectra in heavy-ion collisions at beam energy scan energies.

We calculate the radiative parton energy loss in the deeply inelastic scattering (DIS) off a large nucleus within a generalized high-twist approach. The final gluon radiation spectrum is a convolution of the hard partonic part and the transverse momentum dependent (TMD) quark-gluon correlation function, without the twist expansion in the transverse momentum of the initial gluons used in original high twist approach. The TMD quark-gluon correlation function can be factorized approximately as the product of initial quark distribution and TMD gluon distribution which can be used to define the generalized or TMD jet transport coefficient. The radiation spectrum will recover the Gylassy-Levai-Vitev (GLV) result in the first order of the opacity expansion, under the static scattering center and soft gluon radiation approximation. We also investigate numerically the difference as a result of the soft gluon radiation approximation, under the static scattering center approximation.

It is widely acknowledged that heavy flavor probes are sensitive to the properties of the quark-gluon plasma and are often considered an important tool for the plasma tomography studies. Forward rapidity observables can provide further insight on the dynamics of the medium due to the interplay between the medium size and the differences in the production spectra of heavy quark probes. In this proceedings we present the nuclear modification factor RAA's for B and D mesons, as well as heavy flavor leptons, in the rapidity range −4.0 < y < 4.0 obtained from relativistic Langevin equation with gluon radiation coupled with a (3+1)-dimensional viscous hydrodynamics medium background. We present comparison with experimental data at mid-rapidity as well as predictions for different rapidity ranges.

The large values and constituent-quark-number scaling of the elliptic flow of low- D mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through strong interactions with quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low- charm quarks, we compare the meson spectra and elliptic flow from a hydrodynamic model to experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in the local comoving frame of the medium before freeze-out. mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both the hydrodynamic and transport models can describe the elliptic flow of D0 mesons at GeV/c as measured in Au+Au collisions at sqrts_{NN}} = 200 GeV. Though the experimental data on spectra are consistent with the hydrodynamic result at small p_T\sim 1 GeV/c, they deviate from the hydrodynamic model at high transverse momentum, GeV/c. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Our comparative study indicates that charm quarks only approach local thermal equilibrium at small, even though they acquire sizable elliptic flow comparable to light-quark hadrons at both small and intermediate.

In this proceeding, we review our recent work using deep convolutional neural network (CNN) to identify the nature of the QCD transition in a hybrid modeling of heavy-ion collisions. Within this hybrid model, a viscous hydrodynamic model is coupled with a hadronic cascade “after-burner”. As a binary classification setup, we employ two different types of equations of state (EoS) of the hot medium in the hydrodynamic evolution. The resulting final-state pion spectra in the transverse momentum and azimuthal angle plane are fed to the neural network as the input data in order to distinguish different EoS. To probe the effects of the fluctuations in the event-by-event spectra, we explore different scenarios for the input data and make a comparison in a systematic way. We observe a clear hierarchy in the predictive power when the network is fed with the event-by-event, cascade-coarse-grained and event-fine-averaged spectra. The carefully-trained neural network can extract high-level features from pion spectra to identify the nature of the QCD transition in a realistic simulation scenario.

Using deep convolutional neural network (CNN), the nature of the QCD transition can be identified from the final-state pion spectra from hybrid model simulations of heavy-ion collisions that combines a viscous hydrodynamic model with a hadronic cascade “after-burner”. Two different types of equations of state (EoS) of the medium are used in the hydrodynamic evolution. The resulting spectra in transverse momentum and azimuthal angle are used as the input data to train the neural network to distinguish different EoS. Different scenarios for the input data are studied and compared in a systematic way. A clear hierarchy is observed in the prediction accuracy when using the event-by-event, cascade-coarse-grained and event-fine-averaged spectra as input for the network, which are about 80%, 90% and 99%, respectively. A comparison with the prediction performance by deep neural network (DNN) with only the normalized pion transverse momentum spectra is also made. High-level features of pion spectra captured by a carefully-trained neural network were found to be able to distinguish the nature of the QCD transition even in a simulation scenario which is close to the experiments.

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.

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.

Using deep convolutional neural network (CNN), the nature of the QCD transition can be identified from the final-state pion spectra from hybrid model simulations of heavy-ion collisions that combines a viscous hydrodynamic model with a hadronic cascade “after-burner”. Two different types of equations of state (EoS) of the medium are used in the hydrodynamic evolution. The resulting spectra in transverse momentum and azimuthal angle are used as the input data to train the neural network to distinguish different EoS. Different scenarios for the input data are studied and compared in a systematic way. A clear hierarchy is observed in the prediction accuracy when using the event-by-event, cascade-coarse-grained and event-fine-averaged spectra as input for the network, which are about 80%, 90% and 99%, respectively. A comparison with the prediction performance by deep neural network (DNN) with only the normalized pion transverse momentum spectra is also made. High-level features of pion spectra captured by a carefully-trained neural network were found to be able to distinguish the nature of the QCD transition even in a simulation scenario which is close to the experiments.

We revisit radiative parton energy loss in deeply inelastic scattering (DIS) off a large nucleus within the perturbative QCD approach. We calculate the gluon radiation spectra induced by double parton scattering in DIS without collinear expansion in the transverse momentum of initial gluons as in the original high-twist approach. The final radiative gluon spectrum can be expressed in terms of the convolution of hard partonic parts and unintegrated or transverse momentum dependent (TMD) quark-gluon correlations. The TMD quark-gluon correlation can be factorized approximately as a product of initial quark distribution and TMD gluon distribution which can be used to define the generalized or TMD jet transport coefficient. Under the static scattering center and soft radiative gluon approximation, we recover the result by Gylassy-Levai-Vitev in the first order of the opacity expansion. The difference as a result of the soft radiative gluon approximation is investigated numerically under the static scattering center approximation.

We study the longitudinal decorrelations of elliptic, triangular and quadrangular flows in heavy-ion collisions at the LHC and RHIC energies. The event-by-event CLVisc (3+1)-dimensional hydrodynamics model, combined with the fully fluctuating AMPT initial conditions, is utilized to simulate the space-time evolution of the strongly-coupled quark-gluon plasma. Detailed analysis is performed for the longitudinal decorrelations of flow vectors, flow magnitudes and flow orientations. We find strong correlations between final-state longitudinal decorrelations of anisotropic flows and initial-state longitudinal structures and collision geometry: while the decorrelation of elliptic flow shows a non- monotonic centrality dependence due to initial elliptic geometry, typically the longitudinal flow decorrelations are larger in lower energy and less central collisions where the mean lengths of the string structure are shorter in the initial states.

We study the longitudinal decorrelations of elliptic, triangular and quadrangular flows in heavy-ion collisions at the LHC and RHIC energies. The event-by-event CLVisc (3+1)-dimensional hydrodynamics model, combined with the fully fluctuating AMPT initial conditions, is utilized to simulate the space-time evolution of the strongly-coupled quark-gluon plasma. Detailed analysis is performed for the longitudinal decorrelations of flow vectors, flow magnitudes and flow orientations. We find strong correlations between final-state longitudinal decorrelations of anisotropic flows and initial-state longitudinal structures and collision geometry: while the decorrelation of elliptic flow shows a non- monotonic centrality dependence due to initial elliptic geometry, typically the longitudinal flow decorrelations are larger in lower energy and less central collisions where the mean lengths of the string structure are shorter in the initial states.

We revisit radiative parton energy loss in deeply inelastic scattering (DIS) off a large nucleus within the perturbative QCD approach. We calculate the gluon radiation spectra induced by double parton scattering in DIS without collinear expansion in the transverse momentum of initial gluons as in the original high-twist approach. The final radiative gluon spectrum can be expressed in terms of the convolution of hard partonic parts and unintegrated or transverse momentum dependent (TMD) quark-gluon correlations. The TMD quark-gluon correlation can be factorized approximately as a product of initial quark distribution and TMD gluon distribution which can be used to define the generalized or TMD jet transport coefficient. Under the static scattering center and soft radiative gluon approximation, we recover the result by Gylassy-Levai-Vitev in the first order of the opacity expansion. The difference as a result of the soft radiative gluon approximation is investigated numerically under the static scattering center approximation.

We report on broadly based systematic investigations of the modeling components for open heavy-flavor diffusion and energy loss in strongly interacting matter in their application to heavy-flavor observables in high-energy heavy–ion collisions, conducted within an EMMI Rapid Reaction Task Force framework. Initial spectra including cold-nuclear-matter effects, a wide variety of space-time evolution models, heavy-flavor transport coefficients, and hadronization mechanisms are scrutinized in an effort to quantify pertinent uncertainties in the calculations of nuclear modification factors and elliptic flow of open heavy-flavor particles in nuclear collisions. We develop procedures for error assessments and criteria for common model components to improve quantitative estimates for the (low-momentum) heavy-flavor diffusion coefficient as a long-wavelength characteristic of QCD matter as a function of temperature, and for energy loss coefficients of high-momentum heavy-flavor particles.

We revisit the calculation of multiple parton scattering of a heavy quark in nuclei within the framework of recently improved high-twist factorization formalism, in which gauge invariance is ensured by a delicate setup of the initial partons' transverse momenta. We derive a new result for medium modified heavy quark fragmentation functions in deeply inelastic scattering. It is consistent with the previous calculation of light quark energy loss in the massless limit, but leads to a new correction term in the heavy quark case, which vanishes in the soft gluon radiation limit. We show numerically the significance of the new correction term in the calculation of heavy quark energy loss as compared to previous studies and with soft gluon radiation approximation.

We develop a general formalism for the quantum kinetics of chiral fermions in a background electromagnetic field based on a semiclassical expansion of covariant Wigner functions in the Planck constant. We demonstrate to any order of that only the time-component of the Wigner function is independent while other components are explicit derivative. We further demonstrate to any order of that a system of quantum kinetic equations for multiple-components of Wigner functions can be reduced to one chiral kinetic equation involving only the single-component distribution function. These are remarkable properties of the quantum kinetics of chiral fermions and will significantly simplify the description and simulation of chiral effects in heavy ion collisions and Dirac/Weyl semimetals. We present the unintegrated chiral kinetic equations in four-momenta up to O(2) and the integrated ones in three-momenta up to O. We find that some singular terms emerge in the integration over the time component of the four-momentum, which result in a new source term contributing to the chiral anomaly, in contrast to the well-known scenario of the Berry phase term. Finally we rewrite our results in any Lorentz frame with a reference four-velocity and show how the non-trivial transformation of the distribution function in different frames emerges in a natural way.

The Linear Boltzmann Transport (LBT) model coupled to hydrodynamical background is extended to include transport of both light partons and heavy quarks through the quark–gluon plasma (QGP) in high-energy heavy-ion collisions. The LBT model includes both elastic and inelastic medium-interaction of both primary jet shower partons and thermal recoil partons within perturbative QCD (pQCD). It is shown to simultaneously describe the experimental data on heavy and light flavor hadron suppression in high-energy heavy-ion collisions for different centralities at RHIC and LHC energies. More detailed investigations within the LBT model illustrate the importance of both initial parton spectra and the shapes of fragmentation functions on the difference between the nuclear modifications of light and heavy flavor hadrons. The dependence of the jet quenching parameter qˆ on medium temperature and jet flavor is quantitatively extracted.

We revisit the calculation of multiple parton scattering of a heavy quark in nuclei within the framework of recently improved high-twist factorization formalism, in which gauge invariance is ensured by a delicate setup of the initial partons' transverse momenta. We derive a new result for medium modified heavy quark fragmentation functions in deeply inelastic scattering. It is consistent with the previous calculation of light quark energy loss in the massless limit, but leads to a new correction term in the heavy quark case, which vanishes in the soft gluon radiation limit. We show numerically the significance of the new correction term in the calculation of heavy quark energy loss as compared to previous studies and with soft gluon radiation approximation.

We report on broadly based systematic investigations of the modeling components for open heavy-flavor diffusion and energy loss in strongly interacting matter in their application to heavy-flavor observables in high-energy heavy–ion collisions, conducted within an EMMI Rapid Reaction Task Force framework. Initial spectra including cold-nuclear-matter effects, a wide variety of space-time evolution models, heavy-flavor transport coefficients, and hadronization mechanisms are scrutinized in an effort to quantify pertinent uncertainties in the calculations of nuclear modification factors and elliptic flow of open heavy-flavor particles in nuclear collisions. We develop procedures for error assessments and criteria for common model components to improve quantitative estimates for the (low-momentum) heavy-flavor diffusion coefficient as a long-wavelength characteristic of QCD matter as a function of temperature, and for energy loss coefficients of high-momentum heavy-flavor particles.

We develop a general formalism for the quantum kinetics of chiral fermions in a background electromagnetic field based on a semiclassical expansion of covariant Wigner functions in the Planck constant. We demonstrate to any order of that only the time-component of the Wigner function is independent while other components are explicit derivative. We further demonstrate to any order of that a system of quantum kinetic equations for multiple-components of Wigner functions can be reduced to one chiral kinetic equation involving only the single-component distribution function. These are remarkable properties of the quantum kinetics of chiral fermions and will significantly simplify the description and simulation of chiral effects in heavy ion collisions and Dirac/Weyl semimetals. We present the unintegrated chiral kinetic equations in four-momenta up to O(2) and the integrated ones in three-momenta up to O. We find that some singular terms emerge in the integration over the time component of the four-momentum, which result in a new source term contributing to the chiral anomaly, in contrast to the well-known scenario of the Berry phase term. Finally we rewrite our results in any Lorentz frame with a reference four-velocity and show how the non-trivial transformation of the distribution function in different frames emerges in a natural way.

The Linear Boltzmann Transport (LBT) model coupled to hydrodynamical background is extended to include transport of both light partons and heavy quarks through the quark–gluon plasma (QGP) in high-energy heavy-ion collisions. The LBT model includes both elastic and inelastic medium-interaction of both primary jet shower partons and thermal recoil partons within perturbative QCD (pQCD). It is shown to simultaneously describe the experimental data on heavy and light flavor hadron suppression in high-energy heavy-ion collisions for different centralities at RHIC and LHC energies. More detailed investigations within the LBT model illustrate the importance of both initial parton spectra and the shapes of fragmentation functions on the difference between the nuclear modifications of light and heavy flavor hadrons. The dependence of the jet quenching parameter qˆ on medium temperature and jet flavor is quantitatively extracted.

We establish a linear Boltzmann transport (LBT) model coupled to hydrodynamical background to study hard parton evolution in heavy-ion collisions. Both elastic and inelastic scatterings are included in our calculations; and heavy and light flavor partons are treated on the same footing. Within this LBT model, we provide good descriptions of heavy and light hadron suppression and anisotropic flow in heavy-ion collisions. Angular correlation functions between heavy and light flavor hadrons are studied for the first time and shown able to quantify not only the amount of heavy quark energy loss, but also how the parton energy is re-distributed in parton showers.

We establish a linear Boltzmann transport (LBT) model coupled to hydrodynamical background to study hard parton evolution in heavy-ion collisions. Both elastic and inelastic scatterings are included in our calculations; and heavy and light flavor partons are treated on the same footing. Within this LBT model, we provide good descriptions of heavy and light hadron suppression and anisotropic flow in heavy-ion collisions. Angular correlation functions between heavy and light flavor hadrons are studied for the first time and shown able to quantify not only the amount of heavy quark energy loss, but also how the parton energy is re-distributed in parton showers.

Heavy-ion collisions at very high colliding energies are expected to produce a quark-gluon plasma (QGP) at the highest temperature obtainable in a laboratory setting. Experimental studies of these reactions can provide an unprecedented range of information on properties of the QGP at high temperatures. We report theoretical investigations of the physics perspectives of heavy-ion collisions at a future high-energy collider. These include initial parton production, collective expansion of the dense medium, jet quenching, heavy-quark transport, dissociation and regeneration of quarkonia, photon and dilepton production. We illustrate the potential of future experimental studies of the initial particle production and formation of QGP at the highest temperature to provide constraints on properties of strongly interaction matter.

Predictions made in Albacete et al. [Int. J. Mod. Phys. E 22 (2013) 1330007] prior to the LHC p+Pb run at sNN = 5 TeV are compared to currently available data. Some predictions shown here have been updated by including the same experimental cuts as the data. Some additional predictions are also presented, especially for quarkonia, that were provided to the experiments before the data were made public but were too late for the original publication.

A linearized Boltzmann transport (LBT) model coupled with hydrodynamical background is established to describe the evolution of jet shower partons and medium excitations in high energy heavy-ion collisions. We extend the LBT model to include both elastic and inelastic processes for light and heavy partons in the quark-gluon plasma. A hybrid model of fragmentation and coalescence is developed for the hadronization of heavy quarks. Within this framework, we investigate how heavy flavor observables depend on various ingredients, such as different energy loss and hadronization mechanisms, the momentum and temperature dependences of the transport coefficients, and the radial flow of the expanding fireball. Our model calculations show good descriptions of the D meson suppression and elliptic flow observed at the Larege Hadron Collider and the Relativistic Heavy-Ion Collider. The prediction for the Pb-Pb collisions at sNN=5.02 TeV is provided.

A sizable contribution to heavy-quark production in high-energy hadronic and nuclear collisions comes from heavy quark-antiquark pair production from gluon splitting during the parton shower evolution. We investigate the effect of gluon-medium interaction on open heavy flavor spectra in ultrarelativistic heavy-ion collisions. The interaction of hard gluons and heavy quarks with the hot QCD medium is simulated by utilizing a Langevin transport model that simultaneously incorporates contributions from collisional and radiative processes. It is found that, while the gluon splitting channel has quite an important contribution to the single D-meson production cross section, its influence on the final heavy-meson nuclear modification turns out to be quite modest because the average lifetime of hard gluons is short before splitting into heavy-quark-antiquark pairs during the evolution and propagation of the parton shower.

Heavy-ion collisions at very high colliding energies are expected to produce a quark-gluon plasma (QGP) at the highest temperature obtainable in a laboratory setting. Experimental studies of these reactions can provide an unprecedented range of information on properties of the QGP at high temperatures. We report theoretical investigations of the physics perspectives of heavy-ion collisions at a future high-energy collider. These include initial parton production, collective expansion of the dense medium, jet quenching, heavy-quark transport, dissociation and regeneration of quarkonia, photon and dilepton production. We illustrate the potential of future experimental studies of the initial particle production and formation of QGP at the highest temperature to provide constraints on properties of strongly interaction matter.

A linearized Boltzmann transport (LBT) model coupled with hydrodynamical background is established to describe the evolution of jet shower partons and medium excitations in high energy heavy-ion collisions. We extend the LBT model to include both elastic and inelastic processes for light and heavy partons in the quark-gluon plasma. A hybrid model of fragmentation and coalescence is developed for the hadronization of heavy quarks. Within this framework, we investigate how heavy flavor observables depend on various ingredients, such as different energy loss and hadronization mechanisms, the momentum and temperature dependences of the transport coefficients, and the radial flow of the expanding fireball. Our model calculations show good descriptions of the D meson suppression and elliptic flow observed at the Larege Hadron Collider and the Relativistic Heavy-Ion Collider. The prediction for the Pb-Pb collisions at sNN=5.02 TeV is provided.

A sizable contribution to heavy-quark production in high-energy hadronic and nuclear collisions comes from heavy quark-antiquark pair production from gluon splitting during the parton shower evolution. We investigate the effect of gluon-medium interaction on open heavy flavor spectra in ultrarelativistic heavy-ion collisions. The interaction of hard gluons and heavy quarks with the hot QCD medium is simulated by utilizing a Langevin transport model that simultaneously incorporates contributions from collisional and radiative processes. It is found that, while the gluon splitting channel has quite an important contribution to the single D-meson production cross section, its influence on the final heavy-meson nuclear modification turns out to be quite modest because the average lifetime of hard gluons is short before splitting into heavy-quark-antiquark pairs during the evolution and propagation of the parton shower.