Advances in High Energy Physics

Volume 2014 (2014), Article ID 946934, 4 pages

http://dx.doi.org/10.1155/2014/946934

## Nuclear Stopping in Central Au+Au Collisions at RHIC Energies

College of Physics and Electronic Engineering, Guangxi Teachers Education University, Nanning, Guangxi 530001, China

Received 2 May 2014; Accepted 2 June 2014; Published 17 June 2014

Academic Editor: Fu-Hu Liu

Copyright © 2014 Ying Yuan and Nana Guan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.

#### Abstract

Nuclear stopping in central Au+Au collisions at relativistic heavy-ion collider (RHIC) energies is studied in the framework of a cascade mode and the modified ultrarelativistic quantum molecular dynamics (UrQMD) transport model. In the modified mode, the mean field potentials of both formed and “preformed” hadrons (from string fragmentation) are considered. It is found that the nuclear stopping is increasingly influenced by the mean-field potentials in the projectile and target regions with the increase of the reaction energy. In the central region, the calculations of the cascade model considering the modifying factor can describe the experimental data of the PHOBOS collaboration.

#### 1. Introduction

The study of strongly interacting matter at extreme temperatures and densities is provided a chance by heavy-ion collisions at ultrarelativistic energies [1–5]. In particular, the study of transport phenomena is majorly important to the understanding of many fundamental properties [6, 7]. The quark-gluon plasma (QGP) is expected to form at ultrarelativistic energies. Many microscopic transport models have been used to study the properties. The motivation is to provide a satisfied description of final-state hadron production for accurate and comprehensive simulations used in modern particle and nuclear physics experiments [8].

The main purpose of this work is to extract the information on nuclear stopping by comparison of the pseudorapidity distributions of charged particles from a transport-model simulation with data. This goal can be achieved by studying the pseudorapidity distribution of charged particles from central Au+Au collisions at relativistic heavy ion collider (RHIC) energies, within a transport model, the ultrarelativistic quantum molecular dynamics (UrQMD) model. And the modifying factor is considered in this cascade mode. This method is advantageous to directly compare existing data, and it can describe the experimental data very well.

#### 2. Ultrarelativistic Quantum Molecular Dynamics Transport Model

The UrQMD model is a microscopic many-body transport approach and can be applied to study proton-proton (pp), proton-nucleus (pA), and nucleus-nucleus (AA) interactions over an energy range from the heavy-ion synchrotron (SIS) to RHIC. This transport model is based on the covariant propagation of color strings, constituent quarks, and diquarks (as string ends) accompanied by mesonic and baryonic degree of freedom [9]. In present model, the subhadronic degrees of freedom enter via the introduction of a formation time for hadrons produced in the fragmentation of strings [10–12], which are dominant at the early stage of heavy ion collisions (HICs) with high super proton synchrotron (SPS) and RHIC energies.

##### 2.1. The Mean-Field Treatments

The UrQMD model is based on parallel principles as the quantum molecular dynamics (QMD) model; hadrons are represented by Gaussian wave packets in phase space and the phase space of hadron is propagated according to Hamilton’s equation of motion [13],

Here, and are the coordinate and momentum of the hadron , respectively. The points on the coordinate and momentum represent the coordinate and momentum for the time derivative. The Hamiltonian consists of the kinetic energy and the effective interaction potential energy , In the standard UrQMD model, the potential energy includes the two-body and three-body Skyrme-, Yukawa-, Coulomb-, and Pauli-terms [13, 14],

In the modified version of UrQMD (based on the version 3.3), the following two terms are further added: (1) the density-dependent symmetry potential term and (2) the momentum-dependent term [15]. In this work, the soft momentum dependent (SM) equation of state (EoS) is used, which is described in [16]. The Yukawa-, Pauli-, and symmetry-potentials of baryons become negligible, while the Skyrme- and the momentum-dependent parts of potentials still influence the whole dynamical process of HICs above the 2 GeV [17]. During the formation time, the “preformed” particles (string fragments that will be projected onto hadron states later on) are usually treated to be free-streaming, while reduced cross sections are only included for leading hadrons. Recently, the mean-field potentials for both formed and “preformed” particles are considered for a better understanding of Hanbury-Brown-Twiss (HBT) time-related puzzle [18] and nuclear stopping [19]. In this paper, the pseudorapidity distribution about charged particles at RHIC energies will be shown. And the relativistic effect on the relative distance and the relative momentum and a covariance-related reduced factor used for the update of potentials [17, 20] are considered in calculations.

##### 2.2. The Modifying Factor for the UrQMD Calculations

It is known that the calculation results of the UrQMD model about the (pseudo)rapidity distribution of charged particles are somewhat larger than that of the experimental data [9, 19, 21, 22]. So, in this work, the modifying factor for the UrQMD calculations is led in, and it is marked as which multiplies directly the pseudorapidity distribution of charged particles as the renormalization. At RHIC energies, this factor increases with the increase of the reaction energy.

#### 3. Nuclear Stopping at RHIC Energies

At RHIC energies, as the pseudorapidity distribution of charged particles in the transverse direction has not been provided by experimental physicists, we study the nuclear stopping with the longitudinal pseudorapidity distribution. Figures 1, 2, and 3 depict the pseudorapidity distributions of charged particles for central Au+Au collisions (<6% of total cross section ) at RHIC, respectively. The experimental data (circles) of charged particles produced in Au+Au collisions at , 130, and 200 GeV are taken from [23, 24]. The curves are our calculated results with different treatments shown in the panels. In the calculations, besides a cascade mode shown in the panel, we also show the results with potentials of both formed and “preformed” hadrons (“pf-part and f-B SM-EoS”) in the panel. Cross sections used in the model are not modified by the nuclear medium in this energy region.

In Figure 1 ( GeV), it is found that the shape of the pseudorapidity distributions of measured charged particles is similar to the cascade calculations at GeV, but the results have only a half of the experimental data. The calculations with potentials give a Gaussian-like distribution, and they cannot describe the pseudorapidity distribution of charged particles completely. When the modifying factor is imported in the cascade calculations, it can describe well the experimental data. At this center-of-mass energy, the value of is 1.72 with per degree of freedom () being 1.04.

In Figure 2 ( GeV), the curves and scatters are similar to Figure 1. But the difference is that the cascade calculations cannot describe the experimental data, even if the modifying factor is imported in the cascade calculations. When the value of is 1.9 with being 0.33, the cascade results are according to the experimental data at . We define this region as central region. The left of central region is called target region, and the right is called projectile region. At this energy, the reaction principle of the projectile (or target) region and central region is very different. So in this work, we consider that the calculations with potentials in the regions of and are moved horizontally to both sides, respectively. In other words, the calculations of pseudorapidity distribution in the target region are moved to left side, and the corresponding values of reduce one. In a similar way, the calculations in the projectile region are moved to right side, and the corresponding values of are plus one. In this figure, it can be seen that the calculations of the UrQMD model in view of this method are in keeping with the experimental data as well.

In Figure 3 ( GeV), the curves and scatters are similar to Figure 2, and the processed method is also the same as Figure 2. But the value of is 1.95 with being 0.52, and the cascade results are according to the experimental data at . The processing method is similar to Figure 2. In this figure, the calculations of pseudorapidity distribution in the region of are moved to left side, and the corresponding values of reduce one. In a similar way, the calculations in the region of are moved to right side, and the corresponding values of are plus one. In this figure, it can be seen that the calculations of the UrQMD model in view of this method are in keeping with the experimental data as well.

From these three figures, it is known that the value of modifying factor increases with the increase of reaction energy. At RHIC energies, with the growth of reaction energy, the regions are enlarging which is affected by the potentials of both formed and “preformed” hadrons. At GeV, since the reaction energy is not large enough, the impact of the mean-field potentials is not necessary to be considered in this reaction energy. In other words, there is no obvious difference between the projectile (target) and central regions. When the reaction energy is taken GeV, a lot of energy is gathered, and the QGP has been formed in the central region. In this region, the calculations of cascade UrQMD model can describe the experimental data very well. A fraction of energy is left in the projectile or target region, so we need to consider the impact of the mean-field potentials, like in Figure 2. At GeV, it is similar to the case in Figure 2. The only different thing is that the central region is enlarged comparing with the Figure 2. When the QGP has been formed, the pseudorapidity distribution of charged particles will not be affected by the mean-field potentials.

#### 4. Summary and Outlook

In summary, we have presented the pseudorapidity distribution of charged particles for central Au+Au reactions at RHIC energies. The UrQMD transport model (version 3.3) with the cascade mode and the modified mode (considering the potentials of both formed and “preformed” hadrons) has been used in all calculations. Based on this model we have investigated the effect of mean-field potentials on nuclear stopping. It is found that nuclear stopping is increasingly influenced by the mean-field potentials in the projectile and target regions with the increase of reaction energy. Since most of the incident energy is concentrated in the central region, the calculations of the cascade model considering the modifying factor can describe the experimental data in the projectile and target region. In conclusion, our combined calculations with the cascade mode and the modified mode in the different regions can give the agreement between calculation results and the experimental data. But the intrinsic reaction mechanism has yet to be studied in depth.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors acknowledge support by the Cluster-Computing Center of School of Science (C3S2) of Huzhou Teachers College for Scientific Computing. This work is supported by the Natural Science Foundation of Guangxi Zhuangzu Autonomous Regions of China under Grant no. 2012GXNSFBA053011 and the Starting Foundation of Scientific Research of the Doctor of Guangxi Teachers Education University of China.

#### References

- C. Alt, T. Anticic, B. Baatar et al., “Energy dependence of Λ and
*Ξ*production in central Pb + Pb collisions at 20A,30A,40A,80A, and 158A GeV measured at the CERN Super Proton Synchrotron,”*Physical Review C: Nuclear Physics*, vol. 78, no. 3, Article ID 034918, 2008. View at Publisher · View at Google Scholar · View at Scopus - J.-X. Sun, F.-H. Liu, and E.-Q. Wang, “Pseudorapidity distributions of charged particles and contributions of leading nucleons in Cu-Cu collisions at high energies,”
*Chinese Physics Letters*, vol. 27, no. 3, Article ID 032503, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. A. Rahim, S. Fakhraddin, E.-Q. Wang, F.-H. Liu, and J.-X. Sun, “Singly and doubly charged projectile fragments in nucleus-emulsion collisions at dubna energy in the framework of the multi-source model,”
*Chinese Physics Letters*, vol. 28, no. 8, Article ID 082501, 2011. View at Publisher · View at Google Scholar · View at Scopus - B.-C. Li and M. Huang, “Strongly coupled matter near phase transition,”
*Journal of Physics G: Nuclear and Particle Physics*, vol. 36, no. 6, Article ID 064062, 2009. View at Publisher · View at Google Scholar · View at Scopus - F.-H. Liu, “Anisotropic emission of charged mesons and structure characteristic of emission source in heavy ion collisions at 1-2A GeV,”
*Chinese Physics B*, vol. 17, no. 3, pp. 883–895, 2008. View at Publisher · View at Google Scholar · View at Scopus - B.-C. Li, Y.-Y. Fu, L.-L. Wang, and F.-H. Liu, “Dependence of elliptic flows on transverse momentum and number of participants in Au + Au collisions at 200 GeV,”
*Journal of Physics G: Nuclear and Particle Physics*, vol. 40, no. 2, Article ID 025104, 2013. View at Publisher · View at Google Scholar · View at Scopus - C. Escano-Rodriguez, D. Durand, A. Chbihi, and J. D. Frankland, “A study of nuclear stopping in central symmetric nuclear collisions at intermediate energies,” http://arxiv.org/abs/nucl-ex/0503007.
- K. Abdel-Waged, N. Felemban, and V. V. Uzhinskii, “Geant4 hadronic cascade models analysis of proton and charged pion transverse momentum spectra from p + Cu and Pb collisions at 3, 8, and 15 GeV/c,”
*Physical Review C: Nuclear Physics*, vol. 84, no. 1, Article ID 014905, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. Petersen, M. Bleicher, S. A. Bass, and H. Stocker, “UrQMD-2.3—changes and comparisons,” http://arxiv.org/abs/0805.0567.
- B. Andersson, G. Gustafson, and B. Nilsson-Almqvist, “A model for low-pT hadronic reactions with generalizations to hadron-nucleus and nucleus-nucleus collisions,”
*Nuclear Physics B*, vol. 281, no. 1-2, pp. 289–309, 1987. View at Google Scholar · View at Scopus - B. Nilsson-Almqvist and E. Stenlund, “Interactions between hadrons and nuclei: the Lund Monte Carlo, FRITIOF version 1.6,”
*Computer Physics Communications*, vol. 43, no. 3, pp. 387–397, 1987. View at Google Scholar · View at Scopus - T. Sjöstrand, “High-energy-physics event generation with PYTHIA 5.7 and JETSET 7.4,”
*Computer Physics Communications*, vol. 82, no. 1, pp. 74–89, 1994. View at Google Scholar · View at Scopus - S. A. Bass, M. Belkacem, M. Bleicher et al., “Microscopic models for ultrarelativistic heavy ion collisions,”
*Progress in Particle and Nuclear Physics*, vol. 41, pp. 255–369, 1998. View at Google Scholar · View at Scopus - M. Bleicher, E. Zabrodin, C. Spieles et al., “Relativistic hadron-hadron collisions in the ultra-relativistic quantum molecular dynamics model,”
*Journal of Physics G: Nuclear and Particle Physics*, vol. 25, no. 9, pp. 1859–1896, 1999. View at Publisher · View at Google Scholar · View at Scopus - S. A. Bass, C. Hartnack, H. Stöcker, and W. Greiner, “Azimuthal correlations of pions in relativistic heavy-ion collisions at 1 GeV/nucleon,”
*Physical Review C*, vol. 51, no. 6, pp. 3343–3356, 1995. View at Publisher · View at Google Scholar · View at Scopus - Q. Li, Z. Li, S. Soff, M. Bleicher, and H. Stöcker, “Probing the equation of state with pions,”
*Journal of Physics G: Nuclear and Particle Physics*, vol. 32, no. 2, pp. 151–164, 2006. View at Publisher · View at Google Scholar · View at Scopus - Q. Li and M. Bleicher, “A model comparison of resonance lifetime modifications, a soft equation of state and non-Gaussian effects on
*π*-*π*Correlations at FAIR/AGS energies,”*Journal of Physics G: Nuclear and Particle Physics*, vol. 36, no. 1, Article ID 015111, 2009. View at Publisher · View at Google Scholar · View at Scopus - Q. Li, M. Bleicher, and H. Stöcker, “The effect of pre-formed hadron potentials on the dynamics of heavy ion collisions and the HBT puzzle,”
*Physics Letters B: Nuclear, Elementary Particle and High-Energy Physics*, vol. 659, no. 3, pp. 525–530, 2008. View at Publisher · View at Google Scholar · View at Scopus - Y. Yuan, Q. F. Li, Z. X. Li, and F. H. Liu, “Transport model study of nuclear stopping in heavy-ion collisions over the energy range from 0.09A to 160A GeV,”
*Physical Review C*, vol. 81, no. 3, Article ID 034913, 2010. View at Google Scholar - M. Isse, A. Ohnishi, N. Otuka, P. K. Sahu, and Y. Nara, “Mean-field effects on collective flows in high-energy heavy-ion collisions from AGS to SPS energies,”
*Physical Review C: Nuclear Physics*, vol. 72, no. 6, Article ID 064908, 2005. View at Publisher · View at Google Scholar · View at Scopus - W. Reisdorf, “Dynamics of multifragmentation in heavy ion collisions,”
*Progress of Theoretical Physics*, no. 140, pp. 111–133, 2000. View at Google Scholar · View at Scopus - K. Dey and B. Bhattacharjee, “Rapidity dependence of the produced particles at FAIR energies,”
*Physical Review C*, vol. 89, Article ID 054910, 2014. View at Google Scholar - B. B. Back and PHOBOS Collaboration, “Charged-particle pseudorapidity distributions in Au + Au collisions at $\sqrt{{S}_{NN}=62.4}$ GeV,”
*Physical Review C*, vol. 74, no. 2, Article ID 021901, 2006. View at Google Scholar - B. B. Back and PHOBOS Collaboration, “Significance of the fragmentation region in ultrarelativistic heavy-ion collisions,”
*Physical Review Letters*, vol. 91, no. 5, Article ID 052303, 2003. View at Google Scholar