arXiv:2009.05402v2 [nucl-th] 20 Feb 2021 aoaina h N eaiitcHayInCollider at Ion collisions Pb+Pb Heavy in then Relativistic and [6] BNL (RHIC) the at laboration LC nCR 7,rsetvl.Tepouto of production The respectively. [7], Λ Collider CERN Hadron 3 in Large at (LHC) Collaboration ALICE by TeV ences ossso antihyperon a of consists pernucleus, bevd agn rmAS[1 pt HC[,20] colli- [6, various involving RHIC energies, to collision up [7] [21] LHC AGS and from ranging observed, (anti-)nuclei the- light related of the production and is, Hence on mechanism approaches production understood. oretical fully detailed not Their however, nucleon- and [8]. -nucleon nucleon between strength action at collisions Au+Au oprdwt orsodn omltrebd (anti- three-body normal )nuclei corresponding with compared ∗ coa- the or [9–14] [15–20]. method model lescence thermal statistical the either p neo h as ieieadbnigeeg between energy binding ( and lifetime ( differ- mass, hypertriton the funda- the measuring the of precisely testing ence by theorem e.g., CPT physics, con- mental fun- in a exploring and problems attracts antimatter damental it studying in Hence interest (anti-)nuclei stant [1–3]. light between (anti-)hypernuceli difference and and the behaviour investigate same hyper- to the whether have study antihypernuclei to and opportunity nuclei collisions. unique ion a heavy creates high-energy It in densities temperatures energy high and extreme of conditions under produced Λ 3 orsodn uhr [email protected] Author: Corresponding n a and , ( H )[,5 nA+uadP+bcliinsystems. collision Pb+Pb and Au+Au in 5] [4, H) h nihprrtn( anti-hypertriton The h xsec of existence The yence n hi niyence r copiously are antihypernuclei their and Hypernuclei Λ 3 3 Λ rdcin o rdcinof production for Predictions ( H )hsdsic etrsi ev-o collisions heavy-ion in features distinct has H) 3 e( He 3 Λ )hv encridoti h rmwrsof frameworks the in out carried been have H) 1 h-e She Zhi-Lei isobars aao uC,A+uadP+bcliin rmRI,LC and LHC, RHIC, from t collisions to Pb+Pb insensitive and be Au+Au Cu+Cu, to of found data is production (anti-)nuclei and olsec aaees n tagns ouainfacto isobaric population in strangeness produced and parameters, coalescence oe r rsne nterslsfrcomparison. for results the in presented are model pc olsec oe n the and model coalescence space 3 Λ 3 2 nttt fGohsc n emtc,CiaUiest of University China Geomatics, and Geophysics of Institute e and He) h rdcinof production The )adiscrepniganti-hypertriton corresponding its and H) colo ahmtc n hsc,CiaUiest fGeos of University China Physics, and Mathematics of School 3 .INTRODUCTION I. nttt fPril hsc,CnrlCiaNra Univer Normal China Central Physics, Particle of Institute 44 96 Λ 3 Ru+ ( H √ 3 n ( H s 1 hc a endsoee in discovered been had which , , Λ 3 NN Λ 3 2 44 96 ) h ihetbudantihy- bound lightest the H), agChen Gang , )i ev-o ecin are reactions heavy-ion in H) 3 uand Ru ) u otedffrn inter- different the to due H), 0 e ySA Col- STAR by GeV 200 = Λ 3 44 96 and H Ru+ 40 96 Zr+ olsosat collisions 44 96 a , Λ 2 3 Λ , uand Ru ,a elas well as H, 40 96 ∗ a-e Zhou Dai-Mei , rat Zr √ PACIAE s Λ 3 NN H √ 40 96 2 = s Zr+ NN and oe ihcia antceet h il,yedratio, yield yield, The effect. magnetic chiral with model 3 . H, 76 0 e,uigtednmclycntandphase- constrained dynamically the using GeV, 200 = 40 96 √ rcliin r rdce.Te(anti-)hypertriton The predicted. are collisions Zr 3 3 Λ s H, 3 H in Zheng Liang , NN including the to signal insensitive CME is the collisions interactions. of final-state isobaric two difference the relative between the col- each while in final-state separation lision, the charge the that reduce demonstrated can interactions 28] [27, CME-related Refs. the multiphase study physics. a to from [31], obtained model(AMPT) conditions separa- transport initial charge the to CME-induced tion additional an introduced hs-pc olsec ( coalescence phase-space √ iiy( pidity and td h rdcinof production the study nipratcneuneo h M.Ma CME. the is of separation consequence Charge important an collisions. ion heavy plasma(QGP) high-energy gluon in quark the of properties tromagnetic discussed size system collision of also gap can the it above. fill and to [22–26], used (CME) be Effect Magnetic Chiral e.I,asotsmayi given. is summary short a IV, sec. of properties and production com- Λ to the expect investigate we 3 study, and this pare In systems. collision baric ( model cascade hadron and insses uhas such systems, sion 208 ieacnieitouto othe to introduction concise a vide effect. magnetic chiral the involving a ftesse iefrncesncesinteractions nucleus-nucleus for size system between the of gap vr h eetioa rga ossigof consisting program isobar recent the ever, osmlt h M,in introduced CME, is the [27] simulate separation to charge dipole three-flavor h o HCeeg of energy RHIC top the uain o rdcinadpoete of properties and production for culations and 3 nisobaric in e and He, ( H s nti ae,tepouto ftefia tt hadrons, state final the of production the paper, this In h M ffc a eelsm oooia n elec- and topological some reveal can effect CME The h ae sognzda olw:I et I epro- we II, sect. In follows: as organized is paper The 0 GeV 200 = Pb+ NN DCPC 3 Λ 40 96 )in H) ecia antceet.Experimental effects. magnetic chiral he f(ni)yetio n (anti-)nuclei and (anti-)hypertriton of r 0 e,i aoe osac h rsneof presence the search to favored is GeV, 200 = Zr+ 208 | η 3 63 | eaesuidi eta olsosof collisions central in studied are He p esine,Whn407,China. Wuhan,430074, Geosciences, bcliin.Oecnseta hr xssa exists there that see can One collisions. Pb Cu+ 40 96 oe.Sc I otisornmrclcal- numerical our contains III Sec. model. < , 44 96 2 rcliin ttetpRI nrisof energies RHIC top the at collisions Zr p ine,Whn407,China. Wuhan,430074, ciences, iy ua 309 China. 430079, Wuhan sity, iLn Xie Yi-Long , Ru+ .) hn h yaial constrained dynamically the Then, 0.5). ,and Λ, , 63 44 96 h eut of results the uand Cu Ru+ 44 96 uand Ru 63 ,aesmltdb h parton the by simulated are Λ, Λ 3 44 96 √ DCPC Cu+ H( 197 44 96 uand Ru s 2 Ru+ NN Λ 3 PACIAE 40 96 Au+ ogG Xu Hong-Ge , PACIAE+DCPC )cutri hs w iso- two these in cluster H) 63 Zr+ oe 3]i ple to applied is [33] model ) 0 e ihmidra- with GeV 200 = Cu, 44 96 197 uand Ru 40 96 3] n ninitial an and [32], ) ucliin.How- collisions. Au PACIAE systems collision Zr 197 40 96 Zr+ Au+ 40 96 2 Λ 3 tal. et 44 96 197 ( H Zr+ ihCME with 40 96 Ru+ Zr u and Au, Λ 3 40 96 [27–30] ) In H). rat Zr 44 96 Ru 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 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-) 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. Previous works of (anti- 200 GeV, since their mean number of participating nu- )nuclei and (anti-)hypernuclei production in different cleons ( Npart ) are quite similar 140. Table I shows collision systems, e.g., pp [33, 36, 37], Cu+Cu [20, 38], the correspondingh i integrated yields∼ (dN/dy) of p and p Au+Au [39–42] and Pb+Pb [43, 44] interactions, have with η < 0.1 and 0.35 < p < 1.2 GeV/c, as well as | | T been studied using the same framework. Λ and Λ within η < 0.5 and 0.5 < pT < 8.0 GeV/c, According to the quantum statistical mechanics, one respectively. p and| |p take into account of contributions can estimate the yield of a single particle in the six- from primordial Λ decays. For comparison, the STAR 3 experimental data of 20-40% Au+Au collisions [46, 47]

are also presented. It can be seen from Tab. I that the 4

10 96 96 RuRu/ZrZr collisions at s =200 GeV yields of particles (p, p, Λ, and Λ) in 44Ru+44Ru col- NN

96 96 3

PAC IAE model lision are the same as those of Zr+ Zr collisions at 10 40 40 ) 2 c

Ru+Ru, 0-15% Zr+Zr, 0-15% the specified collision centrality. Moreover, the results 2 -2 of PACIAE model with CME for the two isobaric nuclear 10

1 p x collisions are well consistent with the measured STAR 10 data for Au+Au collisions. The same fitted parameters 0

p of K = 3.0, parj(1) = 0.13, parj(2) = 0.65, and parj(3) 10 dy) ((GeV) dy) T -1

= 0.44 are chosen for these two isobaric nuclear collision 10 dp systems. T

-2 p

10 (2

-3 /

10 N

TABLE I: The integrated yield dN/dy of particles (p, p, Λ 2 d 96 96 96 96 -4

10 and Λ) in 0-15% centrality 44Ru+44Ru and 40Zr+40Zr col- STAR data

NN PACIAE CME Au+Au, 20-40% lisions of √s = 200 GeV by with , as com- -5 pared to 20-40% Au+Au collisions in STAR experimental 10 (a) data [46, 47]. Here p (p) are inclusive of contributions from

1.4 primordial Λ (Λ) decays. (b)

1.2 Particle PACIAE STAR type Ru+Ru(0-15%) Zr+Zr(0-15%) Au+Au(20-40%) 1.0 Ru+Ru ratio Ru+Ru /

0.8 Npart 139.5 1.4 139.5 1.4 142.4 5.3 p ratio = 1

h i ± ± ± p p 11.13 0.03 11.12 0.03 11.85 1.15 0.6

± ± ± Zr+Zr p 9.60 0.02 9.59 0.01 9.33 0.91 ± ± ± 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Λ 5.57 0.02 5.55 0.01 5.70 0.55

p (GeV/c)

± ± ± T Λ 4.51 0.01 4.50 0.01 4.53 0.34 ± ± ± FIG. 1: (a) The transverse momentum spectrum of parti- 96 96 96 96 Figure 1(a) presents the transverse momentum distri- cles ( p, p, Λ, Λ) in midrapidity 44Ru+44Ru and 40Zr+40Zr butions of p (p) and Λ (Λ) (open symbols) in 0-15% collisions at √sNN = 200 GeV. The open symbols show the 96 96 96 96 44Ru+44Ru and 40Zr+40Zr collisions at √sNN = 200 results of PACIAE model with CME, and the solid symbols GeV calculated by the PACIAE model with CME. The show the results from STAR data [46, 47]. For clarity the spectra data are divided by powers of 10. (b) The yield ra- STAR experimental data for 20-40% Au+Au collisions 96 96 tios of particles (p, p, Λ, Λ) produced in 44Ru+44Ru collisions taken from Refs. [46, 47] are shown by the solid symbols. 96 96 It can be seen that the transverse momentum spectrum to that in 40Zr+40Zr collisions. of particles p (p) and Λ (Λ) simulated by PACIAE model with CME are compatible with the STAR data within uncertainties. Besides, figure 1(b) shows the distribution 96 96 TABLE II: The integrated yields dN/dy of (anti-)particles p of the invariant yield ratios of particles for 44Ru+44Ru 3 3 3 3 3 3 96 96 (p), Λ (Λ), and (anti-)nuclei ΛH (ΛH), He ( He), H ( H) and 40Zr+40Zr collisions at √sNN = 200 GeV, as a func- calculated by PACIAE+DCPC model with CME in 0-10% 96 96 96 96 tion of pT . Obviously, one can see from figure 1(b) that 44Ru+44Ru and 40Zr+40Zr collisions of √sNN = 200 GeV there is no significant difference for transverse momen- with η < 0.5. Here p (p) productions from Λ (Λ) feed down tum spectra of (anti-)particles between the two isobaric contribution| | is excluded in the coalescence procedure. nuclear collisions, except the fluctuation at higher pT . Nucleus type Ru+Ru(0-10%) Zr+Zr(0-10%) In the following, we generate 4.0 108 most central 96 96 96× 96 Npart 151.8 1.4 151.8 1.4 (0-10%) events for 44Ru+44Ru and 40Zr+40Zr collisions h i ± ± p 8.03 0.02 8.02 0.02 at √sNN = 200 GeV using the PACIAE model with CME, ± ± respectively. These (anti-)nucleons and (anti-)hyperons p 6.51 0.01 6.50 0.03 ± ± produced within PACIAE with CME are used as input Λ 7.04 0.01 7.03 0.01 ± ± for the DCPC model. The proton productions from Λ Λ 5.64 0.01 5.65 0.01 feed down contribution is excluded in the coalescence 3 −5 ± ± ΛH (10 ) 6.50 0.06 6.45 0.07 procedure. Then, we obtain the integrated yields dN/dy 3 −5 ± ± ΛH (10 ) 3.15 0.03 3.12 0.03 of light (anti-)nuclei and (anti-)hypertriton nuclei with − ± ± 3He (10 5) 8.57 0.09 8.38 0.08 η < 0.5 and pT < 3.0 GeV/c for the most central bin 3 −5 ± ± |of| 0-10%, respectively. Here we choose the parameter He (10 ) 4.41 0.04 4.32 0.04 3 3 3 −5 ± ± ∆m =1.53 MeV for ΛH ( H), and ∆m =2.13 MeV for H (10 ) 8.61 0.08 8.51 0.09 Λ − ± ± 3He (3He) and 3H (3H). 3H (10 5) 4.79 0.04 4.73 0.04 ± ± Table II presents the integrated yields dN/dy of (anti- 4

)hyperons and (anti-)hypertriton (Λ, Λ, 3 H, 3 H), as to particles (p/p, Λ/Λ, 3 H/3 H, 3He/3He, 3H/3H), and Λ Λ Λ Λ well as (anti-)nuclei (p, p, 3He,3 He, 3H,3 H) calculated the mixed ratios (Λ/p,Λ/p, 3 H/3He,3 H/3He, 3 H/3H, Λ Λ Λ by the PACIAE+DCPC model with CME in most cen- 3 3 ΛH/ H) calculated by PACIAE+DCPC model with CME tral (0-10%) 96Ru+96Ru and 96Zr+96Zr collisions at 96 96 96 96 44 44 40 40 in 44Ru+44Ru and 40Zr+40Zr collisions at √sNN = √sNN = 200 GeV, respectively. It can be seen that 200 GeV. One can see from Table III that the yield ra- the yields of (anti-)hypertriton, (anti-)tritium, and (anti- tios of 3 H/3 H, 3He/3He, and 3H/3H are the same within 96 96 96 96 Λ Λ )helium-3 nuclei in central 44Ru+44Ru and 40Zr+40Zr the error range, although their yields are not the same as collisions at √sNN = 200 GeV from PACIAE+DCPC sim- shown in Table II. And the mixed ratio values of (anti- − ulations, are all at the order of 10 5. However, the yields )hypernuclei to (anti-)nuclei(3 H/3He, 3 H/3He, 3 H/3H, Λ Λ Λ of (anti-)hypernuclei are less than that of correspond- 3 H/3H) are also the same in the range of uncertainty in ing (anti-)nuclei with the equal numbers. The Λ central isobaric Ru+Ru and Zr+Zr at √sNN = 200 GeV. yields of (anti-)hypernuclei and (anti-)nuclei in isobaric 96 96 96 96 Besides, the yield ratio results of antiparticles to parti- 44Ru+44Ru and 40Zr+40Zr collisions are the same within cles and their mixed ratios simulated by PACIAE+DCPC the range of uncertainty, indicating that the production model with CME are found to be in agreement with the of (anti-)hypernuclei and light (anti-)nuclei between the above theoretical interpretation within uncertainties. two isobaric nuclei-nuclei collisions is insensitive to the Fig. 2 and Table III show that the ratios of anti-nuclei difference in charge. to nuclei (3 H/3 H, 3He/3He, 3H/3H) are less than 1, In order to understand the fundamental properties of Λ Λ meaning that the yields of antiparticles is less than that antimatter in nuclear collisions, we provide a system- of corresponding particles; similarly, the mixed ratio val- atic investigation to the yield ratios of different (anti- ues indicate that the yields of the (anti-)hypertriton are )nuclei and (anti-)hypernuclei, which are deeply related less than that of (anti-)nuclei. Our simulation results are to the fractions of constituent nucleons in the naive coa- consistent with the STAR data of Cu+Cu [47–49] and lescence framework [11, 15]. For instance, the yield ratio Au+Au [6, 46, 47] collisions at √s = 200 GeV. of 3 H/3 H should be proportional to (p/p)(n/n)(Λ/Λ), NN Λ Λ which is approximate to (p/p)2(Λ/Λ), i.e,

CuCu RuRu(ZrZr) AuAu PbPb 3 H pnΛ p Λ Λ 2 - 1 3 3 3 10

H 3 = ( ) . (7) He H ΛH pnΛ ≃ p Λ 3 3 3

H He H

PACIAE+DCPC - 2

10

RuRu, ZrZr, 200 G eV, 0-10%

1.6

CuCu, AuAu, 200 G eV, 0-10%

s RuRu/ZrZr col l i sions at =200 GeV

NN

- 3 PbPb, 2.76 T eV, 0-10% ) 10

1.4 3

PACIAE+DCPC /c 2

Ru+Ru, 0-10% Zr+Zr, 0-10%

1.2

- 4

10 (GeV 3

1.0 B

3

He

ST AR Au+Au, 0-12% - 5

10

0.8 3

He ST AR Au+Au, 0-12%

3

H ALICE Pb+Pb, 2.76 T eV Ratio

0.6

- 6

10

0.4

60 80 100 120 140 160 180 200 220

STAR

Mass Number A

0.2

C u+C u, 0-10% Au+Au, 0-80%

Au+Au, 20-40% 0.0 3 3 FIG. 3: The coalescence parameters √B3 of ΛH(ΛH), as well 3 3 3 3 96 96 as H( H) and He( He) are compared in isobars 44Ru+44Ru

3 3 3 3 3 3 3 96 96 H H H H H p He H and 40Zr+40Zr collisions at √sNN = 200 GeV with PA-

3 3 3 3 3 3 3 p p H H p He He H He H CIAE+DCPC model with CME. It is also compared with the results from different collisions of Cu+Cu [38], Au+Au [42], and Pb+Pb [43] collisions. The open symbols denote the FIG. 2: The ratios and mixed ratios of (anti-)matter form results computed by PACIAE+DCPC. The solid points take PACIAE+DCPC model with CME (open symbols) in 0-10% from STAR [52] and ALICE [7]. The error bars show statis- Ru+Ru and Zr+Zr collisions, compared with Cu+Cu and tical uncertainties. Au+Au collisions. The data are taken from STAR [6, 46– 49]. The vertical lines and error boxes show statistical and systematic errors, respectively. In nuclear collisions, the invariant yields for produc- tion of (anti-)hypernuclei and (anti-)nuclei can be related Table III represents the yield ratios of antiparticles to the primordial yields of (anti-)nucleons in the coales- 5

96 96 96 96 TABLE III: The (anti-)nucleus ratios from the PACIAE+DCPC model with CME in central (0-10%) 44Ru+44Ru and 40Zr+40Zr collisions at √sNN = 200 GeV. The top section of the table shows the three ratios of anti-nucleus to nucleus, followed by the mixed ratios of (anti-)nucleus to (anti-)nucleus. The ratios between proton, anti-proton, hyperon, and anti-hyperon are shown at the bottom. STAR data are taken from Cu+Cu and Au+Au collisions at √sNN = 200 GeV [6, 46–49], respectively.

Nucleus STAR PACIAE+DCPC ratio Cu+Cu Au+Au Ru+Ru Zr+Zr 3 H/3 H 0.49 0.18 0.07 0.48 0.02 0.48 0.01 Λ Λ − ± ± ± ± 3He/3He 0.46 0.17 0.45 0.02 0.04 0.51 0.01 0.51 0.01 ± ± ± ± ± 3H/3H 0.56 0.01 0.56 0.01 − − ± ± 3 H/3He 0.89 0.28 0.13 0.71 0.02 0.72 0.01 Λ − ± ± ± ± 3 H/3He 0.82 0.16 0.12 0.76 0.02 0.77 0.01 Λ − ± ± ± ± 3 H/3H 0.66 0.02 0.66 0.01 Λ − − ± ± 3 H/3H 0.75 0.02 0.76 0.02 Λ − − ± ± p/p 0.80 0.04 0.79 0.11 0.81 0.01 0.81 0.01 ± ± ± ± Λ/Λ 0.82 0.12 0.80 0.10 0.80 0.01 0.80 0.01 ± ± ± ± Λ/p 0.84 0.09 0.88 0.01 0.88 0.01 ± − ± ± Λ/p 0.83 0.08 0.87 0.01 0.87 0.01 ± − ± ±

cence framework [50, 51] by Eq. (8) ALICE [7] are also presented in Fig. 3.

3 3 d NA d NP A

≈ CuCu RuRu(ZrZr) AuAu PbPb EA 3 BA(EP 3 ) , (8) d PA d PP

t t

s s s s

2.0

3 where Ed3N/d3p stands for the invariant yields of nucle- 3 3 3 ons or light (anti-)nuclei and (anti-)hypernuclei, and A PACIAE+DCPC

RuRu,ZrZr, 200G eV, 0-10% is the atomic mass number, respectively. BA represents

CuCu,AuAu, 200G eV, 0-10% the coalescence parameters, which relates to the freeze- 1.5 1−A PbPb, 2.76T eV, 0-10% out correlation volume, i.e., BA Vf . pA,pp denote )

∝ t 3

their momentum, with pA = App assumed. S (

3 1.0 Fig. 3 presents coalescence parameters √B3 of H 3

Λ S (3 H), as well as 3He (3He) and 3H (3H) in Ru+Ru Λ and Zr+Zr collisions. Meanwhile, the results from dif-

0.5

ferent collision systems of Cu+Cu [20], Au+Au [42], s

ST AR Au+Au, 0-80% and Pb+Pb [43] collisions are compared as a function 3 s ST AR Au+Au, 0-80%

3

√ s of atomic mass number A. One can see that B3 cal- ALICE Pb+Pb, 2.76T eV, 0-10%

3

PACIAE DCPC CME 0.0 culated by + model with is constant ALICE Pb+Pb, 2.76T eV, 0-10% s within the error range from Cu+Cu (A = 63) to Au+Au 3 (A = 197) collisions at top RHIC energy. The coales- 60 80 100 120 140 160 180 200 220 cence parameter is found to drop from RHIC energy to PbPb collisions at the LHC energy. This decreasing can Mass Number A be understood as the correlation volume of the QGP t fireball at LHC becomes larger, as indicated by the pion FIG. 4: Comparison of strangeness population factor s3(s3) 96 96 96 96 HBT measurement [42, 52]. in isobars 44Ru+44Ru and 40Zr+40Zr collisions at √sNN = Specifically, with respect to 0-10% Ru+Ru or Zr+Zr 200 GeV by PACIAE+DCPC model with CME. It is also collisions, the values of √B for nuclei 3 H, 3He, 3H and compared with the results of Cu+Cu [20], Au+Au [42], and 3 Λ Pb+Pb [43] collisions. The open symbols denote the results their corresponding anti-nuclei 3 H, 3He, 3H are about Λ computed by PACIAE+DCPC, and the solid points denote (3.63 0.46) 10−4, (3.90 0.42) 10−4, (3.79 0.39) data from STAR [6] and ALICE [7]. Error bars and error 10−4 ±and (3×.46 0.44) ±10−4,× (3.81 0.40)± 10−×4, boxes denote statistical and systematic errors, respectively. (3.75 0.41) 10± −4, respectively.× It is± clear that× the ± × value of √B3 of (anti-)hypertriton is the same as that of The strangeness population factor s3, should be about (anti-)nuclei within the uncertainties. One can also find one in the coalescence model for particle production, as that the negative (hyper-)nuclei are the same as that of first suggested in Ref. [21]. It is a possible tool to study positive (hyper-)nuclei. Meanwhile, the experiment data the nature of a quark-gluon plasma created in high en- of 0-12% Au+Au in STAR [52] and 0-10% Pb+Pb from ergy nuclear collisions [53], due to its sensitivity to the 6 local baryon-strangeness correlation [54, 55]. This factor with η < 0.5 and pT < 3.0 GeV/c, respectively. We typically is written as predict| | the yield, yield ratio, coalescence parameters, and strangeness population factor of (anti-)hypertriton 3 3 3 3 3 3 3 3 s3 = ( H p)/( He Λ), (9) ( H, H) and (anti-)nuclei ( H, H, He, He) in iso- Λ × × Λ Λ baric Ru+Ru and Zr+Zr collisions. Then the chiral 3 which can be straightforwardly extended to H expressed magnetic effect on the generation of (anti-)hypertriton as (3 H, 3 H) and (anti-)nuclei (3H, 3H, 3He, 3He) in high- Λ Λ energy collisions are studied. It is found that there is st = (3 H p)/(3H Λ). (10) 3 Λ × × no clear difference for the generation and properties of (anti-)hypertriton (3 H, 3 H) and light (anti-)nuclei (3H, In Fig. 4, we compare the values of strangeness Λ Λ 3 3 3 t t H, He, He) in isobaric Ru+Ru and Zr+Zr collision population factor s3(s3) and s3(s3) calculated by PA- CIAE+DCPC model with CME in central Ru+Ru and systems, although these two collision systems have dif- Zr+Zr collisions. Besides, the values of Cu+Cu [20], ferent CME. Au+Au [42], and Pb+Pb [43] collisions varying with In addition, the coalescence parameters √B3 and t mass number A are also presented in this figure. It is the strangeness population factor s3(s3) that produce t t 3 H(3 H), as well as 3H(3H) and 3He(3He) in 96Ru+96Ru shown that the values of s3(s3) and s (s ) for three- Λ Λ 44 44 3 3 96 96 body coalescence are all the same within the error range and 40Zr+40Zr collisions are compared with those of as A increases from 63 to 197 in central (0-10%) Cu+Cu, Cu+Cu, Au+Au, and Pb+Pb collisions at √sNN = 200 Ru+Ru(Zr+Zr) to Au+Au collisions at RHIC energy. GeV by the PACIAE + DCPC model. One can see that t t t the √B3 and s3(s ) are all the same within the error But the values of s3(s3) and s (s ) decrease in Pb+Pb 3 3 3 range as atomic mass number A increases from 63 to collisions at √sNN =2.76 TeV. This may be interpreted as the (anti-)Λ particles freezeout earlier than (anti- 197 in Cu+Cu, Ru+Ru(Zr+Zr) to Au+Au collisions at t t )nucleons but their relative freeze-out time is closer at RHIC energy. But the values of s3(s3) and s3(s3) de- LHC than at RHIC [18]. crease in Pb+Pb collisions at √sNN = 2.76 TeV. This t t may be interpreted as the Λ(Λ) particles freezeout earlier Numerically, the present values of s3, s3 and s3, s3 are 0.87 0.02, 0.82 0.02 and 0.94 0.02, 0.83 0.02 in than (anti-)nucleons but their relative freeze-out time is 0-10%± Ru+Ru collisions,± and 0.88 ±0.01, 0.83 ±0.01 and closer at LHC than at RHIC. The experimental data of 0.94 0.02, 0.83 0.02 in 0-10%± Zr+Zr collisions,± re- Cu+Cu, Au+Au and Pb+Pb collisions from RHIC, LHC ± ± t t are included in the comparison, which shows that our spectively. Meanwhile, the values of s3(s3) and s (s ) 3 3 simulation results for the √B and s (st ) of 3 H (3 H), shown in the Fig. 4 for Au+Au and Pb+Pb collisions by 3 3 3 Λ Λ 3 3 3 3 PACIAE+DCPC model are in agreement with the corre- He ( He), H ( H) are consistent with the experimental sponding available data from STAR [6] and ALICE [7] results within the error range. within uncertainties. ACKNOWLEDGMENT IV. CONCLUSION This work was supported by NSFC(11475149, 11775094, In this paper, we use the PACIAE model with CME 11905188), as well as supported by the high-performance and DCPC model to simulate production of 3 H and 3 H, Λ Λ computing platform of China University of Geosciences. 3 3 3 3 96 96 as well as H, H, He, and He in isobaric 44Ru+44Ru The authors thank Dr. Feng-Xian Liu for fruitful dis- 96 96 and 40Zr+40Zr central collisions at √sNN = 200 GeV cussions.

[1] J.H. Chen, D. Keane, Y.G. Ma et al., Phys. Rep. 760, 754, 360(2016). 1(2018). [8] E. Botta, T. Bressani, and G. Garbarino, Eur. Phys. J. [2] Y.G. Ma, J.H. Chen, and L. Xue, Front. Phys. 7, A 48, 41(2012). 637(2012). [9] V. Topor Pop and S. Das Gupta, Phys. Rev. C 81, [3] P. Braun-Munzinger and B. D¨onigus, Nucl. Phys. A 987, 054911(2010). 144(2019). [10] A. Andronic, P. Braun-Munzinger, J. Stachel et al., [4] J. Adam et al, (STAR Collaboration), Nat. Phys. 16, Phys. Lett. B 697, 203(2011). 409(2020). [11] J. Cleymans, S. Kabana, I. Kraus et al., Phys. Rev. C [5] S. Acharya et al, (ALICE Collaboration), Phys. Lett. B 84, 054916(2011). 797, 134905(2019). [12] S. Pal and W. Greiner, Phys. Rev. C 87, 054905(2013). [6] B.I. Abelev et al., (STAR Collaboration), Science 328, [13] S. Chatterjee and B. Mohanty, Phys. Rev. C 90, 034908 58(2010). (2014). [7] J. Adam et al., (ALICE Collaboration), Phys. Lett. B [14] A. Andronic, P. Braun-Munzinger, K. Redlich et al., Na- 7

ture 561, 321(2018). [36] J.L. Wang, D.K. Li, H.J. Li et al., Int. J. Mod. Phys. E [15] L. Xue, Y.G. Ma, J. H. Chen et al., Phys. Rev. C 85, 23(12), 1450088(2014). 064912(2012). [37] N.A. Ragab, Z.L. She and G. Chen, arXiv:1911.11919. [16] L.L. Zhu, C. M. Ko, and X.J. Yin, Phys. Rev. C 92, [38] F.X. Liu, G. Chen, Z.L. She et al., Eur. Phys. J. A 55(9), 064911(2015). 160(2019). [17] N. Shah, Y.G. Ma, J.H. Chen et al., Phys. Lett. B 754, [39] G. Chen, Y.L. Yan, D.S. Li et al., Phys. Rev. C 86, 6(2016). 054910(2012). [18] K.J. Sun and L.W. Chen, Phys. Rev. C 93, 064909 [40] G. Chen, H. Chen, J. Wu et al., Phys. Rev. C 88, (2016). 034908(2013). [19] P. Liu, J.H. Chen, Y.G. Ma, and S. Zhang, Nucl. Sci. [41] G. Chen, H. Chen, J.L. Wang et al., J. Phys. G: Nucl. Technol. 28, 55(2017). Part. Phys. 41, 115102(2014). [20] F.X. Liu, G. Chen, Z.L. She et al, Phys. Rev. C 99, [42] Z.J. Dong, Q.Y. Wang, G. Chen et al., Eur. Phys. J. A 034904(2019). 54, 144(2018). [21] T. Armstrong et al., (E864 Collaboration), Phys. Rev. [43] Z.L. She, G. Chen, H.G. Xu et al., Eur. Phys. J. A 52, C 70, 024902(2004). 93(2016). [22] Y.F. Sun and C.M. Ko, Phys. Rev. C 98, 014911(2018). [44] Z.L. She, G. Chen, D.M. Zhou et al., arXiv:1909.07070. [23] H.J. Xu, X.B. Wang, H.L. Li et al., Phys. Rev. Lett. [45] H. Nemura, Y. Suzuki, Y. Fujiwara et al., Prog. Theor. 121, 022301(2018). Phys. 103, 929(2000). [24] X.L. Zhao, G.L. Ma and Y.G. Ma, Phys. Rev. C 99, [46] B.I. Abelev et al., (STAR Collaboration), Phys. Rev. C 034903(2019). 79, 034909(2009). [25] J. Adam et al., (STAR collaboration), arXiv:1911.00596. [47] G. Agakishiev et al., (STAR Collaboration), Phys. Rev. [26] J. Zhao, F.Q. Wang. Prog. Part. Nucl. Phys. 107, Lett.108, 072301(2012). 200(2019). [48] M. M. Aggarwal et al., (STAR Collaboration), Phys. [27] G. L. Ma and B. Zhang, Phys. Lett. B 700, 39(2011). Rev. C 83, 034910(2011). [28] W.T. Deng, X.G. Huang, G.L. Ma and G. Wang, Phys. [49] J. Zhou, Light (anti-)nuclei production in the STAR ex- Rev. C 97, 044901(2018). periment at RHIC, Ph.D. thesis, Rice University, 2009 [29] L. Huang, C. W. Ma, and G. L. Ma, Phys. Rev. C 97, (unpublished). 034909(2018). [50] R. Scheibl, U. Heinz, Phys. Rev. C 59, 1585(1999). [30] L. Huang, M.W. Nie, and G.L. Ma, Phys. Rev. C 101, [51] H.H. Gutbrod, A. Sandoval, P.J. Johansen et al., Phys. 024916(2020). Rev. Lett. 37, 667(1976). [31] Z. W. Lin, C. M. Ko, B. A. Li, B. Zhang, and S. Pal, [52] B. I. Abelev et al.,(STAR collaboration), arXiv:0909 Phys. Rev. C 72, 064901(2005). .0566. [32] B.H. Sa, D.M. Zhou, Y.L. Yan et al., Comput. Phys. [53] S. Zhang, J.H. Chen, H. Crawford et al., Phys. Lett. B Commun. 183, 333(2012). 684, 224(2010). [33] Y. L. Yan, G. Chen, X. M. Li et al., Phys. Rev. C 85, [54] V. Koch, A. Majumder, J. Randrup, Phys. Rev. Lett. 024907(2012). 95, 182301(2005). [34] T. Sj¨ostrand, S. Mrenna, and P. Skands, J. High Energy [55] J. Steinheimer, K. Gudima, A. Botvina et al., Phys. Lett. Phys. 05, 026(2006). B 714, 85(2012). [35] B.L. Combridge, J. Kripfgang, and J. Ranft, Phys. Lett. B 70, 234(1977).