From Full stopping to transparency in a holographic model of heavy ion collisions

We numerically simulate planar shock wave collisions in anti-de Sitter space as a model for heavy ion collisions of large nuclei. We uncover a crossover between two different dynamical regimes as a function of the collision energy. At low energies the shocks first stop and then explode in a manner approximately described by hydrodynamics, in close similarity with the Landau model. At high energies the receding fragments move outwards at the speed of light, with a region of negative energy density and negative longitudinal pressure trailing behind them. The rapidity distribution of the energy density at late times around midrapidity is not approximately boost invariant but Gaussian, albeit with a width that increases with the collision energy.

Initial state of ultrarelativistic heavy ion collision

A model for energy, pressure, and flow velocity distributions at the beginning of ultrarelativistic heavy ion collisions is presented, which can be used as an initial condition for hydrodynamic calculations. Our model takes into account baryon recoil for both target and projectile, arising from the acceleration of partons in an effective field F mu nu produced in the collision. The typical field strength (string tension) for RHIC energies is about 512 GeV/fm, which allows us to talk about string ropes. The results show that a quark-gluon plasma forms a tilted disk, such that the direction of the largest pressure gradient stays in the reaction plane, but deviates from both the beam and the usual transverse flow directions. Such initial conditions may lead to the creation of antiflow or third flow component [L. P. Csernai and D. Rhrich, Phys. Rev. Lett. B 458, 454 (1999)].

Statistical hadronization of supercooled quark-gluon plasma

The fast simultaneous hadronization and chemical freeze-out of supercooled quark-gluon plasma, created in relativistic heavy ion collisions, can lead to the reheating of the expanding matter and to the change in a collective flow profile. We use the assumption of statistical nature of the hadronization process, and study quantitatively the freeze-out in the framework of hydrodynamical Bjorken model with different simple quark-gluon plasma equations of state.

Matching stages of heavy-ion collision models

Heavy-ion reactions and other collective dynamical processes are frequently described by different theoretical approaches for the different stages of the process, like initial equilibration stage, intermediate locally equilibrated fluid dynamical stage, and final freeze-out stage. For the last stage, the best known is the Cooper-Frye description used to generate the phase space distribution of emitted, noninteracting particles from a fluid dynamical expansion or explosion, assuming a final ideal gas distribution, or (less frequently) an out-of-equilibrium distribution. In this work we do not want to replace the Cooper-Frye description, but rather clarify the ways of using it and how to choose the parameters of the distribution and, eventually, how to choose the form of the phase space distribution used in the Cooper-Frye formula. Moreover, the Cooper-Frye formula is used in connection with the freeze-out problem, while the discussion of transition between different stages of the collision is applicable to other transitions also. More recently, hadronization and molecular dynamics models have been matched to the end of a fluid dynamical stage to describe hadronization and freeze-out. The stages of the model description can be matched to each other on space-time hypersurfaces (just like through the frequently used freeze-out hypersurface). This work presents a generalized description of how to match the stages of the description of a reaction to each other, extending the methodology used at freeze-out, in simple covariant form which is easily applicable in its simplest version for most applications.

Covariant description of kinetic freeze-out through a finite spacelike layer

The problem of freeze-out (FO) in relativistic heavy-ion reactions is addressed. We develop and analyze an idealized one-dimensional model of FO in a finite layer, based on the covariant FO probability. The resulting post FO phase-space distributions are discussed for different FO probabilities and layer thicknesses.

New method for measuring longitudinal fuctuations and directed flow in ultrarelativistic heavy ion reactions

It has been shown in recent ALICE@LHC measurements that the odd flow harmonics, in particular, a directed flow v1, occurred to be weak and dominated by random fluctuations. In this work we propose a new method, which makes the measurements more sensitive to the flow patterns showing global collective symmetries. We demonstrate how the longitudinal center of mass rapidity fluctuations can be identified, and then the collective flow analysis can be performed in the event-by-event center of mass frame. Such a method can be very effective in separating the flow patterns originating from random fluctuations, and the flow patterns originating from the global symmetry of the initial state.

Antikaon potential in hot and dense matter

The antikaon optical potential in hot and dense nuclear matter is studied within the framework of a coupled-channel self-consistent calculation taking, as bare meson-baryon interaction, the meson-exchange potential of the Jlich group. Typical conditions found in heavy-ion collisions at GSI are explored. As in the case of zero temperature, the angular momentum components larger than L=0 contribute significantly to the finite temperature antikaon optical potential at finite momentum. It is found that the particular treatment of the medium effects has a strong influence on the behavior of the antikaon potential with temperature. Our self-consistent model, in which antikaons and pions are dressed in the medium, gives a moderately temperature dependent antikaon potential which remains attractive at GSI temperatures, contrary to what one finds if only nuclear Pauli blocking effects are included.

Holographic isotropization linearized

The holographic isotropization of a highly anisotropic, homogeneous, strongly coupled, non-Abelian plasma was simplified in ref. [1] by linearizing Einstein"s equations around the final, equilibrium state. This approximation reproduces the expectation value of the boundary stress tensor with a 20% accuracy. Here we elaborate on these results and extend them to observables that are directly sensitive to the bulk interior, focusing for simplicity on the entropy production on the event horizon. We also consider next-to-leading-order corrections and show that the leading terms alone provide a better description of the isotropization process for the states that are furthest from equilibrium.

Pauli distorted double folded potential

The enhancement in the production of even-Z nuclei observed in nuclear fission has also been observed in fragments produced from heavy ion collsions. Beams of 40Ar, 40Cl, and 40Ca at 25 MeV/nucleon were impinged on 58Fe and 58Ni targets. The resulting fragments were detected using the MSU 4pi detector array, which had additional silicon detectors for better isotopic resolution. Comparison of the ratios of yields for each element showed enhancement of even-Z fragment production. The enhancement was more pronounced for reactions with a greater difference in the N/Z of the compound system. However, this effect was less for systems that were more neutron rich. The average N/Z for fragments also displayed an odd-even effect with a lower average N/Z for the even-Z fragments. This is related to the greater availability of neutron-poor isotopes for even-Z nuclei

Double folding with a density-dependent effective interaction and it s analytical approximation

The real part of the optical potential for heavy ion elastic scattering is obtained by double folding of the nuclear densities with a density-dependent nucleon-nucleon effective interaction which was successful in describing the binding, size, and nucleon separation energies in spherical nuclei. A simple analytical form is found to differ from the resulting potential considerably less than 1% all through the important region. This analytical potential is used so that only few points of the folding need to be computed. With an imaginary part of the Woods-Saxon type, this potential predicts the elastic scattering angular distribution in very good agreement with experimental data, and little renormalization (unity in most cases) is needed.

Cherenkov mesons as in-medium quark energy loss

We recently showed that a heavy quark moving su ciently fast through a quark-gluon plasma may lose energy by Cherenkov-radiating mesons [1]. Here we review our previous holographic calculation of the energy loss in N = 4 Super Yang-Mills and extend it to longitudinal vector mesons and scalar mesons. We also discuss phenomenological implications for heavy-ion collision experiments. Although the Cherenkov energy loss is an O(1=Nc) effect, a ballpark estimate yields a value of dE/dx for Nc = 3 which is comparable to that of other mechanisms.

Cherenkov mesons as in-medium quark energy loss

We recently showed that a heavy quark moving su ciently fast through a quark-gluon plasma may lose energy by Cherenkov-radiating mesons [1]. Here we review our previous holographic calculation of the energy loss in N = 4 Super Yang-Mills and extend it to longitudinal vector mesons and scalar mesons. We also discuss phenomenological implications for heavy-ion collision experiments. Although the Cherenkov energy loss is an O(1=Nc) effect, a ballpark estimate yields a value of dE/dx for Nc = 3 which is comparable to that of other mechanisms.

Cherenkov mesons as in-medium quark energy loss

We recently showed that a heavy quark moving sufficiently fast through a quark-gluon plasma may lose energy by Cherenkov-radiating mesons [1]. Here we review our previous holographic calculation of the energy loss in N=4 Super Yang-Mills and extend it to longitudinal vector mesons and scalar mesons. We also discuss phenomenological implications for heavy-ion collision experiments. Although the Cherenkov energy loss is an O(1/Nc) effect, a ballpark estimate yields a value of dE/dx for Nc=3 which is comparable to that of other mechanisms.