Source author record

Che Ming Ko

Che Ming Ko 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

68works
11topics
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

68 published item(s)

preprint2026arXiv

Role of $Σ^*(1385)$ on $Λ$ hyperon polarization in relativistic heavy ion collisions

The effect of $Σ^*(1385)$ baryon resonance on the time evolution of the $Λ$ hyperon polarization in hadronic matter is studied using a kinetic approach. This approach explicitly includes the production of the $Σ^*$ resonance from the $Λ-π$ and $Σ(1192)-π$ scatterings as well as its decay into $Λ+π$ or $Σ+π$. The resulting coupled kinetic equations governing the time evolution of $Λ$, $Σ$ and $Σ^*$ numbers and polarizations are solved for Au-Au collisions at $\sqrt{s_{NN}}=7.7$ GeV and 20-50\% centrality, using initial values determined by thermal yields and the thermal vorticity at chemical freeze-out temperature. As the hadronic matter expands and cools, the $Λ$ polarization is found to increase slightly during early times and then decreases very slowly afterwards, while the $Σ$ polarization remains nearly constant and the $Σ^*$ polarization continuously decreases. Including feed-down contributions to the $Λ$ polarization from the decays of partially polarized $Σ^0$, $Σ^*$, and $Ξ(1322)$ hyperons, where the $Ξ$ polarization is obtained by solving coupled kinetic equations for the $Ξ$ and $Ξ^*(1532)$ system, the resulting $Λ$ polarization becomes smaller and decreases over time. In both cases, however, the time variation of the $Λ$ polarization is sufficiently small to support the assumption of an early freeze-out of $Λ$ spin degree of freedom in relativistic heavy ion collisions.

preprint2022arXiv

Charged pion production from Au + Au collisions at $\sqrt{s_{NN}}=2.4$ GeV in the Relativistic Vlasov-Uehling-Uhlenbeck model

Using the isospin-dependent relativistic Vlasov-Uehling-Uhlenbeck (RVUU) model, we study charged pion ($π^\pm$) production in Au+Au collisions at $\sqrt{s_{NN}}=$ 2.4 GeV. By fitting the density dependence of the $Δ$ resonance production cross section in nuclear medium to reproduce the experimental $π^\pm$ multiplicities measured by the HADES Collaboration, we obtain a good description of the rapidity distributions and transverse momentum spectra of $π^\pm$ in collisions at various centralities. Some shortcomings in the description of $π^{\pm}$ production may indicate the need for including the strong potential on $π^\pm$ in RVUU, which is at present absent. We also calculate the proton rapidity distribution in the most central collisions and compare with the coalescence invariant proton rapidity distribution extracted from preliminary HADES data.

preprint2022arXiv

Spinodal Enhancement of Light Nuclei Yield Ratio in Relativistic Heavy Ion Collisions

Using a relativistic transport model to describe the evolution of the quantum chromodynamic matter produced in Au+Au collisions at $\sqrt{s_{NN}}=3-200$ GeV, we study the effect of a first-order phase transition in the equation of state of this matter on the yield ratio $N_tN_p/ N_d^2$ ($tp/d^2$) of produced proton ($p$), deuteron ($d$), and triton ($t$). We find that the large density inhomogeneities generated by the spinodal instability during the first-order phase transition can survive the fast expansion of the subsequent hadronic matter and lead to an enhanced $tp/d^2$ in central collisions at $\sqrt{s_{NN}}=3-5$ GeV as seen in the experiments by the STAR Collaboration and the E864 Collaboration. However, this enhancement subsides with increasing collision centrality, and the resulting almost flat centrality dependence of $tp/d^2$ at $\sqrt{s_{NN}}=3$ GeV can also be used as a signal for the first-order phase transition.

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.

preprint2022arXiv

Unveiling the dynamics of nucleosynthesis in relativistic heavy-ion collisions

Like nucleosynthesis during the early universe, light nuclei are also produced in relativistic heavy-ion collisions. Although the deuteron ($d$) yields in these collisions can be well described by the statistical hadronization model (SHM), which assumes that particle yields are fixed at a common chemical freezeout near the phase boundary between the quark-gluon plasma and the hadron gas, the recently measured triton ($^3\text{H}$) yields in Au+Au collisions at $\sqrt{s_{NN}}=7.7-200$ GeV are overestimated systematically by this model. Here, we develop a comprehensive kinetic approach to study the effects of hadronic re-scatterings, such as $πNN\leftrightarrowπd$ and $πNNN\leftrightarrowπ^3\text{H}~(^3\text{He})$, on $d$, $^3\text{H}$, and $^3\text{He}$ production in these collisions. We find that these reactions have little effects on the deuteron yield but reduce the $^3\text{H}$ and $^3\text{He}$ yields by about a factor of 1.8 from their initial values given by the SHM. This finding helps resolve the overestimation of triton production in the SHM and provides the evidence for hadronic re-scattering effects on nucleosynthesis in relativistic heavy-ion collisions.

preprint2021arXiv

Angular Momentum Eigenstates of the Isotropic 3-D Harmonic Oscillator: Phase-Space Distributions and Coalescence Probabilities

The isotropic 3-dimensional harmonic oscillator potential can serve as an approximate description of many systems in atomic, solid state, nuclear, and particle physics. In particular, the question of 2 particles binding (or coalescing) into angular momentum eigenstates in such a potential has interesting applications. We compute the probabilities for coalescence of two distinguishable, non-relativistic particles into such a bound state, where the initial particles are represented by generic wave packets of given average positions and momenta. We use a phase-space formulation and hence need the Wigner distribution functions of angular momentum eigenstates in isotropic 3-dimensional harmonic oscillators. These distribution functions have been discussed in the literature before but we utilize an alternative approach to obtain these functions. Along the way, we derive a general formula that expands angular momentum eigenstates in terms of products of 1-dimensional harmonic oscillator eigenstates.

preprint2021arXiv

Evolution of $Λ$ polarization in the hadronic phase of heavy-ion collisions

Using the AMPT + MUSIC+UrQMD hybrid model, we study the global and local spin polarizations of $Λ$ hyperons as functions of the freeze-out temperature of the spin degree of freedom in the hadronic phase of Au+Au collisions at $\sqrt{s_{NN}}=19.6$ GeV. Including contributions from both the thermal vorticity and thermal shear of the hadronic matter, we find that, with the spin freeze-out temperature dropping from the hadronization temperature of 160 MeV to 110 MeV at the kinetic freeze-out, both the global and local spin polarizations of $Λ$ hyperons due to the thermal vorticity decrease by a factor of two, while those due to the thermal shear decrease quickly and become negligibly small at 140 MeV. Our results suggest the importance of understanding the dynamical evolution of the spin degree of freedom in the hadronic stage in relativistic heavy-ion collisions.

preprint2020arXiv

Charmed hadron production in an improved quark coalescence model

We study the production of charmed hadrons $D^{0}$ and $Λ_c^+$ in relativistic heavy-ion collisions using an improved quark coalescence model. In particular, we extend the usual coalescence model by letting a produced hadron to have the same velocity as the center-of-mass velocity of coalesced constituent quarks during hadronization to take into account the effect of collective flow in produced quark-gluon plasma. This results in a shift of charmed resonances of higher masses to larger transverse momenta ($p_T^{}$). Requiring all charm quarks of very low $p_T^{}$ to be converted to hadrons via coalescence and letting charm quarks not undergoing coalescence to hadronize by independent fragmentation, we obtain a good description of the measured yield ratio $Λ_c^+/D^0$ as a function of $p_T^{}$ in $\text{Au} + \text{Au}$ collisions at $\sqrt{s_{NN}}^{}=200$~GeV by the STAR Collaboration at the Relativistic Heavy Ion Collider.

preprint2020arXiv

Constraining the in-medium nucleon-nucleon cross section from the width of nuclear giant dipole resonance

We develop a new lattice Hamiltonian method for solving the Boltzmann-Uehling-Uhlenbeck (BUU) equation. Adopting the stochastic approach to treat the collision term and using the GPU parallel computing to carry out the calculations allows for a rather high accuracy in evaluating the collision term, especially its Pauli blocking, leading thus to a new level of precision in solving the BUU equation. Applying this lattice BUU method to study the width of giant dipole resonance (GDR) in nuclei, where the accurate treatment of the collision term is crucial, we find that the obtained GDR width of $^{208}{\rm Pb}$ shows a strong dependence on the in-medium nucleon-nucleon cross section $σ_{\rm NN}^*$. A very large medium reduction of $σ_{\rm NN}^*$ is needed to reproduce the measured value of the GDR width of $^{208}{\rm Pb}$ at the Research Center for Nuclear Physics in Osaka, Japan.

preprint2020arXiv

Local spin polarizations in relativistic heavy ion collisions

Based on a generalized side-jump formalism for massless chiral fermions, which naturally takes into account the spin-orbit coupling in the scattering of two chiral fermions and the chiral vortical effect in a rotating chiral fermion matter, we have developed a covariant and total angular momentum conserved chiral transport model to study both the global and local polarizations of this matter. For a system of massless quarks of random spin orientations and finite vorticity in a box, we have demonstrated that the model can exactly conserve the total angular momentum of the system and dynamically generate the quark spin polarization expected from a thermally equilibrated quark matter. Using this model to study the spin polarization in relativistic heavy-ion collision, we have found that the local quark spin polarizations depend strongly on the reference frame where they are evaluated as a result of the nontrivial axial charge distribution caused by the chiral vortical effect. We have further shown that because of the anomalous orbital or side-jump contribution to the quark spin polarization, the local quark polarizations calculated in the medium rest frame are qualitatively consistent with the local polarizations of Lambda hyperons measured in experiments.

preprint2020arXiv

Number of constituent quark scaling of elliptic flows in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}=$ 5.02 TeV

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.

preprint2020arXiv

Probing the partonic degrees of freedom in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}=$ 5.02 TeV

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.

preprint2020arXiv

Properties of strange quark stars with isovector interactions

We study the properties of strange quark stars by employing a 3-flavor Nambu-Jona-Lasinio model with both scalar-isovector and vector-isovector interactions. Using the constraint on the vector-isoscalar interaction strength obtained from the elliptic flow splitting between particles and their antiparticles in relativistic heavy-ion collisions, we investigate the dependence of the properties of strange quark stars on the vector-isovector and the scalar-isovector interactions, and compare the results with the state-of-art astrophysical constraints on the compact star radius and mass as well as its tidal deformability from the GW170817 event. Results from our study reinforce the prospect of using both heavy-ion collisions and astrophysical observations to provide constraints on the isovector coupling strength in quark matter and thus the quark matter equation of state as well as the QCD phase structure at finite isospin chemical potentials.

preprint2020arXiv

Spin Polarizations in a Covariant Angular-Momentum-Conserved Chiral Transport Model

Using a covariant and angular-momentum-conserved chiral transport model, which takes into account the spin-orbit interactions of chiral fermions in their scatterings via the side jumps, we study the quark spin polarization in quark matter. For a system of rotating and unpolarized massless quarks in an expanding box, we find that side jumps can dynamically polarize the quark spin and result in a final quark spin polarization consistent with that of thermally equilibrated massless quarks in a self-consistent vorticity field. For the quark matter produced in noncentral relativistic heavy ion collisions, we find that in the medium rest frame both the quark local spin polarizations in the direction perpendicular to the reaction plane and along the longitudinal beam direction show an azimuthal angle dependence in the transverse plane similar to those observed in experiments for the Lambda hyperon.

preprint2020arXiv

Suppression of light nuclei production in collisions of small systems at the Large Hadron Collider

We show that the recently observed suppression of the yield ratio of deuteron to proton and of helium-3 to proton in p+p collisions compared to those in p+Pb or Pb+Pb collisions by the ALICE Collaboration at the Large Hadron Collider (LHC) can be explained if light nuclei are produced from the coalescence of nucleons at the kinetic freeze-out of these collisions. This suppression is attributed to the non-negligible sizes of deuteron and helium-3 compared to the size of the nucleon emission source in collisions of small systems, which reduces the overlap of their internal wave functions with those of nucleons. The same model is also used to study the production of triton and hypertriton in heavy-ion collisions at the LHC. Compared to helium-3 in events of low charged particle multiplicity, the triton is less suppressed due to its smaller size and the hypertriton is even more suppressed as a result of its much larger size.

preprint2019arXiv

Probing QCD critical fluctuations from the yield ratio of strange hadrons in relativistic heavy-ion collisions

By analyzing the available data on strange hadrons in central Pb+Pb collisions from the NA49 Collaboration at the Super Proton Synchrotron (SPS) and in central Au+Au collisions from the STAR Collaboration at the Relativistic Heavy-Ion Collider (RHIC) in a wide collision energy range from $\sqrt{s_{\rm NN}}$ = 6.3 GeV to 200 GeV, we find a possible non-monotonic behavior in the ratio $\mathcal{O}_\text{K-$Ξ$-$ϕ$-$Λ$}$= $\frac{N(K^+)N(Ξ^-)}{N(ϕ)N(Λ)}$ of $K^+$, $Ξ^-$, $ϕ$, and $Λ$ yields as a function of $\sqrt{s_{\rm NN}}$. Based on the quark coalescence model, which can take into account the effect of quark density fluctuations on hadron production, a possible non-monotonic behavior in the dependence of the strange quark density fluctuation on $\sqrt{s_{NN}}$ is obtained. This is in contrast to the coalescence model that does not include quark density fluctuations and also to the statistical hadronization model as both fail to describe even qualitatively the collision energy dependence of the ratio $\mathcal{O}_\text{K-$Ξ$-$ϕ$-$Λ$}$. Our findings thus suggest that the signal and location of a possible critical endpoint in the QCD phase diagram, which is expected to result in large quark density fluctuations, can be found in the on-going Bean Energy Scan program at RHIC.

preprint2016arXiv

Jet Fragmentation via Recombination of Parton Showers

We propose to model hadronization of parton showers in QCD jets through a hybrid approach involving quark recombination and string fragmentation. This is achieved by allowing gluons at the end of the perturbative shower evolution to undergo a non-perturbative splitting into quark and antiquark pairs, then applying a Monte-Carlo version of instantaneous quark recombination, and finally subjecting remnant quarks (those which have not found a recombination partner) to Lund string fragmentation. When applied to parton showers from the PYTHIA Monte Carlo event generator, the final hadron spectra from our calculation compare quite well to PYTHIA jets that have been hadronized with the default Lund string fragmentation. Our new approach opens up the possibility to generalize hadronization to jets embedded in a quark gluon plasma.

preprint2015arXiv

Jet Hadronization via Recombination of Parton Showers in Vacuum and in Medium

We have studied the hadronization of jet parton showers based on the quark recombination model. This is achieved by letting gluons at the end of the perturbative shower evolution undergo a non-perturbative splitting into quark and antiquark pairs, then applying a Monte-Carlo version of instantaneous quark recombination, and finally subjecting remnant quarks (those which have not found a recombination partner) to Lund string fragmentation. When applied to parton showers from the PYTHIA Monte Carlo generator, the final hadron spectra from our calculation compare quite well to PYTHIA jets that have been hadronized with the default Lund string fragmentation. Modeling the quark gluon plasma produced in heavy ion collisions by a blast wave model, we have further studied medium effects on the hadronization of jet shower partons by also including their recombination with the thermal partons from the quark gluon plasma. We find that the latter leads to a significant enhancement of intermediate transverse momentum pions and protons at both RHIC and LHC. Our results thus suggest that medium modification of jet fragmentation provides a plausible explanation for the enhanced production of intermediate transverse momentum hadrons observed in experiments.

preprint2015arXiv

Light (anti-)nuclei production and flow in relativistic heavy-ion collisions

Using the coalescence model based on the phase-space distributions of protons, neutrons, Lambdas and their antiparticles from a multiphase transport (AMPT) model, we study the production of deuteron, triton, helium 3, hypertriton, hyperhelium 3 and their antinuclei in Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV. The resulting transverse momentum spectra, elliptic flows and coalescence parameters for these nuclei are presented and compared with available experimental data. We also show the constituent number scaled elliptic flows of these nuclei and discuss their implications.

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

Spinodal Instabilities of Baryon-Rich Quark-gluon Plasma in the PNJL Model

Using the Polyakov-Nambu-Jona-Lasinia (PNJL) model, we study the spinodal instability of a baryon-rich quark-gluon plasma in the linear response theory. We find that the spinodal unstable region in the temperature and density plane shrinks with increasing wave number of the unstable mode and is also reduced if the effect of Polyakov loop is not included. In the small wave number or long wavelength limit, the spinodal boundaries in both cases of with and without the Polyakov loop coincide with those determined from the isothermal spinodal instability in the thermodynamic approach. Also, the vector interactions among quarks is found to suppress unstable modes of all wave numbers. Moreover, the growth rate of unstable modes initially increases with the wave number but is reduced when the wave number becomes large. Including the collisional effect from quark scattering via the linearized Boltzmann equation, we further find that it decreases the growth rate of unstable modes of all wave numbers. Relevance of these results to relativistic heavy ion collisions is discussed.

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

Heavy quark correlations and the effective volume for quarkonia production

Using the Boltzmann transport approach, we study the effective volume of a correlated heavy quark pair in a partonic medium based on their collision rate. We find that the effective volume is finite and depends sensitively on the momentum of the heavy quark and the temperature of the medium. Generally, it increases linearly with time $t$ at the very beginning and the increase then becomes slower due to multiple scattering, and finally it increases linearly with respect to $t^{3/2}$. We further find that the colliding heavy quark pair has an effective temperature similar to that of the medium even though their initial transverse momentum spectra are far from thermal equilibrium.

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.

preprint2014arXiv

Probing isospin- and momentum-dependent nuclear effective interactions in neutron-rich matter

The single-particle potentials for nucleons and hyperons in neutron-rich matter generally depends on the density and isospin asymmetry of the medium as well as the momentum and isospin of the particle. It further depends on the temperature of the matter if the latter is in thermal equilibrium. We review here the extension of a Gogny-type isospin- and momentum-dependent interaction in several aspects made in recent years and their applications in studying intermediate-energy heavy ion collisions, thermal properties of asymmetric nuclear matter and properties of neutron stars. The importance of the isospin- and momentum-dependence of the single-particle potential, especially the momentum dependence of the isovector potential, is clearly revealed throughout these studies.

preprint2014arXiv

Spin Asymmetry of J/psi in Peripheral Pb+Pb Collisions at LHC

By generalizing the statistical model for particle production to the spin degree of freedom of initially produced J/psi, we study the spin projection J_y of J/psi perpendicular to the reaction plane in peripheral heavy ion collisions at the LHC energy that leads to a strong, albeit of short duration, magnetic field. We find that for J/psis produced directly from charm and anticharm quarks in the color singlet state, like that in the Color-Singlet Model, their yield in the presence of the magnetic field is larger for J_y=0 than for J_y=1 or -1. This leads to a spin asymmetry of finally produced J/ψeven after including their final-state scattering in the produced quark-gluon plasma.

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

Energy dependence of pion in-medium effects on π^-/π^+ ratio in heavy-ion collisions

Within the framework of a thermal model with its parameters fitted to the results from an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model, we study the pion in-medium effect on the charged-pion ratio in heavy-ion collisions at various energies. We find that due to the cancellation between the effects from pion-nucleon s-wave and p-wave interactions in nuclear medium, the π^-/π^+ ratio generally decreases after including the pion in-medium effect. The effect is larger at lower collision energies as a result of narrower pion spectral functions at lower temperatures.

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.

preprint2013arXiv

Shear viscosity of neutron-rich nucleonic matter near its liquid-gas phase transition

Within a relaxation time approach using free nucleon-nucleon cross sections modified by the in-medium nucleon masses that are determined from an isospin- and momentum-dependent effective nucleon-nucleon interaction, we investigate the specific shear viscosity ($η/s$) of neutron-rich nucleonic matter near its liquid-gas phase transition. It is found that as the nucleonic matter is heated at fixed pressure or compressed at fixed temperature, its specific shear viscosity shows a valley shape in the temperature or density dependence, with the minimum located at the boundary of the phase transition. Moreover, the value of $η/s$ drops suddenly at the first-order liquid-gas phase transition temperature, reaching as low as $4\sim5$ times the KSS bound of $\hbar/4π$. However, it varies smoothly for the second-order liquid-gas phase transition. Effects of the isospin degree of freedom and the nuclear symmetry energy on the value of $η/s$ are also discussed.

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.

preprint2012arXiv

Contributions of hyperon-hyperon scattering to subthreshold cascade production in heavy ion collisions

Using a gauged flavor SU(3)-invariant hadronic Lagrangian, we calculate the cross sections for the strangeness-exchange reactions YY to NΞ(Y=Λ, Σ) in the Born approximation. These cross sections are then used in the Relativistic Vlasov-Uehling-Uhlenbeck (RVUU) transport model to study Ξproduction in Ar+KCl collisions at incident energy of 1.76A GeV and impact parameter b=3.5 fm. We find that including the contributions of hyperon-hyperon scattering channels strongly enhances the yield of Ξ, leading to the abundance ratio Ξ^{-}/(Λ+Σ^{0})=3.38E-3, which is essentially consistent with the recently measured value of $(5.6 \pm 1.2_{-1.7}^{+1.8})\times 10^{-3}$ by the HADES collaboration at GSI.

preprint2012arXiv

Effects of hadronic potentials on elliptic flows in relativistic heavy ion collisions

Within the framework of a multiphase transport (AMPT) model that includes both initial partonic and final hadronic interactions, we show that including mean-field potentials in the hadronic phase leads to a splitting of the elliptic flows of particles and their antiparticles, providing thus a plausible explanation of the different elliptic flows between $p$ and ${\bar p}$, $K^+$ and $K^-$, and $π^+$ and $π^-$ observed in recent Beam Energy Scan (BES) program at the Relativistic Heavy-Ion Collider (RHIC).

preprint2012arXiv

Jet Fragmentation via Recombination of Parton Showers

We study hadron production in jets by applying quark recombination to jet shower partons. With the jet showers obtained from PYTHIA and augmented by additional non-perturbative effects, we compute hadron spectra in e+ + e-collisions at sqrt(s)=200 GeV. Including contributions from resonance decays, we find that the resulting transverse momentum spectra for pions, kaons, and protons reproduce reasonably those from the string fragmentation as implemented in PYTHIA.

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

Exotics from Heavy Ion Collisions

Discriminating hadronic molecular and multi-quark states is a long standing problem in hadronic physics. We propose here to utilize relativistic heavy ion collisions to resolve this problem, as exotic hadron yields are expected to be strongly affected by their structures. Using the coalescence model, we find that the exotic hadron yield relative to the statistical model result is typically an order of magnitude smaller for a compact multi-quark state, and larger by a factor of two or more for a loosely bound hadronic molecule. We further find that some of the newly proposed heavy exotic states could be produced and realistically measured at RHIC and LHC.

preprint2011arXiv

Higher-order anisotropic flows and dihadron correlations in Pb-Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV in a multiphase transport model

Using a multiphase transport model that includes both initial partonic and final hadronic scatterings, we have studied higher-order anisotropic flows as well as dihadron correlations as functions of pseudorapidity and azimuthal angular differences between trigger and associated particles in Pb-Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV. With parameters in the model determined previously from fitting the measured multiplicity density of mid-pseudorapidity charged particles in central collisions and their elliptic flow in mid-central collisions, the calculated higher-order anisotropic flows from the two-particle cumulant method reproduce approximately those measured by the ALICE Collaboration, except at small centralities where they are slightly overestimated. Similar to experimental results, the two-dimensional dihadron correlations at most central collisions show a ridge structure at the near side and a broad structure at the away side. The short- and long-range dihadron azimuthal correlations, corresponding to small and large pseudorapidity differences, respectively, are studied for triggering particles with different transverse momenta and are found to be qualitatively consistent with experimental results from the CMS Collaboration. The relation between the short-range and long-range dihadron correlations with that induced by back-to-back jet pairs produced from initial hard collisions is also discussed.

preprint2011arXiv

Imprints of Nuclear Symmetry Energy on Properties of Neutron Stars

Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$-mode frequency and damping time of gravitational wave emission of neutron stars.

preprint2011arXiv

Multi-quark hadrons from Heavy Ion Collisions

Identifying hadronic molecular states and/or hadrons with multi-quark components either with or without exotic quantum numbers is a long standing challenge in hadronic physics. We suggest that studying the production of these hadrons in relativistic heavy ion collisions offer a promising resolution to this problem as yields of exotic hadrons are expected to be strongly affected by their structures. Using the coalescence model for hadron production, we find that compared to the case of a non-exotic hadron with normal quark numbers, the yield of an exotic hadron is typically an order of magnitude smaller when it is a compact multi-quark state and a factor of two or more larger when it is a loosely bound hadronic molecule. We further find that due to the appreciable numbers of charm and bottom quarks produced in heavy ion collisions at RHIC and even larger numbers expected at LHC, some of the newly proposed heavy exotic states could be produced and realistically measured in these experiments.

preprint2011arXiv

Pb-Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV in a multiphase transport model

The multiplicity and elliptic flow of charged particles produced in Pb-Pb collisions at center of mass energy $\sqrt{s_{NN}}=2.76$ TeV from the Large Hadron Collider are studied in a multiphase transport (AMPT) model. With the standard parameters in the HIJING model, which is used as initial conditions for subsequent partonic and hadronic scatterings in the AMPT model, the resulting multiplicity of final charged particles at mid-pseudorapidity is consistent with the experimental data measured by the ALICE Collaboration. This value is, however, increased by about 25% if the final-state partonic and hadronic scatterings are turned off. Because of final-state scatterings, particular those among partons, the final elliptic flow of charged hadrons is also consistent with the ALICE data if a smaller but more isotropic parton scattering cross section than previously used in the AMPT model for describing the charged hadron elliptic flow in heavy ion collisions at the Relativistic Heavy Ion Collider is used. The resulting transverse momentum spectra of charged particles as well as the centrality dependence of their multiplicity density and the elliptic flow are also in reasonable agreement with the ALICE data. Furthermore, the multiplicities, transverse momentum spectra and elliptic flows of identified hadrons such as protons, kaons and pions are predicted.

preprint2011arXiv

Studying Exotic Hadrons in Heavy Ion Collisions

We investigate the possibilities of using measurements in present and future experiments on heavy ion collisions to answer some longstanding problems in hadronic physics, namely identifying hadronic molecular states and exotic hadrons with multiquark components. The yields of a selected set of exotic hadron candidates in relativistic heavy ion collisions are discussed in the coalescence model in comparison with the statistical model. We find that the yield of a hadron is typically an order of magnitude smaller when it is a compact multiquark state, compared to that of an excited hadronic state with normal quark numbers. We also find that some loosely bound hadronic molecules are formed more abundantly than the statistical model prediction by a factor of two or more. Moreover, due to the significant numbers of charm and bottom quarks produced at RHIC and even larger numbers expected at LHC, some of the proposed heavy exotic hadrons could be produced with sufficient abundance for detection, making it possible to study these new exotic hadrons in heavy ion collisions.

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 effect of triangular flow on di-hadron azimuthal correlations in relativistic heavy ion collisions

Using the AMPT model for relativistic heavy ion collisions, we have studied the di-hadron azimuthal angular correlations triggered by emitted jets in Au+Au collisions at center of mass energy $\sqrt{s_{NN}}=200$ GeV and impact parameter $b=8$ fm. A double-peak structure for the associated particles at the away side of trigger particles is obtained after subtracting background correlations due to the elliptic flow. Both the near-side peak and the away-side double peaks in the azimuthal angular correlations are, however, significantly suppressed (enhanced) in events of small (large) triangular flow, which are present as a result of fluctuations in the initial collision geometry. After subtraction of background correlations due to the triangular flow, the away-side double peaks change into a single peak with broad shoulders on both sides. The away side of the di-hadron correlations becomes essentially a single peak after further subtraction of higher-order flows.

preprint2011arXiv

Transport Model Studies of the Baryon-Rich Quark-Gluon Plasma formed in Heavy Ion Collisions

Heavy ion collisions in the low energy run at Relativistic Heavy Ion Collider (RHIC) and future Facility for Antiproton and Ion Research (FAIR) in Germany are expected to produce a quark-gluon plasma that has a finite baryon chemical potential, allowing thus the possibility to study the location of the critical endpoint in the QCD phase diagram. In this talk, using a multiphase transport model, that includes interactions in both initial partonic and final hadronic matters and the transition between these two phases of matter, we discuss the effects of partonic interactions on observables such as the elliptic flow that have played essential roles in studying the properties of the net baryon free quark-gluon plasma produced in heavy ion collisions at higher energies at RHIC. Also, we study the effect of density fluctuations due to a first-order transition between the quark-gluon plasma and hadronic matter on fluctuations of hadron mean transverse momentum and produced deuteron number as well as on two-pion correlations. We further discuss the possibility of studying the mechanism of charm energy loss in the baryon-rich quark-gluon plasma and the properties of phi mesons in hot-dense matter produced in these collisions.

preprint2011arXiv

Triangular flow in heavy ion collisions in a multiphase transport model

We have obtained a new set of parameters in a multiphase transport (AMPT) model that are able to describe both the charged particle multiplicity density and elliptic flow measured in Au+Au collisions at center of mass energy $\sqrt{s_{NN}}=200$ GeV at the Relativistic Heavy Ion Collider (RHIC), although they still give somewhat softer transverse momentum spectra. We then use the model to predict the triangular flow due to fluctuations in the initial collision geometry and study its effect relative to those from other harmonic components of anisotropic flows on the di-hadron azimuthal correlations in both central and mid-central collisions.

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

Analytical relations between nuclear symmetry energy and single-nucleon potentials in isospin asymmetric nuclear matter

Using the Hugenholtz-Van Hove theorem, we derive general expressions for the quadratic and quartic symmetry energies in terms of single-nucleon potentials in isospin asymmetric nuclear matter. These analytical relations are useful for gaining deeper insights into the microscopic origins of the uncertainties in our knowledge on nuclear symmetry energies especially at supra-saturation densities. As examples, the formalism is applied to two model single-nucleon potentials that are widely used in transport model simulations of heavy-ion reactions.

preprint2010arXiv

Density matrix expansion for the MDI interaction

By assuming that the isospin- and momentum-dependent MDI interaction has a form similar to the Gogny-like effective two-body interaction with a Yukawa finite-range term and the momentum dependence only originates from the finite-range exchange interaction, we determine its parameters by comparing the predicted potential energy density functional in uniform nuclear matter with what has been usually given and used extensively in transport models for studying isospin effects in intermediate-energy heavy-ion collisions as well as in investigating the properties of hot asymmetric nuclear matter and neutron star matter. We then use the density matrix expansion to derive from the resulting finite-range exchange interaction an effective Skyrme-like zero-range interaction with density-dependent parameters. As an application, we study the transition density and pressure at the inner edge of neutron star crusts using the stability conditions derived from the linearized Vlasov equation for the neutron star matter.

preprint2010arXiv

Density slope of the nuclear symmetry energy from the neutron skin thickness of heavy nuclei

Expressing explicitly the parameters of the standard Skyrme interaction in terms of the macroscopic properties of asymmetric nuclear matter, we show in the Skyrme-Hartree-Fock approach that unambiguous correlations exist between observables of finite nuclei and nuclear matter properties. We find that existing data on neutron skin thickness $Δr_{np}$ of Sn isotopes give an important constraint on the symmetry energy $E_{sym}({ρ_{0}})$ and its density slope $L$ at saturation density ${ρ_{0}}$. Combining these constraints with those from recent analyses of isospin diffusion and double neutron/proton ratio in heavy-ion collisions at intermediate energies leads to a more stringent limit on $L$ approximately independent of $E_{sym}({ρ_{0}})$. The implication of these new constraints on the $Δr_{np}$ of $^{208}$Pb as well as the core-crust transition density and pressure in neutron stars is discussed.

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.

preprint2010arXiv

Energy release from hadron-quark phase transition in neutron stars and the axial $w$-mode of gravitational waves

Describing the hyperonic and quark phases of neutron stars with an isospin- and momentum-dependent effective interaction for the baryon octet and the MIT bag model, respectively, and using the Gibbs conditions to construct the mixed phase, we study the energy release due to the hadron-quark phase transition. Moreover, the frequency and damping time of the first axial $w$-mode of gravitational waves are studied for both hyperonic and hybrid stars. We find that the energy release is much more sensitive to the bag constant than the density dependence of the nuclear symmetry energy. Also, the frequency of the $w$-mode is found to be significantly different with or without the hadron-quark phase transition and depends strongly on the value of the bag constant. Effects of the density dependence of the nuclear symmetry energy become, however, important for large values of the bag constant that lead to higher hadron-quark transition densities.

preprint2010arXiv

Isospin- and momentum-dependent effective interactions for the baryon octet and the properties of hybrid stars

The isospin- and momentum-dependent MDI interaction, which has been extensively used in intermediate-energy heavy-ion reactions to study the properties of asymmetric nuclear matter, is extended to include the nucleon-hyperon and hyperon-hyperon interactions by assuming same density, momentum and isospin dependence as for the nucleon-nucleon interaction. The parameters in these interactions are determined from the empirical hyperon single-particle potentials in symmetric nuclear matter at saturation density. The extended MDI interaction is then used to study in the mean-field approximation the equation of state of hypernuclear matter and also the properties of hybrid stars by including the phase transition from the hypernuclear matter to the quark matter at high densities. In particular, the effects of attractive and repulsive $Σ$N interactions and different values of symmetry energies on the hybrid star properties are investigated.

preprint2010arXiv

Isospin-dependent pion in-medium effects on charged pion ratio in heavy ion collisions

Using results from the chiral perturbation theory for the s-wave interaction and the delta-resonance model for the p-wave interaction of pions with nucleons, we have evaluated the spectral functions of pions in asymmetric nuclear matter with unequal proton and neutron densities. We find that in hot dense neutron-rich matter the strength of the spectral function of positively charged pion at low energies is somewhat larger than that of negatively charged pion. In a thermal model, this isospin-dependent effect slightly reduces the ratio of negatively charged to positively charged pions that are produced in heavy ion collisions induced by radioactive beams. Relevance of our results to the determination of the nuclear symmetry energy from measured ratio of negatively to positively charged pions produced in heavy ion collisions is discussed.

preprint2010arXiv

Medium effects on charged pion ratio in heavy ion collisions

We have recently studied in the delta-resonance--nucleon-hole model the dependence of the pion spectral function in hot dense asymmetric nuclear matter on the charge of the pion due to the pion p-wave interaction in nuclear medium. In a thermal model, this isospin-dependent effect enhances the ratio of negatively charged to positively charged pions in neutron-rich nuclear matter, and the effect is comparable to that due to the uncertainties in the theoretically predicted stiffness of nuclear symmetry energy at high densities. This effect is, however, reversed if we also take into account the s-wave interaction of the pion in nuclear medium as given by chiral perturbation theory, resulting instead in a slightly reduced ratio of negatively charged to positively charged pions. Relevance of our results to the determination of the nuclear symmetry energy from the ratio of negatively to positively charged pions produced in heavy ion collisions is discussed.

preprint2010arXiv

Transition density and pressure in hot neutron stars

Using the momentum-dependent MDI effective interaction for nucleons, we have studied the transition density and pressure at the boundary between the inner crust and liquid core of hot neutron stars. We find that their values are larger in neutrino-trapped neutron stars than in neutrino-free neutron stars. Furthermore, both are found to decrease with increasing temperature of a neutron star as well as increasing slope parameter of the nuclear symmetry energy, except that the transition pressure in neutrino-trapped neutron stars for the case of small symmetry energy slope parameter first increases and then decreases with increasing temperature. We have also studied the effect of the nuclear symmetry energy on the critical temperature above which the inner crust in a hot neutron star disappears and found that with increasing value of the symmetry energy slope parameter, the critical temperature decreases slightly in neutrino-trapped neutron stars but first decreases and then increases in neutrino-free neutron stars.

preprint2009arXiv

Incompressibility of asymmetric nuclear matter

The incompressibility $K_sat(δ)$ of isospin asymmetric nuclear matter at its saturation density. Our results show that in the expansion of $K_sat(δ)$ in powers of isospin asymmetry $δ$, i.e., $K_sat(δ)$=K_{0}+K_{sat,2}δ^{2}+K_{sat,4}δ^{4}+O(δ^{6})$, the magnitude of the 4th-order K_{sat,4} parameter is generally small. The 2nd-order K_{sat,2} parameter thus essentially characterizes the isospin dependence of the incompressibility of asymmetric nuclear matter at saturation density. Furthermore, the K_{sat,2} can be expressed as K_{sat,2}=K_{sym}-6L-J_{0}/{K_{0}L in terms of the slope parameter $L$ and the curvature parameter $K_{\mathrm{sym}}$ of the symmetry energy and the third-order derivative parameter $J_0$ of the energy of symmetric nuclear matter at saturation density, and we find the higher order $J_0$ contribution to K_{sat,2} generally cannot be neglected. Also, we have found a linear correlation between K_{sym} and $L$ as well as between $J_{0}/K_{0}$ and $K_{0}$. Using these correlations together with the empirical constraints on $K_{0}$ and $L$, the nuclear symmetry energy $E_sym(ρ_{0})$ at normal nuclear density, and the nucleon effective mass, we have obtained an estimated value of K_{sat,2}=-370 +- 120 MeV for the 2nd-order parameter in the isospin asymmetry expansion of the incompressibility of asymmetric nuclear matter at its saturation density.

preprint2009arXiv

Transition Density and Pressure at the Inner Edge of Neutron Star Crusts

Using the nuclear symmetry energy that has been recently constrained by the isospin diffusion data in intermediate-energy heavy ion collisions, we have studied the transition density and pressure at the inner edge of neutron star crusts, and they are found to be 0.040 fm$^{-3}$ $\leq ρ_{t}\leq 0.065$ fm$^{-3}$ and 0.01 MeV/fm$^{3}$ $\leq P_{t}\leq 0.26$ MeV/fm$^{3}$, respectively, in both the dynamical and thermodynamical approaches. We have also found that the widely used parabolic approximation to the equation of state of asymmetric nuclear matter gives significantly higher values of core-crust transition density and pressure, especially for stiff symmetry energies. With these newly determined transition density and pressure, we have obtained an improved relation between the mass and radius of neutron stars.

preprint2008arXiv

Recent Progress and New Challenges in Isospin Physics with Heavy-Ion Reactions

The ultimate goal of studying isospin physics via heavy-ion reactions with neutron-rich, stable and/or radioactive nuclei is to explore the isospin dependence of in-medium nuclear effective interactions and the equation of state of neutron-rich nuclear matter, particularly the isospin-dependent term in the equation of state, i.e., the density dependence of the symmetry energy. Because of its great importance for understanding many phenomena in both nuclear physics and astrophysics, the study of the density dependence of the nuclear symmetry energy has been the main focus of the intermediate-energy heavy-ion physics community during the last decade, and significant progress has been achieved both experimentally and theoretically. In particular, a number of phenomena or observables have been identified as sensitive probes to the density dependence of the nuclear symmetry energy. Experimental studies have confirmed some of these interesting isospin-dependent effects and allowed us to constrain relatively stringently the symmetry energy at sub-saturation densities. The impacts of this constrained density dependence of the symmetry energy on the properties of neutron stars have also been studied, and they were found to be very useful for the astrophysical community. With new opportunities provided by the various radioactive beam facilities being constructed around the world, the study of isospin physics is expected to remain one of the forefront research areas in nuclear physics. In this report, we review the major progress achieved during the last decade in isospin physics with heavy ion reactions and discuss future challenges to the most important issues in this field.

preprint2005arXiv

A Transport Model for Nuclear Reactions Induced by Radioactive Beams

Major ingredients of an isospin and momentum dependent transport model for nuclear reactions induced by radioactive beams are outlined. Within the IBUU04 version of this model we study several experimental probes of the equation of state of neutron-rich matter, especially the density dependence of the nuclear symmetry energy. Comparing with the recent experimental data from NSCL/MSU on isospin diffusion, we found a nuclear symmetry energy of $% E_{sym}(ρ)\approx 31.6(ρ/ρ_{0})^{1.05}$ at subnormal densities. Predictions on several observables sensitive to the density dependence of the symmetry energy at supranormal densities accessible at GSI and the planned Rare Isotope Accelerator (RIA) are also made.

preprint2003arXiv

Charm production from proton-proton collisions

We evaluate the cross sections for charmed hadron production from proton-proton reactions $pp\to\bar D^0pΛ_c^+$ and $pp\to\bar D^{*0}pΛ_c^+$ using a hadronic Lagrangian. With empirical coupling constants and cutoff parameters in the form factors, sum of their cross sections at center-of-mass energy of 11.5 GeV is about 1 $μ$b and is comparable to measured inclusive cross section for charmed hadron production from proton-proton reactions. The cross section decreases to about 1 nb at 40 MeV above threshold.

preprint1997arXiv

Isospin Physics in Heavy-Ion Collisions at Intermediate Energies

In nuclear collisions induced by stable or radioactive neutron-rich nuclei a transient state of nuclear matter with an appreciable isospin asymmetry as well as thermal and compressional excitation can be created. This offers the possibility to study the properties of nuclear matter in the region between symmetric nuclear matter and pure neutron matter. In this review, we discuss recent theoretical studies of the equation of state of isospin-asymmetric nuclear matter and its relations to the properties of neutron stars and radioactive nuclei. Chemical and mechanical instabilities as well as the liquid-gas phase transition in asymmetric nuclear matter are investigated. The in-medium nucleon-nucleon cross sections at different isospin states are reviewed as they affect significantly the dynamics of heavy ion collisions induced by radioactive beams. We then discuss an isospin-dependent transport model, which includes different mean-field potentials and cross sections for the proton and neutron, and its application to these reactions. Furthermore, we review the comparisons between theoretical predictions and available experimental data. In particular, we discuss the study of nuclear stopping in terms of isospin equilibration, the dependence of nuclear collective flow and balance energy on the isospin-dependent nuclear equation of state and cross sections, the isospin dependence of total nuclear reaction cross sections, and the role of isospin in preequilibrium nucleon emissions and subthreshold pion production.

preprint1996arXiv

Excitation functions in central Au+Au collisions from SIS/GSI to AGS/Brookhaven

Using the relativistic transport model (ART), we predict the energy dependence of the stopping power, maximum baryon and energy densities, the population of resonance matter as well as the strength of the transverse and radial flow for central Au+Au reactions at beam momentum from 2 to 12 GeV/c available at Brookhaven's AGS. The maximum baryon and energy densities are further compared to the predictions of relativistic hydrodynamics assuming the formation of shock waves. We also discuss the Fermi-Landau scaling of the pion multiplicity in these reactions.

preprint1996arXiv

Pion Transparency in 500 MeV C(pi,pi') Reactions?

The question whether there is a pion transparency in 500 MeV $C(π, π^{'})$ scatterings is studied using a semiclassical, hadronic transport model. The double differential cross sections of this reaction measured at LAMPF can be largely accounted for, if one uses energy-dependent, anisotropic angular distributions which are fitted to pion-nucleon scattering data for the decay of $Δ(1232)$ and $N^*(1440)$ resonances. The remaining discrepancy between the data and the calculation sets a limit on effects of more exotic processes.