Chen-Chen Guo· Jun Su· Long Zhu

Abstract The isospin fractionations in 124Sn, 107Sn ＋120Sn at 600 MeV/nucleon, and 136Xe, 124Xe ＋ 208Pb at 1000 MeV/nucleon are investigated by the isospin-dependent quantum molecular dynamics model coupled with the statistical code GEMINI. The yield ratio as a function of the binding energy difference for light mirror nuclei 3H/3He, 7Li/7Be, 11B/11C, and 15N/15O is applied to estimate the ratio between neutrons and protons in the gas of the fragmenting system. By comparing the estimated values resulting from the simulations with and without the GEMINI code, it was found that the secondary decay distorts the signal of the isospin fractionation. To minimize the secondary decay effects, the yield ratio of the light mirror nuclei 3H/3He as well as its double ratio between two systems with different isospin asymmetries of the projectiles is recommended as robust isospin observables.

Keywords Quantum molecular dynamics model (QMD) ·Isospin fractionation · Secondary decay effect

1 Introduction

Investigations on fractionation have stimulated significant theoretical and experimental efforts in many branches of science and technology [1,2]. In recent years, there has been increasing interest in nuclear isospin fractionation,which indicates the separation of neutron-rich gas during a liquid-gas phase transition in nuclear matter [3–6]. It was reported in the early 1980s although the concept of isospin fractionation was not applied. For example, it was found that fragment emission amplified the ratio of free neutrons to protons relative to that in the total system [7]. In 2000,the concept of isospin fractionation was used to describe the measurement of the isotopic distributions of light particles and intermediate mass fragments in the collisions between tin isotopes at 50 MeV/nucleon [8]. Isospin fractionation is an important phenomenon in heavy ion collisions in the energy region from the Fermi energy to gigaelectron volt. It has been proposed that the observable associated with the isospin fractionation is sensitive to the isospin dependence of the nuclear equation-of-state in isospin-asymmetric nuclear matter [9–11]. The differential isospin fractionation is also considered to be significantly affected by the isospin splitting of the nucleon effective mass [12, 13].

Several dynamical effects on the isospin fractionation have been studied. For example, the effect of Coulomb interaction on the isospin fractionation has received attention [14–16]. The relationship between fragment chemical composition and the expected isospin fractionation has also been explored [17]. However, studies on the secondary decay effects remain few. Isospin fractionation occurs during the fragmenting process in heavy ion collisions that form excited fragments [18–22]. The decay of these excited fragments changes the isospin asymmetry of the gas (free nucleons) and liquid (fragments) phases[23, 24]. Thus, the secondary decay effects should be considered when comparing model calculations with experimental data. In our previous work, the dynamic properties, secondary decay effects, and isospin dependence of projectile fragmentations in124Sn,107Sn ＋120Sn collisions at 600 MeV/nucleon were studied by an isospindependent quantum molecular dynamics (IQMD) model coupled with the statistical code GEMINI [25–27]. In this work, the isospin fractionation and its secondary decay effects are studied using a momentum-dependent potential version of the IQMD model [28, 29] coupled with the GEMINI code.The method is introduced briefly in Sect. 2.The results and discussions are presented in Sect. 3.Finally, a summary is given in Sect. 4.

2 Theoretical method

The theoretical description is provided in detail in Ref.[25]. Here, the main parts of the model are emphasized.The version of the IQMD model is the IQMD-BNU (Beijing Normal University).

The nuclear interaction is described by the potential energy density of asymmetric nuclear matter with density ρ and asymmetry δ:

where ρ0＝0.16 fm-3is the normal density; p and p′are the momenta of the nucleon;and fτ（r，p）is the phase-space density, with τ ＝1／2 for the neutrons and τ ＝-1／2 for the protons.The parameters α=- 75.86 MeV,β ＝166.43 MeV, γ ＝1.226, Csp= 38.06 MeV, x ＝0, Cm= -88.21 MeV,and Λ=664.86 MeV/c are temperature-independent,while the phase-space density is temperature-dependent.The parameter γi= 0.35, 0.75, and 2 will be applied to the corresponding slope of the symmetry energy at normal density of 45.12, 67.96, and 139.32 MeV, respectively.

Binary nucleon–nucleon (NN) collisions depend on the differential cross section of NN collisions. The cross section of NN collisions in free space was taken from

The evolution by the IQMD model is stopped, and the GEMINI code is switched on when the excitation energy of the prefragments is less than a specified value Estop = 3 MeV/nucleon. In our previous work (e.g., Ref. [25]), the influence of Estop on various observables was studied. It was found that the multiplicity of the final fragments does not depend on the chosen Estop.As shown in Ref. [26],the average fragment multiplicity as a function of the excitation energy manifests the onset of multifragmentation,which occurs at excitation energies of approximately 3 MeV/nucleon. Thus, Estop ＝3 MeV/nucleon was chosen in this work. In the GEMINI code, only the light-particle evaporation mode was considered.The excitation energy of the fragments can be calculated as

3 Results and discussion

To study the isospin fractionation in the projectile fragmentations, it is helpful to compare the isospins between the gas phase (fragments of Z＜3) and liquid phase(intermediate mass fragments and residues).For each reaction, more than 100,000 events were considered to ensure sufficient statistical accuracy. The measured free neutrons and protons came not only from the fragmenting stage, but also from the pre-equilibrium and secondary decay stages. The method proposed by Xu et al. was applied to estimate the ratio between the neutrons and protons emitted by the fragmenting source [8]. The yield ratio M(Z, N)/M(N, Z) as a function of the binding energy difference △B for light mirror nuclei is given by

where T is the temperature of the fragmenting system, and fZ，N（T）／fN，Z（T） is the secondary decay factor. Figure 1 shows the analysis of the yield ratios of four pairs of mirror nuclei:3H/3He,7Li/7Be,11B/11C, and15N/15O, obtained from the projectile fragmentations in124Sn,107Sn ＋120Sn at 600 MeV/nucleon with b = 5 fm. For the IQMD calculations,there was no secondary decay;hence,no correction Fig.1 The yield ratios of the mirror nuclei produced in the projectile fragmentations of124Sn,107Sn ＋120Sn at 600 MeV/nucleon with b ＝5 fm was required (fZ，N（T）／fN，Z（T）＝1). Fitting the yield ratios from the IQMD simulations with Eq. (5), we obtained ρn／ρp= 1.62 for124Sn, and 1.23 for107Sn. The larger ρn／ρpvalues than the projectile N/Z indicate the isospin fractionation. For calculations using IQMD＋GEMINI, the yield ratios display a trend similar to Eq. (5). Only three pairs of mirror nuclei were analyzed because of the scarce production of15N and15O.Fitting the yield ratios obtained from IQMD＋GEMINI using Eq. (5), the ρn／ρp= 1.32 for124Sn, and 1.03 for107Sn. In simulations with secondary decay,the extracted values of ρn／ρpwere smaller than the projectile N/Z(1.48 for124Sn,and 1.14 for107Sn)and those without secondary decay.This indicates that the secondary decay distorts the ρn／ρpvalues. In other words, fitting the final fragments using Eq. (5),the extracted ρn／ρpdoes not reflect the isospin of the gas phase in the fragmentation.We also found that the momentum-dependent potential energy affects the yield ratio M(Z, N)/M(N, Z), and the extraction of ρn／ρp.

Using ρn／ρpas the isospin observable, the symmetry energy dependence of the isospin fractionation in the projectile fragmentations of124Sn,107Sn ＋120Sn at 600 MeV/nucleon as well as136Xe,124Xe ＋208Pb at 1000 MeV/nucleon was studied. The results are shown in Figs. 2 and 3. It was found that the stiff symmetry energy results in a lower ρn／ρpvalue than the soft one.Because of the smaller isospin asymmetries, this symmetry energy dependence is relatively weak for107Sn and124Xe. The projectile-like source was heated but not compressed during the preequilibrium stage. It is supposed that the nuclear gas (light particles) was steamed out from the periphery of the source,where the density is lower than the normal density.The smaller symmetry energy (γi＝2) at subnormal Fig.3 (Color online)Same as Fig. 2 but for136Xe,124Xe ＋208Pb at 1000 MeV/nucleon density provides a smaller repulsive force for neutrons and an attractive force for protons in the neutron-rich system.This leads to a smaller ρn／ρpin the gas phase (Figs.2 and 3), leaving a neutron-rich liquid phase.

Figures 2 and 3 depict the impact parameter dependence of the isospin fractionation. The impact parameter dependence of the ρn／ρpvalues is weak within the uncertainties for124Sn and136Xe. However, for107Sn and124Xe, the general trend is a decreasing ρn／ρpwith increasing impact parameter. The107Sn is the neutron-deficient nucleus. It is known that124Xe is the stable nucleus, whose neutron number is the lowest in stable Xe nuclei. The isospin fractionation of the neutron-deficient system leads to a smaller value of ρn／ρpin the gas phase than that of the initial N/Z value. The anomalous ρn／ρpvalues may be related to the isospin diffusion between the projectile and target.

The ratio between mirror nuclei yields is an isospin observable. Figure 4 shows the yields and ratios of two pairs of mirror nuclei,3H/3He,7Li/7Be, produced in the projectile fragmentations of124Sn ＋120Sn at 600 MeV/nucleon. The calculations with and without the GEMINI code were compared. The yields of the prefragments (by IQMD) contain information on the isospin fractionation.However, some information will be lost during the secondary decay stage where the binding energies rule. The binding energies depend on the symmetry energy at the normal density rather than at a subnormal density. In the decay process,two factors affect the yields of3H and3He.First, the heavier nuclear decay becomes3H and3He.Second, excited3H and3He decay into lighter nuclei,reducing3H and3He. The calculation results suggest that the latter plays a dominant role.This is the reason why the calculated value of IQMD＋GEMINI is lower than that of IQMD,as shown in Fig.4.We discussed this phenomenon in Fig. 8(d) of Ref. [25]. In the figure, it is shown that the yields of the3H,3He, and7Be nuclei are reduced by secondary decay. Owing to the stronger bond for the7Li nucleus, its yields do not change significantly after the secondary decay. As a result, the secondary decay effects of the yield ratio are weak for3H/3He, but strong for7Li/7Be (as shown in Fig. 4c, d). In addition, as mentioned in Ref.[15], the momentum-dependent potential energy also affects the yield of light particles (3H,3He,7Li,7Be) as well as the yield ratio.

The yield ratios of the mirror nuclei3H/3He were used to study the dependence of the projectile fragmentation on symmetry energy.Two soft(γi=0.35 and 0.75)and stiff(γi= 2) symmetry energies were applied in the IQMD＋GEMINI model to perform the collisions of124Sn,107Sn ＋120Sn at 600 MeV/nucleon. In Fig. 5, the yield ratio calculations of3H/3He are shown as a function of the normalized bound charge Zbound/Zpas well as the impact parameter b. The Zbound is defined as the sum of the fragment charges of all fragments with Zi[34]. The bound charge monotonically increased with increasing impact parameter and was applied for event sorting in the experiment. The symmetry energy dependence and impact parameter dependence of ρn／ρp(Figs. 2 and 3) were also determined for the yield ratios of3H/3He. First, the stiff symmetry energy resulted in smaller yield ratios of3H/3He compared with the soft cases. The symmetry energy dependence is weak for107Sn, but relatively strong for124Sn. Second, owing to the coexistence of isospin diffusion and isospin fractionation, the yield ratios of3H/3He for107Sn decreased with increasing b, but those for124Sn are similar for all b.Compared with the isospin observable ρn／ρp,the yield ratios of3H/3He have two advantages,i.e.,smaller uncertainty and weak secondary decay effect. The yield ratios of3H/3He for both107Sn and124Sn decreased with increasing normalized bound charge Zbound/Zp.However, the difference between these two fragmenting systems is more evident in peripheral collisions. We can still see the effect of different symmetry energies within the range of error bar coverage.

The secondary decay effects can be further minimized by using the double ratio of3H/3He ratios between two systems with different isospin asymmetry of projectiles.The double yield ratio (DR) is defined as

where the subscripts 124 and 107 refer to the124Sn ＋120Sn and107Sn ＋120Sn collisions,respectively.For136Xe,124Xe＋208Pb collisions at 1000 MeV/nucleon,a similar DR was defined,but with subscripts of 136 and 124.Figures 6 and 7 show the calculated DR for projectile fragmentations in124Sn,107Sn ＋120Sn collisions at 600 MeV/nucleon, and136Xe,124Xe ＋208Pb collisions at 1000 MeV/nucleon.The DR values increased with increasing b and Zbound/Zp. A comparison of the DR values obtained from the calculations with different symmetry energies indicates that DR can be used as an isospin observable. The isospin observable DR allows one to constrain the symmetry energy by the peripheral collision at hundreds of MeV/nucleon.Using the bound charge in event sorting, the symmetry energy dependence is stronger in the region of a large bound charge. The symmetry energy also affects the error bar coverage.

To further test our theory,136Xe,124Xe ＋208Pb collisions are shown in Fig. 7. As mentioned above, isospin diffusion leads to a weak target dependence of the fragmenting source during peripheral collision. Thus, it is expected that the DR value will increase (i.e., a larger difference between the neutron-rich and neutron-deficient systems) with increasing impact parameter. Because the number of events is proportional to b2, the peripheral events are greater than the central events. The stiff symmetry energy results in a smaller ρn／ρpvalue in the gas phase compared with the soft case.

4 Conclusion

The isospin fractionations in124Sn,107Sn ＋120Sn at 600 MeV/nucleon, and136Xe,124Xe ＋208Pb at 1000 MeV/nucleon were investigated by the IQMD model coupled with the statistical code GEMINI. The yield ratio as a function of the binding energy difference for light mirror nuclei3H/3He,7Li/7Be,11B/11C, and15N/15O was applied to estimate the ratio between neutrons and protons in the gas of the fragmenting system.By comparing the estimated values obtained from the calculations with and without the GEMINI code, it was found that the secondary decay distorts the signal of the isospin fractionation. Using the neutron/proton ratio in the gas as the isospin observable,the symmetry energy dependence of the isospin fractionation was studied. The symmetry energy dependence was weak for107Sn and124Xe fragmentations, but evident for124Sn and136Xe fragmentations.

Because of the secondary decay effects, the neutron/proton ratio in the gas is not a robust isospin observable.Instead,the yield ratio of the light mirror nuclei3H/3He as well as its double ratio between two systems with different isospin asymmetries of the projectiles was considered.As a result of the small binding energy difference,the yield ratio of3H/3He was slightly distorted by the secondary decay.Distortions caused by secondary decay can be minimized by using the double ratio between two systems with different isospin asymmetries of the projectiles. The DR can be used as an isospin observable. Stiff symmetry energy generally leads to a smaller value of the gas, in contrast to the soft case. It is beneficial to obtain high-precision data because the counting for peripheral collisions is larger than for central collisions.