Directed flow of isospin sensitive fragments within a modified clusterization algorithm in heavy-ion collisions

Directed flow of isospin sensitive fragments within a modified clusterization algorithm in heavy-ion collisions 1 2 first-author 2 corresp Corresponding author. E-mail address: ygma@sinap.ac.cn Corresponding author. E-mail address: ygma@sinap.ac.cn ygma@sinap.ac.cn Amity Institute of Nuclear Science and Technology, Amity University, Noida, Uttar Pradesh-201303, India Amity Institute of Nuclear Science and Technology, Amity University, Noida, Uttar Pradesh-201303, India Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China By extending the minimum spanning tree (MST) clusterization algorithm for the binding energy cut, the isospin asymmetry dependence of directed flow for isospin sensitive isobar pairs (neutrons-protons, 3 H- 3 He) is studied from low towards high incident energies. The modified clusterization method (MSTB) has the advantage to identify the fragments at quite early time. It enhances (reduces) the production of free nucleons (fragments) over MST method. The directed flow of isobaric pair 3 H- 3 He is more sensitive towards isospin asymmetry caused by MSTB than isobaric pair n - p . This sensitivity becomes quite strong towards the high incident energy and neutron-rich reaction system. In conclusion, the inclusion of binding energy in clusterization method for the flow studies has been uniquely important for understanding the isospin physics, especially for high density behavior of symmetry energy. Binding energy Isospin Physics Symmetry energy Directed flow 1 Introduction 1 With the availability of radioactive beam facilities, isospin physics has attracted the whole of the nuclear community in last decade. In intermediate energy heavy-ion collisions, the dynamical models such as Quantum Molecular Dynamical (QMD) [ 1 1 , 2 2 , 3 3 , 4 4 , 5 5 , 6 6 , 7 7 , 8 8 , 9 9 , 10 10 , 11 11 , 12 12 , 13 13 , 14 14 , 15 15 ] , Boltzmann-Uehling-Uhlenbeck (BUU) [ 16 16 , 17 17 , 18 18 , 19 19 , 20 20 , 21 21 , 22 22 ] and statistical models [ 23 23 , 24 24 , 25 25 , 26 26 , 27 27 ] were used to get the rich information about isospin physics. The dynamical models are of unique importance due to their ability to follow the time evolution from initialization through compression and expansion to equilibrium state. It is also well known that dynamical models do not simulate the fragments. In order to simulate the fragments, one needs the secondary algorithm for coupling with dynamical models. In the efforts to reproduce the experimental data, many secondary algorithms have been developed to study the nuclei near the drip line [ 11 11 , 12 12 , 13 13 , 15 15 , 28 28 , 29 29 ] , but very few for nuclei away from drip line [ 5 5 , 6 6 , 8 8 , 10 10 , 15 15 ] . The most commonly and widely used algorithm depends on the spatial and momentum coordinates of the nucleons, is known as minimum spanning tree (MST) algorithm [ 12 12 , 15 15 ] . According to this method, two nucleons undergo the cluster formation if the relative distance (| R i – R j |) and relative momentum (| P i – P j |) between the nucleons was less than 3.5–4 fm and 250–268 MeV/ c , respectively. These parameters can be obtained by fitting the experimental data for some of the global observables such as multiplicity of intermediate mass fragments [ 12 12 , 15 15 ] with theoretical results. Recently, MST method was further extended to isospin-dependent MST, in which the cut on the momentum space was kept same, but the cut on the spatial coordinates was constraint on the basis of type of particles. The distance between the different kind of particles was taken as follow: <math id="M1"><mrow><mo>|</mo><msubsup><mi>R</mi><mi>i</mi><mi>p</mi></msubsup><mo>−</mo><msubsup><mi>R</mi><mi>j</mi><mi>p</mi></msubsup><mo>|</mo></mrow></math> = 3 fm, <math id="M2"><mrow><mo>|</mo><msubsup><mi>R</mi><mi>i</mi><mi>n</mi></msubsup><mo>−</mo><msubsup><mi>R</mi><mi>j</mi><mi>p</mi></msubsup><mo>|</mo><mo>=</mo><mo>|</mo><msubsup><mi>R</mi><mi>i</mi><mi>n</mi></msubsup><mo>−</mo><msubsup><mi>R</mi><mi>j</mi><mi>n</mi></msubsup><mo>|</mo></mrow></math> = 6 fm [ 14 14 ] . Since the MST method only depends on the constraints from position and momentum, it seems to be worried about the stability of fragments due to the formation of artificial weakly bound fragments. To avoid from this problem, more complicated methods like Stimulated Annealing Clusterization Algorithm (SACA) [ 13 13 , 28 28 ] , Early Cluster Recognization Algorithm (ECRA) [ 29 29 ] were also developed. These methods were found quite successful, but, very complicated. Moreover, due to the choice of parameters such as cooling parameters, iterations procedure, and choice of minimal can follow the totally different configuration of fragments. Due to the limitation of sharp minimal, in the mildly excited or asymmetric systems, the scope of these two methods was found to restrict at some points. The possible best method for avoiding from the artificial formation of fragments is found a way to constrain the fragments by using an average binding energy cut of 4 MeV/nucleon [ 30 30 ] . One can further improve the method by using the realistic binding energy rather than 4 MeV/nucleon [ 30 30 ] . This method was found to be as simple as MST and found to reproduce the experimental data just like the complicated methods SACA and ECRA etc. The extensive studies for the multiplicity of different kind of fragments with different clusterization methods have been performed [ 12 12 , 13 13 , 14 14 , 15 15 , 28 28 , 30 30 ] , while flow parameters are poorly known with this type of methods [ 1 1 , 8 8 , 11 11 , 12 12 , 31 31 , 32 32 ] . In last two decades, the collective flow of nucleons with non-central collisions was found as a useful tool for extracting the nuclear equation of state (NEOS) of symmetric nuclear matter [ 31 31 , 32 32 ] . With the passage of time, efforts have been done to extract the information of isospin physics in term of symmetry energy using the directed as well as elliptic flow [ 8 8 , 21 21 , 32 32 ] . The studies showed that the experimental results of directed flow can be reproduced by (1) changing the mean field or nucleon-nucleon ( NN ) cross sections in the transport model (2) using the nucleon phase space for calculations [ 31 31 , 33 33 , 34 34 ] . Most of the studies were found to obey the hard equation of state [ 33 33 , 34 34 , 35 35 ] , which opposed the findings with other observables in the literature [ 12 12 , 19 19 , 28 28 , 31 31 , 32 32 , 36 36 , 37 37 , 38 38 ] . In recent study, the soft equation of state was predicted [ 39 39 ] , but, the comparison was made between the directed flow extracted from phase space of nucleons and the directed flow of fragments having Z≤2 [ 40 40 ] . Within the BUU calculations [ 21 21 ] , the importance of directed flow of isospin sensitive fragments n , p , 3 H and 3 He over the directed flow of nucleon phase space has already been put forwarded. The elliptic flow extracted with the simple MST method was also found to be quite sensitive towards the determination of symmetry energy at supra-saturation densities [ 8 8 ] . The results indicated that the directed as well as elliptic flow can be a useful tool to understand the isospin physics or symmetry energy. The main problem with the study was that isospin physics with directed and elliptic flow is either studied from the nucleon phase space or by using the simple MST method. It is worth mentioning here that one cannot guarantee about the stability of fragments formed with simple MST method. Due to the worry about the stability of fragments, it is the prime need of the present time to study the importance of Isospin physics from directed and elliptic flow with the stable fragments. In the present study, first the stable fragments are formed by applying the binding energy cut on the pre-clusters formed with the simple MST method. Afterward, the sensitivity of the isospin sensitive fragments n , p , 3 H and 3 He directed flow is studied from low towards the high incident energy with the variation in isospin asymmetry of the reaction systems. The modified clusterization method is explained in Section-2, the results and discussions are presented in Section-3, followed by the conclusion in Section-4. 2 Methodology 1 In the present work, IQMD model is used, which is discussed in detail in our recent publications [ 5 5 , 6 6 , 10 10 ] , originally developed by Hartnack and Co-workers [ 11 11 , 12 12 , 13 13 ] . The model is modified by the authors for the density dependence of symmetry energy, having form: <math display="block" id="M3"><mrow><msub><mi>E</mi><mrow><mtext>Sym</mtext></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>C</mi><mrow><mi>s</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow><mn>2</mn></mfrac><msup><mrow><mo stretchy="false">(</mo><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>+</mo><mfrac><mrow><msub><mi>C</mi><mrow><mi>s</mi><mo>,</mo><mi>p</mi></mrow></msub></mrow><mn>2</mn></mfrac><msup><mrow><mo stretchy="false">(</mo><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><msub><mi>γ</mi><mi>i</mi></msub></mrow></msup></mrow></math> with the parameters of C s,k = 25 MeV and C s,p = 35.2 MeV. When we set γ i =0.5 and 1.5, respectively, it corresponds to the soft and stiff symmetry energy [ 5 5 , 6 6 ] . Minimum Spanning Tree method with Binding Energy Check (MSTB) is used. This method is a modified version of the normal MST and old MSTB method. The difference between old MSTB and this MSTB is the inclusion of energy from momentum dependent interactions as well as symmetry energy along with Skyrme interactions. The procedure is as follow: the phase space obtained from IQMD is analyzed with simple MST method and pre-clusters are sort out. Since we are not aware about the stability of pre-clusters formed at this stage, the pre-clusters formed from the simple MST are now subjected to the binding energy condition as follow: <math id="M4"><mrow><mfrac><mn>1</mn><mrow><msub><mi>N</mi><mi>f</mi></msub></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>a</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mi>f</mi></msub></mrow></munderover><mrow><mrow><mo>[</mo><mrow><msqrt><mrow><msup><mrow><mo stretchy="false">(</mo><msub><mi>p</mi><mi>α</mi></msub><mo>−</mo><msub><mi>P</mi><mrow><msub><mi>N</mi><mi>f</mi></msub></mrow></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mi>α</mi><mn>2</mn></msubsup></mrow></msqrt><mo>−</mo><msub><mi>m</mi><mi>α</mi></msub><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>β</mi><mo>≠</mo><mi>α</mi></mrow><mrow><msub><mi>N</mi><mi>f</mi></msub></mrow></munderover><mrow><msub><mi>V</mi><mrow><mi>α</mi><mi>β</mi></mrow></msub></mrow></mstyle></mrow><mo>]</mo></mrow></mrow></mstyle><mo><</mo><mo>−</mo><msub><mi>E</mi><mrow><mtext>Bind</mtext></mrow></msub><mo>.</mo></mrow></math> Here, we take E Bind = 4.0 MeV/nucleon if N f ≥3 and E Bind = 0 otherwise. In this equation, N f is the number of nucleons in a fragment, <math id="M5"><mrow><msub><mi>P</mi><mrow><msub><mi>N</mi><mi>f</mi></msub></mrow></msub></mrow></math> is the average momentum of the nucleons bound in the fragment. The requirement of a minimum binding energy excludes loosely bound fragments which will decay later. The realistic value of E Bind changes slightly the fragment multiplicity at intermediate times, but has no influence on the qualitative behavior and on the asymptotic results. However, if by using the realistic binding energy, one searches for the most bound configuration, the results will get affected. At the present time, we had just focused on the bound configurations and hence the average cut of binding energy –4 MeV/nucleon is justified. 3 Results and discussions 1 Several thousands of event for the reactions of 112 Sn+ 112 Sn and 124 Sn+ 124 Sn between the incident energy 50 and 600 MeV/nucleon for the nearly central collisions using the IQMD model coupled with MST and MSTB algorithms are simulated. The soft momentum dependent equation of state with soft symmetry energy and isospin-energy dependent cross sections is employed. There are two methods in the literature used to calculate the directed flow. In the first case, the directed flow is extracted from the mid-rapidity slope of the P x / A versus rapidity distribution Y c.m. / Y beam plots. The rapidity distribution is calculated as follow: <math id="M6"><mrow><mi>Y</mi><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>ln</mi><mfrac><mrow><mi>E</mi><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>p</mi><mi>z</mi></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow><mrow><mi>E</mi><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>−</mo><msub><mi>p</mi><mi>z</mi></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></mfrac><mo>,</mo></mrow></math> where E ( i ) and p z ( i ) are the total energy and longitudinal momentum of i th particle, respectively. The second method is to study the directed transverse in-plane flow <math id="M7"><mrow><msubsup><mi>P</mi><mi>x</mi><mrow><mtext>dir</mtext></mrow></msubsup></mrow></math> is as follow: <math id="M8"><mrow><mrow><mo>〈</mo><mrow><msubsup><mi>P</mi><mi>x</mi><mrow><mtext>dir</mtext></mrow></msubsup></mrow><mo>〉</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>A</mi></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mi>i</mi><mi>A</mi></munderover><mrow><mtext>sign</mtext><mo>{</mo><mi>Y</mi><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>}</mo></mrow></mstyle><msub><mi>P</mi><mi>x</mi></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>,</mo></mrow></math> where Y ( i ) is the rapidity distribution as discussed above and P x ( i ) is the transverse momentum of the i th particle in x-direction. This <math id="M9"><mrow><msubsup><mi>P</mi><mi>x</mi><mrow><mtext>dir</mtext></mrow></msubsup></mrow></math> is defined over entire rapidity region and therefore expected to present an easier way of measuring the in-plane flow rather than complicated P x / A plots. In order to check the stability of fragments and effect of isospin physics, in Fig.1 Fig.1 , the time evolution of free nucleons, light charged particles (LCPs) and intermediate mass fragments (IMFs) for neutron-poor and neutron-rich reaction systems at 50 MeV/nucleon with MST and MSTB algorithms is displayed. In the high density phase, the MSTB algorithm does not find any fragments with reasonable binding energy and hence most of the particles are free, while, there were a lot of artificial fragments with MST. After the high density phase is over, it starts recognizing the fragments, which are real, bound and stable. In all the cases, MSTB helps to identify the fragments quite early. From the figure, it is clear that MSTB enhances the production of free particles and reduces the production of LCPs and IMFs. Moreover, with increase in the size of the fragment, MST takes less time to match with the MSTB results and also the difference between MST and MSTB results goes on decreasing throughout the time evolution. When the results are compared for neutron-poor and neutron-rich reaction systems, it is found that enhanced production of fragments takes place for the more neutron-rich system. The enhanced production is quite effective in free nucleons and LCPs compared to IMFs. It is also shown by us recently that free nucleons and LCPs are more sensitive towards the symmetry energy compared to IMFs [ 10 10 ] . In this regard, the sensitivity of free nucleons and LCPs with MSTB towards the neutron-rich system can also play an important role in the prediction of symmetry energy, which is a topic of separate discussion. After assessing the sensitivity, for the further study, the isospin sensitive free nucleons and LCPs, namely, n , p , 3 H and 3 He are used. Using the first method, we have plotted in Fig.2 Fig.2 , the rapidity distribution of P x / A for neutrons, protons and 3 H particles at 100 MeV/nucleon for the central collisions of 124 Sn reaction system. The mid-rapidity slope of P x / A represents the directed flow. The flow is found to be affected only for the heavier fragments, neutrons and protons flow is least affected. This is due to the isotropic (anisotropic) distribution of free nucleons (LCPs) [ 37 37 ] . Naturally, yield has to be independent of the type of distribution, while momentum of particles is supposed to depend strongly on the type of distribution. One of the examples is that different behavior of flow is obtained in earlier studies for different distribution regions [ 42 42 , 43 43 ] . The sensitivity of flow with MSTB towards the isospin physics at higher incident energy is presented in the following sections. In order to study the role of MSTB at higher incident energies and towards isospin physics, in Fig.3 Fig.3 , the isospin asymmetry dependence of more reliable flow quantity <math id="M10"><mrow><msubsup><mi>P</mi><mi>x</mi><mrow><mtext>dir</mtext></mrow></msubsup></mrow></math> is plotted for isobaric pairs n , p ( 3 H, 3 He) with MST and MSTB algorithm from 50 to 600 MeV/nucleon. In Fig.3 Fig.3 , with increase in the incident energy, isospin asymmetry dependence of directed flow for n, p is found to affect more by the method of clusterization in comparison with the type of particles. In all panels, protons (neutrons) have more (less) positive directed flow with MST as well as MSTB method. This is due to the contribution of Coulomb repulsion in protons production. The directed flow of protons as well as neutrons becomes more positive with MSTB over MST towards the high incident energy. This is true because there was increased content of unstable fragments with MST towards the high incident energy. The MSTB will break these unstable fragments into free particles (as per the condition specified earlier) and hence dominance of collisions will take place. Moreover, higher positive flow of neutrons towards the high incident energy can act as a tool for isospin dependent cross section. Furthermore, these findings are suggesting a shift in the balance energy with different kind of fragments. In Fig.4 Fig.4 , just like n and p , the isospin asymmetry dependence of directed flow for 3 H, 3 He is found to be affected more by the method of clusterization compared to type of particles towards the high incident energy. Up to 100 MeV/nucleon, 3 He ( 3 H) particles have more (less) positive directed flow, but, at and after 200 MeV/nucleon, 3 H ( 3 He) particles have more (less) positive flow with MST as well as MSTB methods. This is due to the dominance of mean field (collisions) at and below (at and above) 100 (200) MeV/nucleon. When the mean field (collisions) dominates, Coulomb interactions (isospin dependent cross sections) make the flow more positive for proton rich (neutron-rich) 3 He ( 3 H) particles. Towards the higher incident energies, the flow of 3 H as well as 3 He becomes more positive with MST over MSTB method. As with the inclusion of binding energy restrictions in clusterization algorithm, most of the 3 H and 3 He were found unstable. Especially, this is true for 3 He and hence directed flow is decreased with MSTB, which has resulted in the increase of the transverse directed flow for neutrons and protons (as discussed before in Fig.3 Fig.3 ). In brief, the directed flow of isospin sensitive particles n , p , 3 H, 3 He is affected by the different method of clusterization. The role of different clusterzation method is dominating for the isobaric pair ( 3 H, 3 He) compared to ( n , p ). Moreover, the behavior of directed flow for ( n , p ) pair is smooth with incident energy, but, the transition is observed for the directed flow of ( 3 H, 3 He) pair with incident energy. The isobaric pair ( n , p ) and ( 3 H, 3 He) has opposite sensitivity with the type of method of clusterization. These points suggest that the methods of clusterization are going to play an important role in understanding the isospin physics and hence need to be handled carefully during the determination of symmetry energy also. 4 Conclusion 1 The modified clusterization method has a time saving advantage in terms of providing the real, bound and stable fragments at quite early time. It has the capability to obtain the results similar to the complicated secondary algorithms SACA and ECRA. Moreover, the complicated methods have the drawback over MSTB of time consuming and consists a lot of adjustable parameters. The directed flow of isobaric pair 3 H- 3 He in comparison to n - p is found to be strongly dependent on the isospin asymmetry caused by MSTB with the increase of incident energy as well as isospin of the reaction systems. This indicates it as a robust approach for understanding the high density behavior of isospin physics or symmetry energy. References Gautam S , Kumari R , Puri R K . Phys Rev C , 2012 , 86 : 034607 . Phys Rev C Jain A , Kumar S , Puri R K . Phys Rev C , 2012 , 85 : 064608 . Phys Rev C Ma Y G and Shen W Q . Phys Rev C , 1995 , 51 : 710 - 715 . Phys Rev C Ma Y G , Wei Y B , Shen W Q , et al . Phys Rev C , 2006 , 73 : 014604 . Phys Rev C Kumar S , Ma Y G , Zhang G Q , et al . Phys Rev C , 2011 , 84 : 044620 . Phys Rev C Kumar S , Ma Y G , Zhang G Q , et al . Phys Rev C , 2012 , 85 : 024620 . Phys Rev C Yan T Z , Ma Y G , Cai X Z , et al . Phys Lett B , 2006 , 638 : 50 - 54 . Phys Lett B Russotto P , Wu P Z , Zoric M , et al . Phys Lett B , 2011 , 697 : 471 - 476 . Phys Lett B Wei Y B , Ma Y G , Shen W Q , et al . Phys Lett B , 2004 , 586 : 225 - 231 . Phys Lett B Kumar S and Ma Y G . Nucl Phys A , 2013 , 898 : 59 - 77 . Nucl Phys A Hartnack C , Puri R K , Aichelin J , et al . Eur Phys J A , 1998 , 1 : 151 - 169 . Eur Phys J A Aichelin J . Phys Rep , 1991 , 202 : 233 - 360 . Phys Rep Puri R K and Aichelin J . J Comput Phys , 2000 , 162 : 245 - 266 . J Comput Phys Zhang Y , Li Z X , Zhou C S , et al . Phys Rev C , 2012 , 85 : 051602 . Phys Rev C Singh J and Puri R K . Phys Rev C , 2000 , 62 : 054602 . Phys Rev C Ma Y G , Shen W Q , Zhu Z Y . Phys Rev C , 1995 , 51 : 1029 - 1032 . Phys Rev C Coupland D D S , Lynch W G , Tsang M B , et al . Phys Rev C , 2011 , 84 : 054603 . Phys Rev C Ma Y G , Shen W Q , Feng J , et al . Phys Rev C , 1993 , 48 : R1492 - R1496 . Phys Rev C Tao C , Ma Y G , Zhang G Q , et al . Nucl Sci Tech , 2013 , 24 : 030502 . Nucl Sci Tech Famiano M , Liu T , Lynch W G , et al . Phys Rev Lett , 2006 , 97 : 052701 . Phys Rev Lett Li B A . Phys Rev Lett , 2002 , 88 : 192701 . Phys Rev Lett Li B A , Chen L W , Ma H R , et al . Phys Rev C , 2007 , 76 : 051601(R) . Phys Rev C Ogul R , Botvina A S , Atav U , et al . Phys Rev C , 2011 , 83 : 024608 . Phys Rev C Zhang G Q , Cao X G , Fu Y , et al . Nucl Sci Tech , 2012 , 23 : 61 - 64 . Nucl Sci Tech Ma Y G . Phys Rev Lett , 1999 , 83 : 3617 - 3620 . Phys Rev Lett Fang D Q , Ma Y G , Zhong C et al . J Phys G , 2007 , 34 : 2173 - 2181 . J Phys G Ma Y G , Su Q M , Shen W Q , et al . Phys Rev C , 1999 , 60 : 024607 . Phys Rev C Vermani Y K and Puri R K , Euro phys Lett , 2009 , 85 : 62001 . Euro phys Lett Dorso C O and Randrup J . Phys Lett B , 1993 , 301 : 328 - 333 . Phys Lett B Goyal S and Puri R K . Phys Rev C , 2011 , 83 : 047601 . Phys Rev C Danielewicz P , Lacey R , Lynch W G . Science , 2002 , 298 : 1592 - 1596 . Science Reisdorf W , Leifels Y , Andronic A , et al . Nucl Phys A , 2012 , 876 : 1 - 60 . Nucl Phys A Molitoris J J , Hahn D , Stocker H . Nucl Phys A , 1986 , 447 : 13 - 26 . Nucl Phys A Sood A D and Puri R K . Phys Rev C , 2006 , 73 : 067602 . Phys Rev C Sood A D and Puri R K . Phys Rev C , 2004 , 70 : 034611 . Phys Rev C Hartnack C , Oeschler H , Aichelin J . Phys Rev Lett , 2006 , 96 : 012302 . Phys Rev Lett Vermani Y K and Puri R K . Nucl Phys A , 2010 , 847 : 243 - 252 . Nucl Phys A Zhang G Q , Ma Y G , Cao X G , et al . Phys Rev C , 2011 , 84 : 034612 . Phys Rev C Gautam S , Chugh R , Sood A D , et al . J Phys G , 2010 , 37 : 085102 . J Phys G Pak R , Li B A , Benenson W , et al . Phys Rev Lett , 1997 , 78 : 1026 - 1029 . Phys Rev Lett Kumar S , Kumar S , Puri R K . Phys Rev C , 2010 , 81 : 014601 . Phys Rev C Han L X , Ma G L , Ma Y G , et al . Phys Rev C , 2011 , 84 : 064907 . Phys Rev C Wang J , Ma Y G , Shen W Q , et al . Nucl Sci Tech , 2013 , 24 : 030501 . Nucl Sci Tech 10.13538/j.1001-8042/nst.2013.05.009 Supported by the Chinese Academy of Sciences Support program for young international scientists (No.2010Y2JB02) , National Natural Science Foundation of China (NSFC) projects (Nos.11220101005, 11035009 and 10979074) and the Knowledge Innovation project of the Chinese Academy of Sciences (No.KJCX2-E W-N01) . 2013-06-27 2013-10-01 2021-04-11 2013 Nuclear Science and Techniques 5 24