Researcher profile

Che Ming Ko

Che Ming Ko contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 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.