No Arabic abstract
We study cluster and hypernuclei production in heavy-ion collisions at relativistic energies employing the Parton-Hadron-Quantum-Molecular-Dynamics (PHQMD) approach, a microscopic n-body transport model based on the QMD propagation of the baryonic degrees of freedom with density dependent 2-body potential interactions. All other ingredients of PHQMD, including the collision integral and the treatment of the quark-gluon plasma (QGP) phase, are adopted from the Parton-Hadron-String Dynamics (PHSD) approach. In PHQMD the cluster formation occurs dynamically, caused by the interactions. The clusters are recognized by the Minimum Spanning Tree (MST) algorithm. We present the PHQMD results for cluster and hypernuclei formation in comparison with the available experimental data at AGS, SPS, RHIC-BES and RHIC fixed target energies. We also provide predictions on cluster production for the upcoming FAIR and NICA experiments. PHQMD allows to study the time evolution of formed clusters and the origin of their production, which helps to understand how such weakly bound objects are formed and survive in the rather dense and hot environment created in heavy-ion collisions. It offers therefore an explanation of the ice in the fire puzzle.
Cluster and hypernuclei production in heavy-ion collisions is presently under active experimental and theoretical investigation. Since clusters are weekly bound objects, their production is very sensitive to the dynamical evolution of the system and its interactions. The theoretical description of cluster formation is related to the n-body problem. Here we present the novel n-body dynamical transport approach PHQMD (Parton-Hadron-Quantum-Molecular Dynamics) which is designed to provide a microscopic description of nuclear cluster and hypernucleus formation as well as of general particle production in heavy-ion reactions at relativistic energies. In difference to the coalescence or statistical models, often used for the cluster formation, in PHQMD clusters are formed dynamically due to the interactions between baryons described on a basis of Quantum Molecular Dynamics (QMD)which allows to propagate the n-body Wigner density and n-body correlations in phase-space, essential for the cluster formation. The clusters are identified by the MST (Minimum Spanning Tree) or the SACA (Simulated Annealing Cluster Algorithm) algorithm which finds the most bound configuration of nucleons and clusters. Collisions among hadrons as well as Quark-Gluon-Plasma formation and parton dynamics in PHQMD are treated in the same way as in the established PHSD (Parton-Hadron-String Dynamics)transport approach. In order to verify our approach with respect to the general dynamics we present here the first PHQMD results for general bulk observables such as rapidity distributions and transverse mass spectra for hadrons ($pi, K, bar K, p, bar p, Lambda, bar Lambda$) from SIS to RHIC energies. We find a good description of the bulk dynamics which allows us to proceed with the results on cluster production, including hypernuclei.
We develop a combined hydro-kinetic approach which incorporates a hydrodynamical expansion of the systems formed in textit{A}+textit{A} collisions and their dynamical decoupling described by escape probabilities. The method corresponds to a generalized relaxation time ($tau_{text{rel}}$) approximation for the Boltzmann equation applied to inhomogeneous expanding systems; at small $tau_{text{rel}}$ it also allows one to catch the viscous effects in hadronic component - hadron-resonance gas. We demonstrate how the approximation of sudden freeze-out can be obtained within this dynamical picture of continuous emission and find that hypersurfaces, corresponding to a sharp freeze-out limit, are momentum dependent. The pion $m_{T}$ spectra are computed in the developed hydro-kinetic model, and compared with those obtained from ideal hydrodynamics with the Cooper-Frye isothermal prescription. Our results indicate that there does not exist a universal freeze-out temperature for pions with different momenta, and support an earlier decoupling of higher $p_{T}$ particles. By performing numerical simulations for various initial conditions and equations of state we identify several characteristic features of the bulk QCD matter evolution preferred in view of the current analysis of heavy ion collisions at RHIC energies.
We review integrated dynamical approaches to describe heavy ion reaction as a whole at ultrarelativistic energies. Since final observables result from all the history of the reaction, it is important to describe all the stages of the reaction to obtain the properties of the quark gluon plasma from experimental data. As an example of these approaches, we develop an integrated dynamical model, which is composed of a fully (3+1) dimensional ideal hydrodynamic model with the state-of-the-art equation of state based on lattice QCD, and subsequent hadronic cascade in the late stage. Initial conditions are obtained employing Monte Car
The experimental data on hadron yields and ratios in central lead-lead and gold-gold collisions at 158 AGeV/$c$ (SPS) and $sqrt{s} = 130$ AGeV (RHIC), respectively, are analysed within a two-source statistical model of an ideal hadron gas. A comparison with the standard thermal model is given. The two sources, which can reach the chemical and thermal equilibrium separately and may have different temperatures, particle and strangeness densities, and other thermodynamic characteristics, represent the expanding system of colliding heavy ions, where the hot central fireball is embedded in a larger but cooler fireball. The volume of the central source increases with rising bombarding energy. Results of the two-source model fit to RHIC experimental data at midrapidity coincide with the results of the one-source thermal model fit, indicating the formation of an extended fireball, which is three times larger than the corresponding core at SPS.
We study charm production in ultra-relativistic heavy-ion collisions by using the Parton-Hadron-String Dynamics (PHSD) transport approach. The initial charm quarks are produced by the PYTHIA event generator tuned to fit the transverse momentum spectrum and rapidity distribution of charm quarks from Fixed-Order Next-to-Leading Logarithm (FONLL) calculations. The produced charm quarks scatter in the quark-gluon plasma (QGP) with the off-shell partons whose masses and widths are given by the Dynamical Quasi-Particle Model (DQPM), which reproduces the lattice QCD equation-of-state in thermal equilibrium. The relevant cross sections are calculated in a consistent way by employing the effective propagators and couplings from the DQPM. Close to the critical energy density of the phase transition, the charm quarks are hadronized into $D$ mesons through coalescence and/or fragmentation. The hadronized $D$ mesons then interact with the various hadrons in the hadronic phase with cross sections calculated in an effective lagrangian approach with heavy-quark spin symmetry. The nuclear modification factor $R_{AA}$ and the elliptic flow $v_2$ of $D^0$ mesons from PHSD are compared with the experimental data from the STAR Collaboration for Au+Au collisions at $sqrt{s_{NN}}$ =200 GeV and to the ALICE data for Pb+Pb collisions at $sqrt{s_{NN}}$ =2.76 TeV. We find that in the PHSD the energy loss of $D$ mesons at high $p_T$ can be dominantly attributed to partonic scattering while the actual shape of $R_{AA}$ versus $p_T$ reflects the heavy-quark hadronization scenario, i.e. coalescence versus fragmentation. Also the hadronic rescattering is important for the $R_{AA}$ at low $p_T$ and enhances the $D$-meson elliptic flow $v_2$.