Source author record

Taesoo Song

Taesoo Song appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

In-medium effects in $ϕ$ meson production in heavy-ion collisions from subthreshold to relativistic energies

We investigate the hidden strange $ϕ$ meson production in heavy-ion collisions from subthreshold ($E_{kin}\approx 1$ A GeV) to relativistic ($E_{kin}\approx 21$ A TeV) energies as well as its coupling to the open strange mesons (kaons, antikaons) and their productions. Our study is based on the off-shell microscopic Parton-Hadron-String Dynamics (PHSD) transport approach which is applicable for the dynamical description of strongly interacting hadronic and partonic degrees-of-freedom created in heavy-ion collisions. Implementing novel meson-baryon and meson-hyperon production channels for $ϕ$ mesons, calculated within a T-matrix coupled channel approach based on the extended SU(6) chiral effective Lagrangian model, along with the collisional broadening of the $ϕ$-meson in-medium spectral function, we find a substantial enhancement of $ϕ$ meson production in heavy-ion collisions, especially at sub- and near-thresholds. This allows to describe the experimentally observed strong enhancement of the $ϕ/K^-$ ratio at low energies without including hypothetical decays of heavy baryonic resonances to $ϕ$ as in alternative approaches. Moreover, we show that in spite of a stronger contribution from enhanced $ϕ$ to $K^-$ production, the majority of the experimental data for different A+A systems at low energies favor the scenario with in-medium modifications of the kaon and antikaon properties in the hot and dense environment. Moreover, we study the influence of the final state interactions of $K, \bar K$ mesons on the reconstruction of $ϕ$'s by the by invariant mass method.

preprint2022arXiv

Transport Model Comparison Studies of Intermediate-Energy Heavy-Ion Collisions

Transport models are the main method to obtain physics information from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. Calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, and more realistic calculations of heavy-ion collisions. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. We review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known as BUU and QMD type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. We often indicate improved strategies in performing transport simulations and thus provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

preprint2021arXiv

Comparison of heavy ion transport simulations: Ag+Ag collisions at Elab = 1.58 AGeV

We compare the microscopic transport models UrQMD, PHSD, PHQMD, and SMASH to make predictions for the upcoming Ag+Ag data at $E_\mathrm{lab}=1.58$~$A$GeV ($\sqrt{s_\mathrm{NN}}=2.55$~GeV) by the HADES collaboration. We study multiplicities, spectra and effective source temperatures of protons, $π^{\pm,0}$, $K^\pm$, the $η$, $Λ+Σ^0$ and the $Ξ^-$ within these models. Despite variations in the detailed implementation of the dynamics in the different models, the employed transport approaches all show consistent multiplicities of the bulk of investigated hadrons. The main differences are in the $Ξ^-$ production, which is treated differently between UrQMD/SMASH on one side employing high mass resonance states with explicit decays to $\mathrm{Resonance}\rightarrow Ξ+K+K$ in contrast to PHSD/PHQMD which account only non-resonant $Ξ$ production channels. A comparison of the spectra, summarized by effective source temperatures, shows that all models provide similar source temperatures around $T_\mathrm{source}=80-95$~MeV, and show substantial radial flow on the order of $\langle v_T\rangle=0.22c-0.3c$ even for such a small system.

preprint2020arXiv

$J/ψ$ near $T_c$

We calculate the mass shift and thermal decay width of the $J/ψ$ near the QCD transition temperature $T_c$ by imposing two independent constraints on these variables that can be obtained first by solving the Schrödinger equation and second from the QCD sum rule approach. While the real part of the potential is determined by comparing the QCD sum rule result for charmonium and the D meson to that from the potential model result, the imaginary potential is taken to be proportional to the perturbative form multiplied by a constant factor, which in turn can be determined by applying the two independent constraints. The result shows that the binding energy and the thermal width becomes similar in magnitude at around $T=1.09T_c$, above which the sum rule analysis also becomes unstable, strongly suggesting that the $J/ψ$ will melt slightly above $T_c$.

preprint2020arXiv

D meson mass and heavy quark potential at finite temperature

Based on the observation that the heavy quark-antiquark potential value at infinity corresponds to twice the D meson mass, we constrain the asymptotic value of the heavy quark potential in a hot medium through a QCD sum rule calculation of the D meson at finite temperature. We find that to correctly reproduce the QCD sum rule results as well as a recent model calculation for the D meson mass near the critical temperature, the heavy quark potential should be composed mostly of the free energy with an addition of a small but nontrivial fraction of the internal energy. Combined with a previous study comparing potential model results for the $J/ψ$ to a QCD sum rule calculation, we conclude that the composition of the effective heavy quark potential should depend on the inter-quark distance. Namely, the potential is dominated by the free energy at short distance, while at larger separation, it has a fraction of about 20% of internal energy.

preprint2020arXiv

Exploring non-equilibrium quark-gluon plasma effects on charm transport coefficients

In this article we investigate how the drag coefficient $A$ and $\hat{q}$, the transverse momentum transfer by unit length, of charm quarks are modified if the QGP is not in complete thermal equilibrium using the dynamical quasi-particle model (DQPM) which reproduces both, the equation-of-state of the QGP and the spatial diffusion coefficient of heavy quarks as predicted by lattice QCD calculations. We study three cases: a) the QGP has an anisotropic momentum distribution of the partons which leads to an anisotropic pressure b) the QGP partons have higher or lower kinetic energies as compared to the thermal expectation value, and c) the QGP partons have larger or smaller pole masses of their spectral function as compared to the pole mass from the DQPM at the QGP temperature. In the last two cases we adjust the number density of partons to obtain the same energy density as in an equilibrated QGP. In the first scenario we find that if the transverse pressure exceeds the longitudinal one for small heavy quark momenta $A$ becomes larger and $\hat{q}$ smaller as compared to an isotropic pressure. For heavy quarks with large momentum both, $A$ and $\hat{q}$ , approach unity. If the partons have less kinetic energy or a smaller pole mass as compared to a system in equilibrium charm quarks lose more energy. In the former case $\hat{q}$ decreases whereas in the latter case it increases for charm quark with a low or intermediate transverse momentum. Thus each non-equilibrium scenario affects $A$ and $\hat{q}$ of charm quarks in a different way. The modifications in our scenarios are of the order 20-50\% at temperatures relevant for heavy ion reactions. These modifications have to be considered if one wants to determine these coefficients by comparing heavy ion data with theoretical predictions from viscous hydrodynamics or Langevin equations.

preprint2020arXiv

Open charm and dileptons from relativistic heavy-ion collisions

We study the dynamics of open charm production and the dilepton radiation of the semi-leptonic decays of correlated $D\bar D$ pairs versus the quark-gluon plasma (QGP) radiation and hadronic sources in relativistic heavy-ion collisions. Our study is based on the Parton-Hadron-String Dynamics (PHSD) transport approach employing a non-perturbative QCD description of the strongly interacting quark-gluon plasma (sQGP) in terms of dynamical quasiparticles and the EoS based on lattice QCD. We compare the PHSD results for charm observables with the calculations from BAMPS (Boltzmann Approach to Multi-Parton Scatterings) which is based on perturbative QCD with massless partons and interaction cross sections calculated in leading order of the QCD coupling. We compare the $p_T$ dependence of the ratio $R_{AA}$ of $D$-mesons in $A+A$ over $p+p$ collisions scaled by the number of binary collisions $N_{bin}$ as well as the elliptic flow $v_2$ of $D$-mesons calculated within the PHSD and BAMPS at LHC energies. In other study, based on the PHSD calculations we find that the dileptons from correlated $D-$meson semi-leptonic decays dominate the 'thermal' radiation from the QGP in central Pb+Pb collisions at the intermediate masses ($1.2 < M < 3$ GeV) for higher invariant energies However, for invariant energies $\sqrt{s_{NN}} < 40$ GeV the QGP radiation overshines the contribution from $D,{\bar D}$ decays such that one should observe a rather clear signal from the partonic dilepton radiation. This finding provides promising perspectives to measure the QGP radiation in the dilepton experiments at RHIC BES and the future FAIR/NICA facilities.

preprint2020arXiv

Traces of non-equilibrium effects, initial condition, bulk dynamics and elementary collisions in the charm observables

We study how heavy quark dynamics is affected by the nature of the bulk evolution of the QCD matter, the initial condition of the system, and the treatment of elementary interactions between heavy quarks and the surrounding medium. For the same initial condition and the same QGP expansion scenario we discuss the consequences of the assumption of a local equilibrium by comparing the consequences for the nuclear modification factor and the elliptic flows of charm quarks. For this purpose we employ the parton-hadron-string dynamics (PHSD), which is an off-shell microscopic transport approach, as well as the linearized-Boltzmann (LB) scheme obtained by coarse graining the PHSD bulk and assuming local equilibrium for the interactions of the charm quarks with the bulk. The $R_{\rm AA}$ of charm quarks stemming from the later LB approach is also compared to a genuine fluid dynamics evolution initiated by the coarse grained PHSD, which allows to further assess the consequences of reducing the full n-body dynamics. We then proceed to a systematic comparison of PHSD with MC@HQ, another transport model for heavy flavors which also relies on LB approach. In particular, we investigate the consequences for the nuclear modification factor of charm quarks if we vary separately the initial heavy quark distribution function in matter, the expansion dynamics of the QGP and the elementary interactions of heavy quarks of these models. We find that the results for both models vary significantly depending on the details of the calculation. However, both models achieve very similar predictions for key heavy quark observables for certain combinations of initial condition, bulk evolution and interactions. We conclude that this ambiguity limits our ability to determine the different properties of the system based on the current set of observables.

preprint2016arXiv

Charm production in Pb+Pb collisions at the Large Hadron Collider energy

We study charm production in Pb+Pb collisions at $\sqrt{s_{\rm NN}}=$2.76 TeV in the Parton-Hadron-String-Dynamics transport approach and the charm dynamics in the partonic and hadronic medium. The charm quarks are produced through initial binary nucleon-nucleon collisions by using the PYTHIA event generator taking into account the (anti-)shadowing incorporated in the EPS09 package. The produced charm quarks interact with off-shell massive partons in the quark-gluon plasma and are hadronized into $D$ mesons through coalescence or fragmentation close to the critical energy density, and then interact with hadrons in the final hadronic stage with scattering cross sections calculated in an effective Lagrangian approach with heavy-quark spin symmetry. The PHSD results show a reasonable $R_{\rm AA}$ and elliptic flow of $D$ mesons in comparison to the experimental data for Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV from the ALICE Collaboration. We also study the effect of temperature-dependent off-shell charm quarks in relativistic heavy-ion collisions. We find that the scattering cross sections are only moderately affected by off-shell charm degrees of freedom. However, the position of the peak of $R_{\rm AA}$ for $D$ mesons depends on the strength of the scalar partonic forces which also have an impact on the $D$ meson elliptic flow. The comparison with experimental data on the $R_{\rm AA}$ suggests that the repulsive force is weaker for off-shell charm quarks as compared to that for light quarks. Furthermore, the effects from radiative charm energy loss appear to be low compared to the collisional energy loss up to transverse momenta of $\sim$ 15 GeV/c.

preprint2016arXiv

The hadronization time of heavy quark in nuclear matter

We study the hadronization time of heavy quark in nuclear matter by using the coalescence model and the spatial diffusion constant of heavy quark from lattice Quantum Chromodynamic calculations, assuming that the main interaction of heavy quark at the critical temperature is hadronization. It is found that the hadronization time of heavy quark is about 3 fm/c for $2πT_c D_s=6$, if a heavy quark is combined with the nearest light antiquark in coordinate space without any correlation between momentum of heavy quark and that of light antiquark which form a heavy meson. However, the hadronization time reduces to 0.6-1.2 fm/c for charm and 0.4-0.9 fm/c for bottom, depending on heavy meson radius, in the presence of momentum correlation. Considering the interspace between quarks and antiquarks at the critical temperature, it seems that the hadronization of heavy quark does not happen instantaneously but gradually for a considerable time, if started from the thermal distribution of quarks and antiquarks.

preprint2016arXiv

Understanding transport simulations of heavy-ion collisions at 100 and 400 AMeV: Comparison of heavy ion transport codes under controlled conditions

Transport simulations are very valuable for extracting physics information from heavy-ion collision experiments. With the emergence of many different transport codes in recent years, it becomes important to estimate their robustness in extracting physics information from experiments. We report on the results of a transport code comparison project. 18 commonly used transport codes were included in this comparison: 9 Boltzmann-Uehling-Uhlenbeck-type codes and 9 Quantum-Molecular-Dynamics-type codes. These codes have been required to simulate Au+Au collisions using the same physics input for mean fields and for in-medium nucleon-nucleon cross sections, as well as the same initialization set-up, the impact parameter, and other calculational parameters at 100 and 400 AMeV incident energy. Among the codes we compare one-body observables such as rapidity and transverse flow distributions. We also monitor non-observables such as the initialization of the internal states of colliding nuclei and their stability, the collision rates and the Pauli blocking. We find that not completely identical initializations constitute partly for different evolutions. Different strategies to determine the collision probabilities, and to enforce the Pauli blocking, also produce considerably different results. There is a substantial spread in the predictions for the observables, which is much smaller at the higher incident energy. We quantify the uncertainties in the collective flow resulting from the simulation alone as about $30\%$ at 100 AMeV and $13\%$ at 400 AMeV, respectively. We propose further steps within the code comparison project to test the different aspects of transport simulations in a box calculation of infinite nuclear matter. This should, in particular, improve the robustness of transport model predictions at lower incident energies where abundant amounts of data are available.

preprint2015arXiv

Heavy quark diffusion from coherent color fields in relativistic heavy-ion collisions

The diffusion coefficients of heavy quarks from the coherent color electromagnetic fields which are generated in the early stage of relativistic heavy-ion collisions are calculated at midrapidity, and compared with those obtained from collisions within a thermalized quark-gluon plasma. The coherent color fields are modeled such that they are initially longitudinal and then become isotropic. We found that the diffusion coefficients from the coherent color fields are larger than those from collisions except for very fast heavy quarks, and the color fields are less effective for heavy-quark energy loss. The importance of coherent color fields for heavy-quark diffusion decreases as energy density decreases.

preprint2015arXiv

Mean-field potential effects on particle and antiparticle elliptic flows in the beam-energy scan program at RHIC

The elliptic flow splitting between particles and their antiparticles has recently been observed by the STAR Collaboration in the beam-energy scan program at the Relativistic Heavy Ion Collider. In studies based on transport models, we have found that this splitting can be explained by the different mean-field potentials acting on particles and their antiparticles in the produced baryon-rich matter. In particular, we have shown that the experimentally measured relative elliptic flow difference can help constrain the vector coupling constant in the Nambu-Jona-Lasinio model used in describing the partonic stage of heavy-ion collisions. This information is useful for locating the critical point in the QCD phase diagram and thus understanding the phase structure of QCD.

preprint2015arXiv

Quarkonium formation time in relativistic heavy-ion collisions

We calculate the quarkonium formation time in relativistic heavy-ion collisions from the space-time correlator of heavy quark vector currents in a hydrodynamics background with the initial nonequilibrium stage expanding only in the longitudinal direction. Using in-medium quarkonia properties determined with the heavy quark potential taken to be the free energy from lattice calculations and the fact that quarkonia can only be formed below their dissociation temperatures due to color screening, we find that $Υ$(1S), $Υ$(2S), $Υ$(3S), $J/ψ$ and $ψ^\prime$ are formed, respectively, at 1.2, 6.6, 8.8, 5.8, and 11.0 fm/c after the quark pair are produced in central Au+Au collisions at the top energy of Relativistic Heavy Ion Collider (RHIC), and these times become shorter in semi-central collisions. We further show, as an example, that including the effect of formation time enhances appreciably the survivability of $Υ$(1S) in the produced hot dense matter.

preprint2015arXiv

Tomography of the Quark-Gluon-Plasma by Charm Quarks

We study charm production in ultra-relativistic heavy-ion collisions by using the Parton-Hadron-String Dynamics (PHSD) transport approach. The initial charm quarks are produced by the Pythia event generator tuned to fit the transverse momentum spectrum and rapidity distribution of charm quarks from Fixed-Order Next-to-Leading Logarithm (FONLL) calculations. The produced charm quarks scatter in the quark-gluon plasma (QGP) with the off-shell partons whose masses and widths are given by the Dynamical Quasi-Particle Model (DQPM) which reproduces the lattice QCD equation-of-state in thermal equilibrium. The relevant cross section are calculated in a consistent way by employing the effective propagators and couplings from the DQPM. Close to the critical energy density of the phase transition, the charm quarks are hadronized into $D$ mesons through coalescence and/or fragmentation depending on transverse momentum. The hadronized $D$ mesons then interact with the various hadrons in the hadronic phase with cross sections calculated in an effective lagrangian approach with heavy-quark spin symmetry. Finally, the nuclear modification factor $\rm R_{AA}$ and the elliptic flow $v_2$ of $D^0$ mesons from PHSD are compared with the experimental data from the STAR Collaboration for Au+Au collisions at $\sqrt{s_{\rm NN}}$ =200 GeV. We find that in the PHSD the energy loss of $D$ mesons at high $p_T$ can be dominantly attributed to partonic scattering while the actual shape of $\rm R_{AA}$ versus $p_T$ reflects the heavy quark hadronization scenario, i.e. coalescence versus fragmentation. Also the hadronic rescattering is important for the $\rm R_{AA}$ at low $p_T$ and enhances the $D$-meson elliptic flow $v_2$.

preprint2014arXiv

Charmonia formation in quark-gluon plasma

Using the color evaporation model, the cross section for charmonium production in p+p collision is calculated in quark-gluon plasma. The threshold energy for open charms is given by the free energy potential from lattice calculations, the initial charm quark pairs by the Pythia simulations, and their time evolution by solving the Langevin equation. It is found that the threshold energy which decreases with temperature reduces the cross section while the invariant mass of charm pair which decreases by collisions enhances it. As a result, charmonia production is suppressed by $30$\sim$50\%$ while $J/ψ$ production is similar or enhanced compared to in vacuum.

preprint2014arXiv

Double ratio of charmonia in p+Pb collisions at sqrt(s_NN)=5.02 TeV

Based on a kinetic description of J/psi dissociation and production in an expanding quark-gluon plasma that is described by a 2+1 dimensional ideal hydrodynamics, we have studied the double ratio R_{p Pb}(psi')/R_{p Pb}(J/psi) of charmonia in p+Pb collisions at sqrt(s_NN)=5.02 TeV by including not only the cold nuclear matter effects but also the hot nuclear matter effects. We find that the double ratio of prompt charmonia is significantly suppressed in the most central collisions as a result of the hot nuclear matter effects.

preprint2014arXiv

Effects of medium modification of pion production threshold in heavy ion collisions and the nuclear symmetry energy

Using the relativistic Vlasov-Uehling-Uhlenbeck (RVUU) equation based on mean fields from the nonlinear relativistic models, we study the effect of medium modification of pion production threshold on the total pion yield and the $π^-/π^+$ ratio in Au+Au collisions. We find that the in-medium threshold effect enhances both the total pion yield and the $π^-/π^+$ ratio, compared to those without this effect. Furthermore, including the medium modification of the pion production threshold in asymmetric nuclear matter leads to a larger $π^-/π^+$ ratio for the $NLρδ$ model with a stiffer symmetry energy than the $NLρ$ model with a softer symmetry energy, opposite to that found without the in-medium threshold effect. Experimental data from the FOPI Collaboration are reproduced after including a density-dependent cross section for $Δ$ baryon production from nucleon-nucleon collisions, which suppresses the total pion yield but hardly changes the $π^-/π^+$ ratio. The large errors in the experimentally measured $π^-/π^+$ ratio prevent, however, the distinction between the predictions from the $NLρ$ and $NLρδ$ models.

preprint2014arXiv

Partonic mean-field effects on matter and antimatter elliptic flows

Using a partonic transport model based on the Nambu-Jona-Lasinio model, we study the effect of scalar and vector mean fields on the elliptic flows of quarks and antiquarks in relativistic heavy ion collisions in Au+Au collisions at $\sqrt{s_{\rm NN}}=7.7~$GeV and impact parameter $b=8~{\rm fm}$ that leads to the production of a baryon-rich matter. Although the scalar mean field, which is attractive for both quarks and antiquarks, reduces both their elliptic flows, the vector mean field, which is repulsive for quarks and attractive for antiquarks, leads to a splitting of their elliptic flows, and this effect increases with the strength of the vector coupling in the baryon-rich quark matter. Converting quarks and antiquarks at hadronization to hadrons via the quark coalescence model, we further study the dependence of the transverse momentum integrated relative elliptic flow differences between protons and antiprotons, lambda and anti-lambdas, and positively and negatively charged kaons on the strength of the quark vector coupling. These results are then compared with the experimental data measured by the STAR Collaboration in the Beam Energy Scan program at the Relativistic Heavy Ion Collider.

preprint2013arXiv

Elliptic flow splitting as a probe of the QCD phase structure at finite baryon chemical potential

Using a partonic transport model based on the 3-flavor Nambu-Jona-Lasinio model and a relativistic hadronic transport model to describe, respectively, the evolution of the initial partonic and the final hadronic phase of heavy-ion collisions at energies carried out in the Beam-Energy Scan program of the Relativistic Heavy Ion Collider, we have studied the effects of both the partonic and hadronic mean-field potentials on the elliptic flow of particles relative to that of their antiparticles. We find that to reproduce the measured relative elliptic flow differences between nucleons and antinucleons as well as between kaons and antikaons requires a vector coupling constant as large as 0.5 to 1.1 times the scalar coupling constant in the Nambu-Jona-Lasinio model. Implications of our results in understanding the QCD phase structure at finite baryon chemical potential are discussed.

preprint2013arXiv

Free energy versus internal energy potential for heavy quark systems at finite temperature

Using the QCD sum rule with its operator product expansion reliably estimated from lattice calculations for the pressure and energy density of the quark-gluon plasma, we calculate the strength of the wave function of J/psi at the origin and find that it decreases with temperature above the transition temperature. This result is shown to follow exactly that obtained from the solution of the Schrödingier equation for a charm and anticharm quark pair using the free energy from lattice calculations as the potential and is in sharp contrast to that using the deeper potential associated with the internal energy, which shows an enhanced strength of the J/psi wave function at the origin. Our result thus has resolved the long-standing question of whether the free energy potential or the internal energy potential should be used in analyzing the spectrum of heavy quark systems at finite temperature.

preprint2013arXiv

Gluon dissociation of J/psi beyond the dipole approximation

Using a nonrelativistic potential model, we calculate the cross section for the leading-order gluon dissociation of J/psi by including the full gluon wave function. We find that the resulting cross section as a function of gluon energy is reduced by about a factor of three at its maximum value compared to that calculated in the dipole approximation that is usually adopted in theoretical studies. The effect of the reduced cross section on the J/psi dissociation width at finite temperature is also discussed.

preprint2013arXiv

Hot medium effects on J/psi production in p+Pb collisions at sqrt{s_{NN}}=5.02 TeV

Based on a kinetic description of J/psi dissociation and production in an expanding quark-gluon plasma that is described by a 2+1 dimensional ideal hydrodynamics, we have studied the hot medium effects on J/psi production in p+Pb collisions at sqrt{s_{NN}}=5.02 TeV. Including also the cold nuclear matter effects, we are able to reproduce recent experimental results on the nuclear modification factor R_{pPb}(J/psi) measured by the ALICE Collaboration. We have also made predictions for the R_{pPb} of J/psi and the double ratio R_{pPb}^{pro}(psi')/R_{pPb}^{pro}(J/psi) of prompt quarkonia produced in the most central 10% p+Pb collisions. We find that different from the cold nuclear matter effects, the R_{pPb}(J/psi) is slightly smaller than that in the minimum bias collisions, and the double ratio is significantly less than one at backward rapidity.

preprint2012arXiv

Charmonium production from nonequilibrium charm and anticharm quarks in quark-gluon plasma

Parameterizing the charm and anticharm quark momentum distributions by the Tsallis distribution, we study the nonequilibrium effect on the charmonium production rate in a quark-gluon plasma up to the next-to-leading order in perturbative QCD. We find that nonequilibrium charm and anticharm quarks suppress the charmonium production rate compared to that from equilibrated ones. We further show that the suppression factor calculated with the charm quark relaxation time, which has been frequently used in the literature, is close to our results.

preprint2011arXiv

Bottomonia suppression in heavy-ion collisions

Using the two-component model that includes both initial production from nucleon-nucleon hard scattering and regeneration from produced quark-gluon plasma (QGP), we study the effect of medium modifications of the binding energies and radii of bottomonia on their production in heavy-ion collisions. We find that the contribution to bottomonia production from regeneration is small and the inclusion of medium effects is generally helpful for understanding the observed suppression of bottomonia production in experiments carried out at both the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC).

preprint2011arXiv

Charmonium production in heavy-ion collisions from SPS to LHC

Using the two-component model that includes charmonium production from initial nucleon-nucleon hard scattering and regeneration in the produced quark-gluon plasma, we study $J/ψ$ production in heavy ion collisions at SPS, RHIC and LHC. For the expansion dynamics of produced hot dense matter, we use a schematic viscous hydrodynamic model with the specific shear viscosity taken, respectively, to be twice and ten times the lower bound of $1/4π$ suggested by the Ads/CFT correspondence in the quark-gluon plasma and in the hadron gas. For the initial dissociation and the subsequent thermal decay of charmonia in the hot dense matter, we use the screened Cornell potential to describe the properties of charmonia and the perturbative QCD to calculate their dissociation cross sections. Including regeneration of charmonia in the quark-gluon plasma via a kinetic equation with in-medium chamonium decay widths, we obtain a good description of measured nuclear modification factors at SPS and RHIC. A reasonable description of measured nuclear modification factor of $J/ψ$s of high transverse momenta in heavy ion collisions at LHC is also obtained if we include the contribution to the $J/ψ$ yield from the decay of bottomed hadrons.

preprint2011arXiv

The effect of initial fluctuations on bottomonia suppression in relativistic heavy-ion collisions

Using the screened Cornell potential and the next-to-leading order perturbative QCD to determine, respectively, the properties of bottomonia and their dissociation cross sections in a quark-gluon plasma, we study in a 2+1 ideal hydrodynamics the effect of initial fluctuations on bottomonia production in relativistic heavy-ion collisions. We find that while initial fluctuations hardly affect the yield of the 1S ground state bottomonium, their effect on that of excited bottomonium states is not small. Compared to the case with smooth initial conditions, the survival probability of excited bottomonia is reduced at low transverse momentum and increased at high transverse momentum. The observed suppression of the excited bottomonia relative to the ground state bottomonium by the Compact Muon Solenoid (CMS) collaborations at an average transverse momentum can, however, be described at present with both smooth and fluctuating initial conditions.

preprint2011arXiv

The flow of heavy flavor in hydrodynamics

The flow of charm is calculated in 2+1 ideal hydrodynamics by introducing the charge of $c\bar{c}$ pair assuming that the number of $c\bar{c}$ pairs is conserved in relativistic heavy-ion collisions. It is found that the mean radial flow velocity of charm quarks is smaller than that of bulk matter by 10$\sim$15 \% and the measured $v_2$ of heavy-flavor electrons is reproduced up to $p_T^e=$ 1.5 GeV/c in Au+Au collision at RHIC. The same flow is applied to regenerated $J/ψ$ and its $v_2$ is discussed.

preprint2010arXiv

$J/ψ$ production and elliptic flow in relativistic heavy-ion collisions

Using a two-component model for charmonium production, which includes contributions from both the initial hard nucleon-nucleon scattering and from the regeneration in the quark-gluon plasma, we study the nuclear modification factor $R_{AA}$ and elliptic flow $v_2$ of $J/ψ$ in relativistic heavy ion collisions. For the expansion dynamics of produced hot dense matter, we introduce a schematic fireball model with its transverse acceleration determined from the pressure gradient inside the fireball and azimuthally anisotropic expansion parameterized to reproduce measured $v_2$ of light hadrons. We assume that light hadrons freeze out at the temperature of 120 MeV while charmonia at 160 MeV, similar to the kinetic and chemical freeze-out temperatures in the statistical model, respectively. For the properties of charmonia in the quark-gluon plasma, we use the screening mass between their charm and anticharm quarks and their dissociation cross sections given by the perturbative QCD (pQCD) in the leading order and up to the next-to-leading order, respectively. For the relaxation time of charm and anticharm quarks in the quark-gluon plasma, we also use the one calculated in the leading order of pQCD. Modeling the effect of higher-order corrections in pQCD by introducing multiplicative factors to the dissociation cross section of charmonia and the elastic scattering cross sections of charm and anticharm quarks, we find that this effect is small for the $R_{AA}$ of $J/ψ$ as they suppress the number of initially produced $J/ψ$ but enhance the number of regenerated ones. The higher-order corrections increase, however, the $v_2$ of $J/ψ$. Our results suggest that the $v_2$ of $J/ψ$ can play an important role in discriminating between $J/ψ$ production from the initial hard collisions and from the regeneration in the quark-gluon plasma.

preprint2010arXiv

$R_{AA}$ of $J/ψ$ near mid-rapidity in heavy ion collisions at $\sqrt{s_{NN}}=200$ GeV

We build up a model to reproduce the experimentally measured $R_{AA}$ of $J/ψ$ near midrapidty in Au+Au collision at $\sqrt{s_{NN}}=200$ GeV. The model takes into account the $J/ψ$ suppression from the quark-gluon plasma and hadron gas as well as the nuclear absorption of primordial charmonia and the regeneration effects at the hadronization stage, and hence is a generalization of the two component model introduced by Grandchamp and Rapp. The improvements in this work are twofold; the addition of the initial local temperature profile and a consistent use of QCD NLO formula for both the dissociation cross section in the hadron gas and the thermal decay widths in the quark-gluon plasma for the charmonium states. The initial local temperature profile is determined from the assumption that the local entropy density is proportional to a formula involving the number densities of the number of participants and of the binary collisions that reproduces the multiplicities of charged particles at chemical freeze-out. The initial local temperature profile brings about a kink in the $R_{AA}$ curve due to the initial melting of $J/ψ$. The initially formed fireball, composed of weakly interacting quarks and gluons with thermal masses that are extracted from lattice QCD, follows an isentropic expansion with cylindrical symmetry. The fit reproduces well the Au+Au as well as the Cu+Cu data. The same method is applied to predict the $R_{AA}$ expected from the Pb+Pb collision at LHC energy.

preprint2010arXiv

Dilepton production in a schematic causal viscous hydrodynamics

Assuming that in the hot dense matter produced in relativistic heavy-ion collisions, the energy density, entropy density, and pressure as well as the azimuthal and space-time rapidity components of the shear tensor are uniform in the direction transversal to the reaction plane, we derive a set of schematic equations from the Isreal-Stewart causal viscous hydrodynamics. These equations are then used to describe the evolution dynamics of relativistic heavy-ion collisions by taking the shear viscosity to entropy density ratio of $1/4π$ for the initial quark-gluon plasma (QGP) phase and of ten times this value for the later hadron-gas (HG) phase. Using the production rate evaluated with particle distributions that take into account the viscous effect, we study dilepton production in central heavy-ion collisions. Compared with results from the ideal hydrodynamics, we find that although the dilepton invariant mass spectra from the two approaches are similar, the transverse momentum spectra are significantly enhanced at high transverse momenta by the viscous effect. We also study the transverse momentum dependence of dileptons produced from QGP for a fixed transverse mass which is essentially absent in the ideal hydrodynamics, and find that this so-called transverse mass scaling is violated in the viscous hydrodynamics, particularly at high transverse momenta.