arxiv:nucl-th/ v3 20 Dec 2005

Transcription

1 A Multi-Phase Transort Model for Relativistic Heav Ion Collisions Zi-Wei Lin Phsics Deartment, The Ohio State Universit, Columbus, Ohio 4321 Che Ming Ko Cclotron Institute and Phsics Deartment, Texas A&M Universit, College Station, Texas arxiv:nucl-th/41111v3 2 Dec 25 Bao-An Li and Bin Zhang Deartment of Chemistr and Phsics, Arkansas State Universit, State Universit, Arkansas Subrata Pal Deartment of Phsics, Michigan State Universit, East Lansing, Michigan (Dated: Februar 8, 28) We describe in detail how the different comonents of a multi-hase transort (AMPT) model, that uses the Heav Ion Jet Interaction Generator (HIJING) for generating the initial conditions, Zhang s Parton Cascade (ZPC) for modeling artonic scatterings, the Lund string fragmentation model or a quark coalescence model for hadronization, and A Relativistic Transort (ART) model for treating hadronic scatterings, are imroved and combined to give a coherent descrition of the dnamics of relativistic heav ion collisions. We also exlain the wa arameters in the model are determined, and discuss the sensitivit of redicted results to hsical inut in the model. Comarisons of these results to exerimental data, mainl from heav ion collisions at the Relativistic Heav Ion Collider (RHIC), are then made in order to extract information on the roerties of the hot dense matter formed in these collisions. PACS numbers: q, Mh, 24.1.Lx I. INTRODUCTION Colliding heav ions at relativistic energies makes it ossible to subject nuclear matter to the extreme condition of large comression, leading to energ densities that can exceed that for roducing a lasma of deconfined quarks and gluons, that is believed to have existed during the first microsecond after the Big Bang. Exeriments at RHIC at the Brookhaven National Laborator with center-of-mass energ u to s NN = 2 GeV in Au+Au collisions thus rovide the oortunit to stud the roerties of this so-called quark-gluon lasma (QGP). At the future Large Hadron Collider (LHC) at CERN, which will allow Pb+Pb collisions at s NN = 5.5 TeV, the roduced quark-gluon lasma will have an even higher temerature and a nearl vanishing net baron chemical otential. Man observables have been measured at RHIC, such as the raidit distributions of various articles and their transverse momentum sectra u to ver high transverse momentum, the centralit deendence of these observables, the ellitic flows of various articles, as well as both identical and non-identical two-article correlations. To understand these extensive exerimental results, man theoretical models have been introduced. It ranges from thermal models [1, 2, 3, 4] based on the assumtion Present address: 31 Sarkman Drive, VBRH E-39, The Universit of Alabama in Huntsville, Huntsville, AL of global thermal and chemical equilibrium to hdrodnamic models [5, 6, 7, 8, 9, 1, 11] based onl on the assumtion of local thermal equilibrium, and to transort models [12, 13, 14, 15, 16, 17, 18, 19, 2, 21, 22, 23, 24, 25, 26] that treat non-equilibrium dnamics exlicitl. The thermal models have been ver successful in accounting for the ield of various articles and their ratios, while the hdrodnamic models are articularl useful for understanding the collective behavior of low transverse momentum articles such as the ellitic flow [8, 9, 1, 11]. Since transort models treat chemical and thermal freeze-out dnamicall, the are also natural and owerful tools for studing the Hanbur-Brown-Twiss interferometr of hadrons. For hard rocesses that involve large momentum transfer, aroaches based on the erturbative quantum chromodnamics (QCD) using arton distribution functions in the colliding nuclei have been used [27, 28]. Also, the classical Yang-Mills theor has been develoed to address the evolution of arton distribution functions in nuclei at ultra-relativistic energies [29, 3, 31] and used to stud the hadron raidit distribution and its centralit deendence at RHIC [32, 33, 34]. These roblems have also been studied in the QCD based final-state saturation model [35, 36, 37]. Although studies based on the QCD [38] have shown that thermalization could be achieved in collisions of ver large nuclei and/or at extremel high energ, even though the strong couling constant at the saturation scale is asmtoticall small, the dense matter created in heav ion collisions at RHIC ma, however, not achieve full thermal or chemical equilibrium as a re-

2 2 sult of its finite volume and energ. To address such non-equilibrium man-bod dnamics, we have develoed a multi-hase transort (AMPT) model that includes both initial artonic and final hadronic interactions and the transition between these two hases of matter [39, 4, 41, 42, 43, 44, 45, 46, 47, 48, 49, 5]. The AMPT model is constructed to describe nuclear collisions ranging from + A to A + A sstems at center-of-mass energies from about s NN = 5 GeV u to 55 GeV at LHC, where strings and minijets dominate the initial energ roduction and effects from final-state interactions are imortant. For the initial conditions, the AMPT model uses the hard minijet artons and soft strings from the HIJING model. The ZPC arton cascade is then used to describe scatterings among artons, which is followed b a hadronization rocess that is based on the Lund string fragmentation model or b a quark coalescence model. The latter is introduced for an extended AMPT model with string melting in which hadrons, that would have been roduced from string fragmentation, are converted instead to their valence quarks and antiquarks. Scatterings among the resulting hadrons are described b the ART hadronic transort model. With arameters, such as those in the string fragmentation, fixed b the exerimental data from heav ion collisions at SPS, the AMPT model has been able to describe reasonabl man of the exerimental observations at RHIC. In this aer, we give a detailed descrition of the different comonents of the AMPT model, discuss the arameters in the model, show the sensitivit of its results to the inut to the model, and comare its redictions with exerimental data. The aer is organized as follows. In Section II, we describe the different comonents of the AMPT model: the HIJING model and string melting, the ZPC model, the Lund string fragmentation model, and the quark coalescence model used for the scenario of string melting, and the extended ART model. Tests of the AMPT model against data from and reactions are given in Section III. Results from the AMPT model for heav ion collisions at SPS energies are discussed in Section IV for hadron raidit distributions and transverse momentum sectra, baron stoing, and antiroton roduction. In Section V, we show results at RHIC for hadron raidit distributions and transverse momentum sectra, article ratios, baron and antibaron roduction, and the roduction of multistrange barons as well as J/ψ. We further show results from the AMPT model with string melting on hadron ellitic flows and two-ion interferometr at RHIC. In Section VI, we resent the redictions from AMPT for hadron raidit and transverse momentum distributions in Pb+Pb collisions at the LHC energ. Discussions on ossible future imrovements of the AMPT model are resented in Section VII, and a summar is finall given in Section VIII. II. THE AMPT MODEL The AMPT model consists of four main comonents: the initial conditions, artonic interactions, the conversion from the artonic to the hadronic matter, and hadronic interactions. The initial conditions, which include the satial and momentum distributions of minijet artons and soft string excitations, are obtained from the HIJING model [51, 52, 53, 54]. Currentl, the AMPT model uses the HIJING model version [55], which does not include baron junctions [56]. Scatterings among artons are modeled b Zhang s arton cascade (ZPC) [18], which at resent includes onl two-bod scatterings with cross sections obtained from the QCD with screening masses. In the default AMPT model [39, 4, 41, 42, 43, 44, 46, 47, 49], artons are recombined with their arent strings when the sto interacting, and the resulting strings are converted to hadrons using the Lund string fragmentation model [57, 58, 59]. In the AMPT model with string melting [45, 48, 5], a quark coalescence model is used instead to combine artons into hadrons. The dnamics of the subsequent hadronic matter is described b a hadronic cascade, which is based on the ART model [14, 25] and extended to include additional reaction channels that are imortant at high energies. These channels include the formation and deca of K resonance and antibaron resonances, and baronantibaron roduction from mesons and their inverse reactions of annihilation. Final results from the AMPT model are obtained after hadronic interactions are terminated at a cutoff time (t cut ) when observables under stud are considered to be stable, i.e., when further hadronic interactions after t cut will not significantl affect these observables. We note that two-bod artonic scatterings at all ossible times have been included because the algorithm of ZPC, which roagates artons directl to the time when the next collision occurs, is fundamentall different from the fixed time ste method used in the ART model. In Figs. 1 and 2, we show, resectivel, the schematic structures of the default AMPT model [39, 4, 41, 42, 43, 44, 46, 47] and the AMPT model with string melting [45, 48, 5] described above. The full source code of the AMPT model in the Fortran 77 language and instructions for users are available online at the OSCAR website [6] and also at the EPAPS website [61]. The default AMPT model is named as version 1.x, and the AMPT model with string melting is named as version 2., where value of the integer extension x or increases whenever the source code is modified. Current versions of the AMPT model is 1.11 for the default model and 2.11 for the string melting model, resectivel. In the following, we exlain in detail each of the above four comonents of the AMPT model and the wa the are combined to describe relativistic heav ion collisions.

3 3 FIG. 1: (Color online) Illustration of the structure of the default AMPT model. FIG. 2: (Color online) Illustration of the structure of the AMPT model with string melting. A. Initial conditions 1. The default AMPT model In the default AMPT model, initial conditions for heav ion collisions at RHIC are obtained from the HI- JING model [51, 52, 53, 54]. In this model, the radial densit rofiles of the two colliding nuclei are taken to have Woods-Saxon shaes, and multile scatterings among incoming nucleons are treated in the eikonal formalism. Particle roduction from two colliding nucleons is described in terms of a hard and a soft comonent. The hard comonent involves rocesses in which the momentum transfer is larger than a cutoff momentum and is evaluated b the QCD using the arton distribution function in a nucleus. These hard rocesses lead to the roduction of energetic minijet artons and are treated via the PYTHIA rogram. The soft comonent, on the other hand, takes into account non-erturbative rocesses with momentum transfer below and is modeled b the formation of strings. The excited strings are assumed to deca indeendentl according to the Lund JETSET fragmentation model. From the and total cross sections and the ratio of σ el /σ tot in the energ range 2 < s < 18 GeV, it has been found that the exerimental data can be fitted with a nucleon-nucleon soft cross section σ s (s) = 57 mb at high energies and = 2 GeV/c [51]. The indeendence of these two arameters on the colliding energ is due to the use of the Duke-Owens set 1 for the arton distribution function [62] in the nucleon. With different arton distribution functions, an energ-deendent ma be needed to fit the same and data [63, 64]. We note that since the number of hard collisions in an A+A collision roughl scales as A 4/3 and grows fast with colliding energ while the number of strings roughl scales as A, minijet roduction becomes more imortant when the energ of heav ion collisions increases [51, 65]. Because of nuclear shadowing, both quark [66] and gluon [67] distribution functions in nuclei are different from the simle suerosition of their distributions in a nucleon. This effect has been included in the HIJING model via the following imact-arameter-deendent but Q 2 (and flavor)-indeendent arameterization [52]: R A (x, r) fa a (x, Q 2, r) Af N a (x, Q 2 ) = ln 1/6 A(x 3 1.2x x) [ α A (r) 1.8(A1/3 1) ] x e x2 /.1,(1) ln(a + 1) where x is the light-cone momentum fraction of arton a, and f a is the arton distribution function. The imactarameter deendence of the nuclear shadowing effect is controlled b α A (r) =.133(A 1/3 1) 1 r 2 /RA 2, (2) with r denoting the transverse distance of an interacting nucleon from the center of the nucleus with radius R A = 1.2A 1/3. Note that there is a modified HIJING model which uses a different arameterization for the nuclear shadowing that is also flavor-deendent [63]. To take into account the Lorentz boost effect, we have introduced a formation time for minijet artons that de-

4 4 ends on their four momenta [68]. Secificall, the formation time for each arton in the default AMPT model is taken to have a Lorentzian distribution with a half width t f = E/m 2 T, where E and m T are the arton energ and transverse mass, resectivel. Initial ositions of formed minijet artons are calculated from those of their arent nucleons using straight-line trajectories. 2. The AMPT model with string melting Although the artonic art in the default AMPT model includes onl minijets from the HIJING model, its energ densit can be ver high in heav ion collisions at RHIC. As shown in Fig. 3 for the time evolutions of the energ and number densities of artons and hadrons in the central cell of central (b = fm) Au+Au collisions at s NN = 2 GeV in the center-of-mass frame, the artonic energ densit during the first few fm/c s of the collision is more than an order of magnitude higher than the critical energ densit ( 1 GeV/fm 3 ) for the QCD hase transition, similar to that redicted b the high densit QCD aroach [69]. The shar increase in energ and number densities at about 3 fm/c is due to the exclusion of energies that are associated with the excited strings in the artonic stage. Keeing strings in the high energ densit region [7] thus underestimates the artonic effect in these collisions. We note that the central cell in the above calculation is chosen to have a transverse radius of 1 fm and a longitudinal dimension between.5t and.5t, where time t starts when the two nuclei are full overlaed in the longitudinal direction. To model the above effect in high energ densit regions, we extend the AMPT model to include the string melting mechanism [45, 48, 5], i.e, all excited strings, that are not rojectile and target nucleons without an interactions, are converted to artons according to the flavor and sin structures of their valence quarks. In articular, a meson is converted to a quark and an antiquark, while a baron is first converted to a quark and a diquark with weights according to relations from the SU(6) quark model [71], and the diquark is then decomosed into two quarks. The quark and diquark masses are taken to be the same as in the PYTHIA rogram [59], e.g. m u = 5.6 MeV/c 2, m d = 9.9 MeV/c 2, and m s = 199 MeV/c 2. We further assume that the above two-bod decomosition is isotroic in the rest frame of the arent hadron or diquark, and the resulting artons do not undergo scatterings until after a formation time given b t f = E H /m 2 T,H, with E H and m T,H denoting the energ and transverse mass of the arent hadron. Similar to minijet artons in the default AMPT model, initial ositions of the artons from melted strings are calculated from those of their arent hadrons using straight-line trajectories. The above formation time for artons is introduced to reresent the time needed for their roduction from strong color fields. Although we consider hadrons be- n (fm 3 ), ε (GeV/fm 3 ) Au+Au (s 1/2 =2A GeV,b=) arton stage hadron stage t (fm/c) FIG. 3: Energ and number densities of minijet artons and formed hadrons in the central cell as functions of time for central (b= fm) Au+Au collisions at s NN = 2 GeV from the default AMPT model, where the energ stored in the excited strings is absent in the arton stage and is released onl when hadrons are formed. fore string melting as a convenient ste in modeling the string melting rocess, choosing a formation time that deends on the momentum of the arent hadron ensures that artons from the melting of same hadron would have the same formation time. The advantage of this choice is that the AMPT model with string melting reduces to HIJING results in the absence of artonic and hadronic interactions as these artons would then find each other as closest artners at the same freeze-out time and thus coalesce back to the original hadron. We note that the tical string fragmentation time of about 1 fm/c is not alied to the melting of strings as the fragmentation rocess involved here is considered just as an intermediate ste in modeling arton roduction from the energ field of the strings in an environment of high energ densit. B. Parton cascade In the transort aroach, interactions among artons are described b equations of motion for their Wigner distribution functions, which describe semi-classicall their densit distributions in hase sace. These equations can be aroximatel written as the following Boltzmann equations: µ µ f a (x,, t) n ε

5 5 = m n b 1,b 2,,b m c 1,c 2,,c n m i=1 n j=1 ( m (2π) 4 δ 4 bk [ k=1 d 3 bi (2π) 3 2E bi f bi (x, bi, t) d 3 cj (2π) 3 2E cj M m n 2 ) n cl l=1 m δ abq δ 3 ( bq ) + q=1 ] n δ acr δ 3 ( cr ).(3) r=1 In the above, f a (x,, t) is the distribution function of arton te a at time t in the hase sace, and M m n denotes the matrix element of the multi-arton interaction m n. If one considers onl two-bod interactions, these equations reduce to µ µ f(x,, t) σf(x 1, 1, t)f(x 2, 2, t), (4) where σ is the cross section for artonic two-bod scattering, and the integral is evaluated over the momenta of other three artons with the integrand containing factors such as a delta function for momentum conservation. The Boltzmann equations are solved using Zhang s arton cascade (ZPC) [18], in which two artons undergo scattering whenever the aroach each other with a closest distance smaller than σ/π. At resent, ZPC includes onl arton two-bod scattering such as gg gg with cross sections calculated from the QCD. For gluon elastic scattering, the leading-order QCD gives dσ gg dt = 9πα2 s 2s 2 9πα2 s 2 ( 3 ut ( 1 t u 2 s 2 us t 2 st ) u 2 ), (5) where α s is the strong couling constant, and s, t and u are standard Mandelstam variables for elastic scattering of two artons. The second line in the above equation is obtained b keeing onl the leading divergent terms. Since the scattering angle ranges from to π/2 for identical articles, one then has [18] dσ gg dt 9πα2 s 2t 2, (6) if the scattering angle is between and π. The singularit in the total cross section can be regulated b a Debe screening mass µ, leading to dσ gg dt 9πα 2 s 2(t µ 2 ) 2, σ gg = 9πα2 s 2µ µ 2 /s. (7) The screening mass µ is generated b medium effects and is thus related to the arton hase-sace densit. For the artonic sstem exected to be formed in Au+Au collisions at RHIC, the value of µ is on the order of one inverse fermi [18]. For massless artons in a lasma at temerature T, their average colliding energ is s 18T, thus µ < s for µ = 3 fm 1 leads to the requirement T > 141 MeV. Since s > µ 2 generall holds in hot QGP, the following simlified relation between the total arton elastic scattering cross section and the screening mass is used in the ZPC [72] σ gg 9πα2 s 2µ 2. (8) A value of 3 fm 1 for the screening mass µ thus leads to a total cross section of about 3 mb for the elastic scattering between two gluons. B changing the value of the screening mass µ, different cross sections can be obtained, and this will be used in studing the effect of arton cross sections in heav ion collisions at RHIC. This cross section is used in AMPT not onl in the default model, which includes onl scatterings of minijet gluons, but also in the string melting model, which onl includes scatterings of quarks/antiquarks of all flavors. We have therefore neglected in the latter case the difference between the Casimir factors for quarks and gluons. We note that minijet artons roduced from hard scatterings in the HIJING model can lose energ b gluon slitting and transfer their energies to nearb soft strings. In the AMPT model, this so-called jet quenching in the HIJING model is relaced b arton scatterings in ZPC. Since onl two-bod scatterings are included in ZPC, higher-order contributions to the jet energ loss are still missing in the AMPT model. C. Hadronization Two different hadronization mechanisms are used in the AMPT model for the two different initial conditions introduced in Sec. II A. In the default AMPT model, minijets coexist with the remaining art of their arent nucleons, and the together form new excited strings after artonic interactions. Hadronization of these strings are described b the Lund string model. In the AMPT model with string melting, these strings are converted to soft artons, and their hadronization is based on a simle quark coalescence model, similar to that in the ALCOR model [73]. 1. Lund string fragmentation for the default AMPT model Hadron roduction from the minijet artons and soft strings in the default AMPT model is modeled as follows. After minijet artons sto interacting, i.e., after the no longer scatter with other artons, the are combined with their arent strings to form excited strings, which are then converted to hadrons according to the Lund string

6 6 fragmentation model [57, 58]. In the Lund model as imlemented in the JETSET/PYTHIA routine [59], one assumes that a string fragments into quark-antiquark airs with a Gaussian distribution in transverse momentum. A suression factor of.3 is further introduced for the roduction of strange quark-antiquark airs relative to that of light quark-antiquark airs. Hadrons are formed from these quarks and antiquarks b using a smmetric fragmentation function [57, 58]. Secificall, the transverse momentum of a hadron is given b those of its constituent quarks, while its longitudinal momentum is determined b the Lund smmetric fragmentation function [74] f(z) z 1 (1 z) a ex( b m 2 /z), (9) with z denoting the light-cone momentum fraction of the roduced hadron with resect to that of the fragmenting string. The average squared transverse momentum is then given b 2 2 = f(z)d 2 dz f(z)d2 dz zmax z(1 z) a ex( b m 2 /z)dz = b z max (1 z) a ex( b m 2 /z)dz. (1) For massless articles, it reduces to 2 = 1 1 z(1 z)a dz b 1 (1 z)a dz = 1 b(2 + a). (11) Since quark-antiquark air roduction from string fragmentation in the Lund model is based on the Schwinger mechanism [75] for article roduction in strong field, its roduction robabilit is roortional to ex( πm 2 /κ), where κ is the string tension, i.e., the energ in a unit length of string. Due to its large mass, strange quark roduction is suressed b the factor e π(m2 s m2 u )/κ, comared to that of light quarks. Also, the average squared transverse momentum of roduced articles is roortional to the string tension, i.e., 2 κ. Comaring this with Eq. (11), one finds that the two arameters a and b in the Lund fragmentation function are aroximatel related to the string tension b κ 1 b(2 + a). (12) After roduction from string fragmentation, hadrons are given an additional roer formation time of.7 fm/c [76]. Positions of formed hadrons are then calculated from those of their arent strings b following straightline trajectories. 2. Quark coalescence for the AMPT model with string melting After artons in the string melting scenario sto interacting, we model their hadronization via a simle quark coalescence model b combining two nearest artons into a meson and three nearest quarks (antiquarks) into a baron (antibaron). Since the invariant mass of combined artons forms a continuous sectrum instead of a discrete one, it is generall imossible to conserve 4- momentum when artons are coalesced into a hadron. At resent, we choose to conserve the three-momentum during coalescence and determine the hadron secies according to the flavor and invariant mass of coalescing artons [77]. For seudo-scalar and vector mesons with same flavor comosition, the meson with mass closer to the invariant mass of coalescing quark and antiquark air is formed. E.g., whether a π or a ρ is formed from the coalescence of a air of ū and d quarks deends on if the invariance mass of the quarks is closer to the π mass or the centroid of ρ mass. The same criterion alies to the formation of octet and deculet barons that have same flavor comosition. It is more comlicated to treat the formation robabilities of flavor-diagonal mesons such as π and η in the seudo-scalar meson octet, and ρ and ω in the vector meson octet. Neglecting the mixing of η meson with the s s state, we take the following aroach for these flavor-diagonal mesons within the SU(2) flavor sace. For π formation from a uū or d d air, the robabilit P π is determined from the average of the numbers of formed π + and π mesons divided b the total number of uū and d d airs. Thus the total number of π, n π, is determined b aling the robabilit P π to each uū or d d air. The robabilit P ρ for forming a ρ meson from a uū or d d air is determined b a similar rocedure. After sorting all uū or d d airs according to their invariant masses, the lightest n π airs are assigned as π mesons. The rest of uū or d d airs form ρ mesons according to the robabilit of P ρ /(1 P π ), and the remaining airs form ω and η mesons with equal robabilities. The above quark coalescence model includes the formation of all mesons and barons listed in the HIJING rogram [54] excet η, Σ and Ξ, which are not resent in our hadronic transort model as well as K S and K L states. The resulting hadrons are given an additional formation time of.7 fm/c in their rest frame before the are allowed to scatter with other hadrons during the hadron cascade. As artons freeze out dnamicall at different times in the arton cascade, hadron formation from their coalescence thus occurs at different times, leading to the aearance of a coexisting hase of artons and hadrons during hadronization. D. Hadron cascade In the AMPT model, the following hadrons with all ossible charges are exlicitl included: π, ρ, ω, η, K, K, and φ for mesons; N,, N (144), N (1535), Λ, Σ, Ξ, and Ω for barons and corresonding antibarons. Man other higher resonances are taken into account imlicitl as intermediate states in scatterings between the above articles [14, 25]. Interactions among these

7 7 hadrons and corresonding inverse reactions are included as discussed in following subsections. 1. The ART model Hadron cascade in the AMPT model is based on the ART model [14, 78], which is a relativistic transort model originall develoed for heav ion collisions at AGS energies. The ART model includes baron-baron, baron-meson, and meson-meson elastic and inelastic scatterings. It treats exlicitl the isosin degrees of freedom for most article secies and their interactions, making it suitable for studing isosin effects in heav ion collisions [79]. Since it includes mean-field otentials for nucleons and kaons, the ART model can also be used for studing the effect due to the hadronic equation of state. Resonances such as ρ and are formed from ion-ion and ion-nucleon scattering, resectivel, with cross sections given b the standard Breit-Wigner form, and the also deca according to their resective widths. In all calculations resented in this stud, the masses and widths of resonances are taken to be their values in the vacuum, i.e., effects due to ossible modifications in dense hadronic matter [8] are neglected. Also, we have turned off the otentials in the AMPT model as their effects are much less imortant than scatterings in high energ heav ion collisions such as at SPS and RHIC. For baron-baron scatterings, the ART model includes the following inelastic channels: NN N( N ), NN ( N (144)), NN NN(πρω), (N ) NN, and N (144) NN (1535). In the above, N denotes either N (144) or N (1535), and the smbol ( N ) denotes a or an N, Also included are reaction channels relevant for kaon roduction, i.e., (N N )(N N ) (N )(ΛΣ)K. Details on their cross sections and the momentum deendence of resonance widths can be found in the original ART model [14]. For meson-baron scatterings, the ART model includes the following reaction channels for the formation and deca of resonances: πn ( N (144) N (1535)), and ηn N (1535). There are also elastic scatterings such as (πρ)(n N ) (πρ)(n N ). As an examle, the cross section for the elastic scattering of π N is evaluated b including heavier baron resonances with masses u to 2. GeV/c 2 as intermediate states using the Breit- Wigner form but neglecting interferences between the amlitudes from different resonances [14]. The ART model further includes inelastic reaction channels such as πn (πρη) and kaon roduction channels such as (πρωη)(n N ) K(ΛΣ). Kaon elastic scatterings with nucleons and baron resonances are included with a constant cross section of 1 mb [14]. Antikaon elastic scatterings with nucleons and inelastic channels, such as K(N N ) π(λσ), are included [81] using arameterized exerimental data [82]. Also included are kaon roduction channels involving three-bod final states, (πρω)(n N ) K KN [81]. Because of the difficult associated with the three-bod kinematics, the inverse kaon annihilation reactions of the above channels are neglected. For meson-meson interactions, the ART model includes both elastic and inelastic ππ interactions, with the elastic cross section consisting of ρ meson formation and the remaining art treated as elastic scattering. Kaon roduction from inelastic scatterings of light mesons is included via the the reactions (πη)(πη) K K and (ρω)(ρω) K K. Kaon or antikaon elastic scatterings with mesons in the SU(2) multilets excet the ion are included using a constant cross section of 1 mb [14], while the kaon-ion elastic scattering is modeled through the K resonance [42]. 2. Exlicit inclusion of K mesons The original ART model [14] includes the K resonance imlicitl through elastic πk scattering with the standard Breit-Wigner form for the cross section [83], i.e., σ πk = 6 mb/ [ 1 + 4( ] s m K ) 2 /Γ 2 K. Since the K meson not onl enhances elastic scattering between ion and kaon but also adds to strange article roduction through its addition to the strangeness degenerac, which becomes imortant when the hadronic matter is highl excited, K and K are included exlicitl in the hadronic hase of the AMPT model [42]. In addition to its formation from πk scattering and its deca, elastic scatterings of K with (ρωη) are included using same constant cross section of 1 mb as those for kaons. Inelastic reaction channels of (πη)(ρω) K K or K K and πk K (ρω) are also included [84]. 3. Baron-antibaron annihilation and roduction In heav ion collisions at or above SPS energies, antibaron roduction becomes significant and needs to be treated exlicitl during hadron cascade. The AMPT model initiall onl includes N N annihilation [4]. It was later extended to include (N N )( N N ) annihilation and also the inverse reactions of baron-antibaron air roduction from mesons [42]. The total cross section for annihilation is known emiricall, and the data has been arameterized as [85] σ = 67 mb/.7 lab, (13) where lab in GeV/c is the roton momentum in the rest frame of the antiroton. Following Ref. [14], a maximum cross section of 4 mb is imosed at low lab. Using hase-sace considerations [86], the branching ratios of annihilation to different multi-ion states are determined according to [86] M n ( s) = C [ 1 6π 2 ( s m π ) 3 ] n (4n 4)!(2n 1) (2n 1)! 2 (3n 4)!,(14)

8 8 where s is the center-of-mass energ of the roton and antiroton, and n is the number of ions in the final state. The dominant final states at moderate energies then involve several ions. For examle, the branching ratios at lab = 4 GeV/c are.33,.161,.36,.286,.151,.49 and.11, resectivel, for n from 3 to 9. For baron-antibaron annihilation channels involving baron resonances or N, their annihilation cross sections and branching ratios are taken to be the same as for annihilation at same center-of-mass energ. To include the inverse reactions that roduce baronantibaron airs during hadron cascade, which currentl onl treats scattering of two articles, we have further assumed that the final state of three ions is equivalent to a πρ state; the four-ion final state is equivalent to ρρ and πω with equal robabilities; the five-ion state is equivalent to ρω; and the six-ion state is equivalent to ωω. The cross sections for baron-antibaron air roduction from two mesons are then obtained from detailed balance relations. As shown later in the aer (Fig. 18), the above aroximate treatment of antibaron annihilation and roduction via two-meson states gives a satisfactor descrition of measured antiroton ield in central Pb+Pb collisions at SPS. σ (mb) Multistrange baron roduction from strangeness-exchange reactions KN Λπ KN Σπ KΛ Ξπ KΣ Ξπ KΞ Ωπ s s (GeV) FIG. 4: (Color online) Isosin-averaged cross sections for (multi)strange baron roduction as functions of center-ofmass energ of interacting antikaon and nucleon/heron. Productions of multistrange barons such as Ξ and Ω are included in the AMPT model [44] through the strangeness-exchange reactions such as K(ΛΣ) πξ, KΞ πω, (15) Since there is no exerimental information on their cross sections, we have assumed that the matrix elements for K(ΛΣ) πξ and KΞ πω are the same as that for the reaction KN πσ [87] at the same amount of energ above corresonding thresholds. The isosin-averaged cross section for the reaction KN πσ can be related to the cross sections for the reactions K Σ π and K n Σ π, which are known emiricall and have been arameterized in Ref. [82], b σ KN Σπ = 3 2 (σ K Σ π + σ K n Σ π ). (16) In Fig. 4, this cross section as well as the isosin-averaged cross sections for other multistrange baron roduction reactions are shown as functions of the center-of-mass energ above the threshold values of the interacting antikaon and nucleon/heron. The cross sections for the inverse multistrange baron destruction reactions are then determined b detailed balance relations. We note that these cross sections are comarable to those redicted b the couled-channel calculations based on the SU(3) invariant hadronic Lagrangian with emirical masses and couling constants [88]. With these cross sections, the ART transort model is able to describe the measured Ξ roduction in heav ion collisions at the AGS energies [89]. We note that, for strange barons Λ, Ξ, Ω and their anti-articles, onl their interactions with mesons have been included, while their annihilations b barons have not been included in the AMPT model at resent. 5. φ meson roduction and scattering The AMPT model also includes φ meson formation from and deca to kaon-antikaon air with the formation cross section given b the standard Breit-Wigner form [47]. Inelastic scatterings of the hi meson include baron-baron channels, (N N )(N N ) φnn, and meson-baron channels, (πρ)(n N ) φ(n N ), where the cross sections for the forward-going reactions are taken from the one-boson-exchange model [9]. The meson-baron channels also include K(ΛΣ) φn with cross section taken from a kaon-exchange model [91]. Phi meson scatterings with mesons such as π, ρ, K and φ have been studied before, and the total collisional width was found to be less than 35 MeV/c 2 [92]. A recent calculation based on the Hidden Local Smmetr Lagrangian [93] shows, however, that the collisional rates of φ with seudo-scalar (π, K) and vector (ρ, ω, K, φ) mesons are areciabl larger. Assuming that the matrix elements are indeendent of center-of-mass energ, we have included all these ossible reactions, i.e., φ(πρω) (KK )( K K ), and φ(kk ) (πρω)(kk ),

9 9 with cross sections determined from the artial collisional widths given in Ref. [93]. The cross section for the elastic scattering of the φ meson with a nucleon is set to 8 mb while the φ meson elastic cross section with a meson is set to 5 mb. The value of 8 mb is the φn total cross section [91] estimated from the φ meson hotoroduction data [94] and thus reresents the uer bound on the φ meson elastic cross section with a nucleon; quark counting then gives 5 mb as the uer bound on the φ meson elastic cross section with a meson. Note that, in most calculations of our revious stud on φ meson roductions at SPS and RHIC [47], the elastic cross section for φ meson scattering with a nucleon was taken to be.56 mb as extracted in Ref. [95] using the vector meson dominance model and the older φ meson hotoroduction data, while based on results of Ref. [93] the φk elastic cross section was extracted to be about 2 mb [47], which was then used as the φ meson elastic cross section with other mesons [96]. 6. Other extensions Other extensions of the ART model have also been made in the hadronic hase of the AMPT model. Antibaron resonances such as and N have been included exlicitl with their formations, decas, and scatterings analogous to those of baron resonances [44, 48]. Also, inelastic meson-meson collisions such as ππ ρρ have been added, and elastic scatterings between π and (ρωη) have been included with cross sections taken to be 2 mb. To address chemical equilibration of η mesons, which affects the height (or the λ arameter) of the correlation functions in two-ion interferometr, inelastic scatterings of η meson with other mesons have also been included with a constant cross section of 5 mb [48], which is roughl in line with recent theoretical redictions based on the Hidden Local Smmetr Lagrangian [97]. III. RESULTS FOR AND COLLISIONS The default AMPT results with no ocorn mechanism (see discussions on the ocorn mechanism in Sec. IV B) for and collisions are essentiall the same as the results from the HIJING model. In this section, we comare the results from the default AMPT model with or without the ocorn mechanism against available data from and collisions [98], where HIJING values for the a and b arameters, a =.5 and b =.9 GeV 2, are used in Eq. (9) for the Lund fragmentation function. A. Raidit distributions In Fig. 5, the AMPT results on charged article seudoraidit distribution are comared with the UA5 data for collisions at s = 2 GeV [99]. Since the AMPT results given b the solid or the long-dashed curves reresent all inelastic events with no trigger conditions, the should be comared with the UA5 inelastic data. The AMPT non-single-diffractive (NSD) results have included the UA5 NSD trigger b requiring events to have at least one charged article each in both ends of the seudoraidit intervals 2 < η < 5.6 [99]. We see that the AMPT model with or without the ocorn mechanism agrees reasonabl with both NSD and inelastic data. The kaon raidit distribution from the AMPT model is comared with the UA5 NSD data in Fig. 6, and both results with or without the ocorn mechanism also agree reasonabl with the data. dn ch /dη UA5 inelastic AMPT AMPT, no ocorn UA5 NSD AMPT NSD AMPT NSD, no ocorn η FIG. 5: (Color online) Pseudoraidit distributions of charged articles for collisions at s = 2 GeV. For collisions at P lab = 4 GeV/c, the raidit distributions of ions, kaons, rotons and antirotons are shown, resectivel, in Fig. 7, Fig. 8, Fig. 9, and Fig. 1. Curves with circles are measured cross sections from the LEBC-EHS Collaboration [1], while for comarison the AMPT results have been scaled u b the inelastic cross section at this energ (σ inel = 32 mb). We see that similar descritions of the charged article data are obtained with and without the ocorn mechanism. Including the ocorn mechanism gives, however, a better descrition of measured ion, kaon and antiroton ields and the shae of the roton raidit distribution. We have thus included the ocorn mechanism in all AMPT calculations. B. Transverse momentum sectra The transverse momentum sectra of charged ions and rotons for collisions at P lab = 24 GeV/c from

10 1 dn/d.2.1 (K + +K )/2 UA5 NSD AMPT NSD AMPT NSD, no ocorn dn/d x σ inel K K + K+ dσ/d data K dσ/d data AMPT AMPT, no ocorn CMS FIG. 6: (Color online) Same as Fig. 5 for kaons. FIG. 8: (Color online) Same as Fig. 7 for kaons. dn/d x σ inel π π + π + dσ/d data π dσ/d data AMPT AMPT, no ocorn dn/d x σ inel dσ/d data AMPT AMPT, no ocorn FIG. 7: (Color online) Raidit distributions of ions for collisions at P lab = 4 GeV/c. Circles are data from the LEBC-EHS Collaboration [1] CMS FIG. 9: (Color online) Same as Fig. 7 for rotons. the AMPT model are shown in Fig. 11. The are seen to reroduce reasonabl well the exerimental data from Ref. [11]. At the Tevatron energ of s = 1.8 TeV, results from the AMPT model for the transverse momentum sectra of ions, kaons and antirotons are shown in Fig. 12 and comared with data from the E735 Collaboration [12]. The AMPT NSD results have included the trigger b requiring events to have at least one charged article each in both ends of the seudoraidit intervals 3 < η < 4.5 [13]. We see that measured momentum sectra excet for antirotons are reroduced reasonabl well b the AMPT model. Note that the E735 T sectrum data shown in Fig. 12 have been averaged over raidit from weighing each track b the raidit range of the sectrometer [12], thus the are for d 2 N/(2π T d T d) instead for d 2 N/(2π T d T dη) even though the E735 sectrometer covers the accetance of.36 < η < 1..

11 dn/d x σ inel 1.5 _ dσ/d data AMPT AMPT, no ocorn d 2 N/(2π T d T d) (c 2 /GeV 2 ) K + +K _ (/5) π + +π E735 NSD AMPT NSD CMS T (GeV/c) dn/(2π T d T ) (full hase sace) (c 2 /GeV 2 ) FIG. 1: (Color online) Same as Fig. 7 for antirotons π (/1) data AMPT π + (x1) T (GeV/c) FIG. 11: (Color online) Transverse momentum sectra of charged ions and rotons for collisions at P lab = 24 GeV/c with data from Ref. [11]. C. Energ deendence Fig. 13 shows the energ deendence of dn ch /dη at η = for and collisions together with UA5 and CDF data at the Tevatron [99, 14]. The AMPT NSD results above s = 1 TeV have included the CDF NSD trigger b requiring events to have at least one charged article each in both ends of the seudoraidit intervals of 3.2 < η < 5.9 [14], and the AMPT NSD results FIG. 12: (Color online) Transverse momentum sectra of charged ions, charged kaons and antirotons from collisions at s = 1.8 TeV. below s = 1 TeV have included the UA5 NSD trigger b requiring events to have at least one charged article each in both ends of the seudoraidit intervals of 2 < η < 5.6 [99]. The agreement with the Tevatron data is reasonable. dn ch /dη (η=) UA5 NSD CDF NSD AMPT NSD s (GeV) FIG. 13: Energ deendence of dn ch /dη at midseudoraidit for and collisions. The energ deendence of the full hase sace K + /π + ratio for collisions is shown in Fig. 14 together with

12 12 the data comiled in Fig. 7 of Ref. [15]. It is seen that the AMPT model reroduces the data reasonabl well. K + /π + (full hase sace).1.5 data AMPT s (GeV) FIG. 14: Energ deendence of the K + /π + ratio in full hase sace for collisions. The energ deendence of the mean transverse momenta of ions, kaons and antirotons are shown in Fig. 15, where filled circles reresent the NSD data from the E735 Collaboration [16] and oen circles reresent data from the CERN Intersecting Storage Rings (ISR) [17]. Note that the AMPT NSD results shown in Fig. 15 have included the E735 NSD trigger [13] b requiring events to have at least one charged article each in both ends of the seudoraidit intervals 3 < η < 4.5. However, the trigger condition for the CERN ISR data is different. From the comarison with the AMPT inelastic results (solid curves), we see that NSD triggers have larger effects at lower energies. IV. RESULTS AT SPS ENERGIES To make redictions for heav ion collisions at RHIC, we first use the AMPT model to stud heav ion collisions at SPS. In articular, arameters in the AMPT model are determined b fitting the exerimental data from central Pb+Pb collisions at the laborator energ of 158A GeV, corresonding to a center-of-mass energ of about s NN = 17.3 GeV. A. Raidit distributions In Fig. 16, we show the raidit distributions of negativel charged articles (uer left anel), net-rotons and antirotons (uer right anel), charged ions (lower <P T > at mid raidit (GeV/c) E735 data CERN ISR data AMPT AMPT w. E735 NSD _ (π + +π )/ s (GeV) (K + +K )/2 FIG. 15: (Color online) Energ-deendence of the mean transverse momenta of ions, kaons and antirotons for and collisions. dn/d dn/d AMPT using HIJING a & b h π + NA49 (5%) dn/d NA49 (5%) NA49 rel π K + 3 dn/d 2 1 _ K FIG. 16: (Color online) Raidit distributions in central (b 3 fm) Pb+Pb collisions at s NN = 17.3 GeV. Circles and squares are exerimental data, while dashed curves are results from the AMPT model using default a and b arameters as in the HIJING model. left anel), and charged kaons (lower right anel) in central (b 3 fm) Pb+Pb collisions at SPS, obtained from the AMPT model with the default values of a =.5 and b =.9 GeV 2 in the HIJING model for the string fragmentation function. Comared with exerimental data for 5% most central Pb+Pb collisions from the NA49

13 13 Collaboration [15, 18, 19], it is seen that the AMPT model with these a and b arameters under-redicts both the negativel charged article [39] and kaon multilicities. dn/d dn/d AMPT using modifed a & b h π + NA49 (5%) dn/d NA49 (5%) NA49 rel π K + 3 dn/d 2 1 _ K FIG. 17: (Color online) Same as Fig. 16 but with solid curves for results from the AMPT model using modified a and b arameters in the Lund string fragmentation function. To increase the article multilicit in the AMPT model, we var the a and b arameters, and find that the exerimental data can be reasonabl reroduced with a = 2.2 and b =.5 GeV 2 as shown in Fig. 17. As seen from Eq. (9), a larger value of a corresonds to a softer fragmentation function, i.e., it leads to a smaller average transverse momentum for roduced hadrons and thus increases the article multilicit. The modified arameters also affect the string tension. According to Eq. (12), it is now 7% larger, and this leads to a strangeness suression factor of.33 instead of the default value of.3. Since the HIJING model with default a and b arameters reroduces the charged article multilicities in and collisions, the AMPT model with these arameters can robabl describe eriheral heav ion collisions. The a and b values in the AMPT model is thus exected to deend on the centralit of heav ion collisions, but this has not been studied. Since the robabilit for minijet roduction in the HI- JING model is ver small in collisions at SPS energies (onl about 4 minijet artons for a central Pb+Pb event on average), the artonic stage in the default AMPT model does not la an role for most observables in these collisions. Final-state hadronic scatterings are, however, imortant, and their effects are illustrated in Fig. 18. It is seen that final-state interactions reduce the ion and antiroton ields, but increase strangeness roduction, e.g., the sum of K and Λ as well as the sum of K + and dn/d dn/d AMPT AMPT, no FSI h π + NA49 (5%) dn/d NA49 (5%) NA49 rel π K + 3 dn/d 2 1 K FIG. 18: (Color online) Same as Fig. 17 but with dashed curves for results from the AMPT model without artonic and hadronic final-state interactions. Λ are both increased b about 2%. Overall, the kaon ield in the default AMPT model is larger than that in the HIJING model b about 4% due to combined effects of modified Lund string fragmentation and final-state interactions. B. Baron stoing and antiroton roduction Understanding baron stoing in relativistic heav ion collisions is imortant as it is related to the total energ deosited in the roduced hot dense matter during the collisions. The observed relativel large baron stoing in these collisions has led to the novel suggestion of gluon junction transort in an initial excited baron and its subsequent deca into a slowl moving baron and three leading beam mesons [56]. In the AMPT model, we have taken instead a more henomenological aroach to baron stoing in relativistic heav ion collisions b including the ocorn mechanism for baron-antibaron air roduction from string fragmentation. The ocorn mechanism introduces two additional baron roduction channels, i.e., the B B and BM B configurations in the Lund fragmentation model, which are controlled b two arameters in the JETSET/PYTHIA routine [59] used in the HIJING model. The first arameter MSTJ(12) is changed from 1 as set in the default HIJING model to 3 in the AMPT model [4, 42] in order to activate the ocorn mechanism, and the second arameter PARJ(5) controls the relative ercentage of the B B and BM B channels. We find that with equal robabilities for the B B and BM B configurations the net-baron raidit distribu-

14 14 tion at SPS can be reroduced as shown in Fig. 17 [42]. Without the ocorn mechanism, as in the default HI- JING model, the net-baron raidit distribution would eak at a larger raidit [39]. We also see from Fig. 18 that the antiroton ield at SPS is sensitive to the antibaron annihilation and roduction channels discussed in Sec. II D3. Without these reactions in the hadronic hase, the antiroton ield is too high comared to reliminar NA49 data. d 2 N/(m T dm T d) (.5<<.5) (arb.unit) C. Transverse momentum sectra K + π + NA44 data AMPT, no FSI AMPT m T m (GeV/c 2 ) FIG. 19: (Color online) Transverse momentum sectra at mid-raidit from the AMPT model with and without including hadron cascade comared with exerimental data from the NA44 collaboration for central Pb+Pb collisions at SPS. For the transverse momentum sectra, results from the default AMPT model are shown b solid curves in Fig. 19 for mid-raidit ions, kaons, and rotons in central (b 3 fm) Pb+Pb collisions at the SPS energ of snn = 17.3 GeV. Comared with exerimental data from NA44, given b solid diamonds, the AMPT model gives a reasonable descrition u to transverse mass of about 1 GeV/c 2 above the article mass. Without including rescatterings in the hadronic hase, the inverse sloe arameters for the transverse mass sectra of kaons and rotons from the AMPT model, given b dashed curves, are significantl reduced due to the absence of transverse flow that is induced b final-state interactions. dn ch /dη dn/d V. RESULTS AT RHIC ENERGIES A. Raidit distributions PHOBOS 13AGeV PHOBOS 2AGeV AMPT 13A GeV AMPT 2A GeV η π + BRAHMS π + BRAHMS π π dn/d dn/d 4 2 BRAHMS BRAHMS K + BRAHMS K + BRAHMS K FIG. 2: (Color online) Raidit distributions for Au+Au collisions at s NN = 13 and 2 GeV. Circles are the PHOBOS data for 6% most central collisions or the BRAHMS data for 5% most central collisions at s NN = 2 GeV, and curves are results from the default AMPT model for b 3 fm. Since the number of strings associated with soft interactions in the HIJING model deends weakl on the colliding energ, and the atomic number of Au is close to that of Pb, we use the same modified arameters in the Lund string fragmentation model used for Pb+Pb collisions at SPS energies for Au+Au collisions at RHIC energies. In Fig. 2, results for central (b 3 fm) Au+Au collisions at center-of-mass energies of s NN = 13 GeV (dashed curves) and 2 GeV (solid curves) are shown together with data from the PHOBOS Collaboration [11, 111, 112] and the BRAHMS Collaboration [113]. We find that measured total charged article seudoraidit distributions at both energies are roughl reroduced. More detailed comarisons on seudoraidit distributions at different centralities and different RHIC energies have been carried out b the BRAHMS Collaboration [114, 115], where results from the default AMPT model are comatible with the data. However, comared to central (to 5%) BRAHMS data, the AMPT model tends to over-redict the height of the raidit distributions of charged ions, kaons, rotons and antirotons. Note that the BRAHMS data on roton and antirotons at = in Fig. 2 have been corrected for feed-down from weak decas. Other models have also been used to stud hadron raidit distributions at RHIC. While results from the HI- K