Pre-Equilibrium Clustering in Production of Spectator Fragments in Collisions of Relativistic Nuclei

: An algorithm of pre-equilibrium clustering of spectator matter based on the construction of the minimum spanning tree (MST) is presented. The algorithm was implemented in the Abrasion-Ablation Monte Carlo for Colliders (AAMCC) model designed to study the characteristics of spectator matter in collisions of relativistic nuclei. Due to accounting for the pre-equilibrium clusters in modelling 208 Pb– 208 Pb collisions at the LHC, the agreement of simulation results with experimental data on the average multiplicities of spectator nucleons was improved. The results of the AAMCC-MST were compared with experimental data on the interactions of 197 Au nuclei in nuclear photoemulsion. Comparison of the yields of spectator nuclei calculated for 16 O– 16 O collisions with the yields measured in interactions of 16 O with light nuclei of photoemulsion made it possible to estimate the effect of MST-clustering in small nuclear systems.


Introduction
The characteristics of the secondary particles formed as a result of collisions of relativistic nuclei, as well as the conditions of formation of quark-gluon plasma, significantly depend on the initial geometry of the collision. One of the methods to determine the geometry of a nucleus-nucleus collision event, including its centrality and reaction plane, is the detection of forward nucleons using Zero Degree Calorimeters (ZDC) [1]. Such nucleons represent spectator matter outside the overlap region of the colliding nuclei, and they are emitted at small angles with respect to the beam axis. The calculations with the Abrasion-Ablation Monte Carlo for Colliders (AAMCC) model [2] demonstrated [3,4] that spectator matter is represented not only by free nucleons but also by nuclear fragments. In addition to the numbers of free spectator neutrons, other characteristics of spectator nucleons and nuclear fragments can be considered to improve the determination of collision centrality [5]. Multiple production of intermediate mass fragments (IMF) with charges 3 ≤ Z IMF ≤ 30 known as multifragmentation was studied in several experiments. In particular, the multifragmentation of 600 MeV/nucleon 197 Au projectiles on Cu target was studied in ALADIN experiments [6]. According to the predictions of the well-known Statistical Model of Multifragmentation SMM [7], a hot thermalized nuclear system (prefragment) with an excitation energy above 3 MeV/nucleon undergoes an explosive decay into three or more fragments of comparable masses. Such less hot fragments may remain excited and their de-excitation occurs later by sequential evaporation of nucleons and light nuclei.
The SMM is based on the assumption that the prefragment is an undivided system in which thermodynamic equilibrium has been established by the time of its decay due to intensive interaction of constituent nucleons [7]. It is assumed that the average density of the decaying system is less than 1/3 of the normal nuclear density. It is obvious that the shape of a prefragment significantly affects the establishment of thermodynamic equilibrium. In particular, in central collisions of identical spherical heavy nuclei, prefragments have narrow crescent-shaped structures mostly arranged along semicircles. As a result, the establishment of the thermodynamic equilibrium in the entire system is hindered by a poor connectivity between nucleons. The connectivity of spectator matter is reduced also in ultracentral collisions (with impact parameter b < 1 fm) of elongated deformed nuclei, such as 238 U, if the long axes of the nuclear ellipsoids are perpendicular to each other. With such orientations called tip-body and side-side [8,9], spectator prefragments from nuclei with long axes perpendicular to the direction of their movement are divided into two independent clusters by the overlapping zone of the nuclei. The possibility of modeling the statistical decay of a prefragment is also questionable in the case of central collisions of light nuclei, in which the spectator matter is represented by a small number of weakly bound nucleons. In the AAMCC model described in our previous publications [2][3][4][5], each nucleusnucleus collision is considered as a two-stage abrasion-ablation process. At the abrasion stage, the removal of all participating nucleons from the primary nuclei is simulated. At the subsequent ablation stage, the decay of both prefragments representing spectator matter from the initial nuclei is simulated. The above-described effects of the loss of connectivity between nucleons in prefragments of specific shape were not taken into account in the previous version of AAMCC. In the present study, the AAMCC model is supplemented by a mechanism of pre-equilibrium clusterization of spectator matter to account for such effects. In this extension of AAMCC called AAMCC-MST, pre-equilibrium clusterization is simulated immediately after the abrasion stage on the basis of the Minimum Spanning Tree (MST) clustering algorithm.
In Section 2 of the present paper, an updated AAMCC-MST model, including the algorithm for the pre-equilibrium clustering modeling based on MST-clusterisation (MSTclustering), is presented. In Section 3, the influence of MST-clustering on various distributions of spectator fragments formed as a result of the collisions of 10.7 GeV/nucleon 197 Au nuclei with 109 Ag nuclei is discussed, and a comparison with experimental data is provided [10]. In Section 4, the impact of such clustering on the average multiplicity of spectator neutrons and protons as functions of the impact parameter in 208 Pb-208 Pb collisions at the LHC at √ s NN = 5.02 TeV is examined, and calculations are compared with preliminary data [11]. In Section 5, the influence of MST-clustering on the probabilities of the formation of various elements as spectator fragments in interactions of 200 GeV/nucleon 16 O nuclei with light nuclei is discussed, and a comparison with data [12] is made. In Section 6, the findings of the present work are summarized.

Description of AAMCC-MST Model
In this paper, the AAMCC model [2,3,5,13] is used to describe the formation of spectator matter in nucleus-nucleus collisions as a two-stage process:

1.
Removal of participating nucleons from the initial nuclei (abrasion) resulting in the formation of an excited prefragments; 2.
De-excitation of prefragments (ablation) by their decay into spectator nucleons and nuclear fragments.
At the abrasion stage, the collision impact parameter and the positions of neutrons and protons in the initial nuclei are sampled using the Glauber Monte Carlo (GMC) model v3.2 [14] according to the corresponding distributions of neutron and proton densities. Nucleons that have experienced collisions with the nucleons of the other nucleus are considered as participants, while others are treated as spectators.
The secondary ablation stage begins with the calculation of the excitation energy of the spectator prefragments from each side. Depending on the relative mass of the prefragment α = A pf /A, where A pf and A are the mass numbers of the prefragment and the initial nucleus, respectively, several methods of calculating the total E or specific = E /A pf excitation energy are available in AAMCC.
The Ericson formula is applied [15] to events with a small number a = A − A pf of nucleons removed from each of the initial nuclei. In this case, the probability distribution of E calculated for each prefragment is determined by a or, in other words, by the numbers of the respective hole states in the potential of the initial nucleus. The Ericson formula [15] is obtained by convoluting energy distributions of single-hole states. Each removed nucleon numbered as a adds excitation energy to the prefragment. On average, this value is equal to the following: where E max = 40 MeV sets the maximum energy of a single hole state in the initial nucleus [16]. In Ref. [16], it is also shown that after the removal of more than 20% of nucleons from initial 208 Pb in peripheral or semi-central collisions with various targets, the average excitation energy per nucleon obtained on the basis of the Ericson formula: significantly exceeds the typical binding energy (∼8 MeV) per prefragment nucleon. This implies a complete disintegration of the prefragment into individual nucleons. This contradicts the data of the ALADIN [6] experiment, which show that spectator matter in peripheral collisions of heavy nuclei is represented not only by free nucleons but also by nuclear fragments. This suggests that with a large number a of removed nucleons the properties of the initial nuclear potential change significantly, and the Ericson formula becomes inapplicable. Therefore, for events with significant numbers of removed nucleons a ∼ A, a phenomenological dependence of on α is used: as proposed by the ALADIN [6] collaboration. The methods (2) and (3) were discussed in detail in Ref. [16]. In the present study, a hybrid method to calculate the correlation between and α is used. The dependence (2) is switched to (3) at the value of α sw determined from the condition of continuity of the expressions (2) and (3) used at low and high , respectively. The 1/A term can be neglected for heavy nuclei, and equating the expressions (2) and (3) provides the following.
The term 1/A cannot be neglected for light nuclei; thus, the resulting equation is solved numerically by the iteration method with the values used in this study ( max = 11.5 MeV and E max = 40 MeV) for heavy projectile nuclei, including 208 Pb and α sw = 0.9287, while α sw = 0.8663 is used for 16 In previous implementations of abrasion-ablation models, particularly in Refs. [2,16,17], spectator matter formed after the ablation stage is considered as a single system with thermodynamic equilibrium established by the time of its decay. It is assumed that equilibrium is achieved due to the intense interaction between all constituent nucleons. However, as shown in Figure 1 for 208 Pb− 208 Pb collisions, such a condition of connectivity is fulfilled in peripheral collisions but violated in central ones. Indeed, in central collisions of identical spherical nuclei, spectator matter has a shape of a narrow crescent and, therefore, looses connectivity due to the low density of nucleons and larger distances between them. This is in contrast to peripheral collisions with spectator nucleons located densely in relatively large prefragments (see Figure 1). It can be expected that in central collisions a separation of spectator matter into clusters occurs before the establishment of thermodynamic equilibrium as a result of pre-equilibrium processes. Similarly, in central collisions of light projectiles with heavy nuclei, a specific ring shape of spectator matter may hinder the establishment of a thermodynamic equilibrium in the entire spectator system. An advanced algorithm based on the construction of a minimum spanning tree (MST) [18] was developed in the present study to take into account the above-described features of the shape of spectator matter resulting in pre-equilibrium clustering. The coordinates of the centers of spectator nucleons in three-dimensional space are used to determine individual nucleons and groups of nucleons that can be attributed to clusters. The MSTclustering algorithm is implemented in the following steps. All prefragment nucleons are represented as vertices of a complete weighted undirected graph for which its edge weights are equal to the moduli of the distances in 3D space between the corresponding vertex nucleons. Then, using the Kruskal algorithm [19], a minimum spanning tree with the minimum possible sum of all the edge weights is found. In the next step, the heavy edges for which their lengths are greater than some critical distance d between the nucleons are removed. This critical distance d is considered as a free parameter of the model. After that, a depth-first search algorithm is used to determine the nucleons connected by light edges that form clusters. A typical result of the MST algorithm is presented in Figure 2. The number of free nucleons and nucleon clusters found by the MST algorithm significantly depends on the parameter d. It is clear, that the application of the MST algorithm to large prefragments formed in peripheral collisions should not result in their separation into many clusters. Therefore, the value of d should not be less than the average distance of 1.6-1.8 fm between the centers of nucleons in nuclei in their ground states. This defines the values 2 < d < 4 fm tested in the present study. In modeling peripheral collisions with d 0 > 3 fm, the model underestimates the numbers of light spectator fragments and overestimates the numbers of heavy fragments. In contrast, in modeling peripheral collisions with d 0 < 2.7 fm spectator matter becomes unstable and emits many nucleons even at very low excitation energy. As a result, d 0 = 2.7 fm was found to be optimal. However, calculations with d 0 = 2.7 fm underestimate the number of free spectator neutrons in central 208 Pb-208 Pb collisions with respect to measurements [11].
It can be expected that the average multiplicity of nucleons and nucleon clusters defined by the MST algorithm for prefragments should increase with the decrease in the their average (bulk) density. The average distance between the groups of nucleons forming different clusters increases in diluted nuclear systems. However, in the AAMCC model based on the Glauber theory, nucleons of colliding nuclei are considered frozen. During a short collision time, any changes of the spatial coordinates of nucleons are neglected, and the radial expansion of the spectator system in the plane transverse to the beam axis is not taken into account. The AAMCC model also neglects the mutual Coulomb repulsion of the prefragment protons. For these reasons, direct estimation of the prefragment density at the pre-equilibrium stage is not possible in AAMCC.
However, the lower limit for the prefragment density at the pre-equilibrium stage can be estimated from the average prefragment density ρ at a later stage, after reaching the thermodynamic equilibrium. This is motivated by the expansion of the prefragment on its way to the thermal equilibrium. Following [20], one can conclude that pre-equilibrium ρ preq and eqilibrium ρ densities correlate, ρ ∼ 0.4 · ρ preq , when they are evaluated from remnant masses for events with given fragment multiplicity. In Ref. [21] break-up densities were calculated from the fractional Coulomb barrier energy, determined by fitting spectra of intermediate mass fragments produced in 3 He-197 Au, 4 He-197 Au, and 14 N-197 Au collisions. It was assumed that for the considered reactions, the heavy thermal source is spherical, and its expansion is radial. It was found [21] that the average density of a hot nuclear system with excitation energy ∼ 2 MeV is similar to the ground state nuclear density ρ 0 ≈ 0.16 fm −3 . At higher , ρ drops significantly [21], and at ≥ 5 MeV it reaches a plateau value of ρ/ρ 0 ∼ 1/3. A similar estimation is adopted by the Statistical Multifragmentation Model (SMM) [7]. There it is assumed that at the moment of an explosive decay of a hot termalized nuclear system with ≥ 3-4 MeV into fragments, its average density is also lower than ρ 0 : ρ/ρ 0 ∼ 1/4-1/3. In Ref. [22], ρ/ρ 0 was extracted by considering the level density parameters required to fit the caloric curves measured in several experiments for various fragmenting systems with 30 < A < 240. A useful compilation of results for ρ/ρ 0 obtained in [21,22] is presented in Ref. [23], where it is discussed that while they agree at 2 < < 3.5 MeV, the break-up densities derived from Coulomb barrier systematics [21] are lower at > 3.5 MeV compared to those from caloric curve data [22]. A moderate decrease in ρ/ρ 0 to ∼ 1/2 at ≥ 5 MeV was reported in Ref. [22]. We attribute the difference between [21,22] to very different techniques used in these papers and consider it as an estimate of systematic uncertainties of the methods.
The dependence of ρ( )/ρ 0 on the excitation energy at the time of pre-equilibrium clusterization is postulated in the AAMCC-MST model taking into account the abovementioned experimental results [21,22]. It is assumed that ρ( )/ρ 0 = 1 at ≤ 0 . At higher excitation energies > 0 , a phenomenological approximation of ρ( )/ρ 0 by a power function is introduced to account for the expansion: where γ = −1.02 ± 0.10 and 0 = 2.17 ± 0.23 MeV are the fit parameters that ensure the continuity of ρ( )/ρ 0 . The dependence (5) is presented in Figure 3 together with the data obtained in Refs. [21,22]. As can be seen from this figure, the results of [21,22] differ from each other at > 3.5 MeV. Our approximation (5) is closer to the data of Ref. [21] at ≥ 3-4 MeV, and they are consistent with the values ρ/ρ 0 ∼1/4-1/3 adopted in the SMM model [6,7], which is used in the present work to simulate subsequent decays of clusters resulting from the MST algorithm.  In the AAMCC-MST model, the phenomenological approximation (5) is used to determine the dependence of the clustering parameter d on the excitation energy at the preequilibrium stage. It is assumed that the excitation energy E of the spectator prefragment is divided between the resulting clusters proportionally to the number of nucleons they contain. In this manner, the of each cluster on average corresponds to the specific excitation energy of the entire prefragment. At the same time, it can be expected that for a given cluster partition, the spatial dimensions of clusters increase with an increase in their excitation energy per nucleon due to an increase in the motion of intracluster nucleons and outward pressure. As a result, the internal density of clusters decreases, approaching the density at which thermodynamic equilibrium is established. Due to a short duration of the pre-equilibrium stage, there is no significant global expansion of the cluster system. Therefore, with an increase in the excitation energy of the system, the distances between the cluster boundaries correlating with the parameter d of the MST algorithm decrease on average. The exact dependence of d on the internal density of clusters is unknown; thus, in the AAMCC-MST model, a phenomenological dependence of d ∝ ρ 1/3 ( ) is assumed. Using this dependence and Equation (5), the clustering parameter d is estimated as follows: where d 0 = 2.7 fm is the clustering parameter at the normal nuclear density ρ 0 and 0 = 2.17 MeV. The decrease in d with the increase in reflects a decrease in the connectivity of the prefragment system with the growth of . All these effects increase the number of free nucleons in central nucleus-nucleus collisions.
Statistical models of the Geant4 [24] toolkit are used in AAMCC-MST to simulate decays of clusters formed at the pre-equilibrium stage. At ≤ 4 MeV, the Weisskopf-Ewing model [25] is used to simulate sequential evaporation of nucleons. For large , the Statistical Multifragmentation Model (SMM) [7] of Geant4 [24] is enabled. A version of the Fermi Break-Up model [26] is used to simulate the decays of excited light nuclei up to 18 O.

Fragmentation of 197 Au in Nuclear Emulsion
Various distributions and correlations between charged secondary fragments were measured in experiments on fragmentation of heavy relativistic nuclei in nuclear emulsion, particularly in Ref. [10]. This technique does not allow registering neutrons and accurately determining the mass numbers of spectator fragments, but it allows determining the charges of all fragments in each fragmentation event. Several charge distributions of fragments measured by the EMU-01/12 collaboration for collisions of 10.7 GeV/nucleon 197 Au nuclei with nuclear photoemulsion (AgBrC) [10] are shown in Figure 4. The results of calculations using the AAMCC model with and without the MST algorithm are presented in the same figure. The dependences of the average maximum fragment charge Z max in each event and the average numbers of hydrogen and helium fragments, N Z=1 and N Z=2 , respectively, are presented as functions of the sum of charges Z bn of fragments with Z ≥ n in each event. The average number of intermediate mass fragments (3 ≤ Z IMF ≤ 30) in the event, denoted by M IMF as a function of Z bound = Z b2 , is of particular interest as an indicator of multifragment decays of spectator matter. The values Z b2 and Z b3 are directly connected with the size of the prefragment and, accordingly, with event centrality.
The dependence of Z max on Z b3 , measured for fragmentation of 600 MeV/nucleon 197 Au on a copper target using the ALADIN detector [6], is also presented in Figure 4. This dependence almost coincides with the one measured in Ref. [10] at significantly higher energy of 197 Au projectiles. This indicates the universality of this dependence. As can be seen from Figure 4, both variants of the AAMCC model describe well the dependence of Z max on Z b3 measured in two different experiments. Since, by definition, in each event, Z max cannot exceed Z b3 , all the points of such dependence are located below the diagonal line. The points near the diagonal at large Z b3 > 60 characterize peripheral events in which all spectator matter, except the lightest fragments, such as neutrons, protons, deuterons, tritium nuclei, 3 He, and 4 He, is represented mainly by a single fragment with Z max ≈ Z b3 . In semi-central events represented by points below the diagonal in the range 30 < Z b3 < 50, the prefragment charge is distributed among several spectator fragments of comparable charge. Both variants of calculation of the dependence of Z max on Z b3 agree with the data in general, but the agreement improves when MST-clustering is enabled.
The results obtained with MST also turn out to be closer to the measured correlation between M IMF and Z bound , which has a characteristic shape with a maximum in semicentral events, see Figure 4. In Ref. [27] such an evolution of the average multiplicity of intermediate mass fragments with an increase in the impact parameter was identified as rise and fall of multifragmentation [27] (see also [6,28]). Both AAMCC variants qualitatively reflect the dependence of M IMF on Z bound but overestimate the multiplicity in central and semi-central events at Z bound < 50. Still, the calculation with MST turns out to be closer to the data at 30 < Z bound < 50.
As observed from Figure 4, the average multiplicities of the spectator hydrogen N Z=1 and helium N Z=2 fragments essentially depend on Z b3 , and, accordingly, on centrality. The general shape of these dependencies is correctly reproduced by the model, with a satisfactory quantitative description of the multiplicity in peripheral events. At the same time, the multiplicity of helium nuclei is underestimated by AAMCC both with and without taking into account pre-equilibrium MST-clustering. The most noticeable difference between N Z=1 calculated with AAMCC-MST and the data is found for central and semicentral collisions. As expected, the inclusion of MST significantly increases the yields of protons and deuterons in central collisions. Such a discrepancy between the AAMCC-MST model and the data may stem from the restriction imposed in measurements [10] on the proton emission angle. The efficiency of a similar cut in calculations depends on the transverse momentum distribution of spectator protons obtained with the AAMCC model and requires further investigation.

Multiplicity of Spectator Nucleons in 208 Pb-208 Pb Collisions
The average multiplicities of spectator neutrons N n and protons N p were measured in 208 Pb-208 Pb collisions at √ s NN = 5.02 TeV depending on the impact parameter b [11].
These preliminary results were obtained by means of the forward neutron (ZN) and proton (ZP) Zero Degree Calorimeters (ZDC) [1] installed in the ALICE experiment at the LHC. The value of b was estimated from the multiplicity of hadrons produced in the overlap region of colliding nuclei. It was measured using other detectors of the ALICE setup by means of various methods for centrality determination [29]. According to Ref. [11], ZN calorimeters provided full acceptance (>99%) for detecting spectator neutrons. However, the reported acceptance of ZP calorimeters for protons was significantly lower. It was defined by the settings of the LHC magnetic field and in general did not exceed 70-75%. In Figure 5, N n and N p calculated by means of the AAMCC model with and without the pre-equilibrium MST-clustering are presented in comparison with the ALICE data [11]. The measured values of N p were divided by 0.7 to account for the acceptance of ZP. As can be seen from Figure 5, the average multiplicities of neutrons calculated by the AAMCC model are lower than those measured at all impact parameters. The inclusion of MST-clustering in the calculations reduces the discrepancy between the theory and experiment in central and semi-central events with b < 10 fm. In such events, the spectator matter at the pre-equilibrium stage, due to its shape, is more susceptible to separation into fragments and free nucleons, as explained in Section 2. As expected, the inclusion of MST-clustering has practically no effect on the multiplicity of nucleons calculated for peripheral events. In general, the application of MST-clustering moves the calculated N n and N p closer to the measured ones, but further tuning of the AAMCC parameters is necessary for an accurate description of the data. with and without MST-clustering (solid and dotted histograms, respectively). ALICE data [11] are represented by circles.
As known [30], in addition to deuterons, tritons, and 3 He produced directly in nucleusnucleus collisions, some additional 2 H, 3 H, and 3 He are formed at a later stage by the coalescence of free neutrons and protons. Among other publications, in Refs. [30][31][32], the formation of deuterons was modeled by the phase space coalescence of protons and neutrons. Within this approach, the probability of creating deuterons in a certain momentum space volume is proportional to the numbers of produced neutrons and protons in this volume and to the coalescence parameter B 2 . As shown [31], the d/p ratio in nucleus-nucleus collisions calculated by the coalescence model remains well below 0.01 at √ s NN > 10 GeV, in full agreement with the corresponding data. Therefore, in the case of spectator neutrons and protons, their conversion to deuterons is also expected to be unlikely, and the primary yields of free nucleons calculated with AAMCC-MST, Figure 5, will not be significantly affected.

Yields of Various Elements from Fragmentation of 16 O
Since in collisions of light nuclei spectator matter consists of a small number of nucleons, its properties should be studied separately. The results of modeling of collisions of 16 O with light nuclei (CNO) in nuclear photoemulsion are presented in Figure 6. In this figure, the multiplicity distribution of spectator α-particles and the probability of production of spectator nuclei of elements He, Li, Be, B, C, and N are presented. Calculations were performed with and without MST-clustering, and the results are compared in this figure with the data on the fragmentation of 16 O in photoemulsion [12]. According to [12], the reduced measured yields of lithium and nitrogen nuclei can be explained by the odd numbers of nucleons in these nuclei. This behavior is considered as an experimental evidence for the formation of α-clusters in primary colliding nuclei and in produced fragments. The consideration of the multiplicity distribution of α-particles presented in Figure 6 does not allow us to make a definite conclusion in favor of selecting MST-clustering in the AAMCC model, because both variants of the model significantly underestimate the probabilities of the formation of two and three spectator α-particles. At the same time, the probabilities of production of certain elements calculated by AAMCC with MST-clustering, which are presented in Figure 6, turn out to be quite close to the data [12] for He and N nuclei. A satisfactory description of the yield of B nuclei is achieved with both calculation options, but the probability of the production of Li nuclei is overestimated by the present model. At the same time, the calculated probabilities of the production of 8 Be+ 4 He+ 4 He and 12 C+ 4 He in fragmentation of 16 O are markedly underestimated by the AAMCC-MST model, as can be understood from both panels of Figure 6. A possible reason for such underestimation may be due to neglecting α-particle clustering [33,34] in initial 16 O nuclei in calculations. Such clustering may increase the probability of breakup channels containing α-particles. It is planned to take into account the cluster structure 16 O in future versions of AAMCC. This will allow predicting the composition of spectator fragments in 16 O- 16 O collisions in future runs [35] at a large hadron collider.

Conclusions
A new version of the Abrasion-Ablation Monte Carlo for Colliders model supplemented with the pre-equilibrium MST-clustering algorithm (AAMCC-MST) is presented. The inclusion of the MST-clustering makes it possible to take into account the peculiarities of the geometry of spectator matter in central collisions of medium and heavy relativistic nuclei. Due to MST-clustering, the probability of separation of spectator matter into excited clusters increases at the pre-equilibrium stage of fragmentation. As a result, the yields of free spectator nucleons in central and semi-central collisions of nuclei also increase. It was shown that the average multiplicities of neutrons and protons measured [11] by the ALICE collaboration as functions of the impact parameter in 208 Pb-208 Pb in collisions at √ s NN = 5.02 TeV are better described by the new version of AAMCC-MST taking into account MST-clustering compared to the version of the AAMCC without MST. The dependence of the clustering parameter on the excitation energy of the system made it possible to improve the description of the charge distributions of spectator fragments measured in experiments on fragmentation of 197 Au in nuclear photoemulsion.
At the same time, the inclusion of MST-clustering in AAMCC does not improve significantly the description of experimental data on fragmentation of relativistic 16 O nuclei on light nuclei in photoemulsion. It is expected that by taking into account the α-clustering structure of 16 O in next versions of AAMCC, the agreement with the experiment will be improved. While the production of 2 H, 3 H, and 3 He via the coalescence of free spectator nucleons does not play the main role, the inclusion of this process in AAMCC-MST is also considered as a possible future development of this model.  Data Availability Statement: Data sharing is not applicable to this article as no new data were created or analysed in this study.