Predictions for production of 3 H and 3 H in isobaric 96Ru+96Ru and 96Zr+96Zr Λ Λ 44 44 40 40 collisions at √sNN = 200 GeV Zhi-Lei She1,2, Gang Chen2,∗, Dai-Mei Zhou3, Liang Zheng2, Yi-Long Xie2, Hong-Ge Xu2 1 Institute of Geophysics and Geomatics, China University of Geosciences, Wuhan,430074, China. 2 School of Mathematics and Physics, China University of Geosciences, Wuhan,430074, China. 3 Institute of Particle Physics, Central China Normal University, Wuhan 430079, China. 3 3 3 3 3 3 The production of ΛH and ΛH, as well as H, H, He, and He are studied in central collisions of 96 96 96 96 isobars 44Ru+44Ru and 40Zr+40Zr at √sNN = 200 GeV, using the dynamically constrained phase- space coalescence model and the PACIAE model with chiral magnetic effect. The yield, yield ratio, coalescence parameters, and strangeness population factor of (anti-)hypertriton and (anti-)nuclei 96 96 96 96 produced in isobaric 44Ru+44Ru and 40Zr+40Zr collisions are predicted. The (anti-)hypertriton and (anti-)nuclei production is found to be insensitive to the chiral magnetic effects. Experimental data of Cu+Cu, Au+Au and Pb+Pb collisions from RHIC, LHC, and the results of PACIAE+DCPC model are presented in the results for comparison. I. INTRODUCTION sion systems, such as 63Cu+63Cu, 197Au+197Au, and 208Pb+208Pb collisions. One can see that there exists a Hypernuclei and their antihypernuclei are copiously gap of the system size for nucleus-nucleus interactions between 63Cu+63Cu and 197Au+197Au collisions. How- produced under conditions of extreme high temperatures 96 96 ever, the recent isobar program consisting of 44Ru+44Ru and energy densities in high-energy heavy ion collisions. 96 96 It creates a unique opportunity to study whether hyper- and 40Zr+40Zr collisions at the top RHIC energies of nuclei and antihypernuclei have the same behaviour and √sNN = 200 GeV, is favored to search the presence of to investigate the difference between light (anti-)nuclei Chiral Magnetic Effect (CME) [22–26], and it can also and (anti-)hypernuceli [1–3]. Hence it attracts a con- be used to fill the gap of collision system size discussed stant interest in studying antimatter and exploring fun- above. damental problems in physics, e.g., testing the funda- The CME effect can reveal some topological and elec- mental CPT theorem by precisely measuring the differ- tromagnetic properties of the quark gluon plasma(QGP) ence of the mass, lifetime and binding energy between in high-energy heavy ion collisions. Charge separation is 3 hypertriton (ΛH) and its corresponding anti-hypertriton an important consequence of the CME. Ma et al. [27–30] (3 H) [4, 5] in Au+Au and Pb+Pb collision systems. Λ introduced an additional CME-induced charge separa- The anti-hypertriton (3 H), the lightest bound antihy- tion to the initial conditions obtained from a multiphase Λ transport model(AMPT) [31], to study the CME-related pernucleus, consists of a antihyperon Λ , a antiproton physics. Refs. [27, 28] demonstrated that the final-state p, and a antineutron n, which had been discovered in interactions can reduce the charge separation in each col- Au+Au collisions at √sNN = 200 GeV by STAR Col- lision, while the relative difference of the CME signal laboration at the BNL Relativistic Heavy Ion Collider between the two isobaric collisions is insensitive to the (RHIC) [6] and then in Pb+Pb collisions at √sNN =2.76 final-state interactions. TeV by ALICE Collaboration at Large Hadron Collider (LHC) in CERN [7], respectively. The production of In this paper, the production of the final state hadrons, 3 H (3 H) has distinct features in heavy-ion collisions including p, p, Λ, and Λ, are simulated by the parton Λ Λ and hadron cascade model (PACIAE) [32], and an initial arXiv:2009.05402v2 [nucl-th] 20 Feb 2021 compared with corresponding normal three-body (anti- )nuclei 3He (3He) and 3H (3H), due to the different inter- three-flavor dipole charge separation [27] is introduced 96 96 96 96 action strength between hyperon-nucleon and nucleon- to simulate the CME, in 44Ru+44Ru and 40Zr+40Zr at nucleon [8]. Their detailed production mechanism is, the top RHIC energy of √sNN = 200 GeV with midra- pidity ( η < 0.5). Then, the dynamically constrained however, not fully understood. Hence the related the- | | oretical approaches on production of light (anti-)nuclei phase-space coalescence (DCPC) model [33] is applied to 3 3 3 3 and H ( H) have been carried out in the frameworks of study the production of ΛH( H) cluster in these two iso- Λ Λ Λ either the statistical thermal method [9–14] or the coa- baric collision systems. In this study, we expect to com- lescence model [15–20]. pare and investigate the production and properties of 3 H (3 H) in 96Ru+96Ru and 96Zr+96Zr collision systems The existence of 3 H (3 H) in heavy-ion reactions are Λ Λ 44 44 40 40 Λ Λ involving the chiral magnetic effect. observed, ranging from AGS [21] up to RHIC [6, 20] and LHC [7] collision energies, involving various colli- The paper is organized as follows: In sect. II, we pro- vide a concise introduction to the PACIAE with CME and DCPC model. Sec. III contains our numerical cal- culations for production and properties of 3 H (3 H). In Λ Λ ∗Corresponding Author: [email protected] sec. IV, a short summary is given. 2 II. MODELS dimension phase space by an integral d~qd~p The PACIAE model [32] is based on PYTHIA 6.4 [34] Y1 = , (2) h3 and is designed for various collision systems ranging from ZH6E proton induced reactions (p+p and p+A), to nuclear where H and E represent the Hamiltonian and energy reactions (A+A). Generally, this entire model has four of the particle, respectively. Similarly, the yield of N main physics stages composed of the parton initiation, particle cluster can also be computed using the following parton rescattering, hadronization, and hadron rescat- integral tering. At the first stage, the nucleus-nucleus collision is de- d~q1d~p1...d~qN d~pN YN = ... 3N . (3) composed into the nucleon-nucleon (NN) collisions ac- H6E h cording to the collision geometry and NN total cross Z Z section. The strings created in the NN collisions will In addition, equation (3) must meet the following con- break up into free partons leading to the formation of straint conditions the deconfined quark-gluon matter. After that, the de- composed partons interact with each other relied on the m0 6 minv 6 m0 + ∆m, (4) 2 2 LO-pQCD parton-parton cross sections [35]. Here, a→K factor is added to account for non-perturbative ~q 6 D , (i = j; i, j =1, 2,...,N). (5) QCD and higher-order corrections. Then, the hadroniza- | ij | 0 6 tion conducts via either the Lund string fragmenta- where tion model [34] or the phenomenological coalescence model [32]. The last step is the hadron rescattering pro- N 2 N 2 1/2 cess happening among the generated hadrons until the m = E ~p , (6) inv i − i hadronic freeze-out. (For more details see Ref. [32]). " i=1 i=1 # X X To study the CME-related physics, an additional CME-induced charge separation mechanism [27], which and Ei, ~pi(i=1,2,. ,N) are the energies and momenta of particles, respectively. m0 and D0 denote the rest mass switches py values for a fraction f of the downward mov- and diameter of light (anti-)nuclei or (anti-)hypernuclei. ing u(d) quarks with those of the upward moving u(d) Here, the radius values R = 1.74, 1.61, 5.0 fm are cho- quarks, is needed to introduce into the initial conditions sen for 3He (3He), 3H (3H), and 3 H (3 H) [19, 21, 45] in in the original PACIAE model [32]. The fraction f can be Λ Λ this simulation, respectively. ∆m represents the allowed described as mass uncertainty, and ~q presents the distance be- | ij | +(−) +(−) tween particles i-th and j-th. The integration in Eq. (3) N↑(↓) N↓(↑) should be replaced by the summation over discrete dis- f = − , (1) N +(−) + N +(−) tributions, as a coarse graining process in the transport ↑(↓) ↓(↑) model. where N denotes the number of a given quark, + and represent positive and negative charges, and III. RESULTS AND DISCUSSION are− the moving directions of quarks along the↑ y axis,↓ respectively. As the Ref. [28] mentioned, the initial At first, we can obtain the final state particles in cen- charge separation fractions f are different between iso- tral collisions of isobaric Ru+Ru and Zr+Zr using the bars 96Ru+96Ru and 96Zr+96Zr collisions since they 44 44 40 40 PACIAE model with CME. This simulation works on the have a same nucleon number but the 10% difference in assumption that (anti-)hyperons heavier than Λ (Λ) have proton number. In this work, we introduce an initial already decayed, and the model parameters are fixed on three-flavor quarks(u,d,s) charge separation into origi- the default values given in PYTHIA model, except the nal PACIAE model [32] to investigate the CME-related K factor and string fragmentation parameters parj(1), physics. parj(2), and parj(3). These selected parameters are con- The DCPC model [33] is developed to calculate produc- firmed by roughly fitting production of p (p) and Λ (Λ) tion of light (anti-)nuclei and (anti-)hypernuclei, after 96 96 96 96 in 0-15% 44Ru+44Ru and 40Zr+40Zr collisions to STAR PACIAE the final-state particles produced by model [32] data in 20-40% centrality Au+Au collisions at √sNN = in high energy collisions.