# Hybrid approaches to heavy ion collisions and future perspectives

###### Abstract

We present the current status of hybrid approaches to describe heavy ion collisions and their future challenges and perspectives. First we present a hybrid model combining a Boltzmann transport model of hadronic degrees of freedom in the initial and final state with an optional hydrodynamic evolution during the dense and hot phase. Second, we present a recent extension of the hydrodynamical model to include fluctuations near the phase transition by coupling a chiral field to the hydrodynamic evolution.

-titleHot and Cold Baryonic Matter – HCBM 2010

## 1 Introduction

Heavy ion collisions are an excellent tool to study the properties of dense and hot nuclear matter Scherer:1999qq (). One of the most interesting features to explore in heavy ion collisions is the phase diagram of quantum chromodynamics (QCD). At high temperatures and large baryonic densities the degrees of freedom of strongly interacting matter are supposed to be partonic. There are indications that this state of matter, the quark-gluon plasma, was created in heavy ion collisions at RHIC and behaves as an almost ideal fluid Adams:2005dq (); Back:2004je (); Arsene:2004fa (); Adcox:2004mh (). It is one of the major challenges to learn about the quark-gluon plasma (QGP) from the final hadronic state observed in the detectors. The phase transition of QCD is twofold: besides the transition from deconfined partonic to confined hadronic matter, the chiral symmetry gets spontaneously broken at lower temperatures while the symmetry is restored in the high temperature phase. Lattice QCD studies have revealed that these two aspects form one crossover transition at least for vanishing baryochemical potential at critical temperatures given between MeV Borsanyi:2010bp (); Soldner:2010xk (). Studies of effective low energy models of QCD show that the phase transition is a first order transition at high baryon densities and low temperatures Scavenius:2000qd (). Consequently the first order phase transition line terminates at a critical end point (CEP) in the plane of the phase diagram.

A broad variety of experimental signatures is proposed to study the QCD phase diagram in experiment. One very famous signature is the ’horn’ structure in the excitation function for the kaon to pion ratio :2007fe (). It was attributed to the onset of deconfinement at SPS energies Gazdzicki:1997hm (). Other proposed signatures cover aspects of the chiral phase transition, a suggested critical point Stephanov:1998dy (); Stephanov:1999zu () and the first order phase transition Mishustin:1998eq (). It is, however, not clear how much of a potential signal is developed in a realistic scenario of a heavy ion collision. The fast dynamics present during the creation and expansion of the fireball of strongly interacting matter makes it difficult to defer results from thermodynamic calculations directly. A thorough theoretical understanding of phase transitions in the environment of a heavy ion collision is necessary to make profound predictions for the running low energy RHIC program and the upcoming FAIR project.

We present results from a hybrid model approach to heavy ion collisions Petersen:2008dd (); Andrade:2006yh (); Teaney:2001av (); Hirano:2007gc () and give an outlook on chiral fluid dynamic models which allow for the explicit propagation of critical fluctuations through the phase transition.

## 2 The hybrid approach

The hybrid approach describes a heavy ion collision in three steps, the inital state, the hot and dense stage with the phase transition and the final interactions after hadronization. By fixing the inital state and the freeze-out of particles from the hydrodynamically propagated densities one can focus on the effects of the phase transition and viscosity in the intermediate state. The model presented here is an extention of the standard transport UrQMD model by the optional use of an intermediate hydrodynamic evolution Petersen:2008dd (). This hybrid approach is realized in the latest version 3.3 of UrQMD and can be downloaded from www.urqmd.org.

### 2.1 The initial state

The inital collisions are obtained from the Ultra-relativistic Quantum Molecular Dynamics model Bass:1998ca (); Bleicher:1999xi (). When the Lorentz-contracted nuclei have passed through each other all nucleons have potentially interacted at least once. This is the earliest time to assume local thermalization, which is needed for the validity of a hydrodynamic description of the expansion of the created matter. The transition time to the hydrodynamic evolution is

(1) |

where is the radius of the nuclei and is the Lorentz gamma factor of the two colliding nuclei in the center of mass frame. To map the particles from the UrQMD to the hydrodynamic fields they are described by Gaussian distributions with a width of fm. The net-baryon density, the energy density and the initial velocity profiles are then transformed onto a -dimensional space-grid to initialize the hydrodynamic evolution. Thereby the event-by-event fluctuations in the initial conditions are taken into account. The spectators are not considered in the hydrodynamic evolution but propagated further in the cascade.

### 2.2 The equation of state in the hydrodynamic evolution

The most important input to the hydrodynamic expansion of the system besides the initial conditions is the equation of state. It strongly influences the dynamics of the systems and naturally describes the phase transition in local equilibrium. In the framework of the hybrid model different equations of state can be implemented and their influence directly seen as the inital and final state remain unchanged. The hadron gas equation of state, for example, allows for a direct comparision of the underlying dynamics of the hydrodynamic treatment versus the transport sitmulation of the hot and dense region, since it includes the same degrees of freedom as in UrQMD and does not exhibit a phase transition. It can be used for baseline calculations to explore the influence of vanishing viscosity and local thermalisation. Other available equations of state are the Bag Model and the chiral equation of state. The Bag Model equation of state has a strong first order phase transition for all baryochemical potentials with a large latent heat. The chiral equation of state is obtained from a chiral hadronic Lagrangian that includes all baryons from the lowest flavor-SU() octet and the multiplets of scalar, pseudo-scalar, vector and axial-vector mesons. Recently a deconfinement equation of state has been derived by including partonic degrees of freedom in the chiral Lagrangian Dexheimer:2009hi (), very much in the same manner as the Polyakov-loop extention of quark meson models. It has been successfully implemented to include investigations on the confinement-deconfinement transition to the framework of the hybrid model Steinheimer:2009nn ().

The hydrodynamic expansion is performed by a full ()-dimensional SHASTA code Rischke:1995ir (); Rischke:1995mt () solving for energy and momentum conservation of an ideal fluid.

### 2.3 Freeze-out

When the system has expanded and diluted to a point where the hydrodynamical prescription becomes questionable and local equilibration can not be assumed any longer, a transition to a Boltzmann transport description is performed. In order to transfer the hydrodynamic densities to particles the Cooper-Frye freeze-out equation

(2) |

is employed. It connects the boosted phase-space distribution and the momentum-space distribution along the hypersurface . For the present calculations a gradual transition is chosen, which corresponds approximately to an iso-eigentime transition. Transverse slices of thickness fm are frozen out when the energy density in each cell of that slice has dropped below a critical energy density . It is important to note that for the correct transition the degrees of freedom of both sides, the hydrodynamic and the UrQMD Bravina:1998it (), need to be the same. Therefore, the last step in the hydrodynamic evolution is performed with an equation of state of a hadronic gas. The Monte Carlo sampling implementation of the Cooper-Frye equation conserves energy,baryon number, electric charge and strangeness in each single event.

## 3 Results from the hybrid model

We present results from the hybrid approach to heavy ion collisions, particle multiplicities and ratios of strange particles to pions.

### 3.1 Particle multiplicities and mean transverse mass excitation function

To study the dependence of the equation of state, viscous effects and nonequilibrium dynamics the hybrid model results can be compared to pure UrQMD 2.3 calculations. In Fig. (1) the multiplicities (top) and the mean transverse mass excitation function (bottom) are shown for pions (left) and kaons (right) compared to data from most central Au+Au/ Pb+Pb collisions at intermediate center-of-mass energies GeV. While the yield of pions is reduced in the hybrid model calculation compared to the transport calculation due to entropy conservation in the hydrodynamic expansion, as seen in the upper left plot of Fig. (1), it does only show a weak dependance on the equation of state. In the kaon multiplicities, upper right plot of Fig. (1), one clearly sees that the UrQMD 2.3 underpredicts the kaon yields in the whole energy range. In the Cooper-Frye freeze-out prescription the strange particles are produced according to their thermal distribution which leads to an enhanced yield in the hybrid model. Except for the Bag Model equation of state one sees a clear excess over the pure transport calculation, and a splitting between different equations of state. Since the transverse mass spectra are (Fig. (1, bottom) more sensitive to the pressure in the transverse plan it shows a more pronounced dependence on the equation of state. For both pions and kaons the chiral and the hadron gas equation of state give similar results. Only the Bag Model equation of state yields lower pion mean transverse masses at higher energies as it is expected for a first order phase transition. Again in the Bag Model equation of state the values of the mean transverse masses for kaons, lower right plot of Fig. (1), are even below the nonequilibrium dynamics of the transport model.

### 3.2 The kaon to pion ratio

The strangeness production is enhanced in the hybrid model due to the local thermal equilibrium asssumed for the hydrodynamic expansion Petersen:2009zi (). We can, therefore, expect to obtain different results for the kaon to pion and strange baryon to pion ratios. These are shown for central Au+Au/Pb+Pb collisions at center-of-mass energies GeV in Fig. (2). There are two different freeze-out criteria of (solid line) and (dashed line). The gray line is a pure UrQMD 2.3 calculation that underestimates the data in all four cases for higher energies. This is due to an over-production of pion yields compared to the strange particle production. As expected the hybrid model gives a better description of the data for energies above GeV. At lower energies the assumption of local thermal equilibrium might not be justified and nonequilibrium propagation seems to be necessary.

## 4 Critical fluctuations in chiral hydrodynamics

While the hybrid model gives an excellent framework to study the effects of an (local) equilibrium versus a full nonequlibrium transport propagation, it is not possible to propagate fluctuations directly in the hydrodynamic evolution. In models of chiral hydrodynamics Mishustin:1998yc (); Paech:2003fe () critical fluctuations at the phase transition can be studied directly and in non-equilibrium. For this purpose a field theoretical model that exhibits a phase transition is coupled to a background hydrodynamic evolution to simulate the expansion and cooling of the fireball in the course of a heavy ion collision. An effective model of QCD that exhibits a chiral phase transition and a critical point is the quark meson model with constituent quarks Scavenius:2000qd ()

(3) | |||

with the constituent quark field , the coupling between the quarks and the chiral fields . The potential is given by

(4) |

where the vacuum expectation of the sigma field is MeV and of the pion field is . with MeV is the explicit symmetry breaking term. Then . yields a sigma mass MeV. In order to have zero potential energy in the ground state a term is added.

### 4.1 The field equation of motion

The order parameter of chiral symmetry breaking is the sigma field. From thermodynamic considerations it is expected to fluctuate largely at the critical point. Our purpose is, however, to propagate the chiral fields in nonequilibrium and study the effects of the phase transition. It is expected that critical fluctuations at the critical point are diminished while nonequilibrium fluctuations at the first order phase transition are enhanced. We will in the following concentrate on the evolution of the sigma field and neglect fluctuations of the pion fields around their vacuum expectation value. The sigma field interacts with the quarks and antiquarks, which we assume to be in local thermal equilibrium. Due to this interaction the sigma field is damped and according to Langevin dynamics gets random kicks from the heat bath. The classical equation of motion for the chiral fields is

(5) |

where is the Fermi-Dirac distribution, the degeneracy factor and the energy of the quarks. The mass of the quarks is dynamically generated by a nonzero value of the sigma field. For the damping coefficient we take the value Biro:1997va (). Then the noise term is, in Markovian approximation,

(6) | |||||

(7) |

For the hydrodynamic evolution of the quark fluid the equation of state is needed. The pressure can be obtained from the thermodynamic relations in mean-field approximation

(8) |

The strength of the phase transition at can be tuned by changing the strength of the coupling . It is a crossover for small coupling and a first order phase transition for larger couplings, we use . Here, the effective potential has two degenerate minima allowing for the study of nucleation and spinodal decomposition Mishustin:1998eq (); Chomaz:2003dz (). For an intermediate coupling of a critical point with one very flat minimum can be found. Energy density and pressure are related by the thermodynamic relations

(9) | |||||

(10) |

Besides the limited degrees of freedom in this model, the main difference to the chiral equation of state used in the hybrid model calculations above is that the field value is not the minimum of the thermodynamic potential. It is explicitely propagated and, thus, the energy density and the pressure explicitely depend on the actual value of the sigma field. In the numerical simulations below the equation of motion of the sigma field (5) is solved by a staggered leap-frog algorithm. Initially we assume the sigma field to be in equilibrium with the fluid in order to be sensitive to the effects of the phase transition only.

### 4.2 The hdyrodynamic expansion

The equations of relativistic hydrodynamics of energy and momentum conservation are are

(11) |

where is the energy-momentum tensor of an ideal fluid. The source term reflects the energy dissipation from the field into the fluid. By using the explicit form of the equation of motion (5) we derive for the source term

(12) |

The inclusion of the source term guarantees energy and momentum conservation of the entire system. The fluid dynamic part is again propagated by a SHASTA code. The energy density is initiated as the equilibrium energy density at , ellipsoidal in x-y-plane and homogenous in z-direction, smoothed by a Woods Saxon type distribution.

## 5 Numerical results

The crucial quantity to look at is the intensity of the sigma fluctuations. It is given by

(13) |

This quantity can be considered as the number of sigma particles produced from the excited sigma field, once the nonlinearities can be neglected in the equation of motion of the sigma field. In order to take into account the mass change of the sigma particle in a hot environment the average temperature in the hot region is calculated. With the sigma mass

(14) |

the energy of the modes is

(15) |

The time evolution of (13) is shown in Fig. (3) for a first order phase transition and in Fig. (4) for a critical point. In a first order phase transition the intensity of sigma fluctuations between and fm is much larger than at a critical point. This is due to nonequilibrium effects at the first order phase transition where the effective potential has to two degenerate minima. Since we initialized the sigma field in equilibrium the increase in the sigma fluctuations is due to the phase transition. In the beginning of the evolution the sigma field is in the minimum which is stable above the critical temperatures . It becomes the unstable minimum for temperatures below the critical temperature . The barrier, which separates the two minima, however, prevents the sigma field from relaxing to the new equilibrium value. The damping to forces it further to remain in the local but unstable minimum. This effect is known as supercooling. The noise term turnes out to be too small to make nucleation an effective mechanism for relaxation. It is more likely that the field will remain in its unstable configuration almost until the barrier disappears and it can relax via spinodal decomposition Chomaz:2003dz ().

## 6 Summary

We have presented results from hybrid models to describe heavy ion collisions. The fully integrated Boltzmann and ()-hydrodynamic model is well established to test the effects of a nonequilibrium versus a (local) equilibrium evolution and of viscosity. By switching between different equations of state and comparing to the original UrQMD 2.3 transport model the effects of a the phase transition have been studied. It is interesting to mention that the assumption of local thermal equilibrium in the hybrid model gives a good description of the ’horn’ structure in ratios without invoking a phase transition. The implementation of a chiral fluid dynamic model, which combines the explicit propagation of the order parameter of the chiral phase transition and a fluid dynamic expansion of a quark fluid, enables us to study the phase transition in heavy ion collisions beyond the equation of state. We have included damping and noise term in the nonequilibrium dynamics of the chiral field in order to study the enhancement of the fluctuations at the phase transition and the relaxation to equilibrium. We have found that the intensity of fluctuations of the sigma field is significantly larger at the first order phase transition than at the critical point where large fluctuations in equilibrium would be expected.

## Acknowledgements

This work was supported by the Hessian LOEWE initiative through the Helmholtz International Center for FAIR (HIC for FAIR). The computational resources were provided by the Frankfurt Center for Scientific Computing (CSC). The authors thank Igor Mishustin and Hannah Peterson for fruitful and interesting discussions.

## References

- (1) S. Scherer et al., Prog. Part. Nucl. Phys. 42 (1999) 279.
- (2) J. Adams et al. [STAR Collaboration], Nucl. Phys. A 757, 102 (2005)
- (3) B. B. Back et al., Nucl. Phys. A 757, 28 (2005)
- (4) I. Arsene et al. [BRAHMS Collaboration], Nucl. Phys. A 757, 1 (2005)
- (5) K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 757, 184 (2005)
- (6) S. Borsanyi, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, C. Ratti and K. K. Szabo [Wuppertal-Budapest Collaboration], JHEP 1009 (2010) 073
- (7) W. Soldner [for the HotQCD collaboration], PoS LATTICE2010 (2010) 215
- (8) O. Scavenius, A. Mocsy, I. N. Mishustin and D. H. Rischke, Phys. Rev. C 64 (2001) 045202
- (9) C. Alt et al. [NA49 Collaboration], Phys. Rev. C 77 (2008) 024903
- (10) M. Gazdzicki, J. Phys. G 23 (1997) 1881 [arXiv:nucl-th/9706036].
- (11) M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. Lett. 81 (1998) 4816
- (12) M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. D 60 (1999) 114028
- (13) I. N. Mishustin, Phys. Rev. Lett. 82 (1999) 4779
- (14) H. Petersen, J. Steinheimer, G. Burau, M. Bleicher and H. Stocker, Phys. Rev. C 78 (2008) 044901
- (15) R. Andrade, F. Grassi, Y. Hama, T. Kodama and O. . J. Socolowski, Phys. Rev. Lett. 97 (2006) 202302 [arXiv:nucl-th/0608067].
- (16) D. Teaney, J. Lauret and E. V. Shuryak, RHIC,” arXiv:nucl-th/0110037.
- (17) T. Hirano, Prog. Theor. Phys. Suppl. 168 (2007) 347 [arXiv:0704.1699 [nucl-th]].
- (18) S. A. Bass et al., Prog. Part. Nucl. Phys. 41, 255 (1998)
- (19) M. Bleicher et al., J. Phys. G 25, 1859 (1999)
- (20) V. A. Dexheimer and S. Schramm, Phys. Rev. C 81 (2010) 045201
- (21) J. Steinheimer, V. Dexheimer, H. Petersen, M. Bleicher, S. Schramm and H. Stoecker, Phys. Rev. C 81 (2010) 044913
- (22) D. H. Rischke, S. Bernard and J. A. Maruhn, Nucl. Phys. A 595 (1995) 346
- (23) D. H. Rischke, Y. Pursun and J. A. Maruhn, Nucl. Phys. A 595 (1995) 383 [Erratum-ibid. A 596 (1996) 717]
- (24) L. V. Bravina et al., J. Phys. G 25 (1999) 351 [arXiv:nucl-th/9810036].
- (25) H. Petersen, M. Mitrovski, T. Schuster and M. Bleicher, Phys. Rev. C 80 (2009) 054910
- (26) S. V. Afanasiev et al. [The NA49 Collaboration], Phys. Rev. C 66, 054902 (2002)
- (27) L. Ahle et al. [E866 Collaboration and E917 Collaboration], Phys. Lett. B 476, 1 (2000)
- (28) J. L. Klay et al. [E-0895 Collaboration], Phys. Rev. C 68 (2003) 054905
- (29) M. K. Mitrovski et al. [NA49 Collaboration], J. Phys. G 32 (2006) S43
- (30) I. N. Mishustin and O. Scavenius, Phys. Rev. Lett. 83 (1999) 3134
- (31) K. Paech, H. Stoecker and A. Dumitru, Phys. Rev. C 68 (2003) 044907
- (32) T. S. Biro and C. Greiner, Phys. Rev. Lett. 79 (1997) 3138
- (33) P. Chomaz, M. Colonna and J. Randrup, Phys. Rept. 389 (2004) 263.