Search

The general intuition that heavier partons suffer weaker energy loss inside a quark-gluon plasma (QGP) medium is critically re-examined. Within a linear Boltzmann transport model that includes both perturbative and non-perturbative interactions between heavy quarks and the QGP, we find that while the radiative energy loss is suppressed by the parton mass, heavier partons can experience stronger non-perturbative scatterings with the medium. Their competition may result in less energy loss of bottom quarks than charm quarks at low transverse momentum ($p_\mathrm{T}$) but an inverse order at high $p_\mathrm{T}$. Our model calculation provides a good description of the weaker nuclear modification on bottom particles than charm particles at low $p_\mathrm{T}$ observed at the RHIC experiment, and predicts an opposite hierarchy at high $p_\mathrm{T}$. A larger momentum space transport coefficient ($\hat{q}$) and a smaller spatial diffusion coefficient ($D_\mathrm{s}$) are found for bottom quarks than for charm quarks., Comment: 9 pages, 6 figures

In a quasi-particle model of QCD matter at finite temperature with thermal masses for quarks and gluons from hard thermal loops, the equation of state (EOS) can be described by an effective temperature dependence of the strong coupling $g(T)$. Assuming the same effective coupling between the exchanged gluon and thermal partons, the EOS can also be related to parton energy loss.} Based on the quasi-particle linear Boltzmann transport (QLBT) model coupled to a (3+1)-dimensional viscous hydrodynamic model of the quark-gluon plasma (QGP) evolution and a hybrid fragmentation-coalescence model for heavy quark hadronization, we perform a Bayesian analysis of the experimental data on $D$ meson suppression $R_{\rm AA}$ and anisotropy $v_2$ at RHIC and the LHC. We achieve a simultaneous constraint on the QGP EOS and the heavy quark transport coefficient, both consistent with the lattice QCD results., Comment: 8 pages, 5 figures

We investigate the hydrodynamic helicity polarization of $\Lambda$ hyperons, defined as the projection of the spin polarization vector along the directions of particle momenta, at RHIC-BES energies by utilizing the relativistic (3+1)D CLVisc hydrodynamics framework with SMASH initial conditions. As opposed to local spin polarization at high energy collisions, our hydrodynamic simulations demonstrate that the helicity polarization induced by the kinetic vorticity dominates over other contributions at intermediate and low collision energies. Our findings provide an opportunity to probe the fine structure of local kinetic vorticity as a function of azimuthal angle at intermediate and low collision energies by mapping our predictions to the future measurements in experiments., Comment: 8 pages, 3 figures

The quenching of light and heavy flavor hadrons in relativistic heavy-ion collisions probes the color and flavor dependences of parton energy loss through a color-deconfined quark-gluon plasma (QGP), and thus reveals the properties of QCD matter at extremely high density and temperature. By combining a next-to-leading order perturbative QCD calculation of parton production, a general ansatz of parton energy loss functions and parton fragmentation functions, we calculate the nuclear modification of various hadron species -- charged hadrons, $D$ mesons and $B$-decayed $J/\psi$ -- over a wide transverse momentum regime. Comparing our calculations to the experimental data using the Bayesian statistical analysis, we perform a first simultaneous extraction of the energy loss functions of gluons ($g$), light quarks ($q$), charm quarks ($c$) and bottom quarks ($b$) inside the QGP. We find that the average parton energy loss at high energies follows the expected hierarchy of $\langle \Delta E_g \rangle > \langle \Delta E_q \rangle \sim \langle \Delta E_c \rangle > \langle \Delta E_b \rangle$, while the parton energy loss distribution can further test the QCD calculations of parton interaction with the dense nuclear matter. We also find that the reduction of experimental uncertainties can significantly improve the precision of the extracted parton energy loss functions inside the QGP., Comment: 7 pages, 5 figures

We report our recent study on the global and local polarization of $\Lambda$ hyperons in Au+Au collisions at RHIC-BES energies within the (3+1)-dimensional CLVisc hydrodynamics framework. We present our numerical results for the global polarization as the function of collision energies and the local polarization along the beam direction as functions of azimuthal angle in $20-50$% centrality at $\sqrt{s_{NN}}$=7.7 GeV Au+Au collision energy. We have discussed the effects of initial conditions, Spin Hall effect and baryon diffusion., Comment: 4 pages, 3 figures. Contribution to the proceedings of Quark Matter 2023 (Houston, TX, 3-9 Sep. 2023)

The attractive interaction between $J/\psi$ and $\psi(3770)$ has to be strong enough to form the molecular resonant structure, i.e., $X(6900)$. We argue that since $\psi(3770)$ decays predominantly into a $D\bar D$ pair, the interactions between $J/\psi$ and $\psi(3770)$ may be significantly enhanced due to the three point $D\bar D$ loop diagram. The enhancement comes from the anomalous threshold located at $t=-1.288$GeV$^2$, which effect propagates into the $s$-channel partial wave amplitude in the vicinity of $\sqrt{s}\simeq 6.9$GeV. This effect may be helpful in forming the $X(6900)$ peak., Comment: 6 pages, 5 figures

The analyses of the LHCb data on $X(6900)$ found in di-$J/\psi$ and $J/\psi \psi(3686)$ systems are performed using a momentum-dependent {Flatt\'{e}}-like parameterization. The use of the pole counting rule and spectral density function sum rule give consistent conclusions that $X(6900)$ may not be a molecule of $J/\psi \psi(3686)$. Nevertheless it is still possible that $X(6900)$ be a molecule of higher states, such as $J/\psi\psi(3770)$, $\chi_{c0}\chi_{c2}$, etc., Comment: 7pages with 9 figures

The quenching of light and heavy flavor hadrons in relativistic heavy-ion collisions probes the color and flavor dependences of parton energy loss through a color-deconfined quark-gluon plasma (QGP), and thus reveals the properties of QCD matter at extremely high density and temperature. By combining a next-to-leading order perturbative QCD calculation of parton production, a general ansatz of parton energy loss functions and parton fragmentation functions, we calculate the nuclear modification of various hadron species -- charged hadrons, $D$ mesons and $B$-decayed $J/\psi$ -- over a wide transverse momentum regime. Comparing our calculations to the experimental data using the Bayesian statistical analysis, we perform a first simultaneous extraction of the energy loss functions of gluons ($g$), light quarks ($q$), charm quarks ($c$) and bottom quarks ($b$) inside the QGP. We find that the average parton energy loss at high energies follows the expected hierarchy of $\langle \Delta E_g \rangle > \langle \Delta E_q \rangle \sim \langle \Delta E_c \rangle > \langle \Delta E_b \rangle$, while the parton energy loss distribution can further test the QCD calculations of parton interaction with the dense nuclear matter. We also find that the reduction of experimental uncertainties can significantly improve the precision of the extracted parton energy loss functions inside the QGP., Comment: 7 pages, 5 figures

We perform a systematic study on the local and global spin polarization of $\Lambda$ and $\overline{\Lambda}$ hyperons in relativistic heavy-ion collisions at beam energy scan energies via the (3+1)-dimensional CLVisc hydrodynamics model with AMPT and SMASH initial conditions. Following the quantum kinetic theory, we decompose the polarization vector as the parts induced by thermal vorticity, shear tensor and the spin Hall effect (SHE). We find that the polarization induced by SHE and the total polarization strongly depends on the initial conditions. At $7.7$GeV, SHE gives a sizeable contribution and even flips the sign of the local polarization along the beam direction for AMPT initial condition, which is not observed for SMASH initial condition. Meanwhile, the local polarization along the out-of-plane direction induced by SHE with AMPT initial condition does not always increase with decreasing collision energies. Next, we find that the polarization along the beam direction is sensitive to the baryon diffusion coefficient, but the local polarization along the out-of-plane direction is not. Our results for the global polarization of $\Lambda$ and $\overline{\Lambda}$ agree well with the STAR data. Interestingly, the global polarization of $\overline{\Lambda}$ is not always larger than that of $\Lambda$ due to various competing effects. Our findings are helpful for understanding the polarization phenomenon and the detailed structure of quark-gluon plasma in relativistic heavy-ion collisions., Comment: 13 pages, 8 figures

Jet quenching has long been regarded as one of the key signatures for the formation of quark-gluon plasma in heavy-ion collisions. Despite significant efforts, the separate identification of quark and gluon jet quenching has remained as a challenge. Here we show that $J/\psi$ in high transverse momentum ($p_T$) region provides a uniquely sensitive probe of in-medium gluon energy loss since its production at high $p_T$ is particularly dominated by gluon fragmentation. Such gluon-dominance is first demonstrated for the baseline of proton-proton collisions within the framework of leading power NRQCD factorization formalism. We then use the linear Boltzmann transport model combined with hydrodynamics for the simulation of jet-medium interaction in nucleus-nucleus collisions. The satisfactory description of experimental data on both nuclear modification factor $R_{AA}$ and elliptic flow $v_2$ reveals, for the first time, that the gluon jet quenching is the driving force for high $p_T$ $J/\psi$ suppression. This novel finding is further confirmed, in a robust and model-independent way, by the data-driven Bayesian analyses of relevant experimental measurements, from which we also obtain the first quantitative extraction of the gluon energy loss distribution in the quark-gluon plasma., Comment: 6 pages, 5 figures

High-energy heavy-ion collisions at the Relativistic Heavy-Ion Collider and the Large Hadron Collider have evolved from qualitative understanding into precise extraction of the properties of the Quantum Chromodynamics medium at extremely high temperature. Jet quenching has offered unique insights into the transport properties of the Quark-Gluon Plasma (QGP) created in these energetic collisions. Apart from medium modification on jets, jet-induced medium excitation constitutes another crucial aspect of jet-QGP interaction, and is indispensable in understanding the soft components of jets. We review recent theoretical and phenomenological developments on medium response to jet energy loss, including an overview of both weakly and strongly coupled approaches for describing the thermalization and propagation of energy deposition from jets, effects of medium response on jet observables, and exploration of its unique signatures in jet-hadron correlations. Jet-induced medium excitation is shown to be an essential component in probing the in-medium dynamics of jets and a critical step towards precise extraction of the QGP properties., Comment: 25 pages, 9 figures, prepared for Ann. Rev. Nucl. Part. Sci

The interactions of quarks and gluons are strong at non-perturbative region. The equation of state (EoS) of a strongly-interacting quantum chromodynamics (QCD) medium can only be studied using the first-principle lattice QCD calculations. However, the complicated QCD EoS can be reproduced using simple statistical formula by treating the medium as a free parton gas whose fundamental degree of freedoms are dressed quarks and gluons called quasi-particles, with temperature-dependent masses. We use deep neural network and auto differentiation to solve this variational problem in which the masses of quasi gluons, up/down and strange quarks are three unknown functions, whose forms are represented by deep neural network. We reproduce the QCD EoS using these machine learned quasi-particle masses, and calculate the shear viscosity over entropy density ($\eta/s$) as a function of temperature of the hot QCD matter.

Jet quenching has long been regarded as one of the key signatures for the formation of quark-gluon plasma in heavy-ion collisions. Despite significant efforts, the separate identification of quark and gluon jet quenching has remained as a challenge. Here we show that $J/\psi$ in high transverse momentum ($p_T$) region provides a uniquely sensitive probe of in-medium gluon energy loss since its production at high $p_T$ is particularly dominated by gluon fragmentation. Such gluon-dominance is first demonstrated for the baseline of proton-proton collisions within the framework of leading power NRQCD factorization formalism. We then use the linear Boltzmann transport model combined with hydrodynamics for the simulation of jet-medium interaction in nucleus-nucleus collisions. The satisfactory description of experimental data on both nuclear modification factor $R_{AA}$ and elliptic flow $v_2$ reveals, for the first time, that the gluon jet quenching is the driving force for high $p_T$ $J/\psi$ suppression. This novel finding is further confirmed, in a robust and model-independent way, by the data-driven Bayesian analyses of relevant experimental measurements, from which we also obtain the first quantitative extraction of the gluon energy loss distribution in the quark-gluon plasma., Comment: 6 pages, 5 figures

While perturbative QCD is sufficient for understanding the color, mass and energy dependences of parton energy loss and jet quenching at large transverse momentum in heavy-ion collisions, a simultaneous description of heavy flavor nuclear modification factor $R_\mathrm{AA}$ and elliptic flow coefficient $v_2$ at low and intermediate $p_\mathrm{T}$ still remains a challenge due to the effects from non-perturbative interactions. In this work, we extend the linear Boltzmann transport model by implementing a generalized Cornell-type potential that incorporates both short-range Yukawa interaction and long-range color confining interaction between heavy quarks and the QGP medium. Combining our new approach for heavy-quark-QGP interaction with a (3+1)-dimensional hydrodynamic model CLVisc for the QGP evolution and a hybrid fragmentation-coalescence model for heavy quark hadronization, we obtain a satisfactory description of heavy meson $R_\mathrm{AA}$ and $v_2$ from low to intermediate to high $p_\mathrm{T}$ observed at both RHIC and the LHC. By model-data-comparison, we extract for the first time the in-medium heavy quark potential from open heavy flavor measurements; the result is in agreement with the lattice QCD calculation. Our study indicates that while jet quenching at high $p_\mathrm{T}$ is dominated by perturbative QCD interactions, non-perturbative interactions are indispensible for understanding heavy flavor quenching and flow at low and intermediate $p_\mathrm{T}$., Comment: 9 pages, 5 figures

We present a unique signal of jet-induced medium excitations: the enhancement of baryon-to-meson ratios around the quenched jets. To illustrate this, we study jet-particle correlations and the distributions of jet-induced identified particles with respect to the jet direction in Pb+Pb collisions at the LHC via a multi-phase transport model. We find a strong enhancement of baryon-to-meson ratios for associated particles at intermediate transverse momentum around the triggered jets in Pb+Pb collisions relative to p+p collisions, due to the coalescence of jet-excited medium partons. Since the lost energy from jets can diffuse to large angles, such baryon-to-meson-ratio enhancement is more pronounced for larger relative distance from the jet axis. We argue that the experimental confirmation of the enhancement of jet-induced baryon-to-meson ratios around the jets will provide an unambiguous evidence for the medium response to jet quenching in heavy-ion collisions., Comment: 6 pages, 3 figures

We explore the system size dependence of heavy-quark-QGP interaction by studying the heavy flavor meson suppression and elliptic flow in Pb-Pb, Xe-Xe, Ar-Ar and O-O collisions at the LHC. The space-time evolution of the QGP is simulated using a (3+1)-dimensional viscous hydrodynamic model, while the heavy-quark-QGP interaction is described by an improved Langevin approach that includes both collisional and radiative energy loss inside a thermal medium. Within this framework, we provides a reasonable description of the $D$ meson suppression and flow coefficients in Pb-Pb collisions, as well as predictions for both $D$ and $B$ meson observables in other collision systems yet to be measured. We find a clear hierarchy for the heavy meson suppression with respect to the size of the colliding nuclei, while their elliptic flow coefficient relies on both the system size and the geometric anisotropy of the QGP. Sizable suppression and flow are predicted for both $D$ and $B$ mesons in O-O collisions, which serve as a crucial bridge of jet quenching between large and small collision systems. Scaling behaviors between different collision systems are shown for heavy meson suppression factor and the bulk-eccentricity-rescaled heavy meson elliptic flow as functions of the number of participant nucleons in heavy-ion collisions., Comment: 10 pages, 9 figures

We perform the first study on asymmetric longitudinal decorrelations of elliptic, triangular and quadrangular flows in proton-nucleus collisions at the LHC and RHIC energies. To measure the longitudinal flow decorrelations for asymmetric collision systems, we propose a new set of rapidity-asymmetric flow decorrelation functions. Our event-by-event hydrodynamic calculations show that the flow decorrelations in proton-going direction are larger than those in nucleus-going direction. We also find that proton-nucleus collisions at RHIC have larger longitudinal flow decorrelation effects than those at the LHC. Our study opens a new window to probe the longitudinal properties and the origin of flows in relativistic nuclear collisions., Comment: 6 pages, 3 figures

Jet-medium interaction involves two important effects: jet energy loss and medium response. The search for jet-induced medium excitations is one of the hot topics in jet quenching study in relativistic nuclear collisions. In this work, we perform a systematic study on how the lost energy from hard jets evolves with the bulk medium and redistributes in the final state of heavy-ion collisions via a multi-phase transport model. In particular, the ($\Delta \eta, \Delta \phi$) distribution of charged particles with respect to the jet axis and jet shape function are studied for various Pb+Pb collision centralities and for different transverse momentum intervals of charged particles. Our numerical result shows a strong enhancement of soft particles at large angles for Pb+Pb collisions relative to p+p collisions at the LHC, qualitatively consistent with recent CMS data. This indicates that a significant fraction of the lost energy from hard jets is carried by soft particles at large angles away from the jet axis., Comment: 10 pages, 9 figures

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 $R_\mathrm{AA}$ and the elliptic flow $v_2$ 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 $R_\mathrm{AA}$ and $v_2$, we extract the heavy quark transport parameter $\hat{q}$ and diffusion coefficient $D_\mathrm{s}$ in the temperature range of $1-4~T_\mathrm{c}$, and compare them with the lattice QCD results and other phenomenological studies., Comment: 13 pages, 9 figures

Recent experiments have observed large anisotropic collective flows in high multiplicity proton-lead collisions at the Large Hadron Collider (LHC), which indicates the possible formation of mini quark-gluon plasma (QGP) in small collision systems. However, no jet quenching has been confirmed in such small systems so far. To understand this intriguing result, the system size scan experiments have been proposed to bridge the gap between large and small systems. In this work, we perform a systematic study on both heavy and light flavor jet quenching in different collision systems at the LHC energies. Using our state-of-the-art jet quenching model, which combines the next-to-leading-order perturbative QCD framework, a linear Boltzmann transport model and the (3+1)-dimensional viscous hydrodynamics simulation, we provide a good description of nuclear modification factor $R_{\rm AA}$ for charged hadrons and $D$ mesons in central and mid-central Pb+Pb and Xe+Xe collisions measured by CMS collaboration. We further predict the transverse momentum and centrality dependences of $R_{AA}$ for charged hadrons, $D$ and $B$ mesons in Pb+Pb, Xe+Xe, Ar+Ar and O+O collisions at the LHC energies. Our numerical results show a clear system size dependence for both light and heavy flavor hadron $R_{AA}$ across different collision systems. Sizable jet quenching effect is obtained for both heavy and light flavor hadrons in central O+O collisions at the LHC energies. Our study provides a significant bridge for jet quenching from large to small systems, and should be helpful for finding the smallest QGP droplet and the disappearance of QGP in relativistic nuclear collisions., Comment: 12 pages, 8 figures

To study the bulk properties of the quark-gluon-plasma (QGP) produced at the 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 equation for baryon diffusion with the NEOS-BQS equation of state, fluctuating initial conditions from 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 increases mildly with the collision energy due to larger radial flow. We further compute the multiple-particle cumulant ratio $v_2\{4\}/v_2\{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., Comment: 10 pages, 6 figures

We briefly summarize our recent study on the number of constituent quark (NCQ) scaling of hadron elliptic flows in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}=$ 5.02 TeV. With the inclusion of hadron production via the quark coalescence model at intermediate $p_T$, the viscous hydrodynamics at low $p_T$, and jet fragmentation at high $p_T$, our $Hydro-Coal-Frag$ model provides a nice description of the $p_T$-spectra and differential elliptic flow $v_2(p_T)$ of pions, kaons and protons over the $p_T$ range from 0 to 6 GeV. Our results demonstrate that including the quark coalescence is essential for reproducing the observed approximate NCQ scaling of hadron $v_2$ at intermediate $p_T$ in experiments, indicating strongly the existence of partonic degrees of freedom and the formation of quark-gluon plasma in high multiplicity p--Pb collisions at the LHC., Comment: 4 pages, 2 figures, QM2019 proceedings

Using a model based on the Color Glass Condensate framework and the dilute-dense factorization, we systematically study the azimuthal angular correlations between a heavy flavor meson and a light reference particle in proton-nucleus collisions. The obtained second harmonic coefficients (also known as the elliptic flows) for $J/\psi$ and $D^0$ agree with recent experimental data from the LHC. We also provide predictions for the elliptic flows of $\Upsilon$ and $B$ meson, which can be measured in the near future at the LHC. This work can shed light on the physics origin of the collectivity phenomenon in the collisions of small systems., Comment: 14 pages, 7 figures. This work is an extension of our earlier publication 10.1103/PhysRevLett.122.172302 (arXiv:1901.10320); v2 with minor updates

Within an advanced Langevin-hydrodynamics framework coupled to a hybrid fragmentation-coalescence hadronization model, we study heavy flavor quenching and flow in relativistic heavy-ion collisions. We investigate how the initial heavy quark spectrum, the energy loss and hadronization mechanisms of heavy quarks in medium, the evolution profile of pre-equilibrium stage, the flow of medium and the temperature dependence of heavy quark diffusion coefficient influence the suppression and elliptic flow of heavy mesons at RHIC and the LHC. Our result shows that different modeling of initial conditions, pre-equilibrium evolution and in-medium interaction can individually yield about 10-40% uncertainties in D meson suppression and flow at low transverse momentum. We also find that a proper combination of collisional versus radiative energy loss, coalescence versus fragmentation in hadronization, and the inclusion of medium flow are the most important factors for describing the suppression and elliptic flow of heavy mesons., Comment: 9 pages, 7 figures; to be published in Chinese Physics C

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 $R_\text{AA}$'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., Comment: 4 pages, 4 figures, contribution to the proceedings of Quark Matter - XXVIII international conference on ultrarelativistic heavy-ion collisions (Wuhan, China, 4-9 November, 2019)

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 (GLV) 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., Comment: 33 pages in RevTeX with 30 figures, final version appeared in PRD with additional typos corrected

Heavy flavor probes are sensitive to the properties of the quark gluon plasma (QGP) produced in relativistic heavy-ion collisions. A huge amount of effort has been devoted to studying different aspects of the heavy-ion collisions using heavy flavor particles. In this work, we study the dynamics of heavy quark transport in the QGP medium using the rapidity dependence of heavy flavor observables. We calculate the nuclear modification of $\text{B}$ and $\text{D}$ meson spectra as well as spectra of leptons from heavy flavor decays in the rapidity range $[-4.0,4.0]$. We use an implementation of the improved Langevin equation with gluon radiation on top of a (3+1)-dimensional relativistic viscous hydrodynamical background for several collision setups. We find that the rapidity dependence of the heavy quark modification is determined by the interplay between the smaller size of the medium, which affects the path length of the heavy quarks, and the softer heavy quark initial production spectrum. We compare our results with available experimental data and present predictions for open heavy flavor meson $R_\text{AA}$ at finite rapidity., Comment: 7 pages, 5 figures

We have developed a numerical framework for a full solution of the relativistic Boltzmann equations for the quark-gluon matter using the multiple Graphics Processing Units (GPUs) on distributed clusters. Including all the $2 \to 2$ scattering processes of 3-flavor quarks and gluons, we compute the time evolution of distribution functions in both coordinate and momentum spaces for the cases of pure gluons, quarks and the mixture of quarks and gluons. By introducing a symmetrical sampling method on GPUs which ensures the particle number conservation, our framework is able to perform the space-time evolution of quark-gluon system towards thermal equilibrium with high performance. We also observe that the gluons naturally accumulate in the soft region at the early time, which may indicate the gluon condensation., Comment: 32 pages, 12 figures

We have developed a numerical framework for a full solution of the relativistic Boltzmann equations for the quark-gluon matter using the multiple Graphics Processing Units (GPUs) on distributed clusters. Including all the $2 \to 2$ scattering processes of 3-flavor quarks and gluons, we compute the time evolution of distribution functions in both coordinate and momentum spaces for the cases of pure gluons, quarks and the mixture of quarks and gluons. By introducing a symmetrical sampling method on GPUs which ensures the particle number conservation, our framework is able to perform the space-time evolution of quark-gluon system towards thermal equilibrium with high performance. We also observe that the gluons naturally accumulate in the soft region at the early time, which may indicate the gluon condensation., Comment: 32 pages, 12 figures

We investigate the role of partonic degrees of freedom in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}=$ 5.02 TeV carried out at the Large Hadron Collider (LHC) by studying the production and collective flow of identified hadrons at intermediate $p_T$ via the coalescence of soft and hard partons generated from viscous hydrodynamics (VISH2+1) and the energy loss model (LBT), respectively. We find that combining the intermediate $p_T$ hadrons from the coalescence with the low $p_T$ hadrons from hydrodynamics and high $p_T$ hadrons from the jet fragmentation, our Hydro-Coal-Frag model provides a nice description of the measured $p_T$-spectra and the differential elliptic flow $v_2(p_T)$ of pions, kaons and protons over the $p_T$ range from 0 to 6 GeV. We further demonstrate the necessity of including the quark coalescence contribution to reproduce the experimentally observed approximate number of constituent quark scaling of hadron $v_2$ at intermediate $p_T$. Our results thus indicate the importance of partonic degrees of freedom and also hint at the possible formation of quark-gluon plasma in high multiplicity p+Pb collisions at the LHC., Comment: 6 pages, 3 figures, published version

We develop for charmed hadron production in relativistic heavy-ion collisions a comprehensive coalescence model that includes an extensive set of $s$ and $p$-wave hadronic states as well as the strict energy-momentum conservation, which ensures the boost invariance of the coalescence probability and the thermal limit of the produced hadron spectrum. By combining our hadronization scheme with an advanced Langevin-hydrodynamics model that incorporates both elastic and inelastic energy loss of heavy quarks inside the dynamical quark-gluon plasma, we obtain a successful description of the $p_\mathrm{T}$-integrated and differential $\Lambda_c/D^0$ and $D_s/D^0$ ratios measured at RHIC and the LHC. We find that including the effect of radial flow of the medium is essential for describing the enhanced $\Lambda_c/D^0$ ratio observed in relativistic heavy-ion collisions. We also find that the puzzling larger $\Lambda_c/D^0$ ratio observed in Au+Au collisions at RHIC than in Pb+Pb collisions at the LHC is due to the interplay between the effects of the QGP radial flow and the charm quark transverse momentum spectrum at hadronization. Our study further suggests that charmed hadrons have larger sizes in medium than in vacuum., Comment: 6 pages, 5 figures

With our coupled jet-fluid model, we study the nuclear modifications of full jets and jet structures for single inclusive jets and $\gamma$-jets in Pb+Pb collisions at $5.02$~ATeV and $2.76$~ATeV. The in-medium evolution of full jet shower is described by a set of coupled transport equations including the effects of collisional energy loss, transverse momentum broadening and medium-induced splitting process. The dynamical evolution of bulk medium is simulated by solving relativistic hydrodynamic equation with source term which accounts for the energy and momentum deposited by hard jet shower to soft medium. Our study demonstrates that the hydrodynamic medium response to jet propagation significantly enhances the broadening of jet shape at large angles and is essential for the cone-size dependence of jet energy loss and nuclear modification factor of inclusive jet production. It is also found that the nuclear modification pattern of jet shape is sensitive to jet energy but has weak dependence on the flavor of the parton that initiates the jet. Our result can naturally explain the different nuclear modification patterns of jet shape functions for single inclusive jet and $\gamma$-jet events as observed by the CMS Collaboration, and can be tested in the future by measuring the jet shape function over a wider range of jet energies in heavy-ion collisions., Comment: 10 pages, 8 figures. To be published in PLB

Relativistic heavy-ion experiments have observed similar quenching effects for (prompt) $D$ mesons compared to charged hadrons for transverse momenta larger than 6-8~GeV, which remains a mystery since heavy quarks typically lose less energies in quark-gluon plasma than light quarks and gluons. Recent measurements of the nuclear modification factors of $B$ mesons and $B$-decayed $D$ mesons by the CMS Collaboration provide a unique opportunity to study the flavor hierarchy of jet quenching. Using a linear Boltzmann transport model combined with hydrodynamics simulation, we study the energy loss and nuclear modification for heavy and light flavor jets in high-energy nuclear collisions. By consistently taking into account both quark and gluon contributions to light and heavy flavor hadron productions within a next-to-leading order perturbative QCD framework, we obtain, for the first time, a satisfactory description of the experimental data on the nuclear modification factors for charged hadrons, $D$ mesons, $B$ mesons and $B$-decayed $D$ mesons simultaneously over a wide range of transverse momenta (8-300~GeV). This presents a solid solution to the flavor puzzle of jet quenching and constitutes a significant step towards the precision study of jet-medium interaction. Our study predicts that at transverse momenta larger than 30-40~GeV, $B$ mesons also exhibit similar suppression effects to charged hadrons and $D$ mesons, which may be tested by future measurements., Comment: 6 pages, 5 figures

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 (GLV) 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., Comment: 33 pages in RevTeX with 30 figures, final version appeared in PRD with additional typos corrected

Heavy flavor probes are sensitive to the properties of the quark gluon plasma (QGP) produced in relativistic heavy-ion collisions. A huge amount of effort has been devoted to studying different aspects of the heavy-ion collisions using heavy flavor particles. In this work, we study the dynamics of heavy quark transport in the QGP medium using the rapidity dependence of heavy flavor observables. We calculate the nuclear modification of $\text{B}$ and $\text{D}$ meson spectra as well as spectra of leptons from heavy flavor decays in the rapidity range $[-4.0,4.0]$. We use an implementation of the improved Langevin equation with gluon radiation on top of a (3+1)-dimensional relativistic viscous hydrodynamical background for several collision setups. We find that the rapidity dependence of the heavy quark modification is determined by the interplay between the smaller size of the medium, which affects the path length of the heavy quarks, and the softer heavy quark initial production spectrum. We compare our results with available experimental data and present predictions for open heavy flavor meson $R_\text{AA}$ at finite rapidity., Comment: 7 pages, 5 figures

Heavy quarks serve as ideal probes of the QGP properties produced in energetic nuclear collisions, and provide a unique opportunity to study the mass effects on parton energy loss. We develop a multi-stage approach for heavy quark evolution inside the QGP, in which heavy quarks first undergo a rare-scattering multiple-emission evolution at momenta large compared to their mass (sensitive only to the transverse diffusion coefficient $\hat{q}$ ), and then evolve through a single-scattering induced emission (Gunion-Bertsch) stage at momenta comparable to their mass [sensitive to not only $\hat{q}$, but also the longitudinal drag $\hat{e}$ and diffusion $\hat{e}_2$ coefficients]. This multi-stage approach is coupled to a (2+1)-D viscous hydrodynamic model for a quantitative investigation of charm vs. beauty quark energy loss inside the QGP. Based on this approach, we find that drag induced radiation has a considerable impact on the energy loss of intermediate $p_T$ massive beauty quarks. This effect increases the suppression of B mesons and narrows the difference between the $R_{AA}$ of B and D mesons. Our results are consistent with the experimental data at the LHC and contribute to a more quantitative understanding of the transverse momentum dependence of the mass hierarchy of parton energy loss inside the QGP., Comment: 5 pages, 2 figures, Latex, Proceedings of the International Conference on Hard and Electromagnetic Probes of High-Energy Nuclear Collisions, 30 September - 5 October 2018, Aix-Les-Bains, Savoie, France

Using the dilute-dense factorization in the Color Glass Condensate framework, we investigate the azimuthal angular correlation between a heavy quarkonium and a charged light hadron in proton-nucleus collisions. We extract the second harmonic $v_2$, commonly known as the elliptic flow, with the light hadron as the reference. This particular azimuthal angular correlation between a heavy meson and a light hadron has first been measured at the LHC recently. The experimental results indicate that the elliptic flows for heavy-flavor mesons ($J/\psi$ and $D^0$) are almost as large as those for light hadrons. Our calculation demonstrates that this result can be naturally interpreted as an initial state effect due to the interaction between the incoming partons from the proton and the dense gluons inside the target nucleus. Since the heavy quarkonium $v_2$ exhibits a weak mass dependence according to our calculation, we predict that the heavy quarkonium $\Upsilon$ should have a similar elliptic flow as compared to that of the $J/\psi$, which can be tested in future measurements., Comment: 6 pages, 3 figures;v2: a few minor typo corrected

We study the suppressions of high transverse momentum single hadron and dihadron productions in high-energy heavy-ion collisions based on the framework of a next-to-leading-order perturbative QCD parton model combined with the higher-twist energy loss formalism.Our model can provide a consistant description for the nuclear modification factors of single hadron and dihadron productions in central and non-central nucleus-nucleus collisions at RHIC and the LHC energies. We quantitatively extract the value of jet quenching parameter $\hat q$ via a global $\chi^2$ analysis, and obtain ${\hat{q}}/{T^3} = 4.1 \sim 4.4$ at $T = 378$~MeV at RHIC and ${\hat{q}}/{T^3} = 2.6 \sim 3.3$ at $T = 486$~MeV at the LHC, which are consistent with the results from JET Collaboration. We also provide the predictions for the nuclear modification factors of dihadron productions in Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV and in Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}$ = 5.44 TeV.

Heavy flavor probes are sensitive to the properties of the quark gluon plasma (QGP) produced in relativistic heavy-ion collisions. A huge amount of effort has been devoted to studying different aspects of the heavy-ion collisions using heavy flavor particles. In this work, we study the dynamics of heavy quark transport in the QGP medium using the rapidity dependence of heavy flavor observables. We calculate the nuclear modification of $\text{B}$ and $\text{D}$ meson spectra as well as spectra of leptons from heavy flavor decays in the rapidity range $[-4.0,4.0]$. We use an implementation of the improved Langevin equation with gluon radiation on top of a (3+1)-dimensional relativistic viscous hydrodynamical background for several collision setups. We find that the rapidity dependence of the heavy quark modification is determined by the interplay between the smaller size of the medium, which affects the path length of the heavy quarks, and the softer heavy quark initial production spectrum. We compare our results with available experimental data and present predictions for open heavy flavor meson $R_\text{AA}$ at finite rapidity., Comment: 7 pages, 5 figures

We study the medium-induced gluon emission process experienced by a hard jet parton propagating through the dense nuclear matter in the framework of deep inelastic scattering off a large nucleus. We work beyond the collinear rescattering expansion and the soft gluon emission limit, and derive a closed formula for the medium-induced single gluon emission spectrum from a heavy or light quark jet interacting with the dense nuclear medium via transverse and longitudinal scatterings. Without performing the collinear rescattering expansion, the medium-induced gluon emission spectrum is controlled by the full distribution of the differential elastic scattering rates between the propagating partons and the medium constituents. We further show that if one utilizes heavy static scattering centers for the traversed nuclear matter and takes the soft gluon emission limit, our result can reduce to the first order in opacity Djordjevic-Gyulassy-Levai-Vitev formula., Comment: 21 pages, 14 figures

We perform a systematic study on the decorrelation of anisotropic flows along the pseudorapidity in relativistic heavy-ion collisions at the LHC and RHIC energies. The dynamical evolution of the QGP fireball is simulated via the CLVisc (ideal) (3+1)-dimensional hydrodynamics model, with the fully fluctuating initial condition from A-Multi-Phase-Transport (AMPT) model. Detailed analysis is performed on the longitudinal decorrelations of elliptic, triangular and quadrangular flows in terms of flow vectors, flow magnitudes and flow orientations (event planes). It is found that pure flow magnitudes have smaller longitudinal decorrelation than pure flow orientations, and the decorrelation of flow vectors is a combined effect of both flow magnitudes and orientations. The longitudinal decorrelation of elliptic flow has a strong and non-monotonic centrality dependence due to the initial elliptic collision geometry: smallest decorrelation in mid-central collisions. In contrast, the decorrelations of triangular and quadrangular flows have weak centrality dependence, slightly larger decorrelations in more peripheral collisions. Our numerical results for Pb+Pb collisions at the LHC are in good agreement with the ATLAS data, while our RHIC results predict much larger longitudinal decorrelations as compared to the LHC. We further analyze the longitudinal structures of the AMPT initial conditions and find that the final-state longitudinal decorrelation effects are strongly correlated with the lengths of the initial string structures in the AMPT model. The decorrelation effects are typically larger at lower collision energies and in more peripheral collisions due to shorter lengths of the string structures in the initial states., Comment: 14 pages, 14 figures

We study the medium-induced gluon emission from a hard quark jet traversing the dense nuclear matter within the framework of deep inelastic scattering off a large nucleus. We extend the previous work and compute the single gluon emission spectrum including both transverse and longitudinal momentum exchanges between the hard jet parton and the medium constituents. On the other hand, with only transverse scattering and using static scattering centers for the traversed medium, our induced gluon emission spectrum in the soft gluon limit reduces to the Gyulassy-Levai-Vitev one-rescattering-one-emission formula., Comment: 16 pages, 13 figures

In this paper, we study the production of isolated-photon plus a jet in $pp$ and $PbPb$ collisions, which can be used as an important probe to the jet transport property in quark gluon plasma created in heavy ion collisions. Normally, there are two types of observables associated with the production of isolated-photon plus a jet, namely, the azimuthal angular correlation and the transverse momentum imbalance. To understand both observables in the full kinematical region, we need to employ the perturbative QCD calculation, which takes into account the hard splitting of partons, together with the Sudakov resummation formalism, which resums soft gluon splittings. Furthermore, by introducing energy-loss into the system, we calculate the enhancement of the momentum imbalance distribution for $AA$ as compared to $pp$ collisions and make predictions for future unfolded experimental data. In addition, in order to extract the jet transport coefficient more precisely in our numerical calculation, we also distinguish quark jets from gluon jets, since they interact with quark gluon plasma with different strengths. This work provides a reliable theoretical tool for the calculation of the gamma-jet correlation, which can lead us to a more precise extraction of the jet transport coefficient in relativistic heavy-ion collisions., Comment: 11 pages, 6 figures

It is argued that heavy-quarks traversing a Quark Gluon Plasma, undergo a multi-stage modification process, similar to the case of light flavors. Such an approach is applied to heavy-quark energy loss, which includes a rare-scattering, multiple emission formalism at momenta large compared to the mass (sensitive only to the transverse diffusion coefficient $\hat{q}$), and a single scattering induced emission formalism (Gunion-Bertsch) at momenta comparable to the mass of the quark [sensitive to $\hat{q}$, and the longitudinal drag ($\hat{e}$) and diffusion ($\hat{e}_2$) coefficients]. The application of such a multi-stage approach within a 2+1D viscous fluid-dynamical simulation leads to simultaneous agreement with experimental data for the nuclear modification factor of both $B$ and $D$ mesons, as measured by the CMS collaboration at the LHC. The extracted transport coefficients are found to be consistent with those for light flavors., Comment: 6 pages, 3 figures

The nuclear modification of jet splitting in relativistic heavy-ion collisions at RHIC and the LHC energies is studied based on the higher twist formalism. Assuming coherent energy loss for the two splitted subjets, a non-monotonic jet energy dependence is found for the nuclear modification of jet splitting function: strongest modification at intermediate jet energies whereas weaker modification for larger or smaller jet energies. Combined with the smaller size and lower density of the QGP medium at RHIC than at the LHC, this explains the CMS-STAR groomed jet puzzle -- strong nuclear modification of the momentum sharing $z_g$ distribution at the LHC whereas no obvious modification of the $z_g$ distribution at RHIC. In contrast, the observed nuclear modification pattern of the groomed jet $z_g$ distribution cannot be explained solely by independent energy loss of the two subjets. Our result may be tested in future measurements of groomed jets with lower jet energies at the LHC and larger jet energies at RHIC, for different angular separations between the two subjets., Comment: 10 pages, 12 figures

We present here a possible resolution of the jet quenching parameter $\hat{q}$ puzzle by exploring the momentum fraction $x$ and scale $Q^2$ dependence of $\hat{q}$, in addition to temperature dependence. We explore the momentum broadening of the jet due to Glauber gluon exchange with the quark-gluon plasma (QGP). This has led us to define the momentum fraction $x$ for the partons in the QGP. We also show here, for the first time, possible forms of the parton distribution function (PDF) of the QGP. The QGP-PDF as input in $\hat{q}$ has been used to fit the PHENIX and CMS data simultaneously. We observe that the scale evolution of the QGP-PDF and the energy of the jet are the missing ingredients responsible for the enhancement of $\hat{q}$ at the same temperature in RHIC collisions compared to LHC collisions., Comment: 4 pages, 8 figures, talk at conference Quark-Matter 2017

We present the first attempt to extract the "$x$" dependence of the parton distribution function (PDF) of the quark gluon plasma (QGP). In the absence of knowledge regarding the mass of a QGP constituent, we define the new variable $x_N$, the nucleonic momentum fraction, which represents the ratio of the momentum of the parton to that of a self-contained section of the plasma that has the mass of a nucleon. Calculations are compared to data for single hadron suppression in terms of the nuclear modification factor $R_{AA}$ and the azimuthal anisotropy parameter $v_{2}$, as a function of transverse momentum $p_{T}$, centrality and energy of the collision. It is demonstrated that the scale evolution of the QGP-PDF is responsible for the reduction in normalization of the jet transport coefficient $\hat{q}$, between fits to Relativistic Heavy-Ion Collider (RHIC) and Large Hadron Collider (LHC) data; a puzzle, first discovered by the JET collaboration. Best fits to data are suggestive of the presence of quasi-particles with wide dispersion relations in the QGP., Comment: 5 pages, 4 figures, RevTex4.1

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 $\hat{q}$ on medium temperature and jet flavor is quantitatively extracted., Comment: 6 pages, 6 figures

We study the nuclear modifications of full jets and their structures in relativistic heavy-ion collisions including the effect of hydrodynamic medium response to jet quenching. To study the evolutions of the full jet shower and the traversed medium with energy and momentum exchanges between them, we formulate a coupled jet-fluid model consisting of a set of jet transport equations and relativistic hydrodynamics equations with source terms. In our model, the full jet shower interacts with the medium and gets modified via collisional and radiative processes during the propagation. Meanwhile, the energy and momentum are deposited from the jet shower to the medium and then evolve with the medium hydrodynamically. The full jet defined by a cone size in the final state includes the jet shower and the particles produced from jet-induced flow. We apply our model to calculate the full jet energy loss and the nuclear modifications of jet rate and shape in Pb+Pb collisions at $2.76{\rm A~TeV}$. It is found that the inclusion of jet-induced flow contribution leads to stronger jet-cone size dependence for jet energy loss and jet suppression. Jet-induced flow also has a significant contribution to jet shape function and dominates at large angles away from the jet axis., Comment: 13 pages, 6 figures

We study the nuclear modifications of full jets and their structures in relativistic heavy-ion collisions including the effect of hydrodynamic medium response to jet quenching. To study the evolutions of the full jet shower and the traversed medium with energy and momentum exchanges between them, we formulate a coupled jet-fluid model consisting of a set of jet transport equations and relativistic hydrodynamics equations with source terms. In our model, the full jet shower interacts with the medium and gets modified via collisional and radiative processes during the propagation. Meanwhile, the energy and momentum are deposited from the jet shower to the medium and then evolve with the medium hydrodynamically. The full jet defined by a cone size in the final state includes the jet shower and the particles produced from jet-induced flow. We apply our model to calculate the full jet energy loss and the nuclear modifications of jet rate and shape in Pb+Pb collisions at $2.76{\rm A~TeV}$. It is found that the inclusion of jet-induced flow contribution leads to stronger jet-cone size dependence for jet energy loss and jet suppression. Jet-induced flow also has a significant contribution to jet shape function and dominates at large angles away from the jet axis., Comment: 13 pages, 6 figures