arxiv: v1 [hep-ph] 23 Jan 2018

Transcription

1 Baryon-antibaryon dynamics in relativistic heavy-ion collisions E. Seifert and W. Cassing Institut für Theoretische Physik, Universität Gießen, Germany (Dated: January 4, 8) arxiv:8.7557v [he-h] 3 Jan 8 The dynamics of baryon-antibaryon annihilation and reroduction (B B 3M) is studied within the Parton-Hadron-String Dynamics (PHSD) transort aroach for Pb+Pb and Au+Au collisions as a function of centrality from lower Suer Proton Synchrotron (SPS) u to Large Hadron Collider (LHC) energies on the basis of the quark rearrangement model (QRM). At Relativistic Heavy-Ion Collider (RHIC) energies we find a small net reduction of baryon-antibaryon (B B) airs while for the LHC energy of s NN =.76 GeV a small net enhancement is found relative to calculations without annihilation (and reroduction) channels. Accordingly, the sizeable difference between data and statistical calculations in Pb+Pb collisions at s NN =.76 TeV for roton and antiroton yields [9], where a deviation of.7 σ was claimed by the ALICE Collaboration, should not be attributed to a net antiroton annihilation. This is in line with the observation that no substantial deviation between the data and statistical hadronization model (SHM) calculations is seen for antihyerons, since according to the PHSD analysis the antihyerons should be modified by the same amount as antirotons. As the PHSD results for article ratios are in line with the ALICE data (within error bars) this might oint towards a deviation from statistical equilibrium in the hadronization (at least for rotons/antirotons). Furthermore, we find that the B B 3M reactions are more effective at lower SPS energies where a net suression for antirotons and antihyerons u to a factor of.5 can be extracted from the PHSD calculations for central Au+Au collisions. PACS: 4..-i; 4..Cn; 4..Jv; q; q I. INTRODUCTION Relativistic and ultra-relativistic heavy-ion collisions offer the unique ossibility to study a new hase of matter, i.e. a quark-gluon lasma (QGP), as well as ossibly the hase boundary between the hadronic and artonic hase. Lattice Quantum-Chromo-Dynamics (lqcd) calculations suggest that at vanishing baryon chemical otential (µ B =) there is a crossover hase transition from hadronic to artonic degrees of freedom [ 6] for the deconfinement hase transition as well as for the restoration of chiral symmetry. However, at some finite baryon chemical otential the crossover might turn to a firstorder hase transition imlying a critical endoint in the QCD hase diagram [7]. Since lattice calculations so far suffer from the fermion-sign roblem, no first rincile information on the hase boundary can be extracted from lqcd at large µ B, whereas at low µ B Taylor exansions of the thermodynamic otential (in owers of µ B /T ) rovide an alternative solution. The studies in Refs. [8, 9] show that for heavy-ion reactions at Relativistic Heavy- Ion Collider (RHIC) and Large Hadron Collider (LHC) energies the hase boundary is a crossover and the critical temerature for deconfinement T c is ractically the same as at µ B =. Due to the high energy densities reached in Au+Au (Pb+P collisions at RHIC and LHC energies as well as strong artonic interactions the final hadron yields turn out exerimentally to be close to thermal and chemical equilibrium as described by a grand-canonical ensemble of non-interacting hadronic states (with excluded volume corrections) [ 7]. In fact, the thermal analysis of hadron yields at midraidity show a high degree of thermalization [8], however, a sizeable difference between data and statistical calculations o u in Pb+Pb collisions at s NN =.76 TeV for roton and antiroton yields [9], where a deviation of.7 σ is obtained [8]. It has been argued in Refs. [9, ] that this deviation might be due to final state hadronic B B annihilation after chemical freezeout. On the other hand such a reduction was not seen in the relative yields of strange baryons/antibaryons to ions []. In Refs. [] the UrQMD transort model [, 3] has been emloyed as an hadronic afterburner to evaluate the final-state interactions and in articular the effects from B B annihilation after chemical freezeout, however, the backward channels had been discarded thus violating detailed balance [4]. This issue has been further addressed in Ref. [5] in a simlified model for the sace-time evolution but incororating detailed balance for the chemical reactions. In the latter study it was found that a net B B reduction by annihilation in central Pb+Pb collisions at s NN =.76 GeV of 4% might result thus coing aroximately with the exerimental observation in Ref. [9]. A more refined aroach - incororating detailed balance - has been roosed in Ref. [6] which solves chemical rate equations on to of + hydrodynamic evolution. At RHIC energies the authors reort a reduction of B B airs by about 5- %; results for LHC energies from this model are not known to the authors. Nevertheless, the imact of B B annihilation and reroduction by the inverse many-body channels should be calculated on a fully microscoic basis including detailed balance. A first ste in this direction has been taken in Ref. [7] where the three-body fusion of nonstrange seudoscalar and vector mesons to B B airs has been incororated in the Hadron-String Dynamics (HSD) transort aroach [8] that referentially describes the hadronic hase and rovides results close to the UrQMD transort

2 model [, 3] for SPS energies as demonstrated in Refs. [9, 3]. In Ref. [7] the matrix element squared for baryon-antibaryon annihilation has been extracted from the exerimental data on annihilation and the threebody meson channels have been determined on the basis of detailed balance. It was found that in central collisions of heavy nuclei at SPS energies the annihilation of antinucleons is almost comensated by the inverse recreation channels. A recent extension of the model has been resented in Ref. [3] within the Parton-Hadron- String Dynamics (PHSD) aroach [3 34] where the full strangeness sector has been included for the 3 reactions. The resulting model (denoted by PHSD4.) has been alied to central Pb+Pb collisions in the SPS energy regime and it was found again that B B annihilation and reroduction comensate each other to a large extend. We recall that the PHSD transort aroach [3 34] suerseeds the HSD aroach by a coule of asects that become essential with increasing bombarding energy: the formation of an initial artonic hase with quark and gluon quasiarticle roerties that are fitted to lattice QCD results in thermodynamic equilibrium A dynamical hadronization scheme on the basis of covariant transition rates Inclusion of further hadronic reactions in the strangeness sector with full baryon-antibaryon symmetry Inclusion of essential asects of chiral symmetry restoration in the hadronic hase [35]. Whereas the latter develoments are imortant for the lower SPS energy regime to account for the strangeness enhancement seen exerimentally in heavy-ion collisions, the formation of a artonic hase is mandatory to understand the hysics at higher SPS, RHIC and LHC energies. This has been demonstrated in a coule of PHSD studies in the ast for heavy-ion reactions from s NN = 4 GeV to.76 TeV [36 39]. Since multistrange baryons and antibaryons at to SPS energies and above no longer stem from string fragmentation (as in HSD [7]) but referentially from hadronization at energy densities around.5 GeV/fm 3 the issue of three-meson fusion reactions for the formation of baryon-antibaryon (B B) airs and the annihilation of B B airs to multile mesons has to be investigated (in addition to Ref. [3]) at RHIC and LHC energies. This work is organized as follows: In Sec. we recaitulate shortly the ingredients of PHSD and the quark rearrangement model for baryon-antibaryon annihilation and recreation (B B 3M) in the version 4. [3]. In Sec. 3 we resent results for antibaryons and multi-strange baryons from PHSD simulations for central Pb+Pb (Au+Au) collisions at RHIC and LHC energies in comarison to exerimental data and then focus on the centrality deendence of baryons and antibaryons. We will comare simulations using the baryon-antibaryon annihilation and formation only (B B) annihilation without the 3 channels. Global excitation functions for mesons, baryons and antibaryons will be rovided in Sec. 4 as well as excitation functions for the imact of final-state interactions and in articular the B B 3M reactions. We conclude our study with a summary in Sec. 5. II. REMINDER OF THE PHSD TRANSPORT APPROACH The Parton-Hadron-String Dynamics (PHSD) is a microscoic covariant transort aroach whose formulation is based on the Kadanoff-Baym equations [4 43] for Green s functions in hase-sace reresentation in first order gradient exansion beyond the quasiarticle aroximation [44]. The PHSD transort aroach describes in a consistent manner the whole time evolution of a relativistic heavy-ion collision as it incororates a hadronic and a artonic hase as well as dynamical transitions between the resective degrees-offreedom. The roerties of the quarks, antiquarks and gluons in the QGP hase are described by the Dynamical Quasi-Particle Model (DQPM) [45, 46], whose three arameters are fixed to reroduce the lqcd equation-ofstate at vanishing baryon chemical otential and which is based on effective roagators for the artons. The quarks and gluons have finite masses as well as widths that are given, resectively, by the real and imaginary arts of the retarded self energies resulting from twoarticle-irreducible diagrams of the effective full roagators. In PHSD simulations of nucleus-nucleus collisions color neutral strings (described by the FRITIOF Lund model [47]) are formed from the initial hard nucleon-nucleon scatterings. These strings decay into rehadrons with a formation time of.8 fm/c in which they do not interact. The string ends are identified with leading hadrons and may interact instantly with reduced cross sections according to the constituent quark model [8]. In case that the local energy density surasses the critical value of ɛ c.5 GeV/fm 3 the rehadrons dissolve into colored effective quarks, antiquarks and gluons given by the DQPM at the given local energy density. These artons then roagate in their self-generated mean field and interact via quasi elastic collisions between quarks, antiquarks and gluons. Additionally, q q airs may annihilate into a gluon and gluons may decay to q q airs. As the system exands the local energy density will dro close to or below ɛ c and will start to hadronize into off-shell mesons and baryons. During hadronization the energy, three-momentum, and quantum numbers are

3 3 conserved in each event [3]. In the hadronic hase the articles interact with each other (as in HSD) via elastic and inelastic collisions satisfying the detailed balance relations. The cross sections are taken from exeriments or effective models. A. Quark rearrangement model for B + B roduction and annihilation The quark rearrangement model (in the context of B B annihilation and reroduction) goes back to Ref. [7] and has been extended to the SU(3) flavor sector recently [3]. It is based on the exerimental observation of a dominant annihilation of into (on averag 5 ions at invariant energies.3 GeV s 4 GeV. Now the final number of 5 ions may be interreted as an initial annihilation into πρρ with the ρ mesons decaying subsequently into two ions each. The channel ππρ then leads to 4 final ions, the channel πωρ to 6 final ions, the channel ρωρ to 7 final ions etc. Accordingly, the baryon-antibaryon annihilation in the first ste is a 3 reaction with a conserved number of quarks and antiquarks. This is the basic assumtion of the quark rearrangement model which is illustrated in Fig. of Ref. [3]. By allowing the mesons M i to be any member of the or nonets one can describe an arbitrary B B annihilation and recreation by rearranging the quark and antiquark content, where B is a member of the baryon octet or deculet. This aroach gives a realistic descrition for annihilation and we assume that for other baryon-antibaryon airs than a similar annihilation attern holds. Since there are no measurements of annihilation cross sections other than n and this is our best guess at resent. B. 3 reactions in kinetic theory The treatment of 3 reactions can be incororated in the collision term in kinetic theory on a channel-bychannel basis by emloying detailed balance as formulated in Refs. [7, 3]. The matrix elements squared M c B B ( s) for a channel c, which - aart from hasesace integrals - determine the B B annihilation rate, are assumed to be the same for all flavor channels c. We recall that more than 5 individual mass channels are incororated and the 3-body hase-sace integrals have to be evaluated for each of these channels as a function of invariant energy s [3]. The technical solution to this roblem and detailed tests of the algorithm in a finite box with eriodic boundary conditions have been described in Ref. [3]. The matrix element squared M B B( s) has been extracted from the exerimental data on annihilation as a function of the relative momentum and it has been assumed that the roduct of relative velocity (in the center-of-mass system) v rel and annihilation cross section σ ann for other flavor channels is the same as for the channel. In this work we discard an exlicit strangeness suression factor which was found in Ref. [3] (Fig. 9) to have a minor imact on the actual results for relativistic heavy-ion collisions. Within these assumtions the kinetic aroach has no additional free arameter when incororating the 3 reactions in addition to the and channels [36]. The actual solution of the covariant transition rates in case of heavyion collisions is erformed by Monte-Carlo emloying the in-cell method as in Refs. [7, 3]. III. PHSD SIMULATIONS FOR HEAVY-ION COLLISIONS In this section we show the influence of the B B 3M reactions on heavy-ion collisions in articular at RHIC and LHC energies in extension of the calculations at SPS energies in Ref. [3]. Before coming to the actual results for hadron sectra we comare in Fig. the reaction rates for the total baryon-antibaryon annihilation and formation from PHSD in 5% central Pb+Pb (Au+Au) collisions at s NN = GeV (, and.76 TeV ( integrated over raidity. The solid blue lines denote the rates for B B annihilation when discarding the reroduction channels; the red solid lines stand for the B B annihilation rate when including the backward channels whereas the dashed lines dislay the reroduction rate in the latter case. Whereas the meson-fusion rate dominates at early times at the LHC energy over the B B annihilation rate ( the situation is inverse at the to RHIC energy. Without regeneration of B B airs (blue solid lines) the annihilation rates are lower than in case of B B reroduction which is, however, an unhysical limit and dislayed only for orientation. The exlicit deendence of ratios versus s NN will be discussed in Sec. IV. A. Hadron transverse mass sectra at RHIC and LHC energies We recall that raidity and transverse momentum sectra of antibaryons from PHSD in central Pb+Pb collisions at SPS energies have been shown in Ref. [3] in comarison to the available data. We here continue with PHSD results for antibaryons and mesons in 5% central Pb+Pb (Au+Au) collisions at the to RHIC energy ( s NN = GeV) and the LHC energy of s NN =.76 TeV. In Fig. we dislay the calculated transverse momentum sectra for rotons, ositive and negative ions as well as for kaons and antikaons in comarison to the data from the PHENIX Collaboration [48]. Whereas the hadron sectra are quite well described at lower transverse momenta there is a deficit at high T for all hadron secies in the PHSD calculations. We note that the hadron formation at the to RHIC energy at midraidity roceeds essentially by hadronization, i.e. by dynamical coalescence, which imlies that the quarks and antiquarks at hadronization have softer transverse momenta

4 4 dn/dt [/fm] Au+Au, s NN = GeV -5% central B B<->3M: --> B B<->3M: <-- B B ->3M: --> 8 Pb+Pb, s 6 NN =.76 TeV 4-5% central t [fm] Figure. (Color onlin The reaction rates of the B B 3M reactions (solid lin as a function of time in 5% central Au+Au collisions at s NN = GeV and Pb+Pb collisions at s NN =.76 TeV integrated over raidity. The solid blue lines denote the rates for B B annihilation when discarding the reroduction channels; the red solid lines stand for the B B annihilation rate when including the backward channels whereas the dashed lines dislay the reroduction rate in the latter case. in PHSD than in exeriment. The total hadron densities at midraidity are only marginally affected by the underestimated high T tail and we may conclude that the hadron densities within PHSD are sufficiently realistic such that rather solid results for the annihilation and fusion rates should emerge. The full red lines show the sectra from calculations with the 3 reactions included while the dashed lines corresond to calculations with the 3 reactions discarded. Since there are almost no differences between the lines we can conclude again that the 3 reactions have ractically no imact on baryon and meson sectra (cf. Ref. [3] for SPS energies). In Fig. 3 we show the same hadron T sectra at midraidity for 5% central Pb+Pb collisions at s NN =.76 TeV in comarison to the data from the ALICE Collaboration [49 5]. In this case the descrition of the data is rather good (excet for rotons) and again there is no visible imact of the 3 reactions on these transverse momentum sectra. We note in assing that the flow coefficients v n (for n =, 3, 4, 5) from PHSD for this system are also in a very good agreement with the exerimental measurements as shown in Ref. [53]. Thus we may state that the densities of the most abundant hadrons aear to be well under control in PHSD in articular at the LHC energy. We continue with the antibaryon transverse momentum /(π T )d N/(d T dy) [GeV - ] -3 Au+Au, s NN = GeV -5% central, y <.5 PHENIX w/ B B<->3M π+ -3 π - -3 K K T [GeV] Figure. (Color onlin The transverse momentum sectra for rotons, ositive and negative ions as well as for kaons and antikaons from PHSD at midraidity in comarison to the data from the PHENIX Collaboration [48] for 5% central Au+Au collisions at s NN = GeV. The full red lines show the results of calculations with the 3 reactions included while dashed lines corresond to calculations with the 3 reactions discarded. sectra at midraidity for to RHIC and LHC energies, which are dislayed in Figs. 4 and 5, resectively, in comarison with the data from the PHENIX, STAR and ALICE Collaborations [48, 5, 5, 54, 55]. Here again the low momentum sectra for, Λ + Σ, Ξ, and Ω + Ω + are roughly described at low momenta, however, the high T tails are missed considerably at s NN = GeV in Fig. 4 while they look somewhat better at the LHC energy in Fig. 5. We note that again there is no visible imact of the 3 reactions on these transverse momentum sectra which imlies that there is no sizeable

5 5 /(π T )d N/(d T dy) [GeV - ] Pb+Pb, s NN =.76 TeV -5% central, y <.5 ALICE w/o B B<->3M w/ B B<->3M 3 π + 3 π - K + K T [GeV] Figure 3. (Color onlin The transverse momentum sectra for rotons, ositive and negative ions as well as for kaons and antikaons from PHSD at midraidity in comarison to the data from the ALICE Collaboration [49 5] for 5% central Pb+Pb collisions at s NN =.76 TeV. The full red lines show the results of calculations with the 3 reactions included while dashed lines corresond to calculations with the 3 reactions discarded. /(π T )d N/(d T dy) [GeV - ] Au+Au, s NN = GeV -5% central, y <.5 PHENIX STAR w/o B B<->3M w/ B B<->3M Λ + Σ Ξ - Ω - + Ω T [GeV] Figure 4. (Color onlin The transverse momentum sectra for, Λ + Σ, Ξ, and Ω + Ω + from PHSD at midraidity in comarison to the data from the PHENIX and STAR Collaborations [48, 54, 55] for 5% central Au+Au collisions at snn = GeV. The full red lines show the results of calculations with the 3 reactions included while dashed lines corresond to calculations with the 3 reactions discarded. net B B annihilation at these ultra-relativistic energies. B. Centrality deendence of baryons and antibaryons at RHIC and LHC We continue with T integrated raidity densities for baryons and antibaryons as a function of centrality in terms of the number of articiating nucleons N art which is calculated within PHSD. Fig. 6 shows the raidity density of baryons and antibaryons from PHSD at midraidity in comarison to data from the PHENIX and STAR Collaborations [48, 54, 55] for 5% central Au+Au collisions at s NN = GeV. When discarding the 3 reactions (blue dotted lines) the exerimental data are slightly overestimated (excet for Λ + Σ ), while calculations with only B B annihilation (dashed red lines) show a slight tendency to underestimate the data. The results from PHSD calculations with the 3 reactions included are dislayed by the black solid lines and lie in between the other limits. This oints towards a small net B B annihilation at the to RHIC energy for

6 6 /(π T )d N/(d T dy) [GeV - ] Pb+Pb, s NN =.76 TeV -5% central, y <.5 ALICE w/o B B<->3M w/ B B<->3M Λ + Σ Ξ - Ω - + Ω + dn/dy y= 5 5 PHENIX STAR w/o B B<->3M w/ B B ->3M w/ B B<->3M 5 Au+Au, s NN = GeV 6 Λ+Σ Ξ Ω - + Ω N art g) Λ+ Σ N art f) T [GeV] Figure 5. (Color onlin The transverse momentum sectra for, Λ + Σ, Ξ, and Ω + Ω + from PHSD at midraidity in comarison to the data from the ALICE Collaboration [5, 5] for 5% central Pb+Pb collisions at s NN =.76 TeV. The full red lines show the results of calculations with the 3 reactions included while dashed lines corresond to calculations with the 3 reactions discarded. all baryons/antibaryons considered. We will quantify this net annihilation in Sec. IV. The situation is somewhat different at LHC energies. Fig. 7 shows the raidity density of baryons and antibaryons from PHSD at midraidity in comarison to data from the ALICE Collaboration [49 5] for 5% central Pb+Pb collisions at s NN =.76 TeV. The blue dotted lines dislay the calculated results when discarding the 3 reactions and the dashed red lines corresond to calculations with only B B annihilation. The results from PHSD calculations with the 3 reactions included are dislayed by the black solid lines and lie in all cases slightly Figure 6. (Color onlin The raidity density of baryons and antibaryons from PHSD at midraidity in comarison to data from the PHENIX and STAR Collaborations [48, 54, 55] for 5% central Au+Au collisions at s NN = GeV. The dashed red lines show the results of calculations with only B B annihilation, the solid black lines show results with the 3 reactions included while dotted blue lines corresond to calculations with the 3 reactions discarded. above the other limits indicating a net B B roduction at the LHC instead of an absortion. The calculations with only B B annihilation (red dashed lin underestimate the exerimental data (excet for Ξ and Ω ). In articular the,, Λ + Σ, and Λ + Σ multilicities are (within error bars) in line with exerimental observation at all centralities (when including the B B 3M channels) contrary to the results of the SHM quoted in Ref. [9]. On the other hand the Ξ, Ξ, Ω and Ω baryons are slightly overestimated in more central collisions when including the B B 3M channels. We attribute these results to a deviation from statistical equilibrium in the hadronization incororated in PHSD.

7 dn/dy y= ALICE w/o B B<->3M w/ B B ->3M w/ B B<->3M Λ+Σ Ξ - Ω - Pb+Pb, s NN =.76 TeV N art g) Λ+ Σ Ω N art Figure 7. (Color onlin The raidity density of baryons and antibaryons from PHSD at midraidity in comarison to the data from the ALICE Collaboration [49 5] for 5% central Pb+Pb collisions at s NN =.76 TeV. The dashed red lines show the results of calculations with only B B annihilation, the solid black lines show results with the 3 reactions included while dotted blue lines corresond to calculations with the 3 reactions discarded. IV. EXCITATION FUNCTIONS In this Sec. we will quantify the net effect of the B B 3M channels for central Pb+Pb (Au+Au) collisions as a function of the bombarding energy or s NN, resectively, including the revious results from Ref. [3]. A. Hadron yields at midraidity f) h) dn/dy y= 3-3 π + π - K + K - AGS SPS BES RHIC LHC Λ+Σ Ξ - Λ+ Σ w/ B B<->3M w/o B B<->3M s NN [GeV] Figure 8. (Color onlin The midraidity yields of mesons (, baryons ( and antibaryons ( from PHSD as a function of the invariant energy s NN for central heavy-ion collisions in comarison with the exerimental data taken from Refs. [48 5, 54 68]. The solid lines refer to calculations including the B B 3M channels while the dashed lines dislay calculations without these channels. The article yields from PHSD are connected by lines to draw the eye although exerimental data and calculations do not always corresond to the same centrality selection (and system) for different bombarding energies. 7 In Fig. 8 we first show the erformance of PHSD4. with resect to hadron roduction (at midraidity) in central Pb+Pb (Au+Au) collisions from s NN = 3.5 GeV to.76 TeV, i.e. by roughly 3 orders of magnitude in invariant energy. The solid lines refer to calculations including the B B 3M channels while the dashed lines dislay calculations without these channels. The arti- cle yields at midraidity from PHSD are connected by lines to draw the eye although exerimental data (taken from Refs. [48 5, 54 68]) and calculations do not always corresond to the same centrality selection (and system) for different bombarding energies. However, for given snn data and calculation corresond to the same centrality and collision system. From Fig. 8 we see that

8 8 [dn/dy y= (w/ <->3)] / [dn/dy y= (w/o <->3)] % central Pb+Pb/Au+Au Λ + Σ Ω + s NN [GeV] Figure 9. (Color onlin Ratios of -5% central midraidity yields from calculations with the full B B 3M reactions to calculations without them for the antibaryons as a function of the invariant energy s NN. [dn/dy y= (w/ <->3)] / [dn/dy y= (only ->3)] % central Pb+Pb/Au+Au Λ + Σ Ω + s NN [GeV] Figure. (Color onlin Ratios of -5% central midraidity yields from PHSD calculations with the full B B 3M reactions to calculations with only annihilation as a function of the invariant energy s NN. PHSD essentially reroduces the exerimental observations for ions, kaons and antikaons in the whole energy range. We recall that at AGS and SPS energies this is essentially due to the incororation of chiral symmetry restoration (cf. Refs. [35, 36]). The same holds true for the baryon and antibaryon excitation functions excet for the energy regime GeV < s NN < GeV where PHSD underestimates the baryons and antibaryons. The reason for this discreancy is resently not understood. However, by comaring the hadron yields from calculations with (solid lines) and without (dashed lines) the B B 3M channels we find no essential differences by eye. B. Quantitative imact of many-body reactions In this Subsec. we will quantify the effect of the B B 3M channels and B B 3M channels in 5% central Pb+Pb collisions for 3.5 GeV s NN.76 TeV. To this end we show in Fig. 9 the ratio of the antibaryons, λ + Σ, Ξ and Ω (at midraidity) from PHSD calculations including the B B 3M channels to calculations without them. At low s NN 3.5 we observe a sizeable net annihilation of antirotons and antihyerons by about a factor of two which is essentially due to the fact that here the nucleon density is very large comared to the antinucleon density. Practically the same holds for the strangeness S = ± sector while the net suression of Ξ is only %. For Ω s there is no net suression within error bars which results from the statistical errors of both calculations. With increasing invariant energy the net annihilation of antirotons and antihyerons disaears at s NN GeV, i.e. at the to SPS and lower RHIC energies. For s NN = 3 GeV and GeV we find a small net annihilation for, Λ+ Σ, and Ξ which turns to a small enhancement at the LHC energy of s NN =.76 TeV as quoted before. This is in contrast to the results of the model calculations in Refs. [, 5]. The small net suression of antirotons at the to RHIC energy, however, is in line with the results from Ref. [6] which also incororate detailed balance for the annihi-

9 9 lation channels. We interret the tiny enhancement of antibaryons at the LHC energy to result from the huge meson abundances which in hase-sace are slightly overoulated in PHSD relative to baryon-antibaryon airs at hadronization. In order to investigate the effect of the B B annihilation channels we show in Fig. the ratio of the antibaryons, Λ + Σ, Ξ and Ω (at midraidity) from PHSD calculations including the B B 3M channels to calculations with only the annihilation channels for the same reactions as in Fig. 9. Although this ratio is an unhysical quantity it allows to shed light on the relative imortance of the annihilation channels. For all antibaryons in Fig. this ratio is larger than unity which imlies that the back-reactions have some imact on the final antibaryon multilicities. This effect is most ronounced at lower SPS energies, where the baryon densities are large comared to the antibaryon densities, and dros below 5% enhancement for invariant energies above about GeV (within error bars). At to RHIC and LHC energies these modifications are below the % level since baryon and antibaryon densities are comarable and all elastic and inelastic channels are equal for time reversed states. Only the relative weight of baryons to mesons changes slightly. V. SUMMARY In this study we have emloyed the extended quark rearrangement model (QRM) for baryon-antibaryon annihilation (B B 3M) from Ref. [3] - incororated in PHSD4. - for the hadron roduction in heavy-ion collisions at ultra-relativistic energies. We recall - using simulations in a box with eriodic boundary conditions - that the numerical imlementation of the quark rearrangement model including the strangeness sector satisfies the detailed balance 3 relations on a channelby-channel basis as well as differentially as a function of the invariant energy s [3]. It has been found that the effects from the (B B 3M) reaction channels on the meson, baryon and antibaryon sectra is only moderate, although nonzero. At the to RHIC energy we find a small net suression of B B airs relative to calculations without these channels whereas at the LHC energy of s NN =.76 TeV there is even a net enhancement of B B airs which we attribute to the higher meson densities. The PHSD net antibaryon enhancement is in contrast to the results of the model calculations in Refs. [, 5] at the LHC energy whereas the small net suression of antirotons at the to RHIC energy is in line with the results from Ref. [6] which also incororate detailed balance for the annihilation channels. Accordingly, the sizeable difference between data and statistical calculations in Pb+Pb collisions at s NN =.76 TeV for roton and antiroton yields [9], where a deviation of.7 σ was obtained [8], should not be attributed to a net antiroton annihilation. On the other hand, no substantial deviation between the data and SHM calculations is seen for antihyerons [], which according to the PHSD analysis should be modified by the same amount as antirotons (cf. Fig. 9). To summarize our results for the LHC energy we show in Fig. the article ratios at midraidity from % central Pb+Pb collisions in comarison to the data from the ALICE Collaboration. Since the PHSD results are in line with data (within error bars) this oints towards a ossible deviation from statistical equilibrium in the hadronization (at least for rotons and antirotons). Yield dn/dy 3 π + π - ALICE w/ B B<->3M K + K - Λ Λ PHSD4. Pb+Pb, s NN =.76 TeV -% central, y <.5 Ξ - Ω - Ω + π + π - K + K - Λ Λ Ξ - Ω - Ω + Figure. (Color onlin The midraidity yields from PHSD4. including the B B 3M reactions (solid lines) for -% central Pb+Pb collisions at s NN =.76 TeV. The data oints are taken from Ref. [8]. Furthermore, we find that the B B 3M reactions are more imortant at lower SPS or FAIR/NICA energies where a net suression for antirotons and antihyerons u to a factor of.5 is seen in the PHSD calculations (cf. Fig. 9). In this energy regime further data on antibaryons (also with multile strangeness) will be needed with high statistics as a function of centrality and system size; a task well suited for the ucoming facilities at FAIR and NICA. So far the baryon-antibaryon dynamics is not sufficiently understood and oen for further uzzles. ACKNOWLEDGMENTS The authors acknowledge insiring discussions with E. L. Bratkovskaya, P. Moreau, A. Palmese, T. Steinert and R. Stock. We thank the Helmholtz International Center for FAIR (HIC for FAIR), the Helmholtz Graduate School for Hadron and Ion Research (HGS-HIR, and the Helmholtz Research School for Quark Matter Studies in Heavy-Ion Collisions (H-QM) for suort. The comutational resources have been rovided by the Center for Scientific Comuting (CSC) in the framework of the Landes-Offensive zur Entwicklung Wissenschaftlichökonomischer Exzellenz (LOEWE) and the Green IT Cube at FAIR.

10 [] Y. Aoki et al., Phys. Lett. B 643, 46 (6). [] S. Borsanyi et al., JHEP 9, 73 (); JHEP, 77 (); JHEP 8, 6 (). [3] S. Borsanyi et al., Phys. Lett B 73, 99 (4); Phys. Rev. D 9, 455 (5). [4] P. Petreczky [HotQCD Collaboration], PoS LATTICE, 69 (); AIP Conf. Proc. 5, 3 (3). [5] H.-T. Ding, F. Karsch, and S. Mukherjee, Int. J. Mod. Phys. E 4, 537 (5). [6] A. Bazavov et al., Phys. Rev. D 9, 9453 (4). [7] P. Senger et al., Lect. Notes Phys. 84, 68 (). [8] A. Bazavov et al., Phys. Rev. D 95, 5454 (7). [9] F. Karsch, PoS CPOD3, 46 (3). [] F. Becattini and R.Fries, Landolt-Boernstein 3, 8 (). [] Braun-Munzinger, K.Redlich, and J.Stachel, in Quark- Gluon Plasma 3, eds. H.C.Hwa and X.N.Wang, World Scientific 4,.49. [] F. Becattini, M. Gadzdicki, A. Keranen, J. Manninen, and R. Stock, Phys. Rev. C 69,495 (4). [3] A. Andronic et al., Nucl. Phys. A 77, 67 (6). [4] J. Rafelski and M. Danos, Phys. Lett. B 97, 67 (98). [5] A. Tounsi and K. Redlich, J. Phys. G 8, 95 (). [6] P. Braun-Munzinger, V. Koch, T. Schäfer, and J. Stachel, Phys. Ret. 6, 76 (6). [7] A. Andronic, D. Blaschke, P. Braun-Munzinger, J. Cleymans, K. Fukushima et al., Nucl. Phys. A 837, 65 (). [8] A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, e-print: arxiv:7.945 [9] ALICE Collaboration, B. Abelev et al., Phys. Rev. C 88, 449 (3). [] F. Becattini et al., Phys. Rev. C 85, 449 (); Phys. Rev. Lett., 83; Phys. Lett. B 764, 4 (7). [] ALICE Collaboration, J. Adam et al., Nature Phys. 3, 535 (7). [] S. A. Bass et al., Prog. Part. Nucl. Phys. 4, 55 (998). [3] M. Bleicher et al., J. Phys. G 5, 859 (999). [4] F. Becattini, M. Bleicher, J. Steinheimer, and R. Stock, arxiv: [5] Y. Pan and S. Pratt, Phys.Rev. C 89, 449 (4). [6] P. Huovinen and J. I. Kausta, Phys. Rev. C 69, 49 (4). [7] W. Cassing, Nucl. Phys. A 7, 68 (). [8] W. Cassing and E. L. Bratkovskaya, Phys. Re. 38, 65 (999). [9] H. Weber, E. L. Bratkovskaya, W. Cassing, and H. Stöcker, Phys. Rev. C 67, 494 (3). [3] E. L. Bratkovskaya, M. Bleicher, M. Reiter, S. Soff, H. Stöcker, M. van Leeuwen, S. A. Bass, and W. Cassing, Phys. Rev. C 69, 5497 (4). [3] E. Seifert and W. Cassing, arxiv:7.665 [3] W. Cassing and E. L. Bratkovskaya, Phys. Rev. C 78, 3499 (8). [33] W. Cassing and E. Bratkovskaya, Nucl. Phys. A 83, 5 (9). [34] E. L. Bratkovskaya, W. Cassing, V. P. Konchakovski and O. Linnyk, Nucl. Phys. A 856, 6 (). [35] W. Cassing, A. Palmese, P. Moreau, and E. L. Bratkovskaya, Phys. Rev. C 93, 49 (6). [36] A. Palmese, W. Cassing, E. Seifert, T. Steinert, P. Moreau and E. L. Bratkovskaya, Phys. Rev. C 94, 449 (6). [37] O. Linnyk, E. L. Bratkovskaya, and W. Cassing, Prog. Part. Nucl. Phys. 87, 5 (6). [38] O. Linnyk et al., Phys. Rev. C 88, 3494 (3); Phys. Rev. C 9, 5494 (5). [39] V. P. Konchakovski, E. L. Bratkovskaya, W. Cassing, V. D. Toneev, S. A. Voloshin and V. Voronyuk, Phys. Rev. C 85, 449 (). [4] L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, New York, 96). [4] J. S. Schwinger, J. Math. Phys., 47 (96). [4] L. V. Keldysh, Sov. Phys. JETP, 8 (965). [43] W. Botermans and R. Malfliet, Phys. Ret. 98, 5 (99). [44] W. Cassing and S. Juchem, Nucl. Phys. A 665, 377 (); Nucl. Phys. A 67, 47 (). [45] A. Peshier, Phys. Rev. D 7, 346 (4). [46] A. Peshier, J. Phys. G 3, S37 (5). [47] B. Nilsson-Almqvist and E. Stenlund, Comut. Phys. Commun. 43, 387 (987). [48] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 69, 3499 (4). [49] K. Aamodt et al. [ALICE Collaboration], Phys. Rev. Lett. 6, 33 (). [5] B. Abelev et al. [ALICE Collaboration], Phys. Lett. B 7, 5 (3); Phys. Rev. C 88, 449 (3); Phys. Rev. Lett., 3 (3); Phys. Lett. B 78, 6 (4). [5] E. Abbas et al. [ALICE Collaboration], Phys. Lett. B 76, 6 (3); Phys. Rev. Lett., 3 (3). [5] J. Adam et al. [ALICE Collaboration], Phys. Lett. B 754, 373 (6). [53] V. P. Konchakovski, W. Cassing, and V.D. Toneev, J. Phys. G 4, 556 (5). [54] G. Agakishiev et al. [STAR Collaboration], Phys. Rev. Lett. 8, 73 (). [55] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 98, 63 (7). [56] K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 89, 93 (); Phys. Rev. C 69, 494 (4). [57] C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 89, 93 (). [58] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 9, 83 (4). [59] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. C 96, no. 4, 4494 (7). [6] L. Ahle et al. [E8 Collaboration], Phys. Rev. Lett. 8, 65 (998). [6] L. Ahle et al. [E866 and E97 Collaborations], Phys. Lett. B 476, (); Phys. Lett. B 49, 53 (). [6] B. B. Back et al. [E97 Collaboration], Phys. Rev. Lett. 87, 43 (). [63] J. L. Klay et al. [E-895 Collaboration], Phys. Rev. C 68, 5495 (3). [64] P. Chung et al. [E895 Collaboration], Phys. Rev. Lett. 9, 3 (3). [65] S. V. Afanasiev et al. [NA49 Collaboration], Phys. Rev. C 66, 549 (). [66] C. Alt et al. [NA49 Collaboration], Phys. Rev. Lett. 94, 93 (5); Phys. Rev. C 73, 449 (6); Phys. Rev. C 77, 493 (8); Phys. Rev. C 78, 3498

11 (8). [67] T. Anticic et al. [NA49 Collaboration], Phys. Rev. Lett. 93, 3 (4); Phys. Rev. C 8, 3496 (9); Phys. Rev. C 83, 49 (). [68] M. K. Mitrovski et al. [NA49 Collaboration], J. Phys. G 3, S43 (6).