Solving the puzzle of high light (anti)-nuclei production in ultra-relativistic heavy ion collisions

Tim Neidig,∗ Kai Gallmeister,† and Carsten Greiner Institut f¨urTheoretische Physik, Johann Wolfgang Goethe-Universit¨at, Max-von-Laue-Strasse 1, 60438 Frankfurt am Main, Germany

Marcus Bleicher Institut f¨urTheoretische Physik, Johann Wolfgang Goethe-Universit¨at, Max-von-Laue-Strasse 1, 60438 Frankfurt am Main, Germany and Helmholtz Research Academy Hesse for FAIR (HFHF), GSI Helmholtz Center, Campus Frankfurt, Max-von-Laue-Straße 12, 60438 Frankfurt am Main, Germany

Volodymyr Vovchenko Nuclear Science Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA (Dated: September 2, 2021) The creation of loosely bound objects in heavy ion collisions, e.g. light clusters, near the phase transition temperature (Tch ≈ 155 MeV) has been a puzzling observation that seems to be at odds with nucleosynthesis suggesting that deuterons and other clusters are formed only below a temperature T ≈ 0.1 − 1 MeV. We solve this puzzle by showing that the light cluster abundancies in heavy ion reactions stay approximately constant from chemical freeze-out to kinetic freeze-out. To this aim we develop an extensive network of coupled reaction rate equations including stable hadrons and hadronic resonances to describe the temporal evolution of the abundancies of light (anti-)(hyper-)nuclei in the late hadronic environment of an ultrarelativistic heavy ion collision. It is demonstrated that the chemical equilibration of the light nuclei occurs on a very short timescale as a consequence of the strong production and dissociation processes. However, because of the partial chemical equilibrium of the stable hadrons, including the nucleon feeding from ∆ resonances, the abundancies of the light nuclei stay nearly constant during the evolution and cooling of the hadronic phase. This solves the longstanding contradiction between the thermal fits and the late stage coalescence (and the Big Bang nucleosynthesis) and explains why the observed light cluster yields are compatible with both a high chemical production temperature and a late state emission as modelled by coalescence. We also note in passing that the abundancies of the light clusters in the present approach are in excellent agreement with those measured by ALICE at LHC.

The production of light nuclei in an expanding dissociation of such light nuclei [11–14]. The two latter hadronic medium has been the driving force of the Big descriptions have certain shortcomings, for instance, the Bang nucleosynthesis that has led to the creation of the energy conservation is not obeyed in the coalescence ap- light atomic nuclei in the universe. Recently, similar con- proach whereas the former treatment of nuclei as point- ditions were recreated in ultra-relativistic heavy ion col- like particles in the kinetic approach may be question- lisions to allow for a deeper understanding of light clus- able, but recent works are going beyond this treatment ter production under controlled collisions. Especially, e.g. [13]. At first glance, both approaches seem incom- the yields of light nuclei like deuterons, tritons, helium-3 patible with each other, however, as we will show below and helium-4, their anti-particles, and also hyper-tritons, this is not the case (see also [6] for an alternative discus- have been measured by the ALICE collaboration at LHC sion of this topic). arXiv:2108.13151v2 [hep-ph] 1 Sep 2021 [1–3]. It was found that their abundances are in a re- Our explanation is based on the assumption of detailed markable agreement with the predictions of the statisti- balance – the law of mass action. Recently, this princi- cal hadronization model, characterized by the chemical ple was used to show how the final yields of light nuclei freeze-out temperature of T = 155 MeV and nearly van- ch are determined by the initial (chemical) freeze-out mul- ishing net baryon density [4,5]. This was a surprising re- tiplicities of the stable hadrons, akin with the nuclear sult, because one would usually expect (in line with Big equivalent of the Saha ionization equation of cosmology Bang nucleosynthesis) that the production of such loosely [15]. In [16], this idea of detailed balance resp. a Saha bound states with binding energies of O(2 MeV) would equation has also been studied in a picture, where the only happen at much lower on the order of usual hadron prescription has been extended by so the binding energy. The idea of late, i.e. low temperature called Hagedorn resonances [16]. cluster formation is usually realized via the coalescence of nucleons in the final state [6–10] or through kinetic trans- In a heavy-ion setup one uses the framework of partial port approach describing the continuous production and chemical equilibrium (PCE) [17] to describe the hadronic environment after the chemical freeze-out. The total 2 abundances of stable hadrons are fixed by introducing thermal and chemical equilibrium. The contributions of non-equilibrium chemical potentials [15, 18], while the resonance R decays and regenerations to the change rate yields of resonances change [19, 20], being determined of the multiplicities of the decay products, dNX /dt and by the relative equilibrium of decays and regenerations. dNY /dt are given by the rhs. of Eq. (3) with the opposite The light nuclei are assumed to be in relative chemical sign. equilibrium in this thermal bath. Considering the nuclear reactions, catalized by some In this letter we relax the assumption of instanta- particle (e.g. pionic [11]), similar expressions can be writ- neous equilibration of nuclear reactions, as well as res- ten down, based on the knowledge of the basic 2 → n onance decays and regenerations, by means of a set of process. coupled (reaction) rate equations describing the chemi- The decay rates used are thermal averaged values over inelastic cal composition of stable hadrons (pions, kaons, (anti- experimentally measured cross sections (e.g. σπ+d ) or )nucleons, Λs), hadronic resonances and light (strange) resonance decay widths taken from Particle Data Ta- (anti-)nuclei. Such a description of the chemical evo- bles [28], lution of the degrees of freedom and especially also of rare probes is well known since a long time in relativistic σA+X→aN+X vrel α = NX heavy ion collisions (e.g. (multi-) strange baryons [21– A+X→aN+X V (4)

24]) and in cosmology (e.g. nucleosynthesis, see [25] and =α ˜A+X→aN+X NX . references therein). The particles of special interest in the rate equations Here we introduced the numerical factorα ˜ to stress the are the light nuclei d(= H2), t(= H3), He3, He4 and the explicit dependence of α of the catalyzing mesons, whose 3 multiplicities are also time dependent, as they are part strange light nucleus HΛ (and their antiparticles). In the mesonic sector, π, K and K are considered as stable of the network. particles, while additionally the resonances ρ, ω and K∗ For the light nuclei we approximate the cross sections and K∗ are included. The baryonic sector is covered by summing over their constituents, e.g. σt+π→3N+π ≈ 2σ −σ . Since experimental data for the by nucleons N(= n, p), the ∆(1232) resonance and the d+π→2N+π p+π→p+π 3 strange Λ baryon as the dominating species. Also here all cross section for π/K +HΛ are not available, the values of antiparticles are taken into account. For simplicity, we 90 mb resp. 30 mb are used, based on the reaction π/K + 3 omit higher lying ∆ and N ∗ resonances [26]. In order to H . compensate this omission, an effective degeneracy for the In order to calculate the multiplicities of particles cre- ∆(1232) resonances is introduced, which has been chosen ated in heavy ion collisions, we incorporate the expan- eff sion and the cooling of the fireball in the hadronic phase, to be twice the fundamental degeneracy, g∆ = 2g∆, in order to fit the initial multiplicities in accordance to the through a time dependence of the volume, V (t), and the statistical hadronisation model [5]. temperature, T (t). For a given (expanding) volume V In the rate equation approach two different types of the temperature T is obtained by in the framework of reactions are considered: first the catalized break up (and the PCE, so that the abundance of each stable hadron fusion) reactions for the nuclei species (including the resonance contribution) are con- served, and by conserving the total entropy of the system A + X a · N + X (1) as well as the net baryon number and net strangeness. These relations provide a set of non-linear equations, and secondly the decay and formation of unstable reso- which can be solved for any T ≤ Tch to obtain V (T ) nances and µi(T )(i = S, B, N, π etc.). The initial conditions at Tch namely V (Tch) and the baryon and strangeness R X + Y. (2) chemical potential are obtained from the experimental As all the reactions are of the type (2 ↔ n), the collision data [1–3,5, 29, 30]. Finally, the time dependence of the rates can be expressed in terms of binary cross sections. volume is given by the parametrization that incorporates In the case of resonance decays, the change in the mul- longitudinal and transverse expansion [31], tiplicity of resonance R with time t [27] may be expressed 2 2 t t⊥ + t as e.g. for two decay particles V (t) = Vch 2 2 (5) tch t⊥ + tch dNR XY = −αR X+Y (NR − cR NX NY ) (3) with t⊥ = 6.5 fm/c and tch = 9 fm/c. In the present work, dt all calculations assume an initial temperature Tch = where α can be interpreted as the decay rate of 155 MeV for the chemical freeze-out. R X+Y R → X + Y , but also as a scaled cross section for X + In total, the network contains 23 rate equations, given eq eq eq XY . Y → R. The factor cR = NR /(NX NY ), given by the by the number of stable and unstable particles and an- eq ratio of the multiplicities in equilibrium, Ni , is dictated tiparticles. In the following, we will only list the most eq by the detailed balance. Any Ni is calculated at full important ones. 3

As a first example, the rate equation for the deuteron is given as 10 1

dNd X N 2 2 = − α˜d+x 2N+x Nx(Nd − cd NN ) (6) 10 3 dt x=π,K,K 10 5 N eqN eq··· XY ... . A B with the generalized definition cAB... = N eqN eq··· and multiplicities X Y 10 7 d XX X2 He3 using shortcuts like c(··· ) = c(··· ), while the rate equation He4 H3 for the nucleons is given by 10 9 0.07 0.08 0.09 0.10 0.11 0.12 0.13 0.14 0.15 T [GeV]

dNN Nπ X FIG. 1. The multiplicity of the light nuclei as function =α ˜∆ N+π (N∆ − c∆ NN Nπ) + AX dt of the (decreasing) temperature T . Solid lines represent the X=d,t,He3,He4 results of the rate equations, while dashed curves show the X AX N AX result of the Saha equation. The colored bands represent the × α˜X+x A +x Nx(NX − cX NN ) X experimental data (ALICE) [1–3] with their uncertainty. x=π,K,K X N 2Λ 2 + 2˜α N (N 3 − c 3 N N ). H3 +x NNΛ+x x H H N Λ Λ Λ Λ x=π,K,K (7) in the whole temperature range, supporting our previous findings.

Here AX is the nucleon content of a nucleus X. The Next we investigate the equilibration time of light rate equations for pions, being a catalyzing particle in nuclei starting from an equilibrated hadron gas with- many equations, is itself quite simple, since its affected out clusters. We also demonstrate that cluster (re- only by the resonances ∆, ρ and ω, )equilibration is robust with respect to the thermal value when the recombination reactions start not immediately dNπ Nπ after the chemical freeze-out but at lower temperatures =α ˜∆ Nπ (N∆ − c∆ NN Nπ) dt (nucleosynthesis starts only at T = 90, 100, 120 and Nπ 155 MeV). For this purpose we look at the behavior of +α ˜ (N∆ − c NN Nπ) ∆ Nπ ∆ (8) the deuteron to proton ratio obtained in these different π2 2 + 2˜α (Nρ − c N ) ρ 2π ρ π scenarios as shown in Fig.2 as function of time. Qualita- π3 3 tively the same results were found for the other light nu- + 3˜αω 3π (Nω − cω Nπ) . clei, which are indeed approaching the equilibrium value The chemical processes described via Eqs. (7) and (8), slightly faster. and also the others conserve the total number of stable hadrons, either direct or carried in a resonance or light 3.0 Equilibrium nuclei. When the rates are indeed large compared to 3*Equilibrium the expansion rate, then the abundancies of the stable 2.5 Nucleosynthesis after 90 MeV d p Nucleosynthesis after 100 MeV N N

Nucleosynthesis after 120 MeV o hadrons fulfill PCE. i 2.0

t Nucleosynthesis at 155 MeV a r In Fig.1 we show the results for the yields of the light d

e 1.5 s i nuclei using the rate equations. These results are com- l a

m 1.0

pared to the ones obtained in the Saha equation limit [15], r o as well as to the experimental data of the ALICE Col- n 0.5 laboration [1–3]. The calculation starts at the chemi- 0.0 cal freeze-out, Tch = 155 MeV, and covers temperatures 10 12 14 16 18 fm down to T = 70 MeV. The yields of light nuclei ex- t [ c ] hibit only a very small temperature dependence [18], FIG. 2. Ratio of deuterons to protons normalised to the with a minor exception of H3 , which is attributed to Λ same ratio when starting from equilibrium at Tch = 155 MeV the missing hyperon resonances [15]. Generally, a very for different initial conditions. good agreement with the experimental data is obtained for the whole temperature range. This indicates that one can not distinguish the emission of light clusters from the It is notable that one reaches the (T = 155 MeV)- chemical freeze-out hypersurface from an emission dur- equilibrium ratio even if one starts with the potential ing the later stages of the reaction by looking at the mul- regeneration reactions only at 120 MeV. For lower tem- tiplicities alone. Let us remark in passing that the rate peratures like 100 MeV or 90 MeV only 88 % or 70 % of equations results stay within 5-10 % of the Saha equation the equilibrium value is reached. The typical kinetic 4 freeze-out is expected to take place between 100 MeV <

Tkin < 120 MeV [2, 19] in central collisions, so that one is 2.5 d inverse rates fit tempted to take the conclusion that the freeze-out condi- He3 inverse rates ] c

m fit f 2.0 [ He4 inverse rates tions for the light nuclei are close to the (T = 155 MeV)-

e fit m equilibrium and dictated by the principle of PCE. i t

n 1.5 o

We also want to stress the connection between the an- i t a r nihilation rates α for each reaction and the equilibration b i l

i 1.0 u time of a particle species. For fixed volume and a given q temperature the hadrons and resonances are in full ther- e 0.5 mal and chemical equilibrium and the abundancies of the 0.08 0.09 0.10 0.11 0.12 0.13 0.14 0.15 light nuclei are set to zero. As a particular example we T [GeV] show the normalised ratio of deuterons as a function of time. The results are shown in Fig.3. By performing an FIG. 4. Equilibration time of different light nuclei for dif- ferent temperatures in a box (V = 4000 fm3). The solid line represents the inverse of the sum of the rates and the dots are

1.2 obtained by exponential fits to curves from Fig.3 and alike.

1.0 u q d e d N N

o 0.8 i t

a equations approach as follows: r

d 0.6 e s i l dNN 5 a +=α ˜ (−N N + cπ N 5) m 0.4 N N π r N+N 5π NN

o dt

n Equilibrium 0.2 105 MeV dN 5 130 MeV N π 5 +=α ˜ (−NN NN + c Nπ) (10) 155 MeV dt N+N 5π NN 0.0 10 12 14 16 18 20 fm dNπ 5 t [ c ] π 5 += 5˜α (NN NN − c Nπ) . dt N+N 5π NN FIG. 3. Time evolution of the number of deuterons normal- In principle, the cross section for these reactions can ized to the equilibrium value with V = 4000 fm3 and start- ing from zero initial yield for three different temperatures: be taken from the data on inelastic p+p scattering. Here, however, we vary its thermal averaged value in a broad T = Tch = 155 MeV (orange), T = 130 MeV (green) and T = 105 MeV (red) range spanning 0-100 mb. Typically 50 mb [18] is the ap- propriate value for the inelastic/annihilation cross sec- tion. The corresponding effect of the deuteron yields at exponential fit to these different temperature branches, decreasing temperatures T < Tch is depicted in Fig.5. one can extract the equilibration time of the deuterons. The results indicate that baryon annihilation leads to For Fig.4 this has also be done for He 3 and He4. By a suppression of the deuteron number at intermediate comparing these results to the inverse of the sum of the stages of the expansion, but these modifications are small rates, those are in perfect agreement with the formula (5 %) at the kinetic freeze-out at T = 100 MeV. The in- Eq. (9) for a light nucleus A, fluence of the baryon annihilation is somewhat stronger for the heavier nuclei. Here, the modifications are up 1 X 25%, but still in the range of the experimental uncer- = α˜ Nx = αA. (9) eq A+x aN+x τA tainty. Hence, our main conclusions obtained without x=π,K,K incorporating NN annihilation are unaffected. It follows from Fig.4 that the equilibration time of the In summary, we have shown that the light nuclei abun- light nuclei are very short and in the range 0.5-1.5 fm/c. dances equilibrate quickly towards their partial chemi- So far with our network, the amount of stable hadrons cal equilibrium yields in an expanding hadronic environ- is strictly conserved. Finally, we now would like to study ment. Because of the PCE, the abundance of the light the potential role of baryon-antibaryon annihilations on nuclei show only a minor temperature dependence when the abundances of light nuclei. The possible relevance decreasing the decoupling temperature. This explains of baryon annihilation has been raised since a long time why the late stage coalescence results in similar light nu- [32, 33], and has recently gained considerable attention at clei yields as a thermal analysis with Tch = 155 MeV, SPS, RHIC, and LHC energies [34–39]. These processes albeit coalescence happens at a much lower temperature. can be easily implemented in our network. Let us for To this aim, we employed a set of 23 coupled rate instance take the process N + N 5π [18]. We take 5π equations among stable hadrons, resonances, and light as a mean value of annihilating pions, it could in principle (strange)(anti-)nuclei. The corresponding chemical equi- vary in a range 3-7π [23]. This reaction is not part of the libration times are are found to be very short, i.e. between original PCE but can be incorporated within the rate 0.5 fm/c close to the chemical freeze-out and 1-1.5 fm/c 5

ACKNOWLEDGEMENTS

1.15 Rate Equation without N + N 5

N + N 5 = 10 mb

N + N 5 = 30 mb = 50 mb V.V. thanks V. Koch and D. Oliinychenko for fruit- 1.10 N + N 5 = 70 mb N + N 5 ful discussions. K.G. was supported by the Bundesmin- N + N 5 = 100 mb ) ) c T

T 1.05

( isterium f¨ur Bildung und Forschung (BMBF), grant ( d d N N No. 3313040033. We jointly acknowledge support by the 1.00 Deutsche Forschungsgemeinschaft (DFG) through the CRC-TR 211 ’Strong-interaction matter under extreme 0.95 conditions’. V.V. acknowledges the support through the 0.08 0.10 0.12 0.14 T [GeV] Feodor Lynen Program of the Alexander von Humboldt foundation, the U.S. Department of Energy, Office of Sci- FIG. 5. Deuteron yield at various temperatures normalised ence, Office of Nuclear Physics, under contract number to the equilibrium value at Tch = 155 MeV in a system de- DE-AC02-05CH11231231, and within the framework of scribed by the network of rate equations with (strong) baryon- the Beam Energy Scan Theory (BEST) Topical Collab- antibaryon annihilation reactions N + N for various values of oration. the corresponding cross section.

∗ at lower temperatures around T ∼ 100 MeV of the (po- Corresponding author; [email protected] kin † Present address: Institut f¨ur Theoretische Physik, tential sequential) kinetic freeze-out of the hadronic fire- Justus-Liebig-Universit¨at,Heinrich-Buff-Ring 16, 35392 ball. Consequently, the abundances stay very close to Gießen, Germany their ”initial” values given by the statistical model at [1] J. Adam et al. (ALICE), Phys. Lett. B 754, 360 (2016), the chemical freeze-out and may be accurately described arXiv:1506.08453 [nucl-ex]. by a Saha equation [15]. [2] J. Adam et al. (ALICE), Phys. Rev. C 93, 024917 (2016), arXiv:1506.08951 [nucl-ex]. This implies that the light nuclei (e.g. deuterons) do [3] S. Acharya et al. (ALICE), Phys. Rev. C 101, 044906 not necessarily have to be formed at the chemical freeze- (2020), arXiv:1910.14401 [nucl-ex]. out: Even if the nucleosynthesis only starts at tempera- [4] A. Andronic, P. Braun-Munzinger, J. Stachel, tures as low as 120 MeV, the statistical model values are and H. Stocker, Phys. Lett. B 697, 203 (2011), essentially reached in the final state. If the nucleosynthe- arXiv:1010.2995 [nucl-th]. sis would start only at 90 MeV one still obtains 70 % of [5] A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, Nature 561, 321 (2018), arXiv:1710.09425 the initial equilibrium yield. [nucl-th]. We have also considered the effect of annihilations be- [6] S. Mrowczynski, Acta Phys. Polon. B 48, 707 (2017), tween nucleons and anti-nucleons and found that it only arXiv:1607.02267 [nucl-th]. leads to a small modification in the total yields of the [7] S. Sombun, K. Tomuang, A. Limphirat, P. Hillmann, C. Herold, J. Steinheimer, Y. Yan, and M. Bleicher, light nuclei. As such, our study emphasizes the validity Phys. Rev. C 99, 014901 (2019), arXiv:1805.11509 [nucl- of the (partial) chemical equilibrium scenario for the ex- th]. panding hadronic fireball created in heavy-ion collisions. [8] W. Zhao, L. Zhu, H. Zheng, C. M. Ko, and H. Song, One may be even come to the conclusion that the re- Phys. Rev. C 98, 054905 (2018), arXiv:1807.02813 [nucl- th]. markable agreement of the experimental yields of light [9] F. Bellini, K. Blum, A. P. Kalweit, and M. Puccio, Phys. nuclei with the yields obtained within the statistical Rev. C 103, 014907 (2021), arXiv:2007.01750 [nucl-th]. hadronization description [5] at the temperature of Tch ≈ [10] S. Gl¨aßel, V. Kireyeu, V. Voronyuk, J. Aichelin, 155 MeV leads strong support to the principle of the par- C. Blume, E. Bratkovskaya, G. Coci, V. Kolesnikov, and tial chemical equilibrium, and hence the validity of the M. Winn, (2021), arXiv:2106.14839 [nucl-th]. Saha equation. [11] D. Oliinychenko, L.-G. Pang, H. Elfner, and V. Koch, Phys. Rev. C 99, 044907 (2019), arXiv:1809.03071 [hep- All results given are based on the law of mass action. ph]. If there are (strong) annihilation reactions, (strong) gen- [12] D. Oliinychenko, C. Shen, and V. Koch, Phys. Rev. C eration processes must also take place according to prin- 103, 034913 (2021), arXiv:2009.01915 [hep-ph]. ciple of detailed balance. It does, however, not tell how [13] K.-J. Sun, R. Wang, C. M. Ko, Y.-G. Ma, and C. Shen, the (re-)generation processes are realized microscopically, (2021), arXiv:2106.12742 [nucl-th]. [14] J. Staudenmaier, D. Oliinychenko, J. M. Torres-Rincon, but they have to occur in time. Of course, the under- and H. Elfner, (2021), arXiv:2106.14287 [hep-ph]. standing of the microscopic process of the generating pro- [15] V. Vovchenko, K. Gallmeister, J. Schaffner-Bielich, cesses is a formidable quantum mechanical task, which is and C. Greiner, Phys. Lett. B 800, 135131 (2020), outside the scope of the present study. arXiv:1903.10024 [hep-ph]. 6

[16] K. Gallmeister and C. Greiner, Eur. Phys. J. A 57, 62 [28] P. A. Zyla et al. (Particle Data Group), PTEP 2020, (2021), arXiv:2007.08258 [hep-ph]. 083C01 (2020). [17] H. Bebie, P. Gerber, J. L. Goity, and H. Leutwyler, Nucl. [29] B. Abelev et al. (ALICE), Phys. Rev. C 88, 044910 Phys. B 378, 95 (1992). (2013), arXiv:1303.0737 [hep-ex]. [18] X. Xu and R. Rapp, Eur. Phys. J. A 55, 68 (2019), [30] B. B. Abelev et al. (ALICE), Phys. Rev. Lett. 111, arXiv:1809.04024 [nucl-th]. 222301 (2013), arXiv:1307.5530 [nucl-ex]. [19] A. Motornenko, V. Vovchenko, C. Greiner, and [31] Y. Pan and S. Pratt, Phys. Rev. C 89, 044911 (2014). H. Stoecker, Phys. Rev. C 102, 024909 (2020), [32] R. Rapp and E. V. Shuryak, Phys. Rev. Lett. 86, 2980 arXiv:1908.11730 [hep-ph]. (2001), arXiv:hep-ph/0008326. [20] B. Tomasik, P. Hillmann, and M. Bleicher, (2021), [33] C. Greiner and S. Leupold, J. Phys. G 27, L95 (2001), arXiv:2107.03830 [nucl-th]. arXiv:nucl-th/0009036. [21] P. Koch, B. Muller, and J. Rafelski, Phys. Rept. 142, [34] J. Steinheimer, J. Aichelin, and M. Bleicher, Phys. Rev. 167 (1986). Lett. 110, 042501 (2013), arXiv:1203.5302 [nucl-th]. [22] H. W. Barz, B. L. Friman, J. Knoll, and H. Schulz, [35] F. Becattini, M. Bleicher, T. Kollegger, T. Schuster, Nucl. Phys. A 484, 661 (1988), [Erratum: Nucl.Phys.A J. Steinheimer, and R. Stock, Phys. Rev. Lett. 111, 492, 663–663 (1989)]. 082302 (2013), arXiv:1212.2431 [nucl-th]. [23] J. Noronha-Hostler, C. Greiner, and I. A. Shovkovy, [36] E. Seifert and W. Cassing, Phys. Rev. C 97, 024913 Phys. Rev. Lett. 100, 252301 (2008), arXiv:0711.0930 (2018), arXiv:1710.00665 [hep-ph]. [nucl-th]. [37] E. Seifert and W. Cassing, Phys. Rev. C 97, 044907 [24] J. Noronha-Hostler, M. Beitel, C. Greiner, and (2018), arXiv:1801.07557 [hep-ph]. I. Shovkovy, Phys. Rev. C 81, 054909 (2010), [38] O. Savchuk, V. Vovchenko, V. Koch, J. Steinheimer, and arXiv:0909.2908 [nucl-th]. H. Stoecker, (2021), arXiv:2106.08239 [hep-ph]. [25] E. W. Kolb and M. S. Turner, The Early Universe, [39] O. Garcia-Montero, J. Staudenmaier, A. Sch¨afer,J. M. Vol. 69 (1990). Torres-Rincon, and H. Elfner, (2021), arXiv:2107.08812 [26] Essentially the same evolution of particle multiplicities is [hep-ph]. obtained when using the full particle list, as implemented [40] V. Vovchenko and H. Stoecker, Comput. Phys. Commun. within the Thermal-FIST package [40] and reported in 244, 295 (2019), arXiv:1901.05249 [nucl-th]. Ref. [41]. [41] V. Vovchenko, “Nuclear clusters in an off-equilibrium [27] The time can be interpreted as the comoving Bjorken thermal model,” (2020), Talk at Mini-Workshop “Origin time. of nuclear clusters in hadronic collisions”.