Decay modes of highly excited nuclei

Decay modes of highly excited nuclei 1 2 first-author 1 corresp dqfang@sinap.ac.cn; dqfang@sinap.ac.cn; dqfang@sinap.ac.cn 1 corresp ygma@sinap.ac.cn ygma@sinap.ac.cn ygma@sinap.ac.cn Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China University of Chinese Academy of Sciences, Beijing 100049, China The deexcitation of single excited 112 Sn nuclei at T = 1–30 MeV is simulated using the isospin-dependent quantum molecular dynamics (IQMD) model and GEMINI model. The fragmentation mechanism, critical behavior, and kinematic characteristics are investigated within these two models. The results show that the IQMD model can be applied to the analysis of fragmentation processes, critical points, and slope temperature extraction. The results of IQMD are generally consistent with experimental 〈 M IMF 〉- Z bound data. However, GEMINI can reproduce the experimental data better than IQMD for isotopic distributions. Excited nuclei IQMD model GEMINI model Deexcitation 1 Introduction 1 Nuclear fragmentation in heavy-ion collisions is one of the most challenging subjects. Initially, the proton-induced reactions were found to produce fragments of all mass numbers [ 1 1 ]. Subsequently, the multifragmentation phenomenon was discovered. This phenomenon can be described as many particles heavier than α being produced at the same time in a single collision at energies above 30 MeV/nucleon [ 2 2 ]. Since then, many experiments and theoretical works have been dedicated to exploring the reaction mechanism. However, the phenomenon is still not fully understood [ 3 3 - 11 11 3-11 ]. In intermediate-energy heavy-ion collision processes, at first, the nuclear system is compressed and heated. Then, the system breaks up into several fragments, including light particles, during the expansion stage. Many theoretical models have been proposed to describe this reaction scenario. A brief introduction to these models can be found in Ref. [ 12 12 ]. One may treat it as a statistical process, which means that before the system disintegrates into lightly charged particles, neutrons, and IMFs (intermediate mass fragments) simultaneously, statistical equilibrium is reached in some subsystems during the expansion. This kind of equilibrium can be maintained until each subsystem becomes deexcited and the whole system reaches the freeze-out stage [ 3 3 ]. The final products are the sum of the decay results from all the subsystems. To further understand this physical process, the isospin-dependent quantum molecular dynamics model (IQMD) [ 13 13 , 14 14 ] and the GEMINI statistical model [ 15 15 ] are adopted to investigate the decay process. The combination of IQMD and GEMINI is applied in many articles to obtain isotopic distributions which are close to experimental data. In this combination, the expansion dynamics are done by IQMD and GEMINI is just used for the decay of the fragments produced by IQMD. J. Hubele et al. [ 16 16 ] also emphasized that it is necessary to have a complete dynamical treatment of the decay process. To compare the decay processes between the two models, it is interesting to study the deexcitation of a single excited nucleus within IQMD and GEMINI. In this way, the influences of the parameters and approximations employed in these models could be reduced [ 17 17 ]. By comparing them with the experimental data, it can be learned how well the IQMD and GEMINI models can describe the deexcitation processes of excited nuclei. This paper is organized as follows: a brief introduction to the IQMD and GEMINI models is given in Sect. 2. In Sect. 3, the calculated results and discussions are presented. Finally, the conclusions are given in Sect. 4. 2 Model description 1 2.1 Dynamical model: IQMD 2 The QMD model is a many-body theory that describes two-nuclei collisions from intermediate to relativistic energies. The IQMD model, which is based on QMD, has been successfully applied to heavy-ion collisions with consideration of the isospin degree-of-freedom for nucleons [ 18 18 , 19 19 ]. The nucleons are represented by Gaussian wave packets with width L ( L = 2.16 fm 2 ). When the nucleons propagate in the mean field, they are governed by the following potential: <math display="block" id="M1"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>U</mi><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><msub><mi>τ</mi><mi>z</mi></msub><mo stretchy="false">)</mo><mo>=</mo></mrow></mtd><mtd columnalign="left"><mrow><mi>α</mi><mo stretchy="false">(</mo><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo><mo>+</mo><mi>β</mi><msup><mrow><mo stretchy="false">(</mo><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo></mrow><mi>γ</mi></msup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msub><mi>τ</mi><mi>z</mi></msub><mo stretchy="false">)</mo><msub><mi>V</mi><mtext>c</mtext></msub></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><msub><mi>C</mi><mrow><mtext>sym</mtext></mrow></msub><mfrac><mrow><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>n</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>p</mtext></msub><mo stretchy="false">)</mo></mrow><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><msub><mi>τ</mi><mi>z</mi></msub><mo>+</mo><msup><mi>U</mi><mrow><mtext>Yuk</mtext></mrow></msup><mo>+</mo><msup><mi>U</mi><mrow><mtext>MDI</mtext></mrow></msup><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where α = -390.1 MeV, β = 320.3 MeV, and γ = 1.14 are the momentum-dependent nuclear equation of state (EOS) parameters, with an incompressibility of K = 200 MeV [ 14 14 ]. ρ , ρ 0 , ρ n , and ρ p are the total, normal, neutron, and proton densities, respectively. τ z (= 1 for neutron or -1 for proton) is the z th component of the isospin degree-of-freedom. V c is the Coulomb potential and C sym is the strength of the symmetry energy, which equals 32 MeV. The last two parameters, U Yuk and U MDI , represent the Yukawa potential and the momentum dependent interaction, respectively. In this work, a single excited nuclear source is simulated at different temperatures. Usually, the initial nucleus is in the state T = 0 ( T is the temperature), corresponding to the ground state. The moment distribution is calculated with the Fermi gas approximation <math id="M2"><mrow><msubsup><mi>P</mi><mtext>F</mtext><mi>i</mi></msubsup><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo stretchy="false">)</mo><mo>=</mo><mi>ℏ</mi><msup><mrow><mo stretchy="false">[</mo><mn>3</mn><msup><mi>π</mi><mn>2</mn></msup><msub><mi>ρ</mi><mi>i</mi></msub><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup></mrow></math> after its initial coordinate is sampled by the Monte Carlo method. ρ i ( i = n or p) is the local density of neutrons or protons, respectively. To study the states of nuclei with T > 0, we replace the previous moment distribution with the Fermi–Dirac distribution, <math display="block" id="M3"><mrow><mi>n</mi><mo stretchy="false">(</mo><msub><mi>e</mi><mtext>k</mtext></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mi>g</mi><mo stretchy="false">(</mo><msub><mi>e</mi><mtext>k</mtext></msub><mo stretchy="false">)</mo></mrow><mrow><msup><mi>e</mi><mrow><mfrac><mrow><msub><mi>e</mi><mtext>k</mtext></msub><mo>−</mo><mi>μ</mi></mrow><mi>T</mi></mfrac></mrow></msup><mo>+</mo><mn>1</mn></mrow></mfrac><mn>,</mn></mrow></math> where <math id="M4"><mrow><msub><mi>e</mi><mtext>k</mtext></msub><mo>=</mo><mfrac><mrow><msup><mi>p</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>m</mi></mrow></mfrac></mrow></math> is the state kinetic energy. p and m represent the momentum and mass of the nucleon, respectively. <math id="M5"><mrow><mi>g</mi><mo stretchy="false">(</mo><msub><mi>e</mi><mtext>k</mtext></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mi>V</mi><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><mn>2</mn><mi>m</mi></mrow><mrow><msup><mi>ℏ</mi><mn>2</mn></msup></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mfrac><mn>3</mn><mn>2</mn></mfrac></mrow></msup><msqrt><mrow><msub><mi>e</mi><mtext>k</mtext></msub></mrow></msqrt></mrow></math> is the state density. <math id="M6"><mrow><mi>V</mi><mo>=</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mi>π</mi><msup><mi>r</mi><mn>3</mn></msup></mrow></math> , with r = r 0 A 1/3 , is the volume of the source, and the source radius r is determined by the radius parameter r 0 and nuclear mass A . μ is the chemical potential, which is determined by the following integral equation: <math display="block" id="M7"><mrow><mfrac><mn>1</mn><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><mn>2</mn><mi>m</mi></mrow><mrow><msup><mi>ℏ</mi><mn>2</mn></msup></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mfrac><mn>3</mn><mn>2</mn></mfrac></mrow></msup><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>∞</mi></msubsup></mrow></mstyle><mfrac><mrow><msqrt><mrow><msub><mi>e</mi><mtext>k</mtext></msub></mrow></msqrt></mrow><mrow><msup><mi>e</mi><mrow><mfrac><mrow><msub><mi>e</mi><mtext>k</mtext></msub><mo>−</mo><msub><mi>μ</mi><mi>i</mi></msub></mrow><mi>T</mi></mfrac></mrow></msup><mo>+</mo><mn>1</mn></mrow></mfrac><mtext>d</mtext><msub><mi>e</mi><mtext>k</mtext></msub><mo>=</mo><msub><mi>ρ</mi><mi>i</mi></msub><mn>.</mn></mrow></math> In this paper, the simulated temperatures are from 1 MeV to 30 MeV and the excited nucleus evolves in the mean field without collisions after the initialization [ 20 20 ]. 2.2 Statistical model: GEMINI 2 GEMINI is a Monte Carlo code written in 1986 that follows the decay of a compound nucleus through a series of sequential binary decays like light-particle evaporation and symmetric fission. All possible binary decays are performed until the resulting products cannot decay any further [ 15 15 ]. There are four basic parameters as the input: fragment mass number A , charge number Z , excitation energy E * , and angular moment L . To study the same excited nuclear source, the excitation energy and angular moment of the hot nucleus are calculated in the IQMD when the initialization is just completed: <math display="block" id="M8"><mrow><msub><mover accent="true"><mi>L</mi><mo>→</mo></mover><mi>i</mi></msub><mo>=</mo><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>i</mi></msub><mo>×</mo><msub><mover accent="true"><mi>P</mi><mo>→</mo></mover><mi>i</mi></msub><mn>,</mn></mrow></math> <math display="block" id="M9"><mrow><msup><mi>E</mi><mo>∗</mo></msup><mo>=</mo><msubsup><mi>E</mi><mrow><mtext>bind</mtext></mrow><mrow><mtext>excited</mtext></mrow></msubsup><mo>−</mo><msubsup><mi>E</mi><mrow><mtext>bind</mtext></mrow><mrow><mtext>ground</mtext></mrow></msubsup><mn>.</mn></mrow></math> <math id="M10"><mrow><msub><mover accent="true"><mi>R</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> and <math id="M11"><mrow><msub><mover accent="true"><mi>P</mi><mo>→</mo></mover><mi>i</mi></msub></mrow></math> are vectors of the coordinates and momentum, respectively, and i represents the i th nucleon of the nucleus [ 21 21 ]. <math id="M12"><mrow><msubsup><mi>E</mi><mrow><mtext>bind</mtext></mrow><mrow><mtext>excited</mtext></mrow></msubsup></mrow></math> is the binding energy of the excited nucleus and it can be calculated from the IQMD model. The <math id="M13"><mrow><msubsup><mi>E</mi><mrow><mtext>bind</mtext></mrow><mrow><mtext>ground</mtext></mrow></msubsup></mrow></math> is obtained from the nuclear data table [ 22 22 ]. The latest atomic mass evaluation can be seen in Ref. [ 23 23 ]. 3 Results and discussion 1 In this section, the fragmentation mechanisms are first analyzed. Then, the critical behavior and slope temperature are compared between the two models. Finally, the comparison of the results between the calculated and experimental data is discussed. In this work, a nuclear source of normal density with 50 protons and 62 neutrons ( 112 Sn) is chosen. The simulation in IQMD is calculated up to 200 fm/ c and the number of simulated events is 1000 for each temperature point. The fragments are constructed by a coalescence model, in which nucleons with relative distances smaller than 3.5 fm and relative momenta smaller than 300 MeV/c are combined into one cluster [ 24 24 ]. 3.1 Fragmentation mechanism 2 Before comparing the fragmentation phenomena of the two models, it has to be checked whether the time-length of 200 fm/c is sufficient for the nucleus decay in IQMD. This question is investigated by looking at the mean multiplicity and production rate of the IMFs, which is defined as the fragments with 3 ≤ Z ≤ 20. From Fig. 1 Fig. 1 , it is found that the 〈 M IMF 〉 remains almost the same after 100 fm/ c at different temperatures. The production rates of the IMFs become almost zero after 100 fm/ c . It indicates that the fragment formation is finished within 100 fm/ c and the evolution time of 200 fm/ c is enough in IQMD. Figure 2 Figure 2 presents the mean multiplicity of the IMFs calculated by IQMD and GEMINI at different temperatures. They have very similar trends at lower temperatures, but the maximum value of the mean IMF multiplicity of IQMD is much bigger than that of GEMINI. The T-value corresponding to the maximum value of the mean IMF multiplicity is also higher for IQMD. In addition, nuclei in IQMD can keep the fragments’ structure, while they are almost completely disassembled into light particles in GEMINI at very high temperatures. The insert graph shows the correlation between the temperature and mean excitation energy calculated by Eq. (5). It is similar to a Fermi gas curve at low excitation energies and shows linear growth at high excitation energies. Thus, the excitation energies of the nuclei used as input for GEMINI at different temperatures are reasonable in this respect. To conveniently compare the fragmentation mechanisms between the two models at different temperatures, three fragmentation phenomena, i.e., spallation, fission, and multifragmentation, are distinguished by the multiplicity M H of the heavy fragments in a single event: M H = 1 for spallation, M H = 2 for fission, and M H = 0 with at least one fragment of Z > 2 for multifragmentation. Fragments with Z > 20 are defined as heavy fragments [ 25 25 ]. It is well known that the binary decay modes that cause the particles evaporated one-by-one from the excited nuclei are adopted in the GEMINI code. Thus, the multifragmentation that is caused by chemical instability cannot be treated with GEMINI. Here, the multifragmentation only concerns the decay results for GEMINI. In Fig. 3 Fig. 3 , the fractions of spallation, fission, and multifragmentation at different temperatures are presented for IQMD and GEMINI. Fission is a rare process that only occurs at lower temperatures. It is interesting to see the opposing trends of the spallation and multifragmentation below T = 10 MeV. There exists an apparent deviation between IQMD and GEMINI in the multifragmentation. The deviation corresponds to the discrepancy in the mean multiplicity of the IMFs at high temperatures, as we can see in Fig. 2 Fig. 2 . In general, the main fragmentation phenomena are similar in IQMD and GEMINI, which indicates that a statistical model and dynamical model may obtain comparable results for the decay process. 3.2 Critical behavior 2 It is a significant topic in nuclear physics to explore the liquid-gas phase transition of finite nuclei that is associated with the nuclear equation of state and clarification of the mechanism when a heated nucleus disassembles [ 26 26 , 27 27 ]. Many signatures have been applied to identify the critical region. The Fisher droplet model is extensively used in the study of multifragmentation and the extracted value of the power-law exponent is reasonable for the critical behavior analysis [ 28 28 , 29 29 ]. We employed a power-law fit for the charge distribution in the range Z = 2–7 for IQMD and Z =1–7 for GEMINI. The effective Fisher-law parameter τ eff is extracted from the following equation: <math display="block" id="M14"><mrow><mtext>d</mtext><mi>N</mi><mo>/</mo><mtext>d</mtext><mi>Z</mi><mo>∼</mo><msup><mi>Z</mi><mrow><mo>−</mo><msub><mi>τ</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></msup><mn>.</mn></mrow></math> The upper and bottom panels of Fig. 4 Fig. 4 show the charge distributions at T = 5, 15, and 25 MeV for GEMINI and IQMD, respectively. The charge distributions become steeper with increasing temperatures and the widths of the charge distributions of GEMINI are narrower than those of IQMD. The fitting curves are consistent with the IQMD data, but the fitting curves for GEMINI are not very good, and the errors of τ eff are quite large, as presented in Fig. 5 Fig. 5 . Thus, the charge distributions of GEMINI are not well described by the power-law dependence. The overall trend for τ eff with temperature is that it first falls to a minimum at about 8 MeV for GEMINI and 10 MeV for IQMD and then rises slowly. The minimum values of τ eff are 1.79 for GEMINI and 2.14 for IQMD, respectively. The Fisher droplet model predicted a reasonable critical exponent ~ 2.21 for the critical behavior [ 28 28 ]. The result of IQMD is closer to the critical exponent. It indicates that the critical behavior is better presented within IQMD model and a possible critical point around T = 10 MeV. The occurrence of the largest fluctuation indicates that the system undergoes a continuous transition. Campi has suggested using the event-by-event distribution of the logarithm of the largest fragment’s size versus the normalized second moment. The normalized second moment is represented by S 2 with the heaviest fragment removed, <math display="block" id="M15"><mrow><msub><mi>S</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><msub><mi>Z</mi><mi>i</mi></msub><mo>≠</mo><msub><mi>Z</mi><mrow><mtext>max</mtext></mrow></msub></mrow></msub><mrow><msubsup><mi>Z</mi><mi>i</mi><mn>2</mn></msubsup><mo>⋅</mo><msub><mi>n</mi><mi>i</mi></msub><mo stretchy="false">(</mo><msub><mi>Z</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow></mstyle></mrow><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><msub><mi>Z</mi><mi>i</mi></msub><mo>≠</mo><msub><mi>Z</mi><mrow><mtext>max</mtext></mrow></msub></mrow></msub><mrow><msub><mi>Z</mi><mi>i</mi></msub><mo>⋅</mo><msub><mi>n</mi><mi>i</mi></msub><mo stretchy="false">(</mo><msub><mi>Z</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow></mstyle></mrow></mfrac><mn>.</mn></mrow></math> Z i is the charge number of the i th fragment and n i ( Z i ) is the multiplicity of fragment Z i . Campi plots are presented in Fig. 6 Fig. 6 and they are instructive for analyzing the critical behavior [ 30 30 ]. There are three selected temperatures for each model. At low temperatures, the upper branch, which represents the liquid phase, is strongly dominant. At T =15 MeV, the lower Z max , which represents the gas phase, is strongly dominant. The transition from the liquid-dominated branch to the vapor branch then occurs, which indicates that the region of maximal fluctuations is in the range of 5 to 15 MeV. To further identify the critical point, we also introduce the quantity γ 2 , which was suggested by Campi: <math display="block" id="M16"><mrow><msub><mi>γ</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>M</mi><mn>2</mn></msub><msub><mi>M</mi><mn>0</mn></msub></mrow><mrow><msubsup><mi>M</mi><mn>1</mn><mn>2</mn></msubsup></mrow></mfrac><mn>,</mn></mrow></math> <math display="block" id="M17"><mrow><msub><mi>M</mi><mi>k</mi></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><msub><mi>Z</mi><mi>i</mi></msub><mo>≠</mo><msub><mi>Z</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munder></mstyle><msubsup><mi>Z</mi><mi>i</mi><mi>k</mi></msubsup><mo>⋅</mo><msub><mi>n</mi><mi>i</mi></msub><mo stretchy="false">(</mo><msub><mi>Z</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mn>.</mn></mrow></math> M 0 , M 1 , M 2 and M k are the zeroth, first, second and k th moments of the charge distribution, respectively [ 31 31 ]. Figure 7 Figure 7 presents the temperature dependence of the average values of γ 2 obtained from the event-by-event analysis. There are well-defined peaks in γ 2 for both models, while the maximum value of γ 2 is larger for IQMD. The maximum values are reached at T = 6 MeV and 8 MeV for GEMINI and IQMD, respectively. They are defined as the critical points where the largest fluctuations in fragment sizes are obtained. However, the peak value of γ 2 is expected to reach at least 2 for the critical behavior. Although a maximum value of γ 2 lower than 2 has been observed in a light system with A = 36 [ 27 27 ], the value extracted by GEMINI in this work is suspect for such heavy system like 112 Sn. Thus, the most probable critical region is T = 8–10 MeV in our calculation. 3.3 Kinematic characteristic 2 Nuclear temperature is a concept derived from the compound nucleus definition[ 32 32 ]. The Coulomb [ 33 33 , 34 34 ], recoil [ 35 35 ], and collective effects [ 36 36 , 37 37 ] have been observed with the slope thermometer, whose value for the temperature is extracted from the distribution of the kinetic energy spectra. Fig. 8 Fig. 8 presents the kinetic energy spectra of 2 H at four temperatures for the two models and the fitting curves are calculated by the surface-emission-type method [ 33 33 ]: <math display="block" id="M18"><mrow><mi>Y</mi><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mo>∝</mo><mfrac><mrow><mo stretchy="false">(</mo><mi>E</mi><mo>−</mo><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mrow><msubsup><mi>T</mi><mrow><mtext>slope</mtext></mrow><mn>2</mn></msubsup></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mrow><mo>−</mo><mfrac><mrow><mo stretchy="false">(</mo><mi>E</mi><mo>−</mo><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mrow><msub><mi>T</mi><mrow><mtext>slope</mtext></mrow></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> where E 0 reflects the Coulomb repulsion. As the temperature or excitation energy increases, the final kinetic energy distribution of 2 H becomes broader, but the peak positions of the spectra that are close to the Coulomb barriers seem to be insensitive to the temperature. This phenomenon can be explained by the fact that the Coulomb barrier is dependent on the charge of the initial emitting source, which is fixed in the simulation [ 38 38 , 39 39 ]. To further investigate the kinetic characteristics, the slope values of the center-of-mass kinetic energy spectra for H and He were extracted, which are shown in Fig. 9 Fig. 9 . One can observe that the slope temperatures are mass dependent for the IQMD model [ 6 6 , 40 40 , 41 41 ]. The extracted T slope values increase with increasing mass numbers for the same elements that are associated with the radial flow effect. Usually, the Coulomb effect can be observed by comparing the T slope values extracted from two fragments with the same mass number but different charge numbers. However, there is no obvious discrepancy between 3 H and 3 He. It may be that the Coulomb effect is not so apparent for there is no compression process as in the collision reaction. The T slope values of the neutron are also extracted, which are almost the same as those of the proton. Generally, the T slope values of 4 He are lower than 3 He in the collision experiments for the abundant secondary emission of 4 He from heavy fragments. The 4 He production in different evolution times were checked and there were no obvious decreases in the momentum distributions, which indicates that this effect is limited in the simulation. The results of GEMINI show similar mass dependent trends, but deviations between the fitting curves and calculated data occur with increasing temperatures, which can be seen in Fig. 8 Fig. 8 . The kinetic energy distributions of 3 He become rather dispersive at high temperatures, thus only several extracted values are given. In fact, the decay nucleus no longer acts like a thermal source for the increasing intensity of the sequential binary decay characteristic. Thus, the extracted T slope values are more reliable in the lower temperature region for GEMINI. It is interesting to see that the T slope values of 2 H are closest to the given temperatures for both models. 3.4 Comparison with experimental data 2 The projectile fragmentation of 124 Sn-on- nat Sn at 600 MeV/nucleon has been studied with the ALADIN forward spectrometer at the GSI Schwerionen Synchrotron [ 42 42 ]. To quantitatively compare it with the experimental 〈 M IMF 〉 versus Z bound projectile data, where Z bound is defined as the sum of the charges of all fragments with Z ≥ 2 in a single event, a series of primary nuclear sources are mixed in the following way to simulate the excited residual nuclei produced in the experiment. First, the average mass of the equilibrated source <math id="M19"><mrow><mover accent="true"><mi>A</mi><mo stretchy="true">¯</mo></mover></mrow></math> has the correlation with the excitation energy per nucleon E x /A suggested by Botvina et al. [ 43 43 ]: <math display="block" id="M20"><mrow><mi>α</mi><mo>≡</mo><mfrac><mrow><mover accent="true"><mi>A</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><msub><mi>A</mi><mn>0</mn></msub></mrow></mfrac><mo>=</mo><mn>1</mn><mo>−</mo><mn>0.001</mn><mo stretchy="false">(</mo><msub><mi>E</mi><mi>x</mi></msub><mo>/</mo><mi>A</mi><mo stretchy="false">)</mo><mo>−</mo><mn>0.015</mn><mo stretchy="false">(</mo><msub><mi>E</mi><mi>x</mi></msub><mo>/</mo><mi>A</mi><msup><mo stretchy="false">)</mo><mn>2</mn></msup><mn>,</mn></mrow></math> where A is assumed to have a Gaussian distribution around <math id="M21"><mrow><mover accent="true"><mi>A</mi><mo stretchy="true">¯</mo></mover></mrow></math> and A 0 is the mass number of the projectile. The width of the distribution is decided by the following formula: <math display="block" id="M22"><mrow><mi>σ</mi><mo>=</mo><msub><mi>σ</mi><mn>0</mn></msub><msub><mi>A</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><msub><mi>c</mi><mn>0</mn></msub><mfrac><mrow><msup><mi>ϵ</mi><mn>2</mn></msup></mrow><mrow><msubsup><mi>ϵ</mi><mrow><mtext>max</mtext></mrow><mn>2</mn></msubsup></mrow></mfrac><mo stretchy="false">)</mo></mrow></math> ϵ = E x /A and ϵ max represents the maximum value of ϵ . The parameters ϵ max , σ 0 , and c 0 are 8.13 MeV, 0.07, and 2, respectively, in the simulation and the details can be seen in Refs. [ 42 42 , 43 43 ]. Then, the distribution of the mass versus excitation energy is generated with the two equations above, as shown in Fig. 10 Fig. 10 . Next, 20 primary nuclear sources are chosen with initial temperatures derived from the correlation of the excitation energy and temperature. As the correlation is not precise, integer temperature values are adopted. The specific mass numbers and mixed proportions of the nuclear sources are presented in Fig. 11 Fig. 11 . The charge numbers of the different sources are obtained with the N/Z ratio remaining the same as the projectile 124 Sn. The parameters of the nuclear sources are identical for the two models. In Fig. 12 Fig. 12 , it can be seen that the simulated results of IQMD are generally consistent with the experimental data, while the results of GEMINI show apparent deviations. However, the rising and falling trends of the GEMINI data are close to the experimental data. Thus, it is reasonable continue exploring the isotopic distribution. As there is no impact parameter in the single nucleus decay process, the yields of the isotopes cannot be compared with the experimental cross sections. However, the peak position and shape of the distribution are important information. Thus, the relative yield is studied, where the isotope with the highest yield is set to 1. In Fig. 13 Fig. 13 , the shapes of the isotopic distributions by GEMINI are very close to the experimental data, while those of IQMD are much wider. To quantitatively compare the parameters of the isotopic distributions between the two models, the average mass numbers and the standard deviations σ of the mass numbers for all isotopes in the events with 10 ≤ Z bound ≤ 40 are presented in Fig. 14 Fig. 14 . The average mass numbers of the isotopes for different charge numbers Z are almost the same for IQMD and GEMINI. The σ is much larger for IQMD, which means a broader isotopic distribution. It may be that the final fragments in the IQMD model do not completely reach the ground state while the experimental fragments are in the ground state. In the GEMINI model, extensive comparisons with data from a large range of compound-nucleus masses have been used to optimize the default parameters. The good overall agreement with the experimental data for the isotopic distributions in the calculations by GEMINI can be attributed to this. 4 Conclusion 1 The deexcitation processes of single excited 112 Sn nuclei are simulated at T = 1 ~ 30 MeV with a dynamical model (IQMD) and a statistical model (GEMINI). The fragmentation mechanisms, critical behavior, and kinematic characteristics are investigated, and some results are compared with experimental data. It is found that IQMD and GEMINI give similar fragmentation phenomena in different temperature regions. However, the power-law fit and the extraction of γ 2 are more reasonable with the IQMD model. Moreover, the extracted slope temperatures of IQMD are more reasonable and closer to the original temperatures than those of GEMINI, which indicates a better kinematic description. The simulation results of IQMD with mixed thermal sources are consistent with experimental 〈 M IMF 〉 - Z bound data of the projectiles, but GEMINI can reproduce the widths of the experimental isotopic distributions better than IQMD. In general, both IQMD and GEMINI could be used to explore the decay properties of excited nuclei like the fragment distributions or fragmentation mechanisms. However, the multfragmentation phenomenon, which is related to the chemical instability, cannot be treated with GEMINI. The critical phenomena and kinematic characteristics are also associated with dynamical behaviors. Thus, they might be better described within the IQMD model. However, GEMINI shows advantages in reproducing isotopic distributions. Thus, the IQMD model can reproduce the decay properties of an excited nucleus rather well, and it is shown that the IQMD model can be used alone in most cases. References S.B. Kaufman , M.W. Weisfield , E.P. Steinberg et al ., Nuclear reactions of 197Au with 11.5- and 300-GeV protons . Phys. Rev. C 14 , 1121 ( 1976 ). doi: 10.1103/PhysRevC.14.1121 http://doi.org/10.1103/PhysRevC.14.1121 Nuclear reactions of 197Au with 11.5- and 300-GeV protons A.I. Warwick , H.H. Wiemann , H.H. Gutbrod et al ., Breakup of spectator residues in relativistic nuclear collisions . Phys. Rev. C 27 , 1083 ( 1983 ). doi: 10.1103/PhysRevC.27.1083 http://doi.org/10.1103/PhysRevC.27.1083 Breakup of spectator residues in relativistic nuclear collisions R. Nebauer , J. Aichelin , Multifragmentation - what the data tell us about the different models . Nucl. Phys. A 681 , 353 ( 2001 ). doi: 10.1016/S0375-9474(00)00539-X http://doi.org/10.1016/S0375-9474(00)00539-X Multifragmentation - what the data tell us about the different models D.Q. Fang , W.Q. Shen , J. Feng et al ., Isospin effect of fragmentation reactions induced by intermediate energy heavy ion and its disapperance . Phys. Rev. C 61 , 044610 ( 2000 ). doi: 10.1103/PhysRevC.61.044610 http://doi.org/10.1103/PhysRevC.61.044610 Isospin effect of fragmentation reactions induced by intermediate energy heavy ion and its disapperance N. Marie , A. Chbihi , J.B. Natowitz et al ., Experimental determination of fragment excitation energies in multifragmentation events . Phys. Rev. C 58 , 256 ( 1998 ). doi: 10.1103/PhysRevC.58.256 http://doi.org/10.1103/PhysRevC.58.256 Experimental determination of fragment excitation energies in multifragmentation events T.T. Ding , C.W. Ma , An improved thermometer for intermediate-mass fragments . Nucl. Sci. Tech 27 , 132 ( 2016 ). doi: 10.1007/s41365-016-0142-2 http://doi.org/10.1007/s41365-016-0142-2 An improved thermometer for intermediate-mass fragments K. Hagel , M. Gonin , R. Wada et al ., Violent collisions and multifragment final states in the 40Ca+40Ca reaction at 35 MeV/nucleon . Phys. Rev. C 50 , 2017 ( 1994 ). doi: 10.1103/PhysRevC.50.2017 http://doi.org/10.1103/PhysRevC.50.2017 Violent collisions and multifragment final states in the 40Ca+40Ca reaction at 35 MeV/nucleon M. Mocko , M.B. Tsang , L. Andronenko et al ., Projectile fragmentation of 40Ca, 48Ca, 58Ni, and 64Ni at 140 MeV/nucleon . Phys. Rev. C 74 , 054612 ( 2006 ). doi: 10.1103/PhysRevC.74.054612 http://doi.org/10.1103/PhysRevC.74.054612 Projectile fragmentation of 40Ca, 48Ca, 58Ni, and 64Ni at 140 MeV/nucleon P. Liu , J.H. Chen , Y.G. Ma et al ., Production of light nuclei and hypernuclei at High Intensity Accelerator Facility energy region . Nucl. Sci. Tech 28 , 55 ( 2017 ). doi: 10.1007/s41365-017-0207-x http://doi.org/10.1007/s41365-017-0207-x Production of light nuclei and hypernuclei at High Intensity Accelerator Facility energy region R. Ogul , N. Buyukcizmeci , A. Ergun et al ., Production of neutron-rich exotic nuclei in projectile fragmentation at Fermi energies . Nucl. Sci. Tech 28 , 18 ( 2016 ). doi: 10.1007/s41365-016-0175-6 http://doi.org/10.1007/s41365-016-0175-6 Production of neutron-rich exotic nuclei in projectile fragmentation at Fermi energies Z.Q. Feng , Nuclear dynamics and particle production near threshold energies in heavy-ion collisions . Nucl. Sci. Tech 29 , 40 ( 2018 ). doi: 10.1007/s41365-018-0379-z http://doi.org/10.1007/s41365-018-0379-z Nuclear dynamics and particle production near threshold energies in heavy-ion collisions C.W Ma and Y.G. Ma , Shannon information entropy in heavy-ion collisions . Prog. Part. Nucl. Phys . 99 , 120 ( 2018 ). doi: 10.1016/j.ppnp.2018.01.002 http://doi.org/10.1016/j.ppnp.2018.01.002 Shannon information entropy in heavy-ion collisions J. Aichelin , “Quantum” molecular dynamics—a dynamical microscopic n-body approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions . Phys. Rep . 202 , 233 ( 1991 ). doi: 10.1016/0370-1573(91)90094-3 http://doi.org/10.1016/0370-1573(91)90094-3 “Quantum” molecular dynamics—a dynamical microscopic n-body approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions C. Hartnack , R.K. Puri , J. Aichelin et al ., Modelling the many-body dynamics of heavy ion collisions: Present status and future perspective . Eur. Phys. J. A 1 , 151 ( 1998 ). doi: 10.1007/s100500050045 http://doi.org/10.1007/s100500050045 Modelling the many-body dynamics of heavy ion collisions: Present status and future perspective R.J. Charity , M.A. McMahan , G.J. Wozniak et al ., Systematics of complex fragment emission in niobium-induced reactions . Nucl. Phys. A 483 , 371 ( 1988 ). doi: 10.1016/0375-9474(88)90542-8 http://doi.org/10.1016/0375-9474(88)90542-8 Systematics of complex fragment emission in niobium-induced reactions J. Hubele , P. Kreutz , V. Lindenstruth et al ., Statistical fragmentation of Au projectiles at E/A=600 MeV . Phys. Rev. C 46 , 1577 ( 1992 ). doi: 10.1103/PhysRevC.46.R1577 http://doi.org/10.1103/PhysRevC.46.R1577 Statistical fragmentation of Au projectiles at E/A=600 MeV W. Müller , M. Begemann-Blaich , J. Aichelin , Deexcitation of single excited nuclei in the QMD model . Phys. Lett. B 298 , ( 1993 ) 27 . doi: 10.1016/0370-2693(93)91700-W http://doi.org/10.1016/0370-2693(93)91700-W Deexcitation of single excited nuclei in the QMD model X.G. Cao , X.Z. Cai , Y.G. Ma et al ., Nucleon-nucleon momentum-correlation function as a probe of the density distribution of valence neutrons in neutron-rich nuclei . Phys. Rev. C 86 , 044620 ( 2012 ). doi: 10.1103/PhysRevC.86.044620 http://doi.org/10.1103/PhysRevC.86.044620 Nucleon-nucleon momentum-correlation function as a probe of the density distribution of valence neutrons in neutron-rich nuclei X.Y. Sun , D.Q. Fang , Y.G. Ma et al ., Neutron/proton ratio of nucleon emissions as a probe of neutron skin . Phys. Lett. B 682 , 396 ( 2010 ). doi: 10.1016/j.physletb.2009.11.031 http://doi.org/10.1016/j.physletb.2009.11.031 Neutron/proton ratio of nucleon emissions as a probe of neutron skin D.Q. Fang , Y.G. Ma , C.L. Zhou , Shear viscosity of hot nuclear matter by the mean free path method . Phys. Rev. C 89 , 047601 ( 2014 ). doi: 10.1103/PhysRevC.89.047601 http://doi.org/10.1103/PhysRevC.89.047601 Shear viscosity of hot nuclear matter by the mean free path method Z.T. Dai , D.Q. Fang , Y.G. Ma et al ., Effect of neutron skin thickness on projectile fragmentation . Phys. Rev. C 91 , 034618 ( 2015 ). doi: 10.1103/PhysRevC.91.034618 http://doi.org/10.1103/PhysRevC.91.034618 Effect of neutron skin thickness on projectile fragmentation G. Audi and A.H. Wapstra , C. Thibault , The Ame2003 atomic mass evaluation: (II) . Tables, graphs and references. Nucl. Phys. A 729 , 337 ( 2003 ). doi: 10.1016/j.nuclphysa.2003.11.003 http://doi.org/10.1016/j.nuclphysa.2003.11.003 The Ame2003 atomic mass evaluation: (II) M. Wang , G. Audi , F.G. Kondev et al ., The AME2016 atomic mass evaluation (II) . Tables, graphs and references. Chinese Physics C 41 , 030003 ( 2017 ). doi: 10.1088/1674-1137/41/3/030003 http://doi.org/10.1088/1674-1137/41/3/030003 The AME2016 atomic mass evaluation (II) T.Z. Yan , Y.G. Ma , X.Z. Cai et al ., Scaling of anisotropic flows and nuclear equation of state in intermediate energy heavy ion collisions . Chin. Phys 16 , 9 ( 2007 ). doi: 10.1088/1009-1963/16/9/031 http://doi.org/10.1088/1009-1963/16/9/031 Scaling of anisotropic flows and nuclear equation of state in intermediate energy heavy ion collisions C. Lewenkopf , J. Dreute , A. Abul-Magd et al ., Fragmentation of gold projectiles with energies of 200-980 MeV/nucleon. II. Multiplicity distributions and correlations . Phys. Rev. C 44 , 1065 ( 1991 ). doi: 10.1103/PhysRevC.44.1065 http://doi.org/10.1103/PhysRevC.44.1065 Fragmentation of gold projectiles with energies of 200-980 MeV/nucleon. II. Multiplicity distributions and correlations J.B. Natowitz , R. Wada , K. Hagel , et al ., Caloric curves and critical behavior in nuclei . Phys. Rev. C 65 , 034618 ( 2002 ). doi: 10.1103/PhysRevC.65.034618 http://doi.org/10.1103/PhysRevC.65.034618 Caloric curves and critical behavior in nuclei Y.G. Ma , J.B. Natowitz , R. Wada et al ., Critical behavior in light nuclear systems: Experimental aspects . Phys. Rev. C 71 , 054606 ( 2005 ). doi: 10.1103/PhysRevC.71.054606 http://doi.org/10.1103/PhysRevC.71.054606 Critical behavior in light nuclear systems: Experimental aspects M.E. Fisher , The theory of equilibrium critical phenomena . Rep. Prog. Phys 30 , 615 ( 1969 ). doi: 10.1088/0034-4885/31/1/508 http://doi.org/10.1088/0034-4885/31/1/508 The theory of equilibrium critical phenomena J.E. Finn , S. Agarwal , A. Bujak et al ., Nuclear Fragment Mass Yields from High-Energy Proton-Nucleus Interactions . Phys. Rev. Lett 49 1321 ( 1982 ). doi: 10.1103/PhysRevLett.49.1321 http://doi.org/10.1103/PhysRevLett.49.1321 Nuclear Fragment Mass Yields from High-Energy Proton-Nucleus Interactions X. Campi , Multifragmentation: nuclei break up like percolation clusters . J. Phys. A 19 , 917 ( 1986 ). doi: 10.1088/0305-4470/19/15/010 http://doi.org/10.1088/0305-4470/19/15/010 Multifragmentation: nuclei break up like percolation clusters X. Campi , H. Krivine , Observables in nuclear fragmentation . Z. Phys. A 344 , 81 ( 1992 ). doi: 10.1007/BF01291024 http://doi.org/10.1007/BF01291024 Observables in nuclear fragmentation N. Bohr , Neutron Capture and Nuclear Constitution . Nature 137 , 351 ( 1936 ). doi: 10.1038/137344a0 http://doi.org/10.1038/137344a0 Neutron Capture and Nuclear Constitution P.M. Milazzo , G. Vannini , M. Azzano et al ., Temperature measurement of fragment emitting systems in Au+Au 35 MeV/nucleon collisions . Phys. Rev. C 58 , 953 ( 1998 ). doi: 10.1103/PhysRevC.58.953 http://doi.org/10.1103/PhysRevC.58.953 Temperature measurement of fragment emitting systems in Au+Au 35 MeV/nucleon collisions T. Odeh , R. Bassini , M. Begemann-Blaich et al ., Fragment Kinetic Energies and Modes of Fragment Formation . Phys. Pev. Lett 84 , 4557 ( 2000 ). doi: 10.1103/PhysRevLett.84.4557 http://doi.org/10.1103/PhysRevLett.84.4557 Fragment Kinetic Energies and Modes of Fragment Formation A.S. Hirsh , A. Bujak , J.E. Finn et al ., Experimental results from high energy proton-nucleus interactions, critical phenomena, and the thermal liquid drop model of fragment production . Phys. Rev. C 29 , 508 ( 1984 ). doi: 10.1103/PhysRevC.29.508 http://doi.org/10.1103/PhysRevC.29.508 Experimental results from high energy proton-nucleus interactions, critical phenomena, and the thermal liquid drop model of fragment production P.J. Siemens and J.O. Rasmussen , Evidence for a Blast Wave from Compressed Nuclear Matter . Phys. Rev. Lett 42 , 880 ( 1979 ). doi: 10.1103/PhysRevLett.42.880 http://doi.org/10.1103/PhysRevLett.42.880 Evidence for a Blast Wave from Compressed Nuclear Matter K.S. Lee , U. Heinz and E. Schnedermann , Search for collective transverse flow using particle transverse momentum spectra in relativistic heavy-ion collisions . Z. Phys. C 48 , 525 ( 1990 ). doi: 10.1007/BF01572035 http://doi.org/10.1007/BF01572035 Search for collective transverse flow using particle transverse momentum spectra in relativistic heavy-ion collisions J. Su , L. Zhu , W.J. Xie et al ., Nuclear temperatures from kinetic characteristics . Phys. Rev. C 85 , 017604 ( 2012 ). doi: 10.1103/PhysRevC.85.017604 http://doi.org/10.1103/PhysRevC.85.017604 Nuclear temperatures from kinetic characteristics C.C. Guo , J. Su , F.S. Zhang , Comparison between nuclear thermometers in central Xe+Sn collision . Nucl. Sci. Tech 24 , 050513 ( 2013 ). doi: 10.13538/j.1001-8042/nst.2013.05.013 http://doi.org/10.13538/j.1001-8042/nst.2013.05.013 Comparison between nuclear thermometers in central Xe+Sn collision X. Liu , W. Lin , M. Huang et al ., Freezeout concept and dynamical transport model in intermediate-energy heavy-ion reactions . Phys. Rev. C 92 , 014623 ( 2015 ). doi: 10.1103/PhysRevC.92.014623 http://doi.org/10.1103/PhysRevC.92.014623 Freezeout concept and dynamical transport model in intermediate-energy heavy-ion reactions X.Q. Liu , M.R. Huang , R. Wada et al ., Symmetry energy extraction from primary fragments in intermediate heavy-ion collisions . Nucl. Sci. Tech . 26 , S20508 ( 2015 ). doi: 10.13538/j.1001-8042/nst.26.S20508 http://doi.org/10.13538/j.1001-8042/nst.26.S20508 Symmetry energy extraction from primary fragments in intermediate heavy-ion collisions R. Ogul , A.S. Botvina , U. Atav et al ., Isospin-dependent multifragmentation of relativistic projectiles . Phys. Rev. C 83 , 024608 ( 2011 ). doi: 10.1103/PhysRevC.83.024608 http://doi.org/10.1103/PhysRevC.83.024608 Isospin-dependent multifragmentation of relativistic projectiles A.S. Botvina , I.N. Mishustin , M. Begemann-Blaich et al ., Multifragmentation of spectators in relativistic heavy-ion reactions . Nuclear Physics A 584 , 737 ( 1995 ). doi: 10.1016/0375-9474(94)00621-S http://doi.org/10.1016/0375-9474(94)00621-S Multifragmentation of spectators in relativistic heavy-ion reactions 10.1007/s41365-018-0427-8 This work was supported by the National Natural Science Foundation of China (Nos. 11421505, 11475244, and 11175231) and the National Key Research and Development Program of China (High Precision Nuclear Physics Experiments). 2018-02-11 2018-03-16 2018-03-19 2018-04-25 2018-06-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0427-8 http://dx.doi.org/10.1007/s41365-018-0427-8 . 2018 Nuclear Science and Techniques 06 29